Bio-CAP: a versatile and highly sensitive technique to purify and characterise regions of non-methylated DNA

Extensive DNA methylome rearrangement during early lamprey embryogenesis

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

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.

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.

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.

Polycomb repressive complex 1 shapes the nucleosome landscape but not accessibility at target genes

Polycomb group (PcG) proteins are transcriptional repressors that play important roles in regulating gene expression during animal development. In vitro experiments have shown that PcG protein complexes can compact chromatin to limit the activity of chromatin remodeling enzymes and access of the transcriptional machinery to DNA. In fitting with these ideas, gene promoters associated with PcG proteins have been reported to be less accessible than other gene promoters. However, it remains largely untested in vivo whether PcG proteins define chromatin accessibility or other chromatin features. To address this important question, we examine the chromatin accessibility and nucleosome landscape at PcG protein-bound promoters in mouse embryonic stem cells using the assay for transposase accessible chromatin (ATAC)-seq. Combined with genetic ablation strategies, we unexpectedly discover that although PcG protein-occupied gene promoters exhibit reduced accessibility, this does not rely on PcG proteins. Instead, the Polycomb repressive complex 1 (PRC1) appears to play a unique role in driving elevated nucleosome occupancy and decreased nucleosomal spacing in Polycomb chromatin domains. Our new genome-scale observations argue, in contrast to the prevailing view, that PcG proteins do not significantly affect chromatin accessibility and highlight an underappreciated complexity in the relationship between chromatin accessibility, the nucleosome landscape, and PcG-mediated transcriptional repression.

FBXL19 recruits CDK-Mediator to CpG islands of developmental genes priming them for activation during lineage commitment

CpG islands are gene regulatory elements associated with the majority of mammalian promoters, yet how they regulate gene expression remains poorly understood. Here, we identify FBXL19 as a CpG island-binding protein in mouse embryonic stem (ES) cells and show that it associates with the CDK-Mediator complex. We discover that FBXL19 recruits CDK-Mediator to CpG island-associated promoters of non-transcribed developmental genes to prime these genes for activation during cell lineage commitment. We further show that recognition of CpG islands by FBXL19 is essential for mouse development. Together this reveals a new CpG island-centric mechanism for CDK-Mediator recruitment to developmental gene promoters in ES cells and a requirement for CDK-Mediator in priming these developmental genes for activation during cell lineage commitment.

The SET1 Complex Selects Actively Transcribed Target Genes via Multivalent Interaction with CpG Island Chromatin

Chromatin modifications and the promoter-associated epigenome are important for the regulation of gene expression. However, the mechanisms by which chromatin-modifying complexes are targeted to the appropriate gene promoters in vertebrates and how they influence gene expression have remained poorly defined. Here, using a combination of live-cell imaging and functional genomics, we discover that the vertebrate SET1 complex is targeted to actively transcribed gene promoters through CFP1, which engages in a form of multivalent chromatin reading that involves recognition of non-methylated DNA and histone H3 lysine 4 trimethylation (H3K4me3). CFP1 defines SET1 complex occupancy on chromatin, and its multivalent interactions are required for the SET1 complex to place H3K4me3. In the absence of CFP1, gene expression is perturbed, suggesting that normal targeting and function of the SET1 complex are central to creating an appropriately functioning vertebrate promoter-associated epigenome.

•

The CFP1/SET1 complex engages in dynamic and stable chromatin-binding events

•

CFP1 uses multivalent chromatin interactions to select active CpG island promoters

•

SET1A occupancy at CpG island promoters is predominately defined by CFP1

•

CFP1 targets SET1 to shape promoter-associated H3K4me3 and gene expression

Biochemical Identification of Nonmethylated DNA by BioCAP-Seq

CpG islands are regions of vertebrate genomes that often function as gene regulatory elements and are associated with most gene promoters. CpG island elements usually contain nonmethylated CpG dinucleotides, while the remainder of the genome is pervasively methylated. We developed a biochemical approach called biotinylated CxxC affinity purification (BioCAP) to unbiasedly isolate regions of the genome that contain nonmethylated CpG dinucleotides. The resulting highly pure nonmethylated DNA is easily analyzed by quantitative PCR to interrogate specific loci or via massively parallel sequencing to yield genome-wide profiles.

MLL-AF4 Spreading Identifies Binding Sites that Are Distinct from Super-Enhancers and that Govern Sensitivity to DOT1L Inhibition in Leukemia

Understanding the underlying molecular mechanisms of defined cancers is crucial for effective personalized therapies. Translocations of the mixed-lineage leukemia (MLL) gene produce fusion proteins such as MLL-AF4 that disrupt epigenetic pathways and cause poor-prognosis leukemias. Here, we find that at a subset of gene targets, MLL-AF4 binding spreads into the gene body and is associated with the spreading of Menin binding, increased transcription, increased H3K79 methylation (H3K79me2/3), a disruption of normal H3K36me3 patterns, and unmethylated CpG regions in the gene body. Compared to other H3K79me2/3 marked genes, MLL-AF4 spreading gene expression is downregulated by inhibitors of the H3K79 methyltransferase DOT1L. This sensitivity mediates synergistic interactions with additional targeted drug treatments. Therefore, epigenetic spreading and enhanced susceptibility to epidrugs provides a potential marker for better understanding combination therapies in humans.

•

MLL-AF4 binding requires an unmethylated CpG (uCpG) island and Menin

•

MLL-AF4 and Menin can spread into the gene body of some targets

•

Spreading targets are highly transcribed and have an aberrant chromatin signature

•

Spreading of MLL-AF4 is a predictor of sensitivity to DOT1L inhibitors

Improved Prediction of Non-methylated Islands in Vertebrates Highlights Different Characteristic Sequence Patterns

Non-methylated islands (NMIs) of DNA are genomic regions that are important for gene regulation and development. A recent study of genome-wide non-methylation data in vertebrates by Long et al. (eLife 2013;2:e00348) has shown that many experimentally identified non-methylated regions do not overlap with classically defined CpG islands which are computationally predicted using simple DNA sequence features. This is especially true in cold-blooded vertebrates such as Danio rerio (zebrafish). In order to investigate how predictive DNA sequence is of a region’s methylation status, we applied a supervised learning approach using a spectrum kernel support vector machine, to see if a more complex model and supervised learning can be used to improve non-methylated island prediction and to understand the sequence properties of these regions. We demonstrate that DNA sequence is highly predictive of methylation status, and that in contrast to existing CpG island prediction methods our method is able to provide more useful predictions of NMIs genome-wide in all vertebrate organisms that were studied. Our results also show that in cold-blooded vertebrates (Anolis carolinensis, Xenopus tropicalis and Danio rerio) where genome-wide classical CpG island predictions consist primarily of false positives, longer primarily AT-rich DNA sequence features are able to identify these regions much more accurately.

DNA methylation is one of the most heavily studied epigenetic modifications, with active research into the topic taking place for over 30 years. Over that time a relationship between non-methylated stretches of genomic DNA and the regulation of nearby genes has been identified, and it was noted that approximately half of all human genes have an unexpectedly high number of unmethylated CpG dinucleotides in their promoters. Nevertheless, until relatively recently we haven’t been able to directly measure methylation gnome-wide, leading to the widespread use of CpG-rich regions (or CpG Islands) as a proxy for unmethylated genomic regions. Fortunately, recently developed experiments finally allow us to look at DNA methylation genome-wide, and have shown that CpG islands do not overlap true non-methylated islands to the extent we would have hoped, especially in cold-blooded vertebrates. We show that by using these experimentally determined non-methylated regions as training data, we are able to improve the prediction of non-methylated regions genome-wide in comparison to classical CpG island prediction methods, in both warm-blooded and cold-blooded vertebrates. Additionally, we show that while CpG content is an important predictor of methylation status in warm-blooded vertebrates, in cold-blooded vertebrates higher predictive accuracy can be achieved using longer stretches of relatively AT-rich DNA sequence.

Successful transmission and transcriptional deployment of a human chromosome via mouse male meiosis

Most human aneuploidies originate maternally, due in part to the presence of highly stringent checkpoints during male meiosis. Indeed, male sterility is common among aneuploid mice used to study chromosomal abnormalities, and male germline transmission of exogenous DNA has been rarely reported. Here we show that, despite aberrant testis architecture, males of the aneuploid Tc1 mouse strain produce viable sperm and transmit human chromosome 21 to create aneuploid offspring. In these offspring, we mapped transcription, transcriptional initiation, enhancer activity, non-methylated DNA, and transcription factor binding in adult tissues. Remarkably, when compared with mice derived from female passage of human chromosome 21, the chromatin condensation during spermatogenesis and the extensive epigenetic reprogramming specific to male germline transmission resulted in almost indistinguishable patterns of transcriptional deployment. Our results reveal an unexpected tolerance of aneuploidy during mammalian spermatogenesis, and the surprisingly robust ability of mouse developmental machinery to accurately deploy an exogenous chromosome, regardless of germline transmission.

Protection of CpG islands from DNA methylation is DNA-encoded and evolutionarily conserved

DNA methylation is a repressive epigenetic modification that covers vertebrate genomes. Regions known as CpG islands (CGIs), which are refractory to DNA methylation, are often associated with gene promoters and play central roles in gene regulation. Yet how CGIs in their normal genomic context evade the DNA methylation machinery and whether these mechanisms are evolutionarily conserved remains enigmatic. To address these fundamental questions we exploited a transchromosomic animal model and genomic approaches to understand how the hypomethylated state is formed in vivo and to discover whether mechanisms governing CGI formation are evolutionarily conserved. Strikingly, insertion of a human chromosome into mouse revealed that promoter-associated CGIs are refractory to DNA methylation regardless of host species, demonstrating that DNA sequence plays a central role in specifying the hypomethylated state through evolutionarily conserved mechanisms. In contrast, elements distal to gene promoters exhibited more variable methylation between host species, uncovering a widespread dependence on nucleotide frequency and occupancy of DNA-binding transcription factors in shaping the DNA methylation landscape away from gene promoters. This was exemplified by young CpG rich lineage-restricted repeat sequences that evaded DNA methylation in the absence of co-evolved mechanisms targeting methylation to these sequences, and species specific DNA binding events that protected against DNA methylation in CpG poor regions. Finally, transplantation of mouse chromosomal fragments into the evolutionarily distant zebrafish uncovered the existence of a mechanistically conserved and DNA-encoded logic which shapes CGI formation across vertebrate species.

Evolution of Epigenetic Regulation in Vertebrate Genomes

Empirical models of sequence evolution have spurred progress in the field of evolutionary genetics for decades. We are now realizing the importance and complexity of the eukaryotic epigenome. While epigenome analysis has been applied to genomes from single-cell eukaryotes to human, comparative analyses are still relatively few and computational algorithms to quantify epigenome evolution remain scarce. Accordingly, a quantitative model of epigenome evolution remains to be established. We review here the comparative epigenomics literature and synthesize its overarching themes. We also suggest one mechanism, transcription factor binding site (TFBS) turnover, which relates sequence evolution to epigenetic conservation or divergence. Lastly, we propose a framework for how the field can move forward to build a coherent quantitative model of epigenome evolution.

BsRADseq: screening DNA methylation in natural populations of non-model species

Epigenetic modifications are expected to occur at a much faster rate than genetic mutations, potentially causing isolated populations to stochastically drift apart, or if they are subjected to different selective regimes, to directionally diverge. A high level of genome‐wide epigenetic divergence between individuals occupying distinct habitats is therefore predicted. Here, we introduce bisulfite‐converted restriction site associated DNA sequencing (bsRADseq), an approach to quantify the level of DNA methylation differentiation across multiple individuals. This reduced representation method is flexible in the extent of DNA sequence interrogated. We showcase its applicability in three natural systems, each comprising individuals adapted to divergent environments: a diploid plant (Heliosperma, Caryophyllaceae), a tetraploid plant (Dactylorhiza, Orchidaceae) and an animal (Gasterosteusaculeatus, Gasterosteidae). We present a robust bioinformatic pipeline, combining tools for RAD locus assembly, SNP calling, bisulfite‐converted read mapping and DNA methylation calling to analyse bsRADseq data with or without a reference genome. Importantly, our approach accurately distinguishes between SNPs and methylation polymorphism (SMPs). Although DNA methylation frequency between different positions of a genome varies widely, we find a surprisingly high consistency in the methylation profile between individuals thriving in divergent ecological conditions, particularly in Heliosperma. This constitutive stability points to significant molecular or developmental constraints acting on DNA methylation variation. Altogether, by combining the flexibility of RADseq with the accuracy of bisulfite sequencing in quantifying DNA methylation, the bsRADseq methodology and our bioinformatic pipeline open up the opportunity for genome‐wide epigenetic investigations of evolutionary and ecological relevance in non‐model species, independent of their genomic features.

Whole-genome bisulfite sequencing maps from multiple human tissues reveal novel CpG islands associated with tissue-specific regulation

CpG islands (CGIs) are one of the most widely studied regulatory features of the human genome, with critical roles in development and disease. Despite such significance and the original epigenetic definition, currently used CGI sets are typically predicted from DNA sequence characteristics. Although CGIs are deeply implicated in practical analyses of DNA methylation, recent studies have shown that such computational annotations suffer from inaccuracies. Here we used whole-genome bisulfite sequencing from 10 diverse human tissues to identify a comprehensive, experimentally obtained, single-base resolution CGI catalog. In addition to the unparalleled annotation precision, our method is free from potential bias due to arbitrary sequence features or probe affinity differences. In addition to clarifying substantial false positives in the widely used University of California Santa Cruz (UCSC) annotations, our study identifies numerous novel epigenetic loci. In particular, we reveal significant impact of transposable elements on the epigenetic regulatory landscape of the human genome and demonstrate ubiquitous presence of transcription initiation at CGIs, including alternative promoters in gene bodies and non-coding RNAs in intergenic regions. Moreover, coordinated DNA methylation and chromatin modifications mark tissue-specific enhancers at novel CGIs. Enrichment of specific transcription factor binding from ChIP-seq supports mechanistic roles of CGIs on the regulation of tissue-specific transcription. The new CGI catalog provides a comprehensive and integrated list of genomic hotspots of epigenetic regulation.

ESCI award lecture: regulation, function and biomarker potential of DNA methylation

Methylation of DNA and modifications of histones have emerged as intricately involved in gene regulation as they cross-talk and respond in multiple ways to the activity of transcription factors. Measuring these epigenome components has become a powerful tool to identify regulatory principles and biomarkers that predict cellular state during development or disease. Here, I will focus on DNA methylation as a reversible epigenetic modification of DNA that has been studied in great detail at the level of the genome. Recent advances in sequencing have identified unexpected dynamics of this modification, which are tightly linked to gene regulation. Understanding how DNA methylation patterns are read and how they contribute to regulation will be critical to interpret and utilize genomic maps of DNA methylation. As these patterns are dynamic during cellular differentiation and perturbed in disease, they present an opportunity to use DNA methylation as a biomarker.

piRBase: a web resource assisting piRNA functional study

piRNAs are a class of small RNAs that is most abundantly expressed in the animal germ line. Presently, substantial research is going on to reveal the functions of piRNAs in the epigenetic and post-transcriptional regulation of transposons and genes. A piRNA database for collection, annotation and structuring of these data will be a valuable contribution to the field, and we have therefore developed the piRBase platform which integrates various piRNA-related high-throughput data. piRBase has the largest collection of piRNAs among existing databases, and contains at present 77 million piRNA sequences from nine organisms. Repeat-derived and gene-derived piRNAs, which possibly participate in the regulation of the corresponding elements, have been given particular attention. Furthermore, epigenetic data and reported piRNA targets were also collected. To our knowledge, this is the first piRNA database that systematically integrates epigenetic and post-transcriptional regulation data to support piRNA functional analysis. We believe that piRBase will contribute to a better understanding of the piRNA functions. Database URL: http://www.regulatoryrna.org/database/piRNA/

Application of histone modification-specific interaction domains as an alternative to antibodies

Post-translational modifications (PTMs) of histones constitute a major chromatin indexing mechanism, and their proper characterization is of highest biological importance. So far, PTM-specific antibodies have been the standard reagent for studying histone PTMs despite caveats such as lot-to-lot variability of specificity and binding affinity. Herein, we successfully employed naturally occurring and engineered histone modification interacting domains for detection and identification of histone PTMs and ChIP-like enrichment of different types of chromatin. Our results demonstrate that histone interacting domains are robust and highly specific reagents that can replace or complement histone modification antibodies. These domains can be produced recombinantly in Escherichia coli at low cost and constant quality. Protein design of reading domains allows for generation of novel specificities, addition of affinity tags, and preparation of PTM binding pocket variants as matching negative controls, which is not possible with antibodies.

DNA methylation profiles define stem cell identity and reveal a tight embryonic-extraembryonic lineage boundary

Embryonic (ES) and epiblast (EpiSC) stem cells are pluripotent but committed to an embryonic lineage fate. Conversely, trophoblast (TS) and extraembryonic endoderm (XEN) stem cells contribute predominantly to tissues of the placenta and yolk sac, respectively. Here we show that each of these four stem cell types is defined by a unique DNA methylation profile. Despite their distinct developmental origin, TS and XEN cells share key epigenomic hallmarks, chiefly characterized by robust DNA methylation of embryo-specific developmental regulators, as well as a subordinate role of 5-hydroxymethylation. We also observe a substantial methylation reinforcement of pre-existing epigenetic repressive marks that specifically occurs in extraembryonic stem cells compared to in vivo tissue, presumably due to continued high Dnmt3b expression levels. These differences establish a major epigenetic barrier between the embryonic and extraembryonic stem cell types. In addition, epigenetic lineage boundaries also separate the two extraembryonic stem cell types by mutual repression of key lineage-specific transcription factors. Thus, global DNA methylation patterns are a defining feature of each stem cell type that underpin lineage commitment and differentiative potency of early embryo-derived stem cells. Our detailed methylation profiles identify a cohort of developmentally regulated sequence elements, such as orphan CpG islands, that will be most valuable to uncover novel transcriptional regulators and pivotal "gatekeeper" genes in pluripotency and lineage differentiation.

Epigenetic conservation at gene regulatory elements revealed by non-methylated DNA profiling in seven vertebrates

Two-thirds of gene promoters in mammals are associated with regions of non-methylated DNA, called CpG islands (CGIs), which counteract the repressive effects of DNA methylation on chromatin. In cold-blooded vertebrates, computational CGI predictions often reside away from gene promoters, suggesting a major divergence in gene promoter architecture across vertebrates. By experimentally identifying non-methylated DNA in the genomes of seven diverse vertebrates, we instead reveal that non-methylated islands (NMIs) of DNA are a central feature of vertebrate gene promoters. Furthermore, NMIs are present at orthologous genes across vast evolutionary distances, revealing a surprising level of conservation in this epigenetic feature. By profiling NMIs in different tissues and developmental stages we uncover a unifying set of features that are central to the function of NMIs in vertebrates. Together these findings demonstrate an ancient logic for NMI usage at gene promoters and reveal an unprecedented level of epigenetic conservation across vertebrate evolution. DOI:http://dx.doi.org/10.7554/eLife.00348.001.

KDM2B links the Polycomb Repressive Complex 1 (PRC1) to recognition of CpG islands

CpG islands (CGIs) are associated with most mammalian gene promoters. A subset of CGIs act as polycomb response elements (PREs) and are recognized by the polycomb silencing systems to regulate expression of genes involved in early development. How CGIs function mechanistically as nucleation sites for polycomb repressive complexes remains unknown. Here we discover that KDM2B (FBXL10) specifically recognizes non-methylated DNA in CGIs and recruits the polycomb repressive complex 1 (PRC1). This contributes to histone H2A lysine 119 ubiquitylation (H2AK119ub1) and gene repression. Unexpectedly, we also find that CGIs are occupied by low levels of PRC1 throughout the genome, suggesting that the KDM2B-PRC1 complex may sample CGI-associated genes for susceptibility to polycomb-mediated silencing. These observations demonstrate an unexpected and direct link between recognition of CGIs by KDM2B and targeting of the polycomb repressive system. This provides the basis for a new model describing the functionality of CGIs as mammalian PREs.DOI:http://dx.doi.org/10.7554/eLife.00205.001.