Recently, two genome-wide association studies (GWASs) of schizophrenia (SCZ) in Han Chinese identified several susceptibility loci. Replication efforts aiming to validate the GWAS findings were made and focused on the top hits. We conducted a more extensive follow-up study in an independent sample of 1471 cases and 1528 matched controls to verify 26 genetic variants by including nine top single-nucleotide polymorphisms (SNPs) that reached genome-wide significance and 17 promising SNPs nominated in the initial discovery phase. rs8073471 in an intron of tubulin-folding cofactor D (TBCD) obtained nominal significance (P<0.01) in single SNP analysis. Logistic regression identified significant interaction between rs3744165 (5'-untranslated region variant of exon 2 of zinc finger protein 750 (ZNF750), and in an intron of TBCD) and rs8073471 (Deviance test P-value=2.77 × 10(-34)). Both SNPs are located at 17q25, an interesting region that has been implicated in SCZ. By using the Genotype-Tissue Expression (GTEx) data set, we implemented an expression quantitative trait loci epistasis analysis to explore the association between the genotype combinations of the two SNPs and gene expression levels in 13 areas of human central nervous system. We observed that rs3744165 × rs8073471 interaction modulated the expression profile of TEAD3 (P=1.87 × 10(-8)), SH3TC2 (P=2.00 × 10(-8)), KCNK9 (P=5.20 × 10(-7)) and PPDPF (P=1.13 × 10(-6)) in postmortem cortex tissue; EFNA1 (P=7.26 × 10(-9)), RNU4ATAC (P=2.32 × 10(-8)) and NUPL2 (P=6.79 × 10(-8)) in cerebellum tissue. To the best of our knowledge, our study is the first one that links TBCD and ZNF750 mutations to SCZ susceptibility and to the transcript levels in human brain tissues. Further efforts are needed to understand the role of those variants in the pathogenesis of SCZ.