Effect of the Glucagon-Like Peptide-1 Receptor Agonists on Autonomic Function in Subjects with Diabetes: A Systematic Review and Meta-Analysis
Article information
Abstract
Background
In addition to the metabolic effects in diabetes, glucagon-like peptide 1 receptor (GLP-1R) agonists lead to a small but substantial increase in heart rate (HR). However, the GLP-1R actions on the autonomic nervous system (ANS) in diabetes remain debated. Therefore, this meta-analysis evaluates the effect of GLP-1R agonist on measures of ANS function in diabetes.
Methods
According to the Cochrane Collaboration and Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) statement, we conducted a meta-analysis considering clinical trials in which the autonomic function was evaluated in diabetic subjects chronically treated with GLP-1R agonists. The outcomes were the change of ANS function measured by heart rate variability (HRV) and cardiac autonomic reflex tests (CARTs).
Results
In the studies enrolled, HR significantly increased after treatment (P<0.001), whereas low frequency/high frequency ratio did not differ (P=0.410); no changes in other measures of HRV were detected. Considering CARTs, only the 30:15 value derived from lying-to-standing test was significantly lower after treatment (P=0.002), but only two studies reported this measurement. No differences in other CARTs outcome were observed.
Conclusion
The meta-analysis confirms the HR increase but seems to exclude an alteration of the sympatho-vagal balance due to chronic treatment with GLP-1R agonists in diabetes, considering the available measures of ANS function.
INTRODUCTION
The diabetic autonomic neuropathy is defined as a heterogeneous category of disorders of the autonomic nervous system (ANS) in individuals with either diabetes mellitus or metabolic derangements of pre-diabetes, when other potential causes have been excluded [1]. In particular, the cardiac autonomic neuropathy (CAN) is the manifestation of an ANS imbalance, due to the impairment of autonomic control of the cardiovascular system [1,2]. CAN affects at least 20% of unselected patients, and this incidence raises up to 65% of subjects with either increasing age or long diabetes duration [2-4]. However, the actual CAN prevalence varies, depending upon diagnostic criteria, patient cohort, and testing modality [2,5]. According to standard cardiac autonomic reflex tests (CARTs), the CAN prevalence is reported about 7% in type 2 diabetes mellitus (T2DM) and increases with diabetes duration by 4.6% to 6% per year [2,5]. It is well known that CAN incidence is influenced by diabetic disease duration, patient age, glycaemic control, and concomitant metabolic syndrome features [6,7]. The CAN diagnosis in people with diabetes is extremely relevant, influencing the prognosis for cardiovascular morbidity and predicting the overall cardiovascular risk [8]. Indeed, CAN is related to silent myocardial ischemia, stroke, postural hypotension, exercise intolerance and enhanced intraoperative instability [9]. Therefore, from a clinical perspective, an individual with diabetes and CAN is at higher risk of mortality and of cardiovascular complications with heavy impact on morbidity and prognosis.
Glucagon-like peptide 1 receptor (GLP-1R) agonists represent a relatively new class of anti-hyperglycemic agents, addressing most of the pathophysiological mechanisms involved in the development of T2DM. The main GLP-1R agonist actions are the stimulation of insulin secretion, the inhibition of glucagon secretion, the delay of the gastric emptying time and the stimulation of neogenesis of insulin-secreting cells [10]. Moreover, GLP-1R agonists show favourable effects on body weight and metabolic profile, with a lower risk of hypoglycaemia. Moreover, the favourable metabolic effects, together with a reduction in blood pressure (BP), contribute to reduce the cardiovascular risk. For all these reasons, the GLP1-R agonist use is constantly increasing to treat T2DM.
Among beneficial GLP-1R agonist effects, heart rate (HR) increase has been observed. Despite the resting HR increase could be considered a safety concern [11], GLP-1R agonist administration is associated to major adverse cardiovascular events reduction, including stroke, cardiovascular and allcause mortality, as suggested by the cardiovascular outcome trials (Supplementary Table 1) [12-18]. This GLP-1R agonist favourable effect on cardio-metabolic health is becoming increasingly evident, although the underlying mechanisms remain largely unknown. Improvement of hypertension [19], endothelial function [20] and a natriuretic GLP-1R agonist-related effect [21] have been proposed. However, while a small but substantial increase in HR by about 3 bpm has been reported by the majority of the available trials [22], the action of GLP-1R agonists on autonomic function in diabetes remains still debated. With this in mind, it is clear that the GLP-1R agonist effects on HR and on ANS need to be reconciled with the favourable cardiovascular outcomes in clinical trials. Thus, we performed this meta-analysis with the main aim to highlight the GLP-1R agonist actions on available measures of autonomic function in diabetic people.
METHODS
This meta-analysis was performed according to the Cochrane Collaboration and Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) statement. To ensure originality and transparency of the review process, the meta-analysis was a priori registered in the International Prospective Register of Systematic Reviews (PROSPERO; registration number CRD42020218063). The literature search was performed until October 10th, 2021 considering the following string: (((autonomic function)) OR ((autonomic dysfunction)) OR (cardiac autonomic neuropathy) AND ((((((((GLP-1 receptor agonist)) OR (GLP-1R agonist)) OR (semaglutide)) OR (liraglutide)) OR (exenatide)) OR (lixisenatide) OR (dulaglutide)) OR (albiglutide). Medline, Embase, and Cochrane databases were considered.
All available GLP-1R agonists were considered potentially eligible, whether applied in clinical trial aiming at evaluating autonomic function before and after chronic administration in participants with diabetes. Thus, each parameter of autonomic function was considered before and after chronic treatment in the study groups. Moreover, when available, mean parameters were compared between study and control groups after treatment.
Endpoints
CARTs represent the gold standard in autonomic testing [2,23]. CARTs involve measuring autonomic responses through changes in HR and BP to provocative physiological manoeuvres. The standard CARTs recommended for diagnosis of CAN include: the deep breathing (DB) test (expiration/inspiration [E/I] ratio), the lying-to-standing (LS) test (30:15), the Valsalva manoeuvre (VM) and the BP response to standing [2]. E/I index from DB test represents the ratio between the 3 maximum and the 3 minimum RR intervals (the intervals between two consecutive eletrocardiogram R waves) in a cycle of expiration and inspiration. In LS test the maximum/minimum 30:15 value is the ratio of the longest RR interval measured between the 25th and 35th beat after the change of posture and the shortest RR interval measured between the 10th and 20th beat. The VM test indicates the ratio between the longest RR interval after exhalation and the shortest RR interval during exhalation. Other approaches in clinical research are currently available to evaluate CAN such as heart rate variability (HRV), baroreflex sensitivity, muscle sympathetic nerve activity, heart sympathetic imaging [24]. Non-invasive and widely used in clinical research, HRV provides key information about autonomic—parasympathetic and sympathetic—modulation of the cardio-vascular system. The measurement of HRV includes two domains. In the frequency domain, the components of the HRV obtained by spectral analysis consist of low frequency (LF) and high frequency (HF) indices. These indices provide information about both the sympathetic and parasympathetic influences on heart. Thus, the LF/HF ratio is the index of sympathetic-parasympathetic balance. Time domain measures of the normal RR intervals mainly include differences between the longest and the shortest RR intervals, the standard deviations of RR interval (SDNN), and the square root of the mean squared difference of successive RR intervals (RMSSD) [24].
Study selection and inclusion criteria
The literature search evaluated all clinical trials with the following inclusion criteria: (1) either interventional or observational; (2) in which the autonomic function was evaluated; (3) in people with diabetes; (4) treated with GLP-1R agonists. In particular, change of ANS function before and after chronic GLP-1R agonist treatment using measures of HRV was considered the main outcome. Treatment was considered chronic when GLP-1R agonist was administered for more than 4 weeks. Both participants with type 1 diabetes mellitus (T1DM) and T2DM were considered eligible. Randomization and presence of controls were not considered as inclusion criteria.
Data collection process and quality
Two authors (C.G. and D.S.) separately performed the literature search, collecting abstracts of each study. Each abstract was evaluated for inclusion criteria and data were extracted from each study considered eligible. C.G. and D.S. performed quality control checks on extracted data. For the literature search, the primary endpoint was the LF/HF ratio after chronic GLP-1R agonist administration. Secondary endpoints extracted were: HR, SDNN, and RSMSSD as measures of HRV and E/I ratio, 30:15, VM as parameters of CARTs.
The risk of bias was assessed independently by two investigators (C.G. and D.S.), using Cochrane risk-of-bias algorithm. In particular, the following criteria were considered for each included trial, (1) randomization method, (2) concealment of allocation, (3) presence or absence of blinding to treatment allocation, (4) presence or absence of blinding of outcome assessment, (5) potential incomplete data biases, and (6) reporting biases.
Data synthesis and analysis
Using the Review Manager (RevMan) 5.4 Software (version 5.4.1, The Nordic Cochrane Centre, The Cochrane Collaboration, Copenhagen, Denmark; 2014), continuous variables were comprehensively evaluated as inverse variance of mean variables. When data were reported in the original manuscript as median or logarithm, they were transformed in mean±standard deviation. Indeed, mean±standard deviation is required for the meta-analytic approach. However, since these parameters could be obtained using different approaches, the meta-analyses were performed using standard mean difference. Considering that studies included in the meta-analysis reported different treatment durations and different time-points in which endpoints were evaluated, we considered the last available time point for each trial (median 24 weeks; range, 12 to 72). The degree of heterogeneity among the studies included in each analysis was examined by inspecting both the scatter in the data points and the overlap in their confidence intervals (CIs), and by performing I2 statistics. The inverse variance with the fixed model was initially chosen, whereas the random effect model was preferred in case of I2 higher than 60%. Sensitivity analyses were performed, considering the type of diabetes and the GLP-1R agonist used, when possible. Values of P<0.05 were considered statistically significant.
Data availability statement
All data generated or analyzed during this study are included in this published article.
RESULTS
Manuscripts (n=596) were identified by the literature search (Fig. 1). After abstract evaluations, 16 studies were considered for the full text analysis (Fig. 1). Among these, ten studies were excluded (reasons reported in Fig. 1) and six studies were enrolled [25-30] in the final analysis (Fig. 1). Table 1 summarizes the characteristics of the studies included in the meta-analysis. Moreover, Table 2 shows the characteristics of the patient cohorts included in the analysis: a total of 182 individuals with diabetes were enrolled (110 men and 72 women), with a mean age of 54.68±10.7 years. Only two studies enrolled patients with type 1 diabetes mellitus.
Primary outcome
The LF/HF ratio did not differ after treatment in the entire sample evaluated (P=0.410), as well as considering exenatide (P=0.290) or liraglutide (P=0.930) treatment, separately (Fig. 2). This lack of LF/HF ratio difference after treatment was maintained also dividing studies according to the type of diabetes considered (standard mean difference, –0.03; 95% CI, –0.26 to 0.20; I2=0%; P=0.080) (Supplementary Fig. 1). Thus, chronic treatment with GLP-1R agonists seems not to affect LF/HF ratio, also after sensitivity analyses. This was further confirmed comparing study and control groups, showing no differences in the LF/HF ratio after treatment (mean difference, 0.18; 95% CI, –0.66 to 1.01; I2=0%; P=0.680) (Supplementary Fig. 2).
Secondary outcomes
HR was reported in all trials, for a total of 182 patients, evaluated before and after chronic treatment. The HR significantly increased after treatment (P<0.001), with a low heterogeneity rate (6%) (Fig. 3). This significant improvement remained considering the use of liraglutide or exenatide separately (Fig. 3). Since Cacciatori et al. [26] did not report a control group, this study was excluded when study and control groups were compared post-treatment. In this analysis, HR was not significantly different between study and control groups (mean difference, –0.86; 95% CI, –3.26 to 1.55; I2=16%; P=0.480) (Supplementary Fig. 3).
After chronic treatment, SDNN did not significantly change (standard mean difference, 0.08; 95% CI, –0.65 to 0.81; P=0.830) (Supplementary Fig. 4). Considering the study-control comparison, significantly lower SDNN values were detected in the study-compared to control groups (P=0.040) (Fig. 4).
RSMSSD was evaluated in six studies, showing no difference after treatment (standard mean difference, 0.57; 95% CI, –0.41 to 1.55; I2=94%; P=0.250) (Supplementary Fig. 5). Similarly, no difference between study and control groups were detected (standard mean difference, –0.38; 95% CI, –1.20 to 0.44; P=0.360) (Supplementary Fig. 6).
The E/I ratio did not change after chronic GLP-1R agonist administration (mean difference, –1.13; 95% CI, –3.64 to 1.38; I2=95%; P=0.380) (Supplementary Fig. 7). In the study vs. control group comparison, the E/I ratio was not significantly different (mean difference, 0.0; 95% CI, –0.004 to 0.004; P=0.999) (Supplementary Fig. 8).
The 30:15 value was significantly lower after treatment (standard mean difference, 0.50; 95% CI, 0.18 to 0.83; P=0.002), although this result was obtained considering only two studies (Fig. 5). Only Hansen et al. [27] reported LS mean difference in the comparison between study and control groups; thus, this analysis was not feasible.
The VM did not significantly change after treatment (mean difference, 0.04; 95% CI, –0.08 to 0.15; P=0.520) (Supplementary Fig. 9).
Finally, the risk of bias was evaluated (Supplementary Fig. 10), showing a good quality of the included studies, in terms of reporting and incomplete biases. On the contrary, a variable quality in terms of blinding was detected, reflecting the difficulty to perform a completely randomized, double-blind clinical trial in this setting.
DISCUSSION
Here, we confirm that chronic GLP-1R agonist administration increases HR, according to a previously available meta-analysis [22]. Systematically, the HR increase after chronic GLP-1R agonist administration is evident in five of six included trials. Indeed, only one study found no significant changes in HR after GLP-1R agonist administration, but, differently from the others, authors used a short-acting GLP-1R agonist, such as exenatide, for much longer time [28]. Thus, we could speculate that the GLP-1R agonist effect on HR could depend on the molecule used and the duration of the administration. Moreover, together with the HR increase, no significant change in other ANS-related parameters is evident in our meta-analysis. This result suggests that the chronic GLP-1R agonist administration may not influence the sympathetic and parasympathetic functions. Thus, we could speculate that the HR modifications induced by GLP-1R agonist is not consequence of sympathetic or parasympathetic stimulation, but other mechanisms should be involved.
GLP-1R agonists are increasingly used in clinical practice in diabetes, considering the wide range of positive effects on glucose homeostasis, body weight, BP, and the low risk of hypoglycaemia. However, an overall GLP-1R agonist effect on ANS is far from being elucidated, even with the meta-analytic approach. Considering each study separately, interesting results could be extracted. Jaiswal et al. [28] did not detect any change in autonomic function after 18 months of treatment with the short-acting GLP-1R agonist exenatide in patients with T2DM, evaluating either the gold-standard CARTs, such as DB and VM or measures of HRV, such as HR, LF/HF, SDNN, RMSSD. In particular, that study was characterized by different follow-up length (18 months) than other studies included in the analysis (range, 12 to 26 weeks). Therefore, we could not exclude that this difference represents a confounding factor, or a determinant of a kind of mitigation of HR increase over time. On the contrary, Kumarathurai et al. [31] reported detrimental reductions in several HRV indices in T2DM patients treated with liraglutide for 12 weeks. In particular, liraglutide decreased SDNN, RMSSD, and HF without changes in LF/HF ratio compared to placebo [31]. This SDNN decrease persisted despite body weight loss and metabolic parameter improvement, suggesting an impairment in vagal activity after treatment [31]. Moreover, in order to elucidate the role of the cardiac sympatho-vagal balance as a possible mediator of the reported HR increase in GLP-1R agonist treatment, Cacciatori et al. [26] performed CARTs and power spectral analysis of HRV in 28 T2DM patients after exenatide extended-release administration. The main result of this study is an expected HR increase together with an unexpected LF/HF decrease, suggesting the existence of a compensatory mechanism [26]. This phenomenon is characterized by a “shift” of the sympatho-vagal balance with reduction of the LF component and no change in parasympathetic tone (HF) [26]. In contrast, Nystrom et al. [30] demonstrated the absence of changes in sympathetic or parasympathetic activity evaluated by HRV in 62 T2DM individuals receiving 1.8 mg liraglutide once daily in comparison with 4 mg glimepiride once daily. Similarly, in the context of T1DM, changes neither in the HRV domains, nor in cardiac vagal tone and cardiac sensitivity to the baroreflex were observed in 19 patients treated with liraglutide for 26 weeks [25]. This result was confirmed by the post-hoc analysis of Hansen et al. [27], showing no changes in CAN measures after liraglutide 1.8 mg once-daily for 24 weeks. In this study, the E/I ratio declined significantly in both liraglutide and placebo groups.
Accordingly, our findings show no difference in the LF/HF ratio after treatment, considering both different molecules (exenatide and/or liraglutide) and both type of diabetes (T1DM and/or T2DM). Moreover, others HRV measures, such as SDNN and RSMSSD, do not change after chronic GLP-1R agonist administration. These results are confirmed considering liraglutide and exenatide separately, suggesting no differences between molecules. LF/HF ratio represents an index widely used in clinical practice for CAN evaluation in diabetic patients, providing information about autonomic—parasympathetic and sympathetic—modulation of the cardiovascular system. Even if considered among the methods of investigation for cardiac autonomic dysfunction in human research studies [24], LF/HF ratio from HRV study confirms to be a measure not accurate and not directly related to sympatho-vagal balance, according to previous studies [32,33]. Among CARTs measurements, representing the gold standard for the diagnosis of CAN, the only parameter that seems to be influenced by GLP-1R agonist administration is the 30:15 ratio, which decreases after treatment. However, the strength of this result is limited by the small number of trials reporting this parameter. Thus, our meta-analysis suggests that chronic GLP-1R agonist treatment does not influence the sympatho-vagal balance in people with diabetes. Hence, the HR increase could depend on different mechanisms.
In animal models, GLP-1 engages GLP-1R in central, peripheral, and ANSs, enhancing the sympathetic nervous system activity, and reducing the parasympathetic nervous system activity [34]. In this regard, Baggio et al. [34] suggested that the GLP-1R agonist-related HR increase is the final effect of direct chronotropic action, which is attenuated by propanolol but not by atropine. Moreover, the in vivo GLP-1R agonist administration induces c-fos expression—a marker of neuronal activity—in the adrenal medulla, activates neurons involved in autonomic control in the brain, and activates tyrosine hydroxylase transcription in brainstem catecholamine neurons [35]. These findings suggest that the central GLP-1 action could be involved in the regulation of the sympathetic pathway [35]. However, the identification of GLP-1R expression in mouse atria [36] and in monkey sinoatrial node [37] raised the additional possibility that GLP-1 may directly modulate HR through the cardiac GLP-1R. Moreover, Berkelaar et al. [38] evaluated 130 healthy participants undergoing hyperglycaemic clamps and acute exposure to GLP-1 during hyperglycaemia, showing a small acute HR increase without an acute decrease in cardiac vagal control, measured by HRV. In this setting, serum insulin was positively associated with HR [38], suggesting that GLP-1 effect could be mediated by GLP1-driven increase in endogenous insulin [39]. Again, other experimental trials suggested the sympathetic nervous system activation after GLP-1R agonist infusion in healthy individuals [40,41]. All these examples demonstrate that the action of GLP-1R agonist on the sympathetic and parasympathetic systems must be both direct and indirect but this should be further studied with properly designed clinical and experimental trials. Indeed, the present meta-analytic approach is not able to reach strong conclusive results, due to the limited number of studies available. In particular, the scarcity of studies on the topic does not allow to perform further subgroup analyses, representing an important limit intrinsic to the meta-analysis itself. Thus, the lack of significant effects of GLP-1R agonist chronic administration on ANS should be considered in light with the limited amount of data available.
In conclusion, the present meta-analysis confirms the HR increase but seems to exclude an alteration of the sympatho-vagal balance due to chronic treatment with GLP-1R agonists in people with diabetes. Indeed, despite the accumulating data linking GLP-1R signalling to autonomic and neuroendocrine responses, the neural pathways underlying these actions are not fully understood. Furthermore, considering some discrepancies in the available preclinical and clinical findings, it is conceivable to suggest possible species-specific patterns of GLP-1R, as well as differences among GLP-1R agonists. More information is needed on the mechanisms through which the GLP-1R agonists administration may affect autonomic activity in individuals with diabetes.
SUPPLEMENTARY MATERIALS
Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2021.0314.
Notes
CONFLICTS OF INTEREST
No potential conflict of interest relevant to this article was reported.
AUTHOR CONTRIBUTIONS
Conception or design: C.G., D.S.
Acquisition, analysis, or interpretation of data: C.G., D.S.
Drafting the work or revising: C.G., D.S., G.B., M.S.
Final approval of the manuscript: C.G., D.S., G.B., C.P., M.S.
FUNDING
None
Acknowledgements
None