A likelihood procedure described previously was used for association analysis26. To correct for relatedness and potential population stratification, genomic control was used27. Within our study, samples were combined using the Mantel–Haenszel model28. Our results were merged with those of the International Schizophrenia Consortium and the Molecular Genetics of Schizophrenia study by computing weighted averages of Z scores (see Methods).