Emergence of an early SARS-CoV-2 epidemic in the United States.

The emergence of the COVID-19 epidemic in the United States (U.S.) went largely undetected due to inadequate testing. New Orleans experienced one of the earliest and fastest accelerating outbreaks, coinciding with Mardi Gras. To gain insight into the emergence of SARS-CoV-2 in the U.S. and how large-scale events accelerate transmission, we sequenced SARS-CoV-2 genomes during the first wave of the COVID-19 epidemic in Louisiana. We show that SARS-CoV-2 in Louisiana had limited diversity compared to other U.S. states and that one introduction of SARS-CoV-2 led to almost all of the early transmission in Louisiana. By analyzing mobility and genomic data, we show that SARS-CoV-2 was already present in New Orleans before Mardi Gras, and the festival dramatically accelerated transmission. Our study provides an understanding of how superspreading during large-scale events played a key role during the early outbreak in the U.S. and can greatly accelerate epidemics.

A genome sequencing system for universal newborn screening, diagnosis, and precision medicine for severe genetic diseases.

Newborn screening (NBS) dramatically improves outcomes in severe childhood disorders by treatment before symptom onset. In many genetic diseases, however, outcomes remain poor because NBS has lagged behind drug development. Rapid whole-genome sequencing (rWGS) is attractive for comprehensive NBS because it concomitantly examines almost all genetic diseases and is gaining acceptance for genetic disease diagnosis in ill newborns. We describe prototypic methods for scalable, parentally consented, feedback-informed NBS and diagnosis of genetic diseases by rWGS and virtual, acute management guidance (NBS-rWGS). Using established criteria and the Delphi method, we reviewed 457 genetic diseases for NBS-rWGS, retaining 388 (85%) with effective treatments. Simulated NBS-rWGS in 454,707 UK Biobank subjects with 29,865 pathogenic or likely pathogenic variants associated with 388 disorders had a true negative rate (specificity) of 99.7% following root cause analysis. In 2,208 critically ill children with suspected genetic disorders and 2,168 of their parents, simulated NBS-rWGS for 388 disorders identified 104 (87%) of 119 diagnoses previously made by rWGS and 15 findings not previously reported (NBS-rWGS negative predictive value 99.6%, true positive rate [sensitivity] 88.8%). Retrospective NBS-rWGS diagnosed 15 children with disorders that had been undetected by conventional NBS. In 43 of the 104 children, had NBS-rWGS-based interventions been started on day of life 5, the Delphi consensus was that symptoms could have been avoided completely in seven critically ill children, mostly in 21, and partially in 13. We invite groups worldwide to refine these NBS-rWGS conditions and join us to prospectively examine clinical utility and cost effectiveness.

A gene-based association test of interactions for maternal-fetal genotypes identifies genes associated with nonsyndromic congenital heart defects.

The risk of congenital heart defects (CHDs) may be influenced by maternal genes, fetal genes, and their interactions. Existing methods commonly test the effects of maternal and fetal variants one-at-a-time and may have reduced statistical power to detect genetic variants with low minor allele frequencies. In this article, we propose a gene-based association test of interactions for maternal-fetal genotypes (GATI-MFG) using a case-mother and control-mother design. GATI-MFG can integrate the effects of multiple variants within a gene or genomic region and evaluate the joint effect of maternal and fetal genotypes while allowing for their interactions. In simulation studies, GATI-MFG had improved statistical power over alternative methods, such as the single-variant test and functional data analysis (FDA) under various disease scenarios. We further applied GATI-MFG to a two-phase genome-wide association study of CHDs for the testing of both common variants and rare variants using 947 CHD case mother-infant pairs and 1306 control mother-infant pairs from the National Birth Defects Prevention Study (NBDPS). After Bonferroni adjustment for 23,035 genes, two genes on chromosome 17, TMEM107 (p = 1.64e-06) and CTC1 (p = 2.0e-06), were identified for significant association with CHD in common variants analysis. Gene TMEM107 regulates ciliogenesis and ciliary protein composition and was found to be associated with heterotaxy. Gene CTC1 plays an essential role in protecting telomeres from degradation, which was suggested to be associated with cardiogenesis. Overall, GATI-MFG outperformed the single-variant test and FDA in the simulations, and the results of application to NBDPS samples are consistent with existing literature supporting the association of TMEM107 and CTC1 with CHDs.

Paternal genetic variants and risk of obstructive heart defects: A parent-of-origin approach.

Previous research on risk factors for obstructive heart defects (OHDs) focused on maternal and infant genetic variants, prenatal environmental exposures, and their potential interaction effects. Less is known about the role of paternal genetic variants or environmental exposures and risk of OHDs. We examined parent-of-origin effects in transmission of alleles in the folate, homocysteine, or transsulfuration pathway genes on OHD occurrence in offspring. We used data on 569 families of liveborn infants with OHDs born between October 1997 and August 2008 from the National Birth Defects Prevention Study to conduct a family-based case-only study. Maternal, paternal, and infant DNA were genotyped using an Illumina Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR), 95% confidence interval (CI), and likelihood ratio tests from log-linear models were used to estimate the parent-of-origin effect of 877 SNPs in 60 candidate genes in the folate, homocysteine, and transsulfuration pathways on the risk of OHDs. Bonferroni correction was applied for multiple testing. We identified 3 SNPs in the transsulfuration pathway and 1 SNP in the folate pathway that were statistically significant after Bonferroni correction. Among infants who inherited paternally-derived copies of the G allele for rs6812588 in the RFC1 gene, the G allele for rs1762430 in the MGMT gene, and the A allele for rs9296695 and rs4712023 in the GSTA3 gene, RRs for OHD were 0.11 (95% CI: 0.04, 0.29, P = 9.16x10-7), 0.30 (95% CI: 0.17, 0.53, P = 9.80x10-6), 0.34 (95% CI: 0.20, 0.57, P = 2.28x10-5), and 0.34 (95% CI: 0.20, 0.58, P = 3.77x10-5), respectively, compared to infants who inherited maternally-derived copies of the same alleles. We observed statistically significant decreased risk of OHDs among infants who inherited paternal gene variants involved in folate and transsulfuration pathways.

Detecting methylation quantitative trait loci using a methylation random field method.

DNA methylation may be regulated by genetic variants within a genomic region, referred to as methylation quantitative trait loci (mQTLs). The changes of methylation levels can further lead to alterations of gene expression, and influence the risk of various complex human diseases. Detecting mQTLs may provide insights into the underlying mechanism of how genotypic variations may influence the disease risk. In this article, we propose a methylation random field (MRF) method to detect mQTLs by testing the association between the methylation level of a CpG site and a set of genetic variants within a genomic region. The proposed MRF has two major advantages over existing approaches. First, it uses a beta distribution to characterize the bimodal and interval properties of the methylation trait at a CpG site. Second, it considers multiple common and rare genetic variants within a genomic region to identify mQTLs. Through simulations, we demonstrated that the MRF had improved power over other existing methods in detecting rare variants of relatively large effect, especially when the sample size is small. We further applied our method to a study of congenital heart defects with 83 cardiac tissue samples and identified two mQTL regions, MRPS10 and PSORS1C1, which were colocalized with expression QTL in cardiac tissue. In conclusion, the proposed MRF is a useful tool to identify novel mQTLs, especially for studies with limited sample sizes.

Random field modeling of multi-trait multi-locus association for detecting methylation quantitative trait loci.

MOTIVATION: CpG sites within the same genomic region often share similar methylation patterns and tend to be co-regulated by multiple genetic variants that may interact with one another. RESULTS: We propose a multi-trait methylation random field (multi-MRF) method to evaluate the joint association between a set of CpG sites and a set of genetic variants. The proposed method has several advantages. First, it is a multi-trait method that allows flexible correlation structures between neighboring CpG sites (e.g. distance-based correlation). Second, it is also a multi-locus method that integrates the effect of multiple common and rare genetic variants. Third, it models the methylation traits with a beta distribution to characterize their bimodal and interval properties. Through simulations, we demonstrated that the proposed method had improved power over some existing methods under various disease scenarios. We further illustrated the proposed method via an application to a study of congenital heart defects (CHDs) with 83 cardiac tissue samples. Our results suggested that gene BACE2, a methylation quantitative trait locus (QTL) candidate, colocalized with expression QTLs in artery tibial and harbored genetic variants with nominal significant associations in two genome-wide association studies of CHD. AVAILABILITY AND IMPLEMENTATION: https://github.com/chenlyu2656/Multi-MRF. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel association of Dandy-Walker malformation with CAPN15 variants expands the phenotype of oculogastrointestinal neurodevelopmental syndrome.

Oculogastrointestinal neurodevelopmental syndrome has been described in seven previously published individuals who harbor biallelic pathogenic variants in the CAPN15 gene. Biallelic missense variants have been reported to demonstrate a phenotype of eye abnormalities and developmental delay, while biallelic loss of function variants exhibit phenotypes including microcephaly and craniofacial abnormalities, cardiac and genitourinary malformations, and abnormal neurologic activity. We report six individuals from three unrelated families harboring biallelic deleterious variants in CAPN15 with phenotypes overlapping those previously described for this disorder. Of the individuals affected, four demonstrate radiographic evidence of the classical triad of Dandy-Walker malformation including hypoplastic vermis, fourth ventricle enlargement, and torcular elevation. Cerebellar anomalies have not been previously reported in association with CAPN15-related disease. Here, we present three unrelated families with findings consistent with oculogastrointestinal neurodevelopmental syndrome and cerebellar pathology including Dandy-Walker malformation. To corroborate these novel clinical findings, we present supporting data from the mouse model suggesting an important role for this protein in normal cerebellar development. Our findings add six molecularly confirmed cases to the literature and additionally establish a new association of Dandy-Walker malformation with biallelic CAPN15 variants, thereby expanding the neurologic spectrum among patients affected by CAPN15-related disease.

Innovative approach to identify multigenomic and environmental interactions associated with birth defects in family-based hybrid designs.

Genes, including those with transgenerational effects, work in concert with behavioral, environmental, and social factors via complex biological networks to determine human health. Understanding complex relationships between causal factors underlying human health is an essential step towards deciphering biological mechanisms. We propose a new analytical framework to investigate the interactions between maternal and offspring genetic variants or their surrogate single nucleotide polymorphisms (SNPs) and environmental factors using family-based hybrid study design. The proposed approach can analyze diverse genetic and environmental factors and accommodate samples from a variety of family units, including case/control-parental triads, and case/control-parental dyads, while minimizing potential bias introduced by population admixture. Comprehensive simulations demonstrated that our innovative approach outperformed the log-linear approach, the best available method for case-control family data. The proposed approach had greater statistical power and was capable to unbiasedly estimate the maternal and child genetic effects and the effects of environmental factors, while controlling the Type I error rate against population stratification. Using our newly developed approach, we analyzed the associations between maternal and fetal SNPs and obstructive and conotruncal heart defects, with adjustment for demographic and lifestyle factors and dietary supplements. Fourteen and 11 fetal SNPs were associated with obstructive and conotruncal heart defects, respectively. Twenty-seven and 17 maternal SNPs were associated with obstructive and conotruncal heart defects, respectively. In addition, maternal body mass index was a significant risk factor for obstructive defects. The proposed approach is a powerful tool for interrogating the etiological mechanism underlying complex traits.

Maternal Smoking and Congenital Heart Defects, National Birth Defects Prevention Study, 1997-2011.

OBJECTIVES: To assess associations between maternal smoking and congenital heart defects (CHDs) in offspring. STUDY DESIGN: We performed a retrospective case-control study using data for cases of CHD (n = 8339) and nonmalformed controls (n = 11 020) from all years (1997-2011) of the National Birth Defects Prevention Study. Maternal self-reported smoking 1 month before through 3 months after conception was evaluated as a binary (none, any) and categorical (light, medium, heavy) exposure. Multivariable logistic regression was used to estimate aOR and 95% CIs. Stratified analyses were performed for septal defects according to maternal age, prepregnancy body mass index, and maternal race/ethnicity. RESULTS: Multiple CHDs displayed modest associations with any level of maternal periconceptional smoking independent of potential confounders; the strongest associations were for aggregated septal defects (OR, 1.5; 95% CI, 1.3-1.7), tricuspid atresia (OR, 1.7; 95% CI, 1.0-2.7), and double outlet right ventricle (DORV) (OR, 1.5; 95% CI, 1.1-2.1). Tricuspid atresia and DORV also displayed dose-response relationships. Among heavy smokers, the highest odds were again observed for tricuspid atresia (aOR 3.0; 95% CI, 1.5-6.1) and DORV (aOR 1.5; 95% CI, 1.1-2.2). Heavy smokers ≥35 years old more frequently had a child with a septal defect when compared with similarly aged nonsmokers (aOR 2.3; 95% CI, 1.4-3.9). CONCLUSIONS: Maternal periconceptional smoking is most strongly associated with septal defects, tricuspid atresia, and DORV; the risk for septal defects is modified by maternal age.

A genome-wide association study of obstructive heart defects among participants in the National Birth Defects Prevention Study.

Obstructive heart defects (OHDs) share common structural lesions in arteries and cardiac valves, accounting for ~25% of all congenital heart defects. OHDs are highly heritable, resulting from interplay among maternal exposures, genetic susceptibilities, and epigenetic phenomena. A genome-wide association study was conducted in National Birth Defects Prevention Study participants (Ndiscovery  = 3978; Nreplication  = 2507), investigating the genetic architecture of OHDs using transmission/disequilibrium tests (TDT) in complete case-parental trios (Ndiscovery_TDT  = 440; Nreplication_TDT  = 275) and case-control analyses separately in infants (Ndiscovery_CCI  = 1635; Nreplication_CCI  = 990) and mothers (case status defined by infant; Ndiscovery_CCM  = 1703; Nreplication_CCM  = 1078). In the TDT analysis, the SLC44A2 single nucleotide polymorphism (SNP) rs2360743 was significantly associated with OHD (pdiscovery  = 4.08 × 10-9 ; preplication  = 2.44 × 10-4 ). A CAPN11 SNP (rs55877192) was suggestively associated with OHD (pdiscovery  = 1.61 × 10-7 ; preplication  = 0.0016). Two other SNPs were suggestively associated (p < 1 × 10-6 ) with OHD in only the discovery sample. In the case-control analyses, no SNPs were genome-wide significant, and, even with relaxed thresholds ( × discovery < 1 × 10-5 and preplication < 0.05), only one SNP (rs188255766) in the infant analysis was associated with OHDs (pdiscovery  = 1.42 × 10-6 ; preplication  = 0.04). Additional SNPs with pdiscovery < 1 × 10-5 were in loci supporting previous findings but did not replicate. Overall, there was modest evidence of an association between rs2360743 and rs55877192 and OHD and some evidence validating previously published findings.

Exome sequencing of child-parent trios with bladder exstrophy: Findings in 26 children.

Bladder exstrophy (BE) is a rare, lower ventral midline defect with the bladder and part of the urethra exposed. The etiology of BE is unknown but thought to be influenced by genetic variation with more recent studies suggesting a role for rare variants. As such, we conducted paired-end exome sequencing in 26 child/mother/father trios. Three children had rare (allele frequency ≤ 0.0001 in several public databases) inherited variants in TSPAN4, one with a loss-of-function variant and two with missense variants. Two children had loss-of-function variants in TUBE1. Four children had rare missense or nonsense variants (one per child) in WNT3, CRKL, MYH9, or LZTR1, genes previously associated with BE. We detected 17 de novo missense variants in 13 children and three de novo loss-of-function variants (AKR1C2, PRRX1, PPM1D) in three children (one per child). We also detected rare compound heterozygous loss-of-function variants in PLCH2 and CLEC4M and rare inherited missense or loss-of-function variants in additional genes applying autosomal recessive (three genes) and X-linked recessive inheritance models (13 genes). Variants in two genes identified may implicate disruption in cell migration (TUBE1) and adhesion (TSPAN4) processes, mechanisms proposed for BE, and provide additional evidence for rare variants in the development of this defect.

A parent-of-origin analysis of paternal genetic variants and increased risk of conotruncal heart defects.

The association between conotruncal heart defects (CTHDs) and maternal genetic and environmental exposures is well studied. However, little is known about paternal genetic or environmental exposures and risk of CTHDs. We assessed the effect of paternal genetic variants in the folate, homocysteine, and transsulfuration pathways on risk of CTHDs in offspring. We utilized National Birth Defects Prevention Study data to conduct a family-based case only study using 616 live-born infants with CTHDs, born October 1997-August 2008. Maternal, paternal and infant DNA was genotyped using an Illumina® Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR) and 95% confidence intervals (CI) from log-linear models determined parent of origin effects for 921 SNPs in 60 candidate genes involved in the folate, homocysteine, and transsulfuration pathways on risk of CTHDs. The risk of CTHD among children who inherited a paternally derived copy of the A allele on GLRX (rs17085159) or the T allele of GLRX (rs12109442) was 0.23 (95%CI: 0.12, 0.42; p = 1.09 × 10-6 ) and 0.27 (95%CI: 0.14, 0.50; p = 2.06 × 10-5 ) times the risk among children who inherited a maternal copy of the same allele. The paternally inherited copy of the GSR (rs7818511) A allele had a 0.31 (95%CI: 0.18, 0.53; p = 9.94 × 10-6 ] risk of CTHD compared to children with the maternal copy of the same allele. The risk of CTHD is less influenced by variants in paternal genes involved in the folate, homocysteine, or transsulfuration pathways than variants in maternal genes in those pathways.

Testing Allele Transmission of an SNP Set Using a Family-Based Generalized Genetic Random Field Method.

Family-based association studies are commonly used in genetic research because they can be robust to population stratification (PS). Recent advances in high-throughput genotyping technologies have produced a massive amount of genomic data in family-based studies. However, current family-based association tests are mainly focused on evaluating individual variants one at a time. In this article, we introduce a family-based generalized genetic random field (FB-GGRF) method to test the joint association between a set of autosomal SNPs (i.e., single-nucleotide polymorphisms) and disease phenotypes. The proposed method is a natural extension of a recently developed GGRF method for population-based case-control studies. It models offspring genotypes conditional on parental genotypes, and, thus, is robust to PS. Through simulations, we presented that under various disease scenarios the FB-GGRF has improved power over a commonly used family-based sequence kernel association test (FB-SKAT). Further, similar to GGRF, the proposed FB-GGRF method is asymptotically well-behaved, and does not require empirical adjustment of the type I error rates. We illustrate the proposed method using a study of congenital heart defects with family trios from the National Birth Defects Prevention Study (NBDPS).

A Three-Way Interaction among Maternal and Fetal Variants Contributing to Congenital Heart Defects.

Congenital heart defects (CHDs) develop through a complex interplay between genetic variants, epigenetic modifications, and maternal environmental exposures. Genetic studies of CHDs have commonly tested single genetic variants for association with CHDs. Less attention has been given to complex gene-by-gene and gene-by-environment interactions. In this study, we applied a recently developed likelihood-ratio Mann-Whitney (LRMW) method to detect joint actions among maternal variants, fetal variants, and maternal environmental exposures, allowing for high-order statistical interactions. All subjects are participants from the National Birth Defect Prevention Study, including 623 mother-offspring pairs with CHD-affected pregnancies and 875 mother-offspring pairs with unaffected pregnancies. Each individual has 872 single nucleotide polymorphisms encoding for critical enzymes in the homocysteine, folate, and trans-sulfuration pathways. By using the LRMW method, three variants (fetal rs625879, maternal rs2169650, and maternal rs8177441) were identified with a joint association to CHD risk (nominal P-value = 1.13e-07). These three variants are located within genes BHMT2, GSTP1, and GPX3, respectively. Further examination indicated that maternal SNP rs2169650 may interact with both fetal SNP rs625879 and maternal SNP rs8177441. Our findings suggest that the risk of CHD may be influenced by both the intragenerational interaction within the maternal genome and the intergenerational interaction between maternal and fetal genomes.

Association between antibiotic use among pregnant women with urinary tract infections in the first trimester and birth defects, National Birth Defects Prevention Study 1997 to 2011.

BACKGROUND: Previous studies noted associations between birth defects and some antibiotics (e.g., nitrofurantoin, sulfonamides) but not others (e.g., penicillins). It is unclear if previous findings were due to antibiotic use, infections, or chance. To control for potential confounding by indication, we examined associations between antibiotic use and birth defects, among women reporting urinary tract infections (UTIs). METHODS: The National Birth Defects Prevention Study is a multi-site, population-based case-control study. Case infants/fetuses have any of over 30 major birth defects and controls are live-born infants without major birth defects. We analyzed pregnancies from 1997 to 2011 to estimate the association between maternally reported periconceptional (month before conception through the third month of pregnancy) use of nitrofurantoin, trimethoprim-sulfamethoxazole, or cephalosporins and specific birth defects, among women with periconceptional UTIs. Women with periconceptional UTIs who reported penicillin use served as the comparator. RESULTS: Periconceptional UTIs were reported by 7.8% (2029/26,068) of case and 6.7% (686/10,198) of control mothers. Most (68.2% of case, 66.6% of control mothers) also reported antibiotic use. Among 608 case and 231 control mothers reporting at least one periconceptional UTI and certain antibiotic use, compared with penicillin, nitrofurantoin use was associated with oral clefts in the offspring (adjusted odds ratio, 1.97 [95% confidence interval, 1.10-3.53]), trimethoprim-sulfamethoxazole use with esophageal atresia (5.31 [1.39-20.24]) and diaphragmatic hernia (5.09 [1.20-21.69]), and cephalosporin use with anorectal atresia/stenosis (5.01 [1.34-18.76]). CONCLUSION: Periconceptional exposure to some antibiotics might increase the risk for certain birth defects. However, because individual birth defects are rare, absolute risks should drive treatment decisions.Birth Defects Research (Part A) 106:940-949, 2016.© 2016 Wiley Periodicals, Inc.

Detecting maternal-fetal genotype interactions associated with conotruncal heart defects: a haplotype-based analysis with penalized logistic regression.

Nonsyndromic congenital heart defects (CHDs) develop during embryogenesis as a result of a complex interplay between environmental exposures, genetics, and epigenetic causes. Genetic factors associated with CHDs may be attributed to either independent effects of maternal or fetal genes, or the intergenerational interactions between maternal and fetal genes. Detecting gene-by-gene interactions underlying complex diseases is a major challenge in genetic research. Detecting maternal-fetal genotype (MFG) interactions and differentiating them from the maternal/fetal main effects has presented additional statistical challenges due to correlations between maternal and fetal genomes. Traditionally, genetic variants are tested separately for maternal/fetal main effects and MFG interactions on a single-locus basis. We conducted a haplotype-based analysis with a penalized logistic regression framework to dissect the genetic effect associated with the development of nonsyndromic conotruncal heart defects (CTD). Our method allows simultaneous model selection and effect estimation, providing a unified framework to differentiate maternal/fetal main effect from the MFG interaction effect. In addition, the method is able to test multiple highly linked SNPs simultaneously with a configuration of haplotypes, which reduces the data dimensionality and the burden of multiple testing. By analyzing a dataset from the National Birth Defects Prevention Study (NBDPS), we identified seven genes (GSTA1, SOD2, MTRR, AHCYL2, GCLC, GSTM3, and RFC1) associated with the development of CTDs. Our findings suggest that MFG interactions between haplotypes in three of seven genes, GCLC, GSTM3, and RFC1, are associated with nonsyndromic conotruncal heart defects.

Obstructive heart defects associated with candidate genes, maternal obesity, and folic acid supplementation.

Right-sided and left-sided obstructive heart defects (OHDs) are subtypes of congenital heart defects, in which the heart valves, arteries, or veins are abnormally narrow or blocked. Previous studies have suggested that the development of OHDs involved a complex interplay between genetic variants and maternal factors. Using the data from 569 OHD case families and 1,644 control families enrolled in the National Birth Defects Prevention Study (NBDPS) between 1997 and 2008, we conducted an analysis to investigate the genetic effects of 877 single nucleotide polymorphisms (SNPs) in 60 candidate genes for association with the risk of OHDs, and their interactions with maternal use of folic acid supplements, and pre-pregnancy obesity. Applying log-linear models based on the hybrid design, we identified a SNP in methylenetetrahydrofolate reductase (MTHFR) gene (C677T polymorphism) with a main genetic effect on the occurrence of OHDs. In addition, multiple SNPs in betaine-homocysteine methyltransferase (BHMT and BHMT2) were also identified to be associated with the occurrence of OHDs through significant main infant genetic effects and interaction effects with maternal use of folic acid supplements. We also identified multiple SNPs in glutamate-cysteine ligase, catalytic subunit (GCLC) and DNA (cytosine-5-)-methyltransferase 3 beta (DNMT3B) that were associated with elevated risk of OHDs among obese women. Our findings suggested that the risk of OHDs was closely related to a combined effect of variations in genes in the folate, homocysteine, or glutathione/transsulfuration pathways, maternal use of folic acid supplements and pre-pregnancy obesity.