  CpG islands (CGI) marked by bivalent chromatin in stem cells are believed to
be more prone to aberrant DNA methylation in tumor cells. The robustness and
genome-wide extent of this instructive program in different cancer types remain
to be determined. To address this issue we developed a user-friendly approach
to integrate the stem cell chromatin signature in customized DNA methylation
analyses. We used publicly available ChIP-sequencing datasets of several human
embryonic stem cell (hESC) lines to determine the extent of bivalent chromatin
genome-wide. We then created annotated lists of high-confidence bivalent,
H3K4me3-only and H3K27me3-only chromatin regions. The main features of bivalent
regions included localization in CGI/promoters, depletion in retroelements and
enrichment in specific histone modifications, including the poorly
characterized H3K23me2 mark. Moreover, bivalent promoters could be classified
in three clusters based on PRC2 and PolII complexes occupancy. Genes with
bivalent promoters of the PRC2-defined cluster displayed the lowest expression
upon differentiation. As proof-of-concept, we assessed the DNA methylation
pattern of eight types of tumors and confirmed that aberrant cancer-associated
DNA hypermethylation preferentially targets CGI characterized by bivalent
chromatin in hESCs. We also found that such aberrant DNA hypermethylation
affected particularly bivalent CGI/promoters associated with genes that tend to
remain repressed upon differentiation. Strikingly, bivalent CGI were the most
affected by aberrant DNA hypermethylation in both CpG Island Methylator
Phenotype-positive (CIMP+) and CIMP-negative tumors, suggesting that, besides
transcriptional silencing in the pre-tumorigenic cells, the bivalent chromatin
signature in hESCs is a key determinant of the instructive program for aberrant
DNA methylation.
