BACKGROUNDMultiple genes might interact to determine the age at onset of bipolar disorder. We investigated gene-gene interactions related to age at onset of bipolar disorder in the Korean population, using genome-wide association study (GWAS) data.METHODSThe study population consisted of 303 patients with bipolar disorder. First, the top 1000 significant single-nucleotide polymorphisms (SNPs) associated with age at onset of bipolar disorder were selected through single SNP analysis by simple linear regression. Subsequently, the QMDR method was used to find gene-gene interactions.RESULTSThe best 10 SNPs from simple regression were located in chromosome 1, 2, 3, 10, 11, 14, 19, and 21. Only five SNPs were found in several genes, such as FOXN3, KIAA1217, OPCML, CAMSAP2, and PTPRS. On QMDR analyses, five pairs of SNPs showed significant interactions with a CVC exceeding 1/5 in a two-locus model. The best interaction was found for the pair of rs60830549 and rs12952733 (CVC = 1/5, P < 1E-07). In three-locus models, four combinations of SNPs showed significant associations with age at onset, with a CVC of >1/5. The best three-locus combination was rs60830549, rs12952733, and rs12952733 (CVC = 2/5, P < 1E-6). The SNPs showing significant interactions were located in the KIAA1217, RBFOX3, SDK2, CYP19A1, NTM, SMYD3, and RBFOX1 genes.CONCLUSIONSOur analysis confirmed genetic interactions influencing the age of onset for bipolar disorder and identified several potential candidate genes. Further exploration of the functions of these promising genes, which may have multiple roles within the neuronal network, is necessary.