To screen the single nucleotide polymorphisms (SNPs) in the coding regions of VANGL and FZD family members related to the plane cell polarity (PCP) signaling pathway in neural tube defects (NTDs) patients, so as to provide theoretical and experimental basis for the prevention and treatment of NTDs by intervening PCP signal transduction. 112 NTDs patients were collected as the case group and 112 craniocerebral trauma patients as control. Afterwards, blood genomic DNA was extracted and sequenced. The distribution of SNP alleles and genotypes between case and control groups was analyzed. Finally, the NTD rat model was constructed, and the effect of SNPs on the expression level of VANGL and FZD genes was verified by qRT-PCR. GC genotype was newly found at VANGL1 c.346G>A, as well as AT genotype in FZD6 c.97A>G. The distribution of VANGL1 c.346g>A allele and genotype was statistically different between the case and control groups (p < 0.05). The newly found genotype GC increased the risk of NTDs (OR=9.918, 95% CI: 1.234%-79.709%). The results of qRT-PCR showed that the expression level of FZD6 in E11 NTD fetuses were significantly increased (p < 0.05), but there was no obvious difference in the expression of VANGL1. We found a new variant of VANGL1 c.346G>A, whose GC genotype might play an important role in the pathogenesis of NTDs. The SNPs of VANGL1 had no significant effect on its expression level, indicating that it may induce NTDs through other ways. FZD6 was significantly overexpressed in NTDs fetuses.