The interaction of endorepellin and neurexin triggers neuroepithelial autophagy and maintains neural tube development.

Heparan sulfate proteoglycan 2 (HSPG2) gene encodes the matrix protein Perlecan, and genetic inactivation of this gene creates mice that are embryonic lethal with severe neural tube defects (NTDs). We discovered rare genetic variants of HSPG2 in 10% cases compared to only 4% in controls among a cohort of 369 NTDs. Endorepellin, a peptide cleaved from the domain V of Perlecan, is known to promote angiogenesis and autophagy in endothelial cells. The roles of enderepellin in neurodevelopment remain unclear so far. Our study revealed that endorepellin can migrate to the neuroepithelial cells and then be recognized and bind with the neuroepithelia receptor neurexin in vivo. Through the endocytic pathway, the interaction of endorepellin and neurexin physiologically triggers autophagy and appropriately modulates the differentiation of neural stem cells into neurons as a blocker, which is necessary for normal neural tube closure. We created knock-in (KI) mouse models with human-derived HSPG2 variants, using sperm-like stem cells that had been genetically edited by CRISPR/Cas9. We realized that any HSPG2 variants that affected the function of endorepellin were considered pathogenic causal variants for human NTDs given that the severe NTD phenotypes exhibited by these KI embryos occurred in a significantly higher response frequency compared to wildtype embryos. Our study provides a paradigm for effectively confirming pathogenic mutations in other genetic diseases. Furthermore, we demonstrated that using autophagy inhibitors at a cellular level can repress neuronal differentiation. Therefore, autophagy agonists may prevent NTDs resulting from failed autophagy maintenance and neuronal over-differentiation caused by deleterious endorepellin variants.

Folate regulation of planar cell polarity pathway and F-actin through folate receptor alpha.

Folate deficiency contribute to neural tube defects (NTDs) which could be rescued by folate supplementation. However, the underlying mechanisms are still not fully understood. Besides, there is considerable controversy concerning the forms of folate used for supplementation. To address this controversy, we prepared culture medium with different forms of folate, folic acid (FA), and 5-methyltetrahydrofolate (5mTHF), at concentrations of 5 µM, 500 nM, 50 nM, and folate free, respectively. Mouse embryonic fibroblasts (MEFs) were treated with different folates continuously for three passages, and cell proliferation and F-actin were monitored. We determined that compared to 5mTHF, FA showed stronger effects on promoting cell proliferation and F-actin formation. We also found that FOLR1 protein level was positively regulated by folate concentration and the non-canonical Wnt/planar cell polarity (PCP) pathway signaling was significantly enriched among different folate conditions in RNA-sequencing analyses. We demonstrated for the first time that FOLR1 could promote the transcription of Vangl2, one of PCP core genes. The transcription of Vangl2 was down-regulated under folate-deficient condition, which resulted in a decrease in PCP activity and F-actin formation. In summary, we identified a distinct advantage of FA in cell proliferation and F-actin formation over 5mTHF, as well as demonstrating that FOLR1 could promote transcription of Vangl2 and provide a new mechanism by which folate deficiency can contribute to the etiology of NTDs.

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.

Functional characterization of age-dependent p16 epimutation reveals biological drivers and therapeutic targets for colorectal cancer.

BACKGROUND: Methylation of the p16 promoter resulting in epigenetic gene silencing-known as p16 epimutation-is frequently found in human colorectal cancer and is also common in normal-appearing colonic mucosa of aging individuals. Thus, to improve clinical care of colorectal cancer (CRC) patients, we explored the role of age-related p16 epimutation in intestinal tumorigenesis. METHODS: We established a mouse model that replicates two common genetic and epigenetic events observed in human CRCs: Apc mutation and p16 epimutation. We conducted long-term survival and histological analysis of tumor development and progression. Colonic epithelial cells and tumors were collected from mice and analyzed by RNA sequencing (RNA-seq), quantitative PCR, and flow cytometry. We performed single-cell RNA sequencing (scRNA-seq) to characterize tumor-infiltrating immune cells throughout tumor progression. We tested whether anti-PD-L1 immunotherapy affects overall survival of tumor-bearing mice and whether inhibition of both epigenetic regulation and immune checkpoint is more efficacious. RESULTS: Mice carrying combined Apc mutation and p16 epimutation had significantly shortened survival and increased tumor growth compared to those with Apc mutation only. Intriguingly, colon tumors with p16 epimutation exhibited an activated interferon pathway, increased expression of programmed death-ligand 1 (Pdl1), and enhanced infiltration of immune cells. scRNA-seq further revealed the presence of Foxp3+ Tregs and γδT17 cells, which contribute to an immunosuppressive tumor microenvironment (TME). Furthermore, we showed that a combined therapy using an inhibitor of DNA methylation and a PD-L1 immune checkpoint inhibitor is more effective for improving survival in tumor-bearing mice than blockade of either pathway alone. CONCLUSIONS: Our study demonstrated that age-dependent p16 epimutation creates a permissive microenvironment for malignant transformation of polyps to colon cancer. Our findings provide a mechanistic rationale for future targeted therapy in patients with p16 epimutation.

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.

Hypermethylation of PI3K-AKT signalling pathway genes is associated with human neural tube defects.

Neural tube defects (NTDs) are a group of common and severe congenital malformations. The PI3K-AKT signalling pathway plays a crucial role in the neural tube development. There is limited evidence concerning any possible association between aberrant methylation in PI3K-AKT signalling pathway genes and NTDs. Therefore, we aimed to investigate potential associations between aberrant methylation of PI3K-AKT pathway genes and NTDs. Methylation studies of PI3K-AKT pathway genes utilizing microarray genome-methylation data derived from neural tissues of ten NTD cases and eight non-malformed controls were performed. Targeted DNA methylation analysis was subsequently performed in an independent cohort of 73 NTD cases and 32 controls to validate the methylation levels of identified genes. siRNAs were used to pull-down the target genes in human embryonic stem cells (hESCs) to examine the effects of the aberrant expression of target genes on neural cells. As a result, 321 differentially hypermethylated CpG sites in the promoter regions of 30 PI3K-AKT pathway genes were identified in the microarray data. In target methylation analysis, CHRM1, FGF19, and ITGA7 were confirmed to be significantly hypermethylated in NTD cases and were associated with increased risk for NTDs. The down-regulation of FGF19, CHRM1, and ITGA7 impaired the formation of rosette-like cell aggregates. The down-regulation of those three genes affected the expression of PAX6, SOX2 and MAP2, implying their influence on the differentiation of neural cells. This study for the first time reported that hypermethylation of PI3K-AKT pathway genes such as CHRM1, FGF19, and ITGA7 is associated with human NTDs.

CIC de novo loss of function variants contribute to cerebral folate deficiency by downregulating FOLR1 expression.

Background Cerebral folate deficiency (CFD) syndrome is characterised by a low concentration of 5-methyltetrahydrofolate in cerebrospinal fluid, while folate levels in plasma and red blood cells are in the low normal range. Mutations in several folate pathway genes, including FOLR1 (folate receptor alpha, FRα), DHFR (dihydrofolate reductase) and PCFT (proton coupled folate transporter) have been previously identified in patients with CFD. Methods In an effort to identify causal mutations for CFD, we performed whole exome sequencing analysis on eight CFD trios and identified eight de novo mutations in seven trios. Results Notably, we found a de novo stop gain mutation in the capicua (CIC) gene. Using 48 sporadic CFD samples as a validation cohort, we identified three additional rare variants in CIC that are putatively deleterious mutations. Functional analysis indicates that CIC binds to an octameric sequence in the promoter regions of folate transport genes: FOLR1, PCFT and reduced folate carrier (Slc19A1; RFC1). The CIC nonsense variant (p.R353X) downregulated FOLR1 expression in HeLa cells as well as in the induced pluripotent stem cell (iPSCs) derived from the original CFD proband. Folate binding assay demonstrated that the p.R353X variant decreased cellular binding of folic acid in cells. Conclusion This study indicates that CIC loss of function variants can contribute to the genetic aetiology of CFD through regulating FOLR1 expression. Our study described the first mutations in a non-folate pathway gene that can contribute to the aetiology of CFD.

Maternal Hypertension-Related Genotypes and Congenital Heart Defects.

BACKGROUND: Maternal hypertension has been associated with congenital heart defect occurrence in several studies. We assessed whether maternal genotypes associated with this condition were also associated with congenital heart defect occurrence. METHODS: We used data from the National Birth Defects Prevention Study to identify non-Hispanic white (NHW) and Hispanic women with (cases) and without (controls) a pregnancy in which a select simple, isolated heart defect was present between 1999 and 2011. We genotyped 29 hypertension-related single nucleotide polymorphisms (SNPs). We conducted logistic regression analyses separately by race/ethnicity to assess the relationship between the presence of any congenital heart defect and each SNP and an overall blood pressure genetic risk score (GRS). All analyses were then repeated to assess 4 separate congenital heart defect subtypes. RESULTS: Four hypertension-related variants were associated with congenital heart defects among NHW women (N = 1,568 with affected pregnancies). For example, 1 intronic variant in ARHGAP2, rs633185, was associated with conotruncal defects (odds ratio [OR]: 1.3, 95% confidence interval [CI]: 1.1-1.6). Additionally, 2 variants were associated with congenital heart defects among Hispanic women (N = 489 with affected pregnancies). The GRS had a significant association with septal defects (OR: 2.1, 95% CI: 1.2-3.5) among NHW women. CONCLUSIONS: We replicated a previously reported association between rs633185 and conotruncal defects. Although additional hypertension-related SNPs were also associated with congenital heart defects, more work is needed to better understand the relationship between genetic risk for maternal hypertension and congenital heart defects occurrence.

FKBP8 variants are risk factors for spina bifida.

Neural tube defects (NTDs) are a group of severe congenital malformations caused by a failure of neural tube closure during early embryonic development. Although extensively investigated, the genetic etiology of NTDs remains poorly understood. FKBP8 is critical for proper mammalian neural tube closure. Fkbp8-/- mouse embryos showed posterior NTDs consistent with a diagnosis of spina bifida (SB). To date, no publication has reported any association between FKBP8 and human NTDs. Using Sanger sequencing on genomic DNA samples from 472 SB and 565 control samples, we identified five rare (MAF ≤ 0.001) deleterious variants in SB patients, while no rare deleterious variant was identified in the controls (P = 0.0191). p.Glu140* affected FKBP8 localization to the mitochondria and created a truncated form of the FKBP8 protein, thus impairing its interaction with BCL2 and ultimately leading to an increase in cellular apoptosis. p.Ser3Leu, p.Lys315Asn and p.Ala292Ser variants decreased FKBP8 protein level. p.Lys315Asn further increased the cellular apoptosis. RNA sequencing on anterior and posterior tissues isolated from Fkbp8-/- and wildtype mice at E9.5 and E10.5 showed that Fkbp8-/- embryos have an abnormal expression profile within tissues harvested at posterior sites, thus leading to a posterior NTD. Moreover, we found that Fkbp8 knockout mouse embryos have abnormal expression of Wnt3a and Nkx2.9 during the early stage of neural tube development, perhaps also contributing to caudal specific NTDs. These findings provide evidence that functional variants of FKBP8 are risk factors for SB, which may involve a novel mechanism by which Fkbp8 mutations specifically cause SB in mice.

DVL mutations identified from human neural tube defects and Dandy-Walker malformation obstruct the Wnt signaling pathway.

Wnt signaling pathways, including the canonical Wnt/ß-catenin pathway, planar cell polarity pathway, and Wnt/Ca2+ signaling pathway, play important roles in neural development during embryonic stages. The DVL genes encode the hub proteins for Wnt signaling pathways. The mutations in DVL2 and DVL3 were identified from patients with neural tube defects (NTDs), but their functions in the pathogenesis of human neural diseases remain elusive. Here, we sequenced the coding regions of three DVL genes in 176 stillborn or miscarried fetuses with NTDs or Dandy-Walker malformation (DWM) and 480 adult controls from a Han Chinese population. Four rare mutations were identified: DVL1 p.R558H, DVL1 p.R606C, DVL2 p.R633W, and DVL3 p.R222Q. To assess the effect of these mutations on NTDs and DWM, various functional analyses such as luciferase reporter assay, stress fiber formation, and in vivo teratogenic assay were performed. The results showed that the DVL2 p.R633W mutation destabilized DVL2 protein and upregulated activities for all three Wnt signalings (Wnt/ß-catenin signaling, Wnt/planar cell polarity signaling, and Wnt/Ca2+ signaling) in mammalian cells. In contrast, DVL1 mutants (DVL1 p.R558H and DVL1 p.R606C) decreased canonical Wnt/ß-catenin signaling but increased the activity of Wnt/Ca2+ signaling, and DVL3 p.R222Q only decreased the activity of Wnt/Ca2+ signaling. We also found that only the DVL2 p.R633W mutant displayed more severe teratogenicity in zebrafish embryos than wild-type DVL2. Our study demonstrates that these four rare DVL mutations, especially DVL2 p.R633W, may contribute to human neural diseases such as NTDs and DWM by obstructing Wnt signaling pathways.

Somatic mutations in planar cell polarity genes in neural tissue from human fetuses with neural tube defects.

Extensive studies that have sought causative mutation(s) for neural tube defects (NTDs) have yielded limited positive findings to date. One possible reason for this is that many studies have been confined to analyses of germline mutations and so may have missed other, non-germline mutations in NTD cases. We hypothesize that somatic mutations of planar polarity pathway (PCP) genes may play a role in the development of NTDs. Torrent™ Personal Genome Machine™ (PGM) sequencing was designed for selected PCP genes in paired DNA samples extracted from the tissues of lesion sites and umbilical cord from 48 cases. Sanger sequencing was used to validate the detected mutations. The source and distribution of the validated mutations in tissues from different germ layers were investigated. Subcellular location, western blotting, and luciferase assays were performed to better understand the effects of the mutations on protein localization, protein level, and pathway signaling. ix somatic mutations were identified and validated, which showed diverse distributions in different tissues. Three somatic mutations were novel/rare: CELSR1 p.Gln2125His, FZD6 p.Gln88Glu, and VANGL1 p.Arg374His. FZD6 p.Gln88Glu caused mislocalization of its protein from the cytoplasm to the nucleus, and disrupted the colocalization of CELSR1 and FZD6. This mutation affected non-canonical WNT signaling in luciferase assays. VANGL1 p.Arg374His impaired the co-localization of CELSR1 and VANGL1, increased the protein levels of VANGL1, and influenced cell migration. In all, 7/48 (14.5%) of the studied NTD cases contained somatic PCP mutations. Somatic mutations in PCP genes (e.g., FZD6 and VANGL1) are associated with human NTDs, and they may occur in different stages and regions during embryonic development, resulting in a varied distribution in fetal tissues/organs.

Rare copy number variations of planar cell polarity genes are associated with human neural tube defects.

Select single-nucleotide variants in planar cell polarity (PCP) genes are associated with increased risk for neural tube defects (NTDs). However, whether copy number variants (CNVs) in PCP genes contribute to NTDs is unknown. Considering that CNVs are implicated in several human developmental disorders, we hypothesized that CNVs in PCP genes may be causative factors to human NTDs. DNA from umbilical cord tissues of NTD-affected fetuses and parental venous blood samples were collected. We performed a quantitative analysis of copy numbers of all exon regions in the VANGL1, VANGL2, CELSR1, SCRIB, DVL2, DVL3, and PTK7 genes using a CNVplex assay. Quantitative real-time PCR (qPCR) was carried out to confirm the results of CNV analysis. As a result, 16 CNVs were identified among the NTDs. Of these CNVs, 5 loci were identified in 11 NTD probands with CNVs involving DVL2 (exons 1-15), VANGL1 (exons 1-7, exon 8), and VANGL2 (exons 5-8, exons 7 and 8). One CNV (DVL2 exons 1-15) was a duplication and the remaining 15 CNVs were deletions. Eleven CNVs were confirmed by qPCR. One de novo CNV in VANGL1 and one DVL2 were detected from two cases. Compared with unaffected control populations in 1000 Genome, ExAC, MARRVEL, DGV, and dbVar databases, the frequencies of de novo deletion in VANGL1 (1.14%) and de novo duplication in DVL2 (0.57%) were significantly higher in our NTD subjects (p < 0.05). This study demonstrates that de novo CNVs in PCP genes, notably deletions in VANGL1 and gains in DVL2, could contribute to the risk of NTDs.

Approaches to studying the genomic architecture of complex birth defects.

Every year nearly 6 percent of children worldwide are born with a serious congenital malformation, resulting in death or lifelong disability. In the United States, birth defects remain one of the leading causes of infant mortality. Among the common structural congenital defects are conditions known as neural tube defects (NTDs). These are a class of malformation of the brain and spinal cord where the neural tube fails to close during the neurulation. Although NTDs remain among the most pervasive and debilitating of all human developmental anomalies, there is insufficient understanding of their etiology. Previous studies have proposed that complex birth defects like NTDs are likely omnigenic, involving interconnected gene regulatory networks with associated signals throughout the genome. Advances in technologies have allowed researchers to more critically investigate regulatory gene networks in ever increasing detail, informing our understanding of the genetic basis of NTDs. Employing a systematic analysis of these complex birth defects using massively parallel DNA sequencing with stringent bioinformatic algorithms, it is possible to approach a greater level of understanding of the genomic architecture underlying NTDs. Herein, we present a brief overview of different approaches undertaken in our laboratory to dissect out the genetics of susceptibility to NTDs. This involves the use of mouse models to identify candidate genes, as well as large scale whole genome/whole exome (WGS/WES) studies to interrogate the genomic landscape of NTDs. The goal of this research is to elucidate the gene-environment interactions contributing to NTDs, thus encouraging global research efforts in their prevention.

Loss of RAD9B impairs early neural development and contributes to the risk for human spina bifida.

DNA damage response (DDR) genes orchestrating the network of DNA repair, cell cycle control, are essential for the rapid proliferation of neural progenitor cells. To date, the potential association between specific DDR genes and the risk of human neural tube defects (NTDs) has not been investigated. Using whole-genome sequencing and targeted sequencing, we identified significant enrichment of rare deleterious RAD9B variants in spina bifida cases compared to controls (8/409 vs. 0/298; p = .0241). Among the eight identified variants, the two frameshift mutants and p.Gln146Glu affected RAD9B nuclear localization. The two frameshift mutants also decreased the protein level of RAD9B. p.Ser354Gly, as well as the two frameshifts, affected the cell proliferation rate. Finally, p.Ser354Gly, p.Ser10Gly, p.Ile112Met, p.Gln146Glu, and the two frameshift variants showed a decreased ability for activating JNK phosphorylation. RAD9B knockdowns in human embryonic stem cells profoundly affected early differentiation through impairing PAX6 and OCT4 expression. RAD9B deficiency impeded in vitro formation of neural organoids, a 3D cell culture model for human neural development. Furthermore, the RNA-seq data revealed that loss of RAD9B dysregulates cell adhesion genes during organoid formation. These results represent the first demonstration of a DDR gene as an NTD risk factor in humans.

Maternal genetic markers for risk of celiac disease and their potential association with neural tube defects in offspring.

BACKGROUND: We examined the association between the maternal genotype for celiac disease-associated variants and risk of neural tube defects (NTDs). METHODS: We conducted a case-control study, using data from the National Birth Defects Prevention Study. We evaluated 667 cases (women with an offspring with NTD) and 743 controls (women with an offspring without a birth defect). We classified women as having low, intermediate, or high risk of celiac disease based on human leukocyte antigen (HLA) variants. We used logistic regression to assess the relationship between HLA celiac risk group (low, intermediate, high) and risk of NTDs. Fifteen non-HLA variants (identified from genome-wide association studies of celiac disease) were individually evaluated and modeled additively. RESULTS: There was no association between HLA celiac risk group and NTDs (intermediate vs. low risk: aOR, 1.0; 95% CI, 0.8-1.3; high vs. low risk: aOR, 0.8; 95% CI, 0.5-1.3). Of the fifteen non-HLA variants, we observed five significant associations after accounting for multiple comparisons. Three negative associations were observed with rs10903122, rs13314993, rs13151961 (aOR range: 0.69-0.81), and two positive associations were observed with rs13003464 and rs11221332 (aOR range: 1.27-1.73). CONCLUSION: If confirmed, our results suggest that the maternal variants related to celiac disease may be involved in the risk of NTDs.

The TFAP2A-IRF6-GRHL3 genetic pathway is conserved in neurulation.

Mutations in IRF6, TFAP2A and GRHL3 cause orofacial clefting syndromes in humans. However, Tfap2a and Grhl3 are also required for neurulation in mice. Here, we found that homeostasis of Irf6 is also required for development of the neural tube and associated structures. Over-expression of Irf6 caused exencephaly, a rostral neural tube defect, through suppression of Tfap2a and Grhl3 expression. Conversely, loss of Irf6 function caused a curly tail and coincided with a reduction of Tfap2a and Grhl3 expression in tail tissues. To test whether Irf6 function in neurulation was conserved, we sequenced samples obtained from human cases of spina bifida and anencephaly. We found two likely disease-causing variants in two samples from patients with spina bifida. Overall, these data suggest that the Tfap2a-Irf6-Grhl3 genetic pathway is shared by two embryologically distinct morphogenetic events that previously were considered independent during mammalian development. In addition, these data suggest new candidates to delineate the genetic architecture of neural tube defects and new therapeutic targets to prevent this common birth defect.

Maternal Lactase Polymorphism (rs4988235) Is Associated with Neural Tube Defects in Offspring in the National Birth Defects Prevention Study.

BACKGROUND: The risk of neural tube defect (NTD)-affected pregnancies is reduced with adequate folic acid intake during early pregnancy. However, NTDs have been observed among offspring of women with adequate folic acid intake. Some of these women are possibly not absorbing enough folic acid. Because lactase deficiency can lead to poor nutrient absorption, we hypothesized that lactase-deficient women will be at increased risk of having offspring with NTDs. OBJECTIVE: We examined the association between maternal rs4988235 (a lactase deficiency genetic marker) and NTDs in offspring. METHODS: We conducted a case-control study using data from the National Birth Defects Prevention Study, United States, 1997-2009, restricting to non-Hispanic white (NHW) and Hispanic women. Cases were women with an offspring with an NTD (n = 378 NHW, 207 Hispanic), and controls were women with an offspring without a birth defect (n = 461 NHW, 165 Hispanic). Analyses were conducted separately by race/ethnicity, using logistic regression. Women with the CC genotype were categorized as being lactase deficient. To assess potential effect modification, analyses were stratified by lactose intake, folic acid supplementation, dietary folate, and diet quality. RESULTS: Among NHW women, the odds of being lactase deficient were greater among cases compared with controls (OR: 1.37; 95% CI: 1.02, 1.82). Among Hispanic women, the odds of being lactase deficient were significantly lower among cases compared with controls (OR: 0.50, 95% CI: 0.33, 0.77). The association differed when stratified by lactose intake in NHW women (higher odds among women who consumed ≥12 g lactose/1000 kcal) and by dietary folate in Hispanic women (opposite direction of associations). The association did not differ when stratified by folic acid supplementation or diet quality. CONCLUSIONS: Our findings suggest that maternal lactase deficiency is associated with NTDs in offspring. However, we observed opposite directions of effect by race/ethnicity that could not be definitively explained.

Dominant negative GPR161 rare variants are risk factors of human spina bifida.

Spina bifida (SB) is a complex disorder of failed neural tube closure during the first month of human gestation, with a suspected etiology involving multiple gene and environmental interactions. GPR161 is a ciliary G-protein coupled receptor that regulates Sonic Hedgehog (Shh) signaling. Gpr161 null and hypomorphic mutations cause neural tube defects (NTDs) in mouse models. Herein we show that several genes involved in Shh and Wnt signaling were differentially expressed in the Gpr161 null embryos using RNA-seq analysis. To determine whether there exists an association between GPR161 and SB in humans, we performed direct Sanger sequencing on the GPR161 gene in a cohort of 384 SB patients and 190 healthy controls. We identified six rare variants of GPR161 in six SB cases, of which two of the variants were novel and did not exist in any databases. Both of these variants were predicted to be damaging by SIFT and/or PolyPhen analysis. The novel GPR161 rare variants mislocalized to the primary cilia, dysregulated Shh and Wnt signaling and inhibited cell proliferation in vitro. Our results demonstrate that GPR161 mutations cause NTDs via dysregulation of Shh and Wnt signaling in mice, and novel rare variants of GPR161 can be risk factors for SB in humans.

Gene variants in the folate pathway are associated with increased levels of folate receptor autoantibodies.

BACKGROUND: Folate receptors (FRs) facilitate embryonic uptake of folates and are important for proper early embryonic development. There is accumulating evidence that blocking FR autoantibodies contribute to developmental diseases. However, genetic factors associated with the expression of FR autoantibodies remain unknown. OBJECTIVE: We investigated the effects of genetic polymorphisms in folate pathway genes on FR autoantibody titers in women. METHODS: We recruited 302 pregnant women in China. The FR antigen-down immunoassay was used to measure levels of FR autoantibodies including human immunoglobulin G (IgG) and immunoglobulin M (IgM) in maternal plasma. Genotypes were identified by matrix-assisted laser desorption/ionization time of flight mass spectrometry and polymerase chain reaction methods. General linear model was used to analyze the effects of genetic variants on FR autoantibody levels. RESULTS: Significant associations were observed between genotypic variations and levels of FR autoantibodies. Plasma levels of FR autoantibodies in women with the TT genotype at MTHFR rs1801133 were significantly higher than those of women with the CC genotype (IgG: ß = 0.62, 95% CI 0.21-1.04; IgM: ß = 0.42, 95% CI 0.12-0.72). For DNMT3A rs7560488, the level of FR autoantibody IgG significantly increased in the TT genotype compared with CC genotype (ß = 0.90, 95% CI 0.20-1.59). For MTHFD2 rs828903, genotype GG was associated with elevated levels of FR autoantibody IgM compared to the AA genotype (ß = 0.60, 95% CI 0.10-1.10). No association was detected between genetic variants of the DHFR gene with FR autoantibodies levels. CONCLUSION: Genetic variations in MTHFR, DNMT3A, and MTHFD2 genes were associated with elevated plasma levels of FR autoantibodies.