Genome-Wide Analysis of DNA Methylation Differences in Muscle and Fat from Monozygotic Twins Discordant for Type 2 Diabetes

Abstract

Background

Monozygotic twins discordant for type 2 diabetes constitute an ideal model to study environmental contributions to type 2 diabetic traits. We aimed to examine whether global DNA methylation differences exist in major glucose metabolic tissues from these twins.

Methodology/Principal Findings

Skeletal muscle (n = 11 pairs) and subcutaneous adipose tissue (n = 5 pairs) biopsies were collected from 53–80 year-old monozygotic twin pairs discordant for type 2 diabetes. DNA methylation was measured by microarrays at 26,850 cytosine-guanine dinucleotide (CpG) sites in the promoters of 14,279 genes. Bisulfite sequencing was applied to validate array data and to quantify methylation of intergenic repetitive DNA sequences. The overall intra-pair variation in DNA methylation was large in repetitive (LINE1, D4Z4 and NBL2) regions compared to gene promoters (standard deviation of intra-pair differences: 10% points vs. 4% points, P<0.001). Increased variation of LINE1 sequence methylation was associated with more phenotypic dissimilarity measured as body mass index (r = 0.77, P = 0.007) and 2-hour plasma glucose (r = 0.66, P = 0.03) whereas the variation in promoter methylation did not associate with phenotypic differences. Validated methylation changes were identified in the promoters of known type 2 diabetes-related genes, including PPARGC1A in muscle (13.9±6.2% vs. 9.0±4.5%, P = 0.03) and HNF4A in adipose tissue (75.2±3.8% vs. 70.5±3.7%, P<0.001) which had increased methylation in type 2 diabetic individuals. A hypothesis-free genome-wide exploration of differential methylation without correction for multiple testing identified 789 and 1,458 CpG sites in skeletal muscle and adipose tissue, respectively. These methylation changes only reached some percentage points, and few sites passed correction for multiple testing.

Conclusions/Significance

Our study suggests that likely acquired DNA methylation changes in skeletal muscle or adipose tissue gene promoters are quantitatively small between type 2 diabetic and non-diabetic twins. The importance of methylation changes in candidate genes such as PPARGC1A and HNF4A should be examined further by replication in larger samples.

About the Authors:

Rasmus Ribel-Madsen

* E-mail: [email protected]

Affiliations Steno Diabetes Center, Gentofte, Denmark, Section of Metabolic Genetics, The Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark

Mario F. Fraga

Affiliation: Department of Immunology and Oncology, National Centre for Biotechnology, CNB-CSIC, Madrid, Spain

Stine Jacobsen

Affiliation: Steno Diabetes Center, Gentofte, Denmark

Jette Bork-Jensen

Affiliation: Steno Diabetes Center, Gentofte, Denmark

Ester Lara

Affiliation: Department of Cancer Epigenetics, Spanish National Cancer Research Centre, Madrid, Spain

Vincenzo Calvanese

Affiliation: Department of Cancer Epigenetics, Spanish National Cancer Research Centre, Madrid, Spain

Agustin F. Fernandez

Affiliation: Cancer Epigenetics Laboratory, Instituto Universitario de Oncologia del Principado de Asturias, HUCA, Universidad de Oviedo, Oviedo, Spain

Martin Friedrichsen

Affiliation: Steno Diabetes Center, Gentofte, Denmark

Birgitte F. Vind

Affiliation: Department of Endocrinology, Odense University Hospital, Odense, Denmark

Kurt Højlund

Affiliation: Department of Endocrinology, Odense University Hospital, Odense, Denmark

Henning Beck-Nielsen

Affiliation: Department of Endocrinology, Odense University Hospital, Odense, Denmark

Manel Esteller

Affiliation: Cancer Epigenetics and Biology Program (PEBC), Bellvitge Biomedical Research Institute (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain

Allan Vaag

Affiliations Steno Diabetes Center, Gentofte, Denmark, Department of Endocrinology, Rigshospitalet, Copenhagen, Denmark

Pernille Poulsen

Affiliation: Novo Nordisk A/S, Søborg, Denmark

Introduction

Type 2 diabetes (T2D) is a heterogeneous and complex disease resulting from a combination of impaired pancreatic insulin secretion and insulin resistance in tissues such as skeletal muscle, adipose tissue, and liver. The etiology of tissue defects causing T2D is multifactorial. Several gene polymorphisms have been identified [1]–[8], which together with environmental factors such as an adverse fetal environment, aging and obesity increase the risk of T2D [9]–[11].

Phenotypic discordance in monozygotic twins is traditionally attributed to environmental factors distinct for each individual. Therefore, paired analyses of monozygotic twins discordant for a disease phenotype provide an excellent tool for examination of the environmental contribution to the disease in question.

Epigenetics is traditionally referred to as heritable changes in gene expression which are not due to any alteration in DNA sequence. Besides being heritable, epigenetics is altered by environmental factors and hence represents a potential mechanism through which the environment may cause phenotypic variation [12]. Two major classes of epigenetic modifications of the chromatin exist: methylation of cytosine, mostly at cytosine-guanine dinucleotides (CpG)s, and histone modifications, notably acetylations and methylations. The gene promoter regions are rich in CpG sites, forming CpG islands [13], and methylation of these cytosines is thought to silence gene transcription [14], [15]. Histone modifications may result in both activation and silencing of genes [16].

Dietary intervention [17], [18] and exercise [19], [20] have been demonstrated to provoke epigenetic modulation in humans. We previously investigated global differences in methylation of repetitive intergenic DNA sequences and acetylation of histones in young and elderly monozygotic twin pairs. Interestingly, each pair of young twins had essentially similar epigenetic markings, whereas intra-pair differences were substantial in elderly twins, suggesting epigenetic changes to accumulate during life [21]. Therefore, epigenetics has been supposed to play a role in the development of age-related diseases [22]. Data on the role of epigenetics in T2D are still sparse, particularly in the glucose metabolic tissues involved in the pathogenesis of insulin resistant T2D. Increased methylation of PPARGC1A, encoding peroxisome proliferator-activated receptor (PPAR)γ co-activator 1α (PGC-1α), was first reported in pancreatic islets from type 2 diabetic individuals [23], and interestingly, similar T2D-related changes were also found in skeletal muscle [24]. In addition, another gene involved in mitochondrial function, PDK4, has been shown to have increased methylation in skeletal muscle from T2D patients [25]. Among the known genes associated with genetic risk of T2D some have been identified with DNA methylation differences between type 2 diabetic and healthy individuals including INS [26] and PDX1 [27] in pancreatic islets and FTO, SLC30A8 and TCF7L2 in leukocytes [28], [29]. Furthermore, a recent epigenome-wide association study of pancreatic islets from type 2 diabetic and non-diabetic deceased individuals found differential methylation in several genes not previously associated with T2D, including NIBAN and CHAC1 which are involved in endoplasmic reticulum stress [30]. Also individuals genetically predisposed to T2D through a family history of T2D have been shown to have altered DNA methylation in skeletal muscle, among others in mitogen-activated protein kinase pathway genes [20].

Applying an epigenome-wide microarray approach in a unique population of elderly monozygotic twins discordant for T2D, the present study aimed to examine whether likely acquired changes in DNA methylation of gene promoters in skeletal muscle and subcutaneous adipose tissue (SAT) associate with T2D. We hypothesized that the genomic identity of the twins would eliminate genetic causes of T2D-related DNA methylation facilitating the identification of acquired changes in a paired study design. On the other hand, the difficulty of obtaining biopsy material from a larger number of the rare T2D-discordant twin pairs is a limitation of the study.

Results

Subject Characteristics

The study population included 12 Danish monozygotic twin pairs discordant for T2D (Table 1). Skeletal muscle and abdominal SAT biopsies were available from 11 and 5 pairs, respectively. The non-diabetic co-twin had normal glucose tolerance in 2 pairs and impaired glucose tolerance in 10 pairs (Figure 1A). Insulin sensitivity, measured as the glucose infusion rate (GIR), was significantly lower in twins with T2D than in their non-diabetic co-twins (Figure 1B). In addition, fasting plasma glucose, blood hemoglobin A1c (HbA1c) and body mass index (BMI) differed significantly between type 2 diabetic and non-diabetic twins (Table 1).

[Figure omitted. See PDF.]

Figure 1. Discordance for oral glucose tolerance and insulin sensitivity in monozygotic twins. A

Plasma glucose 120 min after a 75-g oral glucose load in type 2 diabetic and non-diabetic twins (n = 12 pairs, P<0.001). B Glucose infusion rate during a euglycemic-hyperinsulinemic clamp in type 2 diabetic and non-diabetic twins (n = 9 pairs, P = 0.006). Data are presented as mean±standard error of the mean.

https://doi.org/10.1371/journal.pone.0051302.g001

[Figure omitted. See PDF.]

Table 1. Characteristics of monozygotic twins discordant for type 2 diabetes.

https://doi.org/10.1371/journal.pone.0051302.t001

Global DNA Methylation Patterns

In skeletal muscle and SAT, the fractions of low methylated (<25%) CpG sites on the array were 64% and 65%, respectively, whereas the fractions of highly methylated (>75%) sites were 13% and 15%. Twins from a pair were relatively similar regarding DNA methylation in each tissue type (Figure 2A−B), whereas the intra-individual methylation difference between skeletal muscle and SAT was considerable (Figure 2C). The intra-pair methylation differences were slightly larger in skeletal muscle than in SAT. Twin pair was a significant predictor of overall methylation pattern in both skeletal muscle (P<0.001) and SAT (P<0.001), whereas overall methylation did not differ according to diabetes status (P>0.1). However, in both muscle (P<0.001) and SAT (P = 0.009), twin pair interacted with diabetes status. The variation, expressed as the standard deviation (SD) of absolute intra-pair methylation differences, was significantly greater in repetitive LINE1, D4Z4 and NBL2 DNA sequences than in promoter regions (SD: 10% points vs. 4% points, P<0.001), but there was no correlation between the degree of variation in these regions (r = −0.05, P = 0.9, Figure 3A). The largest variation was found in methylation of LINE1 CpG sites (SD: 13% points), and it correlated positively with the intra-pair differences in BMI (r = 0.77, P = 0.007, Figure 3B) and plasma glucose 2 hours after an oral glucose tolerance test (OGTT, r = 0.66, P = 0.03, Figure 3C), but not significantly with the intra-pair difference in GIR (r = 0.52, P = 0.2, Figure 3D). The intra-pair variation in D4Z4, NBL2 and promoter regions did not correlate with phenotypic differences in BMI, 2-hour glucose and GIR.

[Figure omitted. See PDF.]

Figure 2. Global gene promoter DNA methylation in monozygotic twins.

The DNA methylation was measured as the β-value ranging from 0 (unmethylated) to 1 (completely methylated) at 26,850 CpG sites located in the promoters of 14,279 genes. The plots are shown for a representative twin pair (#3). A Skeletal muscle (r = 0.95, P<0.001). B SAT (r = 0.97, P<0.001). C Comparison of DNA methylation in SAT and skeletal muscle from the non-diabetic twin.

https://doi.org/10.1371/journal.pone.0051302.g002

[Figure omitted. See PDF.]

Figure 3. Intra-twin pair variation of DNA methylation in repetitive DNA sequences.

The variation is given by the standard deviation (SD) of the absolute methylation differences between monozygotic twins in skeletal muscle. A SD of methylation differences on all repetitive versus all promoter CpG sites (r = −0.05, P = 0.9). B SD of methylation differences on LINE1 repetitive CpG sites versus absolute difference in BMI (r = 0.77, P = 0.007). C SD of methylation differences on LINE1 repetitive CpG sites versus absolute difference in 2-hour plasma glucose (r = 0.66, P = 0.03). D SD of methylation differences on LINE1 repetitive CpG sites versus absolute difference in glucose infusion rate (r = 0.52, P = 0.2).

https://doi.org/10.1371/journal.pone.0051302.g003

Differentially Methylated CpG Sites in Type 2 Diabetes

In the 49 known susceptibility genes for mono- or polygenetic T2D represented on the microarray, 136 CpG sites were examined for differential methylation by a candidate approach. Among these sites, 25 sites located in 17 genes (8 sites in muscle and 17 sites in SAT) were differentially methylated between type 2 diabetic and non-diabetic twins (Table 2, Figure S1B,D–E). These candidate genes were CDKN2A, DUSP9, HNF4A, HHEX, KCNQ1, KLF11, PPARGC1A and SLC30A8 in muscle and ADCY5, CAV1, CIDEC, CDKN2A, CDKN2B, DUSP9, HNF4A, IDE, IRS1, KCNQ1, MTNR1B, TSPAN8 and WFS1 in SAT. Two CpG sites in SAT (CDKN2A and HNF4A) remained statistically significant after permutation correction for multiple testing (Padj = 0.02; Table 2, Figure S1B,D–E). An explorative analysis of all 26,850 array CpG sites showed that 789 (3%) sites in muscle and 1,458 (5%) sites in SAT were differentially methylated (P<0.05) in type 2 diabetic compared with non-diabetic twins before correction for multiple testing. These CpG sites were located in 768 genes in muscle and 1,389 genes in SAT. One CpG site in muscle (IL8) and 7 sites in SAT (ZNF668, HSPA2, C8orf31, CD320, SFT2D3, TWIST1, MYO5A) remained statistically significant after permutation correction (Padj<0.001; Table 3). None of the CpG sites in repetitive LINE1, D4Z4 and NBL2 sequences differed significantly in methylation between type 2 diabetic and non-diabetic twins (P>0.1).

[Figure omitted. See PDF.]

Table 2. Differentially methylated type 2 diabetes susceptibility genes.

https://doi.org/10.1371/journal.pone.0051302.t002

[Figure omitted. See PDF.]

Table 3. Differentially methylated genes.

https://doi.org/10.1371/journal.pone.0051302.t003

Validation of Array Data

We have recently validated the Illumina DNA methylation array’s ability to detect methylation differences of approximately 10% points in paired skeletal muscle biopsies from an overfeeding intervention study [18]. In the present study, five genes (three in SAT and two in muscle) with larger absolute methylation differences (>4% points) were selected for validation of the microarray data in three twin pairs. Bisulfite sequencing (BS) was used in SAT for MCF2, HNF4A and FAP, and bisulfite pyrosequencing (PBS) was used in muscle for PPARGC1A and SLC30A8. The methylation differences for these genes were successfully validated (Figure S1).

Pathway Analyses

Ingenuity Pathway Analysis software was used to identify molecular pathways with a significant proportion of genes having CpG sites differentially methylated between type 2 diabetic and non-diabetic twins. The analysis was done on genes with at least one differentially methylated GpC site based on the uncorrected P-values. Nine muscle pathways were identified, including inflammatory (hepatic fibrosis and IL-6), lipid metabolic (PPARα, PPARγ and sphingolipids) and carbohydrate metabolic (pyruvate and propanoate) pathways (Table S2). In SAT, one carbohydrate metabolism and one circadian rhythm signaling pathway had a significant fraction of genes with methylation changes (Table S2).

Discussion

We have previously used a population of twelve monozygotic twin pairs discordant for T2D to recognize a considerable non-genetic contribution to glucose metabolic disturbances in genetically identical individuals [31]. In the present study, we used a similar twin sample for the investigation of genome-wide common acquired changes in gene promoter DNA methylation in the primary insulin responsive tissues, skeletal muscle and SAT. We found that absolute intra-twin pair methylation differences were relatively small, and most of them did not link significantly to T2D after correction for multiple testing. The findings altogether suggest a modest contribution of acquired DNA methylation differences in skeletal muscle or SAT to the non-genetic component of T2D.

It has been demonstrated in a prospective study that DNA methylation changes occur over time, suggesting an influence by environmental or stochastic events [12]. A study of monozygotic and dizygotic twins has provided evidence that such epigenetic changes are more pronounced in non-CpG island than in CpG island DNA regions [32]. In addition, a clustering of DNA methylation changes in families indicated that the susceptibility to a given epigenetic change may have a genetic component [12]. We hypothesized that the use of genetically identical twins with its inherent correction for the genetic influence on both disease phenotype and DNA methylation would represent the ideal design to evaluate acquired DNA methylation changes associated with T2D.

We investigated the methylation percentage of 26,850 gene promoter CpG sites in skeletal muscle and SAT and found that the majority of these sites were clustered into two groups with either a low or a high methylation percentage. This indicates a high concordance of methylation for the DNA copies present in a biopsy. Furthermore, co-twins showed a greater degree of similarity in gene promoter methylation than unrelated individuals. The overall promoter methylation did not depend on T2D status alone, but interestingly it was significantly associated with the interaction between T2D status and twin pair. This suggests that, within each pair, T2D-related DNA methylation differences might exist.

The intra-individual DNA methylation pattern in gene promoters differed vastly between tissues, proposing that DNA methylation changes associated with tissue development greatly exceed changes occurring in the finally differentiated tissue. Thus, our data support the notion of a major role of genetics as well as tissue specificity in determining the methylation pattern of gene promoters in adults. The intra-pair methylation differences were slightly larger in skeletal muscle than in SAT. This phenomenon could reflect a greater environmental influence on DNA methylation in skeletal muscle, but could also be due to a greater diversity of cell types in this tissue.

In contrast to gene promoters, we found a considerable intra-pair variation in methylation of genome-wide repetitive DNA sequences. We have previously demonstrated that differences in DNA methylation of non-coding regions are more prominent in elderly monozygotic twins who have spent less of their lives together than those having shared the environment for a longer period [21]. A novel finding of our present study is that the variation in methylation levels of repetitive LINE1 DNA, which makes up 17% of the human genome [33], was largest in twin pairs being different in BMI and 2-hour plasma glucose. This result provides evidence that epigenetic variation reflects the degree of dissimilarity between phenotypes of monozygotic twins. However, in agreement with a recent study of human pancreatic islet cells [30] we found no specific methylation changes in repetitive DNA sequences between type 2 diabetic and non-diabetic individuals.

Despite being relatively small, a number of the intra-pair methylation differences in gene promoters identified by the microarrays were successfully validated by two different bisulfite sequencing methods. The validation additionally showed that CpG sites close to the array site were likely to be co-methylated supporting previous observations [14], [17]. Although this technical validation does not exclude the modest biological differences to have occurred by chance, it supports the sensitivity of the microarray to detect methylation differences down to a few percentage points. Previous studies have reported large methylation differences between cancer cells and normal cells [15], which may among others be attributed to the fact that the tumor comprises a clone of abnormal cells. In the, by comparison, healthier and non-dividing muscle and SAT tissues of type 2 diabetic individuals, it could be speculated that environmentally determined methylation changes may take place in only some cells. Recent studies of IGF2 [34], HNF4A [35] and PPARGC1A [17], [18] methylation in blood or skeletal muscle from individuals exposed to an adverse intrauterine environment showed that methylation differences between exposed and unexposed individuals were of a similar size to those found in our study. This was also the case for reported T2D-associated methylation differences of PPARGC1A [23], INS [26] and PDK4 [25] in pancreatic islets or skeletal muscle. mRNA expression of these genes was altered in T2D and in some cases correlated with DNA methylation [23], [26]. However, epigenome-wide DNA methylation profiling studies in pancreatic islets [30] and skeletal muscle [18] with expression analysis of a fraction of the differentially methylated genes show that the relationship between small DNA methylation differences and gene expression is complicated. Several factors such as location of the specific DNA methylation site in relation to gene regulatory regions and co-existence of other regulatory mechanisms complicate the study of smaller methylation differences’ biological impact.

Most known T2D susceptibility gene polymorphisms associate primarily with β-cell dysfunction making the endocrine pancreas the likely tissue for T2D-associated DNA methylation changes. A recent genome-wide study of T2D-associated DNA methylation differences in a small sample of pancreatic islet biopsies identified many differentially methylated genes in T2D, but none of the known susceptibility genes were among these [30]. IRS1 [8] and PPARGC1A [2] belong to the small group of T2D susceptibility genes which might mediate their effects through insulin resistance and therefore be relevant in peripheral tissues such as SAT and skeletal muscle. When the known T2D susceptibility genes represented on the array were specifically evaluated, we found that ADCY5, CAV1, CIDEC, CDKN2A, CDKN2B, DUSP9, HNF4A, IDE, IRS1, KCNQ1, MTNR1B, TSPAN8 and WFS1 had at least one CpG site differentially methylated between type 2 diabetic and non-diabetic twins in SAT. In skeletal muscle, CpG sites in CDKN2A, DUSP9, HNF4A, HHEX, KCNQ1, KLF11, PPARGC1A and SLC30A8 were differentially methylated. However, only the maturity-onset diabetes of the young (MODY) gene HNF4A [36] and the T2D gene CDKN2A [37] in SAT had significant permutation adjusted P-values. Importantly, the absolute methylation differences identified were smaller than those predicted to be found at 80% power, particularly in the adipose tissue samples obtained from only five twin pairs. This finding contrasts the methylation differences of ∼20–30% points identified in a similarly sized sample of pancreatic islets from type 2 diabetic and non-diabetic individuals [30] and could be interpreted as the endocrine pancreas being indeed the tissue with most pathological changes in T2D. In addition to the statistical limitation, some of the differences on low-methylated CpG sites such as in IRS1 had absolute methylation differences that are probably too small to be detected reliably by the DNA methylation microarray. The finding of increased HNF4A methylation in tissues from type 2 diabetic individuals is interesting considering the recent similar result in umbilical cord blood leukocytes from newborns with intrauterine growth restriction [35]. Thus, methylation of this gene could represent a mechanism linking early intrauterine programming to development of T2D later in life. Even though the increased methylation of PPARGC1A in type 2 diabetic individuals was not significant after correction for multiple testing, the result corresponds well with previous findings in β-cells [23] and skeletal muscle from type 2 diabetic patients [24]. To this end, we recently reported elevated methylation of the PPARGC1A promoter in skeletal muscle from lean and otherwise healthy young men born with low birth weight. Moreover, five days of high-fat, high-calorie feeding increased PPARGC1A promoter methylation in healthy men with normal birth weight [17]. Thus, increased methylation of PPARGC1A in skeletal muscle seems to be a consistent finding in patients with overt T2D as well as in individuals at risk of developing T2D. The recently discovered T2D susceptibility gene, DUSP9, is another gene which could be important to insulin resistance. It encodes mitogen-activated kinase phosphatase-4 which counteracts stress-induced insulin resistance [38]. Interestingly, this gene had methylation changes of 8−9% points dependent on tissue which was among the largest found on the array.

The combination of small absolute methylation differences and a relatively small sample size made it difficult to identify statistically significant associations between methylation changes and T2D by an explorative approach. We found 789 and 1,458 of 26,850 CpG sites which were differentially methylated in skeletal muscle and SAT, respectively. These approximately 5% of the CpG sites investigated could be chance findings, and the fact that only a few CpG sites remained significant after correction for multiple testing contributes to this interpretation of the results. Genes for which the permutation corrected P-values were significant included IL8 in muscle and TWIST1 in SAT. The inflammatory cytokine IL-8 is expressed in among others skeletal muscle in response to exercise [39], and the transcription factor Twist1 regulates expression of inflammatory cytokines in adipocytes [40].

Molecular pathway analyses were used to examine whether the 768 muscle and 1,458 SAT genes which, based on the explorative analysis, were most likely to have DNA methylation changes belonged to specific functional pathways. Since the methylation differences in the vast majority of these genes were not significant after correction for multiple testing the results should be interpreted cautiously. We found the genes in muscle to be predominantly involved in inflammation, lipid metabolism, and carbohydrate metabolism. Low-grade inflammation is considered to be an important mechanism in insulin resistance [41], and the lipid metabolic regulators PPARα, PPARγ, as well as PGC-1α have been shown to play a role in the regulation of insulin sensitivity [42]–[44]. In SAT, carbohydrate metabolism and circadian rhythm signaling pathways had a significant fraction of genes differentially methylated between type 2 diabetic and non-diabetic individuals. Recent studies have shown that genetic variation in MTNR1B and CRY2, involved in circadian signaling, associate with T2D or related metabolic traits [3], [45]. Interestingly, both of these genes were differentially methylated in SAT from type 2 diabetic and non-diabetic twins.

In conclusion, common skeletal muscle and adipose tissue gene promoter DNA methylation differences between monozygotic twins discordant for T2D were small, but a number of the differences were found in known T2D-related candidate genes. Larger study groups, preferably with longitudinal sample collection and methylation analysis of the entire promoter are needed to replicate the findings in these genes.

Research Design and Methods

Study Participants

Twelve 53–80 year-old monozygotic twin pairs discordant for T2D were recruited through the Danish Twin Registry, University of Southern Denmark. Six pairs had participated in a previous study [46] where discordance for T2D was recognized based on an OGTT. These pairs were reexamined for the present study. The additional six pairs were recruited based on information from the Twin Registry about known T2D. All study participants had provided written informed consent, and the study was approved by the regional Ethical Committee (Southern Denmark, http://komite.regionsyddanmark.dk) and conducted in accordance with the principles of the Helsinki Declaration.

Clinical Examination

Height and weight were measured in light weight clothes for calculation of BMI. Discordance for T2D was verified by a 75-g OGTT. Insulin sensitivity was measured in a subpopulation (n = 9 pairs) by a euglycemic-hyperinsulinemic clamp (40 mU m−2 min−1) and expressed as the mean GIR during the last 30 min of the clamp period.

Biopsies

Biopsies were excised under local anesthesia (lidocaine) from the vastus lateralis muscle (n = 11 pairs) and from abdominal SAT (n = 5 pairs) using a Bergström needle with suction applied. The tissue samples were frozen immediately in liquid nitrogen and stored at −80°C until processed further.

DNA Extraction

Genomic DNA was extracted from the biopsies using a DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany), and DNA concentrations were determined by a Quant-IT PicoGreen dsDNA Kit (Invitrogen, Carlsbad, CA).

DNA Methylation Microarrays

Methylation was assessed at 27,578 CpG sites primarily located close to the transcription start site (TSS) in the promoters of 14,475 genes using 12-sample Infinium HumanMethylation27 Bead Chips (Illumina, San Diego, CA). 600 ng DNA was bisulfite-treated to deaminate unmethylated cytosines to uracil using an EZ DNA Methylation Kit (Zymo Research, Orange, CA). The array was scanned by a BeadArray Reader (Illumina), and intensity data analyzed using GenomeStudio software (version 2011.1, Illumina). Internal array controls verified the hybridization, staining and washing procedures. The methylation status on each CpG site was expressed as the β-value which is the ratio between fluorescent signal from converted and preserved sequence bead types. DNA samples from twin pairs were analyzed on the same bead chip to eliminate a possible influence of batch effect in the subsequent paired statistical analyses. Due to reruns four pairs of muscle samples were separated on different chips. Batch effect in the data set was examined for by unsupervised cluster analysis using MeV software (version 4.5, available at http://www.tm4.org) and scatter plots of principal components 1 and 2. Since only minimal batch effect was found compared to clustering according to twin pair and sex, no further normalization of the β-values was performed. β-values for probes where intensity did not exceed the background level (P≥0.05) were omitted (in average 54 β-values per muscle sample and 23 β-values per SAT sample). All probes were examined for unique genome alignment and for SNPs affecting the CpG site using NCBI human genome build 36 FASTA files and custom PERL (version 5.10.1, available at http://www.perl.org) scripts. For the analysis of unique alignment, two “bisulfite-treated” reference genomes were constructed replacing cytosine by thymine as forward strand reference, and replacing guanine by adenine as reverse strand reference. The unmethylated bead type probe sequences were compared to these reference genomes using the BLAT algorithm [47]. For the SNP analysis, Illumina’s annotation data were used to mark each target region in a reference genome containing dbSNP build 130 information. All corresponding SNP-masked sequences were extracted and screened for SNPs. 494 probes were excluded due to multiple alignments, and additionally 234 were excluded due to SNP in the CpG site. The microarray data have been deposited in the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo, accession number GSE38291) complying with the Minimum Information About a Microarray Experiment (MIAME) guidelines.

Validation and Methylation of Repetitive DNA

Array data were validated in three genes in SAT by BS and in two genes in skeletal muscle by BPS. DNA methylation of LINE1 (interspersed repeat with 8 CpG sites), D4Z4 (tandem repeat with 9 CpG sites) and NBL2 (tandem repeat with 8 CpG sites) sequences was measured by BPS in the muscle samples. The repetitive distribution of these sequences means that the DNA methylation measured is a genome-wide average value.

Bisulfite sequencing.

1 µg DNA was bisulfite-treated as described previously [48]. Oligonucleotide primers (Sigma-Aldrich, St. Louis, MO) were designed using Methyl Primer Express Software (Applied Biosystems, Foster City, CA) to make the PCR amplicons cover approximately five CpG sites, including the site analyzed on the array (Table S1). The PCR product was separated on a 2% agarose gel, and specific DNA bands were cut and purified by a GFX PCR DNA and Gel Band Purification Kit (GE Healthcare, Buckinghamshire, UK). The specific amplicons were cloned in E. coli using the pGEM-T vector system (Promega, Madison, WI). Twelve positive colonies were collected, and plasmid DNA was purified by a Perfectprep Plasmid 96 Vac Kit (Eppendorf, Hamburg, Germany). The plasmid DNA insert was finally sequenced on an ABI 3100 system (Applied Biosystems). The number of methylated and unmethylated clones was counted, and the methylation percentage calculated as the average methylation of all CpG sites in the amplicon.

Bisulfite pyrosequencing.

BPS was performed on bisulfite-treated DNA using a PyroMark MD pyrosequencing system (Biotage, Uppsala, Sweden). Specific pyrosequencing primers were designed to amplify as many CpGs as conditions permitted (1 to 9 sites) using Assay Design Software (version 1.0.6, Biotage, Table S1). The cytosine methylation percentage was evaluated with the Pyro Q-CpG program (version 1.0.9, Biotage).

Statistical Methods

The statistical analyses were performed using R software (version 2.11.1, available at http://www.r-project.org). β-values are shown in percentages and differences given in percentage point change (type 2 diabetic twin − non-diabetic twin). The differences were examined for normality using the Shapiro-Wilk normality test. Only 6% of the CpG sites in adipose tissue and 9% of the CpG sites in skeletal muscle had P-values <0.05 suggestive of non-normally distributed differences, and therefore parametric statistics were applied.

Analysis of twin pair similarity.

An ANOVA including methylation of all microarray CpG sites as response variable and twin pair, diabetes status and their interaction as explanatory variables was used to evaluate the effect of twin pairs on the overall DNA methylation similarity. For correlation analyses the intra-pair variation in promoter or repetitive DNA methylation was expressed as the SD for the absolute differences.

Analysis of differential methylation by candidate gene approach.

All monogenic and confirmed common variety T2D susceptibility genes reviewed by O’Rahilly [1] were examined for differential methylation if represented on the array. In addition, KCNQ1 [6] and KLF11 [7] as well as the newly discovered genes, ADCY5, GCK, GCKR, PROX1, DUSP9, HMGA2, KLF14, TP53INP1 and RBMS1 [3]–[5] were included. PPARGC1A [2] was included due to an a priori hypothesis of altered DNA methylation in type 2 diabetic individuals [23], [24]. Thus, the genes analyzed were ABCC8, ADAMTS9, ADCY5, AGPAT2, AKT2, BSCL2, CAMK1D, CAV1, CDC123, CDKAL1, CDKN2A, CDKN2B, CEL, CIDEC, DUSP9, HHEX, HMGA2, HNF1A, HNF1B, HNF4A, GCK, GCKR, IDE, IGF2BP2, INS, INSR, IRS1, JAZF1, KCNJ11, KCNQ1, KIF11, KLF11, KLF14, LGR5, LMNA, MTNR1B, NEUROD1, NOTCH2, PDX1, PPARG, PPARGC1A, PROX1, RBMS1, SLC30A8, TBC1D4, TCF7L2, TP53INP1, TSPAN8, and WFS1. These 49 genes were represented by 136 probes. Comparisons between type 2 diabetic and non-diabetic twins were performed by paired t-tests. Two-sided P<0.05 was considered statistically significant, and the P-values were corrected for multiple testing, Padj, using the Westfall-Young resampling method [49]. 10,000 simulations with permutations of the sample labels were used to sample the null distribution.

Analysis of differential methylation by explorative approach.

Comparisons between type 2 diabetic and non-diabetic twins were performed for all 26,850 array CpG sites by paired t-tests, and the P-values were corrected for multiple testing as described above.

Statistical power calculations.

The study population size was limited to the maximal number of Danish T2D discordant monozygotic twin pairs that could be recruited in a 10-year period. From previous studies of PPARGC1A, DNA methylation differences of ∼5% points with ∼5% points SD have been identified [17], [23]. Given the 11 pairs of muscle biopsies and a fixed SD of 5% points the study allowed detecting a 5% point methylation difference in a single paired t-test at 80% power. The size of the minimal detectable differences would increase to 10% in the candidate gene analysis (136 t-tests) and 17% in the genome-wide analysis (26,850 t-tests) when applying Bonferroni correction. For the only 5 pairs of adipose tissue biopsies the similar minimal detectable differences would be 8%, 31% and 117% for single CpG site analysis, candidate gene study and hypothesis-free approach, respectively.

Validation of methylation array data.

Differences in DNA methylation levels between three pairs of type 2 diabetic and non-diabetic twins obtained by BS or BPS were analyzed by paired t-test. Based on the a priori hypothesis of replicating array findings one-sided P<0.05 was considered statistically significant.

Molecular pathway analyses.

Differentially methylated genes identified by the explorative approach, without correction for multiple testing, were included in Ingenuity Pathway Analyses (version 7.5, Ingenuity Systems, Redwood City, CA). This software recognized the pathways from the Ingenuity Pathways Analysis library of canonical pathways (81 metabolic and 283 signaling) which contained a significant fraction of differentially methylated genes. Only genes represented on the array were used for the reference pathways. Fischer’s exact test was used to calculate the probability that the association between the differentially methylated genes and the canonical pathway was explained by chance. These P-values were not corrected for multiple testing.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Validation of array data. The average DNA methylation in the gene promoter area surrounding the array site measured by bisulfite sequencing (BS) or bisulfite pyrosequencing (BPS) is indicated in type 2 diabetic (black bars) and non-diabetic (white bars) twins. Three twin pairs were included in the validations. One-sided P-values are shown for the difference between type 2 diabetic and non-diabetic twins measured by BS or BPS. Gene diagrams, including the number of methylated (black) and unmethylated (white) clones for each CpG site in the amplicon, are shown for BS results. A MCF2 +226 base pairs from transcription start site (TSS) in subcutaneous adipose tissue (P<0.001), B HNF4A −521 base pairs from TSS in subcutaneous adipose tissue (P = 0.04), C FAP +80 base pairs from TSS in subcutaneous adipose tissue (P<0.001), D PPARGC1A −383 base pairs from TSS in skeletal muscle (P = 0.006), E SLC30A8 −174 base pairs from TSS in skeletal muscle (P = 0.002). Data are presented as mean±standard error of the mean.

https://doi.org/10.1371/journal.pone.0051302.s001

(PDF)

Table S1.

Primer sequences.

https://doi.org/10.1371/journal.pone.0051302.s002

(DOC)

Table S2.

Molecular pathways with a significant fraction of differentially methylated genes. No overlap denotes the fraction of genes in each pathway which was not represented on the array. The direction of methylation change in type 2 diabetic versus non-diabetic twins is indicated by arrows (↑: increased, ↓: decreased).

https://doi.org/10.1371/journal.pone.0051302.s003

(DOC)

Acknowledgments

We thank the twins for their participation in the study and M. Modest, L. Moreno and R. Alonso for excellent assistance in the laboratory.

Author Contributions

Conceived and designed the experiments: PP MFF AV KH HB-N. Performed the experiments: RR-M AFF MFF EL VC SJ BFV PP AV. Analyzed the data: RR-M JB-J AFF MF. Contributed reagents/materials/analysis tools: MFF AFF ME. Wrote the paper: RR-M PP AV.

Citation: Ribel-Madsen R, Fraga MF, Jacobsen S, Bork-Jensen J, Lara E, Calvanese V, et al. (2012) Genome-Wide Analysis of DNA Methylation Differences in Muscle and Fat from Monozygotic Twins Discordant for Type 2 Diabetes. PLoS ONE 7(12): e51302. https://doi.org/10.1371/journal.pone.0051302

References

1. O’Rahilly S (2009) Human genetics illuminates the paths to metabolic disease. Nature 462: 307–314.

2. Ek J, Andersen G, Urhammer SA, Gaede PH, Drivsholm T, et al. (2001) Mutation analysis of peroxisome proliferator-activated receptor-gamma coactivator-1 (PGC-1) and relationships of identified amino acid polymorphisms to Type II diabetes mellitus. Diabetologia 44: 2220–2226.

3. Dupuis J, Langenberg C, Prokopenko I, Saxena R, Soranzo N, et al. (2010) New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 42: 105–116.

4. Qi L, Cornelis MC, Kraft P, Stanya KJ, Linda Kao WH, et al. (2010) Genetic variants at 2q24 are associated with susceptibility to type 2 diabetes. Hum Mol Genet 19: 2706–2715.

5. Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, et al. (2010) Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet 42: 579–589.

6. Yasuda K, Miyake K, Horikawa Y, Hara K, Osawa H, et al. (2008) Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus. Nat Genet 40: 1092–1097.

7. Neve B, Fernandez-Zapico ME, Ashkenazi-Katalan V, Dina C, Hamid YH, et al. (2005) Role of transcription factor KLF11 and its diabetes-associated gene variants in pancreatic beta cell function. Proc Natl Acad Sci U S A 102: 4807–4812.

8. Rung J, Cauchi S, Albrechtsen A, Shen L, Rocheleau G, et al. (2009) Genetic variant near IRS1 is associated with type 2 diabetes, insulin resistance and hyperinsulinemia. Nat Genet 41: 1110–1115.

9. Hales CN, Barker DJ, Clark PM, Cox LJ, Fall C, et al. (1991) Fetal and infant growth and impaired glucose tolerance at age 64. BMJ 303: 1019–1022.

10. Fink RI, Kolterman OG, Griffin J, Olefsky JM (1983) Mechanisms of insulin resistance in aging. J Clin Invest 71: 1523–1535.

11. de Vegt F, Dekker JM, Jager A, Hienkens E, Kostense PJ, et al. (2001) Relation of impaired fasting and postload glucose with incident type 2 diabetes in a Dutch population: The Hoorn Study. JAMA 285: 2109–2113.

12. Bjornsson HT, Sigurdsson MI, Fallin MD, Irizarry RA, Aspelund T, et al. (2008) Intra-individual change over time in DNA methylation with familial clustering. JAMA 299: 2877–2883.

13. Bird A, Taggart M, Frommer M, Miller OJ, Macleod D (1985) A fraction of the mouse genome that is derived from islands of nonmethylated, CpG-rich DNA. Cell 40: 91–99.

14. Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, et al. (2006) DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet 38: 1378–1385.

15. Calvanese V, Horrillo A, Hmadcha A, Suarez-Alvarez B, Fernandez AF, et al. (2008) Cancer genes hypermethylated in human embryonic stem cells. PLoS ONE 3: e3294.

16. Cheung P, Allis CD, Sassone-Corsi P (2000) Signaling to chromatin through histone modifications. Cell 103: 263–271.

17. Brons C, Jacobsen S, Nilsson E, Rönn T, Jensen CB, et al. (2010) Deoxyribonucleic acid methylation and gene expression of PPARGC1A in human muscle is influenced by high-fat overfeeding in a birth-weight-dependent manner. J Clin Endocrinol Metab 95: 3048–3056.

18. Jacobsen SC, Brons C, Bork-Jensen J, Ribel-Madsen R, Yang B, et al. (2012) Effects of short-term high-fat overfeeding on genome-wide DNA methylation in the skeletal muscle of healthy young men. Diabetologia 55: 3341–3349.

19. Barres R, Yan J, Egan B, Treebak JT, Rasmussen M, et al. (2012) Acute exercise remodels promoter methylation in human skeletal muscle. Cell Metab 15: 405–411.

20. Nitert MD, Dayeh T, Volkov P, Elgzyri T, Hall E, et al.. (2012) Impact of an Exercise Intervention on DNA Methylation in Skeletal Muscle From First-Degree Relatives of Patients With Type 2 Diabetes. Diabetes: In press doi: 10.2337/db11-1653.

21. Fraga MF, Ballestar E, Paz MF, Ropero S, Setien F, et al. (2005) Epigenetic differences arise during the lifetime of monozygotic twins. Proc Natl Acad Sci U S A 102: 10604–10609.

22. Calvanese V, Lara E, Kahn A, Fraga MF (2009) The role of epigenetics in aging and age-related diseases. Ageing Res Rev 8: 268–276.

23. Ling C, Del GS, Lupi R, Rönn T, Granhall C, et al. (2008) Epigenetic regulation of PPARGC1A in human type 2 diabetic islets and effect on insulin secretion. Diabetologia 51: 615–622.

24. Barres R, Osler ME, Yan J, Rune A, Fritz T, et al. (2009) Non-CpG methylation of the PGC-1alpha promoter through DNMT3B controls mitochondrial density. Cell Metab 10: 189–198.

25. Kulkarni SS, Salehzadeh F, Fritz T, Zierath JR, Krook A, Osler ME (2012) Mitochondrial regulators of fatty acid metabolism reflect metabolic dysfunction in type 2 diabetes mellitus. Metabolism 61: 175–185.

26. Yang BT, Dayeh TA, Kirkpatrick CL, Taneera J, Kumar R, et al. (2011) Insulin promoter DNA methylation correlates negatively with insulin gene expression and positively with HbA(1c) levels in human pancreatic islets. Diabetologia 54: 360–367.

27. Yang BT, Dayeh TA, Volkov PA, Kirkpatrick CL, Malmgren S, et al. (2012) Increased DNA methylation and decreased expression of PDX-1 in pancreatic islets from patients with type 2 diabetes. Mol Endocrinol 26: 1203–1212.

28. Bell CG, Finer S, Lindgren CM, Wilson GA, Rakyan VK, et al. (2010) Integrated genetic and epigenetic analysis identifies haplotype-specific methylation in the FTO type 2 diabetes and obesity susceptibility locus. PLoS ONE 5: e14040.

29. Toperoff G, Aran D, Kark JD, Rosenberg M, Dubnikov T, et al. (2012) Genome-wide survey reveals predisposing diabetes type 2-related DNA methylation variations in human peripheral blood. Hum Mol Genet 21: 371–383.

30. Volkmar M, Dedeurwaerder S, Cunha DA, Ndlovu MN, Defrance M, et al. (2012) DNA methylation profiling identifies epigenetic dysregulation in pancreatic islets from type 2 diabetic patients. EMBO J 31: 1405–1426.

31. Vaag A, Henriksen JE, Madsbad S, Holm N, Beck-Nielsen H (1995) Insulin secretion, insulin action, and hepatic glucose production in identical twins discordant for non-insulin-dependent diabetes mellitus. J Clin Invest 95: 690–698.

32. Kaminsky ZA, Tang T, Wang SC, Ptak C, Oh GH, et al. (2009) DNA methylation profiles in monozygotic and dizygotic twins. Nat Genet 41: 240–245.

33. Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, et al. (2001) Initial sequencing and analysis of the human genome. Nature 409: 860–921.

34. Heijmans BT, Tobi EW, Stein AD, Putter H, Blauw GJ, et al. (2008) Persistent epigenetic differences associated with prenatal exposure to famine in humans. Proc Natl Acad Sci U S A 105: 17046–17049.

35. Einstein F, Thompson RF, Bhagat TD, Fazzari MJ, Verma A, et al. (2010) Cytosine methylation dysregulation in neonates following intrauterine growth restriction. PLoS ONE 5: e8887.

36. Yamagata K, Furuta H, Oda N, Kaisaki PJ, Menzel S, et al. (1996) Mutations in the hepatocyte nuclear factor-4alpha gene in maturity-onset diabetes of the young (MODY1). Nature 384: 458–460.

37. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, et al. (2007) Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 316: 1331–1336.

38. Emanuelli B, Eberle D, Suzuki R, Kahn CR (2008) Overexpression of the dual-specificity phosphatase MKP-4/DUSP-9 protects against stress-induced insulin resistance. Proc Natl Acad Sci U S A 105: 3545–3550.

39. Buford TW, Cooke MB, Willoughby DS (2009) Resistance exercise-induced changes of inflammatory gene expression within human skeletal muscle. Eur J Appl Physiol 107: 463–471.

40. Pettersson AT, Laurencikiene J, Mejhert N, Naslund E, Bouloumie A, et al. (2010) A possible inflammatory role of twist1 in human white adipocytes. Diabetes 59: 564–571.

41. Wellen KE, Hotamisligil GS (2005) Inflammation, stress, and diabetes. J Clin Invest 115: 1111–1119.

42. Guerre-Millo M, Gervois P, Raspe E, Madsen L, Poulain P, et al. (2000) Peroxisome proliferator-activated receptor alpha activators improve insulin sensitivity and reduce adiposity. J Biol Chem 275: 16638–16642.

43. Zierath JR, Ryder JW, Doebber T, Woods J, Wu M, et al. (1998) Role of skeletal muscle in thiazolidinedione insulin sensitizer (PPARgamma agonist) action. Endocrinology 139: 5034–5041.

44. Ling C, Poulsen P, Carlsson E, Ridderstrale M, Almgren P, et al. (2004) Multiple environmental and genetic factors influence skeletal muscle PGC-1alpha and PGC-1beta gene expression in twins. J Clin Invest 114: 1518–1526.

45. Bouatia-Naji N, Bonnefond A, Cavalcanti-Proenca C, Sparso T, Holmkvist J, et al. (2009) A variant near MTNR1B is associated with increased fasting plasma glucose levels and type 2 diabetes risk. Nat Genet 41: 89–94.

46. Poulsen P, Grunnet LG, Pilgaard K, Storgaard H, Alibegovic A, et al. (2009) Increased risk of type 2 diabetes in elderly twins. Diabetes 58: 1350–1355.

47. Kent WJ (2002) BLAT–the BLAST-like alignment tool. Genome Res 12: 656–664.

48. Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB (1996) Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci U S A 93: 9821–9826.

49. Westfall PH, Young SS (1993) Resampling-based multiple testing: examples and methods for P-value adjustment. New York: John Wiley & Sons, Inc.