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.