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 ...