Front GenetFront GenetFront. Genet.Frontiers in Genetics1664-8021Frontiers Media S.A.24223580381858210.3389/fgene.2013.00223GeneticsOriginal Research ArticleFolate-related gene variants in Irish families affected by neural tube defectsFisk GreenRidgely1*ByrneJulianne2CriderKrista S.3GallagherMargaret4KoontzDeborah4BerryRobert J.31Carter Consulting, Inc. and Division of Birth Defects and Developmental Disabilities, National Center on Birth Defects and Developmental Disabilities, Centers for Disease Control and PreventionAtlanta, GA, USA2Boyne Research InstituteDrogheda, Ireland3Division of Birth Defects and Developmental Disabilities, National Center on Birth Defects and Developmental Disabilities, Centers for Disease Control and PreventionAtlanta, GA, USA4Newborn Screening and Molecular Biology Branch, Division of Laboratory Science, National Center for Environmental Health, Centers for Disease Control and PreventionAtlanta, GA, USA

Edited by: Zdenka Pausova, The Hospital for Sick Children, Canada

Reviewed by: Alice Hudder, Lake Erie College of Osteopathic Medicine, USA; Chris A. Murgatroyd, Manchester Metropolitan University, UK; Eva Reischl, Helmholtz Zentrum München - German Research Centre for Environmental Health GmbH, Germany

*Correspondence: Ridgely Fisk Green, National Center on Birth Defects and Developmental Disabilities, Centers for Disease Control and Prevention, 1600 Clifton Rd., MS E–86, Atlanta, GA 30333, USA e-mail: grf1@cdc.gov

This article was submitted to Epigenomics and Epigenetics, a section of the journal Frontiers in Genetics.

0611201320134223097201315102013Copyright © 2013 Fisk Green, Byrne, Crider, Gallagher, Koontz and Berry.2013This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Periconceptional folic acid use can often prevent neural tube defects (NTDs). Variants of genes involved in folate metabolism in mothers and children have been associated with occurrence of NTDs. We identified Irish families with individuals affected by neural tube defects. In these families, we observed that neural tube defects and birth defects overall occurred at a higher rate in the maternal lineage compared with the paternal lineage. The goal of this study was to look for evidence for genetic effects that could explain the discrepancy in the occurrence of these birth defects in the maternal vs. paternal lineage. We genotyped blood samples from 322 individuals from NTD-affected Irish families, identified through their membership in spina bifida associations. We looked for differences in distribution in maternal vs. paternal lineages of five genetic polymorphisms: the DHFR 19 bp deletion, MTHFD1 1958G>A, MTHFR 1298A>C, MTHFR 677C>T, and SLC19A1 80A>G. In addition to looking at genotypes individually, we determined the number of genotypes associated with decreased folate metabolism in each relative (“risk genotypes”) and compared the distribution of these genotypes in maternal vs. paternal relatives. Overall, maternal relatives had a higher number of genotypes associated with lower folate metabolism than paternal relatives (p = 0.017). We expected that relatives would share the same risk genotype as the individuals with NTDs and/or their mothers. However, we observed that maternal relatives had an over-abundance of any risk genotype, rather than one specific genotype. The observed genetic effects suggest an epigenetic mechanism in which decreased folate metabolism results in epigenetic alterations related to the increased rate of NTDs and other birth defects seen in the maternal lineage. Future studies on the etiology of NTDs and other birth defects could benefit from including multigenerational extended families, in order to explore potential epigenetic mechanisms.

neural tube defectsfolate metabolismDHFR 19bp deletionMTHFD1 1958G>AMTHFR 1298A>CMTHFR 677C>TSLC19A1 80A>Gmaternal inheritance
Introduction

Neural tube defects (NTDs) are congenital anomalies of the brain and spinal cord that cause significant morbidity and mortality (Sutton et al., 2008). Historically, Ireland has had high rates of NTDs (Coffey, 1983), with rates decreasing significantly in the last 50 years (Botto et al., 2006). Both environmental and genetic factors have been implicated in the etiology of NTDs. The evidence for environmental factors lies in the protective effect of periconceptional folic acid intake (Czeizel and Dudas, 1992). The evidence for genetic factors lies in the increased recurrence risk in siblings of individuals with NTDs, the increased risk for NTDs among other relatives, and the increased risk for birth defects overall among both close and distant relatives (MRC Vitamin Study Research Group, 1991; Deak et al., 2008). Maternal relatives have higher NTD and overall birth defect rates than paternal relatives, consistent with a combination of genetic and epigenetic mechanisms contributing to NTD and overall birth defects risk (Byrne et al., 1996; Byrne, 2008, 2010, 2011; Deak et al., 2008).

The genetic mechanisms associated with NTDs have been addressed through studies on associations between NTD risk and variants in genes involved in folate metabolism, chosen because of the protective effect of folic acid on NTDs (Molloy et al., 2009; Shaw et al., 2009). Methylene tetrahydrofolate reductase (MTHFR) plays an important role in folate metabolism. The 677 T allele results in reduced enzyme activity and decreased folate metabolism (Frosst et al., 1995). The MTHFR 677 TT genotype in infants and mothers has been associated with increased risk for NTDs [reviewed in Botto and Yang (2000)]. Studies in Irish populations found significantly more individuals with NTDs with the TT and CT genotypes and more parents of affected individuals with the TT genotype (Whitehead et al., 1995; Shields et al., 1999; Kirke et al., 2004; Carter et al., 2011). The MTHFR 1298A>C variant also results in reduced enzyme activity, although its role in NTD risk is less clear (van der Put et al., 1998; Botto and Yang, 2000; Wang et al., 2012a). Parle-McDermott et al. (2003) did not find an association between the MTHFR 1298A>C variant and NTDs in an Irish population. Methylene tetrahydrofolate dehydrogenase 1 (MTHFD1) acts in folate metabolism and DNA synthesis, and the 1958G>A variant results in decreased enzyme stability (Christensen et al., 2009). Studies in Irish populations (Brody et al., 2002; Parle-McDermott et al., 2006) identified a significant excess of MTHFD1 AA homozygotes among mothers of NTD cases compared to controls. A recent larger scale study in an Irish population also identified the MTHFD1 A allele as a risk factor for NTDs (Pangilinan et al., 2012). Dihydrofolate reductase (DHFR) also acts in folate metabolism. Studies have yielded conflicting results on the effect of the 19 bp deletion at intron 1 on red blood cell folate concentrations (Kalmbach et al., 2008; Stanislawska-Sachadyn et al., 2008). Researchers from Trinity College in Dublin observed a slight protective effect of the homozygous DHFR deletion genotype in mothers (Parle-McDermott et al., 2007). SLC19A1 plays an important role in folate transport, and some previous studies have reported associations between the 80A>G variant and NTD risk (Chango et al., 2000; De Marco et al., 2003; Wang et al., 2012b). O'Leary et al. (2006) did not find a difference in the frequency of the SLC19A1 80A>G reduced folate carrier in mothers, fathers and affected individuals compared with controls in an Irish population.

Epigenetic mechanisms related to NTDs have not been well explored. The higher rates of NTDs in maternal relatives compared with paternal relatives suggest the presence of epigenetic effects. Folate metabolism genes are of interest because folate metabolism is critical to the production of one-carbon subunits needed for epigenetic modification (Crider et al., 2012). Decreased folate metabolism, seen with certain gene variants, might hinder proper establishment or maintenance of epigenetic modifications. The epigenetic changes of interest would occur in the maternal genome in utero, that is, during the grandmother's pregnancy (Crider et al., 2012). Thus, in contrast to genetic mechanisms, in which the mothers' and affected individuals' genotypes would be most relevant, for epigenetic mechanisms, the maternal and paternal grandparents' genotypes would be important. In the absence of information on grandparents' genotypes, aunts, uncles, and first cousins provide an indirect measure of the grandparents' genotypes.

In this hypothesis-generating study, we investigated possible mechanisms for the excess risk for NTDs and birth defects overall in the maternal lineage of Irish individuals with NTDs. To do this, we evaluated five polymorphisms in genes involved in folate metabolism in the maternal vs. paternal lineage, both individually and with low folate-metabolizing genotypes grouped together. To our knowledge, this study is the first to look for genomic associations in extended relatives in NTD families and the effect of maternal vs. paternal lineage risk. We chose five polymorphic markers for evaluation, which were among the most intensely studied so far: the 19 bp deletion in DHFR, MTHFD1 1958G>A, MTHFR 1298A>C, MTHFR 677C>T, and SLC19A1 80A>G.

Materials and methodsStudy design and participants

The participants in this cross-sectional study were members of Irish families with NTDs. These families had previously participated in epidemiologic studies of the Boyne Research Institute [for instance, (Byrne, 2011)]. The families were identified initially through their membership in the Louth-Meath (northeastern region of the Republic of Ireland) branch of the Irish Association for Spina Bifida and Hydrocephalus, and the Ballymena (Northern Ireland) branch of the (UK) Association for Spina Bifida and Hydrocephalus. Families were also recruited through word of mouth and through radio announcements in the Louth-Meath area.

To be eligible for the study, male and female participants had to be 15 years of age or older, reside in the northeast region of the Republic of Ireland or Northern Ireland, and be either the individual with an NTD (only one NTD in sibship) or a blood relative of the affected individual. In terms of relationship, relatives (defined according to their relationship to the affected individual) had to be one of the following: parents, siblings, nieces, nephews, offspring (of the affected individual), grandparents, uncles, aunts, first cousins and first cousins once removed (the offspring of the first cousins).

Family members were interviewed in phases. Phase I evaluated pregnancy histories and the health of siblings in the original nuclear families (1995–2002). Phase II subjects (uncles and aunts of the affected individuals) were identified by their siblings during Phase I and interviewed between 2000 and 2002. Phase III consisted of interviews with first cousins of the affected individuals (the children of the uncles and aunts) between 2002 and 2009. Between 2005 and 2006, affected individuals and their siblings were interviewed.

For the present study, conducted in 2007, the response rate among eligible, locatable relatives was 49.5%, or 331 participants. From these, 322 gave a blood sample. Of these with a blood sample, 34 were relatives in families in which the individual with an NTD had anencephaly, 41 in which the affected individual had spina bifida occulta, 240 in which the affected individual had open spina bifida and 7 in which the affected individual had encephalocele.

This study was conducted according to the guidelines from the Declaration of Helsinki and all procedures involving human subjects were approved by the Ethics Board of the Boyne Research Institute and the Ethics Committee of the Health Service Executive of the Irish Government's Department of Health and Children. Written informed consent was obtained from all subjects.

Statistical analysis

Data were managed using Microsoft EXCEL©, and statistical analysis was conducted using SAS 9.1 (SAS Institute, Cary, NC, USA), STATA©, and IBM® SPSS® v19. Measures of association were evaluated with t-tests for continuous variables and chi-square and Fisher's exact tests for categorical variables, with the level of statistical significance set at 0.05. First-degree relatives were parents, siblings, and children; second-degree relatives were aunts, uncle, nieces, nephews, and grandparents; and third-degree and higher relatives were cousins and first cousins once removed. Genotypes were analyzed using an additive model, consistent with previous studies. For logistic regression analyses, genotypes were dichotomized using a recessive model, consistent with previous studies. For each individual, we calculated the number of “risk genotypes” present. “Risk genotypes” were those associated with lower folate metabolism in previous studies (DHFR WT/WT, MTHFD1 1958 AA, MTHFR 1298 CC, MTHFR 677 TT, and SLC19A1 80 GG).

DNA methods

DNA was extracted from blood clots using Gentra Puregene Blood kit and Clotspin Baskets per manufacturer's protocol (Qiagen, Valencia, CA). Quantitation was performed by spectrophotometric reading at 260nm and all samples normalized to 10 ng/μl. For PCR reactions, 2 μl genomic DNA (20 ng) was used. Pyrosequencing technology (Qiagen, Valencia, CA) was used to genotype DNA samples for DHFR 19 bp deletion (rs70991108), MTHFD1 1958G>A (rs2236225), MTHFR 1298A>C (rs1801131), MTHFR 677C>T (rs1801133), and SLC19A1 80A>G (rs1051266). PCR and pyrosequencing primers (Table 1) were designed using PSQ Assay Design Software Version 1.0.6. After amplification, the pyrosequencing primer was used to extend the PCR product by addition of nucleotides in a specific sequential order based on the known sequence. The resulting pyrogram of signal peaks was used to determine the genotype. PCR amplification for the DHFR 19 bp deletion was performed using the GC Rich PCR System (Roche, Indianapolis, IN) as previously described (Parle-McDermott et al., 2007). PCR for all other variants was performed in a 20 μl reaction mix containing 0.5 μM each forward and reverse primer and Applied Biosystems GeneAmp® Fast PCR Master Mix. Amplification was performed using the 9800 Fast PCR System according to manufacturer's suggested protocol (Applied Biosystems, Foster City, CA). Pyrosequencing was carried out as described by manufacturer's protocol (Qiagen, Valencia, CA).

Primer sequences and PCR amplicon size for DHFR, MTHFD1, MTHFR, and SLC19A1 variant analysis.

VariantForward primer (5′–3′)aReverse primer (5′–3′)aSequencing primer (5′–3′)bAmplicon size (bp)
DHFR 19 bp deletion (rs70991108)CAAGAACGGGGACCTGCC*ATCCTCTCGCCGGGAGTCGCCCACGGTCGGGGT114,133
MTHFD1 1958G>A (rs2236225)TCATTTTGCTTAGAGCACAGTAGA*CATCCACTTTCCTGGCTTAAATCACAGTAGAGAGTGCCAA63
MTHFR 1298A>C (rs1801131)GGAGCTGCTGAAGATGTGG*TGGTTCTCCCGAGAGGTAAAGAGGAGCTGACCAGTGA85
MTHFR 677C>T (rs1801133)*ATTGGCAGGTTACCCCAAAGGTCGGTGCATGCCTTCACAACGTGATGATGAAATCG156
SLC19A1 80A>G (rs1051266)*GTGGAACCTGGGCCTGACGTAGGGGGTGATGAAGCTCTCCAAAGGTAGCACACGA114

Biotinylated at 5′ end.

Represent forward primers for DHFR, MTHFD1, MTHFR 1298, and reverse primers for MTHFR 677 and SLC19A1.

cMTHFR 677 and SLC19A1 assays on reverse complement strand.

Quality control

Five percent of the samples were selected at random for repeat testing of the five polymorphisms, and the concordance was checked. We also tested each variant for any significant departure from Hardy-Weinberg equilibrium (HWE) by a chi-square goodness-of-fit test to verify data quality. In addition, Mendelian inheritance for all variants was assessed in 13 families with full trios using the PedCheck program (O'Connell and Weeks, 1998). Water and genotype controls (Coriell Cell Repositories, Camden, NJ) were included on each assay plate. The tests for HWE indicated that none of the variants deviated significantly from HWE, supporting high quality genotype data. All replica test results were 100% concordant.

Results

We obtained blood samples from all categories of relatives, including individuals with NTDs (n = 13), offspring (N = 2), fathers (N = 29), mothers (N = 37), siblings (N = 59), nieces and nephews (N = 2), uncles and aunts (N = 66), first cousins (N = 105), and grandparents (N = 2) from 47 families. Paternal relatives (N = 65) and maternal relatives (N = 115) comprised uncles and aunts, first cousins, first cousins once removed, and grandparents, who were related to the affected individual through the father and mother, respectively. The average age of the study participants was 41 years. The overall distribution of each of the five polymorphisms is shown in Table 2. The genotype distributions for the five polymorphisms in our study were similar to those from other studies, with a few exceptions (Table 3).

Characteristics of the participants, northeast region of the Republic of Ireland and Northern Ireland, 1995–2009.

RangeMedian
Age in yearsOverall16–8641.5
Affected individual18–4728.8
1st degree17–7344.7
2nd degree18–8652.1
3rd degree or higher16–5933.5
N%
GenderMales12037.3
Females20262.7
RelationshipAffected individual134.0
Offspring (of affected individual)20.6
Father299.0
Mother3711.5
Sibling5918.3
Niece/Nephew20.6
Uncle/aunt6620.5
First cousin10532.6
First cousin once removed72.2
Grandfather10.3
Grandmother10.3
1st degree12739.4
2nd degree7021.7
3rd degree or higher11234.8
Genotyping resultsDHFR 19 bp deletion
del/del6119.0
WT/del17253.4
WT/WT8927.6
MTHFD1 1958G>A
AA5115.8
AG15648.5
GG11535.7
MTHFR 1298A>C
AA14545.0
AC14745.7
CC309.3
MTHFR 677C>T
CC13441.6
CT15648.5
TT329.9
SLC19A1 80A>G
AA5818.0
AG16852.2
GG9629.8
LinePaternal relatives6536.0
Maternal relatives11564.0

Distribution of genotypes and comparison with other studies.

DHFRDel/DelDel/WTWT/WTP
N%N%N%
This studyIreland, NTD relatives6118.917253.58927.6Ref.
Parle-McDermott et al., 2007Ireland, controls4618.612851.87329.60.88
Kalmbach et al., 2008US, unaffected21617.164651.439631.50.39
This studyIreland, mothers718.92259.5821.6Ref.
Parle-McDermott et al., 2007Ireland, mothers4918.013047.09635.00.25
MTHFD1A/AA/GG/G
This studyIreland, NTD relatives5115.715648.611535.7Ref.
HapMap (CEU)Utah2017.75447.83934.50.90
Parle-McDermott et al., 2006Ireland, controls15019.541153.420927.10.02
This studyIreland, Mothers821.61746.01232.4Ref.
Parle-McDermott et al., 2006Ireland, mothers6527.011848.06225.00.62
MTHFR 1298A/AA/CC/C
This studyIreland, NTD relatives14545.014745.7309.3Ref.
HapMap (CEU)Utah4943.45145.11311.50.79
van der Put et al., 1998Holland, controls17944.418646.2389.40.99
Mills et al., 2008Ireland, controls51949.443941.8928.80.38
MTHFR 677C/CC/TT/T
This studyIreland, NTD relatives13441.615648.6329.8Ref.
HapMap (CEU)Utah5347.05044.2108.80.62
Shields et al., 1999Ireland, controls11447.110844.6208.30.41
Mynett-Johnson et al., 2002Ireland, unaffected volunteers20934.729949.79415.60.02
This studyIreland, affected individuals215.4969.2215.4Ref.
Shields et al., 1999Ireland, affected individuals10137.311943.95118.80.18
This studyIreland, mothers1746.01848.625.4Ref.
Shields et al., 1999Ireland, mothers8036.710849.53013.80.29
SLC19A1A/AA/GG/G
This studyIreland, NTD relatives5818.016852.29629.8Ref.
HapMap (CEU)Utah1816.06255.43228.60.83
O'Leary et al., 2006Ireland, controls10327.017947.09926.00.02
This studyIreland, affected individuals430.8753.9215.3Ref.
O'Leary et al., 2006Ireland, affected individuals6123.013150.07227.00.60
This studyIreland, mothers718.91951.41129.7Ref.
O'Leary et al., 2006Ireland, mothers6525.011945.08030.00.70
This studyIreland, fathers724.11655.2620.7Ref.
O'Leary et al., 2006Ireland, fathers5822.012748.07930.00.58
Analysis of individual polymorphisms

The distributions of each genotype and allele by relationship to the individual with an NTD, degree of relationship to the affected individual, and maternal vs. paternal lineage are shown in Table 4. We did not observe any significant associations between individual genotypes and the affected individuals or mothers, although the higher folate-metabolizing MTHFR 677 CC genotype tended to be found less often in affected individuals (15.4%) than in all other relatives (42.7%) (Fisher's exact test, p = 0.08) (Table 4 and data not shown). We also did not see any significant associations between individual genotypes and maternal vs. paternal lineage (Table 4).

Distribution of genotypes and alleles according to relationship within NTD-affected families and lineage, northeast region of the Republic of Ireland and Northern Ireland, 1995–2009.

DHFRGenotype frequenciesPbAllele frequenciesPb
Del/DelDel/WTWT/WTDelWT
N%N%N%N%N%
Affected individual215.4646.2538.40.741038.51661.50.45
Mother718.92259.5821.60.663648.63851.40.58
Father620.71551.7827.60.972746.63153.40.89
Sibling1322.03559.31118.70.236151.75748.30.13
Uncle/aunt1624.23045.52030.30.306247.07053.00.71
First cousin1312.46158.13129.50.118741.412358.60.13
Other rel.430.8323.0646.20.061142.31557.70.73
1st degree relatives2620.57256.72922.80.3612448.813051.10.23
2nd degree relatives1724.33347.12028.60.366747.97352.10.61
3rd degree or higher relatives1614.36154.53531.20.209341.513158.50.10
Paternal linea913.83756.91929.30.405542.37557.70.65
Maternal linea2420.95547.83631.310344.812755.2
MTHFD1G/GG/AA/AGA
Affected individual215.4861.5323.10.241246.21453.80.14
Mother1232.41746.0821.60.594155.53344.60.40
Father1448.31034.5517.20.263865.52034.50.36
Sibling1627.23254.21118.60.316454.25445.80.16
Uncle/aunt3451.52842.446.10.0039672.73627.30.0008
First cousin3432.45350.51817.10.6812157.68942.40.40
Other rel.323.1861.5215.40.551453.81246.20.52
1st degree relatives4333.96047.22418.90.3614657.510842.50.20
2nd degree relatives3550.03042.957.10.0110071.44028.60.003
3rd degree or higher relatives3531.25851.81917.00.3412855.49642.90.20
Paternal line2335.43655.469.20.248263.14836.90.93
Maternal line4740.95043.51815.614462.68637.4
MTHFR 1298A/AA/CC/CAC
Affected individual646.2646.217.61.001869.2830.80.88
Mother1540.51746.0513.50.614763.52736.50.40
Father1137.91551.8310.30.723763.82136.20.49
Sibling1830.53559.3610.20.047160.24739.80.05
Uncle/aunt3350.02639.4710.60.529269.74030.30.61
First cousin5451.44441.976.70.2215272.45827.60.09
Other rel.861.5430.817.70.492076.9623.10.31
1st degree relatives4434.66853.61511.80.00915661.49838.60.005
2nd degree relatives3550.02840.0710.00.569870.04230.00.53
3rd degree or higher relatives6053.64540.276.20.0516573.75926.30.02
Paternal line3046.32944.569.20.538968.54131.50.27
Maternal line6354.84438.287.017073.96026.1
MTHFR 677C/CC/TT/TCT
Affected individual215.4969.2215.40.091350.01350.00.08
Mother1746.01848.625.40.595270.32229.70.39
Father1137.91551.7310.40.913763.82136.20.73
Sibling2237.33457.635.10.197866.14033.90.95
Uncle/aunt3248.52842.469.10.449269.74030.30.29
First cousin4441.94744.81413.30.3313564.37535.70.56
Genotype frequenciesPbAllele frequenciesPb
MTHFR 677C/CC/TT/TCT
N%N%N%N%N%
Other rel.646.2538.5215.30.601765.4934.60.96
1st degree relatives5240.96752.886.30.1417167.38332.70.72
2nd degree relatives3448.63042.968.50.539870.04230.00.32
3rd degree or higher relatives4641.15044.61614.30.1214263.48236.60.22
Paternal line3350.82640.069.20.389270.83829.20.16
Maternal line4740.95245.21613.914663.58436.5
SLC19A1A/AA/GG/GAG
Affected individual430.8753.9215.30.321557.71142.30.15
Mother718.91951.41129.70.993344.64155.40.93
Father724.11655.2620.70.453051.72848.30.22
Sibling1118.63559.41322.00.345748.36151.70.31
Uncle/aunt1116.72740.92842.40.044937.18362.90.07
First cousin1615.25653.43331.40.668841.912258.10.44
Other rel.215.4861.5323.10.931246.21453.80.83
1st degree relatives2519.77256.73023.60.0912248.013252.00.06
2nd degree relatives1115.73042.92941.40.075237.18862.90.08
3rd degree or higher relatives1816.05952.73531.30.889542.412957.60.67
Paternal line1218.53452.31929.20.405844.67255.40.21
Maternal line1714.85346.14539.18737.814362.2

Paternal line consists of uncles and aunts, grandparents, first cousins and first cousins once removed who are related to the affected individual through the father; maternal line consists of uncles and aunts, grandparents, first cousins and first cousins once removed who are related to the affected individual through the mother.

For specific relatives, each category was compared with all other relatives to calculate P-values. For example, fathers were compared with all other relatives to calculate the P-value. For degree of relationship to the affected individual, each category was compared with a combination of the other two to calculate the P-value. For example, 1st degree relatives were compared with 2nd and 3rd degree relatives combined to calculate the P-value. P-values were calculated using a chi-square test unless one or more cells contained less than 5, in which case Fisher’s exact test was used.

Relatives of individuals with NTDs might be expected to have genotypes associated with increased or decreased risk for NTDs, and the prevalence of these genotypes could increase or decrease with degree of closeness to the affected individual. For MTHFR 1298A>C, a lower percentage of first-degree relatives (34.6%) had the higher folate-metabolizing AA genotype compared with second-degree (50.0%) relatives and third-degree and higher (53.6%) relatives (p = 0.009). Similarly, the A allele was less common in first-degree relatives (61.4%) than in second-degree relatives (70.0%) and third-degree and higher (73.7%) relatives (p = 0.005). When the genotype (AA vs. A/C, C/C) and relatives (first-degree relatives vs. all other relatives) were dichotomized, the odds ratio for the association of the AA genotype with first-degree relatives was 0.49 (0.30–0.78, p = 0.002) (data not shown). Consistent with this, a lower percentage of siblings, 30.5%, had the AA genotype compared to 48.3% among other relatives (p = 0.04, Table 4). We did not observe a consistent trend of genotypes being more or less common in relatives more closely related to the affected individual for any of the other genotypes (Table 4).

We observed some associations with individual categories of relatives. For example, fewer first cousins had the higher folate-metabolizing DHFR del/del genotype compared with all other relatives (12.4 vs. 22.1%, p = 0.11, Table 4 and data not shown), the higher folate-metabolizing MTHFD1 GG genotype was higher in uncles and aunts compared with all other relatives (51.5 vs. 31.6%, p = 0.003) (Table 4), and the lower folate-metabolizing SLC19A1 GG genotype was more common in uncles and aunts than in other relatives (42.4 vs. 26.6%, p = 0.04) (Table 4). However, the meaning of these isolated associations is unclear.

Number of risk genotypes

While the frequencies of individual genotypes did not differ significantly between maternal and paternal relatives (Table 4), the number of risk genotypes, defined as the number of homozygous genotypes associated with lower folate metabolism, did show significant variation between the maternal and paternal lines (p = 0.02) (Table 5). A lower percentage of maternal relatives (20.9%) had no (0) risk genotypes compared with paternal relatives (41.5%), while more maternal relatives (54.8%) than paternal relatives (35.4%) had 1 risk genotype. When the risk genotypes were dichotomized (0, 1–3), the odds ratio for association with the maternal line compared with the paternal line was 2.69 (1.38–5.25, p = 0.004).

Number of risk genotypes by relationship to affected individuals with NTDs and lineage, northeast region of the Republic of Ireland and Northern Ireland, 1995–2009.

Number of risk genotypes/low metabolizer genotypesP-valuea (Fisher’s exact test)
0123
N%N%N%N%
Maternal line2420.96354.82420.943.40.02
Paternal line2741.52335.41218.534.6
Affected individual or mother2244.01122.01530.024.00.004
All other relatives8631.613047.84717.393.3
Mother1746.0718.91232.412.70.005
All other relatives9131.913447.15017.5103.5

P-values compare maternal vs. paternal line, affected individual or mother vs. all other relatives, and mother vs. all other relatives.

The numbers of risk genotypes in affected individuals and mothers (p = 0.004) or mothers alone (p = 0.005), compared with all other relatives were also statistically significant (Table 5). A higher percentage of affected individuals and mothers (30.0 vs. 17.3%) and mothers alone (32.4 vs. 17.5%) had 2 risk genotypes compared with all other relatives, while a lower percentage of affected individuals and mothers (22.0 vs. 47.8%) and mothers alone (18.9 vs. 47.1%) had 1 risk genotype. A higher percentage of affected individuals and mothers (44.0 vs. 31.6%) and mothers alone (46.0 vs. 31.9%) had 0 risk genotypes.

Discussion

In the Irish families in our study, we found that maternal relatives were more likely than paternal relatives to have one or more genotypes associated with decreased folate metabolism. As expected from previous studies, mothers and affected individuals had differential distributions of risk genotypes compared with all other relatives.

The observed increased number of affected relatives in the maternal lineage of families could be due to genetic and/or epigenetic mechanisms. A simple genetic explanation would be that more maternal relatives have alleles associated with increased risk for NTDs simply due to chance, so the likelihood of inheriting these alleles from the maternal side would be increased. With a simple genetic model, we would expect that relatives would share the same risk genotype as the affected individuals and/or mothers. However, we did not observe this, but instead found that maternal relatives had an over-abundance of any risk genotype, rather than one specific genotype. We found that more maternal relatives had one risk genotype, while more paternal relatives had no risk genotypes. However, if this were the main explanation behind the increased number of affected individuals on the maternal side, we would expect that relatives more closely related to the affected individuals would share more genetic risk factors, and thus the prevalence of higher risk genotypes would be greater in first-degree relatives than in more distant ones. However, with the exception of MTHFR 1298, we did not observe consistent patterns of increasing numbers of risk genotypes with closer degree of relationship to the affected individual.

The mother's genotype during pregnancy plays a significant role in NTD risk, best studied in MTHFR (Yan et al., 2012), and the inheritance of these alleles would necessarily come from the maternal line. In this case, mothers and possibly affected individuals might be expected to have an increased prevalence of genotypes associated with greater NTD risk. In our study, fewer affected individuals had the higher folate-metabolizing MTHFR 677 CC genotype, consistent with earlier reports that the homozygous TT genotype is a risk factor for NTDs (Shields et al., 1999; Botto and Yang, 2000). Likewise, affected individuals and mothers, as well as mothers only, were more likely to have 2 risk genotypes. Our study of these associations was limited by our small sample size and was not the main objective of this study. While our findings were consistent with previous literature, these associations alone were not strong enough to explain the increased number of affected relatives on the mother's side, suggesting that other mechanisms are also at work.

Whenever there is preferential inheritance of a condition from the maternal vs. paternal line, epigenetic effects are strongly suspected. The genes included in our study are not located in known imprinted regions. However, folate metabolism is critical to the production of one-carbon subunits needed for epigenetic modification. Certain variants of folate metabolism genes (possibly including those in our study), which are associated with NTD risk, might affect proper imprinting and establishment or maintenance of other epigenetic changes. These epigenetic effects could occur in the maternal genome and would have been initiated when the mother was a fetus herself, but may also have taken place in the genomes of the mother's parents (Crider et al., 2012). If this mechanism were at work, the maternal and paternal grandparents' genotypes would be of importance. Unfortunately, we were unable to obtain samples from many grandparents. However, aunts, uncles, and first cousins provide an indirect measure of the grandparents' genotypes, and significant associations for these relatives were seen for DHFR, MTHFD1, and SLC19A1. The increased number of risk genotypes in maternal vs. paternal relatives is consistent with an epigenetic mechanism.

Our study is the first to look at the genotypes of second- and third-degree or more relatives in families affected by NTDs and to begin to explore the genetic and epigenetic mechanisms behind the increased NTD and overall birth defects risks in these families. However, our study is preliminary and has a number of limitations. The limitations of our study include a relatively small sample size, lack of triads (both parents and offspring) among distant relatives to allow for heritability testing, and inadequate number of samples from grandparents to evaluate epigenetic effects. Additionally, more maternal than paternal relatives participated in the study. However, while this might have affected the stability of the point estimate for the association of risk genotypes with the maternal vs. paternal line, it did not affect the magnitude or direction. The response rate was relatively low for our study, and the participation of families and of relatives within families could be biased and not reflect the true underlying distribution of families with NTDs in Ireland. However, selection for and participation in the study was independent of genotype and thus participants can be considered close to a random sample of relatives. The candidate gene approach used here evaluated a small number of polymorphisms and would have missed associations with other variants in folate metabolism pathway genes. An important limitation of our study (and others) is that we were not able to measure the mother's intake of folate and folic acid during pregnancy. The genes included in this study function in folate metabolism, so that their effect depends on the levels of folate present in the body. For example, high intake of folic acid and folate could mask genetic effects that would be seen at lower intake levels. Lastly, the number of multiple comparisons evaluated in our study likely produced some spurious findings.

Our results are consistent with both genetic and epigenetic mechanisms contributing to the increased NTD and overall birth defects risk in the maternal lineage. The diversity of genetic risk seen in these families suggests that etiology of NTDs is multi-faceted. Future studies will benefit from including other relatives, especially grandparents, in addition to mothers, fathers, and affected individuals, especially if epigenetic mechanisms are to be addressed.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

We are grateful to the families who chose to participate in our studies, to Rebecca Lawler, Sharon McGinty, and the volunteers who helped to conduct this study, to Claymon Laboratories for biochemical tests, to Stanimila Nikolova for technical assistance and to Mary E. Cogswell, Quanhe Yang, and Ramal Moonesinghe for their insightful comments. The Joseph E. and Marjorie B. Jones Foundation of Washington, DC provided partial support to Julianne Byrne. No additional external funding was received for this study.

Disclaimer

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.

ReferencesBottoL. D.LisiA.BowerC.CanfieldM. A.DattaniN.De ViganC. (2006). Trends of selected malformations in relation to folic acid recommendations and fortification: an international assessment. Birth Defects Res. A. Clin. Mol. Teratol. 76, 693705 10.1002/bdra.2030717029289BottoL. D.YangQ. (2000). 5,10-Methylenetetrahydrofolate reductase gene variants and congenital anomalies: a HuGE review. Am. J. Epidemiol. 151, 862877 10.1093/oxfordjournals.aje.a01029010791559BrodyL. C.ConleyM.CoxC.KirkeP. N.McKeeverM. P.MillsJ. L. (2002). A polymorphism, R653Q, in the trifunctional enzyme methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase/formyltetrahydrofolate synthetase is a maternal genetic risk factor for neural tube defects: report of the Birth Defects Research Group. Am. J. Hum. Genet. 71, 12071215 10.1086/34421312384833ByrneJ. (2008). Birth defects in uncles and aunts from Irish families with neural tube defects. Birth Defects Res. A. Clin. Mol. Teratol. 82, 815 10.1002/bdra.2040618044714ByrneJ. (2010). Birth defects among maternal first cousins in Irish families with a neural tube defect. Ir. J. Med. Sci. 179, 375380 10.1007/s11845-009-0381-x19618238ByrneJ. (2011). Three generations of matrilineal excess of birth defects in Irish families with neural tube defects. Ir. J. Med. Sci. 180, 6972 10.1007/s11845-010-0632-x21057889ByrneJ.CamaA.ReillyM.VigliaroloM.LevatoL.BoniL. (1996). Multigeneration maternal transmission in Italian families with neural tube defects. Am. J. Med. Genet. 66, 303310 10.1002/(SICI)1096-8628(19961218)66:3<303::AID-AJMG13>3.3.CO;2-28985492CarterT. C.PangilinanF.TroendleJ. F.MolloyA. M.VanderMeerJ.MitchellA. (2011). Evaluation of 64 candidate single nucleotide polymorphisms as risk factors for neural tube defects in a large Irish study population. Am. J. Med. Genet. A 155, 1421 10.1002/ajmg.a.3375521204206ChangoA.Emery-FillonN.de CourcyG. P.LambertD.PfisterM.RosenblattD. S. (2000). A polymorphism (80G–>A) in the reduced folate carrier gene and its associations with folate status and homocysteinemia. Mol. Genet. Metab. 70, 310315 10.1006/mgme.2000.303410993718ChristensenK. E.RohlicekC. V.AndelfingerG. U.MichaudJ.BigrasJ. L.RichterA. (2009). The MTHFD1 p.Arg653Gln variant alters enzyme function and increases risk for congenital heart defects. Hum. Mutat. 30, 212220 10.1002/humu.2083018767138CoffeyV. P. (1983). Neural tube defects in Dublin 1953–1954 and 1961–1982. Ir. Med. J. 76, 411413 6642965CriderK. S.YangT. P.BerryR. J.BaileyL. B. (2012). Folate and DNA methylation: a review of molecular mechanisms and the evidence for folate's role. Adv. Nutr. 3, 2138 10.3945/an.111.00099222332098CzeizelA. E.DudasI. (1992). Prevention of the first occurrence of neural-tube defects by periconceptional vitamin supplementation. N. Engl. J. Med. 327, 18321835 10.1056/NEJM1992122432726021307234DeakK. L.SiegelD. G.GeorgeT. M.GregoryS.Ashley-KochA.SpeerM. C. (2008). Further evidence for a maternal genetic effect and a sex-influenced effect contributing to risk for human neural tube defects. Birth Defects Res. A. Clin. Mol. Teratol. 82, 662669 10.1002/bdra.2051118937341De MarcoP.CalevoM. G.MoroniA.MerelloE.RasoA.FinnellR. H. (2003). Reduced folate carrier polymorphism (80A–>G) and neural tube defects. Eur. J. Hum. Genet. 11, 245252 10.1038/sj.ejhg.520094612673279FrosstP.BlomH. J.MilosR.GoyetteP.SheppardC. A.MatthewsR. G. (1995). A candidate genetic risk factor for vascular disease: a common mutation in methylenetetrahydrofolate reductase. Nat. Genet. 10, 111113 10.1038/ng0595-1117647779KalmbachR. D.ChoumenkovitchS. F.TroenA. P.JacquesP. F.D'AgostinoR.SelhubJ. (2008). A 19-base pair deletion polymorphism in dihydrofolate reductase is associated with increased unmetabolized folic acid in plasma and decreased red blood cell folate. J. Nutr. 138, 23232327 10.3945/jn.108.09640419022952KirkeP. N.MillsJ. L.MolloyA. M.BrodyL. C.O'LearyV. B.DalyL. (2004). Impact of the MTHFR C677T polymorphism on risk of neural tube defects: case-control study. BMJ 328, 15351536 10.1136/bmj.38036.646030.EE15155469MillsJ. L.MolloyA. M.Parle-McDermottA.TroendleJ. F.BrodyL. C.ConleyM. R. (2008). Folate-related gene polymorphisms as risk factors for cleft lip and cleft palate. Birth Defects Res. A. Clin. Mol. Teratol. 82, 636643 10.1002/bdra.2049118661527MolloyA. M.BrodyL. C.MillsJ. L.ScottJ. M.KirkeP. N. (2009). The search for genetic polymorphisms in the homocysteine/folate pathway that contribute to the etiology of human neural tube defects. Birth Defects Res. A. Clin. Mol. Teratol. 85, 285294 10.1002/bdra.2056619235830MRC Vitamin Study Research Group. (1991). Prevention of neural tube defects: results of the Medical Research Council Vitamin Study. Lancet 338, 131137 10.1016/0140-6736(91)90133-A1677062Mynett-JohnsonL. A.KeenanC.BlackI. L.LivingstoneW. J.LawlerM.RocheH. M. (2002). Thermolabile methylenetetrahydrofolate reductase (C677T): frequency in the Irish population. Ir. J. Med. Sci. 171, 3739 10.1007/BF0316894011993593O'ConnellJ. R.WeeksD. E. (1998). PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am. J. Hum. Genet. 63, 259266 10.1086/3019049634505O'LearyV. B.PangilinanF.CoxC.Parle-McDermottA.ConleyM.MolloyA. M. (2006). Reduced folate carrier polymorphisms and neural tube defect risk. Mol. Genet. Metab. 87, 364369 10.1016/j.ymgme.2005.09.02416343969PangilinanF.MolloyA. M.MillsJ. L.TroendleJ. F.Parle-McDermottA.SignoreC. (2012). Evaluation of common genetic variants in 82 candidate genes as risk factors for neural tube defects. BMC Med. Genet. 13:62 10.1186/1471-2350-13-6222856873Parle-McDermottA.KirkeP. N.MillsJ. L.MolloyA. M.CoxC.O'LearyV. B. (2006). Confirmation of the R653Q polymorphism of the trifunctional C1-synthase enzyme as a maternal risk for neural tube defects in the Irish population. Eur. J. Hum. Genet. 14, 768772 10.1038/sj.ejhg.520160316552426Parle-McDermottA.MillsJ. L.KirkeP. N.O'LearyV. B.SwansonD. A.PangilinanF. (2003). Analysis of the MTHFR 1298A–>C and 677C–>T polymorphisms as risk factors for neural tube defects. J. Hum. Genet. 48, 190193 10.1007/s10038-003-0008-412730722Parle-McDermottA.PangilinanF.MillsJ. L.KirkeP. N.GibneyE. R.TroendleJ. (2007). The 19-bp deletion polymorphism in intron-1 of dihydrofolate reductase (DHFR) may decrease rather than increase risk for spina bifida in the Irish population. Am. J. Med. Genet. A 143A, 11741180 10.1002/ajmg.a.3172517486595ShawG. M.LuW.ZhuH.YangW.BriggsF. B.CarmichaelS. L. (2009). 118 SNPs of folate-related genes and risks of spina bifida and conotruncal heart defects. BMC Med. Genet. 10:49 10.1186/1471-2350-10-4919493349ShieldsD. C.KirkeP. N.MillsJ. L.RamsbottomD.MolloyA. M.BurkeH. (1999). The “thermolabile” variant of methylenetetrahydrofolate reductase and neural tube defects: an evaluation of genetic risk and the relative importance of the genotypes of the embryo and the mother. Am. J. Hum. Genet. 64, 10451055 10.1086/30231010090889Stanislawska-SachadynA.BrownK. S.MitchellL. E.WoodsideJ. V.YoungI. S.ScottJ. M. (2008). An insertion/deletion polymorphism of the dihydrofolate reductase (DHFR) gene is associated with serum and red blood cell folate concentrations in women. Hum. Genet. 123, 289295 10.1007/s00439-008-0475-y18247058SuttonM.DalyL. E.KirkeP. N. (2008). Survival and disability in a cohort of neural tube defect births in Dublin, Ireland. Birth Defects Res. A. Clin. Mol. Teratol. 82, 701709 10.1002/bdra.2049818803309van der PutN. M.GabreelsF.StevensE. M.SmeitinkJ. A.TrijbelsF. J.EskesT. K. (1998). A second common mutation in the methylenetetrahydrofolate reductase gene: an additional risk factor for neural-tube defects? Am. J. Hum. Genet. 62, 10441051 10.1086/3018259545395WangX. W.LuoY. L.WangW.ZhangY.ChenQ.ChengY. L. (2012a). Association between MTHFR A1298C polymorphism and neural tube defect susceptibility: a metaanalysis. Am. J. Obstet. Gynecol. 206, 251e1–7. 10.1016/j.ajog.2011.12.02122265089WangH. G.WangJ. L.ZhangJ.ZhaoL. X.ZhaiG. X.XiangY. Z. (2012b). Reduced folate carrier A80G polymorphism and susceptibility to neural tube defects: a meta-analysis. Gene 510, 180184 10.1016/j.gene.2012.02.02022975209WhiteheadA. S.GallagherP.MillsJ. L.KirkeP. N.BurkeH.MolloyA. M. (1995). A genetic defect in 5,10 methylenetetrahydrofolate reductase in neural tube defects. QJM 88, 763766 8542260YanL.ZhaoL.LongY.ZouP.JiG.GuA. (2012). Association of the maternal MTHFR C677T polymorphism with susceptibility to neural tube defects in offsprings: evidence from 25 case-control studies. PLoS ONE 7:e41689 10.1371/journal.pone.004168923056169