Background Bivalent genes, of which promoters are marked by both H3K4me3 (trimethylation of histone H3 on lysine 4) and H3K27me3 (trimethylation of histone H3 on lysine 27), play critical roles in development and tumorigenesis. Monomethylation on lysine 4 of histone H3 (H3K4me1) is commonly associated with enhancers, but H3K4me1 is also present at promoter regions as an active bimodal or a repressed unimodal pattern. Whether the co-occurrence of H3K4me1 and bivalent marks at promoters plays regulatory role in development is largely unknown. Results We report that in the process of lineage differentiation, bivalent promoters undergo H3K27me3-H3K4me1 transition, the loss of H3K27me3 accompanies by bimodal pattern loss or unimodal pattern enrichment of H3K4me1. More importantly, this transition regulates tissue-specific gene expression to orchestrate the development. Furthermore, knockout of Eed (Embryonic Ectoderm Development) or Suz12 (Suppressor of Zeste 12) in mESCs (mouse embryonic stem cells), the core components of Polycomb repressive complex 2 (PRC2) which catalyzes H3K27 trimethylation, generates an artificial H3K27me3-H3K4me1 transition at partial bivalent promoters, which leads to up-regulation of meso-endoderm related genes and down-regulation of ectoderm related genes, thus could explain the observed neural ectoderm differentiation failure upon retinoic acid (RA) induction. Finally, we find that lysine-specific demethylase 1 (LSD1) interacts with PRC2 and contributes to the H3K27me3-H3K4me1 transition in mESCs. Conclusions These findings suggest that H3K27me3-H3K4me1 transition plays a key role in lineage differentiation by regulating the expression of tissue specific genes, and H3K4me1 pattern in bivalent promoters could be modulated by LSD1 via interacting with PRC2.

N6-methyladenosine (m6A) is the most prevalent internal modification and reversible epitranscriptomic mark in messenger RNAs (mRNAs) and plays essential roles in a variety of biological processes. However, the dynamic distribution patterns of m6A and their significance during mammalian tissue development are poorly understood. Here, we found that based on m6A distribution patterns, protein-coding genes were classified into five groups with significantly distinct biological features and functions. Strikingly, comparison of the m6A methylomes of multiple mammalian tissues between fetal and adult stages revealed dynamic m6A topological transition during mammalian tissue development, and identified large numbers of genes with significant m6A loss in 5′UTRs or m6A gain around stop codons. The genes with m6A loss in 5′UTRs were highly enriched in developmental stage-specific genes, and their m6A topological transitions were strongly associated with gene expression regulation during tissue development. The genes with m6A gain around the stop codons were associated with tissue-specific functions. Our findings revealed the existence of different m6A topologies among protein-coding genes that were associated with distinct characteristics. More importantly, these genes with m6A topological transitions were crucial for tissue development via regulation of gene expression, suggesting the importance of dynamic m6A topological transitions during mammalian tissue development.