Using restriction fragment length polymorphism and pyrosequencing methods, we genotyped two TNFA gene promoter SNPs (ÀG308A, ÀG238A) and analyzed the haplotype structure in 24 Canadian families of primarily Celtic origin. Our results demonstrate that after correction for multiple testing based on simulations of 10 000 replicates of unlinked/unassociated data, there is evidence for association (P ¼ 0.026) of a specific haplotype (À308A, À238G) with schizophrenia and schizophrenia spectrum disorders with a family-based trimmed haplotype linkage disequilibrium test (Trimhap). Stratifying the 22 families with genome scan data by TNFA promoter haplotypes followed by reanalysis of linkage to schizophrenia throughout the genome, we identified few loci that exhibit a considerable increase in LOD/HLOD scores. A locus on chromosome 1q44 (D1S1609) demonstrated a significant increase (P ¼ 0.025) in LOD score from 0.15 to 3.01 with a broad definition of the schizophrenia phenotype and a dominant mode of inheritance. This result replicates a previously reported positive result of linkage of schizophrenia spectrum disorders to this area of the genome. We also illustrated that simulation studies are pivotal in evaluating the significance of results obtained with newer statistical methods, when multiple, but not independent, tests are performed, and when sample stratification is utilized to reduce the impact of heterogeneity or assess the interaction between loci.