Significant impact of the MTHFR polymorphisms and haplotypes on male infertility risk

Abstract

Background: Methylenetetrahydrofolate reductase (MTHFR) converts 5,10-methylene tetrahydrofolate to 5-methyl tetrahydrofolate and affects the activity of cellular cycles participating in nucleotide synthesis, DNA repair, genome stability, maintenance of methyl pool, and gene regulation. Genetically compromised MTHFR activity has been suggested to affect male fertility. The objective of the present study was to find the impact on infertility risk of c.203G>A, c.1298A>C, and c.1793G>A polymorphisms in the MTHFR gene.

Methods: PCR-RFLP and DNA sequencing were used to genotype the common SNPs in the MTHFR gene in 630 infertile and 250 fertile males. Chi-square test was applied for statistical comparison of genotype data. Linkage disequilibrium between the SNPs and the frequency of common haplotypes were assessed using Haploview software. Biochemical levels of total homocysteine (tHcy) and folic acid were measured. Meta-analysis on c.1298A>C polymorphism was performed using data from ten studies, comprising 2734 cases and 2737 controls.

Results: c.203G>A and c.1298A>C were found to be unrelated to infertility risk. c.1793G>A was protective against infertility (P = 0.0008). c.677C>T and c.1793G>A were in significant LD (D' = 0.9). Folic acid and tHcy level did not correlate with male infertility. Pooled estimate on c.1298A>C data from all published studies including our data showed no association of this polymorphism with male infertility (Odds ratio = 1.035, P = 0.56), azoospermia (Odds ratio = 0.97, P = 0.74), or oligoasthenoteratozoospermia (Odds ratio = 0.92, p = 0.29). Eight haplotypes with more than 1% frequency were detected, of which CCGA was protective against infertility (p = 0.02), but the significance of the latter was not seen after applying Bonferroni correction.

Conclusion: Among MTHFR polymorphisms, c.203G>A and c.1298A>C do not affect infertility risk and c.1793G>A is protective against infertility. Haplotype analysis suggested that risk factors on the MTHFR locus do not extend too long on the DNA string.

Figure 1. Linkage disequilibrium plot.

The number in each cell represents the LD parameter D′ (×100). Each cell is color graduated related to the strength of LD between the two markers. The rs numbers are SNP IDs taken from the Ensembl database. rs1801133, CM056008 NOIATAUM_DMCH, CM056008 NOIATAUM_DMCH, and rs3818762 are in strong LD.

Figure 2. Flow diagram showing inclusion and exclusion of the studies for meta-analysis.

Studies retrieved upon literature search were subjected to inclusion/exclusion criteria as detailed in the methods section.

Figure 3. Forest plots.

Meta-analysis on c.1298A>C SNP in infertility (A), oligoasthenoteratozoospermia (B), and azoospermia (C).

Figure 4. Funnel plot.

Plot of precision by log odds ratio using fixed effect model for observed and imputed sets of studies.