Parental arsenic exposure and tissue-specific DNA methylation in Bangladeshi infants with spina bifida.

An emerging hypothesis linking arsenic toxicity involves altered epigenetic mechanisms, such as DNA methylation. In this study, we examined the relationship between parents' arsenic exposure and DNA methylation in tissues obtained from 28 infants with spina bifida from Bangladesh. We analyzed arsenic in parents' toenails using inductively coupled plasma mass spectrometry (ICP-MS). DNA methylation was measured in infants' dural tissue, buccal swabs, and whole blood using the Illumina Infinium MethylationEPIC BeadChip. We performed epigenome-wide association analyses (EWAS) and tested differentially methylated regions (DMRs). In EWAS, DNA methylation at cg24039697 in dural tissue was positively associated (ß = 0.59, p = 7.6 × 10-9) with father's toenail arsenic concentrations, adjusting for covariates. We did not identify any CpG sites related to father's arsenic exposure in the other tissues, or any CpG sites related to mother's arsenic exposure. Gene ontology analysis identified many biological pathways of interest, including the Wnt signaling pathways. We identified several DMRs across the tissues related to arsenic exposure that included probes mapping to genes that have previously been identified in studies of neural tube defects. This study emphasizes the potential impact of arsenic exposure in fathers, often understudied in epidemiological studies, on DNA methylation in a unique neurological tissue specific to spina bifida.

Maternal arsenic exposure modifies associations between arsenic, folate and arsenic metabolism gene variants, and spina bifida risk: A case‒control study in Bangladesh.

BACKGROUND: Spina bifida is a type of neural tube defect (NTD); NTDs are developmental malformations of the spinal cord that result from failure of neural tube closure during embryogenesis and are likely caused by interactions between genetic and environmental factors. Arsenic induces NTDs in animal models, and studies demonstrate that mice with genetic defects related to folate metabolism are more susceptible to arsenic's effects. We sought to determine whether 25 single-nucleotide polymorphisms (SNPs) in genes involved in folate and arsenic metabolism modified the associations between maternal arsenic exposure and risk of spina bifida (a common NTD) among a hospital-based case-control study population in Bangladesh. METHODS: We used data from 262 mothers and 220 infants who participated in a case‒control study at the National Institutes of Neurosciences & Hospital and Dhaka Shishu Hospital in Dhaka, Bangladesh. Neurosurgeons assessed infants using physical examinations, review of imaging, and we collected histories using questionnaires. We assessed arsenic from mothers' toenails using inductively coupled plasma mass spectrometry (ICP-MS), and we genotyped participants using the Illumina Global Screening Array v1.0. We chose candidate genes and SNPs through a review of the literature. We assessed SNP-environment interactions using interaction terms and stratified models, and we assessed gene-environment interactions using interaction sequence/SNP-set kernel association tests (iSKAT). RESULTS: The median toenail arsenic concentration was 0.42 µg/g (interquartile range [IQR]: 0.27-0.86) among mothers of cases and 0.47 µg/g (IQR: 0.30-0.97) among mothers of controls. We found an two SNPs in the infants' AS3MT gene (rs11191454 and rs7085104) and one SNP in mothers' DNMT1 gene (rs2228611) were associated with increased odds of spina bifida in the setting of high arsenic exposure (rs11191454, OR 3.01, 95% CI: 1.28-7.09; rs7085104, OR 2.33, 95% CI: 1.20-4.and rs2228611, OR 2.11, 95% CI: 1.11-4.01), along with significant SNP-arsenic interactions. iSKAT analyses revealed significant interactions between mothers' toenail concentrations and infants' AS3MT and MTR genes (p = 0.02), and mothers' CBS gene (p = 0.05). CONCLUSIONS: Our results support the hypothesis that arsenic increases spina bifida risk via interactions with folate and arsenic metabolic pathways and suggests that individuals in the population who have certain genetic polymorphisms in genes involved with arsenic and folate metabolism may be more susceptible than others to the arsenic teratogenicity.

Risk of meningomyelocele mediated by the common 22q11.2 deletion.

Meningomyelocele is one of the most severe forms of neural tube defects (NTDs) and the most frequent structural birth defect of the central nervous system. We assembled the Spina Bifida Sequencing Consortium to identify causes. Exome and genome sequencing of 715 parent-offspring trios identified six patients with chromosomal 22q11.2 deletions, suggesting a 23-fold increased risk compared with the general population. Furthermore, analysis of a separate 22q11.2 deletion cohort suggested a 12- to 15-fold increased NTD risk of meningomyelocele. The loss of Crkl, one of several neural tube-expressed genes within the minimal deletion interval, was sufficient to replicate NTDs in mice, where both penetrance and expressivity were exacerbated by maternal folate deficiency. Thus, the common 22q11.2 deletion confers substantial meningomyelocele risk, which is partially alleviated by folate supplementation.

Genome-wide analysis of spina bifida risk variants in a case-control study from Bangladesh.

BACKGROUND: Human studies of genetic risk factors for neural tube defects, severe birth defects associated with long-term health consequences in surviving children, have predominantly been restricted to a subset of candidate genes in specific biological pathways including folate metabolism. METHODS: In this study, we investigated the association of genetic variants spanning the genome with risk of spina bifida (i.e., myelomeningocele and meningocele) in a subset of families enrolled from December 2016 through December 2022 in a case-control study in Bangladesh, a population often underrepresented in genetic studies. Saliva DNA samples were analyzed using the Illumina Global Screening Array. We performed genetic association analyses to compare allele frequencies between 112 case and 121 control children, 272 mothers, and 128 trios. RESULTS: In the transmission disequilibrium test analyses with trios only, we identified three novel exonic spina bifida risk loci, including rs140199800 (SULT1C2, p = 1.9 × 10-7), rs45580033 (ASB2, p = 4.2 × 10-10), and rs75426652 (LHPP, p = 7.2 × 10-14), after adjusting for multiple hypothesis testing. Association analyses comparing cases and controls, as well as models that included their mothers, did not identify genome-wide significant variants. CONCLUSIONS: This study identified three novel single nucleotide polymorphisms involved in biological pathways not previously associated with neural tube defects. The study warrants replication in larger groups to validate findings and to inform targeted prevention strategies.

Altered placental immune cell composition and gene expression with isolated fetal spina bifida.

PROBLEM: Fetal spina bifida (SB) is more common in pregnant people with folate deficiency or anomalies of folate metabolism. It is also known that fetuses with SB have a higher risk of low birthweight, a condition that is typically placental-mediated. We therefore hypothesized that fetal SB would associate with altered expression of key placental folate transporters and an increase in Hofbauer cells (HBCs), which are folate-dependent placental macrophages. METHOD OF STUDY: Folate receptor-α (FRα), proton coupled folate receptor (PCFT), and reduced folate carrier (RFC) protein localization and expression (immunohistochemistry) and HBC phenotypes (HBC abundance and folate receptor-ß [FRß] expression; RNA in situ hybridization) were assessed in placentae from fetuses with SB (cases; n = 12) and in term (n = 10) and gestational age (GA) - and maternal body mass index - matched (n = 12) controls without congenital anomalies. RESULTS: Cases had a higher proportion of placental villous cells that were HBCs (6.9% vs. 2.4%, p = .0001) and higher average HBC FRß expression (3.2 mRNA molecules per HBC vs. 2.3, p = .03) than GA-matched controls. HBCs in cases were largely polarized to a regulatory phenotype (median 92.1% of HBCs). In sex-stratified analyses, only male cases had higher HBC levels and HBC FRß expression than GA-matched controls. There were no differences between groups in the total percent of syncytium and stromal cells that were positive for FRα, PCFT, or RFC protein immunolabeling. CONCLUSIONS: HBC abundance and FRß expression by HBCs are increased in placentae of fetuses with SB, suggesting immune-mediated dysregulation in placental phenotype, and could contribute to SB-associated comorbidities.

Fetal spina bifida associates with dysregulation in nutrient-sensitive placental gene networks: Findings from a matched case-control study.

To improve outcomes in fetuses with spina bifida (SB), better understanding is needed of the molecular drivers of SB and its comorbidities. Pregnant people carrying a fetus with isolated SB (cases; n = 12) or a fetus with no congenital anomalies (controls; n = 21) were recruited at Mount Sinai Hospital, Toronto, Ontario, Canada. Clinical data and placental samples were collected. Placental transcriptome was sequenced (Clariom D microarray) and a nutrient-focused gene expression analysis pipeline was applied to determine whether fetal SB associates with placental dysfunction. Of the 391 differentially expressed genes (DEGs) in cases, 11% (n = 42) had at least one nutrient cofactor, including B vitamins (n = 7 genes), iron/heme (n = 6), and zinc (n = 11). Cases had dysregulation in genes not previously known to associate with SB, and in placental genes that have known links to SB but have not been previously identified in the placenta. Cases also had downregulated nutrient transport and upregulated branching angiogenesis and immune/inflammatory processes. Five nutrient-dependent transcription regulators, collectively predicted to target 46% of DEGs in cases, were identified and were most commonly dependent on B vitamins (n = 3) and zinc (n = 2). Placental gene expression changes were most acute in cases with poor growth. Placentae from fetuses with SB have dysregulation in several gene networks, including those that are sensitive to multiple micronutrients beyond the well-known folic acid. An improved understanding of placental phenotype in fetuses with SB may help identify novel mechanisms associated with comorbidities in fetuses with SB, and reveal new targets to improve fetal outcomes in this population.

ARMC5 controls the degradation of most Pol II subunits, and ARMC5 mutation increases neural tube defect risks in mice and humans.

BACKGROUND: Neural tube defects (NTDs) are caused by genetic and environmental factors. ARMC5 is part of a novel ubiquitin ligase specific for POLR2A, the largest subunit of RNA polymerase II (Pol II). RESULTS: We find that ARMC5 knockout mice have increased incidence of NTDs, such as spina bifida and exencephaly. Surprisingly, the absence of ARMC5 causes the accumulation of not only POLR2A but also most of the other 11 Pol II subunits, indicating that the degradation of the whole Pol II complex is compromised. The enlarged Pol II pool does not lead to generalized Pol II stalling or a generalized decrease in mRNA transcription. In neural progenitor cells, ARMC5 knockout only dysregulates 106 genes, some of which are known to be involved in neural tube development. FOLH1, critical in folate uptake and hence neural tube development, is downregulated in the knockout intestine. We also identify nine deleterious mutations in the ARMC5 gene in 511 patients with myelomeningocele, a severe form of spina bifida. These mutations impair the interaction between ARMC5 and Pol II and reduce Pol II ubiquitination. CONCLUSIONS: Mutations in ARMC5 increase the risk of NTDs in mice and humans. ARMC5 is part of an E3 controlling the degradation of all 12 subunits of Pol II under physiological conditions. The Pol II pool size might have effects on NTD pathogenesis, and some of the effects might be via the downregulation of FOLH1. Additional mechanistic work is needed to establish the causal effect of the findings on NTD pathogenesis.

A Shared Pathogenic Mechanism for Valproic Acid and SHROOM3 Knockout in a Brain Organoid Model of Neural Tube Defects.

Neural tube defects (NTDs), including anencephaly and spina bifida, are common major malformations of fetal development resulting from incomplete closure of the neural tube. These conditions lead to either universal death (anencephaly) or severe lifelong complications (spina bifida). Despite hundreds of genetic mouse models of neural tube defect phenotypes, the genetics of human NTDs are poorly understood. Furthermore, pharmaceuticals, such as antiseizure medications, have been found clinically to increase the risk of NTDs when administered during pregnancy. Therefore, a model that recapitulates human neurodevelopment would be of immense benefit to understand the genetics underlying NTDs and identify teratogenic mechanisms. Using our self-organizing single rosette cortical organoid (SOSR-COs) system, we have developed a high-throughput image analysis pipeline for evaluating the SOSR-CO structure for NTD-like phenotypes. Similar to small molecule inhibition of apical constriction, the antiseizure medication valproic acid (VPA), a known cause of NTDs, increases the apical lumen size and apical cell surface area in a dose-responsive manner. GSK3ß and HDAC inhibitors caused similar lumen expansion; however, RNA sequencing suggests VPA does not inhibit GSK3ß at these concentrations. The knockout of SHROOM3, a well-known NTD-related gene, also caused expansion of the lumen, as well as reduced f-actin polarization. The increased lumen sizes were caused by reduced cell apical constriction, suggesting that impingement of this process is a shared mechanism for VPA treatment and SHROOM3-KO, two well-known causes of NTDs. Our system allows the rapid identification of NTD-like phenotypes for both compounds and genetic variants and should prove useful for understanding specific NTD mechanisms and predicting drug teratogenicity.

A non-coding insertional mutation of Grhl2 causes gene over-expression and multiple structural anomalies including cleft palate, spina bifida and encephalocele.

Orofacial clefts, including cleft lip and palate (CL/P) and neural tube defects (NTDs) are among the most common congenital anomalies, but knowledge of the genetic basis of these conditions remains incomplete. The extent to which genetic risk factors are shared between CL/P, NTDs and related anomalies is also unclear. While identification of causative genes has largely focused on coding and loss of function mutations, it is hypothesized that regulatory mutations account for a portion of the unidentified heritability. We found that excess expression of Grainyhead-like 2 (Grhl2) causes not only spinal NTDs in Axial defects (Axd) mice but also multiple additional defects affecting the cranial region. These include orofacial clefts comprising midline cleft lip and palate and abnormalities of the craniofacial bones and frontal and/or basal encephalocele, in which brain tissue herniates through the cranium or into the nasal cavity. To investigate the causative mutation in the Grhl2Axd strain, whole genome sequencing identified an approximately 4 kb LTR retrotransposon insertion that disrupts the non-coding regulatory region, lying approximately 300 base pairs upstream of the 5' UTR. This insertion also lies within a predicted long non-coding RNA, oriented on the reverse strand, which like Grhl2 is over-expressed in Axd (Grhl2Axd) homozygous mutant embryos. Initial analysis of the GRHL2 upstream region in individuals with NTDs or cleft palate revealed rare or novel variants in a small number of cases. We hypothesize that mutations affecting the regulation of GRHL2 may contribute to craniofacial anomalies and NTDs in humans.

The surface ectoderm exhibits spatially heterogenous tension that correlates with YAP localisation during spinal neural tube closure in mouse embryos.

The single cell layer of surface ectoderm (SE) which overlies the closing neural tube (NT) plays a crucial biomechanical role during mammalian NT closure (NTC), challenging previous assumptions that it is only passive to the force-generating neuroepithelium (NE). Failure of NTC leads to congenital malformations known as NT defects (NTDs), including spina bifida (SB) and anencephaly in the spine and brain respectively. In several mouse NTD models, SB is caused by misexpression of SE-specific genes and is associated with disrupted SE mechanics, including loss of rostrocaudal cell elongation believed to be important for successful closure. In this study, we asked how SE mechanics affect NT morphology, and whether the characteristic rostrocaudal cell elongation at the progressing closure site is a response to tension anisotropy in the SE. We show that blocking SE-specific E-cadherin in ex utero mouse embryo culture influences NT morphology, as well as the F-actin cable. Cell border ablation shows that cell shape is not due to tension anisotropy, but that there are regional differences in SE tension. We also find that YAP nuclear translocation reflects regional tension heterogeneity, and that its expression is sensitive to pharmacological reduction of tension. In conclusion, our results confirm that the SE is a biomechanically important tissue for spinal NT morphogenesis and suggest a possible role of spatial regulation of cellular tension which could regulate downstream gene expression via mechanically-sensitive YAP activity.

Dysregulation of Grainyhead-like 3 expression causes widespread developmental defects.

BACKGROUND: The gene encoding the transcription factor, Grainyhead-like 3 (Grhl3), plays critical roles in mammalian development and homeostasis. Grhl3-null embryos exhibit thoraco-lumbo-sacral spina bifida and soft-tissue syndactyly. Additional studies reveal that these embryos also exhibit an epidermal proliferation/differentiation imbalance. This manifests as skin barrier defects resulting in peri-natal lethality and defective wound repair. Despite these extensive analyses of Grhl3 loss-of-function models, the consequences of gain-of-function of this gene have been difficult to achieve. RESULTS: In this study, we generated a novel mouse model that expresses Grhl3 from a transgene integrated in the Rosa26 locus on an endogenous Grhl3-null background. Expression of the transgene rescues both the neurulation and skin barrier defects of the knockout mice, allowing survival into adulthood. Despite this, the mice are not normal, exhibiting a range of phenotypes attributable to dysregulated Grhl3 expression. In mice homozygous for the transgene, we observe a severe Shaker-Waltzer phenotype associated with hearing impairment. Micro-CT scanning of the inner ear revealed profound structural alterations underlying these phenotypes. In addition, these mice exhibit other developmental anomalies including hair loss, digit defects, and epidermal dysmorphogenesis. CONCLUSION: Taken together, these findings indicate that diverse developmental processes display low tolerance to dysregulation of Grhl3.

Genetic Effects of ITPK1 Polymorphisms on the Risk of Neural Tube Defects: a Population-Based Study.

Inositol is closely related to the occurrence of neural tube defects (NTDs). Inositol 1, 3, 4-trisphosphate 5/6-kinase (ITPK1) gene encoded an essential regulatory enzyme ITPK1, which is involved in inositol metabolism and has a critical role in the development of neural tube and axial mesoderm. It had been reported that some polymorphisms of critical genes in inositol pathways, including ITPK1, were associated with NTDs in Chinese pregnant women; however, the association between fetus ITPK1 polymorphisms and NTDs had not been reported. In a high incidence of NTDs region of China, a case-control study was performed to evaluate the association between fetal ITPK1 polymorphisms and NTDs. The ITPK1 polymorphisms were genotyped by iPLEX® Gold assay. Inositol levels in fetus brain tissues were analyzed. Three genetic polymorphisms of fetus ITPK1's, including rs3818175, rs2295394, and rs4586354, were statistically associated with spina bifida (NTD subtypes). A higher risk of spina bifida was associated with genotype GG of rs3818175, genotype CC of rs4586354, and genotype TT of rs2295394 (OR = 2.66, 95% CI [1.17-6.05], P = 0.017; OR = 2.22, 95% CI [1.02-4.80], P = 0.041; and OR = 2.33, 95% CI [1.00-5.48], P = 0.047), when compared with the other wild-type genotypes CC, TT, and CC, respectively. Decreased brain inositol level was found in NTDs fetuses, compared to normal controls. Inositol levels were found to significantly decrease with rs2295394 (CC genotype), rs4586354 (TT genotype), and rs3818175 (GC genotype) (P < 0.05). The polymorphisms of fetus ITPK1 were associated with the incidence of NTDs and might be a genetic risk factor for spina bifida.

CIC missense variants contribute to susceptibility for spina bifida.

Neural tube defects (NTDs) are congenital malformations resulting from abnormal embryonic development of the brain, spine, or spinal column. The genetic etiology of human NTDs remains poorly understood despite intensive investigation. CIC, homolog of the Capicua transcription repressor, has been reported to interact with ataxin-1 (ATXN1) and participate in the pathogenesis of spinocerebellar ataxia type 1. Our previous study demonstrated that CIC loss of function (LoF) variants contributed to the cerebral folate deficiency syndrome by downregulating folate receptor 1 (FOLR1) expression. Given the importance of folate transport in neural tube formation, we hypothesized that CIC variants could contribute to increased risk for NTDs by depressing embryonic folate concentrations. In this study, we examined CIC variants from whole-genome sequencing (WGS) data of 140 isolated spina bifida cases and identified eight missense variants of CIC gene. We tested the pathogenicity of the observed variants through multiple in vitro experiments. We determined that CIC variants decreased the FOLR1 protein level and planar cell polarity (PCP) pathway signaling in a human cell line (HeLa). In a murine cell line (NIH3T3), CIC loss of function variants downregulated PCP signaling. Taken together, this study provides evidence supporting CIC as a risk gene for human NTD.

Eph and Ephrin Variants in Malaysian Neural Tube Defect Families.

Neural tube defects (NTDs) are common birth defects with a complex genetic etiology. Mouse genetic models have indicated a number of candidate genes, of which functional mutations in some have been found in human NTDs, usually in a heterozygous state. This study focuses on Ephs-ephrins as candidate genes of interest owing to growing evidence of the role of this gene family during neural tube closure in mouse models. Eph-ephrin genes were analyzed in 31 Malaysian individuals comprising seven individuals with sporadic spina bifida, 13 parents, one twin-sibling and 10 unrelated controls. Whole exome sequencing analysis and bioinformatic analysis were performed to identify variants in 22 known Eph-ephrin genes. We reported that three out of seven spina bifida probands and three out of thirteen family members carried a variant in either EPHA2 (rs147977279), EPHB6 (rs780569137) or EFNB1 (rs772228172). Analysis of public databases shows that these variants are rare. In exome datasets of the probands and parents of the probands with Eph-ephrin variants, the genotypes of spina bifida-related genes were compared to investigate the probability of the gene-gene interaction in relation to environmental risk factors. We report the presence of Eph-ephrin gene variants that are prevalent in a small cohort of spina bifida patients in Malaysian families.

A SOX3 duplication and lumbosacral spina bifida in three generations.

Chromosomal aneuploidies, microduplications and microdeletions are the most common confirmed genetic causes of spina bifida. Microduplications of Xq27 containing the SOX3 gene have been reported in 11 cases, confirming the existence of an X-chromosomal locus for spina bifida. A three generation kindred reported here with a SOX3 duplication has been identified in one of 17 kindreds with recurrences in the 29 years of the South Carolina Neural Tube Defect Prevention Program. Other recurrences during this time period included siblings with an APAF1 mutation, siblings with a CASP9 mutation, siblings with a microdeletion of 13q, and two sets of siblings with Meckel syndrome who did not have genetic/genomic studies performed.

Whole exome sequencing identifies potential candidate genes for spina bifida derived from mouse models.

Spina bifida (SB) is the second most common nonlethal congenital malformation. The existence of monogenic SB mouse models and human monogenic syndromes with SB features indicate that human SB may be caused by monogenic genes. We hypothesized that whole exome sequencing (WES) allows identification of potential candidate genes by (i) generating a list of 136 candidate genes for SB, and (ii) by unbiased exome-wide analysis. We generated a list of 136 potential candidate genes from three categories and evaluated WES data of 50 unrelated SB cases for likely deleterious variants in 136 potential candidate genes, and for potential SB candidate genes exome-wide. We identified 6 likely deleterious variants in 6 of the 136 potential SB candidate genes in 6 of the 50 SB cases, whereof 4 genes were derived from mouse models, 1 gene was derived from human nonsyndromic SB, and 1 gene was derived from candidate genes known to cause human syndromic SB. In addition, by unbiased exome-wide analysis, we identified 12 genes as potential candidates for SB. Identification of these 18 potential candidate genes in larger SB cohorts will help decide which ones can be considered as novel monogenic causes of human SB.

Folate deficiency disturbs PEG10 methylation modifications in human spina bifida.

BACKGROUND: Paternally expressed gene 10 (PEG10) is believed to be a key imprinted gene involved in placenta formation. However, its role in human folate-related spina bifida (SB) remains unclear. METHODS: The methylation status of the germline differentially methylated region (gDMR) in the PEG10/sarcoglycan epsilon (SGCE) imprinted cluster was compared between SB patients and control samples. Moreover, the influence of ectopic PEG10 expression on apoptosis was assessed to explore the underlying mechanisms related to folate deficiency-induced aberrant gDMR methylation in SB. RESULTS: The case group exhibited a significant increase in the methylation level of the gDMR and a marked reduction in the mRNA and protein expression of PEG10 compared with the control group. A prominent negative correlation was found between the folate level in brain tissue and gDMR methylation status (r = -0.62, P = 0.001). A cell model treated with a demethylating agent showed a significant elevation of PEG10 transcription level, as well as other imprinted genes in this cluster. In addition, the inhibition of PEG10 was found to be accompanied by aberrant activation of apoptosis in SB. CONCLUSIONS: Our findings suggest that disturbed gDMR methylation of the PEG10/SGCE cluster due to folate deficiency is involved in SB through aberrant activation of apoptosis. IMPACT: Disturbed genomic imprinting has been verified to be involved in neural tube defects (NTDs). However, little is known about the effect of ectopic expression of imprinted gene PEG10 on human NTDs. Aberrant methylation status of the germline differentially methylated region (gDMR) of PEG10/SGCE cluster due to folate deficiency has been found to result in the inhibition of PEG10 and has a marked association with an increased occurrence of spina bifida. Inhibited expression of PEG10 partly is found to be related to the abnormal activation of apoptosis in spina bifida.

Systems biology analysis of human genomes points to key pathways conferring spina bifida risk.

Spina bifida (SB) is a debilitating birth defect caused by multiple gene and environment interactions. Though SB shows non-Mendelian inheritance, genetic factors contribute to an estimated 70% of cases. Nevertheless, identifying human mutations conferring SB risk is challenging due to its relative rarity, genetic heterogeneity, incomplete penetrance, and environmental influences that hamper genome-wide association studies approaches to untargeted discovery. Thus, SB genetic studies may suffer from population substructure and/or selection bias introduced by typical candidate gene searches. We report a population based, ancestry-matched whole-genome sequence analysis of SB genetic predisposition using a systems biology strategy to interrogate 298 case-control subject genomes (149 pairs). Genes that were enriched in likely gene disrupting (LGD), rare protein-coding variants were subjected to machine learning analysis to identify genes in which LGD variants occur with a different frequency in cases versus controls and so discriminate between these groups. Those genes with high discriminatory potential for SB significantly enriched pathways pertaining to carbon metabolism, inflammation, innate immunity, cytoskeletal regulation, and essential transcriptional regulation consistent with their having impact on the pathogenesis of human SB. Additionally, an interrogation of conserved noncoding sequences identified robust variant enrichment in regulatory regions of several transcription factors critical to embryonic development. This genome-wide perspective offers an effective approach to the interrogation of coding and noncoding sequence variant contributions to rare complex genetic disorders.

Altered sacral neural crest development in Pax3 spina bifida mutants underlies deficits of bladder innervation and function.

Mouse models of Spina bifida (SB) have been instrumental for identifying genes, developmental processes, and environmental factors that influence neurulation and neural tube closure. Beyond the prominent neural tube defects, other aspects of the nervous system can be affected in SB with significant changes in essential bodily functions such as urination. SB patients frequently experience bladder dysfunction and SB fetuses exhibit reduced density of bladder nerves and smooth muscle although the developmental origins of these deficits have not been determined. The Pax3 Splotch-delayed (Pax3Sp-d) mouse model of SB is one of a very few mouse SB models that survives to late stages of gestation. Through analysis of Pax3Sp-d mutants we sought to define how altered bladder innervation in SB might arise by tracing sacral neural crest (NC) development, pelvic ganglia neuronal differentiation, and assessing bladder nerve fiber density. In Pax3Sp-d/Sp-d fetal mice we observed delayed migration of Sox10+ NC-derived progenitors (NCPs), deficient pelvic ganglia neurogenesis, and reduced density of bladder wall innervation. We further combined NC-specific deletion of Pax3 with the constitutive Pax3Sp-d allele in an effort to generate viable Pax3 mutants to examine later stages of bladder innervation and postnatal bladder function. Neural crest specific deletion of a Pax3 flox allele, using a Sox10-cre driver, in combination with a constitutive Pax3Sp-d mutation produced postnatal viable offspring that exhibited altered bladder function as well as reduced bladder wall innervation and altered connectivity between accessory ganglia at the bladder neck. Combined, the results show that Pax3 plays critical roles within sacral NC that are essential for initiation of neurogenesis and differentiation of autonomic neurons within pelvic ganglia.

Identification of new candidate genes for spina bifida through exome sequencing.

PURPOSE: Neural tube defects are a group of birth defects caused by failure of neural tube closure during development. The etiology of NTD, requiring a complex interaction between environmental and genetic factors, is not well understood. METHODS: We performed whole-exome sequencing (WES) in six trios, with a single affected proband with spina bifida, to identify rare/novel variants as potential causes of the NTD. RESULTS: Our analysis identified four de novo and ten X-linked recessive variants in four of the six probands, all of them in genes previously never implicated in NTD. Among the 14 variants, we ruled out six of them, based on different criteria and pursued the evaluation of eight potential candidates in the following genes: RXRγ, DTX1, COL15A1, ARHGAP36, TKTL1, AMOT, GPR50, and NKRF. The de novo variants where located in the RXRγ, DTX1, and COL15A1 genes while ARHGAP36, TKTL1, AMOT, GPR50, and NKRF carry X-linked recessive variants. This analysis also revealed that four patients presented multiple variants, while we were unable to identify any significant variant in two patients. CONCLUSIONS: Our preliminary conclusion support a major role for the de novo variants with respect to the X-linked recessive variants where the X-linked could represent a contribution to the phenotype in an oligogenic model.