Lithium response in bipolar disorder is associated with focal adhesion and PI3K-Akt networks: a multi-omics replication study

Lithium is the gold standard treatment for bipolar disorder (BD). However, its mechanism of action is incompletely understood, and prediction of treatment outcomes is limited. In our previous multi-omics study of the Pharmacogenomics of Bipolar Disorder (PGBD) sample combining transcriptomic and genomic data, we found that focal adhesion, the extracellular matrix (ECM), and PI3K-Akt signaling networks were associated with response to lithium. In this study, we replicated the results of our previous study using network propagation methods in a genome-wide association study of an independent sample of 2039 patients from the International Consortium on Lithium Genetics (ConLiGen) study. We identified functional enrichment in focal adhesion and PI3K-Akt pathways, but we did not find an association with the ECM pathway. Our results suggest that deficits in the neuronal growth cone and PI3K-Akt signaling, but not in ECM proteins, may influence response to lithium in BD.


INTRODUCTION
Bipolar disorder (BD) is a chronic psychiatric illness that presents with episodes of mania, depression, and sometimes psychosis.Globally, it is the sixth leading cause of medical disability among people from 15 to 44 years old.Patients with BD are at a higher risk of suicide than those with any other psychiatric or medical illness.Some studies report that roughly 50% of patients will attempt suicide, and up to 20% of untreated patients will complete suicide [1], while treatment by lithium reduces that risk significantly [2,3].Unfortunately, misdiagnosis is common and often delays an accurate treatment.Up to 70% of patients are initially misdiagnosed, usually with major depressive disorder.On average, there is a delay of 8 years before the correct diagnosis of BD is made [4].During this time, patients continue to suffer, may be treated with medications that make their illness course worse, and are at risk of suicide.
Lithium is the gold standard treatment for BD [5].Its mechanism of action is still not completely understood [6].Many studies have investigated the neurotrophic effect of lithium.One theory posits that chronic administration of lithium inhibits glycogen synthase kinase 3 (GSK3β), a serine/threonine kinase.This leads to anti-apoptotic effects and improved cell structural stability [7][8][9][10].GSK3β has also been shown to exhibit interactions with many pathways, including phosphorylation of several components of the PI3K/AKT/mTOR signaling network, as well as regulation of transcription for proteins bound to microtubules [11].Another theory involves the phosphoinositol (PI) cycle.In the PI cycle, lithium inhibits inositol monophosphatase, which ultimately downregulates protein kinase C isozymes such as myristoylated alanine-rich C-kinase substrate (MARCKS).MARCKS is an actin-binding protein found in neuronal processes that is implicated in cytoskeletal restructuring.Its downregulation stabilizes the neuronal membrane and results in neurotrophic effects [7,12].A more recent theory proposes that lithium alters the phosphorylation state of collapsin response mediator protein-2 (CRMP2).CRMP2 regulates cytoskeletal organization, particularly in dendritic spines [13,14].Finally, a study using polygenic score modeling has indicated that the cholinergic and glutamatergic pathways may potentially serve as targets for lithium [15].It is possible that lithium exerts its effects through multiple or all of these pathways.A single definitive model remains elusive, but interactions with neuronal cytoskeleton are possibly involved.
Interestingly, there is a range of responses to treatment with lithium.Previous studies have reported that 20-30% of patients with BD are excellent responders, whereas over 40% fail to demonstrate any significant clinical improvement.These patient populations have been shown to differ from each other both phenotypically and genetically [16].A differential response to lithium has been previously demonstrated between induced pluripotent stem cell (iPSC) neurons derived from lithium responders and non-responders.The hyperexcitability of in vitro neurons derived from BD patients was reversed by lithium treatment, but only in those from patients who were lithium responders [17].This finding is also supported by family studies, which found that the relatives of lithium responders were significantly more likely to be lithium responders as well [18,19].These studies imply that patients with BD could be subcategorized based on biological differences which induce a divergent lithium response.There is a great need to better understand these differences in order to identify possible predictors of treatment response.However, dozens of previous candidate-gene association studies, genome-wide association studies (GWAS), and polygenic risk score analyses of lithium response in BD have failed to identify genetic variants with major effects.Given this pressing need to find pharmacogenetic predictors of response, more advanced methods in integrative genomic analysis are necessary [16].
GWAS inherently face several limitations when used in isolation, including the challenge of genetic heterogeneity.In many disease processes with genetic associations, patients may carry diverse combinations of causal variants that impact multiple genes, creating a net effect across a particular pathway.GWAS of BD primarily detect variants of very small effect size consistent with a polygenic mode of transmission.Since each single nucleotide polymorphism (SNP) contributes only a tiny amount to the overall predisposition to BD, enormous sample sizes are required, and it can be difficult to surmise mechanisms of disease.Network approaches seek to address this biological reality by integrating GWAS results with known protein-protein interactions and other molecular networks.New causal genes may be identified by boosting their interactions with products of known causal genes [20,21].
We have recently reported a combined analysis of transcriptomic and GWAS data from the Pharmacogenomics of Bipolar Disorder (PGBD) study [22] of treatment response to lithium.After using network propagation to reprioritize candidate genes from GWAS data, we found significant overlap between both transcriptomic and GWAS results.The joint analysis yielded a 500 gene network significantly enriched in the following Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways: focal adhesion, ECMreceptor interaction, and PI3K-Akt signaling [23].All three pathways play a role in axon growth and neuronal development [24].Consistent with these results, post-mortem studies have found that in BD, neuronal populations may exhibit a decrease in number, size, and/or amount of dendritic spines [13,25].Given that lithium may have downstream effects on these pathways, it is possible that genetic defects in focal adhesion pathways may provide both a mechanism for susceptibility to BD as well as a target for lithium treatment.
In this study, we aimed to replicate the results of our previous multi-omics study on a larger dataset of over 2000 patients from the International Consortium on Lithium Genetics (ConLiGen) [26].We reprioritized GWAS results using network methods to determine overlap with focal adhesion, ECM-receptor interaction, and PI3K-Akt signaling pathways.

METHODS
Summary statistics were downloaded from the NHGRI-EBI GWAS Catalog [27] on 12/12/2022 for study GCST012487 [26].The data resulted from a GWAS of lithium response in 2563 patients at 22 sites participating in the International Consortium on Lithium Genetics (ConLiGen).We utilized the summary statistics from a combined sample of 2039 European ancestry individuals.In the ConLiGen study, data from over 6 million single nucleotide polymorphisms (SNPs) were tested for association with categorical and continuous retrospective ratings of lithium response using the Alda scale [28,29].The Alda scale includes two scores: score A is a 0-10 retrospective rating of lifetime response, while score B captures factors reducing the confidence in score A such as lack of a documented lithium level, etc.In the ConLiGen study, under the continuous phenotype, participants were rated with the Alda A score, and individuals with a B score greater than 4 were excluded.We used the continuous rather than the dichotomous phenotype as a measure of treatment response because genome-wide significant association was detected with the continuous phenotype in the original GWAS.Quality control and statistical analysis methods are described in the original paper.

SNP, gene, and gene-set analysis
We imported the ConLiGen summary statistics into FUMA (Functional Mapping and Annotation of Genome-Wide Association Studies-https:// fuma.ctglab.nl)[30], a web-based platform for annotating, prioritizing, visualizing and interpreting GWAS results.We utilized the SNP2GENE function to map SNPs to genes and conduct SNP, gene-based, and geneset analysis.We used all default settings, except for setting the maximum lead SNP p value to 1 × 10e−5.

Network analysis
We input the ConLiGen summary statistics into NAGA (Network Assisted Genomic Analysis), an online network propagation tool for pathway boosting and interpretation of genome-wide association studies [21].NAGA provided a reprioritized ranked list of 19,781 genes as output.We then entered the top 500 genes with the highest final heat scores into STRING, an online database that generates mapped networks based on protein-protein interactions [31].STRING additionally analyzes for overrepresentation of user-inputted gene lists in established pathways, using the hypergeometric test [32].Using this function, we tested our a priori hypotheses to identify functional enrichment of the NAGA-generated top 500 gene list in the KEGG hsa04510 focal adhesion pathway, KEGG hsa04512 ECM-receptor interaction, and KEGG hsa04151 PI3K-Akt signaling pathway [33].p values were corrected for multiple testing by STRING using the Benjamini-Hochberg procedure [34].
Overlap between the NAGA-generated top 500 gene list and the KEGG pathways was visualized using Cytoscape [35].A hypergeometric test was conducted to test for overrepresentation of the NAGA-generated 500 gene network in the 500 gene network generated in our previous study [23].

SNP, gene, and gene-set analysis
As reported in the original ConLiGen study, the only SNPs that were significant at a genome-wide significance level of 5e-08 were in linkage disequilibrium with the SNP rs74795342 on chromosome 21 (Supplementary Fig. 1).Using FUMA in our gene-wise analysis, no significant genes were found at a significance level of p < 0.05/18314 = 2.730e−6 (Supplementary Fig. 2).No gene-sets were found to be significant either, using p < 0.05 after Bonferroni correction.The most highly associated genes and gene-sets are listed in Supplementary Tables 1 and 2.

Network analysis
We first tested the three a priori pathways that were significant in our previous study, which had examined an independent sample [23].Using the STRING analysis function, the top 500 reprioritized gene list generated by NAGA was found to be significantly enriched in both the KEGG hsa04510 focal adhesion pathway (p = 1.74e−06) and KEGG hsa04151 PI3K-Akt signaling pathway (p = 1.90e−07) (Table 1).Given the goal of replication and the small number of statistical tests, this was considered as a significant replication of our previous results in an independent sample for the focal adhesion and PI3K-Akt pathways.However, the KEGG hsa04512 ECM-receptor interaction pathway was not found to be significantly enriched (Table 1).The overlapping genes in all three networks can be seen in Figs.1-3.
A hypergeometric test found significant overlap (p = 5.699e−07) between the 500 gene network generated by NAGA and the 500 gene network generated by network propagation analysis in our previous study [23].There were 33 genes that were common to both networks.The top 25 reprioritized genes produced by NAGA are listed in Table 2.All top 500 reprioritized NAGA genes are listed in Supplementary Table 4.
After testing the three a priori hypotheses based on previous results, we tested the top 500 NAGA gene list for enrichment in all pathways in STRING.The top 10 KEGG pathways found to be most strongly enriched are found in Supplementary Table 3.These include cancer and growth pathways (such as Pathways in Cancer, Estrogen Signaling Pathway, Ras Signaling Pathway) as well as the dopaminergic synapse pathway.
We additionally used the STRING analysis function to test for functional enrichment of the top 100, 200, 300, and 400 reprioritized gene lists generated by NAGA in all three a priori KEGG pathways.The results agreed with the primary analysis, since all gene lists were significantly enriched in the KEGG hsa04510 focal adhesion pathway and KEGG hsa04151 PI3K-Akt signaling pathways at a level of p < 0.05.Only the top 100 reprioritized gene list was found to be significantly enriched in the KEGG hsa04512 ECM-receptor interaction pathway (p = 0.0050) inconsistent with a robust result.(Supplementary Table 5).

DISCUSSION
In this study, we attempted to our previous results which were from an independent sample [23].We used network methods via NAGA to reprioritize GWAS results from the ConLiGen study on lithium response and used STRING to test three a priori network hypotheses: KEGG focal adhesion, ECM-receptor interaction and PI3K-Akt signaling.Two of these three networks, KEGG focal adhesion and PI3K-Akt signaling, were enriched in our top 500 reprioritized genes.However, we did not find significant enrichment for the ECM-receptor interaction pathway in the 500 gene network.Besides this pathway, we were otherwise able to replicate the results of our previous paper in a larger, independent sample of patients with BD.We found highly significant overlap between the top 500 gene network generated by NAGA in this study and the 500 gene network generated in the previous study, providing further evidence for replication.
Focal adhesions are points of contact between cells and proteins in the ECM.The formation of cell-ECM adhesion structures is initiated by cell surface integrins and driven by local actin polymerization.These structures function to not only mediate cell attachment to ECM, but also mediate transmembrane signaling.Integrin-ECM ligand binding can induce a number of downstream changes affecting cell shape, growth, and proliferation [36].In neurons, specifically, the actin cytoskeleton of growth cones interacts with the ECM to guide axon development and extension [24,37].
We had originally hypothesized that genetic deficits in focal adhesion, ECM, and PI3K-Akt pathways may impair axonal growth in neurons and determine response to lithium.Though one integrin protein was included in our top 500 genes, in general ECM proteins did not overlap with the top 500 gene list (Supplementary Table 4) (Fig. 2), and the pathway was not significant.This result is inconsistent with our previous study.However, it may suggest the possibility that the deficits influencing lithium response may be inherent to the growth cone rather than components of the ECM.This is supported by a number of studies, which have shown that lithium prevents collapse and induces growth of growth cones [38][39][40].
Fig. 2 Overlap between KEGG ECM-receptor interaction and top 500 genes.KEGG hsa04512 pathway for ECM-receptor interaction adapted to illustrate gene overlap.Genes in yellow overlap with the 500 gene NAGA network.
Previously, neurons derived from induced pluripotent stem cells of patients with BD have shown to exhibit hyperexcitability in vitro.This hyperexcitable phenotype was rescued by lithium only in neurons derived from lithium good responders [17].Elevated neuroactivity in BD may induce vulnerability in neurons through impairment of focal adhesion pathways.Chronic elevation of neuroactivity has been shown to dramatically reduce surface expression of integrin β1 in animal models, leading to axonal and dendritic degeneration and eventually cell death [41].
Unsurprisingly, neurons in patients with BD have been shown to be present with smaller size, fewer numbers, and more limited branching.We had previously proposed that in lithium responders, this deficit is caused by deficits in focal adhesion and is rescued by lithium treatment.Furthermore, we proposed that in patients who are not lithium responders, focal adhesion is not dysregulated, and lithium is unable to address the relevant impairments [42][43][44].Our results in this study are consistent with this hypothesis.
After testing our three a priori hypotheses, we conducted exploratory analyses using network methods.We listed the top 10 most significant KEGG pathways that were associated by STRING with the NAGA generated gene list in Supplementary Table 3.These pathways are mostly cancer pathways associated with cell growth and proliferation or pathways of addiction and other dopamine-related processes.Dopamine neurotransmission has previously been associated with response to lithium treatment in BD [45].Genes in associated cancer pathways show some overlap with focal adhesion as well, which suggests the possibility of shared mechanisms (Fig. 1).
Limitations of our study include the relatively small sample size (N = 2039) and the generalizability of the dataset, given that all participants were of European descent.Additionally, data was collected retrospectively.As a result, outcomes may be less accurate in determining response phenotypes [46] which can blur our findings due to false negatives.
This study also demonstrates the utility of network propagation methods, which can add power to GWAS with limited sample sizes.These methods are beneficial in identifying which genes and gene-sets are of interest to a disease process, but future research is still indicated for confirmation [20,21].
In summary, we replicated our previous results reinforcing that genetic deficits in focal adhesion and PI3K-Akt signaling are associated with lithium response in BD patients.We hypothesize, as before, that malformed axonal growth cones result in shorter and less branched axons and susceptibility to BD in a subpopulation of patients who are lithium responders.This is also consistent with the idea that response to lithium results from a disease mechanism distinct from that of lithium non-responders.Furthermore, we propose that lithium rescues disrupted neuronal growth and axon extension processes by addressing deficits in focal adhesion.A better understanding of the pathophysiology of BD and lithium treatment may lead to the future development of drugs similar to lithium, as well as possible clinical predictors for treatment response.
p values corrected for multiple testing using the Benjamini-Hochberg procedure.*Significant at p < 0.05.

Fig. 1
Fig. 1 Overlap between KEGG focal adhesion and top 500 genes.KEGG hsa04510 pathway for focal adhesion adapted to illustrate gene overlap.Genes in yellow overlap with the 500 gene NAGA network.

Fig. 3
Fig. 3 Overlap between KEGG PI3k-Akt and top 500 genes.KEGG hsa04151 pathway for PI3k-Akt signaling adapted to illustrate gene overlap.Genes in yellow overlap with the 500 gene NAGA network.

Table 1 .
Functional enrichment of NAGA top 500 gene list in focal adhesion, ECM, and PI3K-Akt pathways.
a Data does not exist as gene was not evaluated in FUMA.