Biallelic non-productive enhancer-promoter interactions precede imprinted expression of Kcnk9 during mouse neural commitment.

It is only partially understood how constitutive allelic methylation at imprinting control regions (ICRs) interacts with other regulation levels to drive timely parental allele-specific expression along large imprinted domains. The Peg13-Kcnk9 domain is an imprinted domain with important brain functions. To gain insights into its regulation during neural commitment, we performed an integrative analysis of its allele-specific epigenetic, transcriptomic, and cis-spatial organization using a mouse stem cell-based corticogenesis model that recapitulates the control of imprinted gene expression during neurodevelopment. We found that, despite an allelic higher-order chromatin structure associated with the paternally CTCF-bound Peg13 ICR, enhancer-Kcnk9 promoter contacts occurred on both alleles, although they were productive only on the maternal allele. This observation challenges the canonical model in which CTCF binding isolates the enhancer and its target gene on either side and suggests a more nuanced role for allelic CTCF binding at some ICRs.

Influence of pH on indole-dependent heterodimeric interactions between Anopheles gambiae odorant-binding proteins OBP1 and OBP4.

Insect Odorant Binding Proteins (OBPs) constitute important components of their olfactory apparatus, as they are essential for odor recognition. OBPs undergo conformational changes upon pH change, altering their interactions with odorants. Moreover, they can form heterodimers with novel binding characteristics. Anopheles gambiae OBP1 and OBP4 were found capable of forming heterodimers possibly involved in the specific perception of the attractant indole. In order to understand how these OBPs interact in the presence of indole and to investigate the likelihood of a pH-dependent heterodimerization mechanism, the crystal structures of OBP4 at pH 4.6 and 8.5 were determined. Structural comparison to each other and with the OBP4-indole complex (3Q8I, pH 6.85) revealed a flexible N-terminus and conformational changes in the α4-loop-α5 region at acidic pH. Fluorescence competition assays showed a weak binding of indole to OBP4 that becomes further impaired at acidic pH. Additional Molecular Dynamic and Differential Scanning Calorimetry studies displayed that the influence of pH on OBP4 stability is significant compared to the modest effect of indole. Furthermore, OBP1-OBP4 heterodimeric models were generated at pH 4.5, 6.5, and 8.5, and compared concerning their interface energy and cross-correlated motions in the absence and presence of indole. The results indicate that the increase in pH may induce the stabilization of OBP4 by increasing its helicity, thereby enabling indole binding at neutral pH that further stabilizes the protein and possibly promotes the creation of a binding site for OBP1. A decrease in interface stability and loss of correlated motions upon transition to acidic pH may provoke the heterodimeric dissociation allowing indole release. Finally, we propose a potential OBP1-OBP4 heterodimer formation/disruption mechanism induced by pH change and indole binding.

Variable allelic expression of imprinted genes at the Peg13, Trappc9, Ago2 cluster in single neural cells.

Genomic imprinting is an epigenetic process through which genes are expressed in a parent-of-origin specific manner resulting in mono-allelic or strongly biased expression of one allele. For some genes, imprinted expression may be tissue-specific and reliant on CTCF-influenced enhancer-promoter interactions. The Peg13 imprinting cluster is associated with neurodevelopmental disorders and comprises canonical imprinted genes, which are conserved between mouse and human, as well as brain-specific imprinted genes in mouse. The latter consist of Trappc9, Chrac1 and Ago2, which have a maternal allelic expression bias of ∼75% in brain. Findings of such allelic expression biases on the tissue level raise the question of how they are reflected in individual cells and whether there is variability and mosaicism in allelic expression between individual cells of the tissue. Here we show that Trappc9 and Ago2 are not imprinted in hippocampus-derived neural stem cells (neurospheres), while Peg13 retains its strong bias of paternal allele expression. Upon analysis of single neural stem cells and in vitro differentiated neurons, we find not uniform, but variable states of allelic expression, especially for Trappc9 and Ago2. These ranged from mono-allelic paternal to equal bi-allelic to mono-allelic maternal, including biased bi-allelic transcriptional states. Even Peg13 expression deviated from its expected paternal allele bias in a small number of cells. Although the cell populations consisted of a mosaic of cells with different allelic expression states, as a whole they reflected bulk tissue data. Furthermore, in an attempt to identify potential brain-specific regulatory elements across the Trappc9 locus, we demonstrate tissue-specific and general silencer activities, which might contribute to the regulation of its imprinted expression bias.

The Long Non-Coding RNA HOXA-AS2 Promotes Proliferation of Glioma Stem Cells and Modulates Their Inflammation Pathway Mainly through Post-Transcriptional Regulation.

Glioblastomas represent approximatively half of all gliomas and are the most deadly and aggressive form. Their therapeutic resistance and tumor relapse rely on a subpopulation of cells that are called Glioma Stem Cells (GSCs). Here, we investigated the role of the long non-coding RNA HOXA-AS2 in GSC biology using descriptive and functional analyses of glioma samples classified according to their isocitrate dehydrogenase (IDH) gene mutation status, and of GSC lines. We found that HOXA-AS2 is overexpressed only in aggressive (IDHwt) glioma and GSC lines. ShRNA-based depletion of HOXA-AS2 in GSCs decreased cell proliferation and altered the expression of several hundreds of genes. Integrative analysis revealed that these expression changes were not associated with changes in DNA methylation or chromatin signatures at the promoter of the majority of genes deregulated following HOXA-AS2 silencing in GSCs, suggesting a post-transcriptional regulation. In addition, transcription factor binding motif enrichment and correlation analyses indicated that HOXA-AS2 affects, directly or indirectly, the expression of key transcription factors implicated in GCS biology, including E2F8, E2F1, STAT1, and ATF3, thus contributing to GCS aggressiveness by promoting their proliferation and modulating the inflammation pathway.

L1 chimeric transcripts are expressed in healthy brain and their deregulation in glioma follows that of their host locus.

Besides the consequences of retrotransposition, long interspersed element 1 (L1) retrotransposons can affect the host genome through their antisense promoter. In addition to the sense promoter, the evolutionarily recent L1 retrotransposons, which are present in several thousand copies, also possess an anti-sense promoter that can produce L1 chimeric transcripts (LCT) composed of the L1 5' UTR followed by the adjacent genomic sequence. The full extent to which LCT expression occurs in a given tissue and whether disruption of the defense mechanisms that normally control L1 retrotransposons affects their expression and function in cancer cells, remain to be established. By using CLIFinder, a dedicated bioinformatics tool, we found that LCT expression was widespread in normal brain and aggressive glioma samples, and that approximately 17% of recent L1 retrotransposons, from the L1PA1 to L1PA7 subfamilies, were involved in their production. Importantly, the transcriptional activities of the L1 antisense promoters and of their host loci were coupled. Accordingly, we detected LCT-producing L1 retrotransposons mainly in transcriptionally active genes and genomic loci. Moreover, changes in the host genomic locus expression level in glioma were associated with a similar change in LCT expression level, regardless of the L1 promoter methylation status. Our findings support a model in which the host genomic locus transcriptional activity is the main driving force of LCT expression. We hypothesize that this model is more applicable when host gene and LCT are transcribed from the same strand.

An Expanded Survey of the Moth PBP/GOBP Clade in Bombyx mori: New Insight into Expression and Functional Roles.

We studied the expression profile and ontogeny (from the egg stage through the larval stages and pupal stages, to the elderly adult age) of four OBPs from the silkworm moth Bombyx mori. We first showed that male responsiveness to female sex pheromone in the silkworm moth B. mori does not depend on age variation; whereas the expression of BmorPBP1, BmorPBP2, BmorGOBP1, and BmorGOBP2 varies with age. The expression profile analysis revealed that the studied OBPs are expressed in non-olfactory tissues at different developmental stages. In addition, we tested the effect of insecticide exposure on the expression of the four OBPs studied. Exposure to a toxic macrolide insecticide endectocide molecule (abamectin) led to the modulated expression of all four genes in different tissues. The higher expression of OBPs was detected in metabolic tissues, such as the thorax, gut, and fat body. All these data strongly suggest some alternative functions for these proteins other than olfaction. Finally, we carried out ligand docking studies and reported that PBP1 and GOBP2 have the capacity of binding vitamin K1 and multiple different vitamins.

Widespread overexpression from the four DNA hypermethylated HOX clusters in aggressive (IDHwt) glioma is associated with H3K27me3 depletion and alternative promoter usage.

In human, the 39 coding HOX genes and 18 referenced noncoding antisense transcripts are arranged in four genomic clusters named HOXA, B, C, and D. This highly conserved family belongs to the homeobox class of genes that encode transcription factors required for normal development. Therefore, HOX gene deregulation might contribute to the development of many cancer types. Here, we study HOX gene deregulation in adult glioma, a common type of primary brain tumor. We performed extensive molecular analysis of tumor samples, classified according to their isocitrate dehydrogenase (IDH1) gene mutation status, and of glioma stem cells. We found widespread expression of sense and antisense HOX transcripts only in aggressive (IDHwt) glioma samples, although the four HOX clusters displayed DNA hypermethylation. Integrative analysis of expression, DNA methylation, and histone modification signatures along the clusters revealed that HOX gene upregulation relies on canonical and alternative bivalent CpG island promoters that escape hypermethylation. H3K27me3 loss at these promoters emerges as the main cause of widespread HOX gene upregulation in IDHwt glioma cell lines and tumors. Our study provides the first comprehensive description of the epigenetic changes at HOX clusters and their contribution to the transcriptional changes observed in adult glioma. It also identified putative 'master' HOX proteins that might contribute to the tumorigenic potential of glioma stem cells.

TET3 controls the expression of the H3K27me3 demethylase Kdm6b during neural commitment.

The acquisition of cell identity is associated with developmentally regulated changes in the cellular histone methylation signatures. For instance, commitment to neural differentiation relies on the tightly controlled gain or loss of H3K27me3, a hallmark of polycomb-mediated transcriptional gene silencing, at specific gene sets. The KDM6B demethylase, which removes H3K27me3 marks at defined promoters and enhancers, is a key factor in neurogenesis. Therefore, to better understand the epigenetic regulation of neural fate acquisition, it is important to determine how Kdm6b expression is regulated. Here, we investigated the molecular mechanisms involved in the induction of Kdm6b expression upon neural commitment of mouse embryonic stem cells. We found that the increase in Kdm6b expression is linked to a rearrangement between two 3D configurations defined by the promoter contact with two different regions in the Kdm6b locus. This is associated with changes in 5-hydroxymethylcytosine (5hmC) levels at these two regions, and requires a functional ten-eleven-translocation (TET) 3 protein. Altogether, our data support a model whereby Kdm6b induction upon neural commitment relies on an intronic enhancer the activity of which is defined by its TET3-mediated 5-hmC level. This original observation reveals an unexpected interplay between the 5-hmC and H3K27me3 pathways during neural lineage commitment in mammals. It also questions to which extent KDM6B-mediated changes in H3K27me3 level account for the TET-mediated effects on gene expression.

Residual Solvents in Nanomedicine and Lipid-Based Drug Delivery Systems: a Case Study to Better Understand Processes.

PURPOSE: Complexities surrounding the manufacture and quality control of nanomedicines become increasingly apparent. This research article offers a case study to investigate how, at the laboratory scale, various stages of liposome and nanoparticle synthesis affect the amount of residual solvent found in the formulations. The objective is to bring insights on the reliability of each of these processes to provide final products which meet regulatory standards and facilitate identifying possible bottleneck early during the development process. METHODS: The residual solvent at various stages of preparation and purification was measured by headspace gas chromatography. Liposomes were prepared by two different methods with and without solvent. Polymer nanoparticles prepared via nanoprecipitation and purified by ultrafiltration were studied. The effects of purification by size exclusion chromatography and dialysis were also investigated. RESULTS: The complete removal of residual solvent requires processes which go beyond usual preparation methods. CONCLUSIONS: This work might prove valuable as a reference for scientists of different fields to compare their own practices and streamline the translation of nanomedicines into efficacious and safe drug products.

[Argatroban in heparin-induced thrombocytopenia: A survey of prescribing practices and use in France]. / Argatroban et thrombopénie induite par l'héparine : observatoire des modalités de prescription et d'utilisation en France.

OBJECTIVE: A survey of prescribing practices was carried out in France to ensure that argatroban was used appropriately during the first 18months after it obtained marketing authorization for anticoagulation in adults with heparin-induced thrombocytopenia (HIT). METHODS: This observational study was proposed to public and private hospitals with at least 2 orders of argatroban. All patients who received at least one argatroban injection during the study period had to be included. Their demographic characteristics, the pathology causing heparin treatment, the indication of treatment with argatroban, as well as available real-life clinical and biological monitoring data were retrospectively collected. RESULTS: In the 23 participating centers, the drug was prescribed mainly by the following hospital units: surgery and intensive care (79.3%) (of which 18.3% cardiovascular), nephrology/hemodialysis (14.8%) and internal medicine (5.9%). Among the 169 patients included, with median age of 68 years, 118 (69.8%) had renal impairment (creatinine clearance <90mL/min). The HIT probability scoring and diagnostic tests used differed widely between the centers. The reasons for prescribing the drug were mainly suspected HIT (51.5%), declared confirmed HIT (21.3%) or a history of HIT (14.8%). Seventy-three (73) patients (43.2%) had thrombotic symptoms or a systemic reaction initially. The median initial dose prescribed was 0.5µg/kg/min (ranging from 0.05 to 4.42) and doses >2µg/kg/min were used during hemodialysis. More than half of the patients (58.6%) had no clinical complications. Most of the serious adverse reactions were hemorrhagic (11/12). CONCLUSION: This study illustrates the complexity of treatment for HIT and the need to be familiar with and follow guidelines on the management of HIT, especially for susceptible patients treated in intensive care units.

Comprehensive History of CSP Genes: Evolution, Phylogenetic Distribution and Functions.

In this review we present the developmental, histological, evolutionary and functional properties of insect chemosensory proteins (CSPs) in insect species. CSPs are small globular proteins folded like a prism and notoriously known for their complex and arguably obscure function(s), particularly in pheromone olfaction. Here, we focus on direct functional consequences on protein function depending on duplication, expression and RNA editing. The result of our analysis is important for understanding the significance of RNA-editing on functionality of CSP genes, particularly in the brain tissue.

HOX gene cluster (de)regulation in brain: from neurodevelopment to malignant glial tumours.

HOX genes encode a family of evolutionarily conserved homeodomain transcription factors that are crucial both during development and adult life. In humans, 39 HOX genes are arranged in four clusters (HOXA, B, C, and D) in chromosomes 7, 17, 12, and 2, respectively. During embryonic development, particular epigenetic states accompany their expression along the anterior-posterior body axis. This tightly regulated temporal-spatial expression pattern reflects their relative chromosomal localization, and is critical for normal embryonic brain development when HOX genes are mainly expressed in the hindbrain and mostly absent in the forebrain region. Epigenetic marks, mostly polycomb-associated, are dynamically regulated at HOX loci and regulatory regions to ensure the finely tuned HOX activation and repression, highlighting a crucial epigenetic plasticity necessary for homeostatic development. HOX genes are essentially absent in healthy adult brain, whereas they are detected in malignant brain tumours, namely gliomas, where HOX genes display critical roles by regulating several hallmarks of cancer. Here, we review the major mechanisms involved in HOX genes (de)regulation in the brain, from embryonic to adult stages, in physiological and oncologic conditions. We focus particularly on the emerging causes of HOX gene deregulation in glioma, as well as on their functional and clinical implications.

Toward the design of efficient transglycosidases: the case of the GH1 of Thermus thermophilus.

Using the information available in the sequences of well-characterized transglycosidases found in plants, mutations were introduced in the glycoside hydrolase of the bacterium Thermus thermophilus, with the aim of turning it into an efficient transglycosidase. All mutants happen to have fair catalytic efficiencies, being at worst 25 times less efficient than the wild type. Noteworthy, W120F, one of our high transglycosylation yield (≈ 50%) mutants, is only two times less efficient than the wild type. Interestingly, while in the wild type the sidechain of the acid-base is only found able to sample a pair of equivalent conformations during 0.5-µs-long molecular dynamics simulations, its flexibility is much higher in the case of the high transglycosylation yield mutants. Our results thus suggest that engineering the flexibility of the acid-base of a retaining glycoside hydrolase could be a general way to turn it into an efficient transglycosidase.

Transcriptional alterations in glioma result primarily from DNA methylation-independent mechanisms.

In cancer cells, aberrant DNA methylation is commonly associated with transcriptional alterations, including silencing of tumor suppressor genes. However, multiple epigenetic mechanisms, including polycomb repressive marks, contribute to gene deregulation in cancer. To dissect the relative contribution of DNA methylation-dependent and -independent mechanisms to transcriptional alterations at CpG island/promoter-associated genes in cancer, we studied 70 samples of adult glioma, a widespread type of brain tumor, classified according to their isocitrate dehydrogenase (IDH1) mutation status. We found that most transcriptional alterations in tumor samples were DNA methylation-independent. Instead, altered histone H3 trimethylation at lysine 27 (H3K27me3) was the predominant molecular defect at deregulated genes. Our results also suggest that the presence of a bivalent chromatin signature at CpG island promoters in stem cells predisposes not only to hypermethylation, as widely documented, but more generally to all types of transcriptional alterations in transformed cells. In addition, the gene expression strength in healthy brain cells influences the choice between DNA methylation- and H3K27me3-associated silencing in glioma. Highly expressed genes were more likely to be repressed by H3K27me3 than by DNA methylation. Our findings support a model in which altered H3K27me3 dynamics, more specifically defects in the interplay between polycomb protein complexes and the brain-specific transcriptional machinery, is the main cause of transcriptional alteration in glioma cells. Our study provides the first comprehensive description of epigenetic changes in glioma and their relative contribution to transcriptional changes. It may be useful for the design of drugs targeting cancer-related epigenetic defects.

The long non-coding RNA HOTAIR is transcriptionally activated by HOXA9 and is an independent prognostic marker in patients with malignant glioma.

The lncRNA HOTAIR has been implicated in several human cancers. Here, we evaluated the molecular alterations and upstream regulatory mechanisms of HOTAIR in glioma, the most common primary brain tumors, and its clinical relevance. HOTAIR gene expression, methylation, copy-number and prognostic value were investigated in human gliomas integrating data from online datasets and our cohorts. High levels of HOTAIR were associated with higher grades of glioma, particularly IDH wild-type cases. Mechanistically, HOTAIR was overexpressed in a gene dosage-independent manner, while DNA methylation levels of particular CpGs in HOTAIR locus were associated with HOTAIR expression levels in GBM clinical specimens and cell lines. Concordantly, the demethylating agent 5-Aza-2'-deoxycytidine affected HOTAIR transcriptional levels in a cell line-dependent manner. Importantly, HOTAIR was frequently co-expressed with HOXA9 in high-grade gliomas from TCGA, Oncomine, and our Portuguese and French datasets. Integrated in silico analyses, chromatin immunoprecipitation, and qPCR data showed that HOXA9 binds directly to the promoter of HOTAIR. Clinically, GBM patients with high HOTAIR expression had a significantly reduced overall survival, independently of other prognostic variables. In summary, this work reveals HOXA9 as a novel direct regulator of HOTAIR, and establishes HOTAIR as an independent prognostic marker, providing new therapeutic opportunities to treat this highly aggressive cancer.

DNA methylation profiling reveals a pathological signature that contributes to transcriptional defects of CD34+ CD15- cells in early chronic-phase chronic myeloid leukemia.

Despite the high efficiency of tyrosine kinase inhibitors (TKI), some patients with chronic myeloid leukemia (CML) will display residual disease that can become resistant to treatment, indicating intraclonal heterogeneity in chronic-phase CML (CP-CML). To determine the basis of this heterogeneity, we conducted the first exhaustive characterization of the DNA methylation pattern of sorted CP-CML CD34+ CD15- (immature) and CD34- CD15+ (mature) cells at diagnosis (prior to any treatment) and compared it to that of CD34+ CD15- and CD34- CD15+ cells isolated from healthy donors (HD). In both cell types, we identified several hundreds of differentially methylated regions (DMRs) showing DNA methylation changes between CP-CML and HD samples, with only a subset of them in common between CD34+ CD15- and CD34- CD15+ cells. This suggested DNA methylation variability within the same CML clone. We also identified 70 genes that could be aberrantly repressed upon hypermethylation and 171 genes that could be aberrantly expressed upon hypomethylation of some of these DMRs in CP-CML cells, among which 18 and 81, respectively, were in CP-CML CD34+ CD15- cells only. We then validated the DNA methylation and expression defects of selected candidate genes. Specifically, we identified GAS2, a candidate oncogene, as a new example of gene the hypomethylation of which is associated with robust overexpression in CP-CML cells. Altogether, we demonstrated that DNA methylation abnormalities exist at early stages of CML and can affect the transcriptional landscape of malignant cells. These observations could lead to the development of combination treatments with epigenetic drugs and TKI for CP-CML.

CLIFinder: identification of LINE-1 chimeric transcripts in RNA-seq data.

Summary: L1 Chimeric Transcripts (LCTs) are initiated by repeated LINE-1 element antisense promoters and include the L1 5'UTR sequence in antisense orientation followed by the adjacent genomic region. LCTs have been characterized mainly using bioinformatics approaches to query dbEST. To take advantage of NGS data to unravel the transcriptome composition, we developed Chimeric LIne Finder (CLIFinder), a new bioinformatics tool. Using stranded paired-end RNA-seq data, we demonstrated that CLIFinder can identify genome-wide transcribed chimera sequences corresponding to potential LCTs. Moreover, CLIFinder can be adapted to study transcription from other repeat types. Availability and implementation: The code is available at: https://github.com/GReD-Clermont/CLIFinder; and for Galaxy users, it is directly accessible in the tool shed at: https://toolshed.g2.bx.psu.edu/view/clifinder/clifinder/. Contact: catherine.barriere@uca.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Genotyping and Bio-Sensing Chemosensory Proteins in Insects.

Genotyping is the process of determining differences in the genetic make-up of an individual and comparing it to that of another individual. Focus on the family of chemosensory proteins (CSPs) in insects reveals differences at the genomic level across various strains and biotypes, but none at the level of individuals, which could be extremely useful in the biotyping of insect pest species necessary for the agricultural, medical and veterinary industries. Proposed methods of genotyping CSPs include not only restriction enzymatic cleavage and amplification of cleaved polymorphic sequences, but also detection of retroposons in some specific regions of the insect chromosome. Design of biosensors using CSPs addresses tissue-specific RNA mutations in a particular subtype of the protein, which could be used as a marker of specific physiological conditions. Additionally, we refer to the binding properties of CSP proteins tuned to lipids and xenobiotic insecticides for the development of a new generation of biosensor chips, monitoring lipid blood concentration and chemical environmental pollution.

Detection of the alternative lengthening of telomeres pathway in malignant gliomas for improved molecular diagnosis.

Human malignant gliomas exhibit acquisition of either one of two telomere maintenance mechanisms, resulting from either reactivation of telomerase expression or activation of an alternative lengthening of telomeres (ALT) mechanism. In the present study, we analyzed 63 human malignant gliomas for the presence of ALT-specific extrachromosomal circles of telomeric DNA (C-circles) and measured telomerase expression, telomeric DNA content (Telo/Alu method), and telomeric repeat-containing RNAs (TERRA) levels. We also assessed histomolecular markers routinely used in clinical practice. The presence of C-circles significantly correlated with IDH1/2 mutation, MGMT exon 1 methylation, low Ki-67 immunostaining, increased telomeric DNA content, absence of functional ATRX protein and level of HTERT gene expression. In multivariate analysis, we observed a trend to a correlation between elevated TERRA levels and increased survival. Interestingly, the C-circles assay allowed to detect ALT activation in glioblastomas exhibiting wild-type IDH1/2 and ATRX expression. These results suggest that, after the correlations uncovered here have been confirmed on larger numbers of tumors, telomeric markers might be useful in improving diagnosis. They also point out to the utility of using the specific, sensitive and quantitative C-circle and Telo/Alu assays that can work with as few as 30 ng of tumor DNA.

An annotated list of bivalent chromatin regions in human ES cells: a new tool for cancer epigenetic research.

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.