Identification and validation of antidiabetic drug targets
In accordance with the guidelines issued by the World Health Organization (WHO), we included six common classes of antidiabetic drugs (including metformin, sulfonylureas, GLP-1 receptor agonists, SGLT-2 inhibitors, DPP-4 inhibitors, and thiazolidinediones). The DrugBank (https://go.drugbank.com/) and ChEMBL (https://www.ebi.ac.uk/chembl/) databases were then used to identify target genes for six classes of antidiabetic drugs. For drugs with multiple target genes, we usually choose the one that has been widely studied and reported. We labelled two target loci as “ABCC8/KCNJ11” if they overlapped with each other and had the same IVs. Here we identified a total of six classes of target gene information for antidiabetic drugs, including chromosomal locations and gene loci. Specific details have been shown in Table 1. Metformin was excluded because of inconsistent information on metformin target genes obtained from the two databases.
Table 1
Target genes of antidiabetic drugs from DrugBank and ChEMBL databases
Drug class | Grug targets | Encoding genes | Gene location | Included |
Drugbank | ChEMBL | |
Sulfonylureas | ATP-sensitive potassium channel | ABCC8 | ABCC8 | Chr11: 17414045–17498392 | Yes |
KCNJ11 | KCNJ11 | Chr11: 17406795–17410893 |
Thiazolidinediones | Peroxisome proliferator-activated receptor gamma | PPARG | PPARG | Chr3: 12328867–12475843 | Yes |
Glucagon-like peptide-1 (GLP-1) analogues | Glucagon-like peptide 1 receptor | GLLPR | GLLPR | Chr6: 39016557–39059079 | Yes |
Sodium-glucose cotransporter (SGLT2) inhibitor | Sodium/glucose cotransporter 2 | SLC5A2 | SLC5A2 | Chr16: 31494444–31502090 | NO |
Dipeptidyl peptidase-4 (DPP-IV) inhibitor | Dipeptidyl peptidase IV | DPP4 | DPP4 | Chr16:162848755–162930725 | NO |
Metformin | NA | PRKAB1 | 58 genes | NA | NO |
Genetic instruments for drug target gene expression
Since we studied antidiabetic drugs, we used blood glucose as a biomarker. SNPs that were significantly associated with blood glucose and located within ± 100 kb of the target gene were selected as IVs (P < 5 × 10− 8, MAF > 0.01). The threshold for r2 was set at 0.3 in order to eliminate the linkage disequilibrium (LD). The strength of the genetic IVs was tested by calculating the F-value, and weak IVs with F-value < 10 would be excluded. Expressed as \(\:\text{F}\text{=}\frac{{R}^{2}\times\:(N-2)}{1-{R}^{2}}\), \(\:{R}^{2}=\frac{2\times\:{\beta\:}^{2}\times\:EAF\times\:\left(1-EAF\right)}{2\times\:{\beta\:}^{2}\times\:EAF\times\:\left(1-EAF\right)+2\times\:{SE}^{2}\times\:N\times\:EAF\times\:(1-EAF)}\). SNPs associated with blood glucose were provided by the Genome-Wide Association Study (GWAS) (ebi-a-GCST90025986, PMID: 34226706) from the IEU online database, this GWAS containing information on blood glucose of 400,458 European populations.
Positive control analysis and outcomes
The validity of IVs was confirmed by demonstrating the expected impact of antidiabetic drugs target genes on known outcomes. The intended effect of antidiabetic drugs is type 2 diabetes. In addition, thiazolidinediones improve insulin resistance, using the Insulin Resistance (IR) Index as a positive control. Sulfonylureas and GLP-1 receptor agonists promote insulin secretion, so peak insulin level was used as a positive control. Antidiabetic drugs have different effects on body weight, so body mass index (BMI), hip circumference, and waist circumference were chosen as positive controls. Here we have incorporated eight psychiatric disorders, including attention deficit hyperactivity disorder (ADHD), anorexia nervosa, anxiety disorder, autism spectrum disorder (ASD), bipolar disorder, major depressive disorder, obsessive compulsive disorder (OSB), schizophrenia。These genetic association data for positive controls and outcomes are available in the IEU database (https://gwas.mrcieu.ac.uk/).
Two sample MR analysis
Two sample MR method was used to analyze the causal relationship between antidiabetic drugs and psychiatric disorders, and the flow of the analysis was shown in Fig. 1. This approach must be based on three assumptions: (1) IVs are strongly correlated with exposure factors. (2) IVs can only have an effect on outcomes through exposure factors. (3) IVs are not associated with confounders. If multiple IVs were obtained, the inverse variance weighting (IVW) method was used (16). The IVW method is based on the fact that all SNPs are valid variables and independent of each other. IVW estimates the average effect of genetic variants on causality by means of a weighted linear regression model, and the results of IVW are unbiased if there is no horizontal pleiotropy (17). If only single IV was obtained, the Wald ratio test was used. The Wald ratio is the simplest method for MR analyses using single SNP as an IV without the necessity of heterogeneity and horizontal pleiotropy tests (18). IVW and MR Egger methods were used to detect heterogeneity of results. We concentrated on the Cochrane’s Q value, where a considerable heterogeneity is indicated if Q_P < 0.05 (19). To evaluate horizontal pleiotropy, MR Egger intercept method was employed. Significant horizontal pleiotropy is seen when P < 0.05 (20). Furthermore, the MR-PRESSO method was employed to identify any outliers (21). Finally, we applied the leave-one-out technique to a sensitivity analysis: One SNP at a time was eliminated, and the stability of the remaining SNPs was observed (22).
R version 4.3.2 was utilized to conduct all statistical analyses. Two-sample MR and MR-PRESSO packages were used.
Pharmacovigilance data analysis
In this study, data were acquired using Openvigil 2, an open, clean, and standardized pharmacovigilance data tool. The list of sulfonylureas (A10BB) was determined by the Anatomical Therapeutic Chemical (ATC) classification system. Reports of adverse events (AE) caused by sulfonylureas as the primary suspected drug were collected from the first quarter of 2004 to the second quarter of 2024. All reports were identified using the Medical Dictionary for Regulatory Activities (MedDRA) version 24.1 and focused on analyzing psychiatric disorders.
The potential adverse reaction signals of sulfonylureas were calculated using the algorithms including reporting odds ratio (ROR), proportional reporting ratio (PRR), Bayesian Confidence Propagation Neural Network (BCPNN), and Empirical Bayes Geometric Mean (EBGM). The four algorithms are demonstrated in Supplementary Table S1. A positive AE signal was identified when it met the thresholds for all four methods (ROR: n ≥ 3, lower limit of 95% CI > 1; PRR: χ2 ≥ 4, lower limit of 95% CI > 1; EBGM: EBGM05 (EBGM05 denotes the lower bound of 95% CI) > 2; BCPNN: IC025 (IC025 denotes the lower bound of 95% CI) > 0).