Abstract

The aim of this study is to quantify global DNA methylation and investigate the relationship with diabetes status and polymorphisms in MTHFR C677T and NOS3 G894T genes in mixed ancestry subjects from South Africa. Global DNA methylation was measured, and MTHFR rs1801133 and NOS3 rs1799983 polymorphisms were genotyped using high throughput real-time polymerase chain reaction and direct DNA sequencing. Of the 564 participants, 158 (28%) individuals had T2DM of which 97 (17.2%) were screen-detected cases. Another 119 (21.1%) had prediabetes, that is, impaired fasting glucose, impaired glucose tolerance, or the combination of both, and the remainder 287 (50.9%) had normal glucose tolerance. Global DNA methylation was significantly higher in prediabetes and screen-detected diabetes than in normal glucose tolerance (both ) and in screen-detected diabetes compared to known diabetes on treatment (). There was no difference in global DNA methylation between known diabetes on treatment and normal glucose tolerance (). In multivariable linear regression analysis, only NOS3 was associated with increasing global DNA methylation (; 95% CI: 0.286 to 1.560). The association of global DNA methylation with screen-detected diabetes but not treated diabetes suggests that glucose control agents to some extent may be reversing DNA methylation. The association between NOS3 rs1799983 polymorphisms and DNA methylation suggests gene-epigenetic mechanisms through which vascular diabetes complications develop despite adequate metabolic control.

1. Introduction

Type 2 diabetes mellitus (T2DM) results from an interaction of environmental and genetic factors. Despite the initial enthusiasm from the identification of risk loci for T2DM [1, 2], the clinical utility of these genetic markers for T2DM risk prediction and reduction has remained limited. It is becoming increasingly evident that the gene-environment interaction in T2DM and other diseases can be explained in part by epigenetics. Indeed, a number of recent studies using different types of biological tissues from pancreas to peripheral blood mononuclear cells (PBMCs) are showing that DNA methylation patterns are altered in subjects with diabetes [36]. T2DM rates are increasing rapidly in African populations, with traditional diabetes risk factors failing to explain a great deal of these increases. Emerging alternate pathways involve the interplay between epigenetics, microbiome, and the immune system. Epigenetic mechanisms include DNA methylation, lysine methylation, histone methylation, histone phosphorylation, RNA interference (RNAi), and genomic imprinting [7]. DNA methylation is the most extensively investigated epigenetic mechanism; however, DNA methylation profiling of populations from Africa remains to be investigated.

DNA methylation is an epigenetic process characterized by a covalent modification with the addition of a methyl group to the carbon at position 5 of cytosine nucleotides (H5), process catalyzed DNA methyltransferases (DNMTs) [8]. S-Adenosyl-methionine (SAM-CH3) donates the methyl group in a process catalyzed by methylenetetrahydrofolate reductase (MTHFR), an important enzyme in the folate metabolic pathway. Collectively, this folate and methionine cycles form the transmethylation pathway, which is modulated by a number of genetic and environmental factors that in turn affect DNA methylation [9]. In this regard, DNA hypomethylation has been reported in liver biopsies of subjects with type 2 diabetes (T2DM) with low folate levels [10], while MTHFR polymorphisms have been associated with decreased enzyme activity. Two MTHFR polymorphisms, +677C/T (rs1801133) and +1298A/C (rs1801131), have been widely studied and in vitro studies have shown these polymorphisms to decrease the enzyme activity by 35% and 60%, respectively [11, 12], consequently causing DNA hypomethylation. Another gene that has been shown to be affected by folate supplementation is nitric oxide synthase (NOS) [13, 14]. Nitric oxide synthase (eNOS, also referred to as NOS3) is an endothelial enzyme whose function is to synthesize nitric oxide (NO) from L-arginine [15, 16]. Nitric oxide is a unique molecule with diverse physiologic regulatory functions such as smooth muscle relaxation, inhibition of platelet aggregation, immune regulation, neurotransmission, and blood pressure regulation [1719].

Therefore, this study was designed to examine global DNA methylation in individuals with and without diabetes and to evaluate whether polymorphisms in genes involved in DNA methylation and/or folate metabolism, namely, MTHFR and NOS3, have an effect on global DNA methylation in mixed ancestry subjects from South Africa.

2. Materials and Methods

2.1. Ethical Approval of the Study

The study was approved by the Faculty of Health and Wellness Sciences Ethics Committee of the Cape Peninsula University of Technology (CPUT) (NHREC: REC-230408-014) and was conducted according to the code of ethics of the World Medical Association (Declaration of Helsinki, 1975). All participants who were recruited for the study voluntarily signed written consent after the procedures had been fully explained in the language of their choice. Permission to conduct the study was granted by relevant authorities such as city and community authorities as previously described [20].

2.2. Study Design and Population

The present study was cross-sectional by design, involving participants from a mixed ancestry ethnic population group residing in Bellville South Township in Cape Town, South Africa. A detailed description of the survey and procedures conducted in the study are available elsewhere [21, 22]. Eligible participants were those who are older than 20 years, are residing in Bellville South, are of mixed ancestry origin, are not pregnant, and are not acutely ill.

2.3. Clinical Data

Clinical data which has been previously described [20] was collected in the form of a standardized questionnaire. During this time physical examination was conducted with data collection on blood pressure according to World Health Organization (WHO) guidelines [23] using a semiautomatic digital blood pressure monitor (Rossmax PA, USA) on the right arm in sitting position and anthropometric measurements. Body weight was measured to the nearest 0.1 kg with a Sunbeam EB710 digital bathroom scale, which was calibrated and standardized using a weight of known mass. Measurements were recorded with each subject wearing light clothing, without shoes and socks. Waist circumference was determined using a nonelastic tape at the level of the narrowest part of the torso, as seen from the anterior view. All anthropometric measurements were performed three times and their average was used for analysis. Participants with no history of doctor diagnosed diabetes mellitus underwent a 75 g oral glucose tolerance test (OGTT) as recommended by WHO [24].

2.4. Biochemical Analysis

Blood samples were collected from participants after overnight fast and processed as described in [20]. Plasma glucose levels and glycated haemoglobin (HbA1c) in whole blood collected in an EDTA tube were measured, respectively, by enzymatic hexokinase method and the turbidimetric inhibition immunoassay (Cobas 6000, Roche Diagnostics, Germany). Serum insulin was determined by a microparticle enzyme immunoassay (AxSYM, Abbott). High-density lipoprotein cholesterol (HDL-C) and triglycerides (TG) were estimated by enzymatic colorimetric methods in serum (Cobas 6000, Roche Diagnostics). Low-density lipoprotein cholesterol (LDL-C) was calculated using Friedwald’s formula [25]. Body mass index (BMI) was calculated as weight/height height (kg/m2) and waist-hip-ratio (WHR) as waist/hip circumference. Diabetes was based on a history of doctor diagnosis, fasting blood glucose concentration ≥ 7.0 mmol/L (or 126 mg/dL), and/or 2-hour post-OGTT plasma glucose ≥ 11.1 mmol/L (or 200 mg/dL).

2.5. Genotyping and DNA Methylation

Genomic DNA was extracted from whole blood samples collected in EDTA tubes using the salt extraction method and quantified with the NanoDrop ND-1000 instrument (Nanodrop Technologies, Wilmington, USA). Global DNA methylation was measured with the 5mC DNA ELISA kit according to the manufacturer’s instructions (Zymo Research Corp., Irvine, CA, USA). Briefly, 100 ng of genomic DNA was allowed to bind to the enzyme linked immunosorbent assay (ELISA) plate, and the methylated fraction of DNA was detected using a 5-methylcytosine monoclonal antibody and quantified by an ELISA-like reaction and measuring the absorbance at 450 nm. The MTHFR rs1801133 and NOS3 rs1799983 polymorphisms were genotyped using high throughput real-time polymerase chain reaction (RT-PCR) on the BioRad Optica (BioRad, USA) platform using Taqman genotyping assay (Applied Biosystems, USA). Conventional polymerase chain reaction followed by direct DNA sequencing was performed for analytical validation of high throughput genotyping.

2.6. Statistical Analysis

General characteristics of the study participants are summarized as count and percentage for dichotomous traits and median and 25th–75th percentiles for quantitative traits. Group comparisons used chi square and Kruskal-Wallis tests. Monotonous trends in the distribution of characteristics across quarters of global DNA methylation were assessed using the Cochran-Armitage trend test for proportions and Jonckheere-Terpstra trend test for medians. Robust regressions were then used to assess the effect of various traits on global DNA methylation in models accounting for age, gender, smoking, and glucose tolerance status. SNPs were tested for departure from Hardy-Weinberg Equilibrium (HWE) expectation via a chi square goodness of fit test, and their effect on DNA methylation was investigated in linear regressions models, while always assuming a log-additive genetic model. Results corresponding to values below 5% are described as significant. All analyses used the statistical software R (version 3.2.2 (2015-08-14), R Foundation for statistical computing, Vienna, Austria). SNPs analyses used the packages “genetic” and “SNPassoc.”

3. Results

3.1. Global DNA Methylation and Cardiometabolic Profile

The study cohort consisted of 564 participants, including 438 (78%) females. One hundred and fifty-eight (28%) individuals had T2DM of which 97 (17.2%) were screen-detected cases. Another one hundred and nineteen (21.1%) had prediabetes, that is, impaired fasting glucose, impaired glucose tolerance, or the combination of both, and the remainder 287 (50.9%) had normal glucose tolerance. Global DNA methylation (%) was significantly higher in subjects with prediabetes or diabetes when compared to individuals with normoglycemia (). However, no significant differences were observed between subjects with prediabetes and diabetes. On further analysis, higher global DNA methylation in diabetics was driven by those with screen-detected diabetes. As shown in Figure 1, global DNA methylation was significantly higher in subjects with screen-detected diabetes than in those with normoglycemia () and remained significantly higher when compared to individuals with known diabetes (). However, there was no significant difference between known diabetic and normoglycemia subjects (). The characteristics of participants by glucose tolerance status are summarized in Table 1, mostly showing the expected differences.

The baseline characteristics of participants across quarters of global DNA methylation are summarized in Table 2. The proportion of participants with any diabetes increased across increasing quarters of global DNA methylation () and in a linear fashion ( for linear trend). Across quarters of global DNA methylation levels, significant differences were also observed in the distribution of fasting () and 2-hour glucose levels (), with increasing trends across quarters of global DNA methylation, always in linear fashions (both ). Significant correlations were also observed between global DNA methylation and diabetes (; ), body mass index and waist circumference (both, ; ), fasting blood glucose (; ), and after 2-hour blood glucose (; ). In robust linear regression analysis adjusted for age, gender, status for hyperglycemia, and smoking, any diabetes (; ) was associated with global DNA methylation. When participants with diabetes were distinguished into those with screen-detected diabetes or known diabetes, the association remained significant only for screen-detected diabetes (; ) but not for known diabetes. In similar regression analyses, no other characteristics emerged as a significant predictor of DNA methylation, although borderline associations were observed with female gender, gamma-GT, and systolic blood pressure (Table 2). The association of any diabetes with DNA methylation also remained significant after expansion of the basic models to include the following correlates of diabetes status in the sample: BMI, waist circumference, systolic blood pressure, drinking status, HDL-cholesterol, total cholesterol, triglycerides, log-transformed CRP, and log-transformed gamma-GT (; ).

3.2. Genetic Polymorphisms and DNA Methylation

Table 3 shows the genotype distribution of polymorphisms across the quarters of global DNA methylation. None of the single nucleotide polymorphisms (SNPs) deviated from Hardy-Weinberg Equilibrium (HWE) (all ), and their distribution showed neither significant difference (all ) nor linear trends (all for linear trend) across the quarters of global DNA methylation. In linear regression analysis adjusted for age, gender, status for hyperglycemia, and smoking, only T allele of NOS3 was associated with increasing global DNA methylation with a beta coefficient of 0.943 (95% confidence interval: 0.286 to 1.560) assuming a log-additive genetic model. This association remained significant in sensitivity analysis after exclusion of participants with known diabetes. In this subgroup, the beta coefficient was 1.102 (95% CI: 0.401 to 1.802) for the effect of T allele of NOS3 on global DNA methylation.

4. Discussion

In this study, we measured global DNA methylation in PBMCs of mixed ancestry individuals with different levels of glucose tolerance, in parallel with genetic screening of polymorphisms in the MTHFR and NOS3 genes. Our results show that global DNA methylation is increased in both prediabetic and diabetic states, but this increase was more pronounced in those with screen-detected diabetes, even after adjustment for extraneous factors such as age, gender, smoking, and glucose tolerance status. The distribution of the investigated SNPs did not differ across quarters of DNA methylation, but regression analysis under log-additive genetic model assumption revealed NOS3 to be an independent determinant of global DNA methylation.

The strengths of our study include a well-characterized cohort and using the reference standard OGTT and a detailed drug history to carefully assign individuals to glucose tolerance status. Both T2DM and DNA methylation are affected by age, obesity, and genetic factors. In this regard, we also evaluated the effect of genes involved in DNA methylation and/or folate metabolism. We quantified DNA methylation across the whole genome using enzyme linked immunosorbent assay (ELISA), which is cost-effective and easy technique that is amenable for routine analysis in developing countries. The use of robust analytic methods also increases the reliability of our findings. For instance, by implementing robust correlations to eliminate the effect of outliers, we provide better estimates of the measures of association without loss of power [26].

In this study, we show that DNA methylation is significantly increased in subjects with screen-detected diabetes compared to those with diabetes and those on treatment. We believe that the lack of association between known diabetes and global DNA methylation in our sample is likely due to the reversal of DNA methylation changes by glucose controls agents in these individuals. These findings are in contrast with a study that quantified global DNA methylation by bisulfite pyrosequencing and showed a null effect of diabetes medication on the association between DNA methylation and insulin resistance [27]. The difference between Zhao et al. [27] and our study may be attributed to the sample size of these two studies. Their study analyzed global DNA methylation of 17 diabetics on glucose control agents, whereas 61 diabetics on medication were included in our study. Genetic polymorphisms have been shown to be responsible for individual responses to many oral antidiabetic drugs [28, 29], and these polymorphisms may also be subject to epigenetic regulation. It is thought that antidiabetic drugs could influence methylation trends through gene expression mechanisms and may represent a possible confounding factor for the identification of T2DM related epigenetic profiles [3]. Furthermore, our results are in contrast with those reported in liver biopsies of subjects with T2DM where genome-wide analysis of DNA methylation identified 94% hypomethylated CpG sites in these individuals [10]. Although diabetic subjects in this study were similarly obese as in Nilsson et al. [10], the differences in the cell types and liver versus blood leucocytes as well as the different methods for quantifying DNA methylation may explain differences in DNA methylation.

Although global DNA methylation correlated positively with adiposity indices, there was no correlation with age. The lack of association between DNA methylation and age is at variance with a number of studies that have repeatedly shown decreasing DNA methylation with increasing age [3032]. The narrow age range in our study has likely masked any association between DNA methylation and age. We also did not find an association between the MTHFR C677T polymorphism and global DNA methylation, despite previous reports that have shown an association between MTHFR C677T, methyl group generation, and subsequent DNA methylation [33, 34]. Contrary to our prior hypothesis that MTHFR gene polymorphisms will affect global DNA methylation, we observed that global DNA methylation was associated with NOS3 G894T polymorphism. Reduced production of NO is one of the most important contributors to endothelial dysfunction, an initiator of vascular complications [1519]. NOS3 is encoded by the NOS3 gene located on chromosome 7q35-36 and three polymorphisms (G894T, 4b/a, and T786C) in the NOS3 gene are associated with diabetes and diabetes related traits and reduced enzyme activity and consequently bioavailability of NO and endothelial function [19, 35, 36]. It is not clear how NOS3 may influence global DNA methylation, but a study involving obese children and type 1 diabetes with NOS3 polymorphisms demonstrated an improved endothelial function after supplementation with folate [37]. Similarly, the authors obtained differing results with respect to MTHFR as no improvement in endothelial function was evident in those with MTHFR polymorphisms suggesting two independent mechanisms for vascular response to folate for both NOS3 and MTHFR [37]. Folate plays an important role in DNA methylation and several studies have shown aberrant DNA methylation profiles in relation to folate levels [10, 38]. Individuals with T2DM are at an increased risk of both micro- and macrovascular complications and it has previously been shown that glycemic control is not adequate to halt the progression of vascular complications [3941]. Thus, the observed association between NOS3 G894T gene polymorphisms and global DNA methylation suggests gene-epigenetic mechanisms through which vascular complications develop despite adequate glycemic control in individuals with T2DM.

The current study represents a significant contribution in the African setting, where research about the contribution of epigenetics to the growing burden of noncommunicable disease in the region is lacking [42]. In summary, our data indicates that global DNA methylation is significantly associated with screen-detected diabetes compared to known diabetes on treatment, leading us to speculate that glucose control agents to some extent may be reversing DNA methylation. Furthermore, we showed that DNA methylation is independent of MTHFR polymorphism but is influenced by NOS3 gene polymorphism G894T. It is known that epigenetic variations are regulated in a tissue-specific manner; however, the use of peripheral blood to assess DNA methylation is recommended for epidemiologic samples [43]. Furthermore, DNA methylation in blood presents a feasible marker for the diagnosis, prevention, and management of T2DM. To fully elucidate the role of DNA methylation in this population, future studies should explore gene-specific methylation together with global DNA methylation analysis.

Disclosure

Any opinion, finding, and conclusion or recommendation expressed in this material is that of the authors and the MRC does not accept any liability in this regard.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Authors’ Contributions

Tandi E. Matsha made substantial contributions to conception and design, acquisition of data, analysis and interpretation of data, drafting the article, and final approval of the version to be published. Carmen Pheiffer made substantial contributions to conception and design, acquisition of data, revising it critically for important intellectual content, and final approval of the version to be published. Tinashe Mutize made substantial contributions to design, acquisition of data, revising it critically for important intellectual content, and final approval of the version to be published. Rajiv T. Erasmus revised the paper critically for important intellectual content and contributed to final approval of the version to be published. Andre P. Kengne contributed to analysis and interpretation, revising the paper critically for important intellectual content, and final approval of the version to be published.

Acknowledgments

The authors thank the Bellville South (Ward 009) Community for participating in the study. They are also grateful to the Bellville South Community Health Forum for supporting the engagement with the Bellville South Community. This research project was funded by the South African Medical Research Council (MRC) with funds from National Treasury under its Economic Competitiveness and Support Package (MRC-RFA-UFSP-01-2013/VMH Study). Carmen Pheiffer is supported by the South African National Research Foundation Unique (Grant no. 87838).