FormalPara Key Summary Points

Why carry out this study?

 The impact of basal insulin choice on the risks of all-cause and cardiovascular mortality has not been thoroughly investigated in a real-world setting for people with type 2 diabetes (T2D) from the UK.

 A previous real-world study in Finland observed a lower mortality risk associated with insulin detemir (detemir) than with insulin glargine (glargine), which they hypothesised might be explained by a difference in cardiovascular mortality.

What was learned from the study?

 In the present real-world study in the UK, there was a lower risk of all-cause mortality in insulin-naïve people with T2D who initiated detemir than in those who initiated glargine, after adjustment for potential confounding factors.

 There was also a non-significant trend towards a lower risk of cardiovascular mortality associated with detemir versus glargine, again after adjustment for potential confounders.

 Of potential interest, the lower risks of all-cause and cardiovascular mortality associated with detemir in comparison with glargine were most pronounced in a subset of people with body mass index ≥ 30 kg/m2.

Digital Features

This article is published with digital features, including a summary slide, to facilitate understanding of the article. To view digital features for this article go to https://doi.org/10.6084/m9.figshare.14153456.

Introduction

Poor glucose control in diabetes may result in an increased risk of complications, including microvascular and macrovascular disease [1], with an increase in early mortality [2]. In advanced type 2 diabetes (T2D), where non-insulin therapies are not able to maintain glycaemic control, basal insulin therapy is frequently added to achieve optimal glycaemic control with the aim to reduce the risk of both acute and chronic complications [3].

Both insulin detemir (detemir) and insulin glargine (glargine) are long-acting insulin analogues, administered subcutaneously via injection, with a duration of action of up to, or over, 24 h depending on the product and dose [4, 5]. Both insulins are a pharmacological advance compared to human insulins, but differ in their pharmacokinetic/pharmacodynamic profiles, resulting in differing propensity for causing weight gain, hypoglycaemia and periods of hyperglycaemia or variability in blood glucose levels [6, 7]. As these factors may all contribute to cardiovascular risk, it is possible that different basal insulins may affect prognosis differently. However, there is limited knowledge of the differences between the available basal insulins regarding their impact on the risk of cardiovascular disease (CVD) and early mortality.

In randomised controlled trials (RCTs) [8, 9] and short-term observational studies [10, 11], detemir has shown advantages in terms of less weight gain compared with glargine. To assess the effect of a clinical intervention on mortality, many participants would need to be followed for a long time in RCTs, with surrogate or composite endpoints often utilised instead [12, 13]. Furthermore, the more strict inclusion and exclusion criteria for trial enrolment (in comparison with real-world studies) mean that the results of RCTs tend to be less applicable to the general populations seen in clinical practice [14].

A long-term observational study conducted in Finland observed that detemir was associated with a significantly lower all-cause (29% lower) and cardiovascular (36% lower) mortality risk than glargine [11]. To the authors’ knowledge, this has been the only real-world study to date to compare the all-cause and cause-specific mortality risk in people with T2D treated with detemir in comparison with glargine. Moreover, these findings need to be assessed in other patient populations and with use of other data sources with different sources of potential bias to strengthen the evidence base for these observations and test their generalisability across different settings and patient populations [15].

Analyses of real-world data from large databases, such as the Clinical Practice Research Datalink (CPRD), can be useful to assess response to, and outcomes associated with, T2D treatments, and hence can augment current clinical trial research. Therefore, the aim of this register-based cohort study in the UK was to estimate the differences in all-cause and cardiovascular mortality rates between first-time users of the basal insulins detemir and glargine in a real-world population with T2D by utilising the CPRD GOLD database.

Methods

The study protocol was approved by the Independent Scientific Advisory Committee (ISAC) (Medicine and Healthcare Products Regulatory Agency): approved study protocol number 19_058A.

Study Cohort

The study cohort was identified in the CPRD GOLD national database: a large, primary care database that includes data from people registered at general practices across the UK (England, Wales, Scotland and Northern Ireland). The database contains information on people’s demographics, medical diagnoses and symptoms, prescribing information and lifestyle factors, including alcohol consumption and smoking. Individuals in the database can be linked to secondary care and other health datasets, enabling a wider view of patient care [16].

Eligible people in the database were those with T2D, who were insulin-naïve, aged ≥ 40 years and who initiated treatment with either detemir or glargine between 24 June 2004 and 10 May 2019, in accordance with local guidelines (Fig. 1). People with T2D were defined as those with a record of a read code for T2D (Appendix 1) and without a record of a read code for T1D (at any time) in the database (Appendix 2). People who initiated detemir (Appendix 3) and glargine (Appendix 4) were defined as those with a record of a prescription for the appropriate insulin at or after 24 June 2004 (the date on which both drugs were available in the database and thus were available for prescription in the UK). Eligible people who initiated glargine could receive either the 100 units/ml (U100) or 300 units/ml (U300) formulation. Insulin-naïve people were defined as those without a previous record of a prescription for any insulin (Appendix 5) prior to their prescription for detemir or glargine. Eligible people were those who, at least 1 year prior to initiating detemir or glargine, had first registered with the practice and had acceptable research-quality data (as defined by the CPRD [17]) and the practice was of acceptable research quality (as defined by the CPRD [17]).

Fig. 1
figure 1

Study cohort disposition. CPRD Clinical Practice Research Datalink, detemir insulin detemir, glargine insulin glargine, ONS Office for National Statistics, T1D type 1 diabetes, T2D type 2 diabetes

Outcomes

All-cause mortality was defined as mortality from any cause and was identified in the CPRD GOLD database through a registration of mortality date. Cardiovascular mortality was defined as mortality from CVD and was identified through linkage between the CPRD GOLD database and the Office for National Statistics (ONS) death registration data, which includes information on the official date and cause(s) of death [16]. Specifically, mortality from CVD was identified as a record with the 10th revision of the International Statistical Classification of Diseases and Related Health Problems (ICD-10) codes I00-I99 specified as the underlying cause of death (Supplementary Table S1). Linkage was not available for the total study cohort but was available for practices based in England and Wales that had consented to participate in the linkage scheme.

Follow-up Time

The index date was defined as the first prescription date for detemir or glargine. End of follow-up was defined as the earliest date the patient transferred out of the practice or the date of the last collection for the practice in the CPRD GOLD database, or the end of ONS coverage (14 January 2019). Total follow-up time was calculated as the time (in years) from the index date until death or end of follow-up, whichever came first.

Baseline Characteristics

Socioeconomic characteristics were obtained by linking the CPRD GOLD database to small area-level databases (patient and practice level of index of multiple deprivation [IMD2015]). Data for all other baseline characteristics were obtained from the CPRD GOLD database. These included the following characteristics: age, sex (male, female), calendar year of index date (2004–2019), diabetes duration, smoking status (current, former, never) and glycated haemoglobin (HbA1c). Body mass index (BMI) was calculated from weight and height (at age ≥ 18 years). Dyslipidaemia was defined as total cholesterol > 6.2 mmol/l or low-density lipoprotein cholesterol > 6.2 mmol/l or high-density (HDL) cholesterol ≤ 1 mmol/l for men or HDL cholesterol ≤ 1.2 mmol/l for women or triglyceride ≥ 1.7 mmol/l. Hypertension was defined as systolic blood pressure ≥ 140 mmHg or diastolic blood pressure ≥ 90 mmHg. Kidney function was ascertained by estimated glomerular filtration rate (eGFR) calculated according to the Cockcroft-Gault formula: eGFR in ml/min/1.73 m2 = (140 − age in years) × weight in kg × 1.23 for men, 1.03 for women)/serum creatinine in µmol/l. The latest values prior to the index date were used for weight, height, HbA1c, dyslipidaemia and systolic and diastolic blood pressure. History of CVD, prior hypoglycaemia, nephropathy, retinopathy, neuropathy, chronic kidney disease and cancer were defined as a record of selected read codes prior to the index date (Appendix 6). Bolus insulin, non-insulin glucose-lowering medication and CVD-preventive medication were defined as records of selected product codes within 1 year prior to the index date (Appendix 7).

Statistical Analyses

Data import and pre-processing were performed in SAS Enterprise Guide 7.15. Statistical analyses were performed using R version 3.6.0 and the “survival” package version 2.44-1.1. Total basal insulin exposure during follow-up was further examined by plotting the years of cumulative exposure against the time since the index date. Crude mortality rates were calculated as the number of deaths divided by the total follow-up time and reported as rates per 100 person-years. Furthermore, cumulative event curves showed the probability of mortality versus the time since first prescription.

Main Analyses

Cox proportional hazards models were used to calculate crude and adjusted hazards ratios (HRs) with 95% confidence intervals (CIs) for all-cause and cardiovascular mortality in people treated with detemir relative to glargine. The main analyses were intention-to-treat-like analyses where people were allocated to the treatment group of their first prescription during the entire follow-up period. Age was used as the underlying timescale in the Cox proportional hazards models. Both the crude and adjusted models included the calendar year of the index date (2004–2019) as a categorical variable. Furthermore, the following pre-specified covariates that were expected to be the most likely potential confounders (Table 1) were included in the adjusted models: sex (male, female), diabetes duration (0–1, 1–2, 2–5, 5–10, > 10 years), BMI (categorical, based on World Health Organization classes), history of CVD (yes, no), CVD medication at baseline (antihypertensive [yes, no]; lipid lowering [yes, no]; other cardiovascular [yes, no]), sodium-glucose cotransporter-2 inhibitor and/or glucagon-like peptide-1 receptor agonist (time-varying: yes, no); other non-insulin glucose-lowering medications (time-varying: yes, no) and bolus insulin (time-varying: yes, no; data not summarised in Table 1).

Table 1 Baseline characteristics of the study cohort in relation to the model covariates

In addition to the pre-specified covariates, the fully adjusted models considered the following potential baseline covariates: socioeconomic status (level 1, 2, 3, 4 or 5), smoking status (yes, no), HbA1c (continuous), prior hypoglycaemia (yes, no), nephropathy (yes, no), retinopathy (yes, no), neuropathy (yes, no), chronic kidney disease (yes, no), dyslipidaemia (yes, no), hypertension (yes, no), eGFR (continuous) and history of cancer (yes, no). Due to missing data for some of the additional variables considered, potential covariates were selected using a forward one-step procedure and the complete case population (people without missing data for any of the variables). A covariate was selected for inclusion in the fully adjusted model if the HR of detemir versus glargine changed by > 10% on the natural log-scale. The fully adjusted model was then refitted to the total cohort.

The models were verified by plots of Martingale residuals. Since the results from the pre-specified adjusted and fully adjusted models were very similar, only the results from the crude and fully adjusted models are presented in the results section.

Sub-analyses

The main analyses were repeated for the following subsets of people with cardiovascular risk factors at the index date: those with obesity (BMI ≥ 30 kg/m2); those with duration of diabetes ≥ 5 years; those treated with CVD medication and/or with prior CVD; those with HbA1c above the median level for the total cohort.

Sensitivity Analyses

To investigate the sensitivity of assumptions in the statistical models, the main analyses were repeated for four additional ‘sensitivity’ analyses: (1) people were censored at the time that their prescription data indicated a switch to a different basal insulin from their first prescription during follow-up; (2) all people were censored after 8 years of follow-up; (3) only those with an index date post-2006 were included in the analysis; (4) a time-varying exposure model was used, which allowed for interruptions to basal insulin use. In sensitivity analysis 4, a person was assigned to a period of their ‘current prescription’ if at least one prescription was issued within the previous year; otherwise, the period was defined as ‘no current prescription’, and people were censored if they switched to a different basal insulin from their first prescription. The HRs (detemir/glargine) in sensitivity analysis 4 were estimated for periods with the ‘current prescription’.

Results

Study Cohort

The total cohort comprised 12,847 people (detemir, n = 3031, and glargine, n = 9816) (Fig. 1). Of the total cohort, 53.7% (6897 people) were eligible for linkage to the ONS for cardiovascular mortality data (detemir, n = 1833, and glargine, n = 5064). In the total cohort, the median (interquartile range [IQR]) age was 66.94 (57.50–76.09) years (Table 1). Median (IQR) diabetes duration was similar for people in both the detemir and glargine groups, 7.66 (4.28–11.56) vs. 7.61 (4.08–11.81) years, respectively. Median (IQR) follow-up time was 3.86 (1.74–6.72) vs. 3.59 (1.54–6.72) years for the detemir and glargine groups, respectively. After 5 years, the median (IQR) exposed time was 4.64 (1.93–5.00) and 5.00 (2.93–5.00) years for the detemir and glargine groups, respectively (Supplementary Fig. S1). Fewer than 15% of people switched to another basal insulin during follow-up (14.4% and 9.7% with detemir and glargine, respectively).

All-Cause and Cardiovascular Mortality

The cumulative event curves for all-cause mortality are shown in Fig. 2a. Among the 12,847 people included in the total cohort, there were 3231 deaths: 635 with detemir and 2596 with glargine. The cumulative event curves for cardiovascular mortality are shown in Fig. 2b. Of the 6897 people linked to the ONS, 528 died of cardiovascular causes, 117 of whom had been prescribed detemir and 411 glargine.

Fig. 2
figure 2

a All-cause and b cardiovascular mortality for detemir and glargine. CPRD Clinical Practice Research Datalink, detemir insulin detemir, glargine insulin glargine, ONS Office for National Statistics

Main Analyses

For all-cause mortality, the crude and fully adjusted HRs (95% CI) were 0.85 (0.77; 0.92) and 0.86 (0.79; 0.95), respectively, significantly in favour of detemir (Fig. 3a). For cardiovascular mortality, the crude and fully adjusted HRs (95% CI) were 0.82 (0.66; 1.01) and 0.83 (0.67; 1.03), respectively, with the non-significant numerical trends in favour of detemir (Fig. 3b).

Fig. 3
figure 3

Main analyses and sub-analyses for a all-cause and b cardiovascular mortality. Cox proportional hazards models were used to calculate crude and fully adjusted HRs (detemir/glargine), with age as the underlying timescale and the calendar year of the index date as a categorical variable. The (fully) adjusted models included a list of pre-specified covariates and relevant additional covariates that influenced the HR by > 10% on a log scale. The sub-analyses were undertaken in subsets of people with cardiovascular risk factors. BMI body mass index, CI confidence interval, detemir insulin detemir, glargine insulin glargine, HbA1c glycated haemoglobin, HR hazard ratio

Sub-analyses

Across all the subsets of people with cardiovascular risk factors, detemir was associated with a significantly lower risk of all-cause mortality when compared with glargine (Fig. 3a). Results of the sub-analyses for cardiovascular mortality tended to be consistent with the main analyses: a non-significant numerical trend towards a lower risk in people who initiated detemir than in those who initiated glargine (Fig. 3b). These associations were most pronounced in people with a BMI ≥ 30 kg/m2; for all-cause mortality, the crude and fully adjusted HRs (95% CI) were 0.79 (0.69; 0.91) and 0.80 (0.70; 0.92), respectively. For cardiovascular mortality, the crude and fully adjusted HRs (95% CI) were 0.71 (0.52; 0.99) and 0.69 (0.50; 0.96), respectively.

Sensitivity Analyses

Results of the sensitivity analyses are presented in Supplementary Fig. S2. The point estimates for the HRs comparing detemir with glargine for all-cause or cardiovascular mortality were similar across all of the sensitivity analyses. For all-cause mortality, the significance of associations was unchanged across all the sensitivity analyses, with the exception of the adjusted ‘current prescription’ sensitivity analysis 4 (Supplementary Fig. S2).

Discussion

In this UK-based real-world study, there was an association between all-cause mortality and basal insulin choice in insulin-naïve people with T2D; the mortality risk was significantly lower (by 14%) with detemir versus glargine after adjustment for potential confounders. Additionally, there was a non-significant trend towards a 17% lower risk of cardiovascular mortality associated with detemir versus glargine, after adjustment for potential confounders. In general, these findings were consistent in subsets of people with cardiovascular risk factors and in sensitivity analyses that varied methodological assumptions. In sub-analyses, the associations between basal insulin choice and mortality risk were most pronounced in people with obesity; there were significant adjusted risk reductions of 20% and 31% for all-cause and cardiovascular mortality, respectively, with detemir compared with glargine.

A recent population-based study in Finland investigated the differences in all-cause and cause-specific mortality rates between new users of detemir and glargine who had T2D and were aged ≥ 40 years [11]. The results from that study showed the HR (95% CI) for all-cause and cardiovascular mortality for detemir versus glargine to be 0.71 (0.54; 0.93) and 0.64 (0.43; 0.95), respectively [11]. These findings concur with, but are more pronounced than, those in the current study. This observation may be due to the Finnish study including only glargine U100, whereas users of both glargine U100 and U300 formulations were included in the current study. The follow-up time and covariates tested in the current study also differed from those in the Finnish study.

In sub-analyses, the associations between basal insulin choice and mortality risk were most pronounced in people with obesity (BMI ≥ 30 kg/m2). One potential reason underlying associations between lower all-cause or cardiovascular mortality and the initiation of detemir rather than glargine in people with T2D and obesity is the lower weight gain that has been documented with detemir relative to glargine in people with diabetes [8, 9, 18, 19]. Thus, people with diabetes and obesity may be at a comparatively lower risk of subsequent CVD when taking detemir (as compared with glargine). While the present analysis is unable to answer this question directly, the observation that the HR was numerically lowest in people with obesity is of interest and may warrant further investigation. A previous analysis of RCT data found that the magnitude of the difference in weight gain with detemir and other basal insulins increased with patients’ BMI [20]. Thus, it is plausible that this effect might have contributed to the relative cardiovascular mortality benefit.

Another potential contributing factor to the difference in HR is that the two basal insulins differ in their pharmacokinetic/pharmacodynamic profiles. Compared with neutral protoamine Hagedorn (NPH) insulin, which has highly variable absorption, tending to peak 4–5 h after injection, both detemir and glargine have a slower action, less pronounced peaks and a more consistent glucose-lowering effect over time [21, 22].

The present study has limitations, including that we were not able to completely account for confounding by indication. There may have been misclassification due to the assumption that “no information” in the register indicates “no exposure” for some of the variables, and this may not always have been true. The information in primary care records can be incomplete. The CPRD represents a UK real-world clinical practice population with respect to age, sex, BMI and ethnicity,; however, not necessarily with respect to geography and practice size [23]. Detemir may have been chosen (as per guidelines) for people who have concerns about hypoglycaemia or weight gain; therefore, the results may underestimate the benefits. This study included patients treated with either glargine U100 or U300 and it may not be appropriate to combine these products into a single treatment group as they differ at least in terms of their pharmacokinetic profiles [24]. We only had information available on prescriptions and did not know whether people collected the medication from the pharmacy and adhered to it.

This study also has strengths. For instance, the analyses were based on reliable data from a large national register that enabled assessment of a large cohort with a long duration of follow-up. The CPRD Gold database is the largest, most established and most utilised primary care database in the UK [17], and the information stored within the database has been extensively validated (reviewed in [25]). The data were adjusted for several potential confounding factors, and the sensitivity analyses, enabled by the large number of deaths, demonstrated that the HRs for mortality were robust to changes in the analysis methodology. As observational studies face methodological challenges, including confounding by indication, a new-user study design could help to mitigate these challenges by improving the validity of the study findings [26]. Furthermore, our findings are in line with those available from previous studies. Lastly, the total exposure time was relatively high.

Conclusions

Overall, the results from the present real-world study in the UK have shown that there was an association between all-cause mortality and basal insulin choice in insulin-naïve people with T2D; there was a significantly lower mortality risk with detemir versus glargine after adjustment for potential confounders. Additionally, there was a non-significant trend towards a lower risk of cardiovascular mortality with detemir versus glargine, again after adjustment for potential confounders. Of potential interest, the lower mortality risk associated with detemir was most pronounced in a subset of people with obesity.