What this is
- This research investigates the causal relationship between glucagon-like peptide-1 receptor agonists (GLP1RAs) and ten common mental illnesses.
- Using , the study analyzes genetic data to understand potential neuropsychiatric benefits of GLP1RAs.
- The findings suggest a reduced risk of schizophrenia linked to GLP1RA use, with implications for clinical management.
Essence
- Genetically proxied activation of GLP1RAs is associated with a reduced risk of schizophrenia. However, evidence for other mental illnesses remains insufficient.
Key takeaways
- GLP1RAs are linked to a reduced risk of schizophrenia, with an of 0.84 (95% CI = 0.71β0.98). This finding suggests potential neuropsychiatric benefits beyond glycemic control.
- The study found no causal relationship between GLP1RAs and other mental illnesses, including attention deficit hyperactivity disorder, anorexia nervosa, autism spectrum disorder, and post-traumatic stress disorder. This indicates a need for further research.
Caveats
- The conclusions are based on genetic proxies, which may not directly translate to clinical outcomes due to the lifelong cumulative effects of genetic variants.
- The study population was limited to individuals of European ancestry, which may restrict the generalizability of the findings.
Definitions
- Mendelian randomization (MR): A genetic study method that assesses causal relationships between exposure and outcomes using genetic variants.
- Odds Ratio (OR): A measure of association between an exposure and an outcome, indicating the odds of the outcome occurring in the exposed group compared to the unexposed group.
AI simplified
1. Introduction
Glucagon-like peptide-1 (GLP-1) is an incretin hormone that exerts various physiological functions throughout the body, most notably in regulating blood glucose levels and controlling appetite [1]. Glucagon-like peptide-1 receptor agonist (GLP1RA) helps to regulate blood glucose levels by mimicking endogenous GLP-1 [2]. Semaglutide, originally developed for the treatment of diabetes, was approved by the FDA on 4 June 2021 as a weight-loss medication due to its ability to activate glucagon-like peptide-1 receptor (GLP1R) in the central nervous system [3], reducing food intake and promoting weight loss [4]. With advancements in research, GLP-1 has demonstrated promising potential in areas like cardioprotection and neuroprotection [5], generating significant interest in its effects beyond lowering of glucose and weight reduction.
Mental illnesses refer to a group of severe disorders that lead to brain dysfunction, causing significant public health concerns worldwide [6]. Affecting approximately 18% of the global population [7], they contribute to substantial mortality and economic burden [7,8]. Affected individuals experience disturbances in mood, cognition, or behavior. Mood improvement was initially observed in diabetic patients treated with GLP1RA. A large-scale cohort study of 10,690 diabetes patients over 6 to 7 years found that GLP1RA users had a lower incidence of depressive and anxiety disorders compared to non-users [9]. This observation sparked interest in the potential use of GLP1RAs for managing mental illnesses. GLP1Rs are widely distributed in the brain and play a role in regulating neuronal function [5,10,11], synaptic plasticity [12], and neurotransmitter release [13]. Although evidence from preclinical studies and clinical trials remains limited, liraglutide has been shown to improve executive function and memory in patients with major depressive disorder and bipolar disorder [14] as well as memory in individuals with prediabetes or early-stage type 2 diabetes and obesity [15]. In animal models, GLP1RAs have been found to enhance neurogenesis and regulate stress response circuits associated with mood disorders [16]. In addition, GLP1RAs have also demonstrated anti-inflammatory and antioxidant properties [13], which may offer potential benefits for a variety of mental illnesses. Overall, emerging evidence suggests that GLP1RA may have neuropsychiatric benefits. However, there is currently no direct evidence confirming the neuropsychiatric effects of GLP1RA in humans, and high-quality randomized controlled trials (RCTs) are still lacking. Most existing findings come from observational studies, and further research is still in its early stages. Recently, an observational study with a cohort of up to 2.4 million participants systematically examined the association between GLP1RAs and 175 health outcomes [17]. This study stands as the most extensive investigation of GLP1RAs to date. We took particular interest in the section on mental illnesses and found that the study reported that GLP1RA may reduce the risk of alcohol use disorder (AUD), cannabis use disorder (CUD), schizophrenia (SZ), and other mental health outcomes. Establishing a causal link between GLP1RAs and mental illness management could pave the way for new clinical indications. Equally important is the need to assess the potential psychiatric risks associated with long-term GLP1RA use. Therefore, further exploration of the causal relationships and underlying mechanisms between GLP1RAs and psychiatric disorders is warranted.
Mendelian randomization (MR) is a genetic study that assesses causal relationships between exposure and outcome [18]. Observational studies are often subject to biases such as confounding and reverse causation. In contrast, MR prevents these biases by utilizing the random allocation of genetic variants from parents to offspring during conception, thereby offering robust causal inference capabilities [19]. Based on rigorous instrumental variable selection procedure, drug target MR uses genetic variants to simulate the agonistic effects of drug targets, helping to evaluate their long-term exposure impact on diseases [20,21]. By selecting outcomes of interest, it can be used to study the potential efficacy and safety of new drug targets [22,23] as well as to explore the repurposing potential and adverse effects of existing drugs [20]. When RCTs are unavailable, drug target MR can provide high-level human genetic evidence to support drug repositioning and the identification of adverse effects [24]. We have observed that several MR studies have explored the associations between GLP1RA and conditions such as stroke, Alzheimer's disease, osteoarthritis, and cancer [25,26,27,28], contributing to a comprehensive evaluation of GLP1RA. Drug target MR can provide valuable insights into the clinical management of GLP1RA in relation to mental illnesses.
Overall, accumulating evidence suggests an association between GLP1RA and the risk of various mental illnesses, though the exact causal relationship remains unknown. In this context, we hypothesize that GLP1RA use may influence the risk of mental illness. This study aims to investigate the causal relationships between GLP1RA and 10 common mental illnesses: AUD, CUD, SZ, attention deficit and hyperactivity disorder (ADHD), anorexia nervosa (AN), anxiety disorder (ANX), autism spectrum disorder (ASD), bipolar disorder (BD), major depressive disorder (MDD) and post-traumatic stress disorder (PTSD). Specifically, we first perform drug-target MR analyses to examine the causal relationship between GLP1RA use and the aforementioned disorders. Additionally, for the identified causal relationships, we adopt MR mediation analysis to explore the potential mediating roles of glycemic control and weight management, which aids in understanding the underlying biological mechanisms of GLP1RA. This is the first systematic MR study investigating GLP1RA and the risk of mental illnesses. This research not only contributes to exploring the potential indications of GLP1RA but also provides valuable insights for assessing the risks of mental illnesses with their use.
2. Results
2.1. Selection of Genetic Instruments
The instrumental variables used to simulate the pharmacological effects of GLP1R activated in this study all have F-statistics greater than 40, indicating no bias caused by weak instruments. In the positive control analysis, GLP1RA significantly reduces the risk of T2D [OR (95% CI) = 0.79 (0.65β0.96), p = 0.02] and demonstrates an effect in lowering BMI [OR (95% CI) = 0.93 (0.88β0.98), p < 0.01]. These results further validate the effectiveness of the selected instrumental variables (Figure 1).
2.2. Causal Relationship Between GLP1RA and Mental Illnesses
Results for the causal effects of GLP1RA on mental illnesses are shown in Figure 2. In the MR analysis with outcomes data provided by the Psychiatric Genomics Consortium (PGC), the primary method indicated that GLP1RA was associated with an increased risk of AN [OR (95% CI) = 1.33 (1.05β1.68), p = 0.02] and a decreased risk of SZ [OR (95% CI) = 0.84 (0.70β0.99), p = 0.04]. In the FinnGen dataset, GLP1RA was associated with a decreased risk of ADHD [OR (95% CI) = 0.63 (0.44β0.89), p = 0.01], ASD [OR (95% CI) = 0.40 (0.18β0.91), p = 0.03], and PTSD [OR (95% CI) = 0.63 (0.43β0.94), p = 0.02]. MR-Egger, weighted median, weighted mode, and simple mode methods, which are based on different statistical assumptions, provided consistent estimates of effect direction (Supplementary Tables S1 and S2). This consistency enhances the reliability of the causal inference. There was no evidence of heterogeneity or horizontal pleiotropy (Supplementary Tables S3 and S4), indicating that our analysis was robust.
2.3. Meta-Analysis
We then focused on the outcomes that demonstrated at least one statistically significant association. Specifically, we performed a meta-analysis pooling MR results conducted separately in PGC and FinnGen. This helped to increase statistical power and the breadth of analysis while minimizing potential bias (Figure 3).
A causal relationship between GLP1RA and reduced SZ risk was found [OR (95% CI) = 0.84 (0.71β0.98), I2 = 0.0%, Common effects model]. This is considered as the primary finding of our study. However, there is insufficient evidence to support a causal relationship between GLP1RA and the risk of ADHD [OR (95% CI) = 0.86 (0.48β1.54), I2 = 89.6%, Random effects model], AN [OR (95% CI) = 0.84 (0.29β2.48), I2 = 79.1%, Random effects model], ASD [OR (95% CI) = 0.68 (0.30β1.52), I2 = 74.1%, Random effects model], and PTSD [OR (95% CI) = 0.82 (0.52β1.27), I2 = 72.6%, Random effects model].
2.4. Meditation Analysis
After identifying the causal relationship, conducting MR mediation analysis could help to discover underlying biological mechanisms. We primarily considered the potential mediating roles of T2D and BMI (Figure 4). We found no causal relationship between T2D and the risk of SZ [OR (95% CI) = 0.99 (0.95β1.04), p = 0.81], while a causal relationship between BMI and SZ was found [OR (95% CI) = 0.82 (0.75β0.89), p < 0.001]. In the framework of the two-step MR method, this indicates that the effect of GLP1RA in reducing SZ risk is unlikely to be mediated by glycemic control but may be related to BMI. The coefficient product test using the delta method further supports the mediation effect of BMI (p < 0.05). Detailed results are provided in Supplementary Tables S5 and S6.
3. Discussion
The intriguing link between neuroscience and endocrinology holds the potential to lead to groundbreaking therapies for mental illnesses. In recent years, there has been an increasing interest in the effects of GLP1RA beyond their role in regulating blood glucose. Although observed improvements in mood among diabetic patients receiving GLP1RA therapy have offered valuable insights, the ultimate effects and underlying mechanisms remain unclear. In this study, we performed a robust Mendelian randomization analysis and systematically investigated for the first time the causal relationship between GLP1RA and the risk of 10 mental illnesses. We found that genetically proxied GLP1R activation was associated with a reduced risk of schizophrenia and further mediation analysis suggests that this effect might not be related to glucose control but mediated by BMI.
Cognitive impairment is a core feature of schizophrenia [29], with current antipsychotic therapies primarily addressing positive symptoms, while cognitive deficits remain largely unaddressed [30]. Currently, the clinical application of GLP1RAs in schizophrenia management has primarily focused on their role in weight regulation. When co-administered with antipsychotic medications such as lithium, clozapine, and olanzapine, GLP1RAs have demonstrated the ability to mitigate metabolic disturbances and prevent cognitive deterioration [31]. However, research exploring the cognitive-enhancing effects of GLP1RA remains scarce and is predominantly confined to preclinical studies. Preclinical evidence has shown that GLP1RA can substantially improve spatial learning and memory [32,33,34,35,36,37,38]. In murine models, liraglutide has been observed to prevent MK-801-induced schizophrenia-like behaviors [39]. These findings support the potential benefits of GLP1RA for schizophrenia. Notably, GLP1RA has also been shown to enhance cognitive function in models of diabetes and neurodegenerative diseases, such as Alzheimer's disease [40]. This suggests that GLP1RA may have a direct effect on cognitive function independent of glucose normalization. In summary, these results provide compelling evidence that targeting GLP1R may have a positive impact on the clinical management of schizophrenia. This finding is consistent with the conclusions of large-scale observational studies [17]. This MR study suggests that the use of GLP1RA is associated with a reduced risk of schizophrenia, and the evidence from the animal models supports the potential of GLP1R as a therapeutic target for schizophrenia. Despite the genetic evidence, further studies are needed to ultimately confirm these findings and investigate the underlying molecular mechanisms.
Mediation analysis suggests that the mechanism by which GLP1RA reduces the risk of schizophrenia is related to weight management. Some genome-wide cross-trait analysis uncovers a shared genetic architecture between schizophrenia and BMI, indicating a complex biological connection between them [41,42,43]. While the molecular pathways of GLP1RA signaling have been identified, the direct mechanisms by which they enhance cognitive processes remain unclear. Current research suggests that these mechanisms may involve biological processes such as promoting neurogenesis and synaptic plasticity, reducing neuroinflammation, enhancing insulin signaling, preventing neuronal apoptosis, and minimizing oxidative stress [44]. The mechanisms by which GLP1RA exerts potential benefits on schizophrenia through weight management remain unclear, but clues can be found in the pathophysiology of schizophrenia. Weight loss is commonly associated with reduced systemic inflammation, and neuroinflammation is considered a key factor in the pathogenesis of schizophrenia. One possible mechanism is that GLP1RA may reduce inflammation through weight loss, thereby offering protective effects on the nervous system. This hypothesis is plausible. In current clinical practice, obesity has been associated with adverse psychiatric outcomes in individuals with SZ [45,46]. It is important to note that current studies on the cognitive-improving effects of GLP1RA are mainly based on diabetes and Alzheimer's disease models. Future explorations of the mechanisms should focus more on schizophrenia models to better understand the potential therapeutic effects of GLP1RA.
Another contribution of this study is the lack of evidence for a causal relationship between GLP1RA and the risk of other mental disorders. GLP1RAs are demonstrating significant development potential, and clinical evidence reporting the long-term usage risks for mental disorders remains relatively scarce. It is important to evaluate them as comprehensively as possible. In our study, there was insufficient evidence to support a causal relationship between GLP1RA use and other mental illnesses. Although in the analysis with outcomes from PGC, GLP1RA use was associated with an increased risk of AN, and in the analysis with outcomes from FinnGen, GLP1RA use was associated with a reduced risk of ADHD, ASD, and PTSD, the meta-analysis failed to reach statistical significance. These findings will be regarded as "suggestive evidence" and may provide promising directions for future in-depth studies. We noted that larger research interests are primarily focused on anxiety, depression, eating disorders, and addictive disorders. Improvement in anxiety and depression is of particular interest. Preclinical studies have provided preliminary evidence for the antidepressant and anxiolytic properties of GLP1RA [47], and several clinical studies have reported improvements in depressive symptoms following GLP1RA treatment [48,49]. However, there is also evidence reporting the emergence of anxiety and depressive symptoms associated with the use of semaglutide [50,51]. There is even evidence from the FDA Adverse Event Reporting System (FAERS) suggesting that the use of semaglutide may be associated with an increased risk of suicidal ideation [52]. This makes the relationship between GLP1RA use and anxiety and depression unclear. Additionally, our study found suggestive evidence indicating an increased risk of anorexia nervosa. In contrast, liraglutide is considered a safe and effective treatment for bulimia nervosa and binge eating disorder [53]. These findings suggest the need for continued attention to the effects of GLP1RA use in the management of eating disorders. Another promising research direction for GLP1RA is their potential in controlling the risk of addiction disorders. An observational study involving 2.4 million people supports this [17], but as they pointed out in their discussion, "observational studies cannot establish causality". Currently, there is not enough investigation to determine the exact effects of long-term use. In fact, the conclusions of the vast majority of current studies on the above diseases come from observational research, which cannot establish causality. Research on the psychiatric benefits of GLP1RAs is still in its early stages. Therefore, as a result of a MR study, we recommend conducting more large-scale, long-term follow-up, multi-center randomized controlled trials in the future to validate all of our conclusions. This will contribute to drug repositioning of GLP1RAs and the prediction of adverse reactions, ultimately translating into clinical practice. After confirming the causal relationship, further investigation into the mechanisms is needed to distinguish whether the neuropsychiatric effects of GLP1RAs are direct or mediated through the improvement of the original indications (T2D and obesity). This distinction may be important, as it may reveal a unique therapeutic mechanism of GLP1RAs for mental illnesses.
Our study has several strengths. First, to our knowledge, this is the first study to apply drug target MR to systematically investigate the effects of GLP1RA on mental illnesses, offering stronger causal inference compared to traditional observational studies. Second, due to the differences in GLP1R expression between mice and humans, it is challenging to extend preclinical research findings to clinical trials [54]. Our analysis, based on large-scale human GWAS summary-level data, can provide more valuable insights beneficial to clinical practice. Finally, triangulation across multiple databases and the design of positive controls further improved the design of drug target MR. And meta-analysis helps to avoid potential biases.
However, it must be noted that, like other MR studies on drug targets, our conclusions are subject to several limitations. First, unlike drug administration in clinical practice, which typically occurs over a shorter timeframe, genetic variants used as proxies for drug effects represent lifelong cumulative effects of small-scale perturbations to the drug target [55]. This means that the odds ratios presented in this study cannot be directly used to guide clinical practice. In fact, estimates from genetic proxy analyses are generally smaller than those observed in clinical practice [56]. Furthermore, when larger scale tissue-specific eQTL data become available, it will be important to explore the relationship between GLP1RA and the risk of mental illness across more tissues for a more in-depth evaluation. The study population was limited to individuals of European ancestry, so caution is needed when generalizing the findings to other populations.
4. Materials and Methods
4.1. Study Design
In this study, we implemented the following workflow: First, we obtained the drug target gene GLP1R for GLP1RA from DrugBank and ChEMBL, and we then used cis-eQTLs in blood to proxy the pharmacological effects of the target activation. A series of stringent criteria for instrumental variable selection were established. Subsequently, two sample MR analyses were conducted using psychiatric disorder datasets from two independent databases to explore potential causal relationships. Meta-analysis was employed to pool the results and draw final conclusions. Finally, mediation analyses were performed to investigate whether the discovered causal relationships were mediated by glycemic improvement or weight management. The overall study design is illustrated in Figure 5.
4.2. Data Source
Exposure and outcome data were derived from publicly available GWAS studies. To minimize potential bias caused by population diversity, the analysis was restricted to individuals of European ancestry. Table 1 provides further information on the genetic data resources. Further details on the disease endpoints are available in the corresponding publications. Ethics approvals were obtained by the original studies from relevant authorities, and informed consent was provided by all participants.
4.3. Genetic Instrumental Selection
Information regarding the pharmacologically active protein targets and corresponding encoding genes was retrieved from the DrugBank and the ChEMBL databases [71,72]. In this study, GLP1RA was selected as the exposure. Instrumental variables representing the expression of the drug target gene in blood was obtained from cis-eQTLs provided by eQTLGen. Specifically, SNPs significantly associated (p < 5 Γ 10β8) with the target gene GLP1R (chromosome 6, base position: 39,016,557β39,059,079) within a 100 kb upstream and downstream window were selected. Subsequently, clumping was performed with a linkage disequilibrium (LD) threshold of r2 < 0.3 [73,74]. Instrumental strength was ensured when the F-statistic exceeded 10, indicating a robust association with the exposure traits [75].
The decades-long use of GLP1RAs, principally acylated peptides such as liraglutide and semaglutide, for the treatment of obesity and T2D [76] has demonstrated their efficacy and safety. Accordingly, we designed a positive control MR analysis, where their expected effects on the outcomes can validate the effectiveness of the selected instrumental variables [77]. It is worth noting that the instrumental variables could not be directly identified in the summary GWAS data for BMI as the outcome in the MR analysis. To address this issue, we searched for proxies within a 100 kb window around the variants using an LD threshold of r2 > 0.85 to identify proxy SNPs for the analysis [74]. The setting of a positive control ensured the validity of the instrumental variables employed in this study.
4.4. MR and Sensitivity Analysis
We performed a two-sample MR analysis to estimate the association between GLP1RA and mental illnesses. In this study, inverse variance weighting (IVW) and Mendelian Randomization Robust Adjusted Profile Scores (MR.RAPs) were selected as the primary methods. The IVW method assumes balanced pleiotropy and applies a multiplicative random effects model to provide consistent estimates of causal effects [78]. In contrast, the MR.RAPs method utilizes robust adjusted profile scores for statistical inference, though it may still be biased due to unbalanced pleiotropy [79]. When the number of SNPs is less than 2, the Wald ratio is a feasible method. In addition, weighted mode [79], MR-Egger [80], weighted median [81], and simple mode were used as complementary methods.
We initially used mental illnesses data from PGC as the outcome and subsequently replicated the analysis using the summary-level GWAS data from the Finngen Consortium (R12) (https://www.finngen.fi/enβ). The FinnGen study is an extensive genomics initiative that has analyzed over 500,000 samples from the Finnish biobank, examining the correlation between genetic variation and health data to gain insights into disease mechanisms and predispositions [61]. This project involves Finnish research groups, biobanks, and global industry partners. In addition, all MR analyses were performed with sensitivity tests. We used Cochran's Q test [82] to assess heterogeneity and MR-Egger intercept [83] to detect horizontal pleiotropy.
4.5. Meta Analysis
Diseases that showed statistical significance in at least one of the two independent analyses for the same mental illness were further pooled using meta-analysis, which will help integrate causal conclusions across different databases. The results from the meta-analysis will serve as the primary basis for reporting the causal relationship between GLP1RA and specific mental illness outcomes in this study.
4.6. Mediation Analysis
Although the meta-analysis ultimately identified a causal relationship between GLP1RA and mental illness, the underlying mechanism remains unclear. It is necessary to consider whether the discovered causal relationship is related to the effects of GLP1RA on improving glucose levels or reducing weight; if not, it may suggest a direct effect of GLP1RA on the risk of mental illnesses.
To explore this, we conducted a preliminary investigation using multivariable Mendelian randomization [84], incorporating both the two-step MR method and the coefficient product test. The two-step method requires the existence of causal relationships between the exposure, mediator, and outcome in pairs, otherwise, the mediator may not exhibit a true mediation effect [85]. In the positive control, we identified a causal relationship between GLP1RA and both T2D and BMI. The causal relationship between the potential mediator and the target mental illness will be further considered. The selection of instrumental variables and the MR analysis procedure remain as described previously. If there is no causal relationship between the potential mediator and the outcome, it is unlikely that the causal effect from exposure to outcome is mediated by this potential mediator. If identifies causal relationship, we then apply the coefficient product test to further examine the mediation effect [86].
4.7. Statistical Analysis
Within our MR analysis, the results were quantified as odds ratios (ORs) with 95% confidence intervals (CIs), where the OR represents the change in susceptibility to mental illnesses per standard deviation increase in genetically predicted GLP1R expression. Given that IVW and MR.RAPs were the primary estimation models in this study, results were considered statistically significant when both methods reported p < 0.05 with consistent effect direction estimates across multiple analytical approaches. Additionally, conclusions were regarded as robust when no heterogeneity or pleiotropy was found (p > 0.05).
For the meta-analysis, we used I2 as an indicator of heterogeneity between incorporated studies. When it was below 50%, a common effect model was applied; otherwise, a random effects model was employed, which helped address heterogeneity in data analysis. If the OR obtained from the meta-analysis was consistently less than 1 or greater than 1, the exposure was considered associated with a reduction or increase in disease risk, which constitutes the main finding of this study.
The two-step MR method involves conducting MR analyses between exposure, outcome, and potential mediators in pairs. For each independent MR analysis, the statistical analysis methods are consistent with the standards previously described. The standard error used in the coefficient product test is generated by the delta method and statistical significance was determined by the p obtained from this test, with p < 0.05 indicating the presence of a mediation effect. All analyses were based on the "TwoSampleMR" R package (version 0.5.8).
5. Conclusions
This MR study suggests that GLP1RAs are causally associated with a reduced risk of schizophrenia, providing genetic evidence for their use in the clinical management of schizophrenia. This neuropsychiatric benefit may be associated with weight management effects. Our findings may also provide valuable guidance for selecting weight-loss medications for overweight individuals at high risk of schizophrenia.