The MTHFR promoter hypermethylation pattern associated with the A1298C polymorphism influences lipid parameters and glycemic control in diabetic patients

Background Polymorphisms in the gene encoding methylenetetrahydrofolate reductase (MTHFR) have been investigated as risk factors for microvascular complications of diabetes; however, simultaneous analysis of these polymorphisms and the methylation pattern of the gene has never been conducted. The objective of the present study was to evaluate the simultaneous relationship between MTHFR methylation and MTHFR C6TT7 and A1298C polymorphisms with metabolic, inflammatory and oxidative stress parameters related to microvascular complications, diabetic retinopathy (DR) and diabetic nephropathy (DN) in diabetic patients. Methods A total of 107 patients who were diagnosed in the previous 5 to 10 years were recruited and divided into groups with complications (DR and/or DN) or without complications. Methylation analysis of the gene promoter was conducted using the MSP technique, and analysis of the A1298C and C677T polymorphisms was conducted using the restriction fragment length polymorphism (RFLP) assay. Microalbuminuria was determined using urine samples, and other analytes of interest were determined in blood samples using commercial kits. The Mann–Whitney and Chi square statistical tests were used with significance considered at p < 0.05. Results Subjects with a hypermethylated profile and the 1298AA genotype showed the highest levels of blood glucose (p = 0.03), total cholesterol (p = 0.0001) and LDL cholesterol (p = 0.0006). The same profile was associated with higher levels of HbA1c (p = 0.025), glycemia (p = 0.04) and total cholesterol (0.004) in the control group and total cholesterol (p = 0.005) and LDL cholesterol (p = 0.002) in the complications group. Serum creatinine was higher in subjects in the hypermethylated group with the genotype 677CC only in the control group (p = 0.0020). The methylated profile in presence of 677CC + 1298AA and the 677CT/TT +1298AA haplotypes showed higher levels of total cholesterol (p = 0.0024; 0.0031) and LDL cholesterol (p = 0.0060; 0.0125) than 1298AC/CC carriers. The fasting glycemia was higher in hypermethylated profile in the presence of 677CC/1298AA haplotype (p = 0.0077). Conclusion The hypermethylated methylation profile associated with the 1298AA genotype appeared to be connected to higher values of glycemia, total cholesterol and LDL cholesterol.


Introduction
Diabetes mellitus (DM) is a multifactorial disease characterized by chronic hyperglycemia and represents an important public health demand because it causes significant impacts on the quality of life and health management costs of the patients, especially when chronic complications occur. Microvascular complications are common and are related to several factors, especially diabetes duration and poor glycemic control.
Methylenetetrahydrofolate reductase (MTHFR) is an enzyme involved in the re-methylation of homocysteine (Hcy) in methionine cycle and catalyzes the reduction reaction of 5,10-methylenetetrahydrofolate to 5-methylenetetrahydrofolate. Thus, MTHFR function is one of the factors involved in the control of plasma levels of homocysteine. Elevated plasma levels of Hcy have been related to insulin resistance and to chronic complications such as diabetic retinopathy (DR) and diabetic nephropathy (DN) [1][2][3]. MTHFR is involved in the synthesis of the primary universal donor of methyl groups, S-adenosylmethionine (SAM), which is formed from methionine [4]. SAM is a substrate for the DNA methylation process that consists of an epigenetic mechanism that regulates the transcription of several genes and generally inhibits gene expression.
A previous study showed the association between the methylation profile in the MTHFR gene with the occurrence of chronic complications in diabetic patients. Briefly, comparison of the methylation profile of the MTHFR gene was performed in individuals with type 2 diabetes (DM2) who have been diagnosed with DM for 5-10 years. While the unmethylated profile was not identified in the analyzed samples, the partially methylated pattern predominated (~ 60%). The hypermethylated profile was associated with diabetic retinopathy (DR), but not with diabetic nephropathy (DN). Higher levels of total cholesterol and LDL cholesterol were identified in the hypermethylated profile only in DR patients samples [5].
Similarly, other reports, although controversial, describe gene polymorphisms that seem to interfere with the enzymatic activity of MTHFR and predispose individuals to diabetic complications. Studies indicate that the C677T polymorphism in the MTHFR gene, which is associated or not with inadequate glycemic control, may be linked to the progression of diabetic retinopathy [6,7]. In addition, studies point to a significant association between this genotype and the development of nephropathy in diabetic patients [8,9]. The A1298C polymorphism (rs1801131) has also been linked to the development of DN or macroalbuminuria, either alone or in combination with other alleles of interest, including C677T [10,11]. However, the parallel effect of methylation profile and MTHFR polymorphisms on metabolic aspects in diabetic patients has never been analyzed.
Based on these facts, the objective of this work is to investigate in an unprecedented way the simultaneous relationship between MTHFR methylation profile and the C677T and A1298C polymorphisms and metabolic, inflammatory and oxidative stress factors related to microvascular complications in patients with diabetes for 5 to 10 years. photographs were analyzed according to the standards and recommendations of the Action to Control Cardiovascular Risk in Diabetes (ACCORD) and Early Treatment Diabetic Retinopathy Study (ETDRS). The presence of renal impairment was assessed by quantifying albumin in a 24-h urine sample. Microalbuminuria was defined as a urinary albumin concentration > 30 mg/L/24 h.
Clinical variables (age, date of diabetes diagnostic, glycemic control treatment, presence of hypertension, presence of dyslipidemia) were collected using an appropriate form.

Collection of biological samples
For determination of albuminuria, patients received an appropriate container and were instructed to collect all urine for 24 h. They should then take the container to the researchers, ending the 24 h with the last collection in the hospital, at the time of the start of the exam the day before.
The blood was collected by venipuncture and samples were placed in different containers. For biochemical analysis in general, in the presence of clot activator, for determination of HbA1C in the presence of anticoagulant and for determination of glucose in the presence of anticoagulant sodium fluoride. The samples were centrifuged at a speed of 3000 rpm for 10 min at room temperature for separation of serum or plasma and subjected to analysis within 2 h after collection, except the sample for HbA 1 C that was analyzed in hemolyzed whole blood.
For DNA extraction, blood collection was performed by venous puncture in sterile tubes containing 7.2 mg of K 3 EDTA. Blood samples were stored for up to 20 days at − 20 °C until DNA extraction was performed.
Light density lipoprotein cholesterol (LDL) concentration was determined using the Friedewald formula, For determination of malondialdehyde (MDA), the method based on the reaction of thiobarbituric acid with hydroperoxides decomposition products was used. The total antioxidant capacity (TAC) of the serum was evaluated by the reduction of the 2,2-diphenyl-1-picrylhydrazyl radical, which was measured by a colorimetric method in a spectrophotometer.
Urine was used to detect urinary albumin (M.ALB) using a commercial kit according to the guidelines recommended by the manufacturer (Labtest ® , Lagoa Santa, MG, Brazil) in an automated analyzer. For the calculation of albumin excretion in 24 h, the value of the determination of albumin in the urine collected for 24 h (mg/dL) was multiplied by the total volume of urine collected in liters.

DNA extraction for molecular analyzes
Leukocyte DNA was isolated using an protocol adapted from Miller et al. [5,12].

Determination of MTHFR genotypes C677T and A1298C
The presence of polymorphisms was determined using the restriction fragment length polymorphism (RFLP) technique by amplifying a fragment of the MTHFR gene by PCR using primers and temperature cycles described elsewhere [13,14]. The 198 and 256 base pair fragments generated respectively by the amplification reaction of the polymorphisms C677T and A1298C were digested for 3 h with the endonucleases HinfI (C677) and MobII (A1298C) (Biolabs ® ) according to manufacturer's guidelines. The RFLP profiles were generated by separation of DNA fragments in 15% polyacrylamide gel electrophoresis and observed by staining with silver nitrate or using the Gelred (Biotium ® ) intercalating dye followed by excitation by ultraviolet light in transilluminator. The genotypes 677 CC/CT/TT and 1298 AA/AC/CC were identified by the pattern of bands in the polyacrylamide gel according to the expected sizes described in the literature [13,14].

Analysis of the methylation profile in the MTHFR gene
Methylation analysis in the promoter of the MTHFR gene was performed using the methylation specific PCR (MSP) technique as previously described [5]. In short, the extracted leukocyte DNA was converted (500 ng) by sodium bisulfite protocol using EZ DNA Methylation ™ Kit (ZymoResearch) according to the manufacturer's instructions. Methylation-specific PCR reactions were performed using 100 ng of bisulfite-transformed DNA, and specific pairs of primers for methylated and unmethylated targets in two separate amplification reactions. Methylated and unmethylated DNA (Cells-to-CpG ™ Methylated & Unmethylated DNA Control Kit, Life Technologies), which were modified as previously mentioned, and amplified by PCR, were analyzed together with the samples as a control for methylated and non-methylated conditions. Amplified PCR samples were loaded on 3% agarose gels with Gelred (Biotium ® ) and subjected to electrophoresis. The DNA bands were visualized with ultraviolet light. The samples were categorized as methylated (positive only for the methylated condition) or partially methylated (positive for methylated and nonmethylated conditions). Non-methylated samples (positive only for the non-methylated condition) were not found in the population analyzed.

Statistical analysis
The data were compiled in an Excel ® spreadsheet.

Results
A total of 107 diabetic patients, all with diabetes diagnosed within the previous 5 to 10 years, were included in the following two experimental groups: 47 patients with complications (DR and/or DN) and 60 patients without complications (control). Comparison of the clinical and demographic data of the control groups and patients with complications (Table 1) showed that there are differences in diagnosis time that is slightly higher in the group of patients with complications (p = 0.026). The main difference between the groups is in the pharmacological treatment of glycemic control. While in the control group the use of oral hypoglycemic agents predominates (p = 0.006), in the group of patients with complications, the use of insulin associated with oral hypoglycemic agents is significantly higher (p = 0.0002). Initially, there was an association between the genotypes and the clinical condition after evaluating the simultaneous effect of methylation profile and polymorphism in the MTHFR gene in all studied patients, independent of the clinical condition. For this analysis, the patients were divided into two categories according to the methylation profile of the MTHFR promoter: methylated and partially methylated. From the 107 samples analyzed, none presented the unmethylated profile of the MTHFR gene promoter, so the unmethylated category was not included in this study. The relationship of polymorphisms together with the methylation profile and the clinical conditions of the patients were analyzed, taking into account the following four groups: methylated with complications, methylated without complications, partially methylated with complications and partially methylated without complications.

Influence of the methylation profile of MTHFR + polymorphism A1298C
Tables 2 and 3 describe relationships for the A1298C polymorphism. In Table 2, when all individuals were analyzed together, an association between the AA genotype and a hypermethylated promoter with higher glycemic values was identified (p = 0.03). A relationship between this same molecular profile and lipid profile was also observed, showing higher levels of TC (p = 0.0001) and LDL (p = 0.0006). A post hoc test with TC (methylated) was performed. At alpha significance level of 0.05, effect size and power (1-β err prob) were, respectively, 1.44 and 0.998. When the control and complication groups were analyzed separately (Table 3), it was observed that the effect of the genotype + hypermethylated profile appeared for glycemic control parameters fasting glycemia (p = 0.04) and HbA 1 C (p = 0.025) only in the control group. For the lipid profile parameters, while the relationship appeared in both groups for TC (p = 0.0049 and p = 0.0056, respectively), it remained only in the group with complications for LDL cholesterol (p = 0.002).

Influence of the methylation profile of MTHFR + C677T polymorphism
For the C677T polymorphism, the genotypic distribution was fairly balanced between the with complications and without complications phenotypic groups (data not show). For Table 4, analyzing all individuals (complicated + uncomplicated) there were no significant results. In Table 5, analyzing the methylation profile and clinical genotype, only the intersection between CC and TC + TT in the hypermethylated group for creatinine was significant (p = 0.0040). This trend continued when the division of the methylation profile was funneled between the with complications and control groups (p = 0.0020).
In the group with complications, the triglyceride values were higher in the partially methylated group profile in the presence of CT + TT (0.04) genotypes, while the TAC was higher in the same profile but in the presence of the CC genotype (p = 0.01). In addition, the lipid profile was also affected, with the partially methylated group presenting higher levels of HDL (p = 0.0273) in the control group in the presence of CT + TT genotypes.
Influence of the methylation profile of MTHFR + haplotype C677T/A1298C Table 6 shows the results of intersections between the haplotypes C677T and A1298C, divided by the methylation profile (hypermethylated and partially methylated, respectively). For the hypermethylated group, glycemia was significant between 677CC/1298AA and 677CC/1298AC + CC (p = 0.0077), with values from the latter group (mean = 139.9 and median = 129.0) being significantly improved relative to the former (mean = 205.1 and median = 214.0). For total cholesterol (p = 0.0024) and LDL (p = 0.0060), the same trend was observed, with lower values for those with 677CC/1298AC + CC. In the relationship between the       partially methylated individuals and the haplotype no significant difference was found.

Discussion
The main result of this work is the simultaneous effect of the hypermethylation and A1298C polymorphism of the MTHFR gene on lipid parameters and glycemic control in diabetic patients.
The MTHFR gene is located on the short arm of chromosome 1 (1p26.3), has 11 exons and is expressed through the production of three transcripts of different sizes in the human species, suggesting the involvement of complex gene regulation mechanisms [15]. Analysis of the promoter region by computational analysis revealed the presence of two CpG islands with potential binding sites for various transcription factors,    Min.-Max.  suggesting that transcriptional control of the production of this enzyme may undergo methylation control [16]. Few studies have analyzed the effects of the methylation profile of MTHFR and the occurrence of complications in diabetes. A recent study in which samples of healthy, diabetic, early-stage and advanced-stage DN patients were analyzed an association between MTHFR demethylation and DN morbidity was identified in diabetic patients. While the healthy control group did not have demethylated specimens, the diabetic group and the group with early DN showed a higher level of demethylated samples, and the clinical DN group did not present methylated samples [17]. In addition, these results showed an association between increasing levels of homocysteine (Hcy) and demethylation of MTHFR. This effect of the concentration of HCY on the methylation profile of MTHFR has already been demonstrated in in vitro experiments with human muscle cells, which revealed an even higher production of transcripts in the demethylated MTHFR pattern [18]. The authors argue that there may be a compensatory effect causing MTHFR demethylation when HCY is elevated to increase the expression of the enzyme and consequently the conversion of Hcy to methionine. This effect would be evident in the samples of patients with clinical DN, who had high levels of HCY and demethylated profiles [17]. Another study that analysed the relationship between MTHFR methylation and end-stage renal disease, not specifically in DN, has been reported an association between methylation and higher levels of total cholesterol and LDL cholesterol [19]. Diabetic retinopathy and higher levels of total cholesterol and LDL cholesterol has also been associated with the hypermethylated MTHFR pattern in patients with a diagnosis time of diabetes of 5 to 10 years [5], although the relationship has not been identified in samples from patients with DN.

Table 4 Analysis of the relationship between the MTHFR methylation pattern and the C677T polymorphism of MTHFR and biochemical, inflammatory and oxidative stress variables
In this work the objective was to simultaneously analyze the effect of the methylation profile and polymorphisms in the MTHFR gene on biochemical variables taking into consideration the clinical condition of the diabetic patient.
Two polymorphisms in MTHFR, C677T and A1298C, have been most studied because they affect enzymatic activity. The C677T polymorphism, which causes an alanine substitution by valine (Ala222Val), occurs in exon 4 of the MTHFR gene and results in enzyme thermolability and decreased activity. It has been shown that 677TT homozygotes and 677CT heterozygotes have approximately 70% and 35% reduced MTHFR activity, respectively [13,20]. The 677TT genotype is associated with global hypomethylation in the healthy cell genome of cancer patients [21].
The A1298C polymorphism results in the substitution of glutamic acid for alanine (Glu429Ala), and although the homozygote showed a slight reduction in enzymatic activity, it had no effect on Hcy concentration [14]. This polymorphism is located in the regulatory domain of the enzyme responsible for binding to NADPH and S-adenosylmethionine [22].
In this study, while the C677T polymorphism did not appear to affect the distribution of the patients in the with complications and without complications groups, the 1298C allele genotypes (AC + CC) appeared in a higher percentage in the group without complications in diabetes, although without statistical significance, suggesting a protective effect of the C allele. This is the allele that appears to have a minor impact on MTHFR activity, although studies have shown that individuals with the 677CC/1298CC haplotype have lower homocysteine levels than do those with the 677CC/1298AA haplotype [23]. However, the allele 1298C has appeared as a protective factor for other medical conditions, such as congenital heart defects [24]. A possible explanation for this fact, noted by the authors of the cited articles, is the present linkage disequilibrium between the 677T allele, which causes significant reduction in the enzymatic activity and the allele in question, 1298C, which would result in a false protection effect of the mutant allele 1298C. The study further proposes that the 1298C allele may appear in linkage disequilibrium with another yet unrecognized polymorphism of MTHFR. In addition, the 1298C allele was cited as a protective factor for folate deficiency and related to higher levels of this nutrient in comparison with the 1298AA genotype [20,25]. A possible explanation for these findings is related to the location of the polymorphism, which affects the binding site for the negative allosteric effector S-adenosylmethionine, which could therefore interfere with the inhibitory effect of the enzymatic activity [22]. However, there is controversy about these data, and other studies have identified lower levels of folate in carriers of the 1298C allele [26].
In a previous study by our group, a relationship was found between higher levels of total cholesterol and LDL cholesterol in the MTHFR hypermethylated samples of patients with DR [5]. The proposed mechanism for these findings is related to the reduction of MTHFR expression in the hypermethylated profile [18] that could result in higher concentrations of homocysteine, which in turn stimulates hepatic cholesterol synthesis [27].
In the present study, it was demonstrated that individuals with the hypermethylated profile and AA genotype of the A1298C polymorphism had higher levels of TC and LDL and also glycemia and HbA1C than did those with the genotype containing the C allele. Regarding glycemic control, evidenced through fasting