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

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.     

Identification of a DNA Methylome Profile of Esophageal Squamous Cell Carcinoma and Potential Plasma Epigenetic Biomarkers for Early Diagnosis

DNA methylation is a critical epigenetic mechanism involved in key cellular processes. Its deregulation has been linked to many human cancers including esophageal squamous cell carcinoma (ESCC). This study was designed to explore the whole methylation status of ESCC and to identify potential plasma biomarkers for early diagnosis. We used Infinium Methylation 450k array to analyze ESCC tissues (n = 4), paired normal surrounding tissues (n = 4) and normal mucosa from healthy individuals (n = 4), and combined these with gene expression data from the GEO database. One hundred and sixty eight genes had differentially methylated CpG sites in their promoter region and a gene expression pattern inverse to the direction of change in DNA methylation. These genes were involved in several cancer-related pathways. Three genes were validated in additional 42 ESCC tissues and paired normal surrounding tissues. The methylation frequency of EPB41L3, GPX3, and COL14A1 were higher in tumor tissues than in normal surrounding tissues (P<0.017). The higher methylation frequency of EPB41l3 was correlated with large tumor size (P = 0.044) and advanced pT tumor stage (P = 0.001). The higher methylation frequency of GPX3 and COL14A1 were correlated with advanced pN tumor stage (P = 0.001 and P<0.001). The methylation of EPB41L3, GPX3, and COL14A1 genes were only found in ESCC patients'' plasma, but not in normal individuals upon testing 42 ESCC patients and 50 healthy individuals. Diagnostic sensitivity was increased when methylation of any of the 3 genes were counted (64.3% sensitivity and 100% specificity). These differentially methylated genes in plasma may be used as biomarkers for early diagnosis of ESCC.     

Identification of New Differentially Methylated Genes That Have Potential Functional Consequences in Prostate Cancer

Many differentially methylated genes have been identified in prostate cancer (PCa), primarily using candidate gene-based assays. Recently, several global DNA methylation profiles have been reported in PCa, however, each of these has weaknesses in terms of ability to observe global DNA methylation alterations in PCa. We hypothesize that there remains unidentified aberrant DNA methylation in PCa, which may be identified using higher resolution assay methods. We used the newly developed Illumina HumanMethylation450 BeadChip in PCa (n = 19) and adjacent normal tissues (n = 4) and combined these with gene expression data for identifying new DNA methylation that may have functional consequences in PCa development and progression. We also confirmed our methylation results in an independent data set. Two aberrant DNA methylation genes were validated among an additional 56 PCa samples and 55 adjacent normal tissues. A total 28,735 CpG sites showed significant differences in DNA methylation (FDR adjusted P<0.05), defined as a mean methylation difference of at least 20% between PCa and normal samples. Furthermore, a total of 122 genes had more than one differentially methylated CpG site in their promoter region and a gene expression pattern that was inverse to the direction of change in DNA methylation (e.g. decreased expression with increased methylation, and vice-versa). Aberrant DNA methylation of two genes, AOX1 and SPON2, were confirmed via bisulfate sequencing, with most of the respective CpG sites showing significant differences between tumor samples and normal tissues. The AOX1 promoter region showed hypermethylation in 92.6% of 54 tested PCa samples in contrast to only three out of 53 tested normal tissues. This study used a new BeadChip combined with gene expression data in PCa to identify novel differentially methylated CpG sites located within genes. The newly identified differentially methylated genes may be used as biomarkers for PCa diagnosis.     

Heterogeneity in White Blood Cells Has Potential to Confound DNA Methylation Measurements

Epigenetic studies are commonly conducted on DNA from tissue samples. However, tissues are ensembles of cells that may each have their own epigenetic profile, and therefore inter-individual cellular heterogeneity may compromise these studies. Here, we explore the potential for such confounding on DNA methylation measurement outcomes when using DNA from whole blood. DNA methylation was measured using pyrosequencing-based methodology in whole blood (n = 50–179) and in two white blood cell fractions (n = 20), isolated using density gradient centrifugation, in four CGIs (CpG Islands) located in genes HHEX (10 CpG sites assayed), KCNJ11 (8 CpGs), KCNQ1 (4 CpGs) and PM20D1 (7 CpGs). Cellular heterogeneity (variation in proportional white blood cell counts of neutrophils, lymphocytes, monocytes, eosinophils and basophils, counted by an automated cell counter) explained up to 40% (p<0.0001) of the inter-individual variation in whole blood DNA methylation levels in the HHEX CGI, but not a significant proportion of the variation in the other three CGIs tested. DNA methylation levels in the two cell fractions, polymorphonuclear and mononuclear cells, differed significantly in the HHEX CGI; specifically the average absolute difference ranged between 3.4–15.7 percentage points per CpG site. In the other three CGIs tested, methylation levels in the two fractions did not differ significantly, and/or the difference was more moderate. In the examined CGIs, methylation levels were highly correlated between cell fractions. In summary, our analysis detects region-specific differential DNA methylation between white blood cell subtypes, which can confound the outcome of whole blood DNA methylation measurements. Finally, by demonstrating the high correlation between methylation levels in cell fractions, our results suggest a possibility to use a proportional number of a single white blood cell type to correct for this confounding effect in analyses.     

DNA Methylation Profiling across the Spectrum of HPV-Associated Anal Squamous Neoplasia

Background

Changes in host tumor genome DNA methylation patterns are among the molecular alterations associated with HPV-related carcinogenesis. However, there is little known about the epigenetic changes associated specifically with the development of anal squamous cell cancer (SCC). We sought to characterize broad methylation profiles across the spectrum of anal squamous neoplasia.

Methodology/Principal Findings

Twenty-nine formalin-fixed paraffin embedded samples from 24 patients were evaluated and included adjacent histologically normal anal mucosa (NM; n = 3), SCC-in situ (SCC-IS; n = 11) and invasive SCC (n = 15). Thirteen women and 11 men with a median age of 44 years (range 26–81) were included in the study. Using the SFP₁₀ LiPA HPV-typing system, HPV was detected in at least one tissue from all patients with 93% (27/29) being positive for high-risk HPV types and 14 (93%) of 15 invasive SCC tissues testing positive for HPV 16. Bisulfite-modified DNA was interrogated for methylation at 1,505 CpG loci representing 807 genes using the Illumina GoldenGate Methylation Array. When comparing the progression from normal anal mucosa and SCC-IS to invasive SCC, 22 CpG loci representing 20 genes demonstrated significant differential methylation (p<0.01). The majority of differentially methylated gene targets occurred at or close to specific chromosomal locations such as previously described HPV methylation “hotspots” and viral integration sites.

Conclusions

We have identified a panel of differentially methlylated CpG loci across the spectrum of HPV-associated squamous neoplasia of the anus. To our knowledge, this is the first reported application of large-scale high throughput methylation analysis for the study of anal neoplasia. Our findings support further investigations into the role of host-genome methylation in HPV-associated anal carcinogenesis with implications towards enhanced diagnosis and screening strategies.     

The Relationship of DNA Methylation with Age,Gender and Genotype in Twins and Healthy Controls

Cytosine-5 methylation within CpG dinucleotides is a potentially important mechanism of epigenetic influence on human traits and disease. In addition to influences of age and gender, genetic control of DNA methylation levels has recently been described. We used whole blood genomic DNA in a twin set (23 MZ twin-pairs and 23 DZ twin-pairs, N = 92) as well as healthy controls (N = 96) to investigate heritability and relationship with age and gender of selected DNA methylation profiles using readily commercially available GoldenGate bead array technology. Despite the inability to detect meaningful methylation differences in the majority of CpG loci due to tissue type and locus selection issues, we found replicable significant associations of DNA methylation with age and gender. We identified associations of genetically heritable single nucleotide polymorphisms with large differences in DNA methylation levels near the polymorphism (cis effects) as well as associations with much smaller differences in DNA methylation levels elsewhere in the human genome (trans effects). Our results demonstrate the feasibility of array-based approaches in studies of DNA methylation and highlight the vast differences between individual loci. The identification of CpG loci of which DNA methylation levels are under genetic control or are related to age or gender will facilitate further studies into the role of DNA methylation and disease.     

DNA Isolation Method Is a Source of Global DNA Methylation Variability Measured with LUMA. Experimental Analysis and a Systematic Review

In DNA methylation, methyl groups are covalently bound to CpG dinucleotides. However, the assumption that methyl groups are not lost during routine DNA extraction has not been empirically tested. To avoid nonbiological associations in DNA methylation studies, it is essential to account for potential batch effect bias in the assessment of this epigenetic mechanism. Our purpose was to determine if the DNA isolation method is an independent source of variability in methylation status. We quantified Global DNA Methylation (GDM) by luminometric methylation assay (LUMA), comparing the results from 3 different DNA isolation methods. In the controlled analysis (n = 9), GDM differed slightly for the same individual depending on extraction method. In the population analysis (n = 580) there were significant differences in GDM between the 3 DNA isolation methods (medians, 78.1%, 76.5% and 75.1%; p<0.001). A systematic review of published data from LUMA GDM studies that specify DNA extraction methods is concordant with our findings. DNA isolation method is a source of GDM variability measured with LUMA. To avoid possible bias, the method used should be reported and taken into account in future DNA methylation studies.     

Identification of Methylated Genes Associated with Aggressive Bladder Cancer

Approximately 500,000 individuals diagnosed with bladder cancer in the U.S. require routine cystoscopic follow-up to monitor for disease recurrences or progression, resulting in over $2 billion in annual expenditures. Identification of new diagnostic and monitoring strategies are clearly needed, and markers related to DNA methylation alterations hold great promise due to their stability, objective measurement, and known associations with the disease and with its clinical features. To identify novel epigenetic markers of aggressive bladder cancer, we utilized a high-throughput DNA methylation bead-array in two distinct population-based series of incident bladder cancer (n = 73 and n = 264, respectively). We then validated the association between methylation of these candidate loci with tumor grade in a third population (n = 245) through bisulfite pyrosequencing of candidate loci. Array based analyses identified 5 loci for further confirmation with bisulfite pyrosequencing. We identified and confirmed that increased promoter methylation of HOXB2 is significantly and independently associated with invasive bladder cancer and methylation of HOXB2, KRT13 and FRZB together significantly predict high-grade non-invasive disease. Methylation of these genes may be useful as clinical markers of the disease and may point to genes and pathways worthy of additional examination as novel targets for therapeutic treatment.     

Hypomethylation of Serum Blood Clot DNA,but Not Plasma EDTA-Blood Cell Pellet DNA,from Vitamin B12-Deficient Subjects

Vitamin B12, a co-factor in methyl-group transfer, is important in maintaining DNA (deoxycytidine) methylation. Using two independent assays we examined the effect of vitamin B12-deficiency (plasma vitamin B12<148 pmol/L) on DNA methylation in women of childbearing age. Coagulated blood clot DNA from vitamin B12-deficient women had significantly (p<0.001) lower percentage deoxycytidine methylation (3.23±0.66%; n = 248) and greater [3 H]methyl-acceptance (42,859±9,699 cpm; n = 17) than DNA from B12-replete women (4.44±0.18%; n = 128 and 26,049±2,814 cpm; n = 11) [correlation between assays: r = –0.8538; p<0.001; n = 28]. In contrast, uncoagulated EDTA-blood cell pellet DNA from vitamin B12-deficient and B12-replete women exhibited similar percentage methylation (4.45±0.15%; n = 77 vs. 4.47±0.15%; n = 47) and [3 H]methyl-acceptance (27,378±4,094 cpm; n = 17 vs. 26,610±2,292 cpm; n = 11). Therefore, in simultaneously collected paired blood samples, vitamin B12-deficiency was associated with decreased DNA methylation only in coagulated samples. These findings highlight the importance of sample collection methods in epigenetic studies, and the potential impact biological processes can have on DNA methylation during collection.     

Systematic CpT (ApG) Depletion and CpG Excess Are Unique Genomic Signatures of Large DNA Viruses Infecting Invertebrates

Differences in the relative abundance of dinucleotides, if any may provide important clues on host-driven evolution of viruses. We studied dinucleotide frequencies of large DNA viruses infecting vertebrates (n = 105; viruses infecting mammals = 99; viruses infecting aves = 6; viruses infecting reptiles = 1) and invertebrates (n = 88; viruses infecting insects = 84; viruses infecting crustaceans = 4). We have identified systematic depletion of CpT(ApG) dinucleotides and over-representation of CpG dinucleotides as the unique genomic signature of large DNA viruses infecting invertebrates. Detailed investigation of this unique genomic signature suggests the existence of invertebrate host-induced pressures specifically targeting CpT(ApG) and CpG dinucleotides. The depletion of CpT dinucleotides among large DNA viruses infecting invertebrates is at least in part, explained by non-canonical DNA methylation by the infected host. Our findings highlight the role of invertebrate host-related factors in shaping virus evolution and they also provide the necessary framework for future studies on evolution, epigenetics and molecular biology of viruses infecting this group of hosts.     

Fetal DNA Methylation Associates with Early Spontaneous Preterm Birth and Gestational Age

Spontaneous preterm birth (PTB, <37 weeks gestation) is a major public health concern, and children born preterm have a higher risk of morbidity and mortality throughout their lives. Recent studies suggest that fetal DNA methylation of several genes varies across a range of gestational ages (GA), but it is not yet clear if fetal epigenetic changes associate with PTB. The objective of this study is to interrogate methylation patterns across the genome in fetal leukocyte DNA from African Americans with early PTB (24^1/7–34^0/7 weeks; N = 22) or term births (39^0/7–40^6/7weeks; N = 28) and to evaluate the association of each CpG site with PTB and GA. DNA methylation was assessed across the genome with the HumanMethylation450 BeadChip. For each individual sample and CpG site, the proportion of DNA methylation was estimated. The associations between methylation and PTB or GA were evaluated by fitting a separate linear model for each CpG site, adjusting for relevant covariates. Overall, 29 CpG sites associated with PTB (FDR<.05; 5.7×10⁻¹⁰<p<2.9×10⁻⁶) independent of GA. Also, 9637 sites associated with GA (FDR<.05; 9.5×10⁻¹⁶<p<1.0×10⁻³), with 61.8% decreasing in methylation with shorter GA. GA-associated CpG sites were depleted in the CpG islands of their respective genes (p<2.2×10⁻¹⁶). Gene set enrichment analysis (GSEA) supported enrichment of GA-associated CpG sites in genes that play a role in embryonic development as well as the extracellular matrix. Additionally, this study replicated the association of several CpG sites associated with gestational age in other studies (CRHBP, PIK3CD and AVP). Dramatic differences in fetal DNA methylation are evident in fetuses born preterm versus at term, and the patterns established at birth may provide insight into the long-term consequences associated with PTB.     

A Comprehensive Microarray-Based DNA Methylation Study of 367 Hematological Neoplasms

Background

Alterations in the DNA methylation pattern are a hallmark of leukemias and lymphomas. However, most epigenetic studies in hematologic neoplasms (HNs) have focused either on the analysis of few candidate genes or many genes and few HN entities, and comprehensive studies are required.

Methodology/Principal Findings

Here, we report for the first time a microarray-based DNA methylation study of 767 genes in 367 HNs diagnosed with 16 of the most representative B-cell (n = 203), T-cell (n = 30), and myeloid (n = 134) neoplasias, as well as 37 samples from different cell types of the hematopoietic system. Using appropriate controls of B-, T-, or myeloid cellular origin, we identified a total of 220 genes hypermethylated in at least one HN entity. In general, promoter hypermethylation was more frequent in lymphoid malignancies than in myeloid malignancies, being germinal center mature B-cell lymphomas as well as B and T precursor lymphoid neoplasias those entities with highest frequency of gene-associated DNA hypermethylation. We also observed a significant correlation between the number of hypermethylated and hypomethylated genes in several mature B-cell neoplasias, but not in precursor B- and T-cell leukemias. Most of the genes becoming hypermethylated contained promoters with high CpG content, and a significant fraction of them are targets of the polycomb repressor complex. Interestingly, T-cell prolymphocytic leukemias show low levels of DNA hypermethylation and a comparatively large number of hypomethylated genes, many of them showing an increased gene expression.

Conclusions/Significance

We have characterized the DNA methylation profile of a wide range of different HNs entities. As well as identifying genes showing aberrant DNA methylation in certain HN subtypes, we also detected six genes—DBC1, DIO3, FZD9, HS3ST2, MOS, and MYOD1—that were significantly hypermethylated in B-cell, T-cell, and myeloid malignancies. These might therefore play an important role in the development of different HNs.     

Meta-analysis of human methylation data for evidence of sex-specific autosomal patterns

Background

Several individual studies have suggested that autosomal CpG methylation differs by sex both in terms of individual CpG sites and global autosomal CpG methylation. However, these findings have been inconsistent and plagued by spurious associations due to the cross reactivity of CpG probes on commercial microarrays. We collectively analysed 76 published studies (n = 6,795) for sex-associated differences in both autosomal and sex chromosome CpG sites.

Results

Overall autosomal methylation profiles varied substantially by study, and we encountered substantial batch effects. We accounted for these by conducting random effects meta-analysis for individual autosomal CpG methylation associations. After excluding non-specific probes, we found 184 autosomal CpG sites differentially methylated by sex after correction for multiple testing. In line with previous studies, average beta differences were small. Many of the most significantly associated CpG probes were new. Of note was differential CpG methylation in the promoters of genes thought to be involved in spermatogenesis and male fertility, such as SLC9A2, SPESP1, CRISP2, and NUPL1. Pathway analysis revealed overrepresentation of genes differentially methylated by sex in several broad Gene Ontology biological processes, including RNA splicing and DNA repair.

Conclusions

This study represents a comprehensive analysis of sex-specific methylation patterns. We demonstrate the existence of sex-specific methylation profiles and report a large number of novel DNA methylation differences in autosomal CpG sites between sexes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-981) contains supplementary material, which is available to authorized users.     

Gene Expression Profiling in Human Lung Development: An Abundant Resource for Lung Adenocarcinoma Prognosis

A tumor can be viewed as a special “organ” that undergoes aberrant and poorly regulated organogenesis. Progress in cancer prognosis and therapy might be facilitated by re-examining distinctive processes that operate during normal development, to elucidate the intrinsic features of cancer that are significantly obscured by its heterogeneity. The global gene expression signatures of 44 human lung tissues at four development stages from Asian descent and 69 lung adenocarcinoma (ADC) tissue samples from ethnic Chinese patients were profiled using microarrays. All of the genes were classified into 27 distinct groups based on their expression patterns (named as PTN1 to PTN27) during the developmental process. In lung ADC, genes whose expression levels decreased steadily during lung development (genes in PTN1) generally had their expression reactivated, while those with uniformly increasing expression levels (genes in PTN27) had their expression suppressed. The genes in PTN1 contain many n-gene signatures that are of prognostic value for lung ADC. The prognostic relevance of a 12-gene demonstrator for patient survival was characterized in five cohorts of healthy and ADC patients [ADC_CICAMS (n = 69, p = 0.007), ADC_PNAS (n = 125, p = 0.0063), ADC_GSE13213 (n = 117, p = 0.0027), ADC_GSE8894 (n = 62, p = 0.01), and ADC_NCI (n = 282, p = 0.045)] and in four groups of stage I patients [ADC_CICAMS (n = 22, p = 0.017), ADC_PNAS (n = 76, p = 0.018), ADC_GSE13213 (n = 79, p = 0.02), and ADC_qPCR (n = 62, p = 0.006)]. In conclusion, by comparison of gene expression profiles during human lung developmental process and lung ADC progression, we revealed that the genes with a uniformly decreasing expression pattern during lung development are of enormous prognostic value for lung ADC.     

Identification of Highly Methylated Genes across Various Types of B-Cell Non-Hodgkin Lymphoma

Epigenetic alterations of gene expression are important in the development of cancer. In this study, we identified genes which are epigenetically altered in major lymphoma types. We used DNA microarray technology to assess changes in gene expression after treatment of 11 lymphoma cell lines with epigenetic drugs. We identified 233 genes with upregulated expression in treated cell lines and with downregulated expression in B-cell lymphoma patient samples (n = 480) when compared to normal B cells (n = 5). The top 30 genes were further analyzed by methylation specific PCR (MSP) in 18 lymphoma cell lines. Seven of the genes were methylated in more than 70% of the cell lines and were further subjected to quantitative MSP in 37 B-cell lymphoma patient samples (diffuse large B-cell lymphoma (activated B-cell like and germinal center B-cell like subtypes), follicular lymphoma and Burkitt`s lymphoma) and normal B lymphocytes from 10 healthy donors. The promoters of DSP, FZD8, KCNH2, and PPP1R14A were methylated in 28%, 67%, 22%, and 78% of the 36 tumor samples, respectively, but not in control samples. Validation using a second series of healthy donor controls (n = 42; normal B cells, peripheral blood mononuclear cells, bone marrow, tonsils and follicular hyperplasia) and fresh-frozen lymphoma biopsies (n = 25), confirmed the results. The DNA methylation biomarker panel consisting of DSP, FZD8, KCNH2, and PPP1R14A was positive in 89% (54/61) of all lymphomas. Receiver operating characteristic analysis to determine the discriminative power between lymphoma and healthy control samples showed a c-statistic of 0.96, indicating a possible role for the biomarker panel in monitoring of lymphoma patients.     

A Genome-Wide Methylation Study on Essential Hypertension in Young African American Males

Objective

There is emerging evidence from animal studies suggesting a key role for methylation in the pathogenesis of essential hypertension. However, to date, very few studies have investigated the role of methylation in the development of human hypertension, and none has taken a genome-wide approach. Based on the recent studies that highlight the involvement of inflammation in the development of hypertension, we hypothesize that changes in DNA methylation of leukocytes are involved in the pathogenesis of hypertension.

Method & Results

We conducted a genome-wide methylation analysis on 8 hypertensive cases and 8 normotensive age-matched controls aged 14–23 years and performed validation of the most significant CpG sites in 2 genes in an independent sample of 36 hypertensive cases and 60 normotensive controls aged 14–30 years. Validation of the CpG sites in the SULF1 gene was further conducted in a second replication sample of 36 hypertensive cases and 34 controls aged 15.8–40 years. A CpG site in the SULF1 gene showed higher methylation levels in cases than in healthy controls in the genome-wide step (p = 6.2×10⁻⁵), which was confirmed in the validation step (p = 0.011) for subjects ≤30 years old but was not significant for subjects of all ages combined (p = 0.095).

Conclusion

The identification of a difference in a blood leukocyte DNA methylation site between hypertensive cases and normotensive controls suggests that changes in DNA methylation may play an important role in the pathogenesis of hypertension. The age dependency of the effect further suggests complexity of epigenetic regulation in this age-related disease.     

Epigenome-Wide DNA Methylation in Hearing Ability: New Mechanisms for an Old Problem

Epigenetic regulation of gene expression has been shown to change over time and may be associated with environmental exposures in common complex traits. Age-related hearing impairment is a complex disorder, known to be heritable, with heritability estimates of 57–70%. Epigenetic regulation might explain the observed difference in age of onset and magnitude of hearing impairment with age. Epigenetic epidemiology studies using unrelated samples can be limited in their ability to detect small effects, and recent epigenetic findings in twins underscore the power of this well matched study design. We investigated the association between venous blood DNA methylation epigenome-wide and hearing ability. Pure-tone audiometry (PTA) and Illumina HumanMethylation array data were obtained from female twin volunteers enrolled in the TwinsUK register. Two study groups were explored: first, an epigenome-wide association scan (EWAS) was performed in a discovery sample (n = 115 subjects, age range: 47–83 years, Illumina 27 k array), then replication of the top ten associated probes from the discovery EWAS was attempted in a second unrelated sample (n = 203, age range: 41–86 years, Illumina 450 k array). Finally, a set of monozygotic (MZ) twin pairs (n = 21 pairs) within the discovery sample (Illumina 27 k array) was investigated in more detail in an MZ discordance analysis. Hearing ability was strongly associated with DNA methylation levels in the promoter regions of several genes, including TCF25 (cg01161216, p = 6.6×10⁻⁶), FGFR1 (cg15791248, p = 5.7×10⁻⁵) and POLE (cg18877514, p = 6.3×10⁻⁵). Replication of these results in a second sample confirmed the presence of differential methylation at TCF25 (p(replication) = 6×10⁻⁵) and POLE (p(replication) = 0.016). In the MZ discordance analysis, twins'' intrapair difference in hearing ability correlated with DNA methylation differences at ACP6 (cg01377755, r = −0.75, p = 1.2×10⁻⁴) and MEF2D (cg08156349, r = −0.75, p = 1.4×10⁻⁴). Examination of gene expression in skin, suggests an influence of differential methylation on expression, which may account for the variation in hearing ability with age.