The predictive power of polygenic scores (PGSs) for lithium treatment response in bipolar disorder remains limited.

To enhance the prediction of lithium responsiveness by developing a multi-trait PGS (mt-PGS) combining genetic information from multiple phenotypes implicated in lithium response and/or bipolar disorder aetiology.

We analysed data collected from bipolar disorder patients who had received lithium treatment for at least 6 months and participated in the International Consortium on Lithium Genetics study (N = 2367). The Alda scale was used to assess lithium responsiveness, and treatment outcome was defined as continuous total Alda score (0–10) and categorical outcome (favourable ≥7 versus unfavourable response). PGSs were calculated for 60 phenotypes grouped into 5 clinical–biological clusters: psychiatric and behavioural (#23 phenotypes), cardiometabolic (#17), autoimmune/inflammatory (#5), neurocognitive (#8) and renal function (#7). We applied cross-validated machine learning regression approaches in both outcomes within each cluster, with the selected features from each cluster subsequently combined to construct the final mt-PGS models. Model performance was assessed using explained variance (R2) for the continuous outcome, and both McFadden’s pseudo-R2 and standard classification model parameters for the categorical outcome.

mt-PGS explained between 5.70% (continuous outcome) and 9.11% (categorical outcome) of the interindividual variability in lithium responsiveness. Classification accuracy (area under the curve) for the categorical outcome was 68.24% (95% CI: 64.98–71.66), with a Brier score of 0.257. Of the 5 clusters, the PGSs for psychiatric and behavioural phenotypes were most strongly associated with lithium responsiveness, accounting for 3.67–6.52% of its variability.

By integrating PGSs for multiple relevant phenotypes, predictive accuracy for lithium response improved substantially compared with single-trait methods. Future research incorporating larger, more diverse populations and combining genetic scores with clinical data holds promise for further enhancing prediction and advancing clinical implementation.

Accurate diagnosis of bipolar disorder (BPD) is difficult in clinical practice, with an average delay between symptom onset and diagnosis of about 7 years. A depressive episode often precedes the first manic episode, making it difficult to distinguish BPD from unipolar major depressive disorder (MDD).

We use genome-wide association analyses (GWAS) to identify differential genetic factors and to develop predictors based on polygenic risk scores (PRS) that may aid early differential diagnosis.

Based on individual genotypes from case–control cohorts of BPD and MDD shared through the Psychiatric Genomics Consortium, we compile case–case–control cohorts, applying a careful quality control procedure. In a resulting cohort of 51 149 individuals (15 532 BPD patients, 12 920 MDD patients and 22 697 controls), we perform a variety of GWAS and PRS analyses.

Although our GWAS is not well powered to identify genome-wide significant loci, we find significant chip heritability and demonstrate the ability of the resulting PRS to distinguish BPD from MDD, including BPD cases with depressive onset (BPD-D). We replicate our PRS findings in an independent Danish cohort (iPSYCH 2015, N = 25 966). We observe strong genetic correlation between our case–case GWAS and that of case–control BPD.

We find that MDD and BPD, including BPD-D are genetically distinct. Our findings support that controls, MDD and BPD patients primarily lie on a continuum of genetic risk. Future studies with larger and richer samples will likely yield a better understanding of these findings and enable the development of better genetic predictors distinguishing BPD and, importantly, BPD-D from MDD.

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.

To use a combination of transdiagnostic genetic and clinical factors to predict lithium response in patients with bipolar disorder.

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.

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.

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.

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.

To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics.

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.

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.

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.