Our model-selection analysis was designed to characterise the range of phenotypic effects for loci that showed significant association; however, the statistical evidence and effect size for each contributing disorder can vary. Although the best-fit model for the chromosome 3p21 region included all five disorders, interpretation of these results is complicated by evidence from a PGC GWAS mega-analysis of major depressive disorder.17 In the discovery phase of that analysis, which consisted of 9240 major depressive disorder cases and 9519 controls many of whom overlap with samples reported here, the smallest association p value for this region (rs2535629) was 0·00013. However, no association was noted in a replication dataset of 6783 cases and 50 695 controls (p=0·70) for that disorder, and the combined discovery and replication phase p value was 0·0031. Thus, any association between this region and major depressive disorder is unclear.