Long He
Qiuyu Li1,2†
Yilan Huang
Xiaoyan Zhong- 1Department of Pharmacy, The Affiliated Hospital of Southwest Medical University, Luzhou, China
- 2School of Pharmacy, Southwest Medical University, Luzhou, China
Aims: Glucagon-like peptide-1 receptor agonists (GLP-1 RAs) are employed extensively in the management of type 2 diabetes and obesity. However, there is a paucity of real-world data on their safety and tolerability for metabolic and nutritional adverse events in large sample populations. This study aimed to analyse the metabolic and nutritional safety signatures of different GLP-1 RAs by exploring the Food and Drug Administration (FDA) Adverse Event Reporting System (FAERS).
Methods: AEs data were extracted from the FDA Adverse Event Reporting System database for each GLP-1 RA from the time of its launch until the second quarter of 2023. The reported odds ratio (ROR), proportional reporting ratio (PRR), Empirical Bayesian Geometric Mean and Bayesian Confidence Propagation Neural Network were employed to identify AE signals.
Results: A system organ class of metabolism and nutrition disorders was employed to filter AE reports, resulting in the identification of 10,450 reports for exenatide, 2,860 reports for liraglutide, 240 reports for albiglutide, 4,847 reports for dulaglutide, 2,905 reports for semaglutide, 1,089 reports for tirzepatide, and 13 reports for lixisenatide. Semaglutide (ROR, 3.34; 95%CI, 3.22), liraglutide (ROR, 2.78; 95%CI, 2.69), and exenatide (ROR, 2.15; 95%CI, 2.11) were associated with metabolism and nutrition disorders. The number of AE signals detected were as follows: albiglutide (n = 1), lixisenatide (n = 2), tirzepatide (n = 11), exenatide (n = 12), liraglutide (n = 16), semaglutide (n = 20), dulaglutide (n = 22). Dehydration was the most frequent AE contributing to serious outcomes for liraglutide (n = 318, 23.93%), dulaglutide (n = 434, 20.90%), semaglutide (n = 370, 25.10%) and tirzepatide (n = 70, 32.86%). The time to onset (TTO) of AE was statistically different between exenatide and the other GLP-1 RAs (p < 0.001), and the Weibull parameters for dehydration for liraglutide, dulaglutide, and semaglutide analyses all showed an early failure-type profile.
Conclusion: Our study suggests that exenatide, liraglutide, and semaglutide are more susceptible to metabolic and nutritional AEs than other GLP-1 RAs. Liraglutide, dulaglutide, semaglutide, and tirzepaptide’s potential to induce dehydration, necessitates special attention. Despite certain deficiencies, GLP-1 RAs have considerable potential for the treatment of eating disorders.
1 Introduction
Obesity and diabetes constitute formidable challenge to public health, affecting upwards of one billion individuals worldwide (Bailey et al., 2024). The International Diabetes Federation (IDF) Diabetes Atlas reports that the prevalence of diabetes among individuals aged 20–79 stood at 10.5% in 2021, encompassing approximately 536.6 million people, with projections suggesting an escalation to 12.2% (783.2 million people) by 2045. Global healthcare expenditures related to diabetes in 2021 were estimated at 966 billion USD, and forecasted increase to 1,054 billion USD by 2045 (Sun et al., 2022). Glucagon-like peptide-1 receptor agonists (GLP-1 RAs) have emerged as a key therapeutic option in the management of type 2 diabetes mellitus (T2DM), owing to their efficacy in lowering blood glucose, in addition to providing cardiovascular and renal benefits, and ameliorating metabolic syndrome, as supported by a wealth of evidence-based research (Davies et al., 2022). GLP-1RAs increase satiety, reduce food intake, and delay gastric emptying (Drucker, 2018; Novikoff and Müller, 2024). These effects have further expanded the utility of GLP-1RAs in weight loss in obese patients (Tan et al., 2022; Bailey et al., 2023; Camilleri and Acosta, 2024). So far, the United States Food and Drug Administration (FDA) has approved liraglutide in December 2014 and semaglutide in 2020 for the treatment of obesity.
Some concerns have been raised among healthcare professionals over the safety profiles of GLP-1 Ras. There have been reports of suicidal or self-injurious behaviors (SSIBs) to the European Medicines Agency (EMA), although pharmacovigilance reviews have reported no direct link between GLP-1 RAs and SSIBs (Chen et al., 2023; Chen et al., 2024; Zhou et al., 2024). Increased risks of thyroid cancer and cholecystitis have been reported (He et al., 2022; Bezin et al., 2023).
Recent research has focused on the potential use of GLP-1RAs in treating eating disorders, specifically binge eating (McElroy et al., 2018; Aoun et al., 2024; Balantekin et al., 2024). Despite the appetite-suppressing effects of GLP-1 RAs, concerns have been raised regarding their impact on eating disorder symptoms and management (Bartel et al., 2024; Richards and Khalsa, 2024). Pharmacovigilance studies on the metabolic and nutritional adverse events of GLP-1 RAs are still ongoing. The aims of the present study was to identify potential signals of adverse metabolic and nutritional reactions associated with GLP-1 RA use. By doing so, we hope to provide valuable insights into the application and research of GLP-1 RAs in managing eating disorders and other metabolic conditions.
2 Materials and methods
2.1 Data sources
The data extraction and analysis procedure are shown in Figure 1. This retrospective pharmacovigilance study utilized data from the FDA Adverse Event Reporting System (FAERS) (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html), which consists primarily of adverse event reports reported by health professionals, pharmaceutical manufacturers, attorneys, and individual patients, reports of medication errors with medication administration, and product quality complaints. The database consists of seven data files: patient demographics (DEMO), drug information (DRUG), adverse event information (REAC), patient outcome information (OUTC), reported source information (RPSR), drug treatment date information (THER), and indication for drug (INDI) (Shu et al., 2022). Adverse events were classified using standardized Medical Dictionary for Regulatory Activities (MedDRA) terms for consistency and global recognition. GLP-1 RAs were analyzed as primary suspected (PS) from Q1 2004 to Q2 2023, with duplicate reports removed using recommended methods (Zhang et al., 2024). Data were collected for 16,964,230 patients and included 50,659,236 AE reports, focusing on AEs related to the system organ class (SOC) of metabolism and nutrition disorders. It is worth noting that the serious outcomes include death, life-threatening events, hospitalization, disability, congenital anomalies, and other serious consequences.
Figure 1. Flow chart of data extraction and analysis. A detailed description of the flow chart for data extraction and analysis of adverse events for GLP-1 RAs for metabolism and nutrition disorders in the United States Food and Drug Administration Adverse Event Reporting System (FAERS).
2.2 Disproportionate analysis
Disproportionality analysis is a commonly used method for post-marketing surveillance of adverse reaction signals in pharmaceuticals, which includes reporting ratio of ratios (ROR), proportional reporting ratios (PRR), Empirical Bayesian Geometric Mean (EBGM), and Bayesian Confidence Propagation Neural Network (BCPNN). Several studies have shown that Bayesian methods are superior to frequency counting approaches such as PRR, but EBGM and BCPNN are less sensitive and tend to detect fewer signals than ROR and PRR. The study used four approaches to detect adverse reaction signals (Evans et al., 2001; Szarfman et al., 2002). The system organ class (SOC), high level term (HLT) and preferred term (PT) tiers of the target drugs were detected, simultaneously fulfilling the detection conditions of the four approaches (ROR: lower limit of 95% CI > 1, PRR: χ2 ≧4, lower limit of 95% CI > 1, EBGM. EBGM05 > 2, BCPNN: IC025 > 0) were regarded as potential AE signals of GLP-1RAs. Formulas for disproportionate analysis shown in Supplementary Table S1. All data processing and statistical analysis performed using R (version 4.23), Microsoft Excel (version 2021).
2.3 Statistical methods
The study compared whether there were significant differences in proportion of gender, age, and weight across AEs between the serious and non-serious outcomes among the different medications (including albiglutide, exenatide, dulaglutide, liraglutide, semaglutide, and tirzepatide). Proportions were compared using Pearson’s chi-squared (χ2) or Fisher’s exact test, and the Mann-Whitney U test was applied for continuous non-normal distribution data, such as age and weight. Time to onset (TTO) was defined as the interval between the date of AE onset (EVENT_DT in the DEMO file) and the date of the start of drug use (START_DT in the THER file). Once data containing input errors (for example, where the EVENT_DT was earlier than the START_DT), inaccurate date inputs or missing specific data were excluded. Statistical analyses were performed to obtain the median time to onset of the target AEs, quartiles, and parameters of the Weibull test. The test results were analyzed with reference to previous studies in order to obtain the pattern of AE occurrence (Kinoshita et al., 2020; Mazhar et al., 2021). Data were analyzed using SPSS (version 22.0), and Minitab (version 21), and statistical significance was set at p < 0.05.
3 Results
3.1 Descriptive analysis
The AEs for metabolic and nutritional disorders associated with albiglutide (n = 240), exenatide (n = 10,450), dulaglutide (n = 4,847), liraglutide (n = 2,860), semaglutide (n = 2,905), tirzepatide (n = 1,089), and lixisenatide (n = 13), were screened according to their date of launch. Different GLP-1 RAs AE characteristics are summarized in Table 1. More metabolic and nutritional disease adverse reaction reports occurred in female (n = 13,817, 61.67%) compared to male patients (n = 7,684, 34.30%). More patients within a weight distribution of >100 kg (n = 3,523, 15.72%) with a median weight of 95 (IQR 81.6–112) kg. The age distribution of more patients was 45–65 years (n = 7,880, 35.17%), with a median age of 61 (IQR 53–68) years. The vast majority of the AEs were reported in patients from the United States (n = 19,807, 88.41%), and reports of AEs were submitted primarily by consumers (n = 17,775, 79.34%), followed by physicians (n = 2,328, 10.39%). Hospitalization (n = 3,961, 48.37%) and other serious (important medical events) (n = 3,395, 41.46%) were the most frequently reported serious outcomes.
Table 1. Reported characteristics of adverse events in metabolism and nutrition associated with GLP-1RAs.
3.2 Different GLP-1 RA-related signals
Three GLP-1 RAs displayed positive SOC signals in Table 2: semaglutide (ROR, 3.34; 95%CI, 3.22), liraglutide (ROR, 2.78; 95%CI, 2.69), and exenatide (ROR, 2.15; 95%CI, 2.11) showed an association with metabolism and nutrition AEs. In contrast, lixisenatide (ROR, 2.85; 95% CI, 1.66), tirzepatide (ROR, 2.11; 95% CI, 1.99), dulaglutide (ROR, 2.04; 95% CI, 1.99) and albiglutide (ROR, 0.60; 95% CI, 0.53), did not show any association.
Table 2. Signal detection for GLP-1 RAs-associated metabolic and nutritional adverse events.
Decreased appetite was the most frequently reported AE in the FAERS database, including exenatide (n = 6,125, 58.61%), dulaglutide (n = 2,331, 48.09%), liraglutide (n = 1,174, 41.05%), semaglutide (n = 1,377, 47.40%), and tirzepatide (n = 550, 50.51%), with the exception that it is also one of the therapeutic effects of GLP-1 RAs, therefore, decreased appetite was analyzed only at the HLT level (Supplementary Table S2).
We disproportionately analyzed each GLP-1RA separately at the PT level and identified positive signals according to the screening criteria (Figure 2). Albiglutide (n = 240) and lixisenatide (n = 9) were associated with low numbers of metabolic and nutritional AEs, with one and two signals detected, respectively. In contrast, exenatide, dulaglutide, liraglutide, semaglutide, and tirzepatide were associated with 12, 22, 16, 20, and 11 signals, respectively; Including episodes of food craving, hypoglycaemia, appetite disorder, increased appetite, and food aversion for all five drugs. Notably, IC025 > 3 was detected for eight PTs by the BCPNN approach: liraglutide: weight loss poor, food craving, increased appetite; semaglutide: weight loss poor, food aversion, food craving, increased appetite; tirzepatide: increased appetite.
Figure 2. Results of disproportionate analysis of AEs associated with metabolism and nutrition of seven GLP-1 RAs at the PT level. (A) ROR 95%CI, PRR 95%CI, EBGM05, IC025 for exenatide metabolism and nutrition-associated AEs. (B) ROR 95%CI, PRR 95%CI, EBGM05, IC025 for liraglutide metabolism and nutrition-associated AEs. (C) ROR 95%CI, PRR 95%CI, EBGM05, IC025 for dulaglutide metabolism and nutrition-associated AEs. (D) ROR 95%CI, PRR 95%CI, EBGM05, IC025 for semaglutide metabolism and nutrition-associated AEs. (E) ROR 95%CI, PRR 95%CI, EBGM05, IC025 for tirzepatide metabolism and nutrition-associated AEs. (F) ROR 95%CI, PRR 95%CI, EBGM05, IC025 for albiglutide metabolism and nutrition-associated AEs. (G) ROR 95%CI, PRR 95%CI, EBGM05, IC025 for lixisenatide metabolism and nutrition-associated AEs. (H) Overlap relationship between the positive signals of five GLP-1RAs at the pt levels. GLP-1RA, glucagon-like peptide-1 receptor agonist; AE, adverse event; PT, preferred term; ROR, reporting odds ratio; PRR, proportional reporting rate; EBGM05, lower one-sided 95% confidence limit (95% CI) of the empirical Bayes geometric mean; IC025, lower bound of the information component for the Bayesian confidence propagation neural network.
3.3 Severe versus non-severe cases
As shown in Table 3, with the exception of lixisenatide (n = 14), which was not included because the number of AEs was insufficient for analysis, differential analyses of gender, age, and body weight were performed between severe and non-severe cases in patients with AEs of metabolism and nutrition disorders who were treated with the remaining six GLP-1 RAs. The results showed statistically significant differences between two groups, included gender (χ2 = 13.776, p < 0.001), age (60 vs 62 years; p = 0.011), and body weight (94.43 vs 96.6 kg, p = 0.016) for exenatide; gender (χ2 = 25.577, p < 0.001) for liraglutide; gender (χ2 = 37.710, p < 0.001) and body weight (85 vs 96.15 kg, p = 0.003) for dulaglutide; age (61 vs 65 years; p < 0.001) and body weight (88 vs 93 kg, p = 0.018) for semaglutide; gender (χ2 = 11.664, p < 0.001) for tirzepatide.
Table 3. Differences in clinical characteristics of serious and non-serious reports.
A more detailed summary of the results in Supplementary Table S3 indicate that 1, 7, 11, 15, 13, and 3 AEs of the six GLP-1 RAs (albiglutide, exenatide, liraglutide, dulaglutide, semaglutide, and tirzepatide) tended to be reported as severe AEs (p < 0.05), respectively. Dehydration was the AE with the highest number of reports for severe outcomes of GLP-1RAs, including liraglutide (n = 318, 23.93%), dulaglutide (n = 434, 20.90%), semaglutide (n = 370, 25.10%), tirzepatide (n = 70, 32.86%). Of note, all AE outcomes were severe, including dulaglutide: ketoacidosis (n = 58), marasmus (n = 17), hyperglycaemic hyperosmolar nonketotic syndrome (n = 12), diabetic ketosis (n = 4); liraglutide: diabetic ketoacidosis (n = 104), hyperlipasaemia (n = 9), diabetic metabolic decompensation (n = 6), hyperamylasamia (n = 4); albiglutide: diabetic mellitus inadequate control (n = 23); semaglutide: euglycemic diabetic ketoacidosis (n = 32), ketoacidosis (n = 28), and ketosis (n = 8), starvation ketoacidosis (n = 4).
3.4 Time to onset analysis
With the exception of lixisenatide, which was excluded from the analysis due to insufficient validated data, the time to onset of metabolic and nutritional disorders was summarized in Figure 3. Median time to onset for exenatide was the longest at 54 (IQR: 17–152) days, and dulaglutide was the shortest at 5 (IQR: 0–27) days, 15.5 (IQR: 1–76.25) days for liraglutide, 6 (IQR: 1–41) days for albiglutide, 22 (IQR: 6–62) days for semaglutide and 10 (IQR: 1–38) days for tirzepatide. Cumulative distribution curves indicate the time to onset of metabolic and nutritional AEs after treatment with different GLP-1 RAs, and the difference was significant in the exenatide group compared to the other five GLP-1 RAs (p < 0.001).
Figure 3. Cumulative distribution curves demonstrating the onset time of glucagon-like peptide-1 receptor agonist (GLP-1 RA)-related metabolic and nutritional adverse events after treatment with different GLP-1 RAs.
Given that dehydration was the AE in which the most severe outcomes occurred after the use of multiple GLP-1 RAs, we performed an analysis of three GLP-1 RAs with sufficient data on time to onset in the occurrence of dehydration for analysis (Table 4), and the median time to onset was, in descending order, dulaglutide (5 days, IQR 0.5–38.5), liraglutide (20 days, IQR 4.25–64.75), and semaglutide (26 days, IQR 6–61). In the assessment of Weibull 3-parameter analysis, all shape parameters β and its upper limit of 95% CI were <1.
Table 4. Time to onset analysis of dehydration after treatment with three GLP-1RAs.
4 Discussion
We analyzed data on post-marketing AEs reported in the FAERS database related to the use of GLP-1 RA. Specifically, we adopted disproportionality analysis to identify AEs highly associated with GLP-1 RA post-treatment. We focused on AEs with severe outcomes, while also analyzing the time of onset of the six GLP-1 RAs. A Weibull 3-parameter analysis was employed to examine reports of severe outcomes associated with dehydration to identify temporal patterns of time to onset. Our study summarized the metabolic and nutritional class of AEs associated with GLP-1 RAs.
4.1 Signals of adverse events
We observed that metabolic and nutritional adverse reactions were more frequently reported by females (Table 1), which has previously been interpreted as a greater inclination for females to report AEs to FAERS (Becker, 2022; Chen et al., 2024). Semaglutide (ROR = 3.34) had a high risk of inducing metabolic and nutritional AEs and the most pronounced risk of gastrointestinal adverse effects, which may be related to the fact that the SOC tier of semaglutide had the highest disproportionate signal intensity (Liu et al., 2022).
Exenatide (n = 10,450), dulaglutide (n = 4,847), liraglutide (n = 2,860), semaglutide (n = 2,905), and tirzepatide (n = 1,089) had the greatest number of reported AEs, and the highest number of AEs reported on the effects of appetite, as reflected in our analyses of the HLT stratum (Supplementary Table S2). This is consistent with the fact that this class of drugs produces significant appetite suppression in response to central and peripheral GLP-1 R activation and has been actively used in the treatment of T2DM, obesity, and overweight (Shah and Vella, 2014; Coskun et al., 2018; Brunton and Wysham, 2020; Aldawsari et al., 2023; Bailey et al., 2023). Decreased appetite is common in GLP-1 RAs use and has been demonstrated not to be driven by gastric emptying and gastrointestinal adverse effects resulting from treatment with short- or long-acting GLP-1 RAs (Quast et al., 2021). Unfortunately, anorexia produced by the agitation of the central system GLP-1 receptor is undesired in lean T2DM patients (Borner et al., 2022). Disproportionate signals like increased appetite, food craving, hyperphagia, lack of satiety, and appetite disorder were also detected. Appetite disorder was mined in pharmacovigilance of adverse psychiatric events with GLP-1 RA (Chen et al., 2024), we hypothesized that other effects, such as increased appetite and anorexia, could have arisen during agonism of GLP-1 receptor in the nucleus of the solitary tract (NTS). In addition, different weight loss interventions have different effects on appetite, with low-calorie diets and exercise causing an increase in appetite, and medication (GLP-1 RAs) and bariatric surgery causing decreased appetite (Papamargaritis et al., 2024). Consequently, it is necessary to evaluate the positive and negative effects of GLP-1 RAs on appetite in different patients before using them.
Our study identified a potential risk of dehydration with dulaglutide, liraglutide, semaglutide, and tirzepatide, semaglutide emerged as having the highest risk of various fluid reductions (Supplementary Table S2). In animal experiments, it has been demonstrated that the endogenous GLP-1 receptor system is involved in the control of fluid intake and that, unlike food-generated stimuli, both endogenous GLP-1 and exogenous GLP-1 Ras inhibit fluid intake by activating the central nervous system GLP-1 receptor (McKay et al., 2011; McKay et al., 2014), all of which evidence suggests that GLP-1 RAs are associated with dehydration (Filippatos et al., 2014), which is consistent with the results of the two subsequent clinical trials (Winzeler et al., 2020; Winzeler et al., 2021). Meanwhile, the most frequently reported adverse effects of GLP-1 RAs are gastrointestinal, including nausea, vomiting, and diarrhea (Liu et al., 2022; Chong et al., 2024), all of which predispose patients to dehydration. However, current studies of fluid intake inhibition for GLP-1 RAs are based on routinely administered doses, and little is known about indications such as obesity and overweight for taking higher doses.
In previous real-world studies of GLP-1 RAs-induced hypoglycemia, lixisenatide showed the strongest association (Zhao et al., 2021), consistent with our finding (Figure 2). Not only that, but the risk of hyperglycaemic conditions and hypoglycaemic conditions with lixisenatide was also the highest in our study (Supplementary Table S2). Although GLP-1 RAs by themselves have a low risk of hypoglycemia, when used in combination with sulfonylureas and/or insulin, some experts recommend lowering the dose of sulfonylureas and short-acting and low-acting pancreatic islet analogs before or during the use of GLP-1 RAs (Romera et al., 2019; Smits and Van Raalte, 2021).
Corresponding to the signal intensity for weight loss poor of GLP-1 RAs in our study, a network meta-analysis showed a decreasing trend in weight loss efficacy for tirzepatide, semaglutide, and liraglutide (Pan et al., 2024). The weight loss effect produced by liraglutide in previous studies was individualized across subjects, with more significant fat mass reductions in females, and the starting point for measurement of the weight loss effect is recommended to start from baseline (Santini et al., 2023; Schultes et al., 2024). Differences in individual specificity and tolerance to GLP-1 RAs may be the main reason for differences in weight loss (Ruder, 2023), with racial differences also contributing (Dobbie et al., 2024). Meanwhile, the recent observation that previous use of anti-obesity medication is associated with poorer weight loss outcomes with semaglutide further expands our knowledge that GLP-1 RAs produce poor weight loss (Ghusn et al., 2024).
4.2 Severe versus non-severe cases
The most frequently reported adverse effects associated with GLP-1 RAs were gastrointestinal in nature. Practice guides suggest that these effects can be minimized by gradual dose escalation and effective prevention or mitigation during initiation of treatment, generally resulting in mild or moderate adverse events (Brunton and Wysham, 2020). In recent years, the incidence of cases associated with GLP-1 RAs has increased, as well as the rates of serious medical outcomes and healthcare facility admissions, especially for semaglutide and liraglutide (Gaw et al., 2024). In our study involving both severe and non-severe cases, the analysis of patient age and weight suggested a higher risk of severe metabolic and nutritional AEs in patients with younger age and lighter weight. Our study also shows that GLP-1 RA-associated metabolic and nutritional adverse events appear to be defined worse for women. We failed to find relevant studies to corroborate our findings on weight, age or gender, requiring clinical studies to validate our findings.
In the analysis of serious versus non-serious outcomes in cases, 1, 7, 11, 15, 13, and 3 AEs tended to be reported as serious AEs for albiglutide, exenatide, liraglutide, dulaglutide, semaglutide, and tirzepatide, respectively (Supplementary Table S3). Unexpectedly, dehydration produces the highest number of cases with serious outcomes in liraglutide, dulaglutide, semaglutide, and tirzepatide. It is apparent that the inhibition of fluid intake in the CNS produced by the administration of GLP-1 RAs and the loss of fluid induced by gastrointestinal adverse effects combine to contribute to the onset of dehydration. In elderly patients, who are at inherently higher risk of developing diabetes and dehydration (Atciyurt et al., 2024), low fluid intake and hypovolemia are the main etiologic factors contributing to severe dehydration outcomes in elderly patients (Lambert and Carey, 2023). Reduced fluid intake and gastrointestinal adverse effects due to GLP-1 RAs are potential triggers of acute kidney injury (Filippatos et al., 2014; Long et al., 2024), all of which contribute to the significant number of severe outcomes of dehydration following GLP-1 RA use. At the same time, there was no positive signal for dehydration and fluid intake inhibition with exenatide, suggesting that exenatide is the better drug option in patients at risk for dehydration.
In our study, AEs submitted by multiple GLP-1 RAs regarding ketoacidosis were disproportionate and nearly all AEs had a severe outcome when the effect of the combination was not considered (Supplementary Table S3). In a case report of ketoacidosis resulting from an injection of Wegovy 1.7 mg for the treatment of type 1 diabetes mellitus (T1DM), the effects of appetite suppression and inhibition of gluconeogenesis and glycogenolysis produced by GLP-1 RAs contributed to the reduction of the patient’s insulin push and induced the development of diabetic keto acidosis. Moreover, this could become more prevalent with the widespread use of GLP-1 RAs (Widhalm and Pulido, 2023). A meta-analysis of GLP-1 RAs as adjuvant therapy in patients with T1DM also observed that liraglutide was associated with higher odds of ketosis (OR 1.8; 95% CI, 1.1–2.8) (Park et al., 2023). Coupled with the slightly disproportionate reporting of diabetic ketoacidosis associated with GLP-1 RAs when they occur without combination with insulin (Yang et al., 2022), there is a need to be concerned about ketoacidosis as a potential adverse effect when using GLP-1 RAs.
4.3 Time to onset analysis
In TTO analysis, the median time to onset of GLP-1 RAs except exenatide (54 days, IQR 17–152), liraglutide (15.5 days, IQR 1–76.25), dulaglutide (5 days, IQR 0–27), albiglutide (6 days, IQR 1–41), semaglutide (22 days, IQR 6–62) and tirzepatide (10 days, IQR 1–38) all within 30 days, and the difference in onset of action between exenatide and other GLP-1 RAs was statistically significant (p < 0.001). The pharmacokinetic and pharmacodynamic characteristics of different GLP-1 RAs may be relevant (Chen et al., 2024), no studies have characterized the time to onset of metabolic and nutritional adverse effects of GLP-1 RAs, and we look forward to more clinical practice to validate our results in the future. The Weibull parameter is used to predict the time to onset and can serve to refine pharmacologic management thinking for patients in clinical practice (Mazhar et al., 2021; Shu et al., 2022). The shortest median time to onset of dehydrated AEs was for dulaglutide (5 days, IQR 0.5–38.5), followed by liraglutide (20 days, IQR 4.25–64.75) and semaglutide (26 days, IQR 6–61). Dehydration with three GLP-1 RAs showed early failure-type characteristics, suggesting that the risk of developing dehydration decreased progressively over time.
The desire to explore the effects associated with appetite and feeding prompted further differential analyses of whether metabolic and nutritional AEs are among the AEs with appetite disorders (Supplementary Table S4), which revealed that women appeared to be more susceptible to appetite-related side effects with the use of GLP-1 RAs. On the one hand, estradiol in women can mediate a reduction in food intake through multiple pathways (Eckel, 2011), an activational effect of estrogens may cause anorexia nervosa (Asarian and Geary, 2013). On the other hand, the literature reports a much greater prevalence of psychogenic eating disorders in women than in men in the United States, including anorexia nervosa, bulimia nervosa, and binge eating disorder (Hudson et al., 2007; Guerdjikova et al., 2019). Given the remarkable efficacy of GLP-1 RAs in the treatment of binge eating in several trials, although no GLP-1 RA has been approved for the treatment of any eating disorder symptom, GLP-1 RAs remain a highly promising drug for the treatment of bulimia at this time (Bartel et al., 2024). Based on other potential adverse effects of GLP-1 RAs and positive and/or negative effects on eating disorders (Richards and Khalsa, 2024), we suggest that careful assessment of a patient’s appetite status is necessary when considering GLP-1 RAs for the treatment of eating disorders, especially in female patients.
4.4 Limitations
In terms of the data received in the FAERS database, it is a spontaneous reporting system, and involving reports submitted by different reporters introduce the potential for duplicate records and variable quality of information. Despite our effect on data corrections and deletions of duplication, the large number of missing data services link the quality of data and its subsequent findings to varying degrees. Some of the detected disproportionate signals of adverse reactions overlapped somewhat with the clinical manifestations of the therapeutic drug indications, which prevented our study from distinguishing between them. We selected post-market real-world data of individual GLP-1 RAs for analysis, although this approach provides a more accurate description of the characteristics of each GLP-1 RA, it does not allow us to evaluate the overall intensity and magnitude of targeting a specific drug with disproportionate signaling effects. Finally, the FAERS-based disproportionality analysis reflects only an assessment of signal intensity, which is a statistical association that needs to be validated by further studies (Liu et al., 2022).
5 Conclusion
In this study, we quantitatively analyzed the relationship between each GLP-1 RA and metabolic and nutritional AEs based on the FAERS database at multiple levels and from various perspectives. Exenatide, liraglutide, and semaglutide were associated with metabolic and nutritional AEs, and multiple GLP-1 RAs novel AE signals were identified by disproportionate analysis. TTO analysis suggests metabolic and nutritional AEs of exenatide tend to develop later than other GLP-1 RAs. The occurrence of dehydration following the use of dulaglutide, liraglutide, and semaglutide tends to be within a month. However, the serious outcomes associated with the development of dehydration following the use of GLP-1 RAs warrant increased attention.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html.
Author contributions
LH: Data curation, Methodology, Resources, Software, Visualization, Writing–original draft, Writing–review and editing. QL: Data curation, Methodology, Writing–original draft, Writing–review and editing. YY: Data curation, Software, Visualization, Writing–original draft. JL: Data curation, Software, Writing–original draft. WL: Data curation, Software, Writing–original draft. YH: Conceptualization, Methodology, Supervision, Visualization, Writing–original draft, Writing–review and editing. XZ: Data curation, Methodology, Software, Supervision, Visualization, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
Thanks to the United States FDA for providing a free source of data for the study.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1416985/full#supplementary-material
Abbreviations
AEs, Adverse Events; BCPNN, Bayesian Confidence Propagation Neural Network; EBGM, Empirical Bayesian Geometric Mean; EMA, European Medicines Agency; FAERS, FDA Adverse Event Reporting System; FDA, United States Food and Drug Administration; GLP-1 RAs, Glucagon-like peptide-1 receptor agonists; HLT, High level term; IDF, The International Diabetes Federation; IQR, Interquartile Range; MedDRA, Medical Dictionary for Regulatory Activities; NTS, Nucleus of the solitary tract; PRR, Proportional reporting ratios; PT, Preferred term; ROR, Reporting ratio of ratios; SOC, System organ class; SSIBs, Suicidal or self-injurious behaviors; T1DM, Type 1 diabetes mellitus; T2DM, Type 2 diabetes mellitus; TTO, Time to onset.
References
Aldawsari, M., Almadani, F. A., Almuhammadi, N., Algabsani, S., Alamro, Y., and Aldhwayan, M. (2023). The efficacy of GLP-1 analogues on appetite parameters, gastric emptying, food preference and taste among adults with obesity: systematic review of randomized controlled trials. Diabetes Metab. Syndr. Obes. Targets Ther. 16, 575–595. doi:10.2147/DMSO.S387116
Aoun, L., Almardini, S., Saliba, F., Haddadin, F., Mourad, O., Jdaidani, J., et al. (2024). GLP-1 receptor agonists: a novel pharmacotherapy for binge eating (Binge eating disorder and bulimia nervosa)? A systematic review. J. Clin. Transl. Endocrinol. 35, 100333. doi:10.1016/j.jcte.2024.100333
Asarian, L., and Geary, N. (2013). Sex differences in the physiology of eating. Am. J. Physiol. - Regul. Integr. Comp. Physiol. 305, R1215–R1267. doi:10.1152/ajpregu.00446.2012
Atciyurt, K., Heybeli, C., Smith, L., Veronese, N., and Soysal, P. (2024). The prevalence, risk factors and clinical implications of dehydration in older patients: a cross-sectional study. Acta Clin. belg. 79, 12–18. doi:10.1080/17843286.2023.2275922
Bailey, C. J., Flatt, P. R., and Conlon, J. M. (2023). An update on peptide-based therapies for type 2 diabetes and obesity. Peptides 161, 170939. doi:10.1016/j.peptides.2023.170939
Bailey, C. J., Flatt, P. R., and Conlon, J. M. (2024). Recent advances in peptide-based therapies for obesity and type 2 diabetes. Peptides 173, 171149. doi:10.1016/j.peptides.2024.171149
Balantekin, K. N., Kretz, M. J., and Mietlicki-Baase, E. G. (2024). The emerging role of glucagon-like peptide 1 in binge eating. J. Endocrinol. 262, e230405. doi:10.1530/JOE-23-0405
Bartel, S., McElroy, S. L., Levangie, D., and Keshen, A. (2024). Use of glucagon-like peptide-1 receptor agonists in eating disorder populations. Int. J. Eat. Disord. 57, 286–293. doi:10.1002/eat.24109
Becker, R. (2022). Gender and survey participation: an event history analysis of the gender effects of survey participation in a probability-based multi-wave panel study with a sequential mixed-mode design. Methods Data Anal. 16, 30. doi:10.12758/mda.2021.08
Bezin, J., Gouverneur, A., Pénichon, M., Mathieu, C., Garrel, R., Hillaire-Buys, D., et al. (2023). GLP-1 receptor agonists and the risk of thyroid cancer. Diabetes Care 46, 384–390. doi:10.2337/dc22-1148
Borner, T., Tinsley, I. C., Doyle, R. P., Hayes, M. R., and De Jonghe, B. C. (2022). Glucagon-like peptide-1 in diabetes care: can glycaemic control be achieved without nausea and vomiting? Br. J. Pharmacol. 179, 542–556. doi:10.1111/bph.15647
Brunton, S. A., and Wysham, C. H. (2020). GLP-1 receptor agonists in the treatment of type 2 diabetes: role and clinical experience to date. Postgrad. Med. 132, 3–14. doi:10.1080/00325481.2020.1798099
Camilleri, M., and Acosta, A. (2024). Newer pharmacological interventions directed at gut hormones for obesity. Br. J. Pharmacol. 181, 1153–1164. doi:10.1111/bph.16278
Chen, C., Zhou, R., Fu, F., and Xiao, J. (2023). Postmarket safety profile of suicide/self-injury for GLP-1 receptor agonist: a real-world pharmacovigilance analysis. Eur. Psychiatry 66, e99. doi:10.1192/j.eurpsy.2023.2474
Chen, W., Cai, P., Zou, W., and Fu, Z. (2024). Psychiatric adverse events associated with GLP-1 receptor agonists: a real-world pharmacovigilance study based on the FDA Adverse Event Reporting System database. Front. Endocrinol. 15, 1330936. doi:10.3389/fendo.2024.1330936
Chong, K., Chang, J. K. J., and Chuang, L. M. (2024). Recent advances in the treatment of type 2 diabetes mellitus using new drug therapies. Kaohsiung J. Med. Sci. 40, 212–220. doi:10.1002/kjm2.12800
Coskun, T., Sloop, K. W., Loghin, C., Alsina-Fernandez, J., Urva, S., Bokvist, K. B., et al. (2018). LY3298176, a novel dual GIP and GLP-1 receptor agonist for the treatment of type 2 diabetes mellitus: from discovery to clinical proof of concept. Mol. Metab. 18, 3–14. doi:10.1016/j.molmet.2018.09.009
Davies, M. J., Aroda, V. R., Collins, B. S., Gabbay, R. A., Green, J., Maruthur, N. M., et al. (2022). Management of hyperglycemia in type 2 diabetes, 2022. A consensus report by the American diabetes association (ADA) and the European association for the study of diabetes (EASD). Diabetes Care 45, 2753–2786. doi:10.2337/dci22-0034
Dobbie, L. J., Coelho, C., Mgaieth, F., Chauhan, K., Campbell, S., Shuriye, S., et al. (2024). Liraglutide 3.0 mg in the treatment of adults with obesity and prediabetes using real-world UK data: a clinical evaluation of a multi-ethnic population. Clin. Obes. 14, e12649. doi:10.1111/cob.12649
Drucker, D. J. (2018). Mechanisms of action and therapeutic application of glucagon-like peptide-1. Cell Metab. 27, 740–756. doi:10.1016/j.cmet.2018.03.001
Eckel, L. A. (2011). The ovarian hormone estradiol plays a crucial role in the control of food intake in females. Physiol. Behav. 104, 517–524. doi:10.1016/j.physbeh.2011.04.014
Evans, S. J., Waller, P. C., and Davis, S. (2001). Use of proportional reporting ratios (PRRs) for signal generation from spontaneous adverse drug reaction reports. Pharmacoepidemiol. Drug Saf. 10, 483–486. doi:10.1002/pds.677
Filippatos, T. D., Panagiotopoulou, T. V., and Elisaf, M. S. (2014). Adverse effects of GLP-1 receptor agonists. Rev. Diabet. Stud. RDS 11, 202–230. doi:10.1900/RDS.2014.11.202
Gaw, C. E., Hays, H. L., Kemp, C. A., Kistamgari, S., Spiller, H. A., Rine, N. I., et al. (2024). Glucagon-like peptide-1 receptor agonist cases reported to United States poison centers, 2017-2022. J. Med. Toxicol. Off. J. Am. Coll. Med. Toxicol. 20, 193–204. doi:10.1007/s13181-024-00999-x
Ghusn, W., Fansa, S., Anazco, D., Tama, E., Cifuentes, L., Gala, K., et al. (2024). The association between previous use of anti-obesity medication and semaglutide weight loss outcomes. Diabetes Obes. Metab. 26, 2167–2175. doi:10.1111/dom.15523
Guerdjikova, A. I., Mori, N., Casuto, L. S., and McElroy, S. L. (2019). Update on binge eating disorder. Med. Clin. North Am. 103, 669–680. doi:10.1016/j.mcna.2019.02.003
He, L., Wang, J., Ping, F., Yang, N., Huang, J., Li, Y., et al. (2022). Association of glucagon-like peptide-1 receptor agonist use with risk of gallbladder and biliary diseases: a systematic review and meta-analysis of randomized clinical trials. JAMA Intern. Med. 182, 513–519. doi:10.1001/jamainternmed.2022.0338
Hudson, J. I., Hiripi, E., Harrison, G., Pope, J., and Kessler, R. C. (2007). The prevalence and correlates of eating disorders in the national comorbidity survey replication. Biol. Psychiatry 61, 348–358. doi:10.1016/j.biopsych.2006.03.040
Kinoshita, S., Hosomi, K., Yokoyama, S., and Takada, M. (2020). Time-to-onset analysis of amiodarone-associated thyroid dysfunction. J. Clin. Pharm. Ther. 45, 65–71. doi:10.1111/jcpt.13024
Lambert, K., and Carey, S. (2023). Dehydration in geriatrics: consequences and practical guidelines. Curr. Opin. Clin. Nutr. Metab. Care 26, 36–41. doi:10.1097/MCO.0000000000000880
Liu, L., Chen, J., Wang, L., Chen, C., and Chen, L. (2022). Association between different GLP-1 receptor agonists and gastrointestinal adverse reactions: a real-world disproportionality study based on FDA adverse event reporting system database. Front. Endocrinol. 13, 1043789. doi:10.3389/fendo.2022.1043789
Long, B., Pelletier, J., Koyfman, A., and Bridwell, R. E. (2024). GLP-1 agonists: a review for emergency clinicians. Am. J. Emerg. Med. 78, 89–94. doi:10.1016/j.ajem.2024.01.010
Mazhar, F., Battini, V., Gringeri, M., Pozzi, M., Mosini, G., Marran, A. M. N., et al. (2021). The impact of anti-TNFα agents on weight-related changes: new insights from a real-world pharmacovigilance study using the FDA adverse event reporting system (FAERS) database. Expert Opin. Biol. Ther. 21, 1281–1290. doi:10.1080/14712598.2021.1948529
McElroy, S. L., Mori, N., Guerdjikova, A. I., and Keck, P. E. (2018). Would glucagon-like peptide-1 receptor agonists have efficacy in binge eating disorder and bulimia nervosa? A review of the current literature. Med. Hypotheses 111, 90–93. doi:10.1016/j.mehy.2017.12.029
McKay, N. J., Galante, D. L., and Daniels, D. (2014). Endogenous glucagon-like peptide-1 reduces drinking behavior and is differentially engaged by water and food intakes in rats. J. Neurosci. 34, 16417–16423. doi:10.1523/JNEUROSCI.3267-14.2014
McKay, N. J., Kanoski, S. E., Hayes, M. R., and Daniels, D. (2011). Glucagon-like peptide-1 receptor agonists suppress water intake independent of effects on food intake. Am. J. Physiol. Regul. Integr. Comp. Physiol. 301, R1755–R1764. doi:10.1152/ajpregu.00472.2011
Novikoff, A., and Müller, T. D. (2024). Pharmacological advances in incretin-based polyagonism: what we know and what we don’t. Physiol. Bethesda Md 39, 142–156. doi:10.1152/physiol.00032.2023
Pan, X.-H., Tan, B., Chin, Y. H., Lee, E. C. Z., Kong, G., Chong, B., et al. (2024). Efficacy and safety of tirzepatide, GLP-1 receptor agonists, and other weight loss drugs in overweight and obesity: a network meta-analysis. Obes. Silver Spring Md 32, 840–856. doi:10.1002/oby.24002
Papamargaritis, D., le Roux, C. W., Holst, J. J., and Davies, M. J. (2024). New therapies for obesity. Cardiovasc. Res. 119, 2825–2842. doi:10.1093/cvr/cvac176
Park, J., Ntelis, S., Yunasan, E., Downton, K. D., Yip, T. C. F., Munir, K. M., et al. (2023). Glucagon-like peptide 1 analogues as adjunctive therapy for patients with type 1 diabetes: an updated systematic review and meta-analysis. J. Clin. Endocrinol. Metab. 109, 279–292. doi:10.1210/clinem/dgad471
Quast, D. R., Nauck, M. A., Schenker, N., Menge, B. A., Kapitza, C., and Meier, J. J. (2021). Macronutrient intake, appetite, food preferences and exocrine pancreas function after treatment with short- and long-acting glucagon-like peptide-1 receptor agonists in type 2 diabetes. Diabetes Obes. Metab. 23, 2344–2353. doi:10.1111/dom.14477
Richards, J. R., and Khalsa, S. S. (2024). Highway to the danger zone? A cautionary account that GLP-1 receptor agonists may be too effective for unmonitored weight loss. Obes. Rev. 25, e13709. doi:10.1111/obr.13709
Romera, I., Cebrián-Cuenca, A., Álvarez-Guisasola, F., Gomez-Peralta, F., and Reviriego, J. (2019). A review of practical issues on the use of glucagon-like peptide-1 receptor agonists for the management of type 2 diabetes. Diabetes Ther. 10, 5–19. doi:10.1007/s13300-018-0535-9
Ruder, K. (2023). As semaglutide’s popularity soars, rare but serious adverse effects are emerging. JAMA 330, 2140–2142. doi:10.1001/jama.2023.16620
Santini, S., Vionnet, N., Pasquier, J., Gonzalez-Rodriguez, E., Fraga, M., Pitteloud, N., et al. (2023). Marked weight loss on liraglutide 3.0 mg: real-life experience of a Swiss cohort with obesity. Obes. Silver Spring Md 31, 74–82. doi:10.1002/oby.23596
Schultes, B., Timper, K., Cavadini, G., Rüh, J., and Gerber, P. A. (2024). Weight loss and treatment patterns in a real-world population of adults receiving liraglutide 3.0 mg for weight management in routine clinical practice in Switzerland (ADDRESS study). Diabetes Obes. Metab. 26, 1023–1032. doi:10.1111/dom.15403
Shah, M., and Vella, A. (2014). Effects of GLP-1 on appetite and weight. Rev. Endocr. Metab. Disord. 15, 181–187. doi:10.1007/s11154-014-9289-5
Shu, Y., He, X., Wu, P., Liu, Y., Ding, Y., and Zhang, Q. (2022). Gastrointestinal adverse events associated with semaglutide: a pharmacovigilance study based on FDA adverse event reporting system. Front. Public Health 10, 996179. doi:10.3389/fpubh.2022.996179
Smits, M. M., and Van Raalte, D. H. (2021). Safety of semaglutide. Front. Endocrinol. 12, 645563. doi:10.3389/fendo.2021.645563
Sun, H., Saeedi, P., Karuranga, S., Pinkepank, M., Ogurtsova, K., Duncan, B. B., et al. (2022). IDF Diabetes Atlas: global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res. Clin. Pract. 183, 109119. doi:10.1016/j.diabres.2021.109119
Szarfman, A., Machado, S. G., and O’Neill, R. T. (2002). Use of screening algorithms and computer systems to efficiently signal higher-than-expected combinations of drugs and events in the US FDA’s spontaneous reports database. Drug Saf. 25, 381–392. doi:10.2165/00002018-200225060-00001
Tan, Q., Akindehin, S. E., Orsso, C. E., Waldner, R. C., DiMarchi, R. D., Müller, T. D., et al. (2022). Recent advances in incretin-based pharmacotherapies for the treatment of obesity and diabetes. Front. Endocrinol. 13, 838410. doi:10.3389/fendo.2022.838410
Widhalm, C., and Pulido, V. (2023). Euglycemic diabetic keto acidosis in a type 1 diabetic patient after glucose like peptide-1 administration: a case presentation. J. Investig. Med. High. Impact Case Rep. 11, 23247096231206333. doi:10.1177/23247096231206333
Winzeler, B., da Conceição, I., Refardt, J., Sailer, C. O., Dutilh, G., and Christ-Crain, M. (2020). Effects of glucagon-like peptide-1 receptor agonists on fluid intake in healthy volunteers. Endocrine 70, 292–298. doi:10.1007/s12020-020-02394-2
Winzeler, B., Sailer, C. O., Coynel, D., Zanchi, D., Vogt, D. R., Urwyler, S. A., et al. (2021). A randomized controlled trial of the GLP-1 receptor agonist dulaglutide in primary polydipsia. J. Clin. Invest. 131, e151800. doi:10.1172/JCI151800
Yang, Z., Yu, M., Mei, M., Chen, C., Lv, Y., Xiang, L., et al. (2022). The association between GLP-1 receptor agonist and diabetic ketoacidosis in the FDA adverse event reporting system. Nutr. Metab. Cardiovasc. Dis. 32, 504–510. doi:10.1016/j.numecd.2021.10.003
Zhang, M., Xie, W., Li, J., Zheng, J., and Zhou, Y. (2024). Postmarketing safety profile of brexanolone: a pharmacovigilance analysis based on FDA Adverse Event Reporting System (FAERS). Arch. Womens Ment. Health 27, 35–44. doi:10.1007/s00737-023-01378-1
Zhao, Z., Tang, Y., Hu, Y., Zhu, H., Chen, X., and Zhao, B. (2021). Hypoglycemia following the use of glucagon-like peptide-1 receptor agonists: a real-world analysis of post-marketing surveillance data. Ann. Transl. Med. 9, 1482. doi:10.21037/atm-21-4162
Zhou, J., Zheng, Y., Xu, B., Long, S., Zhu, L., Liu, Y., et al. (2024). Exploration of the potential association between GLP-1 receptor agonists and suicidal or self-injurious behaviors: a pharmacovigilance study based on the FDA Adverse Event Reporting System database. BMC Med. 22, 65. doi:10.1186/s12916-024-03274-6
Keywords: GLP-1 receptor agonists, FAERS, adverse events, metabolism and nutrition disorders, pharmacovigilance
Citation: He L, Li Q, Yang Y, Li J, Luo W, Huang Y and Zhong X (2024) Pharmacovigilance study of GLP-1 receptor agonists for metabolic and nutritional adverse events. Front. Pharmacol. 15:1416985. doi: 10.3389/fphar.2024.1416985
Received: 13 April 2024; Accepted: 18 June 2024;
Published: 08 July 2024.
Edited by:
Tin Wui Wong, Universiti Teknologi MARA Puncak Alam, MalaysiaReviewed by:
Harmanjit Singh, Government Medical College and Hospital, IndiaTimothy C. Hardman, Niche Science and Technology Ltd., United Kingdom
Copyright © 2024 He, Li, Yang, Li, Luo, Huang and Zhong. 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: Yilan Huang, aHlsMzE2MDEzMUAxMjYuY29t; Xiaoyan Zhong, MzEyNTU2NjYyQHFxLmNvbQ==
†These authors share first authorship