Schizophrenia polygenic risk scores, clinical variables and genetic pathways as predictors of phenotypic traits of bipolar I disorder.

AIM: We investigated the predictive value of polygenic risk scores (PRS) derived from the schizophrenia GWAS (Trubetskoy et al., 2022) (SCZ3) for phenotypic traits of bipolar disorder type-I (BP-I) in 1878 BP-I cases and 2751 controls from Romania and UK. METHODS: We used PRSice-v2.3.3 and PRS-CS for computing SCZ3-PRS for testing the predictive power of SCZ3-PRS alone and in combination with clinical variables for several BP-I subphenotypes and for pathway analysis. Non-linear predictive models were also used. RESULTS: SCZ3-PRS significantly predicted psychosis, incongruent and congruent psychosis, general age-of-onset (AO) of BP-I, AO-depression, AO-Mania, rapid cycling in univariate regressions. A negative correlation between the number of depressive episodes and psychosis, mainly incongruent and an inverse relationship between increased SCZ3-SNP loading and BP-I-rapid cycling were observed. In random forest models comparing the predictive power of SCZ3-PRS alone and in combination with nine clinical variables, the best predictions were provided by combinations of SCZ3-PRS-CS and clinical variables closely followed by models containing only clinical variables. SCZ3-PRS performed worst. Twenty-two significant pathways underlying psychosis were identified. LIMITATIONS: The combined RO-UK sample had a certain degree of heterogeneity of the BP-I severity: only the RO sample and partially the UK sample included hospitalized BP-I cases. The hospitalization is an indicator of illness severity. Not all UK subjects had complete subphenotype information. CONCLUSION: Our study shows that the SCZ3-PRS have a modest clinical value for predicting phenotypic traits of BP-I. For clinical use their best performance is in combination with clinical variables.

Fine-mapping genomic loci refines bipolar disorder risk genes.

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

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.

Exploring the genetics of lithium response in bipolar disorders.

Background: Lithium (Li) remains the treatment of choice for bipolar disorders (BP). Its mood-stabilizing effects help reduce the long-term burden of mania, depression and suicide risk in patients with BP. It also has been shown to have beneficial effects on disease-associated conditions, including sleep and cardiovascular disorders. However, the individual responses to Li treatment vary within and between diagnostic subtypes of BP (e.g. BP-I and BP-II) according to the clinical presentation. Moreover, long-term Li treatment has been linked to adverse side-effects that are a cause of concern and non-adherence, including the risk of developing chronic medical conditions such as thyroid and renal disease. In recent years, studies by the Consortium on Lithium Genetics (ConLiGen) have uncovered a number of genetic factors that contribute to the variability in Li treatment response in patients with BP. Here, we leveraged the ConLiGen cohort (N=2,064) to investigate the genetic basis of Li effects in BP. For this, we studied how Li response and linked genes associate with the psychiatric symptoms and polygenic load for medical comorbidities, placing particular emphasis on identifying differences between BP-I and BP-II. Results: We found that clinical response to Li treatment, measured with the Alda scale, was associated with a diminished burden of mania, depression, substance and alcohol abuse, psychosis and suicidal ideation in patients with BP-I and, in patients with BP-II, of depression only. Our genetic analyses showed that a stronger clinical response to Li was modestly related to lower polygenic load for diabetes and hypertension in BP-I but not BP-II. Moreover, our results suggested that a number of genes that have been previously linked to Li response variability in BP differentially relate to the psychiatric symptomatology, particularly to the numbers of manic and depressive episodes, and to the polygenic load for comorbid conditions, including diabetes, hypertension and hypothyroidism. Conclusions: Taken together, our findings suggest that the effects of Li on symptomatology and comorbidity in BP are partially modulated by common genetic factors, with differential effects between BP-I and BP-II.

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 2,039 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.

Lack of guidelines and translational knowledge is hindering the implementation of psychiatric genetic counseling and testing within Europe - A multi-professional survey study.

Genetic research has identified a large number of genetic variants, both rare and common, underlying neurodevelopmental disorders (NDD) and major psychiatric disorders. Currently, these findings are being translated into clinical practice. However, there is a lack of knowledge and guidelines for psychiatric genetic testing (PsychGT) and genetic counseling (PsychGC). The European Union-funded COST action EnGagE (CA17130) network was started to investigate the current implementation status of PsychGT and PsychGC across 35 participating European countries. Here, we present the results of a pan-European online survey in which we gathered the opinions, knowledge, and practices of a self-selected sample of professionals involved/interested in the field. We received answers from 181 respondents. The three main occupational categories were genetic counselor (21.0%), clinical geneticist (24.9%), and researcher (25.4%). Of all 181 respondents, 106 provide GC for any psychiatric disorder or NDD, corresponding to 58.6% of the whole group ranging from 43.2% in Central Eastern Europe to 66.1% in Western Europe. Overall, 65.2% of the respondents reported that genetic testing is offered to individuals with NDD, and 26.5% indicated the same for individuals with major psychiatric disorders. Only 22.1% of the respondents indicated that they have guidelines for PsychGT. Pharmacogenetic testing actionable for psychiatric disorders was offered by 15%. Interestingly, when genetic tests are fully covered by national health insurance, more genetic testing is provided for individuals with NDD but not those with major psychiatric disorders. Our qualitative analyses of responses highlight the lack of guidelines and knowledge on utilizing and using genetic tests and education and training as the major obstacles to implementation. Indeed, the existence of psychiatric genetic training courses was confirmed by only 11.6% of respondents. The question on the relevance of up-to-date education and training in psychiatric genetics on everyday related practice was highly relevant. We provide evidence that PsychGC and PsychGT are already in use across European countries, but there is a lack of guidelines and education. Harmonization of practice and development of guidelines for genetic counseling, testing, and training professionals would improve equality and access to quality care for individuals with psychiatric disorders within Europe.

Association of polygenic score and the involvement of cholinergic and glutamatergic pathways with lithium treatment response in patients with bipolar disorder.

Lithium is regarded as the first-line treatment for bipolar disorder (BD), a severe and disabling mental health disorder that affects about 1% of the population worldwide. Nevertheless, lithium is not consistently effective, with only 30% of patients showing a favorable response to treatment. To provide personalized treatment options for bipolar patients, it is essential to identify prediction biomarkers such as polygenic scores. In this study, we developed a polygenic score for lithium treatment response (Li+PGS) in patients with BD. To gain further insights into lithium's possible molecular mechanism of action, we performed a genome-wide gene-based analysis. Using polygenic score modeling, via methods incorporating Bayesian regression and continuous shrinkage priors, Li+PGS was developed in the International Consortium of Lithium Genetics cohort (ConLi+Gen: N = 2367) and replicated in the combined PsyCourse (N = 89) and BipoLife (N = 102) studies. The associations of Li+PGS and lithium treatment response - defined in a continuous ALDA scale and a categorical outcome (good response vs. poor response) were tested using regression models, each adjusted for the covariates: age, sex, and the first four genetic principal components. Statistical significance was determined at P < 0.05. Li+PGS was positively associated with lithium treatment response in the ConLi+Gen cohort, in both the categorical (P = 9.8 × 10-12, R2 = 1.9%) and continuous (P = 6.4 × 10-9, R2 = 2.6%) outcomes. Compared to bipolar patients in the 1st decile of the risk distribution, individuals in the 10th decile had 3.47-fold (95%CI: 2.22-5.47) higher odds of responding favorably to lithium. The results were replicated in the independent cohorts for the categorical treatment outcome (P = 3.9 × 10-4, R2 = 0.9%), but not for the continuous outcome (P = 0.13). Gene-based analyses revealed 36 candidate genes that are enriched in biological pathways controlled by glutamate and acetylcholine. Li+PGS may be useful in the development of pharmacogenomic testing strategies by enabling a classification of bipolar patients according to their response to treatment.

Immunogenetics of lithium response and psychiatric phenotypes in patients with bipolar disorder.

The link between bipolar disorder (BP) and immune dysfunction remains controversial. While epidemiological studies have long suggested an association, recent research has found only limited evidence of such a relationship. To clarify this, we investigated the contributions of immune-relevant genetic factors to the response to lithium (Li) treatment and the clinical presentation of BP. First, we assessed the association of a large collection of immune-related genes (4,925) with Li response, defined by the Retrospective Assessment of the Lithium Response Phenotype Scale (Alda scale), and clinical characteristics in patients with BP from the International Consortium on Lithium Genetics (ConLi+Gen, N = 2,374). Second, we calculated here previously published polygenic scores (PGSs) for immune-related traits and evaluated their associations with Li response and clinical features. We found several genes associated with Li response at p < 1×10- 4 values, including HAS3, CNTNAP5 and NFIB. Network and functional enrichment analyses uncovered an overrepresentation of pathways involved in cell adhesion and intercellular communication, which appear to converge on the well-known Li-induced inhibition of GSK-3ß. We also found various genes associated with BP's age-at-onset, number of mood episodes, and presence of psychosis, substance abuse and/or suicidal ideation at the exploratory threshold. These included RTN4, XKR4, NRXN1, NRG1/3 and GRK5. Additionally, PGS analyses suggested serum FAS, ECP, TRANCE and cytokine ligands, amongst others, might represent potential circulating biomarkers of Li response and clinical presentation. Taken together, our results support the notion of a relatively weak association between immunity and clinically relevant features of BP at the genetic level.

Association of Polygenic Score and the involvement of Cholinergic and Glutamatergic Pathways with Lithium Treatment Response in Patients with Bipolar Disorder.

Lithium is regarded as the first-line treatment for bipolar disorder (BD), a severe and disabling mental disorder that affects about 1% of the population worldwide. Nevertheless, lithium is not consistently effective, with only 30% of patients showing a favorable response to treatment. To provide personalized treatment options for bipolar patients, it is essential to identify prediction biomarkers such as polygenic scores. In this study, we developed a polygenic score for lithium treatment response (Li+PGS) in patients with BD. To gain further insights into lithium's possible molecular mechanism of action, we performed a genome-wide gene-based analysis. Using polygenic score modeling, via methods incorporating Bayesian regression and continuous shrinkage priors, Li+PGS was developed in the International Consortium of Lithium Genetics cohort (ConLi+Gen: N=2,367) and replicated in the combined PsyCourse (N=89) and BipoLife (N=102) studies. The associations of Li+PGS and lithium treatment response - defined in a continuous ALDA scale and a categorical outcome (good response vs. poor response) were tested using regression models, each adjusted for the covariates: age, sex, and the first four genetic principal components. Statistical significance was determined at P<����������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������.

Investigating the phenotypic and genetic associations between personality traits and suicidal behavior across major mental health diagnoses.

Personality traits influence risk for suicidal behavior. We examined phenotype- and genotype-level associations between the Big Five personality traits and suicidal ideation and attempt in major depressive, bipolar and schizoaffective disorder, and schizophrenia patients (N = 3012) using fixed- and random-effects inverse variance-weighted meta-analyses. Suicidal ideations were more likely to be reported by patients with higher neuroticism and lower extraversion phenotypic scores, but showed no significant association with polygenic load for these personality traits. Our findings provide new insights into the association between personality and suicidal behavior across mental illnesses and suggest that the genetic component of personality traits is unlikely to have strong causal effects on suicidal behavior.

Using polygenic scores and clinical data for bipolar disorder patient stratification and lithium response prediction: machine learning approach.

BACKGROUND: Response to lithium in patients with bipolar disorder is associated with clinical and transdiagnostic genetic factors. The predictive combination of these variables might help clinicians better predict which patients will respond to lithium treatment. AIMS: To use a combination of transdiagnostic genetic and clinical factors to predict lithium response in patients with bipolar disorder. METHOD: This study utilised genetic and clinical data (n = 1034) collected as part of the International Consortium on Lithium Genetics (ConLi+Gen) project. Polygenic risk scores (PRS) were computed for schizophrenia and major depressive disorder, and then combined with clinical variables using a cross-validated machine-learning regression approach. Unimodal, multimodal and genetically stratified models were trained and validated using ridge, elastic net and random forest regression on 692 patients with bipolar disorder from ten study sites using leave-site-out cross-validation. All models were then tested on an independent test set of 342 patients. The best performing models were then tested in a classification framework. RESULTS: The best performing linear model explained 5.1% (P = 0.0001) of variance in lithium response and was composed of clinical variables, PRS variables and interaction terms between them. The best performing non-linear model used only clinical variables and explained 8.1% (P = 0.0001) of variance in lithium response. A priori genomic stratification improved non-linear model performance to 13.7% (P = 0.0001) and improved the binary classification of lithium response. This model stratified patients based on their meta-polygenic loadings for major depressive disorder and schizophrenia and was then trained using clinical data. CONCLUSIONS: Using PRS to first stratify patients genetically and then train machine-learning models with clinical predictors led to large improvements in lithium response prediction. When used with other PRS and biological markers in the future this approach may help inform which patients are most likely to respond to lithium treatment.

Genetic risk for psychiatric illness is associated with the number of hospitalizations of bipolar disorder patients.

OBJECTIVES: Bipolar disorder (BD) has a highly heterogeneous clinical course that is characterized by relapses and increased health care utilization in a significant fraction of patients. A thorough understanding of factors influencing illness course is essential for predicting disorder severity and developing targeted therapies. METHODS: We performed polygenic score analyses in four cohorts (N = 954) to test whether the genetic risk for BD, schizophrenia, or major depression is associated with a severe course of BD. We analyzed BD patients with a minimum illness duration of five years. The severity of the disease course was assessed by using the number of hospitalizations in a mental health facility and a composite measure of longitudinal illness severity (OPCRIT item 90). RESULTS: Our analyses showed that higher polygenic scores for BD (ß = 0.11, SE = 0.03, p = 1.17 × 10-3) and schizophrenia (ß = 0.09, SE = 0.03, p = 4.24 × 10-3), but not for major depression, were associated with more hospitalizations. None of the investigated polygenic scores was associated with the composite measure of longitudinal illness severity (OPCRIT item 90). LIMITATIONS: We could not account for non-genetic influences on disease course. Our clinical sample contained more severe cases. CONCLUSIONS: This study demonstrates that the genetic risk burden for psychiatric illness is associated with increased health care utilization, a proxy for disease severity, in BD patients. The findings are in line with previous observations made for patients diagnosed with schizophrenia or major depression. Therefore, in the future psychiatric disorder polygenic scores might become helpful for stratifying patients with high risk of a chronic manifestation and predicting disease course.

Combining schizophrenia and depression polygenic risk scores improves the genetic prediction of lithium response in bipolar disorder patients.

Lithium is the gold standard therapy for Bipolar Disorder (BD) but its effectiveness differs widely between individuals. The molecular mechanisms underlying treatment response heterogeneity are not well understood, and personalized treatment in BD remains elusive. Genetic analyses of the lithium treatment response phenotype may generate novel molecular insights into lithium's therapeutic mechanisms and lead to testable hypotheses to improve BD management and outcomes. We used fixed effect meta-analysis techniques to develop meta-analytic polygenic risk scores (MET-PRS) from combinations of highly correlated psychiatric traits, namely schizophrenia (SCZ), major depression (MD) and bipolar disorder (BD). We compared the effects of cross-disorder MET-PRS and single genetic trait PRS on lithium response. For the PRS analyses, we included clinical data on lithium treatment response and genetic information for n = 2283 BD cases from the International Consortium on Lithium Genetics (ConLi+Gen; www.ConLiGen.org ). Higher SCZ and MD PRSs were associated with poorer lithium treatment response whereas BD-PRS had no association with treatment outcome. The combined MET2-PRS comprising of SCZ and MD variants (MET2-PRS) and a model using SCZ and MD-PRS sequentially improved response prediction, compared to single-disorder PRS or to a combined score using all three traits (MET3-PRS). Patients in the highest decile for MET2-PRS loading had 2.5 times higher odds of being classified as poor responders than patients with the lowest decile MET2-PRS scores. An exploratory functional pathway analysis of top MET2-PRS variants was conducted. Findings may inform the development of future testing strategies for personalized lithium prescribing in BD.

HLA-DRB1 and HLA-DQB1 genetic diversity modulates response to lithium in bipolar affective disorders.

Bipolar affective disorder (BD) is a severe psychiatric illness, for which lithium (Li) is the gold standard for acute and maintenance therapies. The therapeutic response to Li in BD is heterogeneous and reliable biomarkers allowing patients stratification are still needed. A GWAS performed by the International Consortium on Lithium Genetics (ConLiGen) has recently identified genetic markers associated with treatment responses to Li in the human leukocyte antigens (HLA) region. To better understand the molecular mechanisms underlying this association, we have genetically imputed the classical alleles of the HLA region in the European patients of the ConLiGen cohort. We found our best signal for amino-acid variants belonging to the HLA-DRB1*11:01 classical allele, associated with a better response to Li (p < 1 × 10-3; FDR < 0.09 in the recessive model). Alanine or Leucine at position 74 of the HLA-DRB1 heavy chain was associated with a good response while Arginine or Glutamic acid with a poor response. As these variants have been implicated in common inflammatory/autoimmune processes, our findings strongly suggest that HLA-mediated low inflammatory background may contribute to the efficient response to Li in BD patients, while an inflammatory status overriding Li anti-inflammatory properties would favor a weak response.

Prediction of lithium response using genomic data.

Predicting lithium response prior to treatment could both expedite therapy and avoid exposure to side effects. Since lithium responsiveness may be heritable, its predictability based on genomic data is of interest. We thus evaluate the degree to which lithium response can be predicted with a machine learning (ML) approach using genomic data. Using the largest existing genomic dataset in the lithium response literature (n = 2210 across 14 international sites; 29% responders), we evaluated the degree to which lithium response could be predicted based on 47,465 genotyped single nucleotide polymorphisms using a supervised ML approach. Under appropriate cross-validation procedures, lithium response could be predicted to above-chance levels in two constituent sites (Halifax, Cohen's kappa 0.15, 95% confidence interval, CI [0.07, 0.24]; and Würzburg, kappa 0.2 [0.1, 0.3]). Variants with shared importance in these models showed over-representation of postsynaptic membrane related genes. Lithium response was not predictable in the pooled dataset (kappa 0.02 [- 0.01, 0.04]), although non-trivial performance was achieved within a restricted dataset including only those patients followed prospectively (kappa 0.09 [0.04, 0.14]). Genomic classification of lithium response remains a promising but difficult task. Classification performance could potentially be improved by further harmonization of data collection procedures.

Characterisation of age and polarity at onset in bipolar disorder.

BACKGROUND: Studying phenotypic and genetic characteristics of age at onset (AAO) and polarity at onset (PAO) in bipolar disorder can provide new insights into disease pathology and facilitate the development of screening tools. AIMS: To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics. METHOD: Genome-wide association studies (GWASs) and polygenic score (PGS) analyses of AAO (n = 12 977) and PAO (n = 6773) were conducted in patients with bipolar disorder from 34 cohorts and a replication sample (n = 2237). The association of onset with disease characteristics was investigated in two of these cohorts. RESULTS: Earlier AAO was associated with a higher probability of psychotic symptoms, suicidality, lower educational attainment, not living together and fewer episodes. Depressive onset correlated with suicidality and manic onset correlated with delusions and manic episodes. Systematic differences in AAO between cohorts and continents of origin were observed. This was also reflected in single-nucleotide variant-based heritability estimates, with higher heritabilities for stricter onset definitions. Increased PGS for autism spectrum disorder (ß = -0.34 years, s.e. = 0.08), major depression (ß = -0.34 years, s.e. = 0.08), schizophrenia (ß = -0.39 years, s.e. = 0.08), and educational attainment (ß = -0.31 years, s.e. = 0.08) were associated with an earlier AAO. The AAO GWAS identified one significant locus, but this finding did not replicate. Neither GWAS nor PGS analyses yielded significant associations with PAO. CONCLUSIONS: AAO and PAO are associated with indicators of bipolar disorder severity. Individuals with an earlier onset show an increased polygenic liability for a broad spectrum of psychiatric traits. Systematic differences in AAO across cohorts, continents and phenotype definitions introduce significant heterogeneity, affecting analyses.

Association of Attention-Deficit/Hyperactivity Disorder and Depression Polygenic Scores with Lithium Response: A Consortium for Lithium Genetics Study.

Response to lithium varies widely between individuals with bipolar disorder (BD). Polygenic risk scores (PRSs) can uncover pharmacogenomics effects and may help predict drug response. Patients (N = 2,510) with BD were assessed for long-term lithium response in the Consortium on Lithium Genetics using the Retrospective Criteria of Long-Term Treatment Response in Research Subjects with Bipolar Disorder score. PRSs for attention-deficit/hyperactivity disorder (ADHD), major depressive disorder (MDD), and schizophrenia (SCZ) were computed using lassosum and in a model including all three PRSs and other covariates, and the PRS of ADHD (ß = -0.14; 95% confidence interval [CI]: -0.24 to -0.03; p value = 0.010) and MDD (ß = -0.16; 95% CI: -0.27 to -0.04; p value = 0.005) predicted worse quantitative lithium response. A higher SCZ PRS was associated with higher rates of medication nonadherence (OR = 1.61; 95% CI: 1.34-1.93; p value = 2e-7). This study indicates that genetic risk for ADHD and depression may influence lithium treatment response. Interestingly, a higher SCZ PRS was associated with poor adherence, which can negatively impact treatment response. Incorporating genetic risk of ADHD, depression, and SCZ in combination with clinical risk may lead to better clinical care for patients with BD.

Association of polygenic score for major depression with response to lithium in patients with bipolar disorder.

Lithium is a first-line medication for bipolar disorder (BD), but only one in three patients respond optimally to the drug. Since evidence shows a strong clinical and genetic overlap between depression and bipolar disorder, we investigated whether a polygenic susceptibility to major depression is associated with response to lithium treatment in patients with BD. Weighted polygenic scores (PGSs) were computed for major depression (MD) at different GWAS p value thresholds using genetic data obtained from 2586 bipolar patients who received lithium treatment and took part in the Consortium on Lithium Genetics (ConLi+Gen) study. Summary statistics from genome-wide association studies in MD (135,458 cases and 344,901 controls) from the Psychiatric Genomics Consortium (PGC) were used for PGS weighting. Response to lithium treatment was defined by continuous scores and categorical outcome (responders versus non-responders) using measurements on the Alda scale. Associations between PGSs of MD and lithium treatment response were assessed using a linear and binary logistic regression modeling for the continuous and categorical outcomes, respectively. The analysis was performed for the entire cohort, and for European and Asian sub-samples. The PGSs for MD were significantly associated with lithium treatment response in multi-ethnic, European or Asian populations, at various p value thresholds. Bipolar patients with a low polygenic load for MD were more likely to respond well to lithium, compared to those patients with high polygenic load [lowest vs highest PGS quartiles, multi-ethnic sample: OR = 1.54 (95% CI: 1.18-2.01) and European sample: OR = 1.75 (95% CI: 1.30-2.36)]. While our analysis in the Asian sample found equivalent effect size in the same direction: OR = 1.71 (95% CI: 0.61-4.90), this was not statistically significant. Using PGS decile comparison, we found a similar trend of association between a high genetic loading for MD and lower response to lithium. Our findings underscore the genetic contribution to lithium response in BD and support the emerging concept of a lithium-responsive biotype in BD.