Identification of novel candidate risk genes for myelomeningocele within the glucose homeostasis/oxidative stress and folate/one-carbon metabolism networks.

BACKGROUND: Neural tube defects (NTDs) are the second most common complex birth defect, yet, our understanding of the genetic contribution to their development remains incomplete. Two environmental factors associated with NTDs are Folate and One Carbon Metabolism (FOCM) and Glucose Homeostasis and Oxidative Stress (GHOS). Utilizing next-generation sequencing of a large patient cohort, we identify novel candidate genes in these two networks to provide insights into NTD mechanisms. METHODS: Exome sequencing (ES) was performed in 511 patients, born with myelomeningocele, divided between European American and Mexican American ethnicities. Healthy control data from the Genome Aggregation database were ethnically matched and used as controls. Rare, high fidelity, nonsynonymous predicted damaging missense, nonsense, or canonical splice site variants in independently generated candidate gene lists for FOCM and GHOS were identified. We used a gene-based collapsing approach to quantify mutational burden in case and controls, with the control cohort estimated using cumulative allele frequencies assuming Hardy-Weinberg equilibrium. RESULTS: We identified 45 of 837 genes in the FOCM network and 22 of 568 genes in the GHOS network as possible NTD risk genes with p < 0.05. No nominally significant risk genes were shared between ethnicities. Using a novel approach to mutational burden we identify 55 novel NTD risk associations. CONCLUSIONS: We provide a means of utilizing large publicly available sequencing datasets as controls for sequencing projects examining rare disease. This approach confirmed existing risk genes for myelomeningocele and identified possible novel risk genes. Lastly, it suggests possible distinct genetic etiologies for this malformation between different ethnicities.

Burden of rare deleterious variants in WNT signaling genes among 511 myelomeningocele patients.

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.