A correction has been applied to this article in:
Corrigendum: Decreased Risk of Anxiety in Diabetic Patients Receiving Glucagon-Like Peptide-1 Receptor Agonist: A Nationwide, Population-Based Cohort Study
Wen-Hsuan Tsai
Fung-Chang Sung
Lu-Ting Chiu
Ying-Hsiu Shih2
Shu-I Wu- 1Division of Endocrinology and Metabolism, Department of Internal Medicine, Mackay Memorial Hospital, Taipei, Taiwan
- 2Management Office for Health Data (DryLab), Clinical Trial Research Center (CTC), China Medical University Hospital, Taichung, Taiwan
- 3Department of Health Services Administration, College of Public Health, China Medical University, Taichung, Taiwan
- 4Department of Food Nutrition and Health Biotechnology, Asia University, Taichung, Taiwan
- 5Department of Medicine, MacKay Medical College, New Taipei City, Taiwan
- 6Department of Psychiatry, Mackay Memorial Hospital, Taipei, Taiwan
Background: Previous findings on using Glucagon-like peptide-1 receptor agonist (GLP1-RA) as an antidepressant were conflicting, lacking large-scale studies. We used population-based data to investigate depression and anxiety risk in diabetic patients receiving the medication.
Methods: From claims records of the National Health Insurance Research Database (NHIRD) of Taiwan, we identified cohorts of 10,690 GLP1-RA users and 42,766 propensity score-matched patients without GLP1-RA use from patients with diabetes mellitus (DM) diagnosed in 2011–2017, matched by age, gender, index year, occupation, urbanization, comorbidities, and medications. Incidence, hazard ratios (HR) and 95% confidence interval (CI) of depression and/or anxiety were estimated by the end of 2017.
Results: The overall combined incidence of anxiety and/or depression was lower in GLP1-RA users than in non-users (6.80 versus 9.36 per 1,000 person-years), with an adjusted HR adjusted hazard ratio (aHR) of 0.8 (95% CI: 0.67–0.95) after controlling for covariates. The absolute incidence reduction was greater in anxiety (2.13 per 1,000 person-years) than in depression (0.41 per 1,000 person-years). The treatment effectiveness was significant for women. Patients taking GLP1-RA for longer than 180 days had the incidence of anxiety reduced to 2.93 per 1,000 person-years, with an aHR of 0.41 (95%CI: 0.27–0.61), compared to non-users. Dulaglutide could significantly decrease risks of both anxiety and depression.
Conclusion: Patients with DM receiving GLP1-RA therapy have a greater reduction of the risk of anxiety than that of depression. Our findings strengthen previous research that advocated possible anti-depressant or anxiolytic effects of GLP1-RA and may lead to improved treatment adherence among patients with DM.
Introduction
Diabetes mellitus (DM) is a group of endocrinological and metabolic disorders affecting approximately 422 million patients worldwide (Saeedi et al., 2019). Patients with DM are at an elevated risk of developing acute and serious long-term complications (Fowler, 2008; Gregg et al., 2016; Harding et al., 2019). Considerable evidence has also associated DM with increased risk of incident or recurrent depression (Gavard et al., 1993; Lustman et al., 2000; Anderson et al., 2001; Snoek et al., 2015; Chima et al., 2017; Khaledi et al., 2019). An earlier meta-analysis found that the risk of depression was 2-fold higher for patients with DM than individuals without DM (Anderson et al., 2001). Another meta-analysis review reported that poor glycemic control is associated with the occurrence of depression in patients with both type 1 and type 2 DM (Lustman et al., 2000).
Metformin, thiazolinediones (TZDs), glucagon-like peptide-1 receptor agonists (GLP1-RA), and dipeptidyl peptidase 4 (DPP-4) inhibitors are medications commonly prescribed for glycemic and weight control (Pozzi et al., 2019). A meta-analysis of clinical trials found that the use of pioglitazone, a TZD, was associated with significant improvements in depressive symptoms. The effect was more marked in women. The treatment effectiveness of metformin for diabetic patients was not consistent among studies (Moulton et al., 2018). However, a recent study in Saudi Arabia found that metformin could lower the probability of major depressive symptoms by 70% among patients with polycystic ovary syndrome (PCOS), but did not influence the occurrence of anxiety (AlHussain et al., 2020).
GLP1-RAs can cross the blood brain barrier, exerting function at both peripheral and central systems (Pozzi et al., 2019). A meta-analysis advocated that GLP1-RAs could exert antidepressant or anxiolytic effects on reducing the depression rating score of −2.09 (95% CI −2.28 – −1.91, p < 0.001) for diabetic patients (Pozzi et al., 2019). A United Kingdom study evaluating changes of quality of life for patients with GLP1-RA demonstrated that the therapy significantly reduced the Hospital Anxiety and Depression Scale (HADS) scores, compared with insulin-treated patients. (Grant et al., 2011). A placebo-controlled single-blind study, investigating the effect of GLP1-RA on excessive daytime sleepiness with 16 male DM patients, demonstrated a significant reduction in depression scores after the treatment compared to the baseline scores. Nonetheless, the reduction was not significant compared to placebo (Idris et al., 2013). A 26-weeks randomized controlled trial in the Netherlands found that GLP1-RA treatment in 26 patients significantly improved the quality of life measured by the Problem Areas in Diabetes Scale (PAID) questionnaire, compared to 24 patients using standard treatment (de Wit et al., 2014). Further follow-up study revealed that the effectiveness was sustained to 52 weeks, but there was insignificant improvement in scores of Beck Depression Index (BDI) (de Wit et al., 2016). Another 6-month treatment with liraglutide for 19 women with polycystic ovary syndrome without DM revealed no significant association with reduced depression, compared to 17 controls matched by age (Kahal et al., 2019). A pooled analysis on 5,325 individuals with higher BMI without diabetes evaluating the effectiveness of trials using GLP1-RA for 32–160 weeks also found no significant differences in depression (2.1 versus 2.1 events/100 person-years) and anxiety (1.9 versus 1.7 events/100 person-years) (O'Neil et al., 2017).
Most previous studies on GLP1-RA were undertaken with small sample sizes and reported inconsistent findings on the treatment effectiveness related to risks of depression and anxiety. These studies ascertained conditions of depression or anxiety mainly based on screening instruments rather than clinical diagnoses. A large long-term cohort study with sufficient sample size focusing on DM patients is lacking. Therefore, we used the insurance claims data of Taiwan to conduct this study to evaluate risks of anxiety and/or depression between DM patients with and without GLP1-RA, with sufficient large population-based cohorts and follow-up time.
Methods and Materials
Data Source
The National Health Insurance (NHI) of Taiwan is a compulsory health insurance system, launched in 1995 for all residents. For this study, we used the National Health Insurance Research Database (NHIRD) issued by NHI containing claims data of all insured individuals. The claims data provided records of birth date, gender, places of residence, and enrollee category, as well as diagnoses, drug prescriptions, and treatments for emergency and outpatient visits, and hospitalizations from 2000 to 2017. All identifications of insured people in the NHIRD were replaced with surrogate numbers and analyzed anonymously, hence consents were waived. This study was approved by the Institutional Review Board of Mackay Memorial Hospital (20MMHIS431e).
Study Cohorts
GLP1-RA was one of the available treatment options for DM patients, with approval from the NHI given in 2011. We therefore first identified 4,079,299 patients with the diagnosis of DM during the period of 2011–2017 (Figure 1). Patients who were less than 20 years old or had been diagnosed with depression and/or anxiety at baseline were excluded. Patients who were newly diagnosed in 2011–2017 with DM (ICD-10-CM code E08, E09, E10, E11, and E13; ICD-9-CM codes: 250) at least twice were considered as the potential study population. The date of the initial diagnosis with DM was designated as the index date. Patients with DM prescribed with GLP1-RA (Liraglutide, Dulaglutide, Exenatide) in 2011–2017 were included in the GLP1-RA group. For the comparison group, propensity score-matched DM patients without GLP1-RA use were selected and matched at a 1:4 ratio. The propensity score was calculated by logistic regression using covariates of age, gender, index year, occupation, urbanization, comorbidities, and medications.
FIGURE 1
FIGURE 1. Selection scheme for study population.
Outcome
Both cohorts were followed to estimate follow-up person-years for each individual in both groups until incident depression (ICD-10-CM codes: F32.0, F32.1, F32.2, F32.3, F32.4, F32.5, F32.9, F33.0, F33.1, F33.2, F33.3, F33.40, F33.41, F33.42, F34.1, and F43.21; ICD-9-CM codes: 296.20, 296.21, 296.22, 296.23, 296.24, 296.25, 296.26, 296.30, 296.31, 296.32, 296.33, 296.34, 296.35, 296.36, 311, 300.4, 309.0, and 309.1) and/or anxiety (ICD-10-CM code: F41; ICD-9-CM codes: 300) were identified, or the individual withdrew from the insurance, or the end of 2017. Depression and anxiety were defined as having the diagnoses from at least two outpatient visits or one hospitalization to ensure the validity (Chien et al., 2007). The incidence rate of depression or anxiety was estimated by per 1,000 person-years.
Potential Risk Covariates
Variables that might have an association with the development of depression or anxiety included comorbidities, other medications and demographic variables of occupation and urbanization level of residential area. We classified residential areas into four levels, with level 1 as the most urbanized, and level 4 as the least urbanized. We were interested in baseline comorbidities (Read et al., 2017), including hyperlipidemia (ICD-10-CM code: E78; ICD-9-CM codes: 272), hypertension (ICD-10-CM codes: I10, I11, I12, I13, and I15; ICD-9-CM codes: 401, 402, 403, 404, and 405), coronary artery disease (ICD-10-CM codes: I20, I21, I22, I23, I24, and I25; ICD-9-CM codes: 410, 411, 412, 413, 414, and 429), chronic kidney disease (ICD-10-CM codes: N18 and N19; ICD-9-CM codes: 585, 586, and 593), stroke (ICD-10-CM codes: I60, I61, I62, and I63; ICD-9-CM codes: 430, 431, 432, and 433), heart failure (ICD-10-CM code I50; ICD-9-CM codes: 428), chronic obstructive pulmonary disease (COPD) (ICD-10-CM codes: J44, J45, and J47; ICD-9-CM codes: 491, 493, 494, and 496), connective tissue disease (ICD-10-CM codes L93, L94; ICD-9-CM codes: 695 and 701), obesity (ICD-10-CM code E66; ICD-9-CM codes: 278), malignancy (ICD-10-CM codes D49; ICD-9-CM codes: 239), liver cirrhosis (ICD-10-CM code: K74; ICD-9-CM codes: 571), alcoholic liver disease (ICD-10-CM code: K70; ICD-9-CM codes: 571) and adapted Diabetes Complications Severity Index (aDCSI). The aDCSI is a good measure of diabetes severity that uniquely incorporates a wide range of diabetes complications (Young et al., 2008; Chang et al., 2012). The other anti-glycemic agents, including metformin, thiazolinedione (TZD), sodium-glucose co-transporters 2 (SGLT2) inhibitors, acarbose, sulfonylurea, dipeptidyl peptidase 4 (DPP-4) inhibitors, insulin, and compound agents, were also possible confounders to consider in the analyses.
Statistical Analysis
Demographic characteristics and prevalence of comorbidities in the study and the comparison groups were compared and examined with Chi-square tests for categorical variables and t-tests for continuous variables. The Kaplan-Meier method was used to graphically describe the cumulative incidence of depression and anxiety during the 7-years follow-up period. Cox proportional hazards regression analysis was used to determine the GLP1-RA group to compare group hazard ratio (HR) and 95% confidence interval (CI) of depression or anxiety. Demographic variables, comorbidities, and medications were included in multivariable models to estimate the adjusted hazard ratio (aHR): Model 1 adjusted for demographic factors, Model 2 adjusted for Model 1 variables and comorbidities, and Model 3 adjusted for Model 2 variables and medications. We further evaluated the treatment effectiveness by the length of treatment, stratifying the medication courses into 3 periods, 30–90 days, 91–180 days, and >180 days. All statistical analyses were performed using STATA version 14.0 (StataCorp) and results with p values less than 0.05 considered as significant.
Results
A total of 10,690 DM patients prescribed GLP1-RA and 42,766 comparisons with non-users were identified from the NHIRD. Both groups had similar distributions of age, gender, occupation, urbanization and comorbidities, with a mean age of 53.33 years (SD 13.04) and 45.06% women (Table 1). Although, the percentage of using metformin was higher in non-GLP1-RA users.
TABLE 1
TABLE 1. Baseline demographic factors and comorbidities of study population according to GLP1-RA use.
The cumulative incidence of anxiety was 2.13% lower in GLP1-RA users than non-users (Log-rank test p < 0.001), whereas that of depression was not significantly different between the 2 groups (Figure 2). The overall incidence of depression and/or anxiety was lower in GLP1-RA users than non-users (6.80 versus 9.36 per 1,000 person-years), with an aHR of 0.8 (95% CI: 0.67–0.95) for users, after controlling for demographic factors, comorbidities and medications (Table 2). The difference of incidence rates between the two groups was greater for anxiety than depression. The aHRs of developing anxiety and depression for the GLP1-RA group, compared to non-users, were 0.78 (95% CI: 0.64–0.95) and 0.94 (95% CI: 0.72–1.23), respectively.
FIGURE 2
FIGURE 2. Cumulative incidence of depression and anxiety in diabetes patients with/without GLP1-RA.
TABLE 2
TABLE 2. Incidence of depression and anxiety in diabetes patients with/without GLP1-RA use.
The beneficial effect on anxiety (Table 3) was specific to patients between 40 and 60 years old, with aHR of 0.73 (95% CI: 0.55–0.96). Compared with men, women in both groups exhibited greater incidences of depression and anxiety (Table 4). However, for women, risks of anxiety were significantly lower in GLP1-RA users than non-users. GLP1-RA users with either kind of comorbidity or that used other hypoglycemia agents showed no significant difference on the risk of depression (Table 4). On the other hand, GLP1-RA users without comorbidities except for hypertension showed significant reductions of anxiety (Table 3). GLP1-RA users with hypertension tended to have lower risk of anxiety, with aHR of 0.73 (95% CI: 0.54–0.97). GLP1-RA users that took metformin, sulfonylurea and oral medication combination had lower risk of anxiety, so as to those that did not use TZD, acarbose, SGLT2 inhibitors, DPP4 inhibitors, or insulin (Table 3).
TABLE 3
TABLE 3. Factors associated with anxiety incidence in diabetes patients with/without GLP1-RA: subgroup analysis.
TABLE 4
TABLE 4. Factors associated with depression incidence in diabetes patients with/without GLP1-RA: subgroup analysis.
The incidence of depression or anxiety decreased with increasing duration of treatment after the initiation of GLP1-RA medication (Table 5). The reduction trends were significant after controlling for all covariates. After taking the medicine for 180 days or longer, rates of incidence of depression or anxiety reduced to 2.19 and 2.93 per 1,000 person-years, respectively. The risk of having any depression or anxiety fell to an aHR of 0.5 (95% CI: 0.36–0.69) in GLP1-RA users, compared to non-users. We also evaluated the effect of different GLP1-RA on anxiety or depression in Tables 6 and 7. Our results revealed that dulaglutide could significantly reduce risks of anxiety and depression, while liraglutide and exenatide showed no significant reductions on risks of either anxiety or depression.
TABLE 5
TABLE 5. Risk of depression and anxiety associated with duration of GLP1-RA use.
TABLE 6
TABLE 6. Cox proportional hazard model estimated GLP1-RA hazard ratio of anxiety among different GLP1-RA types.
TABLE 7
TABLE 7. Cox proportional hazard model estimated GLP1-RA hazard ratio of depression among different GLP1-RA types.
Discussion
To our knowledge, this study represents the largest population-based analyses investigating whether GLP1-RA medication is associated with reduced risks of depression or anxiety in DM patients. GLP1-RA users exhibited a significant risk reduction for anxiety, and a moderate reduction for depression, compared with non-users. This treatment effectiveness on anxiety was observed in female users but not in male users. The effectiveness increased with the duration of medication and the significant risk reduction was observed after 6-months or longer therapy. Further age-specific stratified analyses revealed significant reduction of anxiety in patients between 40 and 60 years. GLP1-RA users taking metformin or sulfonylurea at the same time had decreased risk of anxiety, which was also noted in patients with hypertension. As for the specific effectiveness of each GLP1-RA, dulaglutide use could significantly decrease risks of both anxiety and depression, while liraglutide and exenatide showed no significant effect on reductions of either anxiety or depression.
Our findings on reduced risks of anxiety or depression associated with GLP1-RA use were consistent with some previous findings (Bode et al., 2010; Grant et al., 2011; Idris et al., 2013; de Wit et al., 2014), but contrast with others that showed no significant association between GLP1-RA medication and depression (de Wit et al., 2016; O'Neil et al., 2017; Kahal et al., 2019). However, studies with non-significant findings were either based on small sample sizes or were not focused on DM patients. Although Eren-Yazicioglu et al. (2021) reported that exenatide was associated with higher depressive scores indirectly through its effect on perceived stress, their study participants already had lifetime or current psychiatric diagnosis at baseline, which was different from our study (we excluded those with prior anxiety or depression diagnosis at baseline).
Neuroinflammation might be a critical factor for the onset, deterioration, relapse, and maintenance of depression or anxiety (Kopschina Feltes et al., 2017; Paudel et al., 2018; Woelfer et al., 2019). GLP-1 has been reported to promote productions of anti-inflammatory cytokines in various organs, including the adipose tissue, the pancreas, and the brain (Dobrian et al., 2011; Lee et al., 2012; Darsalia et al., 2014; Augestad et al., 2020; Reed et al., 2020). Kim et al. (2020) summarized that GLP1-RA may ameliorate depression by reducing neuroinflammation, balancing neurotransmitter homeostasis, promoting neuronal differentiation and neural stem cell proliferation and improving synaptic function. Albeit we failed to find significant reductions in the overall incidence of depression in GLP1-RA users, our subgroup analysis showed that patients using dulaglutide were at a lower risk of depression than non-users. GLP1-RA treatment-associated decrease in the risk of depression or anxiety in DM patients may be driven by combinations of anti-inflammation, better glycemic control, greater weight loss, and reduced concern about weight gain (Bode et al., 2010).
Women are generally more likely than men to develop an anxiety disorder (Leach et al., 2008; McLean and Anderson, 2009). The present study demonstrated that women were also at a higher incidence of anxiety risk than men even after using GLP1-RA. However, female users had a greater benefit from GLP1-RA therapy, than male users. The treatment effectiveness of GLP1-RA might be due to interactions between GLP-1R, the central and peripheral nervous system, estrogen, and the anorexic actions. GLP-1 is secreted from gut enteroendocrine cells and brain preproglucagon (PPG) neurons, which are known as the peripheral and central GLP-1 systems (Brierley et al., 2021). GLP-1 secreted by the intestine releases into the hepatoportal vein. This activates the vagus nerve to generate a neural signal towards the brain stem, such as the nucleus of the solitary tract nucleus (NTS) and the area postrema (AP), which send axons to the hypothalamus to release GLP-1 and activate the receptors. A new signal is then sent towards peripheral tissues through the autonomic nervous system (ANS) to regulate numerous functions (Cabou and Burcelin, 2011). However, research indicated that the intake-inhibitory effects of the GLP-1 RA, exendin-4 and liraglutide, were mediated by activation of GLP-1R expressed on sub-diaphragmatic vagal afferents as well as in the brain (Kanoski et al., 2011). That is, central and peripheral GLP-1 systems suppress eating via independent gut-brain circuits (Brierley et al., 2021). Central injections of GLP1-RA strongly decrease food intake, due to GLP-1R expression targeting the hypothalamic and brainstem inducing anorexic action. (Larsen et al., 1997; McMahon and Wellman, 1998; Hayes et al., 2008; Hayes et al., 2010; Kanoski et al., 2011). The same CNS regions are implicated for the anorexic action of estrogens (Palmer and Gray, 1986; Osterlund et al., 1998; Merchenthaler et al., 2004; Musatov et al., 2007), and may provide neuroanatomical grounds for the interaction between GLP-1 receptor and gender (Richard et al., 2016) in our results. Such a mechanism was also supported by Richard et al. (2016), who reported that women are more sensitive than men to the food reward impact of central GLP-1 receptor activation. In addition, individuals with obesity and overweight are more likely to have anxiety than non-obese persons (Amiri and Behnezhad, 2019). Since GLP1-RA is well known for its effect on weight loss, this might lead to the elimination of anxiety in women. Furthermore, younger DM patients with fewer comorbidities may be more concerned about body weight. They were probably more willing to adhere to the prescription of GLP1-RA and thus the risk of anxiety might be decreased.
Our study demonstrated that the combination of GLP1-RA with metformin exerted an anxiolytic effect, which was consistent with previous preclinical studies (Fan et al., 2019; Ji et al., 2019; Zemdegs et al., 2019; Turan et al., 2021). We also found that GLP1-RA users taking sulfonylurea had lower risk of anxiety, which may be due to the lower dose of sulfonylurea needed after using GLP1-RA, which could decrease the hypoglycemic risk of sulfonylurea. In our study, GLP1-RA users not using TZD also had a reduced risk of anxiety. This might be owing to the lowered side effect of weight gain from not using TZD. The finding that GLP1-RA users not using DPP-4 inhibitors had better reduction of anxiety was somewhat surprising. Numerous studies have reported that patients with mood and anxiety-related disorders are characterized by increased circulating inflammatory cytokines, including interleukin (IL)-1, IL-6, tumor necrosis factor (TNF), their soluble receptors, and acute phase reactants such as C-reactive protein (CRP) (Maes et al., 1992; Maes, 1999; Sluzewska, 1999; Michopoulos et al., 2017). DPP-4 is a novel adipokine secreted from adipose tissue (Lamers et al., 2011) with a pro-inflammatory role (Cordero et al., 1997). Obese patients may have higher DPP-4 levels from adipose tissue, and DPP-4 might induce depression by activating the immune system and promoting proinflammatory cytokine secretion. Previous studies revealed that although very low levels of the DPP-4 inhibitors were found in the brain (Fuchs et al., 2009; Fura et al., 2009), its neuroprotective effects might be more attributed to peripheral functions rather than directly in the central nervous system (Darsalia et al., 2013; Lin and Huang, 2016). On the other hand, GLP-1 RA crossed the BBB successfully (Hunter and Hölscher, 2012; Athauda and Foltynie, 2016). In Taiwan, combination usage of GLP1-RA and DPP-4 inhibitor is not reimbursed by health insurance. That is, either kind of medication has to be paid by the patient himself/herself. This may be another reason why GLP1-RA users not using DPP-4 inhibitors had better reduction of anxiety, since when the patient needs to pay for the other hypoglycemic agent, it may indicate that his or her glucose control is poorer, and may have more anxiety regarding health conditions compared to those who did not combine DPP-4 inhibitor with GLP1-RA.
Our study found different effects from each subtype of GLP1-RA on anxiety or depression, and dulaglutide was the only subtype associated with significant reduction in such risks. The distinct result may be due to differences in molecular structure that affects the potency, duration, and frequency of use (Monami et al., 2009; Buse et al., 2011). Dulaglutide may have better reductions of anxiety or depression due to its weekly injection formulation, when compared with liraglutide (daily injection). Exenatide and liraglutide also showed neuroprotective effect in few clinical studies among patients with Parkinson’s disease, Alzheimer’s disease, ischemia, traumatic brain injury, neuropathies, neurogenesis, or epilepsy, but not in patients with anxiety disorder as in our study. Erbil et al. (2019) suggested that not all GLP1-RAs have the same neuroprotective effect and the effect of GLP1-RAs on reversing neurodegenerative or neuro-destructive processes might be time and dose-dependent. Behavioral assessments in these preclinical studies are limited to the lesion model and may not cover all aspects of cortical function and all cortical layers. Hence, more studies investigating cortical effects of different subtypes of GLP-1 on depression or anxiety at the whole brain level, with longer assessment periods, and complete neuropsychological evaluations or behavioral observations are still needed both for preclinical and clinical studies.
The strength of this study lies in the very large number of patients used, the largest to-date, that focused on anxiety or depression associated with using GLP1-RA for the treatment for DM. We were able to conduct a longer duration follow-up than previous studies. The stratified data analysis was able to evaluate factors associated with treatment effectiveness. The identification of anxiety and depression were based on at least two clinical diagnoses rather than by screening instruments. However, several limitations should be acknowledged for this study. First, the study was observational, rather than a randomized control design, for the population of Taiwan. The generalizability of findings to other population may be restricted. Second, our efforts to match and adjust for possible confounding factors might be biased by unmeasured or unknown confounders not available in the NHIRD. Information on some potential confounders, including body weight, smoking, diet, exercise, stress, or family history, was not available in the NHIRD, the bias associated with unmeasured factors was unclear, despite our efforts to match study groups and controlling for available variables. Third, information on changes in body weight was unavailable in the NHIRD, we were unable to evaluate risks of anxiety or depression related to body weight changes.
Conclusion
This population-based study demonstrated that recipients of GLP1-RA treatment for DM were at a significantly decreased risk of anxiety. However, the effectiveness was significant only for female users and patients using the medicine for longer than 180 days. The treatment effectiveness was moderate in reducing depression risk, especially for patients that used dulaglutide. These findings strengthen previous research advocating the antidepressant or anxiolytic effects of GLP1-RA.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by The Mackay Memorial Hospital Committee Review Board approved the study protocol (IRB 20MMHIS431e approved with exempt review). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
Conceived and designed the analysis: W-HT, S-IW, M-CT, and F-CS. Contributed data or analysis tools: L-TC and Y-HS. Wrote the paper: W-HT, S-IW, and F-CS.
Funding
This study is supported in part by Taiwan Ministry of Health and Welfare Clinical Trial Center (No. MOHW109-TDU-B-212-114004), MOST Clinical Trial Consortium for Stroke (No. MOST 109-2321-B-039-002), China Medical University Hospital (No. DMR-108-154), Tseng-Lien Lin Foundation of Taichung, Taiwan.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alhussain, F., Alruthia, Y., Al-Mandeel, H., Bellahwal, A., Alharbi, F., Almogbel, Y., et al. (2020). Metformin Improves the Depression Symptoms of Women with Polycystic Ovary Syndrome in a Lifestyle Modification Program. Patient Prefer Adherence 14, 737–746. doi:10.2147/PPA.S244273
Amiri, S., and Behnezhad, S. (2019). Obesity and Anxiety Symptoms: a Systematic Review and Meta-Analysis. Neuropsychiatr 33, 72–89. doi:10.1007/s40211-019-0302-9
Anderson, R. J., Freedland, K. E., Clouse, R. E., and Lustman, P. J. (2001). The Prevalence of Comorbid Depression in Adults with Diabetes: a Meta-Analysis. Diabetes Care 24, 1069–1078. doi:10.2337/diacare.24.6.1069
Athauda, D., and Foltynie, T. (2016). The Glucagon-like Peptide 1 (GLP) Receptor as a Therapeutic Target in Parkinson's Disease: Mechanisms of Action. Drug Discov. Today 21, 802–818. doi:10.1016/j.drudis.2016.01.013
Augestad, I. L., Pintana, H., Larsson, M., Krizhanovskii, C., Nyström, T., Klein, T., et al. (2020). Regulation of Glycemia in the Recovery Phase after Stroke Counteracts the Detrimental Effect of Obesity-Induced Type 2 Diabetes on Neurological Recovery. Diabetes 69, 1961–1973. doi:10.2337/db20-0095
Bode, B. W., Testa, M. A., Magwire, M., Hale, P. M., Hammer, M., Blonde, L., et al. (2010). Patient-reported Outcomes Following Treatment with the Human GLP-1 Analogue Liraglutide or Glimepiride in Monotherapy: Results from a Randomized Controlled Trial in Patients with Type 2 Diabetes. Diabetes Obes. Metab. 12, 604–612. doi:10.1111/j.1463-1326.2010.01196.x
Brierley, D. I., Holt, M. K., Singh, A., De Araujo, A., Mcdougle, M., Vergara, M., et al. (2021). Central and Peripheral GLP-1 Systems Independently Suppress Eating. Nat. Metab. 3, 258–273. doi:10.1038/s42255-021-00344-4
Buse, J. B., Garber, A., Rosenstock, J., Schmidt, W. E., Brett, J. H., Videbæk, N., et al. (2011). Liraglutide Treatment Is Associated with a Low Frequency and Magnitude of Antibody Formation with No Apparent Impact on Glycemic Response or Increased Frequency of Adverse Events: Results from the Liraglutide Effect and Action in Diabetes (LEAD) Trials. J. Clin. Endocrinol. Metab. 96, 1695–1702. doi:10.1210/jc.2010-2822
Cabou, C., and Burcelin, R. (2011). GLP-1, the Gut-Brain, and Brain-Periphery Axes. Rev. Diabet Stud. 8, 418–431. doi:10.1900/RDS.2011.8.418
Chang, H. Y., Weiner, J. P., Richards, T. M., Bleich, S. N., and Segal, J. B. (2012). Validating the Adapted Diabetes Complications Severity Index in Claims Data. Am. J. Manag. Care 18, 721–726.
Chien, I. C., Kuo, C. C., Bih, S. H., Chou, Y. J., Lin, C. H., Lee, C. H., et al. (2007). The Prevalence and Incidence of Treated Major Depressive Disorder Among National Health Insurance Enrollees in Taiwan, 1996 to 2003. Can. J. Psychiatry 52, 28–36. doi:10.1177/070674370705200106
Chima, C. C., Salemi, J. L., Wang, M., Mejia De Grubb, M. C., Gonzalez, S. J., and Zoorob, R. J. (2017). Multimorbidity Is Associated with Increased Rates of Depression in Patients Hospitalized with Diabetes Mellitus in the United States. J. Diabetes Complications 31, 1571–1579. doi:10.1016/j.jdiacomp.2017.08.001
Cordero, O. J., Salgado, F. J., Viñuela, J. E., and Nogueira, M. (1997). Interleukin-12 Enhances CD26 Expression and Dipeptidyl Peptidase IV Function on Human Activated Lymphocytes. Immunobiology 197, 522–533. doi:10.1016/s0171-2985(97)80084-8
Darsalia, V., Hua, S., Larsson, M., Mallard, C., Nathanson, D., Nyström, T., et al. (2014). Exendin-4 Reduces Ischemic Brain Injury in normal and Aged Type 2 Diabetic Mice and Promotes Microglial M2 Polarization. PLoS One 9, e103114. doi:10.1371/journal.pone.0103114
Darsalia, V., Ortsäter, H., Olverling, A., Darlöf, E., Wolbert, P., Nyström, T., et al. (2013). The DPP-4 Inhibitor Linagliptin Counteracts Stroke in the normal and Diabetic Mouse Brain: a Comparison with Glimepiride. Diabetes 62, 1289–1296. doi:10.2337/db12-0988
De Wit, H. M., Vervoort, G. M., Jansen, H. J., De Galan, B. E., and Tack, C. J. (2016). Durable Efficacy of Liraglutide in Patients with Type 2 Diabetes and Pronounced Insulin-Associated Weight Gain: 52-week Results from the Effect of Liraglutide on Insulin-Associated wEight GAiN in Patients with Type 2 Diabetes' (ELEGANT) Randomized Controlled Trial. J. Intern. Med. 279, 283–292. doi:10.1111/joim.12447
De Wit, H. M., Vervoort, G. M., Jansen, H. J., De Grauw, W. J., De Galan, B. E., and Tack, C. J. (2014). Liraglutide Reverses Pronounced Insulin-Associated Weight Gain, Improves Glycaemic Control and Decreases Insulin Dose in Patients with Type 2 Diabetes: a 26 Week, Randomised Clinical Trial (ELEGANT). Diabetologia 57, 1812–1819. doi:10.1007/s00125-014-3302-0
Dobrian, A. D., Ma, Q., Lindsay, J. W., Leone, K. A., Ma, K., Coben, J., et al. (2011). Dipeptidyl Peptidase IV Inhibitor Sitagliptin Reduces Local Inflammation in Adipose Tissue and in Pancreatic Islets of Obese Mice. Am. J. Physiol. Endocrinol. Metab. 300, E410–E421. doi:10.1152/ajpendo.00463.2010
Erbil, D., Eren, C. Y., Demirel, C., Küçüker, M. U., Solaroğlu, I., and Eser, H. Y. (2019). GLP-1's Role in Neuroprotection: a Systematic Review. Brain Inj. 33, 734–819. doi:10.1080/02699052.2019.1587000
Eren-Yazicioglu, C. Y., Kara, B., Sancak, S., Uysal, S. P., Yazici, D., Okuroglu, N., et al. (2021). Effect of Exenatide Use on Cognitive and Affective Functioning in Obese Patients with Type 2 Diabetes Mellitus: Exenatide Use Mediates Depressive Scores through Increased Perceived Stress Levels. J. Clin. Psychopharmacol. 41, 428–435. doi:10.1097/JCP.0000000000001409
Fan, J., Li, D., Chen, H. S., Huang, J. G., Xu, J. F., Zhu, W. W., et al. (2019). Metformin Produces Anxiolytic-like Effects in Rats by Facilitating GABAA Receptor Trafficking to Membrane. Br. J. Pharmacol. 176, 297–316. doi:10.1111/bph.14519
Fowler, M. J. (2008). Microvascular and Macrovascular Complications of Diabetes. Clin. Diabetes 26, 77–82. doi:10.2337/diaclin.26.2.77
Fuchs, H., Binder, R., and Greischel, A. (2009). Tissue Distribution of the Novel DPP-4 Inhibitor BI 1356 Is Dominated by Saturable Binding to its Target in Rats. Biopharm. Drug Dispos 30, 229–240. doi:10.1002/bdd.662
Fura, A., Khanna, A., Vyas, V., Koplowitz, B., Chang, S. Y., Caporuscio, C., et al. (2009). Pharmacokinetics of the Dipeptidyl Peptidase 4 Inhibitor Saxagliptin in Rats, Dogs, and Monkeys and Clinical Projections. Drug Metab. Dispos 37, 1164–1171. doi:10.1124/dmd.108.026088
Gavard, J. A., Lustman, P. J., and Clouse, R. E. (1993). Prevalence of Depression in Adults with Diabetes. An Epidemiological Evaluation. Diabetes Care 16, 1167–1178. doi:10.2337/diacare.16.8.1167
Grant, P., Lipscomb, D., and Quin, J. (2011). Psychological and Quality of Life Changes in Patients Using GLP-1 Analogues. J. Diabetes Complications 25, 244–246. doi:10.1016/j.jdiacomp.2011.03.002
Gregg, E. W., Sattar, N., and Ali, M. K. (2016). The Changing Face of Diabetes Complications. Lancet Diabetes Endocrinol. 4, 537–547. doi:10.1016/S2213-8587(16)30010-9
Harding, J. L., Pavkov, M. E., Magliano, D. J., Shaw, J. E., and Gregg, E. W. (2019). Global Trends in Diabetes Complications: a Review of Current Evidence. Diabetologia 62, 3–16. doi:10.1007/s00125-018-4711-2
Hayes, M. R., De Jonghe, B. C., and Kanoski, S. E. (2010). Role of the Glucagon-Like-Peptide-1 Receptor in the Control of Energy Balance. Physiol. Behav. 100, 503–510. doi:10.1016/j.physbeh.2010.02.029
Hayes, M. R., Skibicka, K. P., and Grill, H. J. (2008). Caudal Brainstem Processing Is Sufficient for Behavioral, Sympathetic, and Parasympathetic Responses Driven by Peripheral and Hindbrain Glucagon-Like-Peptide-1 Receptor Stimulation. Endocrinology 149, 4059–4068. doi:10.1210/en.2007-1743
Hunter, K., and Hölscher, C. (2012). Drugs Developed to Treat Diabetes, Liraglutide and Lixisenatide, Cross the Blood Brain Barrier and Enhance Neurogenesis. BMC Neurosci. 13, 33. doi:10.1186/1471-2202-13-33
Idris, I., Abdulla, H., Tilbrook, S., Dean, R., and Ali, N. (2013). Exenatide Improves Excessive Daytime Sleepiness and Wakefulness in Obese Patients with Type 2 Diabetes without Obstructive Sleep Apnoea. J. Sleep Res. 22, 70–75. doi:10.1111/j.1365-2869.2012.01030.x
Ji, S., Wang, L., and Li, L. (2019). Effect of Metformin on Short-Term High-Fat Diet-Induced Weight Gain and Anxiety-like Behavior and the Gut Microbiota. Front. Endocrinol. (Lausanne) 10, 704. doi:10.3389/fendo.2019.00704
Kahal, H., Kilpatrick, E., Rigby, A., Coady, A., and Atkin, S. (2019). The Effects of Treatment with Liraglutide on Quality of Life and Depression in Young Obese Women with PCOS and Controls. Gynecol. Endocrinol. 35, 142–145. doi:10.1080/09513590.2018.1505848
Kanoski, S. E., Fortin, S. M., Arnold, M., Grill, H. J., and Hayes, M. R. (2011). Peripheral and central GLP-1 Receptor Populations Mediate the Anorectic Effects of Peripherally Administered GLP-1 Receptor Agonists, Liraglutide and Exendin-4. Endocrinology 152, 3103–3112. doi:10.1210/en.2011-0174
Khaledi, M., Haghighatdoost, F., Feizi, A., and Aminorroaya, A. (2019). The Prevalence of Comorbid Depression in Patients with Type 2 Diabetes: an Updated Systematic Review and Meta-Analysis on Huge Number of Observational Studies. Acta Diabetol. 56, 631–650. doi:10.1007/s00592-019-01295-9
Kim, Y. K., Kim, O. Y., and Song, J. (2020). Alleviation of Depression by Glucagon-like Peptide 1 through the Regulation of Neuroinflammation, Neurotransmitters, Neurogenesis, and Synaptic Function. Front. Pharmacol. 11, 1270. doi:10.3389/fphar.2020.01270
Kopschina Feltes, P., Doorduin, J., Klein, H. C., Juárez-Orozco, L. E., Dierckx, R. A., Moriguchi-Jeckel, C. M., et al. (2017). Anti-inflammatory Treatment for Major Depressive Disorder: Implications for Patients with an Elevated Immune Profile and Non-responders to Standard Antidepressant Therapy. J. Psychopharmacol. 31, 1149–1165. doi:10.1177/0269881117711708
Lamers, D., Famulla, S., Wronkowitz, N., Hartwig, S., Lehr, S., Ouwens, D. M., et al. (2011). Dipeptidyl Peptidase 4 Is a Novel Adipokine Potentially Linking Obesity to the Metabolic Syndrome. Diabetes 60, 1917–1925. doi:10.2337/db10-1707
Larsen, P. J., Tang-Christensen, M., Holst, J. J., and Orskov, C. (1997). Distribution of Glucagon-like Peptide-1 and Other Preproglucagon-Derived Peptides in the Rat Hypothalamus and Brainstem. Neuroscience 77, 257–270. doi:10.1016/s0306-4522(96)00434-4
Leach, L. S., Christensen, H., Mackinnon, A. J., Windsor, T. D., and Butterworth, P. (2008). Gender Differences in Depression and Anxiety across the Adult Lifespan: the Role of Psychosocial Mediators. Soc. Psychiatry Psychiatr. Epidemiol. 43, 983–998. doi:10.1007/s00127-008-0388-z
Lee, Y. S., Park, M. S., Choung, J. S., Kim, S. S., Oh, H. H., Choi, C. S., et al. (2012). Glucagon-like Peptide-1 Inhibits Adipose Tissue Macrophage Infiltration and Inflammation in an Obese Mouse Model of Diabetes. Diabetologia 55, 2456–2468. doi:10.1007/s00125-012-2592-3
Lin, C. L., and Huang, C. N. (2016). The Neuroprotective Effects of the Anti-diabetic Drug Linagliptin against Aβ-Induced Neurotoxicity. Neural Regen. Res. 11, 236–237. doi:10.4103/1673-5374.177724
Lustman, P. J., Anderson, R. J., Freedland, K. E., De Groot, M., Carney, R. M., and Clouse, R. E. (2000). Depression and Poor Glycemic Control: a Meta-Analytic Review of the Literature. Diabetes Care 23, 934–942. doi:10.2337/diacare.23.7.934
Maes, M., Lambrechts, J., Bosmans, E., Jacobs, J., Suy, E., Vandervorst, C., et al. (1992). Evidence for a Systemic Immune Activation during Depression: Results of Leukocyte Enumeration by Flow Cytometry in Conjunction with Monoclonal Antibody Staining. Psychol. Med. 22, 45–53. doi:10.1017/s0033291700032712
Maes, M. (1999). Major Depression and Activation of the Inflammatory Response System. Adv. Exp. Med. Biol. 461, 25–46. doi:10.1007/978-0-585-37970-8_2
Mclean, C. P., and Anderson, E. R. (2009). Brave Men and Timid Women? A Review of the Gender Differences in Fear and Anxiety. Clin. Psychol. Rev. 29, 496–505. doi:10.1016/j.cpr.2009.05.003
Mcmahon, L. R., and Wellman, P. J. (1998). PVN Infusion of GLP-1-(7-36) Amide Suppresses Feeding but Does Not Induce Aversion or Alter Locomotion in Rats. Am. J. Physiol. 274, R23–R29. doi:10.1152/ajpregu.1998.274.1.R23
Merchenthaler, I., Lane, M. V., Numan, S., and Dellovade, T. L. (2004). Distribution of Estrogen Receptor Alpha and Beta in the Mouse central Nervous System: In Vivo Autoradiographic and Immunocytochemical Analyses. J. Comp. Neurol. 473, 270–291. doi:10.1002/cne.20128
Michopoulos, V., Powers, A., Gillespie, C. F., Ressler, K. J., and Jovanovic, T. (2017). Inflammation in Fear- and Anxiety-Based Disorders: PTSD, GAD, and beyond. Neuropsychopharmacology 42, 254–270. doi:10.1038/npp.2016.146
Monami, M., Marchionni, N., and Mannucci, E. (2009). Glucagon-like Peptide-1 Receptor Agonists in Type 2 Diabetes: a Meta-Analysis of Randomized Clinical Trials. Eur. J. Endocrinol. 160, 909–917. doi:10.1530/EJE-09-0101
Moulton, C. D., Hopkins, C. W. P., Ismail, K., and Stahl, D. (2018). Repositioning of Diabetes Treatments for Depressive Symptoms: A Systematic Review and Meta-Analysis of Clinical Trials. Psychoneuroendocrinology 94, 91–103. doi:10.1016/j.psyneuen.2018.05.010
Musatov, S., Chen, W., Pfaff, D. W., Mobbs, C. V., Yang, X. J., Clegg, D. J., et al. (2007). Silencing of Estrogen Receptor Alpha in the Ventromedial Nucleus of Hypothalamus Leads to Metabolic Syndrome. Proc. Natl. Acad. Sci. U S A. 104, 2501–2506. doi:10.1073/pnas.0610787104
O'neil, P. M., Aroda, V. R., Astrup, A., Kushner, R., Lau, D. C. W., Wadden, T. A., et al. (2017). Neuropsychiatric Safety with Liraglutide 3.0 Mg for Weight Management: Results from Randomized Controlled Phase 2 and 3a Trials. Diabetes Obes. Metab. 19, 1529–1536. doi:10.1111/dom.12963
Osterlund, M., Kuiper, G. G., Gustafsson, J. A., and Hurd, Y. L. (1998). Differential Distribution and Regulation of Estrogen Receptor-Alpha and -beta mRNA within the Female Rat Brain. Brain Res. Mol. Brain Res. 54, 175–180. doi:10.1016/s0169-328x(97)00351-3
Palmer, K., and Gray, J. M. (1986). Central vs. Peripheral Effects of Estrogen on Food Intake and Lipoprotein Lipase Activity in Ovariectomized Rats. Physiol. Behav. 37, 187–189. doi:10.1016/0031-9384(86)90404-x
Paudel, Y. N., Shaikh, M. F., Shah, S., Kumari, Y., and Othman, I. (2018). Role of Inflammation in Epilepsy and Neurobehavioral Comorbidities: Implication for Therapy. Eur. J. Pharmacol. 837, 145–155. doi:10.1016/j.ejphar.2018.08.020
Pozzi, M., Mazhar, F., Peeters, G. G. A. M., Vantaggiato, C., Nobile, M., Clementi, E., et al. (2019). A Systematic Review of the Antidepressant Effects of Glucagon-like Peptide 1 (GLP-1) Functional Agonists: Further Link between Metabolism and Psychopathology: Special Section on "Translational and Neuroscience Studies in Affective Disorders". Section Editor, Maria Nobile MD, PhD. This Section of JAD Focuses on the Relevance of Translational and Neuroscience Studies in Providing a Better Understanding of the Neural Basis of Affective Disorders. The Main Aim Is to Briefly Summaries Relevant Research Findings in Clinical Neuroscience with Particular Regards to Specific Innovative Topics in Mood and Anxiety Disorders. J. Affect Disord. 257. Section Editor, Maria Nobile MD, PhD. This Section of JAD focuses on the relevance of translational and neuroscience studies in providing a better understanding of the neural basis of affective disorders. The main aim is to briefly summaries relevant research findings in clinical neuroscience with particular regards to specific innovative topics in mood and anxiety disorders. doi:10.1016/j.jad.2019.05.044
Read, J. R., Sharpe, L., Modini, M., and Dear, B. F. (2017). Multimorbidity and Depression: A Systematic Review and Meta-Analysis. J. Affect Disord. 221, 36–46. doi:10.1016/j.jad.2017.06.009
Reed, J., Bain, S., and Kanamarlapudi, V. (2020). Recent Advances in Understanding the Role of Glucagon-like Peptide 1. F1000Res 9, 239. doi:10.12688/f1000research.20602.1
Richard, J. E., Anderberg, R. H., López-Ferreras, L., Olandersson, K., and Skibicka, K. P. (2016). Sex and Estrogens Alter the Action of Glucagon-like Peptide-1 on Reward. Biol. Sex. Differ. 7, 6. doi:10.1186/s13293-016-0059-9
Saeedi, P., Petersohn, I., Salpea, P., Malanda, B., Karuranga, S., Unwin, N., et al. (2019). Global and Regional Diabetes Prevalence Estimates for 2019 and Projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9th Edition. Diabetes Res. Clin. Pract. 157, 107843. 9(th) edition. doi:10.1016/j.diabres.2019.107843
Sluzewska, A. (1999). Indicators of Immune Activation in Depressed Patients. Adv. Exp. Med. Biol. 461, 59–73. doi:10.1007/978-0-585-37970-8_4
Snoek, F. J., Bremmer, M. A., and Hermanns, N. (2015). Constructs of Depression and Distress in Diabetes: Time for an Appraisal. Lancet Diabetes Endocrinol. 3, 450–460. doi:10.1016/S2213-8587(15)00135-7
Turan, I., Sayan Özaçmak, H., Erdem, S., Ergenc, M., and Ozacmak, V. H. (2021). Protective Effect of Metformin against Ovariectomy Induced Depressive- and Anxiety-like Behaviours in Rats: Role of Oxidative Stress. Neuroreport 32, 666–671. doi:10.1097/WNR.0000000000001634
Woelfer, M., Kasties, V., Kahlfuss, S., and Walter, M. (2019). The Role of Depressive Subtypes within the Neuroinflammation Hypothesis of Major Depressive Disorder. Neuroscience 403, 93–110. doi:10.1016/j.neuroscience.2018.03.034
Young, B. A., Lin, E., Von Korff, M., Simon, G., Ciechanowski, P., Ludman, E. J., et al. (2008). Diabetes Complications Severity index and Risk of Mortality, Hospitalization, and Healthcare Utilization. Am. J. Manag. Care 14, 15–23.
Keywords: depression, anxiety, glucagon-like peptide-1 receptor agonist, diabetes, neuroendocrine
Citation: Tsai W-H, Sung F-C, Chiu L-T, Shih Y-H, Tsai M-C and Wu S-I (2022) Decreased Risk of Anxiety in Diabetic Patients Receiving Glucagon-like Peptide-1 Receptor Agonist: A Nationwide, Population-Based Cohort Study. Front. Pharmacol. 13:765446. doi: 10.3389/fphar.2022.765446
Received: 27 August 2021; Accepted: 24 January 2022;
Published: 23 February 2022.
Edited by:
Katarzyna Socała, Marie Curie-Sklodowska University, PolandReviewed by:
Hale Yapici Eser, Koç University, TurkeyJens Juul Holst, University of Copenhagen, Denmark
Copyright © 2022 Tsai, Sung, Chiu, Shih, Tsai and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shu-I Wu, t140@mmc.edu.tw
†These authors have contributed equally to this work