Epigenetic marks such as cytosine methylation are important determinants of cellular and whole-body phenotypes. However, the extent of, and reasons for inter-individual differences in cytosine methylation, and their association with phenotypic variation are poorly characterised. Here we present the first genome-wide study of cytosine methylation at single-nucleotide resolution in an animal model of human disease. We used whole-genome bisulfite sequencing in the spontaneously hypertensive rat (SHR), a model of cardiovascular disease, and the Brown Norway (BN) control strain, to define the genetic architecture of cytosine methylation in the mammalian heart and to test for association between methylation and pathophysiological phenotypes. Analysis of 10.6 million CpG dinucleotides identified 77,088 CpGs that were differentially methylated between the strains. In F1 hybrids we found 38,152 CpGs showing allele-specific methylation and 145 regions with parent-of-origin effects on methylation. Cis-linkage explained almost 60% of inter-strain variation in methylation at a subset of loci tested for linkage in a panel of recombinant inbred (RI) strains. Methylation analysis in isolated cardiomyocytes showed that in the majority of cases methylation differences in cardiomyocytes and non-cardiomyocytes were strain-dependent, confirming a strong genetic component for cytosine methylation. We observed preferential nucleotide usage associated with increased and decreased methylation that is remarkably conserved across species, suggesting a common mechanism for germline control of inter-individual variation in CpG methylation. In the RI strain panel, we found significant correlation of CpG methylation and levels of serum chromogranin B (CgB), a proposed biomarker of heart failure, which is evidence for a link between germline DNA sequence variation, CpG methylation differences and pathophysiological phenotypes in the SHR strain. Together, these results will stimulate further investigation of the molecular basis of locally regulated variation in CpG methylation and provide a starting point for understanding the relationship between the genetic control of CpG methylation and disease phenotypes.

Genomics Core Laboratory MRC Clinical Sciences Centre London United Kingdom

Institute of Physiology Academy of Sciences of the Czech Republic Prague Czech Republic; Institute of Biology and Medical Genetics 1st Medical Faculty Charles University Prague Czech Republic

Integrative Genomics and Medicine Group MRC Clinical Sciences Centre London United Kingdom

Physiological Genomics and Medicine Group MRC Clinical Sciences Centre London United Kingdom

Physiological Genomics and Medicine Group MRC Clinical Sciences Centre London United Kingdom; Institute of Clinical Sciences Imperial College London United Kingdom

Physiological Genomics and Medicine Group MRC Clinical Sciences Centre London United Kingdom; National Heart and Lung Institute Imperial College London United Kingdom

Zobrazit více v PubMed

Li E, Bestor TH, Jaenisch R (1992) Targeted mutation of the DNA methyltransferase gene results in embryonic lethality. Cell 69: 915–926. PubMed

Okano M, Bell DW, Haber DA, Li E (1999) DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell 99: 247–257. PubMed

Bird A (2002) DNA methylation patterns and epigenetic memory. Genes Dev 16: 6–21. PubMed

Kerkel K, Spadola A, Yuan E, Kosek J, Jiang L, et al. (2008) Genomic surveys by methylation-sensitive SNP analysis identify sequence-dependent allele-specific DNA methylation. Nat Genet 40: 904–908. PubMed

Schilling E, El Chartouni C, Rehli M (2009) Allele-specific DNA methylation in mouse strains is mainly determined by cis-acting sequences. Genome Res 19: 2028–2035. PubMed PMC

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. PubMed PMC

Schalkwyk LC, Meaburn EL, Smith R, Dempster EL, Jeffries AR, et al. (2010) Allelic skewing of DNA methylation is widespread across the genome. Am J Hum Genet 86: 196–212. PubMed PMC

Shoemaker R, Deng J, Wang W, Zhang K (2010) Allele-specific methylation is prevalent and is contributed by CpG-SNPs in the human genome. Genome Res 20: 883–889. PubMed PMC

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. PubMed PMC

Gertz J, Varley KE, Reddy TE, Bowling KM, Pauli F, et al. (2011) Analysis of DNA methylation in a three-generation family reveals widespread genetic influence on epigenetic regulation. PLoS Genet 7: e1002228. PubMed PMC

Lienert F, Wirbelauer C, Som I, Dean A, Mohn F, et al. (2011) Identification of genetic elements that autonomously determine DNA methylation states. Nat Genet 43: 1091–1097. PubMed

Xie W, Barr CL, Kim A, Yue F, Lee AY, et al. (2012) Base-resolution analyses of sequence and parent-of-origin dependent DNA methylation in the mouse genome. Cell 148: 816–831. PubMed PMC

Schmitz RJ, He Y, Valdes-Lopez O, Khan SM, Joshi T, et al. (2013) Epigenome-wide inheritance of cytosine methylation variants in a recombinant inbred population. Genome Res 23: 1663–1674. PubMed PMC

Hubner N, Wallace CA, Zimdahl H, Petretto E, Schulz H, et al. (2005) Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease. Nat Genet 37: 243–253. PubMed

Petretto E, Mangion J, Dickens NJ, Cook SA, Kumaran MK, et al. (2006) Heritability and tissue specificity of expression quantitative trait loci. PLoS Genet 2: e172. PubMed PMC

Atanur SS, Birol I, Guryev V, Hirst M, Hummel O, et al. (2010) The genome sequence of the spontaneously hypertensive rat: Analysis and functional significance. Genome Res 20: 791–803. PubMed PMC

Simonis M, Atanur SS, Linsen S, Guryev V, Ruzius FP, et al. (2012) Genetic basis of transcriptome differences between the founder strains of the rat HXB/BXH recombinant inbred panel. Genome Biol 13: r31. PubMed PMC

Monti J, Fischer J, Paskas S, Heinig M, Schulz H, et al. (2008) Soluble epoxide hydrolase is a susceptibility factor for heart failure in a rat model of human disease. Nat Genet 40: 529–537. PubMed PMC

Petretto E, Sarwar R, Grieve I, Lu H, Kumaran MK, et al. (2008) Integrated genomic approaches implicate osteoglycin (Ogn) in the regulation of left ventricular mass. Nat Genet 40: 546–552. PubMed PMC

Pravenec M, Churchill PC, Churchill MC, Viklicky O, Kazdova L, et al. (2008) Identification of renal Cd36 as a determinant of blood pressure and risk for hypertension. Nat Genet 40: 952–954. PubMed

Heinig M, Petretto E, Wallace C, Bottolo L, Rotival M, et al. (2010) A trans-acting locus regulates an anti-viral expression network and type 1 diabetes risk. Nature 467: 460–464. PubMed PMC

McDermott-Roe C, Ye J, Ahmed R, Sun XM, Serafin A, et al. (2011) Endonuclease G is a novel determinant of cardiac hypertrophy and mitochondrial function. Nature 478: 114–118. PubMed PMC

Pravenec M, Klir P, Kren V, Zicha J, Kunes J (1989) An analysis of spontaneous hypertension in spontaneously hypertensive rats by means of new recombinant inbred strains. J Hypertens 7: 217–221. PubMed

Aitman TJ, Critser JK, Cuppen E, Dominiczak A, Fernandez-Suarez XM, et al. (2008) Progress and prospects in rat genetics: a community view. Nat Genet 40: 516–522. PubMed

Bailey TL, Boden M, Buske FA, Frith M, Grant CE, et al. (2009) MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res 37: W202–208. PubMed PMC

Alisch RS, Barwick BG, Chopra P, Myrick LK, Satten GA, et al. (2012) Age-associated DNA methylation in pediatric populations. Genome Res 22: 623–632. PubMed PMC

Sati S, Tanwar VS, Kumar KA, Patowary A, Jain V, et al. (2012) High resolution methylome map of rat indicates role of intragenic DNA methylation in identification of coding region. PLoS One 7: e31621. PubMed PMC

Haas J, Frese KS, Park YJ, Keller A, Vogel B, et al. (2013) Alterations in cardiac DNA methylation in human dilated cardiomyopathy. EMBO Mol Med 5: 413–429. PubMed PMC

Movassagh M, Choy MK, Knowles DA, Cordeddu L, Haider S, et al. (2011) Distinct epigenomic features in end-stage failing human hearts. Circulation 124: 2411–2422. PubMed PMC

Rosjo H, Husberg C, Dahl MB, Stridsberg M, Sjaastad I, et al. (2010) Chromogranin B in heart failure: a putative cardiac biomarker expressed in the failing myocardium. Circ Heart Fail 3: 503–511. PubMed

Bell JT, Pai AA, Pickrell JK, Gaffney DJ, Pique-Regi R, et al. (2011) DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines. Genome Biol 12: R10. PubMed PMC

Rao RH (1993) Insulin resistance in spontaneously hypertensive rats. Difference in interpretation based on insulin infusion rate or on plasma insulin in glucose clamp studies. Diabetes 42: 1364–1371. PubMed

Paternostro G, Clarke K, Heath J, Seymour AM, Radda GK (1995) Decreased GLUT-4 mRNA content and insulin-sensitive deoxyglucose uptake show insulin resistance in the hypertensive rat heart. Cardiovasc Res 30: 205–211. PubMed

Aitman TJ, Glazier AM, Wallace CA, Cooper LD, Norsworthy PJ, et al. (1999) Identification of Cd36 (Fat) as an insulin-resistance gene causing defective fatty acid and glucose metabolism in hypertensive rats. Nat Genet 21: 76–83. PubMed

Price AL, Helgason A, Thorleifsson G, McCarroll SA, Kong A, et al. (2011) Single-tissue and cross-tissue heritability of gene expression via identity-by-descent in related or unrelated individuals. PLoS Genet 7: e1001317. PubMed PMC

Shen-Orr SS, Tibshirani R, Khatri P, Bodian DL, Staedtler F, et al. (2010) Cell type-specific gene expression differences in complex tissues. Nat Methods 7: 287–289. PubMed PMC

Siedlecka U, Arora M, Kolettis T, Soppa GK, Lee J, et al. (2008) Effects of clenbuterol on contractility and Ca2+ homeostasis of isolated rat ventricular myocytes. Am J Physiol Heart Circ Physiol 295: H1917–1926. PubMed PMC

Johnson MD, Mueller M, Game L, Aitman TJ (2012) Single nucleotide analysis of Cytosine methylation by whole-genome shotgun bisulfite sequencing. Curr Protoc Mol Biol Chapter 21: Unit21 23. PubMed

Karimi M, Johansson S, Ekstrom TJ (2006) Using LUMA: a Luminometric-based assay for global DNA-methylation. Epigenetics 1: 45–48. PubMed

Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754–1760. PubMed PMC

McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, et al. (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20: 1297–1303. PubMed PMC

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, et al. (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25: 2078–2079. PubMed PMC

Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological 57: 289–300.

Ward JH (1963) Hierarchical grouping to optimize an objective function. J Am Statist Assoc 48: 236–244.

R Development Core Team (2010) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing.

Saar K, Beck A, Bihoreau MT, Birney E, Brocklebank D, et al. (2008) SNP and haplotype mapping for genetic analysis in the rat. Nat Genet 40: 560–566. PubMed PMC

Churchill GA, Doerge RW (1994) Empirical threshold values for quantitative trait mapping. Genetics 138: 963–971. PubMed PMC

Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11: R106. PubMed PMC

Erill I, O'Neill MC (2009) A reexamination of information theory-based methods for DNA-binding site identification. BMC bioinformatics 10: 57. PubMed PMC

Quinlan AR, Hall IM (2010) BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26: 841–842. PubMed PMC

Huang da W, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 44–57. PubMed

Huang da W, Sherman BT, Lempicki RA (2009) Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37: 1–13. PubMed PMC

Supek F, Bosnjak M, Skunca N, Smuc T (2011) REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One 6: e21800. PubMed PMC

Cap analysis of gene expression reveals alternative promoter usage in a rat model of hypertension

Genetic Complementation of ATP Synthase Deficiency Due to Dysfunction of TMEM70 Assembly Factor in Rat

Systems Genetics Approaches in Rat Identify Novel Genes and Gene Networks Associated With Cardiac Conduction