Druggable genome-wide Mendelian randomization for identifying the role of integrated stress response in therapeutic targets of bipolar disorder

For bipolar disorder (BD), the inconsistency of treatment guidelines and the long phases of pharmacological adjustment remain major challenges. BD is known to be comorbid with many medical and psychiatric conditions and they may share inflammatory and stress-related aetiologies, which could give rise to this association. The integrated stress response (ISR) responds to various stress conditions that lead to alterations in cellular homeostasis. However, as a causative mechanism underlying cognitive deficits and neurodegeneration in a broad range of brain disorders, the impact of ISR on BD is understudied. Mendelian randomization has been widely used to repurpose licensed drugs and discover novel therapeutic targets. Thus, we aimed to identify novel therapeutic targets for BD and analyze their pathophysiological mechanisms, using the summary data-based Mendelian Randomization (SMR) and Bayesian colocalization (COLOC) methods to integrate the summary- level data of the GWAS on BD and the expression quantitative trait locus (eQTL) study in blood. We utilized the GWAS data including 41,917 BD cases and 371,549 controls from the Psychiatric Genomics Consortium and the eQTL data from 31,684 participants of predominantly European ancestry from the eQTLGen consortium. The SMR analysis identified the EIF2B5 gene that was associated with BD due to no linkage but pleiotropy or causality. The COLOC analysis strongly suggested that EIF2B5 and the trait of BD were affected by shared causal variants, and thus were colocalized. Utilizing data in EpiGraphDB we find other putative causal BD genes (EIF2AK4 and GSK3B) to prioritize potential alternative drug targets.