Received date: December 21, 2015; Accepted date: April 05, 2016; Published date: April 10, 2016
Citation: Rai A, Jaiswal SK, Sukla KK, Mishra SK, Lakhotia AR, et al. (2016) Association of Genetic Polymorphisms in Genes Involved at the Branch Point of Nucleotide Biosynthesis and Remethylation with Down Syndrome Birth Risk: A Case-Control Study. J Mol Genet Med 10:207. doi:10.4172/1747-0862.1000207
Copyright: © 2016 Rai A, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Molecular and Genetic Medicine
DNA methylation and nucleic acid biosynthesis are two crucial phenomena for normal chromosomes segregation. From our earlier studies, MTHFR 677T individually and in combination with other gene polymorphisms, micronutrient deficiency and hyperhomocysteinemia was shown to be associated with risk in Down syndrome (DS) mothers. Remethylation and nucleic acid biosynthesis pathways are dependent on the activity of Methylenetetrahydrofolate reductase (MTHFR) and Thymidylate synthase (TYMS) respectively, competing for common substrate molecule 5,10-methelenetetrahydrofolate (5,10-MTHF). Role of MTHFD1 1958 G>A (affecting synthesis of 5,10-MTHF), MTHFR 677 C>T (affecting methylation), TYMS 5'UTR 28 bp repeat polymorphism and TYMS 3'UTR 6bp deletion polymorphism (affecting nucleic acid biosynthesis) in a cohort of 200 case mothers and 187 control mothers (also 146 case triads: mother, father, and child) were studied. We observed a significant association of MTHFR 677 C>T in a co-dominant model (p=0.0428) and dominant model (0.0194) as well as TYMS 5'UTR 28 bp repeat polymorphism in a recessive model (p=0.0005) and dominant model (0.0161). Genetic combination analysis revealed a significant additive effect of certain genotypic combinations (especially combination of MTHFR 677T and TYMS 2R alleles with other alleles or genotypes) in increasing risk. Weak linkage disequilibrium (LD) was observed between TYMS 5' and 3' UTR regions polymorphism in LD analysis. Transmission disequilibrium test (TDT) analysis revealed a consistent trend of preferential allele's transmission from parents. We concluded that genetic interaction of remethylation pathway and the nucleic acid metabolic pathway was significantly associated with risk factors for DS childbirth. However, replication studies are required to validate our observation in the population.
Chromosomal basis Down syndrome (DS) due to abnormal segregation of chromosome 21 has been known for many years. The mechanism behind mal-segregation of the chromosome is still unknown. In 90% of the cases, the nondisjunction occurs during maternal meiosis l. Advanced maternal age during conception and reduced recombination were the only well-established known factor for mal-segregation of chromosome 21 . James et al. hypothesized that maternal genetic polymorphism (MTHFR 677T) may cause hypomethylation in the pericentromeric region the chromosome, thereby leading to nondisjunction of chromosome 21 and further predicted that impaired folate metabolism might act as a maternal risk factor for having a DS child .
However, several comparable studies done on populations from different regions of the world have yielded conflicting results; in developing countries, MTHFR 677T was a significant genetic predisposing factor, while in studies from Europe and the Americas it was not an associated risk factor . Similarly, studies on other folate homocysteine pathway genes [4-6] have shown a much heterogeneity, like that of MTHFR 677T in association with Down syndrome.
Folate metabolic pathway is involved in two essential physiological processes: nucleic acid biosynthesis and DNA repair; and remethylation of homocysteine to form methionine. Methylenetetrahydrofolate dehydrogenase (MTHFD1) is a trifunctional gene having methylene tetrahydrofolate dehydrogenase, methenyltetrahydrofolate cyclohydrolase and formyl tetrahydrofolate synthetase activity and is involved in interconversion of intermediates molecules tetrahydrofolates for consistently providing 5,10- methlenetetrahydrofolate. It mediates the interconversion of 5,10- MTHF, 5,10-methenylTHF, and 10-formylTHF. The last two are donor cofactors for denovo purine and pyrimidine biosynthesis .
Also, one MTHFD1 gene function utilizes MTHF as substrate and competes with MTHFR gene function for the same substrate. Further, Thymidylate synthase (TYMS) competes with MTHFR for 5, 10- MTHF as a substrate. The former is involved in the conversion of dUMP to dTMP thus involved in nucleic acid biosynthesis, replication and repair and later is for remethylation of homocysteine to methionine to form S-adenosylmethionine (Universal methyl donar). Therefore, polymorphism affecting the activity of one of these two enzymes may shift pool of 5,10-MTHF towards either DNA synthesis or DNA methylation in the pathway. Studies also have shown that the decrease in activity of TYMS due to genetic polymorphism may be responsible for the inadequate formation of pyrimidine dTTP and thus impairing nucleic acid biosynthesis, DNA replication, DNA repair and chromosomal damage. Micronuclei formation, a marker for chromosomal damage has been shown significantly higher in younger mothers of DS children .
An earlier study on eastern Indian cohort has shown that the mothers having the MTHFR 677T allele (especially TT homozygotes) were susceptible to having DS children . In a recent study on eastern Indian population, a significant increased DNMT3B -579G/-149T allele frequency was observed in DS case mothers indicating an association of altered DNA methylation as a risk for DS .
Similarly, in a very recent study on Eastern Indian Cohort (on seven genetic polymorphisms SLC19A1 or RFC1 G80A, MTHFR C677T and A1298C, MTR A2756G, MTRR A66G, TCN2 C776G and CBS 844 in 68bp) we have shown that MTHFR 677T is an associated genetic risk factor and gene-gene interactions have a significant role in genetic predisposition in DS mothers . Hyperhomocysteinemia (hypHcy), vitamin-B12 and folate deficiency was significantly higher in the studied cohort of DS mothers. We have also shown from our previous findings that the incidence of hypHcy (30%), Vitamin B12 (50%) and folic acid (11%) deficiency in a healthy population of eastern Indian populace .
The impaired Folate-Hcy pathway as an associated risk in DS mothers from our previous and recent findings led us to investigate the role of four polymorphisms (MTHFD1 G1958A, MTHFR C677T, TYMS 5'UTR 28bp repeats and TYMS 3'UTR 6bpdel) in the predisposition of mothers for DS child birth. MTHFR 677T and TYMS 5'UTR 2R/2R independently and in combination with other genetic polymorphisms was significantly associated with risk in DS mothers.
The cases were referred to the Centre for Genetic Disorders (Genetic diagnostic unit) for karyotyping, upon clinical examination in the Pediatrics Department of the Institute of Medical Sciences, Banaras Hindu University. The period of sample recruitment was from 2010 to 2014. A cohort of 200 mothers of a Down syndrome child (MDS) confirmed by karyotyping (includes 146 DS triads) and 187 control mothers was included in the study. The mean age of the MDS was 27.0 ± 4.6 years at the time of the DS childbirth while the age of MDS at the time of sample collection was little less than 28.5 years.
The age of control mothers was 26.0 ± 4.2 years at the time of the child birth. They were collected immediately after delivery. The age of the 146 DS children and fathers included in the study was found to be 11.5 ± 3.5 months and 32.5 ± 4.2 years, respectively. The socioeconomic status of both the case and control mothers was similar. Both groups belong to a low economical status. All patients belonged to Eastern Uttar Pradesh and the contiguous regions of Bihar, Jharkhand and Chhattisgarh. A peripheral blood sample was drawn from mother, father, and child after obtaining a written consent from the parents. Control mothers were recruited from the Outpatient clinic of the Gynecology Department, Banaras Hindu University, Varanasi, India.
We paid extreme caution during sampling from control mothers who had atleast one healthy child and no history of any spontaneous abortions. Mothers having children with any chromosomal and/or congenital disorders were not included in the study as controls. This study was approved by the Institutional Ethical Committee of the Institute of Medical Sciences, Banaras Hindu University, Varanasi, India.
A peripheral blood samples (5 ml) was drawn from DS children, their parents, and control mothers under aseptic condition and in a heparinized syringe. Within a week of sample collection genomic DNA was extracted from the peripheral blood by high salt method, as described previously .
Genotyping for all the polymorphisms were performed using a polymerase chain reaction-restriction fragment length polymorphism method using gene specific primers [9,14,15]. The PCR product and restriction digested products were electrophoresed on agarose gel.
Allelic frequencies of the SNPs examined were tested for Hardy & Weinberg equilibrium using a Chi-square test. The genetic risk among case and control groups was compared and tested by calculating odds ratio with a 95% class interval. A two-tailed p-value of <0.05 was interpreted as statistically significant (using Graphpad software). For the transmission disequilibrium test, we calculated the Chi-square value with one degree of freedom and a p-value <0.05 was taken to be significant.
Haplotype frequencies and linkage disequilibrium (LD) were calculated using the Haploview software (Version 4.2) developed at The Broad Institute (http://www.broadinstitute.org/), which is based on the EM algorithm. The standardized disequilibrium coefficient (D) and correlation coefficient (r2) between these SNPs were calculated to investigate whether certain allelic combinations of the polymorphisms examined altered the genetic risk in MDS.
MTHFD1 G1958A polymorphism: MTHFD1 G1958A genotype and allele frequency distribution in mothers of DS child and control mothers are shown in Table 1. We didn't observe significantly altered distribution of alleles and genotype among case and control mothers in any of the three applied models: recessive (p=0.97), co-dominant (0.43) and dominant models (p=0.53).
|Gene and SNP||Genotype||Case (%) N=200||Control (%) N=187||OR(95% CI)||p value|
|MTHFD1 G1958A (rs2236225)||GG||25(12.80)||30(16.04)||1.00||Referent|
|GA||104(53.30)||94(50.27)||1.33 (0.73 -2.42)||0.4380|
|AA||66(33.90)||63(33.69)||1.457 (0.67 - 2.37)||0.5837|
|Recessive model||-||-||1.01 (0.66 - 1.54)||0.9742|
|Dominant model||-||-||1.25 (0.71 - 2.23)||0.5310|
|A allele frequency||0.40||0.41||1.07(0.80 - 1.43)||0.6876|
|MTHFR C677T (rs1801133)||CC||122(61)||136(72.73)||1.00||Referent|
|CT||72(36)||50(26.74)||1.61 (1.04 - 2.48)||0.0428|
|TT||6(3)||1(0.53)||6.69 (0.80 - 56.37)||0.1043|
|Recessive model||-||-||5.753 (0.69 - 48.27)||0.1508|
|Dominant model||-||-||1.705 (1.11 - 2.62)||0.0194|
|T allele frequency||0.21||0.14||1.65(1.127 - 2.405)||0.0125*|
|TYMS 5’ UTR 28bp tandem repeats (rs 45445694)||3R3R||57(28.5)||76(40.64)||1.00||Referent|
|3R2R||95(47.5)||92(49.2)||1.38 (0.88 - 2.15)||0.1974|
|2R2R||48(24)||19(10.16)||3.37 (1.79 - 6.34)||0.0002*|
|Recessive model||-||-||2.79 (1.57 - 4.96)||0.0005*|
|Dominant model||-||-||1.18 (1.12 - 2.62)||0.0161|
|2R frequency||0.48||0.35||1.72 (1.28 - 2.29)||0.0003*|
|TYMS 3’ UTR 6bp deletion (rs16430)||6bp0/6bp0||68(34)||73(39.04)||1.00||Referent|
|6bp0/6bpdel||99(50)||76(40.64)||1.40 (0.90 - 2.18)||0.1726|
|6bpdel/6bpdel||33(16.50)||38(20.32)||0.93 (0.5264 - 1.65)||0.9244|
|Recessive model||-||-||0.78 (0.46 - 1.30)||0.4014|
|Dominant model||-||-||1.24 (0.82 - 1.88)||0.3558|
|6bpdel allele frequency||0.41||0.406||1.03 (0.77 - 1.37)||0.9213|
Table 1: Association of MTHFD1 G1958A, MTHFR C677T, TYMS 5’ UTR 28 bp tandem repeats and TYMS 3’ UTR 6bp deletion at 1494 polymorphism and Down syndrome risk. *P-Value significant after Bonferroni multiple correction test.
MTHFR C677T polymorphism: Significant difference in distribution of MTHFR 677T allele was observed between case mothers (0.21) and control mothers (0.14) with odds ratio 1.65(95% CI: 1.13-2.41; p=0.0125). The SNP was found to be associated with increased risk in both models co-dominant (odds ratio 1.61, 95% CI: 1.04-2.48; p=0.0428) as well as in dominant (odds ratio 1.705, 95% CI: 1.11-2.62; p=0.0194). Recessive model of association didn't show significant correlation for risk DS with odds ratio 5.753 (95% CI: 0.69-48.27; p= 0.1508).
TYMS 5'UTR 28 bp repeats (2R/3R) polymorphism: TYMS 5'UTR 28 bp repeats (2R/3R) polymorphism was found to be most significantly associated to increase risk in present study. Minor allele (2R) frequency was significantly higher in case mothers (0.48) compared to control mothers (0.35) with odds ratio 1.72 (95% CI: 1.28-2.29; p=0.0003). TYMS 2R2R genotype was significantly higher in the case of mothers compared to control mothers with odds ratio 3.37 (95% CI: 1.79-6.34; p=0.0002). Recessive model, as well as dominant models of association, revealed significant association to increase risk with odds ratio 2.79 (95% CI: 1.57-4.96; p=0.0005) and 1.18(95% CI: 1.12-2.62; p=0.0161) respectively.
TYMS 3'UTR 6 bp deletion polymorphism: TYMS 3'UTR 6 bp deletion polymorphism didn't reveal association either in allele frequency distribution or in genotypic distribution among case and control mothers. All the model of association remained statistically insignificant among case and control mothers for this polymorphism. Minor allele (6bpdel) frequency distribution was almost same among case mothers (0.410) and control mothers (0.406).
Genotypic combination: Following univariate analysis, we did genetic combination analysis of TYMS 5'UTR 28 bp repeats polymorphism (showing most significant association) with other studied polymorphism. In combination with MTHFD1 G1958A, the only combination TYMS 2R2R/MTHFD1 1958GA showed an additive effect of increasing risk in the case mothers compared to control mothers. Genotypic combinations of TYMS 5'UTR 28 bp repeats polymorphism with MTHFR C677T like 3R3R/CT, 3R2R/CT, 3R2R/TT, 2R2R/CC and 2R2R/CT were significantly higher in the case mothers compared to control mothers (Table 2). In case of genetic interaction of TYMS 5'UTR 28 bp repeats (2R/3R) with TYMS 3'UTR 6bp del polymorphism, the genotypic combinations 3R3R/ 6bp0-6bpdel, 3R2R/6bp0-6bp0, 3R2R/6bp0-6bpdel, 2R2R/6bp0-6bp0, and 2R2R/6bp0-6bpdel were turned out significantly higher in the case mothers compared to control mothers.
|TYMS 5’UTR 28 bp repeats (2R/3R) /MTHFR C677T||3R3R/CC||3R3R/CT||3R3R/TT||3R2R/CC||3R2R/CT||3R2R/TT||2R2R/CC||2R2R/CT||2R2R/TT|
|Case mothers (200)||24||31||1||59||32||5||38||10||0|
|Control mothers (187)||47||26||1||72||22||0||17||2||0|
|Odds ratio 95% CI||Referent||2.24 (1.14 - 4.78)||1.20 (0.12 - 32.72)||1.61 (0.88 - 2.93)||2.85 (1.37 - 5.93)||21.33 (1.13 - 402.02)||4.38 (2.06 - 9.31)||9.79 (1.98 - 48.32)||-|
|TYMS 5’UTR 28 bp repeats (2R/3R) /TYMS 3’UTR 6bp0/6bpdel||3R3R/6bp0-6bp0||3R3R/6bp0-6bpdel||3R3R/6bpdel-6bpdel||3R2R/6bp0-6bp0||3R2R/6bp0-6bpdel||3R2R/6bpdel-6bpdel||2R2R/6bp0-6bp0||2R2R/6bp0-6bpdel||2R2R/6bpdel-6bpdel|
|Case mothers (200)||8||33||16||36||45||14||24||20||4|
|Control mothers (187)||25||27||21||39||40||15||9||8||3|
|Odds ratio 95% CI||Referent||3.82 (1.49 - 9.83)||2.38 (0.85 - 6.66)||4.42 (1.15 - 7.21)||3.52 (1.43 - 8.67)||2.92 (0.99 - 8.58)||8.33 (2.76 - 25.16)||7.81 (2.49 - 24.50)||4.16 (0.76 - 22.72)|
|TYMS 5’UTR 28 bp repeats (2R/3R)/MTHFD G1958A||3R3R/GG||3R3R/GA||3R3R/AA||3R2R/GG||3R2R/GA||3R2R/AA||2R2R/GG||2R2R/GA||2R2R/AA|
|Case mothers (200)||9||21||26||13||50||32||8||33||8|
|Control mothers (187)||11||34||29||17||51||26||2||9||8|
|Odds ratio 95% CI||Referent||0.76 (0.27 - 2.13)||1.10 (0.39 - 3.06)||0.94 (0.30 - 2.92)||1.20 (0.46 - 3.14)||1.50 (0.54 - 4.18)||4.89 (0.82 - 29.07)||4.48 (1.42 - 14.13)||1.22 (0.33 - 4.57)|
Table 2: Genotypic combination of TYMS 5’ UTR 28 bp repeat polymorphism with other polymorphism.
Multivariate analysis: When case and control mothers were stratified and compared in terms of the number of variant alleles harbored, we found more than four variant alleles harboring case mothers were significantly more compared to such control mothers with odds ratio 6.52 (95% CI: 2.17 - 19.61; p=0.008) (Table 3).
|Number of risk allele||Case mothers (%) N=200||Control mothers (%) N=187||OR(95% CI)||p value|
Table 3: Incidence of risk genotypic presence among the case and controls mothers.
Linkage dis-equilibrium (LD) analysis
We performed linkage dis-equilibrium (LD) analysis for the two polymorphisms in TYMS gene. D' value was 0.288 with a LOD of 3.89 and an r2 value of 0.041 indicating weak linkage disequilibrium between the two polymorphisms. TYMS 5'UTR 28 bp repeat 2RTYMS 3'UTR 6bp0 haplotype was significantly higher in case mothers (0.346) compared to control mothers (0.251) with p-value 0.0038. On another side, 3R-6bp0 haplotype was significantly higher in control mothers (0.347) compared to case mothers (0.347) with a p-value of 0.0008 (Table 4).
|Haplotypes TYMS 5’UTR 28 bp repeats (2R/3R) and TYMS 3’UTR 6bp0/6bpdel||Frequency||Case ratio (frequency) N=200||Control ratio (frequency) N=187||Chi square- value||p-value|
|2R-6bp0||0.297||92.2/281.8 (0.246)||137.4/262.6 (0.343)||8.725||0.0031|
|3R-6bp0||0.294||129.8/244.2 (0.347)||97.6/302.4 (0.244)||9.905||0.0016|
|3R-6bpdel||0.288||112.2/261.8 (0.300)||110.4/289.6 (0.276)||0.54||0.4622|
|2R-6bpdel||0.122||39.8/334.2 (0.107)||54.6/345.4 (0.137)||1.625||0.2024|
Table 4: Analysis of haplotypes in MDS and control mothers.
Transmission dis-equilibrium test (TDT) analysis
TDT analysis was performed in a total of 146 DS trios families included in the study to see preferential transmission of alleles from either of the heterozygous parents to DS children (Table 5).
|Genotype||Paternal alleles||Maternal alleles||Total alleles|
|MTHFD1 G1958A G A χ2 – Test P-value||23 24 0.021 0.88||29 30 0.017 0.90||107 133 2.86 0.09|
|TYMS 5’UTR 28 bp tandem repeats 3R 2R χ2 – Test P-value||29 20 1.62 0.20||19 30 2.43 0.12||120 94 3.15 0.07|
|TYMS 3’UTR 6 bp deletion 6bp 6bp deletion χ2 – Test P-value||33 21 2.67 0.10||34 27 0.80 0.37||145 93 11.36 0.0008|
Table 5: Transmission disequilibrium test analysis (TDT).
Although not statistically significant, we observed consistent trend for preferential transmission of alleles to DS children. In case of paternal allele transmission, wild-type alleles TYMS 5'UTR 3R and TYMS 3'UTR 6bp0 (except MTHFD1 1958G) were transmitted more compared to respective risk allele, whereas in case of maternal alleles transmission, risk alleles TYMS 5'UTR 2R, MTHFD1 1958A (except TYMS 3'UTR 6bpdel) were preferentially transmitted to DS children, but this preferential transmission was statistically not significant from both the parents. When we evaluated the total transmitted alleles we found a borderline significant transmission of A allele of MTHFD1 polymorphism and 3R allele of TYMS 28 bp tandem repeat polymorphism. We also found a very significant transmission of wild allele of TYMS 3'UTR 6 bp deletion polymorphism in children (Table 5).
MTHFR C677T has been widely studied polymorphism in case of Down syndrome in different populations and in a recent meta-analysis it has been reported to be independent risk factor along with MTRR A66G and RFC1 G80A . In the present study, MTHFR 677T allele (dominant model MTHFR CT+TT vs CC and co-dominant model MTHFR CT vs TT) was observed to be an associated risk factor and was in concordance with our previous findings [9,11].
In the present study MTHFD1 G1958A polymorphism affecting availability of 5,10-methelene THF didn't turn out to be associated risk with respect to DS. There are conflicting results for association of MTHFD1 G1958A polymorphism with risk of DS [17-20]. The combined genotype MTHFD1 1958 AA/RFC1 80GG has also been reported to be associated with risk for DS . Similarly, in our study the combined genotype of MTHFD1 1958GA/TYMS 5'UTR 2R/2R was associated with increased risk in DS mothers (Table 2). Also, in a meta-analysis, the dominant model of MTHFD1 (GA vs. GG) has been reported to be an associated risk factor . MTHFD1 gene polymorphism has been shown to be associated with increased risk in case of neural tube defects, congenital heart defect and unexplained second-trimester pregnancy loss [14,23,24].
In the present study, TYMS 5'UTR 28bp repeat polymorphism independently turned out associated risk factor in DS mothers. But TYMS 3'UTR region polymorphism didn't reveal any association with DS risk. TYMS and MTHFR genes have been suggested to affect drift of pool of 5,10-methyleneTHF either towards remethylation or nuleic acid biosynthesis pathway . TYMS 5'UTR region 28bp repeats polymorphism is associated with gene expression  and 3'UTR 6bpdel polymorphism affect stability of m-RNA in cytoplasm  thus affecting formation of dTTP for nucleic acid biosynthesis. Apparently, interaction of genetic polymorphism of TYMS and MTHFR might be playing crucial role in chromosome damage and mal-segregation of the chromosome and risk for DS becasuse non-disjunction of chromosome 21 in 95% of the cases is during meiosis I during embryonic development of mothers in grandmothers when the cell have high demand of nucleic acid biosynthesis.
TYMS 3'UTR 6bpdel/6bpdel genotype has been reported to be associated with high level of red folate cell and lower homocysteine when compared with individual with TYMS 6bp0/6bp0 genotype . TYMS 6bpdel/6bpdel genotype has been shown associated with reduced risk of certain cancer [15,29,30]. It is associated with reduced risk for spina bifida in a study from Calfornia . TYMS 6bpdel/6bp0 genotype has been reported to associate to have increased risk over TYMS 6bp0/6bp0 in case of conotruncal heart defect in population of California . In our case, TYMS 3'UTR 6bp0/6bpdel polymorphism didn't turn out significantly altered in distribution of alleles and genotypes.
We also looked into gene-gene interactions; the effect of TYMS 5'UTR 28bp repeats polymorphism in combination with MTHFR C677T, MTHFD1 G1958A and TYMS 3'UTR 6bp0/6bpdel polymorphism was analyzed. TYMS 5'UTR 2R allele was found to be independent risk factor and it is obvious from Table 2 that, MTHFR 677T allele has additive effect. MTHFR 677CT was found to interact with all genotype of TYMS 5'UTR polymorphism (even with 3R3R) to have significant association for increasing risk of DS. The interaction of TYMS polymorphisms was not of significant scale with respect to mothers with MTHFR 677TT genotype due to a few sample sizes. Polymorphism in MTHFD1 G1958A was not found to have significant additive effect to TYMS 5'UTR 28bp repeats polymorphism except for one combination TYMS 2R2R/MTHFD1 GA. We observed additive effect of TYMS 3'UTR polymorphism with TYMS 5'UTR polymorphism. Further, weak linkage dis-equilibrium was observed for the two polymorphisms of TYMS in our case. Haplotypes TYMS 5'UTR 2R-TYMS 3'UTR 6bp0 and TYMS 5'UTR 3R- TYMS 3'UTR 6bp0 were significantly higher in case mothers and control mothers respectively. It is obvious from this interaction that TYMS 5'UTR 3R allele and TYMS 3'UTR 6bp0 allele synergistically are protective in control mothers. Linkage disequilibrium has been reported for TYMS 3'UTR 6bpdel allele and TYMS 5'UTR 3R allele [33,34] but weak linkage diseulilibrium was reported for this haplotype revealing protective effect in another study by Zhu et al., . To the best of our knowledge, there is only one study  so for available for association of these two polymorphisms in case of DS. MTHFR 1298AC/TYMS 2R/2R genotypic combination is associated with reduced risk factor for DS in Italian population . In another study, the 2R/2R genotype has been reported to be protective for coronary artery disease  and in drug response for certain cancer treatment [27,36,37] via interaction with folate intake. Further in a study of Artificial neural network analysis by Coppede et al., revealed TYMS 2R/3R genotype is closely connected with a condition of being mothers of DS and TYMS 3R/3R genotype connected with control mothers . TYMS 2R/2R genotype has been reported to increase risk of spina bifida in infants .
The present study also included 146 DS trios families; we performed transmission disequilibrium analysis to see preferential transmission of any allele from heterozygous parents. We didn't observe significant preferential transmission of any of the alleles in DS children from their heterozygous parents. Although not statistically significant, we observed consistent trend for preferential transmission of alleles to DS children. Wild-type alleles TYMS 5'UTR 3R, TYMS 3'UTR 6bp0 and mutant allele i.e., A allele of MTHFD1 G1958A polymorphism (risk allele) were transmitted more to DS children (Table 5). The MTHFD1 1958A genotype was shown to be decreasing the enzymatic efficiency thereby allowing the substrate to enter more nucleotide biosynthesis. Similarly both the wild alleles of the TYMS polymorphisms studied will promote nucleotide synthesis more promptly than mutant alleles. All these conditions favor nucleotide synthesis which is very crucial during embryonic development. By favoring nucleotide biosynthesis, we believe that preferential transmission of TYMS 3'UTR 6bp0 allele and border line preference of an allele of MTHFD1 G1958A, 3R allele of TYMS 5’UTR in DS children helps in selective advantage to the growing fetus. Further studies are required to validate our observations.
Deficiency of folate and Vit-B12 leading to impaired folate metabolism resulted into aberrant DNA methylation, chromosome breakage, increased frequency of micronuclei formation, defective chromosome recombination and aneuploidy has been already reported [39,40]. Association of genetic polymorphism of folate metabolism with genomic instability and mal-segregation chromosome has been reported in human blood cells [8,41-43]. Further, increased frequency of chromosome damage and mal-segregation of chromosome 13 and 21 have been reported in younger aged mothers compared to control mothers revealing that mothers of DS children have tendency of chromosomal damage and mal-segregation not only in germ lines but also in somatic cells . In another study, Coppede 2009 showed correlation of increased frequency micronuclei formation with MTHFR C677T and MTHFR A1298C in lymphocytes of MDS . In the same study, combined genotypes MTHFR 677TT/MTR 2756AA and MTHFR 1298AC/TYMS2R2R were associated as increased and decreased risk for DS respectively. In the present study, MTHFR C677T and TYMS 5'UTR 28bp repeat polymorphisms independently as well as in combination with each other are significantly associated with risk for DS. Interestingly, both the genes are dependent on the same substrate molecule 5,10-methlenetetrahydrofolate. From our earlier and recent work, we have demonstrated that gene-gene interactions (gene polymorphisms of folate-homocysteine) and micronutrient status (folate, and vitamin B12) independently and cumulatively are significantly associated risk factors in mothers for DS baby birth [9-11].
Coppede have shown that TYMS 5'UTR 2R/2R acts as a protective genotype in control mothers and the logic put forward was that this genotype will allow a slight decrement in nucleotide biosynthesis and thereby promoting more remethylation reaction . This, in turn, will allow proper DNA methylation at pericentromeric region and thereby preventing the nondisjunction. But in a population with folatehomocysteine impairment (like the present cohort), a condition of low micronutrient levels and/or genetic predisposition (MTHFR 677T and other gene polymorphisms), will lead to a limited synthesis of 5, 10- methylenetetrahydrofolate. This will restrict the formation of 5- methylenetetrahydrofolate for remethylation of homocysteine resulting into hyperhomocysteinemia and DNA hypomethylation. If the nucleotide biosynthesis is also hindered in such population, it might prove selective disadvantage for the growing fetus. This can be explained as follows, in 90% of DS cases, non-disjunction of chromosome 21 occurs at meiosis l of gametes formation in embryonic development stage of mothers. At this time, there is high requirement of nucleic acid biosynthesis for DNA replication and repair. Also, in developing fetus there is high need of nucleic acid biosynthesis as well as methylation for proper development of fetus.
In summary, we found significant association of TYMS 5'UTR 2R/2R and MTHFR C677T independently and in combination with increased risk of mothers predisposed to DS baby birth. Two other SNPs, TYMS 3' UTR 6bp del/6bp del and MTHFD1 G1958A were not associated as risk factor for DS, however in genetic combination with TYMS (5'UTR 2R/2R), certain combinations were acting as risk factor. We observed an epistatic interaction between genes involved nucleotide biosynthesis and remethylation and these interactions were significantly associated with genetic predisposition in DS mothers. Replication studies are necessary to confirm our results to elucidate the underlining mechanism.
We thank MDS and control mothers for participation in the study. We record words of gratitude to our teacher Professor Rajiva Raman for critical analysis and comments. We are thankful to Centre for Genetic Disorders, BHU for chromosomal analysis. This work forms part of the contribution of the Disease biology thrust area under DBT sponsored Interdisciplinary School of Life Sciences, BHU. Indian Council of Medical Research (ICMR), New Delhi, India is highly acknowledged for fellowship.
The authors declare no conflict of interest.
Conceived and designed the experiments: S.K.J., K.K.S., A.K.R. Performed the experiments: S.K.J. and S.K.M Analyzed the data: K.K.S., S.K.J., and A.K.R. Wrote the manuscript: K.K.S., S.K.J., and A.K.R. Grants for experiments provided by: A.K.R. Clinical samples provided by: A.K. and A.R.L. Collected the samples: S.K.J. Gave critical comments on the manuscript: A.K. and A.R.L.