ResultsA total of 13 DNVs were identified in 8 anencephaly trios by WES, including two loss of function (LoF) variants detected in pLI > 0.9 genes (SPHKAP, c.2629_2633del, and NCOR1, p.Y1907X). Damaging DNVs were identified in 61.5% (8/13) of the anencephalic cases. Independent validation was conducted in an additional 502 NTD cases. Gene inactivation using targeted morpholino antisense oligomers and rescue assays were conducted in zebrafish, and transfection expression in HEK293T cells. Four DNVs in four cases were identified and predicted to alter protein function, including p.R328Q in WD repeat domain phosphoinositide-interacting 1 (WIPI1). Three variants, p.G313R, p.T418M, and p.L406P, in the WIPI1 gene were identified from the independent replication cohort consisting of 502 cases. Functional analysis suggested that the wipi1 p.L406P and p.R328Q variants most likely displayed loss-of-function effects during embryonic development.