Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic Effects on DNA Methylation and Its Potential Relevance for Obesity in Mexican Americans

  • Melanie A. Carless ,

    mcarless@txbiomedgenetics.org

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Hemant Kulkarni,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Mark Z. Kos,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Jac Charlesworth,

    Affiliations Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America, Menzies Research Institute Tasmania, University of Tasmania, Hobart, Tasmania, Australia

  • Juan M. Peralta,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Harald H. H. Göring,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Joanne E. Curran,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Laura Almasy,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Thomas D. Dyer,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Anthony G. Comuzzie,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • Michael C. Mahaney,

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

  • John Blangero

    Affiliation Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas, United States of America

Abstract

Several studies have identified effects of genetic variation on DNA methylation patterns and associated heritability, with research primarily focused on Caucasian individuals. In this paper, we examine the evidence for genetic effects on DNA methylation in a Mexican American cohort, a population burdened by a high prevalence of obesity. Using an Illumina-based platform and following stringent quality control procedures, we assessed a total of 395 CpG sites in peripheral blood samples obtained from 183 Mexican American individuals for evidence of heritability, proximal genetic regulation and association with age, sex and obesity measures (i.e. waist circumference and body mass index). We identified 16 CpG sites (∼4%) that were significantly heritable after Bonferroni correction for multiple testing and 27 CpG sites (∼6.9%) that showed evidence of genetic effects. Six CpG sites (∼2%) were associated with age, primarily exhibiting positive relationships, including CpG sites in two genes that have been implicated in previous genome-wide methylation studies of age (FZD9 and MYOD1). In addition, we identified significant associations between three CpG sites (∼1%) and sex, including DNA methylation in CASP6, a gene that may respond to estradiol treatment, and in HSD17B12, which encodes a sex steroid hormone. Although we did not identify any significant associations between DNA methylation and the obesity measures, several nominally significant results were observed in genes related to adipogenesis, obesity, energy homeostasis and glucose homeostasis (ARHGAP9, CDKN2A, FRZB, HOXA5, JAK3, MEST, NPY, PEG3 and SMARCB1). In conclusion, we were able to replicate several findings from previous studies in our Mexican American cohort, supporting an important role for genetic effects on DNA methylation. In addition, we found a significant influence of age and sex on DNA methylation, and report on trend-level, novel associations between DNA methylation and measures of obesity.

Introduction

DNA methylation is an epigenetic modification that effectively regulates gene expression via gene silencing, and is thus considered an important mechanism in growth and development. Within humans, the DNA methyltransferases, DNMT1, DNMT3a and DNMT3b, regulate DNA methylation, which involves the covalent addition of a methyl group at the C(5) position of cytosine in cytosine-guanine (CpG) dinucleotides [1], [2]. Emerging evidence from both animal and human studies suggest that DNA methylation may be heritable [3][7], although the mechanisms of such heritability are speculative.

Both genetic and epigenetic mechanisms can underlie disease development through gene dysregulation, and in many cases act synergistically. DNA methylation changes, at least at the single gene level, have been implicated in most complex diseases, including cancer [2], [8], [9], autoimmune diseases (lupus, multiple sclerosis) [8], [9], psychiatric disorders (bipolar disorder, schizophrenia, autism and depression) [10], obesity [11][14], diabetes [15][19], atherosclerosis [20], hypertension [21] and cardiovascular disease [22][24].

Epigenetic regulation has not been well studied in Mexican Americans, who show a distinct disease profile relative to Caucasian populations. In particular, the prevalence of obesity, diabetes and metabolic syndrome are all elevated among Mexican Americans [25][27], which may be mediated, in part, by epigenetic changes. Moreover, obesity is strongly associated with cardiovascular disease, the leading cause of death worldwide and thus merits particular attention. For our study sample, two measures of obesity were taken: body mass index (BMI), which has been historically used to define obesity (i.e. greater than 30 kg/m2); and waist circumference, which, as an index for central adiposity, has been shown to be a more accurate indicator of obesity, as well as an independent risk factor for diabetes [28]. To date, several animal and human studies have implicated DNA methylation in obesity [11][14], but have primarily focused on global methylation levels or gene-specific DNA methylation rather than a genome-wide approach, or have focused on relatively small sample sizes. In this paper, we report the findings from our pilot study of genome-wide methylation data that investigates: 1) the inherent genetic contribution to epigenetic modifications in a family-based Mexican American cohort; 2) the relationship between DNA methylation and age and sex; and 3) the relationship of DNA methylation with two measures of obesity, BMI and waist circumference.

Materials and Methods

Ethics Statement

Written consent was obtained for all individuals in this study. Approval for this study and for the consent process was received from the Institutional Review Board at The University of Texas Health Science Center at San Antonio.

Population Structure

The San Antonio Family Heart Study (SAFHS) was initiated in 1991 to study the effects of cardiovascular related illness in Mexican Americans [29]. The initial population consisted of 1,431 individuals within 42 extended families. This pilot study assesses DNA methylation profiles of a subset of these individuals (n = 183 successfully profiled subjects from an initial selection of 188 subjects) who represent large founder lineages, chosen to maximize genetic variation. Many of these individuals are related. Pairwise relationships include 196 pairs of 1st degree relatives, 221 pairs of 2nd degree relatives, 401 pairs of 3rd degree relatives, 320 pairs of 4th degree relatives and 39 pairs of 5th degree relatives. The mean age of this population is 42.15 years (standard deviation (sd): 13.74; range 19–75 years) and 55% of individuals are female. Each individual has been assessed for basic anthropometric measures, including those relevant to obesity, waist circumference (n = 182) and BMI (n = 183). In this study, the prevalence of obesity (defined as BMI≥30 kg/m2) is 59%, with 86% of the sample falling in the “overweight” category (BMI≥25 kg/m2). The average BMI for the cohort is 32.40 kg/m2 (sd: 7.40 kg/m2), and is 31.21 kg/m2 (sd: 6.81 kg/m2) for males and 33.38 kg/m2 (sd: 7.75 kg/m2) for females. The average waist circumference for the cohort is 104.86 cm (sd: 16.14 cm), and is 104.08 cm (sd: 14.91 cm) for males and 105.50 cm (sd: 17.12 cm) for females.

Methylation Profiling

Using the EZ-96 DNA Methylation™ Kit (Zymo Research Corp, Irvine, CA) and following the manufacturers’ instructions, we performed bisulfite conversion of 1 µg DNA obtained from peripheral blood of an initially selected 188 Mexican American individuals following an overnight fast. We used Illumina GoldenGate technology incorporating the Methylation Cancer Panel I (Illumina, San Diego, CA) to investigate the quantitative level of methylation at 1,505 CpG sites across 807 genes; 28.6% of genes contained only one CpG site, 57.3% contained two CpG sites and 14.1% of genes had three or more CpG sites. Although the CpG sites within this panel were chosen on the basis of their association or potential association with cancer, they fell within genes involved in DNA repair, cell cycle control, differentiation and apoptosis making them suitable for study in disorders other than cancer. The GoldenGate assay was carried out according to manufacturers’ instructions. Throughout the manuscript, CpG sites are annotated as GENE_Position_Strand, where GENE designates the gene name, Position designates the position of the CpG site relative to the transcription start site either within the promoter (P) or exon (E), and Strand designates whether the CpG site is on the forward (F) or reverse (R) strand.

Initial analysis was carried out using GenomeStudio, applying background normalization by using the signals from built-in controls to minimize the amount of variation in background signals between arrays. DNA methylation for each CpG site is reported as a beta value, which is a score between 0 (completely unmethylated) and 1 (fully methylated). Built-in sample dependent and independent controls were assessed for each sample to look for any outliers. The raw data generated from GenomeStudio is given in Table S1. We excluded three samples with call rates less than 90% (based on detection p-values≥0.05). We excluded two additional samples whose coefficient of variation (based on beta values) was beyond the 97.5th percentile, leaving a total of 183 samples for statistical analyses. Samples used in this study were also genotyped for approximately one million SNPs (see below) as part of a larger study; sex and familial relationships were verified during the genotyping process. Any CpG sites with a call rate less than 0.97 were excluded from analysis (n = 103). The platform used in this study has been shown to discriminate DNA methylation levels (beta values) that differ as little as 17% [30], as such, we excluded probes that showed <0.17 variation in beta values from analysis (n = 355). To account for the non-normal distribution typical of many CpG sites, we used the Lilliefors test for normality [31] at the 0.1 level after FDR correction for multiple testing; non-normally distributed probes were excluded from analysis (n = 641). Any remaining X-linked probes were excluded from analysis (n = 11). Thus, in total, 395 CpG sites were included in subsequent analyses. Polymorphisms within the probe sequence may affect probe binding and therefore bias DNA methylation levels detected. We therefore tested for the presence of 39,706,715 known single nucleotide polymorphisms (SNPs), derived from the 1000 Genomes Project, in each of the 395 probes analyzed.

High Density SNP Genotyping

Each individual assessed in this pilot study has been genotyped for approximately one million SNP markers using several Illumina genotyping arrays, including the Illumina HumanHap550v3, with HumanExon510Sv1, Human1Mv1 and Human1M-Duov3 (Illumina, San Diego, CA). The Infinium Whole-Genome Genotyping Assay was employed according to manufacturers’ instructions. Genotyping of the full SAFHS cohort included >2,000 samples, which underwent stringent quality control measures, excluding any samples with call rates <90% or that showed discrepancies in either sex or known familial relationships. According to these criteria, all 188 samples that underwent methylation profiling were eligible for genotyping analysis. Data cleaning in the full SAFHS population excluded monomorphic SNPs and those with a call rate less than 95%, those that were monomorphic, those whose minor allele was present in <10 individuals and those with Hardy-Weinberg Equilibrium test statistics of p≤10−4 (calculated using SOLAR [32] taking familial relationships into account). A total of 995,320 SNPs were available for analysis of genetic regulation. Allele frequencies were calculated using maximum likelihood estimates in SOLAR [32] and genotypes were checked for Mendelian consistency using Simwalk [33]. Finally, Merlin [34] was used to impute missing genotypes conditional on relatives’ genotypes with a weighted average of possible genotypes being used when an individuals’ genotype could not be inferred with certainty. Raw genotype data for the entire SAFHS cohort will be made available shortly through the Center for Collaborative Genomic Studies on Mental Disorders.

Statistical Analysis

Our general computer package for statistical genetic analysis, SOLAR [32], was used to estimate heritability, identify genetic variation influencing DNA methylation and determine associations with age, sex and obesity measures. Polygenic regression models of the form: were used in the analyses where t is phenotypic trait observation on the ith individual, m is the estimated mean trait value, b is the regression coefficient associated with covariate×(total k covariates), g is the polygenic effect and e is the residual error. The polygenic effect g was modeled by including a pairwise kinship matrix with the following kinship coefficients: identical twins, 1; parent-offspring or sibling, 0.5; and grandparent-grandchild, avuncular, half-siblings or double first cousins, 0.25. Further, the kinship coefficients for 3rd, 4th, 5th and 6th degree relatives are 0.0078, 0.0020, 0.0005 and 0.0001, respectively. The covariates sex, age, age2 and their interactions were used in all relevant analyses. Regression terms were estimated for each covariate, the statistical significance of which was assessed using a likelihood ratio test by comparing the log-likelihoods of models in which the covariate effect is estimated to the likelihoods of models in which the covariate effects were constrained to zero. Given that methylation profiling was performed on individuals that were chosen to maximize capture of founder genomes from the largest pedigrees and were optimally chosen for heritability inference, we calculated power to detect additive genetic heritability. We estimated that we have 50% power to detect heritability as small as 0.23 and 80% power to detect heritability as low as 0.36. To detect associations whilst correcting for multiple tests (Bonferroni correction) our power is 80% to detect variants that account for 16.7% of the variation. As such, we present results based on both Bonferroni correction and the less stringent false discovery rate (FDR) of 5%. The tail area-based estimation of FDR was computed using the R program “fdrtool” [35].

We used the BSMAP software package [36] to examine if the autosomal CpG sites that were significantly associated with sex were cross-reactive with the X-chromosome. For this analysis, we used the Hg19 reference genome with the following in silico conversions – reference, reverse complement, fully methylated and fully unmethylated. We then aligned the three probe sequences against these four reference genomes.

We employed the “measured genotype” model [37], [38], as implemented in the genome-wide association procedures in SOLAR [39], to investigate the effect of genetic variation on methylation levels at each CpG site. Each of the observed CpG sites was tested as a possible covariate influencing waist circumference or BMI using a regression approach that allowed for non-independence. We also tested for effects of proximal genetic variation, with the assumption that significant findings (correlations) would likely indicate genetic effects. Two CpG sites (PECAM1_E32_R and PECAM1_P135_F) could not be annotated and were excluded from this analysis. To test for genetic effects on DNA methylation levels, we took the measured-genotype approach to analyses and defined proximal SNPs as those within 100 kb of the flanking region around the genomic location of the CpG site (n = 29,862). According to these criteria, all 393 CpG sites were assessed for genetic effects. For determining the statistical significance of a correlated SNP, the nominal p-value was estimated by comparing the log-likelihoods of a model with and without the SNP. Since we anticipated that many of the SNPs included in these analyses are likely to be in linkage disequilibrium, we used the method suggested by Li and Ji [40] to correct for multiple testing. In this method, the number of independent tests is determined using the eigenvalues of the SNP×SNP correlation matrix and then using a Sidak correction for multiple testing. We conducted these analyses separately for each chromosome and then summated all of the independent tests (correction for 4,209 tests). Additionally, we adjusted for an FDR of 5%.

Results and Discussion

We have conducted a pilot study to examine the effects of genetic variation on DNA methylation in Mexican Americans, in which these effects, to our knowledge, have not yet been examined. Further, we assessed associations between DNA methylation and age and sex, as well as obesity measures, in order to better understand epigenetic contributions to health.

To ensure DNA methylation levels were not significantly influenced by variation within the probe sequence, we tested for the presence of 39,706,715 SNPs derived from the 1000 Genomes Project. We identified SNPs in 163 of the 395 probes used for analysis (see Table S2); five probes contained three SNPs, 38 probes contained two SNPs and the remaining 120 probes contained a single SNP, three SNPs overlapped with two probes (total of 208 SNPs). Given that most of these SNPs are rare and therefore unlikely to be present in our cohort at a significant frequency, we have highlighted probes containing SNPs with a minor allele frequency (MAF)>5%. We identified SNPs within 24 probes that had a MAF>5%, two of these probes contained two SNPs each. Within these 24 probes, eight contained SNPs within 5 bp of the cytosine molecule of the CpG site and an additional six contained SNPs within 10 bp of the cytosine molecule. Although we included all probes in our analyses, regardless of the presence of SNPs, in the following text (and shown in Table S2), we have highlighted any results that are likely to be influenced by the inclusion of such probes. Only probes that showed significant heritability or genetic effects contained SNPs with a MAF>5%.

DNA Methylation is Heritable

Of the 395 high quality CpG sites examined in this study, 95 (∼24%) exhibited nominally significant (p<0.05) heritability, 28 (∼7%) were heritable at an FDR of 5% and 16 (∼4%) exhibited statistically significant heritability after Bonferroni correction for multiple testing (p<1.27×10−4; Table S3). Table 1 shows all sites that are significantly heritable. Across all loci, the mean heritability was 14%, but it was 36% for the 95 CpG sites that showed nominal evidence for heritability, 52% for the 28 sites that were heritable at an FDR of 5% and was 59% for the 16 CpG sites that were significantly heritable after Bonferroni correction. One of the CpG sites that was found to be heritable (LTA_E28) was within a probe sequence that contained a SNP with a MAF of 35%, 19 bp from the CpG site (see Table S2). Although a SNP at this distance is unlikely to affect probe binding, alternative strategies to determine the methylation status of this CpG site may be warranted. An additional five probes that were associated with heritability contained one or two SNPs, although the MAF of these was ≤1% and they are therefore unlikely to influence results due to the presence of perhaps only one or two copies in our cohort.

thumbnail
Table 1. Heritability and genetic regulation of DNA methylation.

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

Several recent studies have examined the heritability of DNA methylation, finding it to be highly variable for individual CpG sites. Using the Illumina GoldenGate platform Boks, et al, [41] estimated the mean heritability to be 22% (peripheral blood cells) across 430 tested loci and ∼3% of loci were found to be significantly heritable (applying Bonferroni correction). Our results using this same platform revealed a similar, although slightly higher, percentage of significantly heritable sites. Using the Illumina HumanMethylation27 platform that assesses ∼27,000 CpG sites, mean heritability has previously been estimated at: 18% in peripheral blood cells [42]; 12%, 7% and 5% in human umbilical vein epithelial cells, cord blood mononuclear cells and placenta, respectively [43]; and 3% in brain [44]. Our overall mean heritability (14%) is within the range of these study estimates and was determined based on a slightly larger sample size and family data rather than twin or singleton data. Recently, highly variable DNA methylation patterns were identified across three different neonatal tissues and no compelling evidence for a common set of highly heritable DNA methylation variants was seen [43], suggesting that the tissue type analyzed plays a role in heritability estimates.

Given the overlap in the platform and quality control strategies employed, we compared our heritability findings to those of Boks and colleagues, who examined heritability in a twin dataset using 430 CpG probes [41]. Of the 395 probes that were deemed of suitable quality for heritability analysis in our study, only 232 of these (59%) were assessed in the study by Boks, et al, reflecting differential DNA methylation profiles across populations, which may represent ethnic-specific differences or may be due to small cohort sizes and subsequent underpowered analyses. Of the 16 sites we identified as heritable (after Bonferroni correction), five were also identified as heritable in the previous study, one site was nominally significant, seven sites were not assessed and three sites were not heritable in the previous study. In particular, DNA methylation at NPR2_P618, MET_E333_F, AXL_P223, HPN_P823_F and HOXA5_P1324_F were found to be heritable in both studies. Of the CpG sites that were assessed in both studies, evidence for heritability of DNA methylation at IRF5_E101_F, DSC2_E90_F and HS3ST2_E145_R was found in our study only, which may represent population-specific DNA methylation changes driven by heritable factors. We should also note that both our study, and that by Boks and colleagues, assessed heritability in small cohorts (183 and 82 individuals, respectively) and may therefore be underpowered to accurately identify heritable factors influencing DNA methylation or may be prone to the detection of false positive results. A recent study by Breton, et al, [45] presented contrasting findings for the AXL_P223 locus, suggesting that DNA methylation was not heritable but was highly influenced by environmental factors, however this too was performed in a small sample set of only 36 twins and as such is likely limited in its power to detect heritability. Nonetheless, given the overlap for several CpG sites across our study and that by Boks, et al, the case for DNA methylation being a heritable trait is becoming stronger and the application of heritability analysis to larger cohorts will ultimately be highly informative.

Proximal Genetic Variants Influence DNA Methylation

We studied the association of a total of 29,862 SNPs within the 100 Kb flanking region of 393 CpG sites (two sites were excluded due to inability to be annotated), applying a Sidak correction for 4,209 independent tests. In total, 214 statistically significant correlations were detected across 27 CpG sites (p<1.19×10−5, 6.9%; Table S4); the most significant association for each CpG site is shown in Table 1. All 16 CpG sites that were significantly heritable also showed at least nominal evidence of being influenced by genetic effects, although several sites that showed nominal evidence for heritability were not found to be regulated by genetic variation (Table S3). We should note that four CpG sites influenced by genetic effects were within probes that also contained SNPs with a MAF>5% (ABL2_P459, GDF10_P95, HHIP_P307, LTA_E28; see Table S2), three of these had SNPs within 5 bp of the CpG site. It is therefore possible that the genetic effects seen for these CpG sites are in fact due to linkage disequilibrium between genetic variants influencing probe binding rather than a true biological effect. The SNP contained within the LTA_E28 probe (rs2239704) was directly tested for association with DNA methylation levels and found to be significantly associated (p = 1.09×10−16; see Table S4). An additional eight probes contained between one and three SNPs, although all had a MAF≤3% and are therefore unlikely to influence the results.

Recent genome-wide analyses have indicated that DNA methylation may indeed be tightly regulated by genetic factors. In fact, across studies 4–9% of CpG sites show evidence of genetic effects by SNPs within 50 kb-1 Mb of the site [42], [46], [47]. Applying stringent criteria for the localization of proximal genetic variants (100 kb) and correction for multiple testing (Li and Ji’s correction for 4,209 independent CpG site/SNP association tests; adjusted alpha value 1.18×10−5), we found that DNA methylation levels within 6.9% of CpG sites tested were correlated with one or more SNPs. Although this is well in the range of the studies reported above, it is slightly higher than what was reported in the study by Boks, et al, (4%) using the same DNA methylation platform [41]. This may be explained by the smaller sample size (n = 91) they used, a smaller set of SNPs to assess (∼550,000 compared to ∼1,000,000 in our study) and a larger defined region for regulatory variation (1 Mb versus 100 kb in our study), which increases the number of statistical tests that have to be corrected for, thereby reducing power to detect associations. In fact, Quon and colleagues systematically varied the window size to detect regulation using 10 kb, 50 kb, 100 kb, 500 kb and 1 Mb windows surrounding the CpG site to identify the optimum window for detection of heritable sites [44]. Their results showed that a window size of 50 kb yielded the highest number of heritable methylation loci, although similar numbers were seen for the 100 kb window, which we used in our study.

Once again, we compared our regulatory findings to those presented by Boks, et al, who used the same platform but in a singleton dataset, examining 512 probes [41]. Of the 393 high quality probes assessed in our study, only 155 (39%) were also assessed in the study by Boks and colleagues. Of the 23 CpG sites that showed evidence for genetic effects in our study, seven were previously identified as regulated by genetic variation in the study by Boks, et al, (defined as FDR<0.05; ALOX12_E85, ALOX12_P223, GNMT_P197_F, HPN_P823, MET_E333, NPR2_P618_F and SLC22A3_P634_F). In all cases where the same SNP was used to test for regulatory effects, the direction of effect was identical in both studies.

In total, all of our significantly heritable CpG sites (p<1.27×10−4) showed at least nominal evidence for genetic effects; 10 of the 16 heritable CpG sites (63%) showed strong evidence for regulation by genetic variation. For nominally heritable CpG sites (p<0.05, n = 94 assessed for genetic effects), 21% showed significant evidence for genetic effects and a further 65% showed nominal evidence for genetic effects; only 14% of nominally heritable sites showed no evidence of genetic effects (see Table S3). Although it is tempting to suggest that heritability associated with DNA methylation is derived solely through association with regulatory variants, the fact that not all nominally heritable sites are regulated by such variants suggests that at least a portion of the heritability may be due to other factors. Based on in silico SNP analysis and animal studies, Kaminsky and colleagues suggest that the molecular mechanisms driving heritable epigenetic differences may not be limited to sequence differences [6]. Further, evidence from mouse studies suggest that epigenetic marks may not be completely erased, resulting in transgenerational inheritance of epigenetic modifications [3], [4].

Age and Sex Effects on DNA Methylation Profiles

We also tested for effects of age and sex on DNA methylation profiles and found substantial evidence for both; Table 2 displays CpG sites significantly associated with age and sex after Bonferroni correction. DNA methylation at 88 of the 395 CpG sites (∼22%) showed nominally significant associations with age (p<0.05), 24 CpG sites (∼6%) were significant at an FDR of 5% and six CpG sites (∼2%) were significant after Bonferroni correction for multiple testing (Table S5). Four of the probes analyzed contained one or two SNPs, however these were all at a MAF≤1% and would be unlikely to influence the results (see Table S2). Only one (CARD15_P302) of the six associations that we identified was also identified as significant by Boks and colleagues [41], with one site nominally significant (p<0.05), and the other four not assessed. For sex, 50 CpG sites (∼13%) showed nominally significant associations (p<0.05) and three (∼1%) of these were significantly associated with sex after Bonferroni correction for multiple testing and at an FDR of 5% (Table S6). One of the significantly associated probes contained a SNP with a MAF of 1% (see Table S2), which is unlikely to influence the results. Similar to our results with age, the study by Boks, et al, [41] did not assess any of the three CpG sites that we identified as significantly associated with sex.

thumbnail
Table 2. Significant age and sex associations with DNA methylation.

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

Recent genome-wide studies have identified a number of age related differentially methylated regions with 78–98% of associated CpG sites showing positive associations [42], [48], [49]. Similarly, of the six CpG sites that we found to be associated with age, five showed a positive association (83%); only the CARD15_P302 site showed a negative direction of effect (i.e. decreased methylation was associated with increased age), and this was the same direction of effect that was observed in the study by Boks, et al [41]. Positive associations were seen for two CpG sites within genes that have previously been implicated with age in one or more genome-wide methylation studies (FZD9 and MYOD1) [42], [49], [50]. Of these, methylation within MYOD1 was associated with age in all three previously published studies, as well as in our own. Although two studies were performed in peripheral blood cells, the third assessed DNA methylation that was consistently associated with aging across four brain regions, suggesting that the DNA methylation changes associated with the MYOD1 gene may be general rather than tissue-specific. In addition to replicating results obtained from these genome-wide studies, we have identified new loci, possibly specific to the Mexican American population, whose DNA methylation status may be related to age. Due to the small number of samples assessed in our study however, such findings would need to be verified.

Despite evidence for sexual dimorphism in a number of diseases, a recent genome-wide association study meta-analysis [51] failed to reveal any genome-wide significant common SNP differences between men and women, which suggests that at least some sex-based differences in disease may be driven more so by either rare variants or epigenetic effects. It is clear that DNA methylation plays a sex-specific role during fetal germline progression [52] and a recent genome-wide methylation study identified DNA methylation differences between the two sexes, although these differences were subtle [53]. Similarly, we identified significant associations between sex and DNA methylation at three CpG sites, each of which has been at least loosely implicated in sex-specific functions of disease. We found highly significant DNA methylation changes at CASP6_P201, showing increased methylation in women. Although a sex-specific role has not yet been defined, an early study in rats suggests that caspase-6 is involved in the control of apoptotic processes at estrus following treatment with the sex-hormone 17β-estradiol [54]. Increased methylation at HSD17B12_P97, which encodes a sex steroid hormone, was also seen in females. HSD17B12 has been shown to be involved in the catalytic conversion of oestrone to estradiol and has also been implicated in fatty acid synthesis (reviewed in [55]), which shows sex-specific differences (reviewed in [56]). DNA methylation at RET_seq_54_S260 was increased in males, and although specific sex effects have not been noted, genetic variants within this gene have been shown to have different effects in males and females, and parent of origin effects involving such variants have been noted in Hirschsprung disease [57], [58]. It is important to note that some level of cross-reactivity with sex chromosomes has been reported in ∼6–10% of the autosomal probes used in some of the higher-throughput Illumina platforms (i.e. 27 k and 450 k methylation chips), due to high sequence homology [59], [60]. This may generate spurious signals that overestimate the effect of DNA methylation on sex and as such, DNA methylation changes that are determined by probe-based methods and are associated with sex should be interpreted with caution. Using BSMAP [36], we observed that none of the three significant sex-associated CpG sites identified in our study (CASP6_P201, HSD17B12_P97 or RET_seq_54_S260) could be mapped to the X-chromosome, using the reference, reverse complement, bisulfite converted fully methylated and bisulfite converted fully unmethylated genomes. Further, given some of the plausible biological roles discussed above, we believe that this study provides additional evidence for a potential role of DNA methylation in sex.

Methylation Profiles show Nominal Association with Obesity Measures

To investigate a potential role of DNA methylation in disease, we looked for association with the obesity measures, waist circumference and BMI. In particular, waist circumference, which measures central obesity, is considered a stronger indicator of risk for diabetes and cardiovascular disease than BMI [28]. Table 3 shows the ten most significant associations with waist circumference, although none of these are within the threshold for multiple testing correction (p<1.27×10−4). In total, we identified 27 and 12 nominally significant (p<0.05) associations between DNA methylation and waist circumference and BMI, respectively (Table S7). Several nominally associated CpG sites are in genes that have a known or plausible role in obesity and diabetes.

thumbnail
Table 3. Most highly significant associations between DNA methylation and obesity measures.

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

Methylation levels of a CpG site within MEST, were found to be negatively associated with both waist circumference (nominal p = 1.29×10−2) and BMI (nominal p = 4.41×10−2). MEST is well known to be involved in obesity and glucose metabolism. Gene expression increases have been shown in white adipose tissue of mice with diet-induced and genetically caused obesity/diabetes [61] and in humans, decreased methylation of MEST has been shown in morbidly obese individuals as well as in newborns of mothers who had gestational diabetes [62]. We also identified a negative association between DNA methylation within NPY and waist circumference (nominal p = 2.96×10−2), but not with BMI. NPY is a hypothalamic peptide which regulates feeding behavior and energy homeostasis and has been proposed as a target for anti-obesity drugs [63]. Studies have shown that knockdown of NPY expression in the hypothalamus increases energy expenditure [64] and overexpression within the paraventricular nucleus results in obesity via increased food intake [65]. We identified a positive association between DNA methylation in HOXA5 and BMI (nominal p = 2.64×10−2), but not waist circumference. HOXA5 has been implicated in body fat distribution, where it is more highly expressed in visceral (i.e. abdominal) adipose tissue as opposed to subcutaneous (i.e. gluteal) adipose tissue [66], [67]. Further, HOXA5 is up-regulated in abdominal subcutaneous adipose tissue following fat loss (one year after bariatric surgery) [68]. In addition, nominal associations were seen between either waist circumference or BMI and DNA methylation within several other genes previously implicated in serum triglyceride levels (ARHGAP9 [69]); obesity and adipogenesis (CDKN2A (p16) [70], [71], FRZB [72], JAK3 [73], PEG3 [74]); and glucose homeostasis and insulin resistance (CDKN2A [75][77], SMARCB1 [78]).

Study Limitations

It is important to note the limitations of this study. As both the number of samples (n = 183) and number of CpG sites assessed (n = 395) are small, the power and scope of this study is limited. We have only assayed DNA methylation in one tissue type, peripheral blood, which may limit any generalized conclusions and may not be representative of changes that would be seen in diseased tissue (e.g. adipose tissue). However, given that lipid profiling, glucose measures and cardiovascular-related biomarkers are measured in blood, and that each of these are impacted by obesity, blood represents a good surrogate tissue for these analyses. Many of the results reported here further validate previously reported data and may also offer insight into any ethnic-specific differences.

Conclusions

We examined the effects of genetic variation, sex and age on DNA methylation and assessed the role of DNA methylation in measures of obesity. Other studies have previously investigated the effects of genetic variation, age and sex using an identical or higher-throughput Illumina platform. However, this is the first study that has examined such factors in a Mexican American population. We have demonstrated here, further evidence supporting a role for genetic effects on DNA methylation, by replicating previous findings and presenting new data. Although several of our findings overlap with previous studies, the identification of novel CpG loci that are heritable, regulated by genetic variation, or that are influenced by age and sex may suggest the existence of ethnic-specific factors that drive regulation of DNA methylation. Such novel findings could also suggest that currently published studies, as well as our own, may be underpowered in their ability to detect true associations or have generated false positive results due to small sample sizes. Differences that are seen across studies could also be due to environmental effects, some of which may themselves be ethnic- or population-specific (e.g. different diets, residential areas). We also provide nominal evidence that DNA methylation may be implicated in obesity in Mexican Americans, however additional follow up studies, particularly in a larger population are required to fully assess such implications.

Supporting Information

Table S1.

Raw data for 188 samples and 1,505 CpG sites.

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

(ZIP)

Table S2.

List of all SNPs within probe sequences.

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

(XLSX)

Table S3.

Heritability estimates for 395 CpG sites assessed.

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

(XLSX)

Table S4.

Complete association analysis between all proximal genetic variants (defined as those within 100 kb) and 393 CpG sites.

https://doi.org/10.1371/journal.pone.0073950.s004

(XLSX)

Table S5.

Complete association analysis between age and DNA methylation at 395 CpG sites.

https://doi.org/10.1371/journal.pone.0073950.s005

(XLSX)

Table S6.

Complete association analysis between sex and DNA methylation at 395 CpG sites.

https://doi.org/10.1371/journal.pone.0073950.s006

(XLSX)

Table S7.

Complete association analysis between obesity measures (waist circumference and BMI) and DNA methylation at 395 CpG sites.

https://doi.org/10.1371/journal.pone.0073950.s007

(XLSX)

Author Contributions

Conceived and designed the experiments: MAC JB. Performed the experiments: MAC. Analyzed the data: MAC HK MZK JC JMP HHHG LA TDD MCM JB. Contributed reagents/materials/analysis tools: JEC AGC. Wrote the paper: MAC HK MZK.

References

  1. 1. Moss TJ, Wallrath LL (2007) Connections between epigenetic gene silencing and human disease. Mutat Res 618: 163–174.
  2. 2. Park SY, Kim BH, Kim JH, Cho NY, Choi M, et al. (2007) Methylation profiles of CpG island loci in major types of human cancers. J Korean Med Sci 22: 311–317.
  3. 3. Roemer I, Reik W, Dean W, Klose J (1997) Epigenetic inheritance in the mouse. Curr Biol 7: 277–280.
  4. 4. Morgan HD, Sutherland HG, Martin DI, Whitelaw E (1999) Epigenetic inheritance at the agouti locus in the mouse. Nat Genet 23: 314–318.
  5. 5. 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.
  6. 6. 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.
  7. 7. Heijmans BT, Kremer D, Tobi EW, Boomsma DI, Slagboom PE (2007) Heritable rather than age-related environmental and stochastic factors dominate variation in DNA methylation of the human IGF2/H19 locus. Hum Mol Genet 16: 547–554.
  8. 8. Millington GW (2008) Epigenetics and dermatological disease. Pharmacogenomics 9: 1835–1850.
  9. 9. Wilson AS, Power BE, Molloy PL (2007) DNA hypomethylation and human diseases. Biochim Biophys Acta 1775: 138–162.
  10. 10. van Vliet J, Oates NA, Whitelaw E (2007) Epigenetic mechanisms in the context of complex diseases. Cell Mol Life Sci 64: 1531–1538.
  11. 11. Waterland RA (2009) Is epigenetics an important link between early life events and adult disease? Horm Res 71 Suppl 113–16.
  12. 12. Milagro FI, Gomez-Abellan P, Campion J, Martinez JA, Ordovas JM, et al. (2012) CLOCK, PER2 and BMAL1 DNA methylation: association with obesity and metabolic syndrome characteristics and monounsaturated fat intake. Chronobiol Int 29: 1180–1194.
  13. 13. Fujiki K, Kano F, Shiota K, Murata M (2009) Expression of the peroxisome proliferator activated receptor gamma gene is repressed by DNA methylation in visceral adipose tissue of mouse models of diabetes. BMC Biol 7: 38.
  14. 14. Wang X, Zhu H, Snieder H, Su S, Munn D, et al. (2010) Obesity related methylation changes in DNA of peripheral blood leukocytes. BMC Med 8: 87.
  15. 15. Shao WJ, Tao LY, Gao C, Xie JY, Zhao RQ (2008) Alterations in methylation and expression levels of imprinted genes H19 and Igf2 in the fetuses of diabetic mice. Comp Med 58: 341–346.
  16. 16. Jiang MH, Fei J, Lan MS, Lu ZP, Liu M, et al. (2008) Hypermethylation of hepatic Gck promoter in ageing rats contributes to diabetogenic potential. Diabetologia 51: 1525–1533.
  17. 17. Park JH, Stoffers DA, Nicholls RD, Simmons RA (2008) Development of type 2 diabetes following intrauterine growth retardation in rats is associated with progressive epigenetic silencing of Pdx1. J Clin Invest 118: 2316–2324.
  18. 18. Mackay DJ, Boonen SE, Clayton-Smith J, Goodship J, Hahnemann JM, et al. (2006) A maternal hypomethylation syndrome presenting as transient neonatal diabetes mellitus. Hum Genet 120: 262–269.
  19. 19. Ling C, Del Guerra S, Lupi R, Ronn 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.
  20. 20. Kim J, Kim JY, Song KS, Lee YH, Seo JS, et al. (2007) Epigenetic changes in estrogen receptor beta gene in atherosclerotic cardiovascular tissues and in-vitro vascular senescence. Biochim Biophys Acta 1772: 72–80.
  21. 21. Smolarek I, Wyszko E, Barciszewska AM, Nowak S, Gawronska I, et al. (2010) Global DNA methylation changes in blood of patients with essential hypertension. Med Sci Monit 16: CR149–155.
  22. 22. Stenvinkel P, Karimi M, Johansson S, Axelsson J, Suliman M, et al. (2007) Impact of inflammation on epigenetic DNA methylation - a novel risk factor for cardiovascular disease? J Intern Med 261: 488–499.
  23. 23. Baccarelli A, Wright R, Bollati V, Litonjua A, Zanobetti A, et al. (2010) Ischemic heart disease and stroke in relation to blood DNA methylation. Epidemiology 21: 819–828.
  24. 24. Movassagh M, Choy MK, Goddard M, Bennett MR, Down TA, et al. (2010) Differential DNA methylation correlates with differential expression of angiogenic factors in human heart failure. PLoS One 5: e8564.
  25. 25. Flegal KM, Carroll MD, Ogden CL, Curtin LR (2010) Prevalence and trends in obesity among US adults, 1999–2008. Jama 303: 235–241.
  26. 26. Romero CX, Romero TE, Shlay JC, Ogden LG, Dabelea D (2012) Changing trends in the prevalence and disparities of obesity and other cardiovascular disease risk factors in three racial/ethnic groups of USA adults. Adv Prev Med 2012: 172423.
  27. 27. Ford ES (2005) Prevalence of the metabolic syndrome defined by the International Diabetes Federation among adults in the U.S. Diabetes Care. 28: 2745–2749.
  28. 28. Janiszewski PM, Janssen I, Ross R (2007) Does waist circumference predict diabetes and cardiovascular disease beyond commonly evaluated cardiometabolic risk factors? Diabetes Care 30: 3105–3109.
  29. 29. Mitchell BD, Kammerer CM, Blangero J, Mahaney MC, Rainwater DL, et al. (1996) Genetic and environmental contributions to cardiovascular risk factors in Mexican Americans. The San Antonio Family Heart Study. Circulation 94: 2159–2170.
  30. 30. Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, et al. (2006) High-throughput DNA methylation profiling using universal bead arrays. Genome Res 16: 383–393.
  31. 31. Lilliefors HW (1967) On the Kolmogorov-Smirnov test for normality with mean and variance unknown. Journal of the American Statistical Association 62: 399–402.
  32. 32. Almasy L, Blangero J (1998) Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet 62: 1198–1211.
  33. 33. Sobel E, Lange K (1996) Descent graphs in pedigree analysis: applications to haplotyping, location scores, and marker-sharing statistics. Am J Hum Genet 58: 1323–1337.
  34. 34. Abecasis GR, Cherny SS, Cookson WO, Cardon LR (2002) Merlin–rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 30: 97–101.
  35. 35. Strimmer K (2008) fdrtool: a versatile R package for estimating local and tail area-based false discovery rates. Bioinformatics 24: 1461–1462.
  36. 36. Xi Y, Li W (2009) BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinformatics 10: 232.
  37. 37. Hopper JL, Mathews JD (1982) Extensions to multivariate normal models for pedigree analysis. Ann Hum Genet 46: 373–383.
  38. 38. Boerwinkle E, Chakraborty R, Sing CF (1986) The use of measured genotype information in the analysis of quantitative phenotypes in man. I. Models and analytical methods. Ann Hum Genet 50: 181–194.
  39. 39. Blangero J, Goring HH, Kent JW Jr, Williams JT, Peterson CP, et al. (2005) Quantitative trait nucleotide analysis using Bayesian model selection. Hum Biol 77: 541–559.
  40. 40. Li J, Ji L (2005) Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity 95: 221–227.
  41. 41. Boks MP, Derks EM, Weisenberger DJ, Strengman E, Janson E, et al. (2009) The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PLoS ONE 4: e6767.
  42. 42. Bell JT, Tsai PC, Yang TP, Pidsley R, Nisbet J, et al. (2012) Epigenome-wide scans identify differentially methylated regions for age and age-related phenotypes in a healthy ageing population. PLoS Genet 8: e1002629.
  43. 43. Gordon L, Joo JE, Powell JE, Ollikainen M, Novakovic B, et al. (2012) Neonatal DNA methylation profile in human twins is specified by a complex interplay between intrauterine environmental and genetic factors, subject to tissue-specific influence. Genome Res 22: 1395–1406.
  44. 44. Quon G, Lippert C, Heckerman D, Listgarten J (2013) Patterns of methylation heritability in a genome-wide analysis of four brain regions. Nucleic Acids Res 41: 2095–2104.
  45. 45. Breton CV, Salam MT, Gilliland FD (2011) Heritability and role for the environment in DNA methylation in AXL receptor tyrosine kinase. Epigenetics 6: 895–898.
  46. 46. Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, et al. (2010) Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet 6: e1000952.
  47. 47. Zhang D, Cheng L, Badner JA, Chen C, Chen Q, et al. (2010) Genetic control of individual differences in gene-specific methylation in human brain. Am J Hum Genet 86: 411–419.
  48. 48. Bocklandt S, Lin W, Sehl ME, Sanchez FJ, Sinsheimer JS, et al. (2011) Epigenetic predictor of age. PLoS ONE 6: e14821.
  49. 49. Hernandez DG, Nalls MA, Gibbs JR, Arepalli S, van der Brug M, et al. (2011) Distinct DNA methylation changes highly correlated with chronological age in the human brain. Hum Mol Genet 20: 1164–1172.
  50. 50. Rakyan VK, Down TA, Maslau S, Andrew T, Yang TP, et al. (2010) Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res 20: 434–439.
  51. 51. Boraska V, Jeroncic A, Colonna V, Southam L, Nyholt DR, et al. (2012) Genome-wide meta-analysis of common variant differences between men and women. Hum Mol Genet 21: 4805–4815.
  52. 52. Kobayashi H, Sakurai T, Miura F, Imai M, Mochiduki K, et al.. (2013) High-resolution DNA methylome analysis of primordial germ cells identifies gender-specific reprogramming in mice. Genome Res.
  53. 53. Lam LL, Emberly E, Fraser HB, Neumann SM, Chen E, et al. (2012) Factors underlying variable DNA methylation in a human community cohort. Proc Natl Acad Sci U S A 109 Suppl 217253–17260.
  54. 54. Leblanc V, Dery MC, Shooner C, Asselin E (2003) Opposite regulation of XIAP and Smac/DIABLO in the rat endometrium in response to 17beta-estradiol at estrus. Reprod Biol Endocrinol 1: 59.
  55. 55. Saloniemi T, Jokela H, Strauss L, Pakarinen P, Poutanen M (2012) The diversity of sex steroid action: novel functions of hydroxysteroid (17beta) dehydrogenases as revealed by genetically modified mouse models. J Endocrinol 212: 27–40.
  56. 56. Lohner S, Fekete K, Marosvolgyi T, Decsi T (2013) Gender Differences in the Long-Chain Polyunsaturated Fatty Acid Status: Systematic Review of 51 Publications. Ann Nutr Metab 62: 98–112.
  57. 57. Emison ES, McCallion AS, Kashuk CS, Bush RT, Grice E, et al. (2005) A common sex-dependent mutation in a RET enhancer underlies Hirschsprung disease risk. Nature 434: 857–863.
  58. 58. Jannot AS, Amiel J, Pelet A, Lantieri F, Fernandez RM, et al. (2012) Male and female differential reproductive rate could explain parental transmission asymmetry of mutation origin in Hirschsprung disease. Eur J Hum Genet 20: 917–920.
  59. 59. Chen YA, Choufani S, Ferreira JC, Grafodatskaya D, Butcher DT, et al. (2011) Sequence overlap between autosomal and sex-linked probes on the Illumina HumanMethylation27 microarray. Genomics 97: 214–222.
  60. 60. Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, et al.. (2013) Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics 8.
  61. 61. Takahashi M, Kamei Y, Ezaki O (2005) Mest/Peg1 imprinted gene enlarges adipocytes and is a marker of adipocyte size. Am J Physiol Endocrinol Metab 288: E117–124.
  62. 62. El Hajj N, Pliushch G, Schneider E, Dittrich M, Muller T, et al.. (2012) Metabolic Programming of MEST DNA Methylation by Intrauterine Exposure to Gestational Diabetes Mellitus. Diabetes.
  63. 63. Yulyaningsih E, Zhang L, Herzog H, Sainsbury A (2011) NPY receptors as potential targets for anti-obesity drug development. Br J Pharmacol 163: 1170–1202.
  64. 64. Chao PT, Yang L, Aja S, Moran TH, Bi S (2011) Knockdown of NPY expression in the dorsomedial hypothalamus promotes development of brown adipocytes and prevents diet-induced obesity. Cell Metab 13: 573–583.
  65. 65. Tiesjema B, la Fleur SE, Luijendijk MC, Adan RA (2009) Sustained NPY overexpression in the PVN results in obesity via temporarily increasing food intake. Obesity (Silver Spring) 17: 1448–1450.
  66. 66. Karastergiou K, Fried SK, Xie H, Lee MJ, Divoux A, et al. (2013) Distinct developmental signatures of human abdominal and gluteal subcutaneous adipose tissue depots. J Clin Endocrinol Metab 98: 362–371.
  67. 67. Gesta S, Bluher M, Yamamoto Y, Norris AW, Berndt J, et al. (2006) Evidence for a role of developmental genes in the origin of obesity and body fat distribution. Proc Natl Acad Sci U S A 103: 6676–6681.
  68. 68. Dankel SN, Fadnes DJ, Stavrum AK, Stansberg C, Holdhus R, et al. (2010) Switch from stress response to homeobox transcription factors in adipose tissue after profound fat loss. PLoS ONE 5: e11033.
  69. 69. Haas BE, Horvath S, Pietilainen KH, Cantor RM, Nikkola E, et al. (2012) Adipose co-expression networks across Finns and Mexicans identify novel triglyceride-associated genes. BMC Med Genomics 5: 61.
  70. 70. Barfield RT, Kilaru V, Smith AK, Conneely KN (2012) CpGassoc: an R function for analysis of DNA methylation microarray data. Bioinformatics 28: 1280–1281.
  71. 71. minfi: Analyze Illumina’s 450 k methylation arrays. Available: http://www.bioconductor.org/packages/release/bioc/html/minfi.html. Accessed 2013 Aug 5.
  72. 72. Hu E, Zhu Y, Fredrickson T, Barnes M, Kelsell D, et al. (1998) Tissue restricted expression of two human Frzbs in preadipocytes and pancreas. Biochem Biophys Res Commun 247: 287–293.
  73. 73. Song HY, Kim MR, Lee MJ, Jeon ES, Bae YC, et al. (2007) Oncostatin M decreases adiponectin expression and induces dedifferentiation of adipocytes by JAK3- and MEK-dependent pathways. Int J Biochem Cell Biol 39: 439–449.
  74. 74. Curley JP, Pinnock SB, Dickson SL, Thresher R, Miyoshi N, et al. (2005) Increased body fat in mice with a targeted mutation of the paternally expressed imprinted gene Peg3. Faseb J 19: 1302–1304.
  75. 75. Davis S, Du P, Bilke S, Triche T (2010) Methylumi: handle illumina methylation data (version 2.0.1). Available: www.bioconductor.org/packages/release/bioc/html/methylumi.html. Accessed 2013 Aug 5.
  76. 76. Du P, Kibbe WA, Lin SM (2008) lumi: a pipeline for processing Illumina microarray. Bioinformatics 24: 1547–1548.
  77. 77. Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr AR, James KD, et al.. (2010) Orphan CpG islands identify numerous conserved promoters in the mammalian genome. PLoS Genet 6.
  78. 78. Gresh L, Bourachot B, Reimann A, Guigas B, Fiette L, et al. (2005) The SWI/SNF chromatin-remodeling complex subunit SNF5 is essential for hepatocyte differentiation. Embo J 24: 3313–3324.