A human B cell methylome at 100-base pair resolution - PubMed (original) (raw)
A human B cell methylome at 100-base pair resolution
Tibor A Rauch et al. Proc Natl Acad Sci U S A. 2009.
Abstract
Using a methylated-DNA enrichment technique (methylated CpG island recovery assay, MIRA) in combination with whole-genome tiling arrays, we have characterized by MIRA-chip the entire B cell "methylome" of an individual human at 100-bp resolution. We find that at the chromosome level high CpG methylation density is correlated with subtelomeric regions and Giemsa-light bands (R bands). The majority of the most highly methylated regions that could be identified on the tiling arrays were associated with genes. Approximately 10% of all promoters in B cells were found to be methylated, and this methylation correlates with low gene expression. Notably, apparent exceptions to this correlation were the result of transcription from previously unidentified, unmethylated transcription start sites, suggesting that methylation may control alternate promoter usage. Methylation of intragenic (gene body) sequences was found to correlate with increased, not decreased, transcription, and a methylated region near the 3' end was found in approximately 12% of all genes. The majority of broad regions (10-44 kb) of high methylation were at segmental duplications. Our data provide a valuable resource for the analysis of CpG methylation patterns in a differentiated human cell type and provide new clues regarding the function of mammalian DNA methylation.
Conflict of interest statement
The authors declare no conflict of interest.
Figures
Fig. 1.
Snapshots of MIRA data from whole-genome methylation analysis. (A) The HOXA gene cluster on chromosome 7 was analyzed in detail. The methylation signal is shown plotted along the chromosome as a −log10 P value score (green). Therefore, the minimum number on the y axis is 0 (when P = 1). The P value score was obtained by NimbleScan software and is derived from the Kolmogorov-Smirnov test comparing the log2 ratios (MIRA vs. input) within a 750-bp window centered at each probe and the rest of the data on the array. The HOXA genes and the direction of transcription are indicated. The red boxes show the CpG islands. Methylation patterns at several CpG islands were confirmed by bisulfite sequencing. Open circles represent unmethylated CpG dinucleotides, and black circles represent methylated CpGs. The MIRA signals correlate with the density of methylated CpGs. (B) Analysis of the XIST locus on the active X chromosome. The promoter of the XIST gene shows a strong MRI.
Fig. 2.
Distribution of DNA methylation along the B cell genome. (A) Distribution of MRIs and methylated CpG islands relative to annotated genes. A total of 42,944 MRIs were identified, and their location with respect to genes and the 5′ and 3′ ends of genes, defined as ±500 bp relative to the transcript start or end, was determined. Intragenic MRIs are located from +500 of transcription start to −500 of transcript end. Methylated CpG islands were defined as those CpG islands that overlap with an MRI. (B) Distribution of MRIs along individual chromosomes. The number of MRIs was plotted against the length of the chromosomes.
Fig. 3.
Methylation patterns along human chromosomes. (A) Chromosomal methylation density map (see text for details). Chromosome 7 is shown as an example. All other chromosomes are displayed in
Figs. S4-S7
. (B) Chromosomal methylation density profile of chromosome 7 (Middle) aligned with gene density profile (Top) and density of SINE elements (Bottom). (C) Correlation of methylation density and Giemsa staining. Giemsa staining intensity (from UCSC database) was partitioned into 5 groups (Giemsa negative and 4 groups, increasing from the 25th to the 100th percentile). Giemsa light bands have higher levels of DNA methylation. (D) Correlation of DNA methylation density with SINE and gene density. Density maps were determined using a 5-Mb sliding window with 500-kb step size along each chromosome. The Spearman correlation between DNA methylation density and gene density or SINE density was calculated, and the distribution of the correlations is plotted.
Fig. 4.
Expression of methylated genes is often explained by alternative promoter usage. In this example, data for the PARP12 gene are shown. The methylation profile is in green, and the red box indicates a CpG island. Initial analysis showed that this gene contains a methylated region in the promoter, but the gene is expressed according to expression microarrays. Using 5′ RACE, we identified an alternative transcript that emanates from an unmethylated promoter. RT-PCR analysis of the 2 alternative transcripts is shown in the bottom gel panel. mRNA1 is initiated from promoter 1, whereas mRNA2 is initiated from promoter 2.
Fig. 5.
Intragenic methylation and methylation at the 3′ end of genes. (A) Correlation between intragenic CpG methylation and gene expression levels. Data for gene expression levels in B cells were obtained from the GNF and GEO databases and were partitioned into 10 groups according to expression level (see Materials and Methods). The y axis is the percentage of methylated genes within each group defined as having at least 1 MRI within the gene body. With increased gene expression levels, there is a tendency toward increased intragenic methylation. Only the most highly expressed genes differ from the trend. (B) Identification of an MRI at the 3′ end of genes. The composite profile of methylation densities along all human RefSeq genes shows MRIs at the 5′ end of genes and near the 3′ end. The y axis is the average log2 ratio of all of the genes at each corresponding genome coordinate window. TSS, transcription start site; TES, transcript end site. (C) Methylation at 3′ gene ends and flanking genes. Genes present in tail-to-head orientation in the configuration indicated were analyzed in more detail. Genes in tail-to-head orientation were 1.5-fold more likely to have 3′-end MRIs than singular genes (P < 0.00001).
Similar articles
- Genome-wide profiling of DNA methylation reveals a class of normally methylated CpG island promoters.
Shen L, Kondo Y, Guo Y, Zhang J, Zhang L, Ahmed S, Shu J, Chen X, Waterland RA, Issa JP. Shen L, et al. PLoS Genet. 2007 Oct;3(10):2023-36. doi: 10.1371/journal.pgen.0030181. Epub 2007 Sep 10. PLoS Genet. 2007. PMID: 17967063 Free PMC article. - DNA methylation analysis of chromosome 21 gene promoters at single base pair and single allele resolution.
Zhang Y, Rohde C, Tierling S, Jurkowski TP, Bock C, Santacruz D, Ragozin S, Reinhardt R, Groth M, Walter J, Jeltsch A. Zhang Y, et al. PLoS Genet. 2009 Mar;5(3):e1000438. doi: 10.1371/journal.pgen.1000438. Epub 2009 Mar 27. PLoS Genet. 2009. PMID: 19325872 Free PMC article. - Cell type-specific DNA methylation at intragenic CpG islands in the immune system.
Deaton AM, Webb S, Kerr AR, Illingworth RS, Guy J, Andrews R, Bird A. Deaton AM, et al. Genome Res. 2011 Jul;21(7):1074-86. doi: 10.1101/gr.118703.110. Epub 2011 May 31. Genome Res. 2011. PMID: 21628449 Free PMC article. - DNA methylation profiling using the methylated-CpG island recovery assay (MIRA).
Rauch TA, Pfeifer GP. Rauch TA, et al. Methods. 2010 Nov;52(3):213-7. doi: 10.1016/j.ymeth.2010.03.004. Epub 2010 Mar 19. Methods. 2010. PMID: 20304072 Free PMC article. Review. - DNA methylation patterns in lung carcinomas.
Pfeifer GP, Rauch TA. Pfeifer GP, et al. Semin Cancer Biol. 2009 Jun;19(3):181-7. doi: 10.1016/j.semcancer.2009.02.008. Epub 2009 Feb 20. Semin Cancer Biol. 2009. PMID: 19429482 Free PMC article. Review.
Cited by
- Triphenyl Phosphate Alters Methyltransferase Expression and Induces Genome-Wide Aberrant DNA Methylation in Zebrafish Larvae.
Negi CK, Bláhová L, Phan A, Bajard L, Blaha L. Negi CK, et al. Chem Res Toxicol. 2024 Sep 16;37(9):1549-1561. doi: 10.1021/acs.chemrestox.4c00223. Epub 2024 Aug 29. Chem Res Toxicol. 2024. PMID: 39205618 Free PMC article. - Stress-induced epigenetic effects driven by maternal lactation in dairy cattle: a comethylation network approach.
López-Catalina A, Reverter A, Alexandre PA, Nguyen LT, González-Recio O. López-Catalina A, et al. Epigenetics. 2024 Dec;19(1):2381856. doi: 10.1080/15592294.2024.2381856. Epub 2024 Jul 23. Epigenetics. 2024. PMID: 39044410 Free PMC article. - Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.
Nair VD, Pincas H, Smith GR, Zaslavsky E, Ge Y, Amper MAS, Vasoya M, Chikina M, Sun Y, Raja AN, Mao W, Gay NR, Esser KA, Smith KS, Zhao B, Wiel L, Singh A, Lindholm ME, Amar D, Montgomery S, Snyder MP, Walsh MJ, Sealfon SC; MoTrPAC Study Group. Nair VD, et al. Cell Genom. 2024 Jun 12;4(6):100421. doi: 10.1016/j.xgen.2023.100421. Epub 2024 May 1. Cell Genom. 2024. PMID: 38697122 Free PMC article. - Establishment of paternal methylation imprint at the H19/Igf2 imprinting control region.
Liao J, Song S, Gusscott S, Fu Z, VanderKolk I, Busscher BM, Lau KH, Brind'Amour J, Szabó PE. Liao J, et al. Sci Adv. 2023 Sep 8;9(36):eadi2050. doi: 10.1126/sciadv.adi2050. Epub 2023 Sep 6. Sci Adv. 2023. PMID: 37672574 Free PMC article.
References
- Bird AP. CpG-rich islands and the function of DNA methylation. Nature. 1986;321:209–213. - PubMed
- Suzuki MM, Bird A. DNA methylation landscapes: Provocative insights from epigenomics. Nat Rev Genet. 2008;9:465–476. - PubMed
- Jones PA. The DNA methylation paradox. Trends Genet. 1999;15:34–37. - PubMed
- Miranda TB, Jones PA. DNA methylation: the nuts and bolts of repression. J Cell Physiol. 2007;213:384–390. - PubMed
- Yoder JA, Walsh CP, Bestor TH. Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 1997;13:335–340. - PubMed
Publication types
MeSH terms
Substances
LinkOut - more resources
Full Text Sources