Sixty-eight SNPs with p<1×10−5, representing 16 independent genomic regions, were identified, though none reached the genome-wide significance threshold of p<5×10−8 (Table 1,Figure 1; Table S1). The most significant association was found in rs4988462 on 3p11 (p=3.72×10−7, OR=1.18). This SNP lies within an intron of POU1F1, though the entire 279 kb region of association in linkage disequilibrium (LD) with rs4988462 contains 16 additional SNPs with p<1×10−5 and includes CHMP2B and POU1F1, as well as the microRNA MIR4795. Regional association and forest plots from the top five independent GWAS signals are provided in Figures S4–S8. Eleven of the 68 SNPs were also identified in the original OCD GWAS with p<1×10−5; none of these SNPs were identified in the TS GWAS at p<1×10−5 (15, 16) (Table S1).
