Synthetic CpG islands reveal DNA sequence determinants of chromatin structure

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically shows phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption-the root of the pathogenesis-is similar in the different disease-relevant cell types. This is possible in principle, because all these cell types are subject to effects from the same causative gene, which has the same kind of function (e.g., methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2 and find that the chromatin accessibility changes in neurons are mostly distinct from changes in B or T cells. This is not because the neuronal accessibility changes occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that regulatory elements disrupted in B/T cells do show chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators and suggest that blood-derived episignatures, although useful diagnostically, may not be well suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Crosstalk within and beyond the Polycomb repressive system

The development of multicellular organisms depends on spatiotemporally controlled differentiation of numerous cell types and their maintenance. To generate such diversity based on the invariant genetic information stored in DNA, epigenetic mechanisms, which are heritable changes in gene function that do not involve alterations to the underlying DNA sequence, are required to establish and maintain unique gene expression programs. Polycomb repressive complexes represent a paradigm of epigenetic regulation of developmentally regulated genes, and the roles of these complexes as well as the epigenetic marks they deposit, namely H3K27me3 and H2AK119ub, have been extensively studied. However, an emerging theme from recent studies is that not only the autonomous functions of the Polycomb repressive system, but also crosstalks of Polycomb with other epigenetic modifications, are important for gene regulation. In this review, we summarize how these crosstalk mechanisms have improved our understanding of Polycomb biology and how such knowledge could help with the design of cancer treatments that target the dysregulated epigenome.

Polycomb repressive complex 2 and its core component EZH2: potential targeted therapeutic strategies for head and neck squamous cell carcinoma

The polycomb group (PcG) comprises a set of proteins that exert epigenetic regulatory effects and play crucial roles in diverse biological processes, ranging from pluripotency and development to carcinogenesis. Among these proteins, enhancer of zeste homolog 2 (EZH2) stands out as a catalytic component of polycomb repressive complex 2 (PRC2), which plays a role in regulating the expression of homologous (Hox) genes and initial stages of x chromosome inactivation. In numerous human cancers, including head and neck squamous cell carcinoma (HNSCC), EZH2 is frequently overexpressed or activated and has been identified as a negative prognostic factor. Notably, EZH2 emerges as a significant gene involved in regulating the STAT3/HOTAIR axis, influencing HNSCC proliferation, differentiation, and promoting metastasis by modulating related oncogenes in oral cancer. Currently, various small molecule compounds have been developed as inhibitors specifically targeting EZH2 and have gained approval for treating refractory tumors. In this review, we delve into the epigenetic regulation mediated by EZH2/PRC2 in HNSCC, with a specific focus on exploring the potential roles and mechanisms of EZH2, its crucial contribution to targeted drug therapy, and its association with cancer markers and epithelial-mesenchymal transition. Furthermore, we aim to unravel its potential as a therapeutic strategy for oral squamous cell carcinoma.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically exhibits phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption - the root of the pathogenesis - is similar in the different disease-relevant cell types. This is possible in principle, since all these cell-types are subject to effects from the same causative gene, that has the same kind of function (e.g. methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2, and find that the chromatin accessibility abnormalities in neurons are mostly distinct from those in B or T cells. This is not because the neuronal abnormalities occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that the regions disrupted in B/T cells do exhibit chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators, and suggest that blood-derived "episignatures" may not be well-suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Engineering CpG island DNA methylation in pluripotent cells through synthetic CpG-free ssDNA insertion

Cellular differentiation requires global changes to DNA methylation (DNAme), where it functions to regulate transcription factor, chromatin remodeling activity, and genome interpretation. Here, we describe a simple DNAme engineering approach in pluripotent stem cells (PSCs) that stably extends DNAme across target CpG islands (CGIs). Integration of synthetic CpG-free single-stranded DNA (ssDNA) induces a target CpG island methylation response (CIMR) in multiple PSC lines, Nt2d1 embryonal carcinoma cells, and mouse PSCs but not in highly methylated CpG island hypermethylator phenotype (CIMP)+ cancer lines. MLH1 CIMR DNAme spanned the CGI, was precisely maintained through cellular differentiation, suppressed MLH1 expression, and sensitized derived cardiomyocytes and thymic epithelial cells to cisplatin. Guidelines for CIMR editing are provided, and initial CIMR DNAme is characterized at TP53 and ONECUT1 CGIs. Collectively, this resource facilitates CpG island DNAme engineering in pluripotency and the genesis of novel epigenetic models of development and disease.

Regulation, functions and transmission of bivalent chromatin during mammalian development

Cells differentiate and progress through development guided by a dynamic chromatin landscape that mediates gene expression programmes. During development, mammalian cells display a paradoxical chromatin state: histone modifications associated with gene activation (trimethylated histone H3 Lys4 (H3K4me3)) and with gene repression (trimethylated H3 Lys27 (H3K27me3)) co-occur at promoters of developmental genes. This bivalent chromatin modification state is thought to poise important regulatory genes for expression or repression during cell-lineage specification. In this Review, we discuss recent work that has expanded our understanding of the molecular basis of bivalent chromatin and its contributions to mammalian development. We describe the factors that establish bivalency, especially histone-lysine N-methyltransferase 2B (KMT2B) and Polycomb repressive complex 2 (PRC2), and consider evidence indicating that PRC1 shapes bivalency and may contribute to its transmission between generations. We posit that bivalency is a key feature of germline and embryonic stem cells, as well as other types of stem and progenitor cells. Finally, we discuss the relevance of bivalent chromtin to human development and cancer, and outline avenues of future research.

The sound of silence: Transgene silencing in mammalian cell engineering

To elucidate principles operating in native biological systems and to develop novel biotechnologies, synthetic biology aims to build and integrate synthetic gene circuits within native transcriptional networks. The utility of synthetic gene circuits for cell engineering relies on the ability to control the expression of all constituent transgene components. Transgene silencing, defined as the loss of expression over time, persists as an obstacle for engineering primary cells and stem cells with transgenic cargos. In this review, we highlight the challenge that transgene silencing poses to the robust engineering of mammalian cells, outline potential molecular mechanisms of silencing, and present approaches for preventing transgene silencing. We conclude with a perspective identifying future research directions for improving the performance of synthetic gene circuits.

Spatiotemporal specificity of correlated DNA methylation and gene expression pairs across different human tissues and stages of brain development

DNA methylation (DNAm) that occurs on promoter regions is primarily considered to repress gene expression. Previous studies indicated that DNAm could also show positive correlations with gene expression. Both DNAm and gene expression profiles are known to be tissue- and development-specific. This study aims to investigate how DNAm and gene expression are coordinated across different human tissues and developmental stages, as well as the biological significance of such correlations. By analyzing 2,239 samples with both DNAm and gene expression data in the same human subjects obtained from six published datasets, we evaluated the correlations between gene and CpG pairs (GCPs) at cis-regions and compared significantly correlated GCPs (cGCPs) across different tissues and brains at different age groups. A total of 37,363 cGCPs was identified in the six datasets; approximately 38% of the cGCPs were positively correlated. The majority (>90%) of cGCPs was tissue- or development-specific. We also observed that the correlation direction can be opposite in different tissues and ages. Further analysis highlights the importance of cGCPs for their cellular functions and potential roles in complex traits and human diseases. For instance, the early developmental brain possessed a highly unique set of cGCPs that were associated with neurogenesis and psychiatric disorders. By assessing the epigenetic factors involved in cGCPs, we discovered novel regulatory mechanisms of positive cGCPs distinct from negative cGCPs, which were related to multiple factors, such as H3K27me3, CTCF, and JARD2. The catalogue of cGCPs compiled can be used to guide functional interpretation of genetic and epigenetic studies.

Identification of a polycomb responsive region in human HoxA cluster and its long-range interaction with polycomb enriched genomic regions

Polycomb and Trithorax group proteins (PcG, TrxG) epigenetically regulate developmental genes. These proteins bind with specific DNA elements, the Polycomb Response Element (PRE). Apart from mutations in polycomb/ trithorax proteins, altered cis-elements like PRE underlie the modified function and thus disease etiology. PREs are well studied in Drosophila, while only a few human PREs have been reported. We have identified a polycomb responsive DNA element, hPRE-HoxA3, in the intron of the HoxA3 gene. The hPRE-HoxA3 represses luciferase reporter activity in a PcG-dependent manner. The endogenous hPRE-HoxA3 element recruits PcG proteins and is enriched with repressive H3K27me3 marks, demonstrating that hPRE-HoxA3 is a part of the PcG-dependent gene regulatory network. Furthermore, it interacts with D11-12, the well-known PRE in the human Hox cluster. hPRE-Hox3 is a part of the 3-dimensional chromosomal domain organization as it is involved in the long-range interaction with other PcG enriched regions of Hox A, B, C, and D clusters.

Comparative Transcriptomics and Methylomics Reveal Adaptive Responses of Digestive and Metabolic Genes to Dietary Shift in Giant and Red Pandas

Both the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to the order Carnivora, but have changed their dietary habits to eating bamboo exclusively. The convergent evolution characteristics of their morphology, genome and gut flora have been found in the two pandas. However, the research on the convergent adaptation of their digestion and metabolism to the bamboo diet, mediated by the dietary shift of the two pandas at the gene-expression and epigenetic regulation levels, is still lacking. We therefore used RNA sequencing among five species (two pandas and three non-herbivore mammals) and bisulfite sequencing among three species (two pandas and a carnivore ferret) to sequence key digestion and metabolism tissues (stomach and small intestine). Our results provide evidence that the convergent differentially expressed genes (related to carbohydrate utilization, bile secretion, Lys and Arg metabolism, vitamin B12 utilization and cyanide detoxification) of the two pandas are adaptive responses to the bamboo diet containing low lipids, low Lys and Arg, low vitamin B12 and high cyanide. We also profiled the genome-wide methylome maps of giant panda, red panda and ferret, and the results indicated that the promoter methylation of the two pandas may regulate digestive and metabolic genes to adapt to sudden environmental changes, and then, transmit genetic information to future generations to evolve into bamboo eaters. Taken together, our study provides new insights into the molecular mechanisms of the dietary shift and the adaptation to a strict bamboo diet in both pandas using comparative transcriptomics and methylomics.

Transcription-coupled and epigenome-encoded mechanisms direct H3K4 methylation

Mono-, di-, and trimethylation of histone H3 lysine 4 (H3K4me1/2/3) are associated with transcription, yet it remains controversial whether H3K4me1/2/3 promote or result from transcription. Our previous characterizations of Arabidopsis H3K4 demethylases suggest roles for H3K4me1 in transcription. However, the control of H3K4me1 remains unexplored in Arabidopsis, in which no methyltransferase for H3K4me1 has been identified. Here, we identify three Arabidopsis methyltransferases that direct H3K4me1. Analyses of their genome-wide localization using ChIP-seq and machine learning reveal that one of the enzymes cooperates with the transcription machinery, while the other two are associated with specific histone modifications and DNA sequences. Importantly, these two types of localization patterns are also found for the other H3K4 methyltransferases in Arabidopsis and mice. These results suggest that H3K4me1/2/3 are established and maintained via interplay with transcription as well as inputs from other chromatin features, presumably enabling elaborate gene control.

PRC2, Chromatin Regulation, and Human Disease: Insights From Molecular Structure and Function

Polycomb repressive complex 2 (PRC2) is a multisubunit histone-modifying enzyme complex that mediates methylation of histone H3 lysine 27 (H3K27). Trimethylated H3K27 (H3K27me3) is an epigenetic hallmark of gene silencing. PRC2 plays a crucial role in a plethora of fundamental biological processes, and PRC2 dysregulation has been repeatedly implicated in cancers and developmental disorders. Here, we review the current knowledge on mechanisms of cellular regulation of PRC2 function, particularly regarding H3K27 methylation and chromatin targeting. PRC2-related disease mechanisms are also discussed. The mode of action of PRC2 in gene regulation is summarized, which includes competition between H3K27 methylation and acetylation, crosstalk with transcription machinery, and formation of high-order chromatin structure. Recent progress in the structural biology of PRC2 is highlighted from the aspects of complex assembly, enzyme catalysis, and chromatin recruitment, which together provide valuable insights into PRC2 function in close-to-atomic detail. Future studies on the molecular function and structure of PRC2 in the context of native chromatin and in the presence of other regulators like RNAs will continue to deepen our understanding of the stability and plasticity of developmental transcriptional programs broadly impacted by PRC2.

DNA binding by polycomb-group proteins: searching for the link to CpG islands

Polycomb group proteins predominantly exist in polycomb repressive complexes (PRCs) that cooperate to maintain the repressed state of thousands of cell-type-specific genes. Targeting PRCs to the correct sites in chromatin is essential for their function. However, the mechanisms by which PRCs are recruited to their target genes in mammals are multifactorial and complex. Here we review DNA binding by polycomb group proteins. There is strong evidence that the DNA-binding subunits of PRCs and their DNA-binding activities are required for chromatin binding and CpG targeting in cells. In vitro, CpG-specific binding was observed for truncated proteins externally to the context of their PRCs. Yet, the mere DNA sequence cannot fully explain the subset of CpG islands that are targeted by PRCs in any given cell type. At this time we find very little structural and biophysical evidence to support a model where sequence-specific DNA-binding activity is required or sufficient for the targeting of CpG-dinucleotide sequences by polycomb group proteins while they are within the context of their respective PRCs, either PRC1 or PRC2. We discuss the current knowledge and open questions on how the DNA-binding activities of polycomb group proteins facilitate the targeting of PRCs to chromatin.

Presence of H3K4me3 on Paternally Expressed Genes of the Paternal Genome From Sperm to Implantation

Genomic imprinting, parent-of-origin-specific gene expression, is controlled by differential epigenetic status of the parental chromosomes. While DNA methylation and suppressive histone modifications established during gametogenesis suppress imprinted genes on the inactive allele, how and when the expressed allele gains its active status is not clear. In this study, we asked whether the active histone-3 lysine-4 trimethylation (H3K4me3) marks remain at paternally expressed genes (PEGs) in sperm and embryos before and after fertilization using published data. Here we show that mouse sperm had the active H3K4me3 at more than half of known PEGs, and these genes were present even after fertilization. Using reciprocal cross data, we identified 13 new transient PEGs during zygotic genome activation. Next, we confirmed that the 12 out of the 13 new transient PEGs were associated with the paternal H3K4me3 in sperm. Nine out of the 12 genes were associated with the paternal H3K4me3 in zygotes. Our results show that paternal H3K4me3 marks escape inactivation during the histone-to-protamine transition that occurs during sperm maturation and are present in embryos from early zygotic stages up to implantation.

Intragenic CpG Islands and Their Impact on Gene Regulation

The mammalian genome is depleted in CG dinucleotides, except at protected regions where they cluster as CpG islands (CGIs). CGIs are gene regulatory hubs and serve as transcription initiation sites and are as expected, associated with gene promoters. Advances in genomic annotations demonstrate that a quarter of CGIs are found within genes. Such intragenic regions are repressive environments, so it is surprising that CGIs reside here and even more surprising that some resist repression and are transcriptionally active within a gene. Hence, intragenic CGI positioning within genes is not arbitrary and is instead, selected for. As a wealth of recent studies demonstrate, intragenic CGIs are embedded within genes and consequently, influence ‘host’ gene mRNA isoform length and expand transcriptome diversity.

Distinctive aspects of the placental epigenome and theories as to how they arise

The placenta has a methylome dramatically unlike that of any somatic cell type. Among other distinctions, it features low global DNA methylation, extensive “partially methylated domains” packed in dense heterochromatin and methylation of hundreds of CpG islands important in somatic development. These features attract interest in part because a substantial fraction of human cancers feature the exact same phenomena, suggesting parallels between epigenome formation in placentation and cancer. Placenta also features an expanded set of imprinted genes, some of which come about by distinctive developmental pathways. Recent discoveries, some from far outside the placental field, shed new light on how the unusual placental epigenetic state may arise. Nonetheless, key questions remain unresolved.

DNA sequence and chromatin modifiers cooperate to confer epigenetic bistability at imprinting control regions

Genomic imprinting is regulated by parental-specific DNA methylation of imprinting control regions (ICRs). Despite an identical DNA sequence, ICRs can exist in two distinct epigenetic states that are memorized throughout unlimited cell divisions and reset during germline formation. Here, we systematically study the genetic and epigenetic determinants of this epigenetic bistability. By iterative integration of ICRs and related DNA sequences to an ectopic location in the mouse genome, we first identify the DNA sequence features required for maintenance of epigenetic states in embryonic stem cells. The autonomous regulatory properties of ICRs further enabled us to create DNA-methylation-sensitive reporters and to screen for key components involved in regulating their epigenetic memory. Besides DNMT1, UHRF1 and ZFP57, we identify factors that prevent switching from methylated to unmethylated states and show that two of these candidates, ATF7IP and ZMYM2, are important for the stability of DNA and H3K9 methylation at ICRs in embryonic stem cells.

Sex-specific chromatin remodelling safeguards transcription in germ cells

Stability of the epigenetic landscape underpins maintenance of the cell-type-specific transcriptional profile. As one of the main repressive epigenetic systems, DNA methylation has been shown to be important for long-term gene silencing; its loss leads to ectopic and aberrant transcription in differentiated cells and cancer¹. The developing mouse germ line endures global changes in DNA methylation in the absence of widespread transcriptional activation. Here, using an ultra-low-input native chromatin immunoprecipitation approach, we show that following DNA demethylation the gonadal primordial germ cells undergo remodelling of repressive histone modifications, resulting in a sex-specific signature in mice. We further demonstrate that Polycomb has a central role in transcriptional control in the newly hypomethylated germline genome as the genetic loss of Ezh2 leads to aberrant transcriptional activation, retrotransposon derepression and dramatic loss of developing female germ cells. This sex-specific effect of Ezh2 deletion is explained by the distinct landscape of repressive modifications observed in male and female germ cells. Overall, our study provides insight into the dynamic interplay between repressive chromatin modifications in the context of a developmental reprogramming system.

A CpG island promoter drives the CXXC5 gene expression

CXXC5 is a member of the zinc-finger CXXC family that binds to unmethylated CpG dinucleotides. CXXC5 modulates gene expressions resulting in diverse cellular events mediated by distinct signaling pathways. However, the mechanism responsible for CXXC5 expression remains largely unknown. We found here that of the 14 annotated CXXC5 transcripts with distinct 5' untranslated regions encoding the same protein, transcript variant 2 with the highest expression level among variants represents the main transcript in cell models. The DNA segment in and at the immediate 5'-sequences of the first exon of variant 2 contains a core promoter within which multiple transcription start sites are present. Residing in a region with high G-C nucleotide content and CpG repeats, the core promoter is unmethylated, deficient in nucleosomes, and associated with active RNA polymerase-II. These findings suggest that a CpG island promoter drives CXXC5 expression. Promoter pull-down revealed the association of various transcription factors (TFs) and transcription co-regulatory proteins, as well as proteins involved in histone/chromatin, DNA, and RNA processing with the core promoter. Of the TFs, we verified that ELF1 and MAZ contribute to CXXC5 expression. Moreover, the first exon of variant 2 may contain a G-quadruplex forming region that could modulate CXXC5 expression.

The chromatin, topological and regulatory properties of pluripotency-associated poised enhancers are conserved in vivo

Poised enhancers (PEs) represent a genetically distinct set of distal regulatory elements that control the expression of major developmental genes. Before becoming activated in differentiating cells, PEs are already bookmarked in pluripotent cells with unique chromatin and topological features that could contribute to their privileged regulatory properties. However, since PEs were originally characterized in embryonic stem cells (ESC), it is currently unknown whether PEs are functionally conserved in vivo. Here, we show that the chromatin and 3D structural features of PEs are conserved among mouse pluripotent cells both in vitro and in vivo. We also uncovered that the interactions between PEs and their target genes are globally controlled by the combined action of Polycomb, Trithorax and architectural proteins. Moreover, distal regulatory sequences located close to developmental genes and displaying the typical genetic (i.e. CpG islands) and chromatin (i.e. high accessibility and H3K27me3 levels) features of PEs are commonly found across vertebrates. These putative PEs show high sequence conservation within specific vertebrate clades, with only a few being evolutionary conserved across all vertebrates. Lastly, by genetically disrupting PEs in mouse and chicken embryos, we demonstrate that these regulatory elements play essential roles during the induction of major developmental genes in vivo.

Sequence determinants, function, and evolution of CpG islands

In vertebrates, cytosine-guanine (CpG) dinucleotides are predominantly methylated, with ∼80% of all CpG sites containing 5-methylcytosine (5mC), a repressive mark associated with long-term gene silencing. The exceptions to such a globally hypermethylated state are CpG-rich DNA sequences called CpG islands (CGIs), which are mostly hypomethylated relative to the bulk genome. CGIs overlap promoters from the earliest vertebrates to humans, indicating a concerted evolutionary drive compatible with CGI retention. CGIs are characterised by DNA sequence features that include DNA hypomethylation, elevated CpG and GC content and the presence of transcription factor binding sites. These sequence characteristics are congruous with the recruitment of transcription factors and chromatin modifying enzymes, and transcriptional activation in general. CGIs colocalize with sites of transcriptional initiation in hypermethylated vertebrate genomes, however, a growing body of evidence indicates that CGIs might exert their gene regulatory function in other genomic contexts. In this review, we discuss the diverse regulatory features of CGIs, their functional readout, and the evolutionary implications associated with CGI retention in vertebrates and possibly in invertebrates.

The interplay between DNA and histone methylation: molecular mechanisms and disease implications

Methylation of cytosine in CpG dinucleotides and histone lysine and arginine residues is a chromatin modification that critically contributes to the regulation of genome integrity, replication, and accessibility. A strong correlation exists between the genome-wide distribution of DNA and histone methylation, suggesting an intimate relationship between these epigenetic marks. Indeed, accumulating literature reveals complex mechanisms underlying the molecular crosstalk between DNA and histone methylation. These in vitro and in vivo discoveries are further supported by the finding that genes encoding DNA- and histone-modifying enzymes are often mutated in overlapping human diseases. Here, we summarize recent advances in understanding how DNA and histone methylation cooperate to maintain the cellular epigenomic landscape. We will also discuss the potential implication of these insights for understanding the etiology of, and developing biomarkers and therapies for, human congenital disorders and cancers that are driven by chromatin abnormalities.

Conserved dual-mode gene regulation programs in higher eukaryotes

Recent genomic data analyses have revealed important underlying logics in eukaryotic gene regulation, such as CpG islands (CGIs)-dependent dual-mode gene regulation. In mammals, genes lacking CGIs at their promoters are generally regulated by interconversion between euchromatin and heterochromatin, while genes associated with CGIs constitutively remain as euchromatin. Whether a similar mode of gene regulation exists in non-mammalian species has been unknown. Here, through comparative epigenomic analyses, we demonstrate that the dual-mode gene regulation program is common in various eukaryotes, even in the species lacking CGIs. In cases of vertebrates or plants, we find that genes associated with high methylation level promoters are inactivated by forming heterochromatin and expressed in a context-dependent manner. In contrast, the genes with low methylation level promoters are broadly expressed and remain as euchromatin even when repressed by Polycomb proteins. Furthermore, we show that invertebrate animals lacking DNA methylation, such as fruit flies and nematodes, also have divergence in gene types: some genes are regulated by Polycomb proteins, while others are regulated by heterochromatin formation. Altogether, our study establishes gene type divergence and the resulting dual-mode gene regulation as fundamental features shared in a broad range of higher eukaryotic species.

The Role of H3K4 Trimethylation in CpG Islands Hypermethylation in Cancer

CpG methylation in transposons, exons, introns and intergenic regions is important for long-term silencing, silencing of parasitic sequences and alternative promoters, regulating imprinted gene expression and determining X chromosome inactivation. Promoter CpG islands, although rich in CpG dinucleotides, are unmethylated and remain so during all phases of mammalian embryogenesis and development, except in specific cases. The biological mechanisms that contribute to the maintenance of the unmethylated state of CpG islands remain elusive, but the modification of established DNA methylation patterns is a common feature in all types of tumors and is considered as an event that intrinsically, or in association with genetic lesions, feeds carcinogenesis. In this review, we focus on the latest results describing the role that the levels of H3K4 trimethylation may have in determining the aberrant hypermethylation of CpG islands in tumors.

Rebooting the Epigenomes during Mammalian Early Embryogenesis

Upon fertilization, terminally differentiated gametes are transformed to a totipotent zygote, which gives rise to an embryo. How parental epigenetic memories are inherited and reprogrammed to accommodate parental-to-zygotic transition remains a fundamental question in developmental biology, epigenetics, and stem cell biology. With the rapid advancement of ultra-sensitive or single-cell epigenome analysis methods, unusual principles of epigenetic reprogramming begin to be unveiled. Emerging data reveal that in many species, the parental epigenome undergoes dramatic reprogramming followed by subsequent re-establishment of the embryo epigenome, leading to epigenetic "rebooting." Here, we discuss recent progress in understanding epigenetic reprogramming and their functions during mammalian early development. We also highlight the conserved and species-specific principles underlying diverse regulation of the epigenome in early embryos during evolution.

A Structural Perspective on Gene Repression by Polycomb Repressive Complex 2

Polycomb Repressive Complex 2 (PRC2) is a major repressive chromatin complex formed by the Polycomb Group (PcG) proteins. PRC2 mediates trimethylation of histone H3 lysine 27 (H3K27me3), a hallmark of gene silencing. PRC2 is a key regulator of development, impacting many fundamental biological processes, like stem cell differentiation in mammals and vernalization in plants. Misregulation of PRC2 function is linked to a variety of human cancers and developmental disorders. In correlation with its diverse roles in development, PRC2 displays a high degree of compositional complexity and plasticity. Structural biology research over the past decade has shed light on the molecular mechanisms of the assembly, catalysis, allosteric activation, autoinhibition, chemical inhibition, dimerization and chromatin targeting of various developmentally regulated PRC2 complexes. In addition to these aspects, structure-function analysis is also discussed in connection with disease data in this chapter.

CpG Islands Shape the Epigenome Landscape

Epigenetic modifications and nucleosome positioning play an important role in modulating gene expression. However, how the patterns of epigenetic modifications and nucleosome positioning are established around promoters is not well understood. Here, we have addressed these questions in a series of genome-wide experiments coupled to a novel bioinformatic analysis approach. Our data reveal a clear correlation between CpG density, promoter activity and accumulation of active or repressive histone marks. CGI boundaries define the chromatin promoter regions that will be epigenetically modified. CpG-rich promoters are targeted by histone modifications and histone variants, while CpG-poor promoters are regulated by DNA methylation. CGIs boundaries, but not transcriptional activity, are essential determinants of H2A.Z positioning in vicinity of the promoters, suggesting that the presence of H2A.Z is not related to transcriptional control. Accordingly, H2A.Z depletion has no impact on gene expression of arrested mouse embryonic fibroblasts. Therefore, the underlying DNA sequence, the promoter CpG density and, to a lesser extent, transcriptional activity, are key factors implicated in promoter chromatin architecture.

Promoter CpG Density Predicts Downstream Gene Loss-of-Function Intolerance

The aggregation and joint analysis of large numbers of exome sequences has recently made it possible to derive estimates of intolerance to loss-of-function (LoF) variation for human genes. Here, we demonstrate strong and widespread coupling between genic LoF intolerance and promoter CpG density across the human genome. Genes downstream of the most CpG-rich promoters (top 10% CpG density) have a 67.2% probability of being highly LoF intolerant, using the LOEUF metric from gnomAD. This is in contrast to 7.4% of genes downstream of the most CpG-poor (bottom 10% CpG density) promoters. Combining promoter CpG density with exonic and promoter conservation explains 33.4% of the variation in LOEUF, and the contribution of CpG density exceeds the individual contributions of exonic and promoter conservation. We leverage this to train a simple and easily interpretable predictive model that outperforms other existing predictors and allows us to classify 1,760 genes-which are currently unascertained in gnomAD-as highly LoF intolerant or not. These predictions have the potential to aid in the interpretation of novel variants in the clinical setting. Moreover, our results reveal that high CpG density is not merely a generic feature of human promoters but is preferentially encountered at the promoters of the most selectively constrained genes, calling into question the prevailing view that CpG islands are not subject to selection.

Genome-wide DNA Methylation Signatures Are Determined by DNMT3A/B Sequence Preferences

Cytosine methylation is an important epigenetic mark, but how the distinctive patterns of DNA methylation arise remains elusive. For the first time, we systematically investigated how these patterns can be imparted by the inherent enzymatic preferences of mammalian de novo DNA methyltransferases in vitro and the extent to which this applies in cells. In a biochemical experiment, we subjected a wide variety of DNA sequences to methylation by DNMT3A or DNMT3B and then applied deep bisulfite sequencing to quantitatively determine the sequence preferences for methylation. The data show that DNMT3A prefers CpG and non-CpG sites followed by a 3'-pyrimidine, whereas DNMT3B favors a 3'-purine. Overall, we show that DNMT3A has a sequence preference for a TNC[G/A]CC context, while DNMT3B prefers TAC[G/A]GC. We extended our finding using publicly available data from mouse Dnmt1/3a/3b triple-knockout cells in which reintroduction of either DNMT3A or DNMT3B expression results in the acquisition of the same enzyme specific signature sequences observed in vitro. Furthermore, loss of DNMT3A or DNMT3B in human embryonic stem cells leads to a loss of methylation at the corresponding enzyme specific signatures. Therefore, the global DNA methylation landscape of the mammalian genome can be fundamentally determined by the inherent sequence preference of de novo methyltransferases.

Dynamic Competition of Polycomb and Trithorax in Transcriptional Programming

Predicting regulatory potential from primary DNA sequences or transcription factor binding patterns is not possible. However, the annotation of the genome by chromatin proteins, histone modifications, and differential compaction is largely sufficient to reveal the locations of genes and their differential activity states. The Polycomb Group (PcG) and Trithorax Group (TrxG) proteins are the central players in this cell type-specific chromatin organization. PcG function was originally viewed as being solely repressive and irreversible, as observed at the homeotic loci in flies and mammals. However, it is now clear that modular and reversible PcG function is essential at most developmental genes. Focusing mainly on recent advances, we review evidence for how PcG and TrxG patterns change dynamically during cell type transitions. The ability to implement cell type-specific transcriptional programming with exquisite fidelity is essential for normal development.

High-complexity regions in mammalian genomes are enriched for developmental genes

Motivation

Unique sequence regions are associated with genetic function in vertebrate genomes. However, measuring uniqueness, or absence of long repeats, along a genome is conceptually and computationally difficult. Here we use a variant of the Lempel-Ziv complexity, the match complexity, Cm, and augment it by deriving its null distribution for random sequences. We then apply Cm to the human and mouse genomes to investigate the relationship between sequence complexity and function.

Results

We implemented Cm in the program macle and show through simulation that the newly derived null distribution of Cm is accurate. This allows us to delineate high-complexity regions in the human and mouse genomes. Using our program macle2go, we find that these regions are twofold enriched for genes. Moreover, the genes contained in these regions are more than 10-fold enriched for developmental functions.

Availability and implementation

Source code for macle and macle2go is available from www.github.com/evolbioinf/macle and www.github.com/evolbioinf/macle2go, respectively; Cm browser tracks from guanine.evolbio.mgp.de/complexity.

Supplementary information

Supplementary data are available at Bioinformatics online.

Enhancer DNA methylation: implications for gene regulation

DNA methylation involves the addition of a methyl group to the fifth carbon of the pyrimidine cytosine ring (5-methylcytosine, 5mC). 5mC is widespread in vertebrate genomes where it is predominantly found within CpG dinucleotides. In mammals, 5mC participates in long-term silencing processes such as X-chromosome inactivation, genomic imprinting, somatic silencing of germline genes, and silencing of repetitive DNA elements. The evidence for 5mC as a dynamic gene-regulatory mechanism is mostly limited to specific examples, and is far from being completely understood. Recent work from diverse model systems suggests that 5mC might not always act as a dominant repressive mechanism and that hypermethylated promoters and enhancers can be permissive to transcription in vivo and in vitro. In this review, we discuss the links between 5mC and enhancer activity, and evaluate the role of this biochemical mechanism in various biological contexts.

A report on DNA sequence determinants in gene expression

The biased usage of nucleotides in coding sequence and its correlation with gene expression has been observed in several studies. A complex set of interactions between genes and other components of the expression system determine the amount of proteins produced from coding sequences. It is known that the elongation rate of polypeptide chain is affected by both codon usage bias and specific amino acid compositional constraints. Therefore, it is of interest to review local DNA-sequence elements and other positional as well as combinatorial constraints that play significant role in gene expression.

Identification of a PRC2 Accessory Subunit Required for Subtelomeric H3K27 Methylation in Neurospora crassa

Polycomb repressive complex 2 (PRC2) catalyzes methylation of histone H3 at lysine 27 (H3K27) in genomic regions of most eukaryotes and is critical for maintenance of the associated transcriptional repression. However, the mechanisms that shape the distribution of H3K27 methylation, such as recruitment of PRC2 to chromatin and/or stimulation of PRC2 activity, are unclear. Here, using a forward genetic approach in the model organism Neurospora crassa, we identified two alleles of a gene, NCU04278, encoding an unknown PRC2 accessory subunit (PAS). Loss of PAS resulted in losses of H3K27 methylation concentrated near the chromosome ends and derepression of a subset of associated subtelomeric genes. Immunoprecipitation followed by mass spectrometry confirmed reciprocal interactions between PAS and known PRC2 subunits, and sequence similarity searches demonstrated that PAS is not unique to N. crassa PAS homologs likely influence the distribution of H3K27 methylation and underlying gene repression in a variety of fungal lineages.

Understanding the interplay between CpG island-associated gene promoters and H3K4 methylation

The precise regulation of gene transcription is required to establish and maintain cell type-specific gene expression programs during multicellular development. In addition to transcription factors, chromatin, and its chemical modification, play a central role in regulating gene expression. In vertebrates, DNA is pervasively methylated at CG dinucleotides, a modification that is repressive to transcription. However, approximately 70% of vertebrate gene promoters are associated with DNA elements called CpG islands (CGIs) that are refractory to DNA methylation. CGIs integrate the activity of a range of chromatin-regulating factors that can post-translationally modify histones and modulate gene expression. This is exemplified by the trimethylation of histone H3 at lysine 4 (H3K4me3), which is enriched at CGI-associated gene promoters and correlates with transcriptional activity. Through studying H3K4me3 at CGIs it has become clear that CGIs shape the distribution of H3K4me3 and, in turn, H3K4me3 influences the chromatin landscape at CGIs. Here we will discuss our understanding of the emerging relationship between CGIs, H3K4me3, and gene expression.

Recombination may occur in the absence of transcription in the immunoglobulin heavy chain recombination centre

Developing B cells undergo V(D)J recombination to generate a vast repertoire of Ig molecules. V(D)J recombination is initiated by the RAG1/RAG2 complex in recombination centres (RCs), where gene segments become accessible to the complex. Whether transcription is the causal factor of accessibility or whether it is a side product of other processes that generate accessibility remains a controversial issue. At the IgH locus, V(D)J recombination is controlled by Eμ enhancer, which directs the transcriptional, epigenetic and recombinational events in the IgH RC. Deletion of Eμ enhancer affects both transcription and recombination, making it difficult to conclude if Eμ controls the two processes through the same or different mechanisms. By using a mouse line carrying a CpG-rich sequence upstream of Eμ enhancer and analyzing transcription and recombination at the single-cell level, we found that recombination could occur in the RC in the absence of detectable transcription, suggesting that Eμ controls transcription and recombination through distinct mechanisms. Moreover, while the normally Eμ-dependent transcription and demethylating activities were impaired, recruitment of chromatin remodeling complexes was unaffected. RAG1 was efficiently recruited, thus compensating for the defective transcription-associated recruitment of RAG2, and providing a mechanistic basis for RAG1/RAG2 assembly to initiate V(D)J recombination.

Telomere dysfunction cooperates with epigenetic alterations to impair murine embryonic stem cell fate commitment

The precise relationship between epigenetic alterations and telomere dysfunction is still an extant question. Previously, we showed that eroded telomeres lead to differentiation instability in murine embryonic stem cells (mESCs) via DNA hypomethylation at pluripotency-factor promoters. Here, we uncovered that telomerase reverse transcriptase null (Tert^-/-) mESCs exhibit genome-wide alterations in chromatin accessibility and gene expression during differentiation. These changes were accompanied by an increase of H3K27me3 globally, an altered chromatin landscape at the Pou5f1/Oct4 promoter, and a refractory response to differentiation cues. Inhibition of the Polycomb Repressive Complex 2 (PRC2), an H3K27 tri-methyltransferase, exacerbated the impairment in differentiation and pluripotency gene repression in Tert^-/-mESCs but not wild-type mESCs, whereas inhibition of H3K27me3 demethylation led to a partial rescue of the Tert^-/- phenotype. These data reveal a new interdependent relationship between H3K27me3 and telomere integrity in stem cell lineage commitment that may have implications in aging and cancer.

A Dimeric Structural Scaffold for PRC2-PCL Targeting to CpG Island Chromatin

Diverse accessory subunits are involved in the recruitment of polycomb repressive complex 2 (PRC2) to CpG island (CGI) chromatin. Here we report the crystal structure of a SUZ12-RBBP4 complex bound to fragments of the accessory subunits PHF19 and JARID2. Unexpectedly, this complex adopts a dimeric structural architecture, accounting for PRC2 self-association that has long been implicated. The intrinsic PRC2 dimer is formed via domain swapping involving RBBP4 and the unique C2 domain of SUZ12. MTF2 and PHF19 associate with PRC2 at around the dimer interface and stabilize the dimer. Conversely, AEBP2 binding results in a drastic movement of the C2 domain, disrupting the intrinsic PRC2 dimer. PRC2 dimerization enhances CGI DNA binding by PCLs in pairs in vitro, reminiscent of the widespread phenomenon of transcription factor dimerization in active transcription. Loss of PRC2 dimerization impairs histone H3K27 trimethylation (H3K27me3) on chromatin at developmental gene loci in mouse embryonic stem cells.

DNA G-quadruplex stability, position and chromatin accessibility are associated with CpG island methylation

CpG islands (CGI) are genomic regions associated with gene promoters and involved in gene expression regulation. Despite their high CpG content and unlike bulk genomic DNA methylation pattern, these regions are usually hypomethylated. So far, the mechanisms controlling the CGI methylation patterning remain unclear. G-quadruplex (G4) structures can inhibit DNA methyltransferases 1 enzymatic activity, leading to CGI hypomethylation. Our aim was to analyse the association of G4 forming sequences (G4FS) and CGI methylation as well as to determine the intrinsic and extrinsic characteristics of G4FS that may modulate this phenomenon. Using methylation data from human embryonic stem cells (hESCs) and three hESC-derived populations, we showed that hypomethylated CpGs located inside CGI (CGI/CpG) tend to be associated with highly stable G4FS (Minimum free energy ≤ -30 kcal·mol^-1 ). The association of highly stable G4FS and hypomethylation tend to be stronger when these structures are located at shorter distances (~ < 150 bp) from GCI/CpGs, when G4FS and CpGs are located within open chromatin and G4FS are inside CGI. Moreover, this association is not strongly influenced by the CpG content of CGI. Conversely, highly methylated CGI/CpG tend to be associated with low stability G4FS. Although CpGs inside CGIs without a G4FS tend to be more methylated, high stability G4FS within CGI neighbourhood were associated with decreased methylation. In summary, our data indicate that G4FS may act as protective cis elements against CGI methylation, and this effect seems to be influenced by the G4FS folding potential, its presence within CGI, CpG distance from G4FS and chromatin accessibility.

Non-core Subunits of the PRC2 Complex Are Collectively Required for Its Target-Site Specificity

The Polycomb repressive complex 2 (PRC2) catalyzes H3K27 methylation across the genome, which impacts transcriptional regulation and is critical for establishment of cell identity. Because of its essential function during development and in cancer, understanding the delineation of genome-wide H3K27 methylation patterns has been the focus of intense investigation. PRC2 methylation activity is abundant and dispersed throughout the genome, but the highest activity is specifically directed to a subset of target sites that are stably occupied by the complex and highly enriched for H3K27me3. Here, we show, by systematically knocking out single and multiple non-core subunits of the PRC2 complex in mouse embryonic stem cells, that they each contribute to directing PRC2 activity to target sites. Furthermore, combined knockout of six non-core subunits reveals that, while dispensable for global H3K27 methylation levels, the non-core PRC2 subunits are collectively required for focusing H3K27me3 activity to specific sites in the genome.

Molecular Mechanisms Directing PRC2 Recruitment and H3K27 Methylation

The polycomb repressive complex 2 (PRC2) is a chromatin-associated methyltransferase catalyzing mono-, di-, and trimethylation of lysine 27 on histone H3 (H3K27). This activity is required for normal organismal development and maintenance of gene expression patterns to uphold cell identity. PRC2 function is often deregulated in disease and is a promising candidate for therapeutic targeting in cancer. In this review, we discuss the molecular mechanisms proposed to take part in modulating PRC2 recruitment and shaping H3K27 methylation patterns across the genome. This includes consideration of factors influencing PRC2 residence time on chromatin and PRC2 catalytic activity with a focus on the mechanisms giving rise to regional preferences and differential deposition of H3K27 methylation. We further discuss existing evidence for functional diversity between distinct subsets of PRC2 complexes with the aim of extracting key concepts and highlighting major open questions toward a more complete understanding of PRC2 function.

Chromatin folding and nuclear architecture: PRC1 function in 3D

Embryonic development requires the intricate balance between the expansion and specialisation of defined cell types in time and space. The gene expression programmes that underpin this balance are regulated, in part, by modulating the chemical and structural state of chromatin. Polycomb repressive complexes (PRCs), a family of essential developmental regulators, operate at this level to stabilise or perpetuate a repressed but transcriptionally poised chromatin configuration. This dynamic state is required to control the timely initiation of productive gene transcription during embryonic development. The two major PRCs cooperate to target the genome, but it is PRC1 that appears to be the primary effector that controls gene expression. In this review I will discuss recent findings relating to how PRC1 alters chromatin accessibility, folding and global 3D nuclear organisation to control gene transcription.

CG dinucleotides enhance promoter activity independent of DNA methylation

Most mammalian RNA polymerase II initiation events occur at CpG islands, which are rich in CpGs and devoid of DNA methylation. Despite their relevance for gene regulation, it is unknown to what extent the CpG dinucleotide itself actually contributes to promoter activity. To address this question, we determined the transcriptional activity of a large number of chromosomally integrated promoter constructs and monitored binding of transcription factors assumed to play a role in CpG island activity. This revealed that CpG density significantly improves motif-based prediction of transcription factor binding. Our experiments also show that high CpG density alone is insufficient for transcriptional activity, yet results in increased transcriptional output when combined with particular transcription factor motifs. However, this CpG contribution to promoter activity is independent of DNA methyltransferase activity. Together, this refines our understanding of mammalian promoter regulation as it shows that high CpG density within CpG islands directly contributes to an environment permissive for full transcriptional activity.

Regulation of Tumor Suppressor Gene CDKN2A and Encoded p16-INK4a Protein by Covalent Modifications

CDKN2A is one of the most studied tumor suppressor genes. It encodes the p16-INK4a protein that plays a critical role in the cell cycle progression, differentiation, senescence, and apoptosis. Mutations in CDKN2A or dysregulation of its functional activity are frequently associated with various types of human cancer. As a cyclin-dependent kinase inhibitor, p16-INK4a forms a complex with cyclin-dependent kinases 4/6 (CDK4/6) thereby competing with cyclin D. It is believed that the helix-turn-helix structures in the content of tandem ankyrin repeats in p16-INK4a are required for the protein interaction with CDK4. Until recently, the mechanisms considered to be involved in the regulation of p16-INK4a functions and cancer development have been mutations in DNA, homozygous or heterozygous gene loss, and methylation of CDKN2A promoter region. In this review, we discuss recent findings on the regulation of p16-INK4a by covalent modifications at both transcriptional and post-translational levels.

DNA G-quadruplex structures mold the DNA methylome

Control of DNA methylation level is critical for gene regulation, and the factors that govern hypomethylation at CpG islands (CGIs) are still being uncovered. Here, we provide evidence that G-quadruplex (G4) DNA secondary structures are genomic features that influence methylation at CGIs. We show that the presence of G4 structure is tightly associated with CGI hypomethylation in the human genome. Surprisingly, we find that these G4 sites are enriched for DNA methyltransferase 1 (DNMT1) occupancy, which is consistent with our biophysical observations that DNMT1 exhibits higher binding affinity for G4s as compared to duplex, hemi-methylated, or single-stranded DNA. The biochemical assays also show that the G4 structure itself, rather than sequence, inhibits DNMT1 enzymatic activity. Based on these data, we propose that G4 formation sequesters DNMT1 thereby protecting certain CGIs from methylation and inhibiting local methylation.

Three classes of response elements for human PRC2 and MLL1/2-Trithorax complexes

Polycomb group (PcG) and Trithorax group (TrxG) proteins are essential for maintaining epigenetic memory in both embryonic stem cells and differentiated cells. To date, how they are localized to hundreds of specific target genes within a vertebrate genome had remained elusive. Here, by focusing on short cis-acting DNA elements of single functions, we discovered three classes of response elements in human genome: Polycomb response elements (PREs), Trithorax response elements (TREs) and Polycomb/Trithorax response elements (P/TREs). In particular, the four PREs (PRE14, 29, 39 and 48) are the first set of, to our knowledge, bona fide vertebrate PREs ever discovered, while many previously reported Drosophila or vertebrate PREs are likely P/TREs. We further demonstrated that YY1 and CpG islands are specifically enriched in the four TREs (PRE30, 41, 44 and 55), but not in the PREs. The three classes of response elements as unraveled in this study should guide further global investigation and open new doors for a deeper understanding of PcG and TrxG mechanisms in vertebrates.