Identification of genetic variants that influence bipolar I disorder (BPD-I) through genome-wide association (GWA) studies is limited in Asian populations. The current study aimed to identify novel common variants for BPD-I in an ethnically homogeneous Taiwanese sample using a multi-stage GWA study design.At the discovery stage, 200 BPD-I patients and 200 controls that combined to form 16 pools were genotyped with 1 million markers. Utilizing a newly developed rank-based method, top-ranked markers were selected. After validation with individual genotyping, a fine-mapping association study was conducted to identify associated loci using 240 patients and 240 controls. At the last stage, independent samples were collected (351 cases and 341 controls) for replication.Among the top-ranked markers from the discovery stage, eight genes and 15 individual SNPs were evaluated in the fine-mapping stage. At this stage, rs7619173, which is not in a gene coding region, showed the most significant association (P = 2 ∗ 10(-5)) with BPD-I. Four genes had empirical P-values<0.05, including KCNH7 (P = 0.0047), MYST4 (P = 0.0047), NRXN3 (P = 0.0095), and SEMA3D (P = 0.037). For markers genotyped in replication samples, rs7619173 exhibited a significant association (P(combined) = 2 ∗ 10(-4)) after multiple testing correction, while markers rs11001178 (MYST4) and rs2217887 (NRXN3) showed weak associations (P(combined) = 0.02) with BPD-I.A multi-stage GWA design has the potential to uncover the underlying pathogenesis of a complex trait. Findings in the present study highlight three loci that warrant further investigation for bipolar.