Bipolar disorder (BD) is an overarching diagnostic class defined by the presence of at least one prior manic episode (BD I) or both a prior hypomanic episode and a prior depressive episode (BD II). Traditionally, BD II has been conceptualized as a less severe presentation of BD I, however, extant literature to investigate this claim has been mixed.

We apply genomic structural equation modeling (Genomic SEM) to investigate divergent genetic pathways across BD's two major subtypes using the most recent GWAS summary statistics from the PGC. We begin by identifying divergences in genetic correlations across 98 external traits using a Bonferroni-corrected threshold. We also use a theoretically informed follow-up model to examine the extent to which the genetic variance in each subtype is explained by schizophrenia and major depression. Lastly, transcriptome-wide SEM (T-SEM) was used to identify neuronal gene expression patterns associated with BD subtypes.

BD II was characterized by significantly larger genetic overlap across non-psychiatric medical and internalizing traits (e.g. heart disease, neuroticism, insomnia), while stronger associations for BD I were absent. Consistent with these findings, follow-up modeling revealed a substantial major depression component for BD II. T-SEM results revealed 35 unique genes associated with shared risk across BD subtypes.

Divergent patterns of genetic relationships across external traits provide support for the distinction of the bipolar subtypes. However, our results also challenge the illness severity conceptualization of BD given stronger genetic overlap across BD II and a range of clinically relevant traits and disorders.