Genome-wide DNA methylation analysis for diabetic nephropathy in type 1 diabetes mellitus
© Bell et al; licensee BioMed Central Ltd. 2010
Received: 4 February 2010
Accepted: 5 August 2010
Published: 5 August 2010
Diabetic nephropathy is a serious complication of diabetes mellitus and is associated with considerable morbidity and high mortality. There is increasing evidence to suggest that dysregulation of the epigenome is involved in diabetic nephropathy. We assessed whether epigenetic modification of DNA methylation is associated with diabetic nephropathy in a case-control study of 192 Irish patients with type 1 diabetes mellitus (T1D). Cases had T1D and nephropathy whereas controls had T1D but no evidence of renal disease.
We performed DNA methylation profiling in bisulphite converted DNA from cases and controls using the recently developed Illumina Infinium® HumanMethylation27 BeadChip, that enables the direct investigation of 27,578 individual cytosines at CpG loci throughout the genome, which are focused on the promoter regions of 14,495 genes.
Singular Value Decomposition (SVD) analysis indicated that significant components of DNA methylation variation correlated with patient age, time to onset of diabetic nephropathy, and sex. Adjusting for confounding factors using multivariate Cox-regression analyses, and with a false discovery rate (FDR) of 0.05, we observed 19 CpG sites that demonstrated correlations with time to development of diabetic nephropathy. Of note, this included one CpG site located 18 bp upstream of the transcription start site of UNC13B, a gene in which the first intronic SNP rs13293564 has recently been reported to be associated with diabetic nephropathy.
This high throughput platform was able to successfully interrogate the methylation state of individual cytosines and identified 19 prospective CpG sites associated with risk of diabetic nephropathy. These differences in DNA methylation are worthy of further follow-up in replication studies using larger cohorts of diabetic patients with and without nephropathy.
Diabetic nephropathy is a serious microvascular complication of diabetes and has become the most common cause of end-stage renal disease (ESRD) in many national ESRD registries [1, 2]. Approximately one third of diabetic individuals will develop clinically apparent nephropathy characterised by persistent proteinuria, hypertension and eventual progressive decline in glomerular filtration rate [3, 4]. Whilst chronic hyperglycaemia is integral to the pathogenesis of diabetic nephropathy there is also strong evidence for a genetic susceptibility to this common complication of diabetes [5–7]. In genetically susceptible individuals, prolonged hyperglycaemia leads to chronic metabolic and haemodynamic changes [8, 9] whose effects, including those driven by transforming growth factor beta 1 (TGFβ1), promote structural abnormalities in the kidney such as glomerular basement membrane thickening, podocyte injury, and mesangial matrix expansion, with the later development of irreversible glomerular sclerosis and tubulointerstitial fibrosis .
Genetic studies have found evidence for linkage to nephropathy at a number of chromosome loci including 2q, 3q22 and 19q [11, 12]. While no consistently replicated genetic associations have been identified, various candidate gene associations in diabetic nephropathy in type 1 diabetes (T1D) have been proposed, such as the SNP rs13293564 G/T substitution in intron 1 of UNC13B , and SOD1 SNPs . An association with the ELMO1 gene and diabetic nephropathy was initially identified in patients with type 2 diabetes mellitus  and this association has now also been reported for diabetic nephropathy in T1D . A recent genome-wide association study identified FRMD3 and CARS gene variants as possible genetic susceptibility factors for T1D diabetic nephropathy .
Epigenetic mechanisms allow alteration of genome function without mutating the underlying sequence. They involve the interacting actions of DNA methylation (the addition of a methyl group to the 5th carbon position of cytosine), histone modifications and non-coding RNAs . A number of indirect lines of evidence point to the involvement of epigenetic changes in diabetic nephropathy. Murine models of disease progression displaying temporal variation in gene expression have indicated these supra-sequence devices may be involved in the pathogenesis . Gene expression changes reflect dynamic alterations in gene transcription and also messenger RNA stability, which may be influenced by the epigenetic modification of the genome in response to chronic hyperglycaemic stress. Altered DNA methylation has been additionally implicated in vascular disease [20, 21]. Furthermore, characteristics observed in diabetic nephropathy such as hyperhomocysteinaemia, dyslipidaemia, inflammation and oxidative stress can promote aberrant DNA methylation [22–24].
DNA methylation in mammalian species, except for pluripotent stem cells [25, 26], appears to occur almost exclusively at consecutive CpG dinucleotides, of which there are approximately 28 million located within the human genome . The majority of cytosines in these dinucleotides remain methylated (~80% in Homo sapiens) , which contrasts with CpG islands (CGIs) in which cytosines are generally unmethylated. CGIs are regions that exceed genome average levels of CpG density and co-localise with 72% of gene promoters , notably in those of highly expressed genes.
In order to assess possible changes in CpG methylation in a genome-wide manner we utilised the high throughput Illumina Infinium® Human Methylation27 BeadChip DNA methylation array to investigate association with nephropathy in a case-control study of 192 T1D patients. These probes interrogate the methylation state of 27,578 individual CpGs located predominately in CGIs within proximal promoter regions, between 1.5 kb upstream and 1 kb downstream of the transcription start sites of 14,475 consensus coding sequence genes throughout the genome. Furthermore 110 miRNA promoters are targeted with 254 CpG loci probes . This platform enables the high-throughput array based investigation of individual CpGs at a magnitude higher level than the previous array based Golden-Gate assay .
The development, as well as the progression, of diabetic nephropathy has been linked to the same inflammatory cell processes that are implicated in the progression of T1D itself, such as the activated T cells involved in the destruction of islet β cells . Thus changes in the expression of cytokines and other inflammatory markers from distinct blood cell types make the investigation of methylation changes in DNA derived from blood a plausible option for examining epigenetic alterations in this disease. DNA methylation changes in peripheral blood have additionally been shown to indicate or be surrogate markers of active disease processes .
This study has investigated the possibility to identify DNA methylation biomarkers in peripheral blood cell-derived DNA that may be associated with the pathogenesis of T1D nephropathy. The study was designed to detect associations between the disease and differences in DNA methylation but not their functional relationship. Undertaken in a case control design, those with and without nephropathy were matched by age of onset and duration of diabetes.
Age at Diagnosis
Bisulphite conversion of DNA
Bisulphite conversion of 1 ug genomic DNA was performed using the EZ-96 DNA Methylation-Gold™ Kit (©Zymo Research Corp., USA) in two 96 well plates according to the manufacturer's instructions. Cases and controls were randomly located between the two plates. Effective bisulphite conversion was checked by successful PCR amplification with a pair of primers specific for converted DNA and unsuccessful amplification with a pair of primers for unconverted DNA, of 24 random samples across both plates (details available on request from authors).
Illumina Infinium methylation assay
DNA samples from cases and controls were interrogated utilising the Illumina Infinium® Human Methylation27 BeadChip. This platform detects the methylation status of 27,578 CpG sites by sequencing-based genotyping of bisulphite treated DNA. The chromosomal distribution and other statistics relevant to these CpG sites are provided in Additional File 1: Supplementary Figure S1. This chemical modification of DNA converts only unmethylated cytosines to uracils, thereafter allowing for highly multiplexed genotyping with single site resolution. Bisulphite converted and unconverted (i.e. methylated) sites are simultaneously evaluated by hybridisation of DNA to site-specific probes attached to beads, one for unmethylated and the other for methylated sites, followed by allele specific base extension that includes a fluorescent label. Two different labels are used, and fluorescent signals are specific for either the T (unmethylated) or C (methylated) alleles. Methylation scores represented as β values are generated for each site using BeadStudio 3.2 software (©Illumina, Inc. 2003-2008) and are computed based on the ratio of methylated to methylated plus unmethylated signal outputs. Thus the β values range from 0 (unmethylated) to 1 (fully methylated) on a continuous scale.
This platform for quantitative methylation, essentially an adaptation of the highly successful Illumina Infinium I Whole Genome Genotype SNP genotyping assay used extensively for GWAS [32, 33], is extremely accurate and powerful. Reproducibility estimates are extremely high with an average r2 correlation for β-values of 0.992 for 24 technical replicates . Due to the high density of probes up to 12 separate samples can be assayed on a single BeadChip, thus allowing for high-throughput processing.
The Illumina Infinium methylation assay was performed according to the manufacturer's instructions, and this along with the technical validation of the assay is detailed in Bibikova et al. . Bisulphite-converted DNA was denatured, amplified, fragmented and subsequently hybridised, with use of a specific hybridisation buffer, to the chip arrays (Sentrix positions A-L). Cases and controls were arranged to randomly approximate a 50:50 proportion on each chip. Primer extension was performed, then staining and coating, followed by imaging with the Illumina BeadArray reader. Internal quality controls included assessment of staining, hybridisation, target removal, extension, bisulphite conversion efficiency, dye specificity and additionally negative controls.
Initial array results were visualised using Illumina® BeadStudio 3.2 (© Illumina Inc 2003-2008). All computations and statistical analyses were performed using the R package (R 2.8.1) (http://www.r-project.org)  and Bioconductor .
Quality control and data normalisation
Quality control and normalisation of methylation data was undertaken using the analysis pipeline reported elsewhere . Briefly, samples were monitored for coverage (fraction of CpGs with detectable intensity values above background) and bisulphite conversion efficiency (BSCE) using the controls provided on the Illumina Beadchip. Of the 192 samples, one sample had failed BSCE and an additional four samples had relatively low coverage and BSCE. These samples were removed, leaving a total of 187 samples with approximately 81% global coverage (22,486 CpG sites). Intra-array normalisation was performed by the Illumina BeadStudio 3.2. software. The resulting β-valued data matrix was normalised further using a quantile-normalisation strategy, designed to reduce unwanted inter-array variation. The normalised and raw data are available from GEO (Gene Expression Omnibus, NCBI) under the accession number GSE20067.
Singular Value Decomposition and Significance of Singular Values
Association between CpG-valued methylation profiles and diabetic nephropathy was performed using a Cox-regression model where time to event equals the time between diagnosis of T1D to the onset of nephropathy. Multivariate regressions were performed for each CpG separately and included factors for chip, BSCE, sex, cohort and age at sample draw. To correct for multiple testing we estimated the false discovery rate (FDR) using the q-value framework . We have shown previously that the analytical q-value estimates are similar to those obtained using permutations of sample labels preserving the potential correlative structure between CpGs .
Methylation profile of arrays
DNA methylation analysis
Ilmn CpG ID
3.266 × 10-06
4.164 × 10-06
5.427 × 10-06
Negative regulation of transcription
8.703 × 10-06
1.457 × 10-05
1.644 × 10-05
2.486 × 10-05
Induction of apoptosis
2.798 × 10-05
2.976 × 10-05
3.056 × 10-05
3.441 × 10-05
NLS-bearing substrate import into nucleus
3.543 × 10-05
3.847 × 10-05
4.092 × 10-05
4.786 × 10-05
4.915 × 10-05
Regulation of calcium ion-dependent exocytosis
5.049 × 10-05
5.122 × 10-05
5.430 × 10-05
To investigate whether there was a significant co-directional change of CpGs within 1 kb of each other, the normalised results for each of those proximal pairs within this distance were combined and compared in a case versus control analysis, with two of the genes in Table 2, PPAPR3 (cg02590345 and cg16898420 - separated by 119 bp) and TRPS1 (cg12569516 and cg21312090 - separated by 534 bp) also being in the top ten results of this analysis (data not shown).
This study investigated the methylation state of approximately 14,000 CGIs in promoter regions of genes spread throughout the genome in a cohort of 192 type 1 diabetic patients, half of whom had developed diabetic nephropathy. A Cox-regression model was utilised whereby the effect of a genetic susceptibility has a multiplicative effect over time on the risk of a subject developing nephropathy. This analysis, after accounting for confounding effects, identified an excess of CpG sites which were observed to have a change in methylation status that could be correlated with duration of diabetes prior to onset of nephropathy with a FDR of 0.15. Of these observations 263 CpGs were identified as conferring increased risk and 162 CpGs conferring decreased risk for nephropathy, respectively and co-localised to 421 unique genic regions. However this dataset did not reveal any significant enrichment of particular physiological gene set pathways.
Using a more stringent FDR cut-off of 0.05, a set of 19 CpGs was identified (in 19 unique gene CGIs). One CpG in this list, cg07341907 is located 18 bp upstream of the transcription start site of a previously identified type 1 diabetic nephropathy gene, UNC13B . Levels of methylation for this CpG were shown to be slightly higher in the case group compared with controls. Whilst expression of this gene has been shown to be increased in models of nephropathy and hyperglycaemia  the expression effect of methylation changes can often be difficult to predict . The other CpG assayed within this genic region, cg02096633, is located ~0.5 kb further downstream, 197 bp into the first intron and shows consistently low methylation levels in both cases and controls (β-value averages 0.051 and 0.052 respectively). The susceptibility SNP, rs13293564 resides within the first intron and in Caucasians (from CEU HapMap data) is within a linkage disequilibrium block of 23 kb that includes the CGI. Tregouet et al. could determine no plausible function for the intronic SNP, however did identify five SNPs in strong linkage disequilibrium with this SNP that reside within the plausible promoter region . Two of which (rs10081672 and rs10972333) were found to affect potential transcription factor binding sites. Additionally, the former SNP (rs10081672) creates or abrogates a CpG site depending on which allele is present. Therefore, modification of the promoter region of this gene whether genetically or epigenetically or in combination could influence its expression by affecting transcription factor binding. Tregouet et al. proposed that UNC13B mediates apoptosis in glomerular cells due to hyperglycemia, and they therefore suggested that this association could indicate initiation of nephropathy .
Whilst the 27K Infinium assay enables high-throughput investigation of individual CpG sites with high resolution, and is a considerable improvement on the previous Golden-Gate methylation assay, this methodology still has some limitations. Only approximately 0.1% of the total number of CpGs within the genome are assessed, however the placement of these CpGs are predominately within the CpG Islands with high likelihood of critical effects on promoter activity. As these CpGs can significantly affect or silence expression, changes can be more deleterious than gene sequence variants, therefore increasing the potential power of the assay. Additionally not all CpG sites need to be assessed as correlation of CpG status between sites can extend up to 1 kb . Whist dramatic switching on or off facilitated by DNA methylation may occur within these promoter regions, more subtle and perhaps dynamic effects may however be missed that are occurring in the surrounding CGI shores (regions within 2 kb of the islands) . Changes in gene body methylation that may influence the expression of certain isoforms, by effecting the inclusion or exclusion of certain exons, are also not assayed. Additionally this array does have a bias towards cancer and imprinted genes; with ~200 of these having 3-20 CpG probes per promoter. Increasing knowledge with regards to the most biologically informative CpG sites and future up scaling will undoubtedly improve targeted platforms, although a more in depth investigation of all of the CpGs of the genome would require a differing approach. An enrichment technique such as MeDIP-Seq would gain a representative genome wide picture but this lacks the base pair resolution of this array . Isolated studies utilising full bisulphite sequencing of a very small number of genomes have recently been published , however scaling up to large numbers of cases and controls currently remains prohibitive due to cost.
Like all genome-wide association studies this study was designed to find statistically significant associations, in this case with DNA methylation variation not with genotype, but not to identify the underlying mechanism(s) for the cause or function of the observed variation. Although DNA methylation is tissue-specific, examination of peripheral whole blood was informative in this case in order to determine whether any genome-wide methylation signal changes could be detected in this easily accessible surrogate tissue. Whilst only a small number of significant sites were identified, due to the stringent False Discovery Rate cut-off implemented, this experiment did identify possible biomarkers in this pathogenic process with more easily implementable clinical utility potential. Therefore, this indicates that active disease processes can be identified in the DNA methylation pattern of peripheral blood and be a possible marker of these . The development, as well as the progression, of diabetic nephropathy has been linked to the same inflammatory cell changes that are implicated in the progression of T1D itself, such as the activated T cells involved in the destruction of islet β cells . Furthermore, additional inflammatory-related or other signals may be present in peripheral blood DNA that require isolation of specific cell types to detect and other methylation change signals may only be able to be identified in the affected renal cells themselves.
This study is one of the first to investigate a complex disease trait utilising this high-throughput DNA methylation 27 K array based assay. The effect of length of T1D status before developing or not developing nephropathy was observed to show a correlation with respect to methylation scores of individual CpG loci. By the examination of a cohort of 192 T1D patients this study has confirmed the utility of this approach to genome wide DNA methylation analysis and the potential prospects for larger and subsequent replication studies of epigenomic factors in common diseases.
List of Abbreviations
Bisulphite Conversion Efficiency
End-Stage Renal Disease
False Discovery Rate
Type 1 Diabetes Mellitus
Singular Value Decomposition.
We thank the Renal Unit Fund, Belfast Health and Social Care Trust for funding this project.
CGB and SB were funded by the Wellcome Trust (084071), and AET was supported by a Heller Research Fellowship. SB is a Royal Society-Wolfson Research Merit Award Holder. We also thank Dr Denise Sadlier, University College Dublin, for providing DNA samples from cases and controls from the Republic of Ireland.
We would further like to thank UCL Genomics for excellent technical support.
- Ansell D, Feehally J, Fogarty D, Tomsom CRV, Williams AJ, Warwick G: UK Renal Registry Annual Report. Bistol. 2008, 2009/06/30 ednGoogle Scholar
- U.S Renal Data System: US Renal Data System 2008 Annual Data Report. Atlas of Chronic Kidney Disease and End Stage Renal Disease in the United States. 2008, Bethesda, MD: National Institutes of Health, National Institutes of Diabetes and Digestive and Kidney Diseases, (National Institutes of Health NIoDaDaKDGoogle Scholar
- Ibrahim HN, Hostetter TH: Diabetic nephropathy. J Am Soc Nephrol. 1997, 8: 487-493.PubMedGoogle Scholar
- O'Connor AS, Schelling JR: Diabetes and the kidney. Am J Kidney Dis. 2005, 46: 766-773. 10.1053/j.ajkd.2005.05.032.View ArticlePubMedGoogle Scholar
- Seaquist ER, Goetz FC, Rich S, Barbosa J: Familial clustering of diabetic kidney disease. Evidence for genetic susceptibility to diabetic nephropathy. N Engl J Med. 1989, 320: 1161-1165. 10.1056/NEJM198905043201801.View ArticlePubMedGoogle Scholar
- Quinn M, Angelico MC, Warram JH, Krolewski AS: Familial factors determine the development of diabetic nephropathy in patients with IDDM. Diabetologia. 1996, 39: 940-945. 10.1007/BF00403913.View ArticlePubMedGoogle Scholar
- Ichinose K, Kawasaki E, Eguchi K: Recent advancement of understanding pathogenesis of type 1 diabetes and potential relevance to diabetic nephropathy. Am J Nephrol. 2007, 27: 554-564. 10.1159/000107758.View ArticlePubMedGoogle Scholar
- Weigert C, Sauer U, Brodbeck K, Pfeiffer A, Haring HU, Schleicher ED: AP-1 proteins mediate hyperglycemia-induced activation of the human TGF-beta1 promoter in mesangial cells. J Am Soc Nephrol. 2000, 11: 2007-2016.PubMedGoogle Scholar
- Katavetin P, Miyata T, Inagi R, Tanaka T, Sassa R, Ingelfinger JR, Fujita T, Nangaku M: High glucose blunts vascular endothelial growth factor response to hypoxia via the oxidative stress-regulated hypoxia-inducible factor/hypoxia-responsible element pathway. J Am Soc Nephrol. 2006, 17: 1405-1413. 10.1681/ASN.2005090918.View ArticlePubMedGoogle Scholar
- Navarro-Gonzalez JF, Mora-Fernandez C: The role of inflammatory cytokines in diabetic nephropathy. J Am Soc Nephrol. 2008, 19: 433-442. 10.1681/ASN.2007091048.View ArticlePubMedGoogle Scholar
- He B, Osterholm AM, Hoverfalt A, Forsblom C, Hjorleifsdottir EE, Nilsson AS, Parkkonen M, Pitkaniemi J, Hreidarsson A, Sarti C, et al: Association of genetic variants at 3 q22 with nephropathy in patients with type 1 diabetes mellitus. Am J Hum Genet. 2009, 84: 5-13. 10.1016/j.ajhg.2008.11.012.View ArticlePubMedPubMed CentralGoogle Scholar
- Rogus JJ, Poznik GD, Pezzolesi MG, Smiles AM, Dunn J, Walker W, Wanic K, Moczulski D, Canani L, Araki S, et al: High-density single nucleotide polymorphism genome-wide linkage scan for susceptibility genes for diabetic nephropathy in type 1 diabetes: discordant sibpair approach. Diabetes. 2008, 57: 2519-2526. 10.2337/db07-1086.View ArticlePubMedPubMed CentralGoogle Scholar
- Tregouet DA, Groop PH, McGinn S, Forsblom C, Hadjadj S, Marre M, Parving HH, Tarnow L, Telgmann R, Godefroy T, et al: G/T substitution in intron 1 of the UNC13B gene is associated with increased risk of nephropathy in patients with type 1 diabetes. Diabetes. 2008, 57: 2843-2850. 10.2337/db08-0073.View ArticlePubMedPubMed CentralGoogle Scholar
- Al-Kateb H, Boright AP, Mirea L, Xie X, Sutradhar R, Mowjoodi A, Bharaj B, Liu M, Bucksa JM, Arends VL, et al: Multiple superoxide dismutase 1/splicing factor serine alanine 15 variants are associated with the development and progression of diabetic nephropathy: the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications Genetics study. Diabetes. 2008, 57: 218-228. 10.2337/db07-1059.View ArticlePubMedGoogle Scholar
- Leak TS, Perlegas PS, Smith SG, Keene KL, Hicks PJ, Langefeld CD, Mychaleckyj JC, Rich SS, Kirk JK, Freedman BI, et al: Variants in intron 13 of the ELMO1 gene are associated with diabetic nephropathy in African Americans. Ann Hum Genet. 2009, 73: 152-159. 10.1111/j.1469-1809.2008.00498.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Pezzolesi MG, Katavetin P, Kure M, Poznik GD, Skupien J, Mychaleckyj JC, Rich SS, Warram JH, Krolewski AS: confirmation of genetic associations at ELMO1 in the GoKinD collection supports its role as a susceptibility gene in diabetic nephropathy. Diabetes. 2009, 58: 2698-2702. 10.2337/db09-0641.View ArticlePubMedPubMed CentralGoogle Scholar
- Pezzolesi MG, Poznik GD, Mychaleckyj JC, Paterson AD, Barati MT, Klein JB, Ng DP, Placha G, Canani LH, Bochenski J, et al: Genome-wide association scan for diabetic nephropathy susceptibility genes in type 1 diabetes. Diabetes. 2009, 58: 1403-1410. 10.2337/db08-1514.View ArticlePubMedPubMed CentralGoogle Scholar
- Richards EJ: Inherited epigenetic variation--revisiting soft inheritance. Nat Rev Genet. 2006, 7: 395-401. 10.1038/nrg1834.View ArticlePubMedGoogle Scholar
- Makino H, Miyamoto Y, Sawai K, Mori K, Mukoyama M, Nakao K, Yoshimasa Y, Suga S: Altered gene expression related to glomerulogenesis and podocyte structure in early diabetic nephropathy of db/db mice and its restoration by pioglitazone. Diabetes. 2006, 55: 2747-2756. 10.2337/db05-1683.View ArticlePubMedGoogle Scholar
- Krause B, Sobrevia L, Casanello P: Epigenetics: new concepts of old phenomena in vascular physiology. Curr Vasc Pharmacol. 2009, 7: 513-520. 10.2174/157016109789043883.View ArticlePubMedGoogle Scholar
- Turunen MP, Aavik E, Yla-Herttuala S: Epigenetics and atherosclerosis. Biochim Biophys Acta. 2009, 1790: 886-891.View ArticlePubMedGoogle Scholar
- Ingrosso D, Perna AF: Epigenetics in hyperhomocysteinemic states. A special focus on uremia. Biochim Biophys Acta. 2009, 1790: 892-899.View ArticlePubMedGoogle Scholar
- Ekstrom TJ, Stenvinkel P: The epigenetic conductor: a genomic orchestrator in chronic kidney disease complications?. J Nephrol. 2009, 22: 442-449.PubMedGoogle Scholar
- Lund G, Zaina S: Atherosclerosis risk factors can impose aberrant DNA methylation patterns: a tale of traffic and homocysteine. Curr Opin Lipidol. 2009, 20: 448-449. 10.1097/MOL.0b013e3283309928.View ArticlePubMedGoogle Scholar
- Ramsahoye BH, Biniszkiewicz D, Lyko F, Clark V, Bird AP, Jaenisch R: Non-CpG methylation is prevalent in embryonic stem cells and may be mediated by DNA methyltransferase 3a. Proc Natl Acad Sci USA. 2000, 97: 5237-5242. 10.1073/pnas.97.10.5237.View ArticlePubMedPubMed CentralGoogle Scholar
- Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, Nery JR, Lee L, Ye Z, Ngo QM, et al: Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009, 462: 315-322. 10.1038/nature08514.View ArticlePubMedPubMed CentralGoogle Scholar
- Beck S, Rakyan VK: The methylome: approaches for global DNA methylation profiling. Trends Genet. 2008, 24: 231-237. 10.1016/j.tig.2008.01.006.View ArticlePubMedGoogle Scholar
- Saxonov S, Berg P, Brutlag DL: A genome-wide analysis of CpG dinucleotides in the human genome distinguishes two distinct classes of promoters. Proc Natl Acad Sci USA. 2006, 103: 1412-1417. 10.1073/pnas.0510310103.View ArticlePubMedPubMed CentralGoogle Scholar
- Bibikova M, Le J, Barnes B, Saedinia-Melnyk S, Zhou L, Shen R, Gunderson K: Genome-wide DNA methylation profiling using Infinium assay. Epigenomics. 2009, 1: 177-200. 10.2217/epi.09.14.View ArticlePubMedGoogle Scholar
- Bibikova M, Fan JB: GoldenGate assay for DNA methylation profiling. Methods Mol Biol. 2009, 507: 149-163. full_text.View ArticlePubMedGoogle Scholar
- Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Gayther SA, Apostolidou S, Jones A, Lechner M, Beck S, Jacobs IJ, Widschwendter M: An epigenetic signature in peripheral blood predicts active ovarian cancer. PLoS One. 2009, 4: e8274-10.1371/journal.pone.0008274.View ArticlePubMedPubMed CentralGoogle Scholar
- Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, Jonasdottir A, Sigurdsson A, Baker A, Palsson A, Masson G, et al: A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science. 2007, 316: 1491-1493. 10.1126/science.1142842.View ArticlePubMedGoogle Scholar
- Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, Boutin P, Vincent D, Belisle A, Hadjadj S, et al: A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 2007, 445: 881-885. 10.1038/nature05616.View ArticlePubMedGoogle Scholar
- Ihaka R, Gentleman R: R: A Language for Data Analysis and Graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314. 10.2307/1390807.Google Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.View ArticlePubMedPubMed CentralGoogle Scholar
- Hastie T, Tibshirani R, Friedman JH: The elements of statistical learning: data mining, inference, and prediction: with 200 full-color illustrations. 2001, New York: SpringerGoogle Scholar
- Leek JT, Storey JD: Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007, 3: 1724-1735. 10.1371/journal.pgen.0030161.View ArticlePubMedGoogle Scholar
- Buja A, Eyuboglu N: Remarks on Parallel Analysis. Multivariate Behavioral Research. 1992, 27: 509-540. 10.1207/s15327906mbr2704_2.View ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.View ArticlePubMedPubMed CentralGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.View ArticlePubMedPubMed CentralGoogle Scholar
- Goldenberg NM, Silverman M: Rab34 and its effector munc13-2 constitute a new pathway modulating protein secretion in the cellular response to hyperglycemia. Am J Physiol Cell Physiol. 2009, 297: C1053-1058. 10.1152/ajpcell.00286.2009.View ArticlePubMedGoogle Scholar
- Suzuki MM, Bird A: DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008, 9: 465-476. 10.1038/nrg2341.View ArticlePubMedGoogle Scholar
- Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, Burger M, Burton J, Cox TV, Davies R, Down TA, et al: DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006, 38: 1378-1385. 10.1038/ng1909.View ArticlePubMedPubMed CentralGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, et al: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178-186. 10.1038/ng.298.View ArticlePubMedPubMed CentralGoogle Scholar
- Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, Graf S, Johnson N, Herrero J, Tomazou EM, et al: A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis. Nat Biotechnol. 2008, 26: 779-785. 10.1038/nbt1414.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/3/33/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.