Comparative effectiveness of dulaglutide versus liraglutide in Asian type 2 diabetes patients: a multi-institutional cohort study and meta-analysis
Cardiovascular Diabetology volume 19, Article number: 172 (2020)
Head-to-head comparison of clinical effectiveness between dulaglutide and liraglutide in Asia is limited. This study was aimed to assess the real-world comparative effectiveness of dulaglutide versus liraglutide.
We conducted a retrospective cohort study by utilizing multi-institutional electronic medical records to identify real-world type 2 diabetes patients treated with dulaglutide or liraglutide during 2016–2018 in Taiwan and followed up until 2019. Effectiveness outcomes were assessed at every 3 months in the 1-year follow-up. Propensity score techniques were applied to enhance between-group comparability. Significant differences in changes of effectiveness outcomes between treatment groups during the follow-up were examined and further analyzed using mixed-model repeated-measures approaches.
A total of 1512 subjects receiving dulaglutide and 1513 subjects receiving liraglutide were identified. At 12 months, significant HbA1c changes from baseline were found in both treatments (dulaglutide: − 1.06%, p < 0.001; liraglutide: − 0.83%, p < 0.001), with a significant between-group difference (− 0.23%, 95% confidence interval − 0.38 to − 0.08%, p < 0.01). Both treatments yielded significant declines in weight, alanine aminotransferase level, and estimated glomerular filtration rate from baseline (dulaglutide: − 1.14 kg, − 3.08 U/L and − 2.08 mL/min/1.73 m2, p < 0.01; liraglutide: − 1.64 kg, − 3.65 U/L and − 2.33 mL/min/1.73 m2, p < 0.001), whereas only dulaglutide yielded a significant systolic blood pressure reduction (− 2.47 mmHg, p < 0.001). Between-group differences in changes of weight, blood pressure, and liver and renal functions at 12 months were not statistically significant.
In real-world T2D patients, dulaglutide versus liraglutide was associated with better glycemic control and comparable effects on changes of weight, blood pressure, and liver and renal functions.
Treatment with glucagon-like peptide-1 receptor agonists (GLP-1RAs) is recommended as a second- or third-line therapeutic option for patients with inadequately controlled type 2 diabetes (T2D) . GLP-1RAs have exhibited promising glycemic efficacy, with additional favorable effects of weight loss, blood pressure reduction, and renal function preservation, and pose a low risk of hypoglycemia, all of which contribute to the desirable cardiovascular outcomes [2, 3].
Long-acting agents (liraglutide, dulaglutide, semaglutide) generally offer better glycemic control than a short-acting agent (lixisenatide, exenatide) . Liraglutide and dulaglutide are the two commonly-used GLP-1RAs in East Asia. There are two head-to-head randomized, phase III trials of dulaglutide versus liraglutide. The international AWARD-6 trial  for patients recruited from nine countries and the other trial for Japanese patients  both demonstrated non-inferior HbA1c reduction between dulaglutide and liraglutide users at 26 weeks (between-group difference: − 0.06% in the AWARD-6 trial and − 0.10% in the Japanese trial, pnon-inferiority < 0.001). The Japanese trial further showed more favorable HbA1c reduction with dulaglutide versus liraglutide at 52 weeks (1.39% versus 1.19%; between-group difference: − 0.20%, p = 0.04) . However, these trials had a relatively short follow-up period (26 weeks [5, 6]), small sample size (< 300 patients per treatment arm [5,6,7]), or a lack of generalizability due to the highly-selected populations [5,6,7].
Real-world evidence obtained from comparative effectiveness research can be used to complement evidence from trials through translating the efficacy of interventions in trials to the effectiveness of them in clinical practice among a broader spectrum of patient populations . Evidence regarding the head-to-head comparative effectiveness of GLP-1RAs is critical for supporting clinical decisions and formulating healthcare reimbursement policies in real-world practice. There are few real-world studies, from Spain , Italy , Canada , United States [12, 13], India , and Scandinavian countries  with a limited number of study subjects (25–585 per treatment arm [9, 11,12,13,14]), patients with a specific comorbidity (i.e., solid organ transplant ), a short follow-up period (e.g., 13 weeks ), few clinical effectiveness measures (i.e., HbA1c only  or HbA1c, weight, and blood pressure only [9,10,11, 13, 14]), or specific severe clinical outcome events of interest (i.e., renal replacement therapy, death from renal causes, and hospitalization for renal events ). To date, there are no published real-word comparative effectiveness studies on GLP-1RAs for Western Pacific patients with T2D. We utilized multi-institutional electronic medical records (EMRs) to identify real-world Taiwanese T2D patients receiving GLP-1RAs for the comparative effectiveness of dulaglutide versus liraglutide on glycemic control, weight, blood pressure, and liver and renal functions.
The Chang Gung Research Database (CGRD) was utilized. It comprises de-identified individual EMRs of disease diagnoses, medical visits (outpatient, inpatient, and emergency room), pharmacy records, examination reports, and laboratory data from seven medical institutes throughout Taiwan, covering 1.3 million individuals (about 6% of Taiwan’s total population) . Its validity for real-world pharmacoepidemiological studies is documented elsewhere [16,17,18,19,20].
As shown in Fig. 1, we identified patients diagnosed with T2D (International Classification of Diseases, Ninth Revision [ICD-9] disease code of 250.X0 or 250.X2; ICD-10 disease code of E11) during 2015–2019. Then, we included patients aged 18 years or older and newly initiated on dulaglutide or liraglutide in 2016–2018. The first prescription date of dulaglutide or liraglutide was the index date. All patients were followed up from the index date until the end of 2019. The study subject identification period of 2016–2018 for GLP-1RAs users allowed for at least 1-year follow-up period. Patients were excluded if they were diagnosed with type 1 diabetes or gestation diabetes in the year before the index date, or treated with exenatide during the study period. To obtain sufficient data for patient baseline conditions, all study patients were required to have at least one clinic visit and one HbA1c record in the year before the index date.
Study effectiveness outcomes
The primary outcome was the comparison of dulaglutide versus liraglutide on the HbA1c change from the index date (baseline) to 3, 6, 9, and 12 consecutive months in the 1-year follow-up. We identified HbA1c records within each 3-month interval and the HbA1c value closest to each corresponding assessment time point was used in the analyses. Body weight, systolic blood pressure (SBP), and liver (alanine aminotransferase [ALT]) and renal (estimated glomerular filtration rate [eGFR]) functions were also measured from baseline and every 3 months in the follow-up. We implemented multiple imputations using the Markov chain Monte Carlo method with an expectation maximization algorithm and combined 10 simulations to deal with missing data in the follow-up .
Analytic procedures are outlined in Fig. 1. To minimize potential selection bias and enhance the comparability of study subjects between two treatment groups, we applied a propensity score (PS) matching procedure . Propensity scores of study subjects were estimated using a multivariable logistic regression model based on various patient baseline characteristics listed in Table 1. We used the nearest-neighbor 1:1 PS matching with a caliper of 0.05 on the PS scale with 8-digit greedy matching .
Primary analyses were based on an intention-to-treat (ITT) scenario where the loss to follow-up in the CGRD, death, or end of the 12-month follow-up, whichever came first, was censored. Analyses were divided into two parts. First, we used the paired t-test to estimate changes in clinical effectiveness at 12 months from baseline within each treatment group for assessing the within-group difference, and then used the two-sample t-test to determine the between-group difference in changes of clinical effectiveness at 12 months from baseline. Second, to consider time-varying changes in biomarkers (e.g., HbA1c) that were repeatedly assessed every 3 months during the follow-up, we performed a mixed-model analysis to consider treatment groups, assessment time points, and the interaction of treatment groups and assessment time points as fixed effects and individual patients as a random effect .
A series of sensitivity and subgroup analyses were conducted. First, to account for possible over-estimation of treatment effects in the ITT analyses where non-adherence to treatments was ignored, we performed the as-treated analysis where patients who switched away from or discontinued the use of a study drug were also censored, in addition to the censoring defined in the ITT analyses (Sensitivity 1). Second, to avoid potential confounding from short-term or accidental use of GLP-1RAs, we performed analyses where only stable users were included (Sensitivity 2). Stable users were defined as patients who had at least three consecutive refills of dulaglutide or liraglutide with any gaps between two consecutive refills of less than 90 days . Third, we performed analyses with adjustment for potential healthy user bias (Sensitivity 3). Specifically, the patients who received GLP-1RAs and also used a dipeptidyl peptidase 4 inhibitor (DPP-4i) or sodium glucose cotransporter 2 inhibitor (SGLT-2i) were identified as possible healthy users because the combined use of a GLP-1RAs with a DPP-4i or SGLT-2i is not reimbursed by Taiwan’s National Health Insurance program and patients have to pay out-of-pocket fees. Under this circumstance, patients who are willing to self-pay for more intensive treatments would be likely to be engaged in healthier behaviors. We re-ran the analyses using a subset of patients who did not use a DPP-4i or SGLT-2i in combination with GLP-1RAs to avoid potential healthy user bias. Fourth, to enhance the study generalizability through retaining study cohort patients as many as possible, we applied two PS weighting procedures, inverse probability of treatment weighting (IPTW) and standardized mortality ratio weighting (SMRW)  (Sensitivities 4 and 5). Specifically, the patients at the 5th to 95th percentiles of the distribution of PS were first trimmed to minimize potential residual confounding . Then, for IPTW, dulaglutide users were weighted as the inverse of the estimated PS and liraglutide users were weighted as the inverse of 1 minus the estimated PS. For SMRW, dulaglutide users were given a weight of 1 and liraglutide users were given a weight based on the ratio of the estimated PS to 1 minus the estimated PS.
In subgroup analyses, the procedures that were performed in the primary analyses were applied to examine the treatment effects on study outcomes in subgroups according to a series of patient baseline characteristics, including HbA1c (≥ 9%, < 9%), age (≥ 65 years, < 65 years), eGFR (≥ 60 mL/min/1.73 m2, < 60 mL/min/1.73 m2), ALT (> upper normal limit [UNL], ≤ UNL), and body mass index (≥ 27 kg/m2, < 27 kg/m2). A two-tail p-value of less than 0.05 was considered statistically significant. Data were analyzed using SAS Enterprise Guide, version 7.1 (SAS Institute, Cary, NC, USA).
We further performed a meta-analysis on clinical effectiveness (i.e., HbA1c, weight and SBP) of liraglutide vs. dulaglutide by pooling the results from prior studies and the present study. Two reviewers (Chang and Shao) independently searched studies from the PubMed and Embase from the inception of database to May 31, 2020 that reported the comparison of liraglutide and dulaglutide. The search strategy and key terms were listed in Additional file 1: Appendix Table S1. Effectiveness outcomes abovementioned were measured from 6 and 12 months of follow-up periods. We included both randomized control trials (RCTs) and observational studies without imposing any language restrictions. Data were presented as mean difference with 95% CIs. We conducted the random-effects model meta-analysis using the reverse invariance method. The statistical heterogeneity was assessed by the statistic I2. To minimize the heterogeneity of included studies, we further conducted subgroups analyses for meta-analysis of RCTs or observational studies only. Data were analyzed by Review Manager version 5.3 (Copenhagen: The Nordic Cochrane Centre, The Cochrane Collaboration, 2014).
Baseline characteristics of study patients
1512 and 1513 patients newly-initiated on dulaglutide and liraglutide, respectively, were included for the analyses under PS weighting procedures (IPTW and SMRW; Fig. 1). After PS matching, 983 matched pairs of dulaglutide and liraglutide new users were identified. Additional file 1: Appendix Figure S1 illustrates the kernel density for the PS distributions of the treatment groups. Table 1 shows patients’ characteristics for dulaglutide and liraglutide users before and after the PS matching. All patients’ characteristics were comparable between treatment groups after the PS matching.
Based on the PS-matched sample, there was a statistically significant change in HbA1c at 12 months from baseline in each treatment group (dulaglutide: − 1.06% [standard deviation (SD): 1.70] and liraglutide: − 0.83% [SD: 1.61]), with a significant between-group difference in HbA1c reduction of − 0.23% (95% confidence interval [CI] − 0.38 to − 0.08%) (Table 2). Figure 2a, b show a large HbA1c reduction trend associated with dulaglutide versus liraglutide over the 1-year follow-up.
Body weight loss, blood pressure control, and liver and renal functions
At 12 months, the changes in body weight, SBP, eGFR, and ALT from baseline in dulaglutide users were − 1.14 kg, − 2.47 mmHg, − 2.08 mL/min/1.73 m2, and − 3.08 U/L, respectively, and those in liraglutide users were − 1.64 kg, − 0.56 mmHg, − 2.33 mL/min/1.73 m2, and − 3.65 U/L, respectively (Table 2 and Fig. 2). Among these outcomes, both dulaglutide and liraglutide users had significant reduction in body weight, eGFR, and ALT at 12 months from baseline, while only dulaglutide users had significant SBP reduction. The between-group differences in the changes in body weight, SBP, eGFR, and ALT levels did not reach statistical significance. Similar results were observed in the mixed-model analyses (Table 2).
Sensitivity and subgroup analyses
A series of sensitivity analyses show consistent results (Additional file 1: Appendix Table S2) with the primary analyses (Table 2). The results of subgroup analyses are summarized in Additional file 1: Appendix Table S3 and Figure S2. Generally, there was a consistent benefit of dulaglutide versus liraglutide on HbA1c and SBP across all subgroups. In contrast, there was some heterogeneous treatment effects for other outcomes across subgroups. For example, at 12 months, liraglutide use yielded greater body weight reduction compared to dulaglutide use among patients with age ≥ 65 years old (between-group difference in weight reduction: 1.30 kg, p < 0.05). Among patients with abnormal liver function (ALT > UNL), the ALT levels in both treatment groups significantly declined at 12 months (i.e., dulaglutide: − 19.13 U/L [SD: 36.31] and liraglutide: − 22.79 U/L [41.85], p < 0.001), despite no statistically significant between-group difference in the ALT change (3.65 U/L, 95% CI − 2.45 to 9.76). The baseline HbA1c level appears to be a significant modifier for the comparative effectiveness of dulaglutide versus liraglutide; the interaction terms of treatment group (dulaglutide versus liraglutide) and baseline HbA1c level (≥ 9% versus < 9%) across different study outcomes were statistically significant (Additional file 1: Appendix Figure S2).
Including the present study, we considered seven studies with up to a total of 5817 T2D patients for meta-analysis. Data extraction flow was detailed in Additional file 1: Appendix Figure S3. We excluded one study  with less than 6 months of follow-up and three studies [9, 13, 15] without sufficient outcome data. The characteristics of the included studies are provided in Additional file 1: Appendix Table S4. The use of dulaglutide versus liraglutide resulted in a significant HbA1c reduction; the pooled mean HbA1c difference between dulaglutide and liraglutide was − 0.17% (95% CI − 0.08 to − 0.26%) and − 0.22% (− 0.11 to − 0.32%) over the 6- and 12-month follow-ups, respectively (Additional file 1: Appendix Figure S4 and Fig. 3). However, there were comparable effects on weight loss and SBP reduction between two drugs (weight loss in Additional file 1: Appendix Figures S5 and S6, SBP reduction in Additional file 1: Appendix Figure S7 and S8). A considerable heterogeneity across the studies was noted for analyzing clinical effectiveness of treatments over a 6-month follow-up (I2 = 52% to 86%).
This large, real-world comparative effectiveness study of GLP-1RAs agents in a Taiwanese population with T2D across multiple medical institutions comprehensively evaluated clinical effectiveness of dulaglutide versus liraglutide. We found that dulaglutide versus liraglutide was associated with a greater HbA1c reduction at 12 months. The benefits of reducing body weight, blood pressure, and ALT levels were also found with these GLP-1RAs treatments, although between-group differences in these beneficial effects were not statistically significant. Moreover, across all pre-specified subgroups, there were consistent beneficial effects of dulaglutide versus liraglutide on HbA1c and SBP, some of which were of statistical significance. Treatment effects of dulaglutide versus liraglutide on clinical outcomes appeared to be modified by baseline HbA1c levels.
The AWARD-6 and Japanese trials first showed the non-inferiority of dulaglutide versus liraglutide at 26 weeks with a mean between-group difference in HbA1c change of − 0.06 to − 0.10% [5, 6]. The post hoc analysis of AWARD-6 showed the similar impact of dulaglutide and liraglutide on relative contribution of basal and postprandial hyperglycemia across HbA1c quartiles after 6 months of treatment . With a longer follow-up period, the Japanese trial reported significant HbA1c reduction for dulaglutide versus liraglutide at 52 weeks with a mean between-group difference in HbA1c change of − 0.20% . However, the small number of study participants in the Japanese trial [6, 7] and the selective and homogenous study populations in well-controlled trial settings [5,6,7] may limit the generalizability of study results to real-world diverse populations treated with GLP-1RAs. Our systematic review found that three real-world studies [10,11,12] have been conducted on the comparative effectiveness of dulaglutide versus liraglutide over a follow-up of 6 or 12 months in non-ethnically Chinese populations (Additional file 1: Appendix Table S4), showing a greater glycemic reduction with dulaglutide versus liraglutide [10, 12], except the Canadian study . Compared to these previous studies with certain limitations such as having a more homogeneous population [5,6,7] or a limited number of patients (i.e., 417–1344 study subjects [5,6,7, 10,11,12]), this present large real-world study of GLP-1RAs in Taiwan reveals that the greater benefit of glycemic control was associated with the use of dulaglutide versus liraglutide over a follow-up of 3–12 months (Fig. 2). It is worth noting that the absolute glycemic reduction with dulaglutide and liraglutide found in this study is slightly smaller than that shown in previous trials. More frequent contact and monitoring, more intensive diabetes education, and better medication adherence in the well-controlled trial setting may contribute this discrepancy. Furthermore, patient populations differ among studies. For example, more complex comorbidities in our study population than previous trials’ populations [5,6,7] might have mitigated the glycemic effect of GLP-1RAs therapy. Nevertheless, consistent with previous studies [7, 10, 12], we found a significantly larger HbA1c reduction associated with dulaglutide versus liraglutide. Our meta-analysis that incorporated all existing data (including the present study) further demonstrated that the use of dulaglutide versus liraglutide possessed a superior glycemic control over a 6-month or 12-month follow-up (Fig. 3 and Additional file 1: Appendix Figure S4). The favorable glycemic control of dulaglutide over liraglutide might be explained by better adherence owing to its once-weekly regimen  and its long-lasting drug action attributable to the longer half-life .
Body weight loss
Weight loss with individual GLP-1RAs drugs has been confirmed in T2D patients . The AWARD-6 trial showed significantly larger body weight reduction associated with liraglutide versus dulaglutide (− 0.71 kg) at 26 weeks , but the Japanese trial showed no difference in weight loss between dulaglutide and liraglutide at 26 or 52 weeks [6, 7]. This implies that the weight loss effect of GLP-1RAs therapy may vary with patients’ baseline body weight and could be culture- or population-specific. Specifically, there was a significant weight reduction with GLP-1RAs use in the AWARD-6 trial participants, which predominantly comprised Caucasians with an average baseline weight of 94.1 kg , whereas there was no significant weight loss with GLP-1RAs use in the Japanese trial, which included Japanese patients with an average baseline body weight of 70.9 kg [6, 7]. Moreover, in Japanese populations, it has been reported that females with the treatment of dulaglutide or liraglutide generally had greater weight loss than males . In our study of the ethnically Chinese population with an average baseline weight of 77.7 kg, the weight loss with dulaglutide (− 1.14 kg) or liraglutide (− 1.64 kg) was statistically significant. Although the weight loss of liraglutide compared to that of dulaglutide was greater, the between-treatment difference was not statistically significant. Another explanation for the weaker weight loss effect of GLP-1RAs in this study compared to that in the AWARD-6 trial  is the concomitant use of glucose-lowering agents (GLAs). Our study subjects were also treated with other GLAs which may have a weight gain effect (e.g., around 60% of patients on sulfonylurea which may have weight gain effect ), whereas the trial patients were treated with only metformin in addition to GLP-1RAs. The concomitant use of other GLAs with a potential weight gain effect may mask the weight benefit of GLP-1RAs. Nevertheless, the favorable weight benefit of liraglutide versus dulaglutide may be explained by the molecule size of the drug and the associated mechanism of weight loss. Because the dulaglutide molecule is larger than the liraglutide molecule, less dulaglutide is transported across the blood–brain barrier or through fenestrated capillaries and thus less effects for increasing satiety and inducing nausea in the central nervous system [29, 33]. As a result, dulaglutide may be less effective for weight reduction compared to liraglutide. Based on the meta-analysis for existing studies (including our study) (Additional file 1: Appendix Figures S5 and S6), the effect of weight loss for liraglutide versus dulaglutide was greater over a 6-month follow-up but comparable over a 12-month follow-up.
Blood pressure control
Reduced blood pressure with GLP-1RAs treatment has been documented in recent cardiovascular outcome trials [34, 35] and linked to better cardiovascular outcomes in GLP-1RAs users . In previous head-to-head comparison trials of GLP-1RAs [5,6,7], blood pressure reduction was found in both dulaglutide and liraglutide groups, while the between-group difference in blood pressure change was not statistically significant. In the present study, only the use of dulaglutide was associated with a significant blood pressure reduction, which is consistent with a recent real-world study from Italy . This discrepancy between trials [5,6,7] and real-world studies (our study and the Italian study ) may be explained by more complicated comorbidities and more use of antihypertension agents in the real-world patient populations, which might mask the blood pressure reduction effects of GLP-1RAs therapy. Limited evidence for blood pressure control between dulaglutide and liraglutide (Additional file 1: Appendix Figures S7 and S8) suggests a need for future research.
Previous studies have supported a reduced ALT level associated with liraglutide use compared to either placebo  or no GLAs use . Consistent with these findings, the present study found that both dulaglutide and liraglutide were associated with a significant ALT decline at 12 months, despite no significant difference in the ALT decline between two treatments. In addition, we found that the magnitude of ALT reduction with GLP-1RAs use was greater in the subgroup of patients with abnormal liver function, which was indicated with a significant interaction between treatment groups (dulaglutide versus liraglutide) and ALT levels (> UNL versus ≤ UNL) (Additional file 1: Appendix Figure S2d). Future research is warranted to corroborate this finding to determine whether the liver outcomes associated with GLP-1RAs use vary with individual GLP-1RAs drugs or patients’ underlying liver function.
Although previous clinical studies have revealed either no change in eGFR (i.e., liraglutide versus placebo ) or reduced eGFR declines (eGFR preservation) (i.e., dulaglutide versus insulin or placebo [35, 40] or liraglutide versus placebo ), there has been no head-to-head comparison of dulaglutide versus and liraglutide on renal outcomes. The present study found that a slight decrease in eGFR was associated with both GLP-1RAs treatments. However, our subgroup analyses further showed that among patients with chronic renal impairment (eGFR < 60 mL/min/1.73 m2), a slight increase in eGFR was associated with both GLP-1RAs treatments, and that the interaction between treatment groups (dulaglutide versus liraglutide) and eGFR levels (eGFR ≥ 60 versus < 60 mL/min/1.73 m2) was statistically significant (Additional file 1: Appendix Figure S2e). This study adds supporting evidence of favorable renal effects with GLP-1RAs therapy in a real-world T2D population that had a relatively poor renal function. Future research should be conducted to determine whether the renal benefit of different GLP-1RAs varies by patients’ underlying renal function.
The present study has several strengths compared to previous studies. First, our study cohort was derived from large EMRs across multiple institutions at different levels of hospitals throughout Taiwan that comprised a diverse real-world T2D population along with their individual-level detailed laboratory measurements, which are typically lacking in administrative datasets, to enrich our analyses. Second, this is the largest real-world Asian study on the comparative effectiveness of GLP-1RAs drugs among a T2D population. Third, we considered a wide range of biochemical markers in the analyses. We are thus able to either corroborate the findings from existing studies or provide additional evidence of clinical effects of GLP-1RAs. Lastly, our rigorous analytic procedures with a series of sensitivity and subgroup analyses that varied with different clinical scenarios ensure the robustness of our study findings and their generalizability to diverse real-world T2D populations treated with GLP-1RAs.
This study also has several limitations. First, as a retrospective design study, possible confounding by indication and unmeasured confounding might not have been avoided. However, we implemented the rigorous PS procedures (e.g., matching, weighting) to enhance between-group comparability and the sensitivity analyses to account for patients’ medication use behaviors (i.e., as-treated scenario, stable GLP-1RAs users) and possible healthy user bias to minimize potential confounding and bias commonly seen in retrospective studies. Second, the present study only assessed the clinical effectiveness of GLP-1RAs in terms of short-term clinical biochemical marker changes (e.g., HbA1c, eGFR), while the hard endpoints of treatments (e.g., cardiovascular disease, death, progression to end-stage renal disease) were not measured. To date, there is a lack of direct head-to-head comparative trials of GLP-1RAs on the long-term cardiovascular safety and mortality, but there are indirect comparisons from three network meta-analysis studies [42,43,44]. This suggests that future studies with the long-term follow-up period on hard outcomes among GLP-1RAs are needed. Third, comparison of dulaglutide with semaglutide would be also important as they are both once-weekly injections. However, semaglutide was not available in Taiwan during our study period and thus was not analyzed in this study. Future research for comparative effectiveness of dulaglutide versus semaglutide is needed. Lastly, it is possible that our study patients had visited other medical institutions outside the CGRD system. The lack of continuity in medical visits would affect the completeness of our patient records. However, because our primary analysis was based on study patients with at least one HbA1c value in the year of the follow-up and our sensitivity analysis was also performed by assessing stable GLP-1RAs users, our study patients are most likely to be loyal patients in the CGRD system. We further analyzed the rates of loss to follow-up between treatment groups and found similar rates between two treatment groups (dulaglutide versus liraglutide: 12.9% versus 14.3%). This implies that the potential bias attributable to incomplete follow-up records could be negligible.
In this large real-world T2D population treated with GLP-1RAs, dulaglutide was associated with more favorable glycemic control compared to liraglutide. Weight loss, blood pressure reduction, and improved liver and renal functions after GLP-1RAs treatments were also found. The differences in these clinical outcomes between dulaglutide and liraglutide were comparable. Future research on the comparative effectiveness of among other GLP-1RAs drugs and a comparison of GLP-1RAs with other GLAs is warranted to support the clinical rationale of selecting GLP-1RAs treatment.
Availability of data and materials
Data sharing is not applicable to this study as data management and analysis were performed on a statistics server through remote access in Chang Gung Medical Foundation in Taiwan, for privacy and safety concerns.
Glucagon-like peptide-1 receptor agonists
Type 2 diabetes
Electronic medical records
Chang Gung Research Database
International Classification of Diseases, Ninth Revision
International Classification of Diseases, Tenth Revision
Systolic blood pressure
Estimated glomerular filtration rate
Dipeptidyl peptidase 4 inhibitor
Sodium glucose cotransporter 2 inhibitor
Inverse probability of treatment weighting
Standardized mortality ratio weighting
American Diabetes Association. 9. Pharmacologic approaches to glycemic treatment: standards of medical care in diabetes-2020. Diabetes Care. 2020;43(Suppl 1):S98–s110.
Kristensen SL, Rorth R, Jhund PS, Docherty KF, Sattar N, Preiss D, Kober L, Petrie MC, McMurray JJV. Cardiovascular, mortality, and kidney outcomes with GLP-1 receptor agonists in patients with type 2 diabetes: a systematic review and meta-analysis of cardiovascular outcome trials. Lancet Diabetes Endocrinol. 2019;7(10):776–85.
Sposito AC, Berwanger O, de Carvalho LSF, Saraiva JFK. GLP-1RAs in type 2 diabetes: mechanisms that underlie cardiovascular effects and overview of cardiovascular outcome data. Cardiovasc Diabetol. 2018;17(1):157.
Pinelli NR, Hurren KM. Efficacy and safety of long-acting glucagon-like peptide-1 receptor agonists compared with exenatide twice daily and sitagliptin in type 2 diabetes mellitus: a systematic review and meta-analysis. Ann Pharmacother. 2011;45(7–8):850–60.
Dungan KM, Povedano ST, Forst T, González JGG, Atisso C, Sealls W, Fahrbach JL. Once-weekly dulaglutide versus once-daily liraglutide in metformin-treated patients with type 2 diabetes (AWARD-6): a randomised, open-label, phase 3, non-inferiority trial. Lancet. 2014;384(9951):1349–57.
Miyagawa J, Odawara M, Takamura T, Iwamoto N, Takita Y, Imaoka T. Once-weekly glucagon-like peptide-1 receptor agonist dulaglutide is non-inferior to once-daily liraglutide and superior to placebo in Japanese patients with type 2 diabetes: a 26-week randomized phase III study. Diabetes Obes Metab. 2015;17(10):974–83.
Odawara M, Miyagawa J, Iwamoto N, Takita Y, Imaoka T, Takamura T. Once-weekly glucagon-like peptide-1 receptor agonist dulaglutide significantly decreases glycated haemoglobin compared with once-daily liraglutide in Japanese patients with type 2 diabetes: 52 weeks of treatment in a randomized phase III study. Diabetes Obes Metab. 2016;18(3):249–57.
Seeger JD, Nunes A, Loughlin AM. Using RWE research to extend clinical trials in diabetes: an example with implications for the future. Diabetes Obes Metab. 2020;22(Suppl 3):35–44.
Tofé S, Argüelles I, Mena E, Serra G, Codina M, Urgeles JR, García H, Pereg V. Real-world GLP-1 RA therapy in type 2 diabetes: a long-term effectiveness observational study. Endocrinol Diabetes Metabol. 2019;2(1):e00051.
Morieri ML, Rigato M, Frison V, Simioni N, D'Ambrosio M, Tadiotto F, Paccagnella A, Lapolla A, Avogaro A, Fadini GP. Effectiveness of dulaglutide vs liraglutide and exenatide once-weekly. A real-world study and meta-analysis of observational studies. Metabolism. 2020;106:154190.
Brown RE, Abitbol A, Bajaj HS, Goldenberg R, Khandwala H, Abdel-Salam S, Aronson R. Patient reported outcomes following initiation of glucagon-like peptide-1 receptor agonists in patients with type 2 diabetes in a specialist endocrinology practice of the LMC diabetes registry: the PROGRESS-diabetes study. Diabetes Res Clin Pract. 2019;156:107820.
Mody R, Huang Q, Yu M, Zhao R, Patel H, Grabner M, Landó LF. Adherence, persistence, glycaemic control and costs among patients with type 2 diabetes initiating dulaglutide compared with liraglutide or exenatide once weekly at 12-month follow-up in a real-world setting in the United States. Diabetes Obes Metab. 2019;21(4):920–9.
Singh P, Taufeeq M, Pesavento TE, Washburn K, Walsh D, Meng S. Comparison of the glucagon-like-peptide-1 receptor agonists dulaglutide and liraglutide for the management of diabetes in solid organ transplant: a retrospective study. Diabetes Obes Metab. 2020;22(5):879–84.
Ghosal S, Sinha B. Liraglutide and dulaglutide therapy in addition to SGLT-2 inhibitor and metformin treatment in Indian type 2 diabetics: a real world retrospective observational study. Clin Diabetes Endocrinol. 2018;4:11.
Pasternak B, Wintzell V, Eliasson B, Svensson AM, Franzén S, Gudbjörnsdottir S, Hveem K, Jonasson C, Melbye M, Svanström H, et al. Use of glucagon-like peptide 1 receptor agonists and risk of serious renal events: scandinavian cohort study. Diabetes Care. 2020;43(6):1326–35.
Shao SC, Chan YY, Kao Yang YH, Lin SJ, Hung MJ, Chien RN, Lai CC, Lai EC. The Chang Gung Research Database-a multi-institutional electronic medical records database for real-world epidemiological studies in Taiwan. Pharmacoepidemiol Drug Saf. 2019;28(5):593–600.
Shao SC, Chang KC, Chien RN, Lin SJ, Hung MJ, Chan YY, Kao Yang YH, Lai EC. Effects of sodium-glucose co-transporter-2 inhibitors on serum alanine aminotransferase levels in people with type 2 diabetes: a multi-institutional cohort study. Diabetes Obes Metab. 2020;22(1):128–34.
Shao SC, Chang KC, Hung MJ, Yang NI, Chan YY, Chen HY, Kao Yang YH, Lai EC. Comparative risk evaluation for cardiovascular events associated with dapagliflozin vs. empagliflozin in real-world type 2 diabetes patients: a multi-institutional cohort study. Cardiovasc Diabetol. 2019;18(1):120.
Shao SC, Lin YH, Chang KC, Chan YY, Hung MJ, Kao Yang YH, Lai EC. Sodium glucose co-transporter 2 inhibitors and cardiovascular event protections: how applicable are clinical trials and observational studies to real-world patients? BMJ Open Diabetes Res Care. 2019;7(1):e000742.
Shao SC, Chang KC, Lin SJ, Chien RN, Hung MJ, Chan YY, Kao Yang YH, Lai EC. Favorable pleiotropic effects of sodium glucose cotransporter 2 inhibitors: head-to-head comparisons with dipeptidyl peptidase-4 inhibitors in type 2 diabetes patients. Cardiovasc Diabetol. 2020;19(1):17.
Harel O, Zhou XH. Multiple imputation: review of theory, implementation and software. Stat Med. 2007;26(16):3057–77.
Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivar Behav Res. 2011;46(3):399–424.
Rassen JA, Shelat AA, Myers J, Glynn RJ, Rothman KJ, Schneeweiss S. One-to-many propensity score matching in cohort studies. Pharmacoepidemiol Drug Saf. 2012;21(Suppl 2):69–80.
Detry MA, Ma Y. Analyzing repeated measurements using mixed models. JAMA. 2016;315(4):407–8.
Hou WH, Chang KC, Li CY, Ou HT. Dipeptidyl peptidase-4 inhibitor use is associated with decreased risk of fracture in patients with type 2 diabetes: a population-based cohort study. Br J Clin Pharmacol. 2018;84(9):2029–39.
Desai RJ, Franklin JM. Alternative approaches for confounding adjustment in observational studies using weighting based on the propensity score: a primer for practitioners. BMJ (Clin Res Ed). 2019;367:l5657.
Brookhart MA, Wyss R, Layton JB, Stürmer T. Propensity score methods for confounding control in nonexperimental research. Circ Cardiovasc Qual Outcomes. 2013;6(5):604–11.
Umpierrez G, Pantalone KM, Atisso CM, Landó LF, Patel H. Relative contribution of basal and postprandial hyperglycaemia stratified by HbA1c categories before and after treatment intensification with dulaglutide. Diabetes Obes Metab. 2019;21(6):1365–72.
Lund A, Knop FK, Vilsbøll T. Glucagon-like peptide-1 receptor agonists for the treatment of type 2 diabetes: Differences and similarities. Eur J Internal Med. 2014;25(5):407–14.
Vilsboll T, Christensen M, Junker AE, Knop FK, Gluud LL. Effects of glucagon-like peptide-1 receptor agonists on weight loss: systematic review and meta-analyses of randomised controlled trials. BMJ (Clin Res Ed). 2012;344:d7771.
Onishi Y, Oura T, Matsui A, Matsuura J, Iwamoto N. Analysis of efficacy and safety of dulaglutide 0.75 mg stratified by sex in patients with type 2 diabetes in 2 randomized, controlled phase 3 studies in Japan. Endocr J. 2017;64(5):553–60.
Tran L, Zielinski A, Roach AH, Jende JA, Householder AM, Cole EE, Atway SA, Amornyard M, Accursi ML, Shieh SW, et al. Pharmacologic treatment of type 2 diabetes: oral medications. Ann Pharmacother. 2015;49(5):540–56.
Madsbad S. Review of head-to-head comparisons of glucagon-like peptide-1 receptor agonists. Diabetes Obes Metab. 2016;18(4):317–32.
Marso SP, Daniels GH, Brown-Frandsen K, Kristensen P, Mann JF, Nauck MA, Nissen SE, Pocock S, Poulter NR, Ravn LS, et al. Liraglutide and cardiovascular outcomes in type 2 diabetes. N Engl J Med. 2016;375(4):311–22.
Gerstein HC, Colhoun HM, Dagenais GR, Diaz R, Lakshmanan M, Pais P, Probstfield J, Riesmeyer JS, Riddle MC, Ryden L, et al. Dulaglutide and cardiovascular outcomes in type 2 diabetes (REWIND): a double-blind, randomised placebo-controlled trial. Lancet. 2019;394(10193):121–30.
Sowers JR, Epstein M, Frohlich ED. Diabetes, hypertension, and cardiovascular disease. Hypertension. 2001;37(4):1053–9.
Cusi K, Sattar N, García-Pérez L-E, Pavo I, Yu M, Robertson KE, Karanikas CA, Haupt A. Dulaglutide decreases plasma aminotransferases in people with Type 2 diabetes in a pattern consistent with liver fat reduction: a post hoc analysis of the AWARD programme. Diabet Med. 2018;35(10):1434–9.
Bajaj HS, Brown RE, Bhullar L, Sohi N, Kalra S, Aronson R. SGLT2 inhibitors and incretin agents: associations with alanine aminotransferase activity in type 2 diabetes. Diabetes Metab. 2018;44(6):493–9.
Davies MJ, Bain SC, Atkin SL, Rossing P, Scott D, Shamkhalova MS, Bosch-Traberg H, Syren A, Umpierrez GE. Efficacy and safety of liraglutide versus placebo as add-on to glucose-lowering therapy in patients with type 2 diabetes and moderate renal impairment (LIRA-RENAL): a randomized clinical trial. Diabetes Care. 2016;39(2):222–30.
Tuttle KR, Lakshmanan MC, Rayner B, Busch RS, Zimmermann AG, Woodward DB, Botros FT. Dulaglutide versus insulin glargine in patients with type 2 diabetes and moderate-to-severe chronic kidney disease (AWARD-7): a multicentre, open-label, randomised trial. Lancet Diabetes Endocrinol. 2018;6(8):605–17.
Mann JFE, Ørsted DD, Brown-Frandsen K, Marso SP, Poulter NR, Rasmussen S, Tornøe K, Zinman B, Buse JB. Liraglutide and renal outcomes in type 2 diabetes. N Engl J Med. 2017;377(9):839–48.
Alfayez OM, Almohammed OA, Alkhezi OS, Almutairi AR, Al Yami MS. Indirect comparison of glucagon like peptide-1 receptor agonists regarding cardiovascular safety and mortality in patients with type 2 diabetes mellitus: network meta-analysis. Cardiovasc Diabetol. 2020;19(1):96.
Zheng SL, Roddick AJ, Aghar-Jaffar R, Shun-Shin MJ, Francis D, Oliver N, Meeran K. Association between use of sodium-glucose cotransporter 2 inhibitors, glucagon-like peptide 1 agonists, and dipeptidyl peptidase 4 inhibitors with all-cause mortality in patients with type 2 diabetes: a systematic review and meta-analysis. JAMA. 2018;319(15):1580–91.
Tsapas A, Avgerinos I, Karagiannis T, Malandris K, Manolopoulos A, Andreadis P, Liakos A, Matthews DR, Bekiari E. Comparative effectiveness of glucose-lowering drugs for type 2 diabetes: a systematic review and network meta-analysis. Ann Intern Med. 2020;173(4):278–86.
The authors would like to thank Chang Gung Memorial Hospitals for providing Chang Gung Research Database data. The interpretation and conclusions contained herein do not represent the presentation of Chang Gung Memorial Hospital.
Parts of this study were presented in abstract and poster forms at ISPOR Europe 2019, Copenhagen, Denmark, November 2–6, 2019.
This project was supported by grants from the Ministry of Science and Technology in Taiwan (Grant MOST 108-2320-B-006-048) (Huang-Tz Ou) and from the Chang Gung Medical Foundation (ID: CMRPG3K0101) (Kai-Cheng Chang). Technical/expert advice for this project was supported by the MCDTR from the National Institute of Diabetes and Digestive and Kidney Diseases (Grant number P30DK092926) (Shihchen Kuo). The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Ethics approval and consent to participate
This study was approved by Institutional Review Board of the Chang Gung Medical Foundation (No. 201801493B0C601).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Distribution of kernel density of the propensity score distribution in dulaglutide and liraglutide users before and after propensity score matching. Table S1. Search strategy and key terms for meta-analysis. Table S2. Comparison of clinical effectiveness between liraglutide and dulaglutide (sensitivity analyses). Table S3. Subgroup analyses for comparison of clinical effectiveness changes between dulaglutide and liraglutide at 12 months (based on the propensity-score-matched sample). Figure S2. Changes in clinical effectiveness between dulaglutide and liraglutide at 12 months stratified by patient subgroup (based on the propensity-score-matched sample). Figure S3. Flow chart of selection of studies included in the meta-analysis. Table S4. Summary of existing studies that head-to-head compared dulaglutide and liraglutide. Figure S4. Forest plot of 6-month difference in HbA1c between dulaglutide and liraglutide. Figure S5. Forest plot of 12-month difference in weight loss between dulaglutide and liraglutide. Figure S6. Forest plot of 6-month difference in weight loss between dulaglutide and liraglutide. Figure S7. Forest plot of 12-month difference in systolic blood pressure change between dulaglutide and liraglutide. Figure S8. Forest plot of 6-month difference in systolic blood pressure change between dulaglutide and liraglutide.
About this article
Cite this article
Chang, KC., Shao, SC., Kuo, S. et al. Comparative effectiveness of dulaglutide versus liraglutide in Asian type 2 diabetes patients: a multi-institutional cohort study and meta-analysis. Cardiovasc Diabetol 19, 172 (2020). https://doi.org/10.1186/s12933-020-01148-8