Non-imprinted allele-specific DNA methylation on human autosomes
© Zhang et al.; licensee BioMed Central Ltd. 2009
Received: 17 July 2009
Accepted: 3 December 2009
Published: 3 December 2009
Differential DNA methylation between alleles is well established in imprinted genes and the X chromosomes in females but has rarely been reported at non-imprinted loci on autosomes.
We studied DNA methylation of cytosine-guanine dinucleotide (CpG) islands on chromosome 21 in leukocytes from several healthy individuals and observed novel cases of pronounced differential methylation of alleles. Allele-specific methylation affected complete CpG islands with methylation differences between alleles of up to 85%. The methylation differences between alleles were strongly correlated with the genotypes, excluding a connection to imprinting. We show that allele-specific methylation can lead to allelic repression of the methylated gene copy. Based on our results, allele-specific methylation is likely to affect about 10% of all human genes and to contribute to allele-specific expression and monoallelic gene silencing. Therefore, allele-specific methylation represents an epigenetic pathway of how genetic polymorphisms may lead to phenotypic variability. In most cases, we observed that some, but not all, heterozygous individuals showed allele-specific methylation, suggesting that allele-specific methylation is the outcome of an epigenetic drift, the direction of which is determined by the genetic differences between the alleles. We could show that the tendency to acquire hypermethylation in one allele was inherited.
We observed that larger differences in methylation levels between individuals were often coupled to allele-specific methylation and genetic polymorphisms, suggesting that the inter-individual variability of DNA methylation is strongly influenced by genetic differences. Therefore, genetic differences must be taken into account in future comparative DNA methylation studies.
DNA methylation is a major epigenetic process that plays essential roles in gene expression regulation, development and disease [1–3]. In mammals, differential DNA methylation between alleles occurs in imprinted genes [4, 5] and on the female X chromosomes [6, 7]. So far, there have been few reports about allele-specific methylation (ASM) on autosomes not connected to the parental inheritance of the alleles [8–10]. In imprinting and X-chromosome inactivation, ASM leads to monoallelic expression of genes, and in some cases of non-imprinted, autosomal ASM, a correlation between DNA methylation and allele-specific gene expression has been documented as well [8, 11].
In a previous work, we discovered ASM of three cytosine-guanine dinucleotide (CpG)-rich regions in gene promoters in leukocyte DNA derived from a healthy individual using bisulfite conversion, subcloning and sequencing . Because of the important potential contribution of ASM to gene expression and disease , we initiated a larger survey to find more examples of ASM and to understand if the ASM of these regions relates to gene imprinting or is sequence-dependent. Therefore, we studied the methylation pattern of 16 CpG-rich regions in gene promoters of chromosome 21 in up to 38 individuals by bisulfite conversion, subcloning and sequencing of individual clones. Additionally, we checked the inter-individual DNA methylation difference at these gene promoters with respect to a potential link to age and gender differences.
Allele-specific methylation of amplicons among individuals and correlation with genotype
Genes and amplicons
SNP (strand) dbsnp build 129
A/C (-) rs2297246
A/C (+) rs56270809
A/G (-) rs2838343
A/G (+) rs55860816
del/no del (+) rs41563015
C/G (+) rs25678
Individuals with SNP
Individuals with SNP and ASM
Coupling of ASM and genotype
Always: C hyper- methylated
Always: C hyper- methylated
18 cases: G hyper- methylated in right part of amplicon
Always: G hyper- methylated
Always: del hyper- methylated
Always: C hyper- methylated
1 case: A hyper- methylated in left part of amplicon
Average allelic methylation difference in ASM cases
54% (G hyper- methylated), 64% (A hyper- methylated)
Individuals with SNP but without ASM
Individuals without SNP
Inter-individual variation of DNA methylation
We observed that amplicons with a high methylation variance among individuals often showed a biphasic distribution of methylation comprising clones that were mainly unmethylated or highly methylated. This led to the question of whether the DNA methylation varied from one cell to another or between the two alleles in each cell. In some of these cases, the alleles could be differentiated by the presence of a single nucleotide polymorphism (SNP) in the sequenced region such that the sequencing reads could be sorted according to the allelic origin. In many of these cases we observed a clear correlation of genotype and methylation levels, indicating the presence of ASM (Table 1). To exclude that clonal PCR after bisulfite conversion might have caused these results, we performed individual bisulfite genomic sequencing analysis twice for cases where ASM was observed and always obtained similar results (Additional file 3). In addition, we used digestion with methylation-sensitive restriction enzymes to confirm the ASM (Additional file 3). For amplicons with ASM, we included more individuals in the analysis.
Allele-specific methylation and expression of the first exon of the C21orf81 gene
The methylation pattern of amplicon 23_2 was studied in 19 individuals, including one (called 'N' in Figure 2 and the following) previously analyzed . Altogether, eight individuals had an A/C SNP in this region and seven of them showed strong ASM (Table 1; Figure 2d; Additional file 4). In all ASM cases, the A allele was less methylated than the C allele (Figure 2e). For the 11 homozygous individuals, biphasic DNA methylation was observed in 7 cases.
We identified the correlation between the SNPs in amplicon 23_1 and 23_2 in the genomic DNA sequence of all individuals with ASM. When both SNPs were present, we always observed A-A and C-C haplotypes. In three such cases (individuals 9, 11 and 17), the A alleles in both amplicons were hypomethylated, indicating that pronounced ASM can span a whole CpG island (Figure 2a; Additional file 4). In individuals 3 and N, ASM was observed in one amplicon (either 23_1 or 23_2) and the other amplicon showed biphasic methylation, but did not contain the SNP, suggesting that the whole CpG island has ASM as well in these cases. However, there were also cases in which the methylation states of alleles in these two amplicons were not coupled. In individuals 4, 12 and 20, ASM was observed in 23_2, but there was no obvious methylation difference between alleles in amplicon 23_1. No methylation difference between alleles in both amplicons was observed in one individual (individual 1).
To determine if ASM leads to monoallelic expression of the C21orf81 gene, we prepared RNA from leukocytes of individuals 9 and 11, which showed strong ASM of both amplicons, and converted it into cDNA. Sequencing of the cDNA revealed that the gene was exclusively expressed from the unmethylated A allele, indicating the ASM correlates with allele-specific gene expression (Figure 2a).
We were able to check the methylation pattern of amplicon 23_1 in the parents of individual 17 (male). Like individual 17, the mother (17_M) had an A/C SNP with the C allele methylated. The father (17_P) had an A/A genotype and showed slight biphasic methylation but did not show clones with high methylation characteristic for C alleles (Additional file 4).
Allele-specific methylation of a CpG island in the HSF2BPgene
We investigated the methylation levels of amplicon 262 in homozygous individuals and observed a strong correlation between DNA methylation and genotype - individuals with A/A genotype were always low methylated while individuals with G/G were highly methylated (Figure 3d). The methylation levels of the A and G alleles from heterozygous individuals were similar to those of individuals with A/A or G/G genotypes, respectively (Figure 3d), suggesting that the methylation difference is determined mainly by the polymorphisms between the alleles. A similar observation was made for amplicons 23_1 and 23_2 (Additional file 6).
Allele-specific methylation of an RIPK4intron
Allele-specific methylation of amplicons in the CBR1gene
We previously identified ASM of two amplicons (176_1 and 176_2) of the CBR1 gene in blood derived from one individual . The allele with a deletion in 176_1 (dbsnp rs41563015) and with C at a C/G SNP site in 176_2 (dbsnp rs25678) was hypermethylated when compared with the allele without deletion and with G at the SNP site. Here, we studied the methylation pattern of these two regions in 20 additional individuals (Additional file 8). For amplicon 176_1, none of these individuals had the deletion or any other SNP in the sequence and all showed low methylation. For amplicon 176_2, 7 out of 20 individuals were heterozygous, but both alleles were always unmethylated. This result suggests that the deletion was more important for ASM to occur in the CBR1 gene promoter than the SNP in amplicon 176_2.
Allele-specific methylation on autosomes is well established in imprinted genes, which carry a differentially methylated region that is methylated depending on the parental origin of the allele . The sizes of such regions vary between hundreds to thousands of base pairs and the methylation difference between alleles can be up to 100%, in particular on maternally methylated alleles. Methylation leads to the exclusive expression of the maternal or paternal copy of the gene.
Recently, the allelic methylation difference of several CpG sites next to SNPs in non-imprinted regions of the genome was reported . Here, we focused on whole CpG islands close to the start sites of genes and identified novel regions of ASM, which all contain a large number of CpG sites and show strong differences in average methylation of up to 85% between the alleles. In some cases the entire promoter of a gene was subject to ASM. The allele-specific differences in methylation observed here are much greater than in previously published cases.
In all cases of ASM described here, there was a full correlation between the average allelic methylation level and allelic genotype, meaning that an allele with one particular genotype was always hypermethylated in different individuals. This observation is not compatible with a connection between ASM and imprinting. We conclude that the ASM identified here is due to genetic polymorphisms between two alleles, providing examples of genotype-epigenotype interactions.
In the current study we investigate ASM in leukocytes. We previously showed biphasic methylation of amplicon 232 in HEK293 cells, but in that sample no SNP was available to discriminate the alleles and confirm the presence of ASM (Additional file 9) . We did not find evidence for ASM in any of the other amplicons in HEK293 cells, HEPG2 cells or human fibroblasts in our previous study (Additional file 9) . Although the methylation status of cultured cells may not fully reflect the methylation status of the original tissue, this finding suggests that tissue specific factors play a role in the generation of ASM. There are several potential mechanisms for how genetic differences could specifically trigger the loss or gain of DNA methylation in one allele, such as binding sites of DNA binding proteins that may be present on one allele but not on the other. Predictions of transcription factor binding sites suggested for all amplicons that the different alleles might specifically interact with some transcription factors (Additional file 10), similar to what has been described recently for one case of cancer-associated ASM .
Our data allow for a rough estimation of how prevalent ASM is in the human epigenome. We identified 4 regions of strong ASM selected out of 190 genes studied for DNA methylation in our previous study . If one considers that we found SNPs allowing us to discriminate the alleles in about 20% of all amplicons that were carefully inspected, one may estimate ASM to occur in about 10% of all genes in leukocytes. The finding that about 11% of all amplicons in different tissues showed a biphasic distribution of methylation  leads to a similar estimate. We conclude, therefore, that ASM likely affects many genes in the human genome.
Recently, allele-specific gene expression has been identified to be a widespread phenomenon in human cells [15–18]. The ASM may correlate with allele-specific gene expression [8, 11] and we also identified one example of such correlation in our study. Aberrant monoallelic silencing of genes contributes to cancer and other diseases [19, 20] and allele specific hypermethylation of genes has been associated with cancer as well . Therefore, ASM may represent an important epigenetic pathway connecting genetic polymorphisms and phenotypic variability [20, 21]. In the future, more genome-wide bisulfite sequencing DNA methylation analysis with DNA from different individuals will be required to uncover the full impact of genetic polymorphisms on DNA methylation and gene regulation for development and disease processes.
We observed pronounced methylation difference between alleles spanning whole CpG islands in the non-imprinting loci of human autosomes, and proved a correlation with allelic gene expression level in one case. We further confirmed that the methylation differences between alleles were strongly correlated with genetic polymorphisms, which can further influence the inter-individual variability of DNA methylation. Our results suggest the ASM can affect many genes in the human genome and the genetic differences must be taken into account in future comparative DNA methylation studies.
Materials and methods
Genomic DNA extraction
Blood DNA was collected from 37 apparently healthy individuals, including 36 of caucasian origin, and 1 from Asia (Additional file 2). All individuals were informed about the study and declared consent in written form. Genomic DNA from blood was extracted using a QIAamp DNA Blood Mini Kit (Qiagen, Hilton, Germany).
DNA methylation analysis
DNA methylation analysis was carried out as described . Briefly, 200 to 300 ng genomic DNA was converted by bisulfite and used as template for PCR with primer designed in the region of interest. Information on amplicons selected and primer sequences are listed in Additional file 1. Around 40 clones for each amplicon were picked and sequenced using ABI BigDye Terminator chemistry (BigDye Terminator v3.1 K (Applied Biosystems Inc, Foster City, CA, USA) and 3730 × l ABI 96-capillary sequencer systems equipped with capillaries of 50-cm separation length. In all ASM cases, the presence of the SNP was confirmed by sequencing of unconverted DNA using primers listed in Additional file 1. In cases where ASM was observed, the bisulfite conversion, subcloning and sequencing analysis was performed twice in independent experiments for the same amplicon to exclude clonal PCR amplification. The final results presented here are a combination of the results from both experiments. The separate results are shown in Additional file 3.
We confirmed the methylation difference between alleles by using the methylation-sensitive restriction enzymes HpaII (for amplicons 23_1, 23_2, and 232) and BstUI (for amplicon 262) (New England Biolabs, Ipswich, MA, USA). Genomic DNA (250 ng) was digested at the appropriate temperature (37°C for HpaII, 60°C for BstUI) for 12 h in buffer recommended by the supplier. Afterwards, the PCR products were amplified from the digested genomic DNA using the following thermal cycling parameters: 15 minutes at 95°C, 35 cycles of 45 s at 95°C, 45 s at 50°C and 1 minute at 72°C, and finally 72°C for 5 minutes. The sequences of the PCR products were analyzed by direct sequencing or sequencing of subclones. The PCR products amplified from undigested genomic DNA were used as control. The primer sequences are listed in Additional file 1.
Allele-specific mRNA expression analysis
Total RNA was extracted from leukocytes (QIAamp RNA Blood Mini Kit, Qiagen) and was reverse-transcribed. The cDNA was used as templates for PCR by using the following thermal cycling parameters: 15 minutes at 95°C, 35 cycles of 45 s at 95°C, 35 s at 55°C and 50 s at 72°C and finally 72°C for 5 minutes. The primers were designed to span different exons of genes C21orf81, HSF2BP and RIPK4 and include exonic SNPs (listed in Additional file 1). The RT-PCR product was analyzed by direct sequencing.
Data analysis and presentation
BiQ analyzer software was used to perform quality control and to derive DNA methylation patterns from the sequencing results . The BDPC (Bisulfite sequencing Data Presentation and Compilation) program was used to present the methylation pattern of single amplicons, single alleles and prepare the methylation map of all amplicons studied . All data can be accessed at .
The following additional data are available with the online version of this paper: detailed information on the amplicons analyzed and the correlated genes, genomic positions and primers used for bisulfite genomic sequencing and gene expression analysis (Additional file 1); detailed information for all individuals selected for the analysis (Additional file 2); confirmation of the bisulfite genomic sequencing results (Additional file 3); DNA methylation patterns of amplicons 23_1 and 23_2 in different individuals (Additional file 4); DNA methylation patterns of amplicon 262 in different individuals (Additional file 5); methylation levels of homozygous individuals compared with the allelic methylation levels of heterozygous individuals for amplicons 23_1, 23_2 and 262 (Additional file 6); methylation patterns of amplicon 232 in different individuals (Additional file 7); ASM of amplicons 176_1 and 176_2 (Additional file 8); compilation of the methylation levels of regions with ASM in HEK293 cells, HEPG2 cells and human fibroblasts from  (Additional file 9); allele-specific transcription factor binding sites predicted at the ASM amplicons and mechanism of ASM (Additional file 10).
single nucleotide polymorphism.
This work was supported by the NGF2 program of the Federal minster of research and education BMBF and by the Robert Bosch Foundation (12.5.1366.0005.0). We acknowledge technical assistance by Ms Sandra Becker, and Dr Philipp Rathert for support and discussions.
- Reik W: Stability and flexibility of epigenetic gene regulation in mammalian development. Nature. 2007, 447: 425-432. 10.1038/nature05918.PubMedView ArticleGoogle Scholar
- Klose RJ, Bird AP: Genomic DNA methylation: the mark and its mediators. Trends Biochem Sci. 2006, 31: 89-97. 10.1016/j.tibs.2005.12.008.PubMedView ArticleGoogle Scholar
- Goll MG, Bestor TH: Eukaryotic cytosine methyltransferases. Annu Rev Biochem. 2005, 74: 481-514. 10.1146/annurev.biochem.74.010904.153721.PubMedView ArticleGoogle Scholar
- Reik W, Walter J: Evolution of imprinting mechanisms: the battle of the sexes begins in the zygote. Nat Genet. 2001, 27: 255-256. 10.1038/85804.PubMedView ArticleGoogle Scholar
- Wilkins JF: Genomic imprinting and methylation: epigenetic canalization and conflict. Trends Genet. 2005, 21: 356-365. 10.1016/j.tig.2005.04.005.PubMedView ArticleGoogle Scholar
- Heard E: Recent advances in X-chromosome inactivation. Curr Opin Cell Biol. 2004, 16: 247-255. 10.1016/j.ceb.2004.03.005.PubMedView ArticleGoogle Scholar
- Panning B: X-chromosome inactivation: the molecular basis of silencing. J Biol. 2008, 7: 30-10.1186/jbiol95.PubMedPubMed CentralView ArticleGoogle Scholar
- Kerkel K, Spadola A, Yuan E, Kosek J, Jiang L, Hod E, Li K, Murty VV, Schupf N, Vilain E, Morris M, Haghighi F, Tycko B: Genomic surveys by methylation-sensitive SNP analysis identify sequence-dependent allele-specific DNA methylation. Nat Genet. 2008, 40: 904-908. 10.1038/ng.174.PubMedView ArticleGoogle Scholar
- Boumber YA, Kondo Y, Chen X, Shen L, Guo Y, Tellez C, Estecio MR, Ahmed S, Issa JP: An Sp1/Sp3 binding polymorphism confers methylation protection. PLoS Genet. 2008, 4: e1000162-10.1371/journal.pgen.1000162.PubMedPubMed CentralView ArticleGoogle Scholar
- Fraenkel S, Mostoslavsky R, Novobrantseva TI, Pelanda R, Chaudhuri J, Esposito G, Jung S, Alt FW, Rajewsky K, Cedar H, Bergman Y: Allelic 'choice' governs somatic hypermutation in vivo at the immunoglobulin kappa-chain locus. Nat Immunol. 2007, 8: 715-722. 10.1038/ni1476.PubMedView ArticleGoogle Scholar
- Chan HW, Kurago ZB, Stewart CA, Wilson MJ, Martin MP, Mace BE, Carrington M, Trowsdale J, Lutz CT: DNA methylation maintains allele-specific KIR gene expression in human natural killer cells. J Exp Med. 2003, 197: 245-255. 10.1084/jem.20021127.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Y, Rohde C, Tierling S, Jurkowski TP, Bock C, Santacruz D, Ragozin S, Reinhardt R, Groth M, Walter J, Jeltsch A: DNA methylation analysis of chromosome 21 gene promoters at single base pair and single allele resolution. PLOS Genet. 2009, 5: e1000438-10.1371/journal.pgen.1000438.PubMedPubMed CentralView ArticleGoogle Scholar
- Richards EJ: Inherited epigenetic variation--revisiting soft inheritance. Nat Rev Genet. 2006, 7: 395-401. 10.1038/nrg1834.PubMedView ArticleGoogle Scholar
- Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, Burger M, Burton J, Cox TV, Davies R, Down TA, Haefliger C, Horton R, Howe K, Jackson DK, Kunde J, Koenig C, Liddle J, Niblett D, Otto T, Pettett R, Seemann S, Thompson C, West T, Rogers J, Olek A, Berlin K, Beck S: DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006, 38: 1378-1385. 10.1038/ng1909.PubMedPubMed CentralView ArticleGoogle Scholar
- Gimelbrant A, Hutchinson JN, Thompson BR, Chess A: Widespread monoallelic expression on human autosomes. Science. 2007, 318: 1136-1140. 10.1126/science.1148910.PubMedView ArticleGoogle Scholar
- Serre D, Gurd S, Ge B, Sladek R, Sinnett D, Harmsen E, Bibikova M, Chudin E, Barker DL, Dickinson T, Fan JB, Hudson TJ: Differential allelic expression in the human genome: a robust approach to identify genetic and epigenetic cis-acting mechanisms regulating gene expression. PLoS Genet. 2008, 4: e1000006-10.1371/journal.pgen.1000006.PubMedPubMed CentralView ArticleGoogle Scholar
- Milani L, Lundmark A, Nordlund J, Kiialainen A, Flaegstad T, Jonmundsson G, Kanerva J, Schmiegelow K, Gunderson KL, Lonnerholm G, Syvanen AC: Allele-specific gene expression patterns in primary leukemic cells reveal regulation of gene expression by CpG site methylation. Genome Res. 2009, 19: 1-11. 10.1101/gr.083931.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Palacios R, Gazave E, Goni J, Piedrafita G, Fernando O, Navarro A, Villoslada P: Allele-specific gene expression is widespread across the genome and biological processes. PLoS ONE. 2009, 4: e4150-10.1371/journal.pone.0004150.PubMedPubMed CentralView ArticleGoogle Scholar
- Feinberg AP, Tycko B: The history of cancer epigenetics. Nat Rev Cancer. 2004, 4: 143-153. 10.1038/nrc1279.PubMedView ArticleGoogle Scholar
- Feinberg AP: An epigenetic approach to cancer etiology. Cancer J. 2007, 13: 70-74. 10.1097/PPO.0b013e31803c6e3b.PubMedView ArticleGoogle Scholar
- Ohlsson R, Tycko B, Sapienza C: Monoallelic expression: 'there can only be one'. Trends Genet. 1998, 14: 435-438. 10.1016/S0168-9525(98)01583-2.PubMedView ArticleGoogle Scholar
- Zhang Y, Rohde C, Tierling S, Stamerjohanns H, Reinhardt R, Walter J, Jeltsch A: DNA methylation analysis by bisulfite conversion, cloning, and sequencing of individual clones. Methods Mol Biol. 2009, 507: 177-187. full_text.PubMedView ArticleGoogle Scholar
- Bock C, Reither S, Mikeska T, Paulsen M, Walter J, Lengauer T: BiQ Analyzer: visualization and quality control for DNA methylation data from bisulfite sequencing. Bioinformatics. 2005, 21: 4067-4068. 10.1093/bioinformatics/bti652.PubMedView ArticleGoogle Scholar
- Rohde C, Zhang Y, Jurkowski TP, Stamerjohanns H, Reinhardt R, Jeltsch A: Bisulfite sequencing Data Presentation and Compilation (BDPC) web server - a useful tool for DNA methylation analysis. Nucleic Acids Res. 2008, 36: e34-10.1093/nar/gkn083.PubMedPubMed CentralView ArticleGoogle Scholar
- Database for Non-imprinted Allele-specific DNA Methylation on Human Autosomes. [http://biochem.jacobs-university.de/ASM/]
- Heinemeyer T, Wingender E, Reuter I, Hermjakob H, Kel AE, Kel OV, Ignatieva EV, Ananko EA, Podkolodnaya OA, Kolpakov FA, Podkolodny NL, Kolchanov NA: Databases on transcriptional regulation: TRANSFAC, TRRD and COMPEL. Nucleic Acids Res. 1998, 26: 362-367. 10.1093/nar/26.1.362.PubMedPubMed CentralView ArticleGoogle Scholar
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.