Genes in the noncanonical WNT signaling pathway controlling planar cell polarity have been linked to the neural tube defect myelomeningocele. We hypothesized that some genes in the WNT signaling network have a higher mutational burden in myelomeningocele subjects than in reference subjects in gnomAD. Exome sequencing data from 511 myelomeningocele subjects was obtained in-house and data from 29,940 ethnically matched subjects was provided by version 2 of the publicly available Genome Aggregation Database. To compare mutational burden, we collapsed rare deleterious variants across each of 523 human WNT signaling genes in case and reference populations. Ten WNT signaling genes were disrupted with a higher mutational burden among Mexican American myelomeningocele subjects compared to reference subjects (Fishers exact test, P � 0.05) and seven different genes were disrupted among individuals of European ancestry compared to reference subjects. Gene ontology enrichment analyses indicate that genes disrupted only in the Mexican American population play a role in planar cell polarity whereas genes identified in both populations are important for the regulation of canonical WNT signaling. In summary, evidence for WNT signaling genes that may contribute to myelomeningocele in humans is presented and discussed.

This is an open access article under the terms of the Creative Commons Attribution-NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Myelomeningocele (MMC) affects one in 1000 newborns annually worldwide and each surviving child faces tremendous lifetime medical and caregiving burdens. Both genetic and environmental factors contribute to disease risk but the mechanism is unclear. This study examined 506 MMC subjects for ultra-rare deleterious variants (URDVs, absent in gnomAD v2.1.1 controls that have Combined Annotation Dependent Depletion score ≥ 20) in candidate genes either known to cause abnormal neural tube closure in animals or previously associated with human MMC in the current study cohort. Approximately 70% of the study subjects carried one to nine URDVs among 302 candidate genes. Half of the study subjects carried heterozygous URDVs in multiple genes involved in the structure and/or function of cilium, cytoskeleton, extracellular matrix, WNT signaling, and/or cell migration. Another 20% of the study subjects carried heterozygous URDVs in candidate genes associated with gene transcription regulation, folate metabolism, or glucose metabolism. Presence of URDVs in the candidate genes involving these biological function groups may elevate the risk of developing myelomeningocele in the study cohort.