Arie B. Brinkman

Dynamic Readers for 5-(Hydroxy)Methylcytosine and Its Oxidized Derivatives

Tet proteins oxidize 5-methylcytosine (mC) to generate 5-hydroxymethyl (hmC), 5-formyl (fC), and 5-carboxylcytosine (caC). The exact function of these oxidative cytosine bases remains elusive. We applied quantitative mass-spectrometry-based proteomics to identify readers for mC and hmC in mouse embryonic stem cells (mESC), neuronal progenitor cells (NPC), and adult mouse brain tissue. Readers for these modifications are only partially overlapping, and some readers, such as Rfx proteins, display strong specificity. Interactions are dynamic during differentiation, as for example evidenced by the mESC-specific binding of Klf4 to mC and the NPC-specific binding of Uhrf2 to hmC, suggesting specific biological roles for mC and hmC. Oxidized derivatives of mC recruit distinct transcription regulators as well as a large number of DNA repair proteins in mouse ES cells, implicating the DNA damage response as a major player in active DNA demethylation.

Quantitative comparison of genome-wide DNA methylation mapping technologies

DNA methylation plays a key role in regulating eukaryotic gene expression. Although mitotically heritable and stable over time, patterns of DNA methylation frequently change in response to cell differentiation, disease and environmental influences. Several methods have been developed to map DNA methylation on a genomic scale. Here, we benchmark four of these approaches by analyzing two human embryonic stem cell lines derived from genetically unrelated embryos and a matched pair of colon tumor and adjacent normal colon tissue obtained from the same donor. Our analysis reveals that methylated DNA immunoprecipitation sequencing (MeDIP-seq), methylated DNA capture by affinity purification (MethylCap-seq), reduced representation bisulfite sequencing (RRBS) and the Infinium HumanMethylation27 assay all produce accurate DNA methylation data. However, these methods differ in their ability to detect differentially methylated regions between pairs of samples. We highlight strengths and weaknesses of the four methods and give practical recommendations for the design of epigenomic case-control studies.DNA methylation plays a key role in regulating eukaryotic gene expression. Although mitotically heritable and stable over time, patterns of DNA methylation frequently change in response to cell differentiation, disease and environmental influences. Several methods have been developed to map DNA methylation on a genomic scale. Here, we benchmark four of these approaches by analyzing two human embryonic stem cell lines derived from genetically unrelated embryos and a matched pair of colon tumor and adjacent normal colon tissue obtained from the same donor. Our analysis reveals that methylated DNA immunoprecipitation sequencing (MeDIP-seq), methylated DNA capture by affinity purification (MethylCap-seq), reduced representation bisulfite sequencing (RRBS) and the Infinium HumanMethylation27 assay all produce accurate DNA methylation data. However, these methods differ in their ability to detect differentially methylated regions between pairs of samples. We highlight strengths and weaknesses of the four methods and give practical recommendations for the design of epigenomic case-control studies.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences

We analysed whole genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. 93 protein-coding cancer genes carried likely driver mutations. Some non-coding regions exhibited high mutation frequencies but most have distinctive structural features probably causing elevated mutation rates and do not harbour driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed 12 base substitution and six rearrangement signatures. Three rearrangement signatures, characterised by tandem duplications or deletions, appear associated with defective homologous recombination based DNA repair: one with deficient BRCA1 function; another with deficient BRCA1 or BRCA2 function; the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operative, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

MBD2/NuRD and MBD3/NuRD, Two Distinct Complexes with Different Biochemical and Functional Properties

ABSTRACT The human genome contains a number of methyl CpG binding proteins that translate DNA methylation into a physiological response. To gain insight into the function of MBD2 and MBD3, we first applied protein tagging and mass spectrometry. We show that MBD2 and MBD3 assemble into mutually exclusive distinct Mi-2/NuRD-like complexes, called MBD2/NuRD and MBD3/NuRD. We identified DOC-1, a putative tumor suppressor, as a novel core subunit of MBD2/NuRD as well as MBD3/NuRD. PRMT5 and its cofactor MEP50 were identified as specific MBD2/NuRD interactors. PRMT5 stably and specifically associates with and methylates the RG-rich N terminus of MBD2. Chromatin immunoprecipitation experiments revealed that PRMT5 and MBD2 are recruited to CpG islands in a methylation-dependent manner in vivo and that H4R3, a substrate of PRMT, is methylated at these loci. Our data show that MBD2/NuRD and MBD3/NuRD are distinct protein complexes with different biochemical and functional properties.

Sequential ChIP-bisulfite sequencing enables direct genome-scale investigation of chromatin and DNA methylation cross-talk

Cross-talk between DNA methylation and histone modifications drives the establishment of composite epigenetic signatures and is traditionally studied using correlative rather than direct approaches. Here, we present sequential ChIP-bisulfite-sequencing (ChIP-BS-seq) as an approach to quantitatively assess DNA methylation patterns associated with chromatin modifications or chromatin-associated factors directly. A chromatin-immunoprecipitation (ChIP)-capturing step is used to obtain a restricted representation of the genome occupied by the epigenetic feature of interest, for which a single-base resolution DNA methylation map is then generated. When applied to H3 lysine 27 trimethylation (H3K27me3), we found that H3K27me3 and DNA methylation are compatible throughout most of the genome, except for CpG islands, where these two marks are mutually exclusive. Further ChIP-BS-seq-based analysis in Dnmt triple-knockout (TKO) embryonic stem cells revealed that total loss of CpG methylation is associated with alteration of H3K27me3 levels throughout the genome: H3K27me3 in localized peaks is decreased while broad local enrichments (BLOCs) of H3K27me3 are formed. At an even broader scale, these BLOCs correspond to regions of high DNA methylation in wild-type ES cells, suggesting that DNA methylation prevents H3K27me3 deposition locally and at a megabase scale. Our strategy provides a unique way of investigating global interdependencies between DNA methylation and other chromatin features.

Whole-Genome Bisulfite Sequencing of Two Distinct Interconvertible DNA Methylomes of Mouse Embryonic Stem Cells

The use of two kinase inhibitors (2i) enables derivation of mouse embryonic stem cells (ESCs) in the pluripotent ground state. Using whole-genome bisulfite sequencing (WGBS), we show that male 2i ESCs are globally hypomethylated compared to conventional ESCs maintained in serum. In serum, female ESCs are hypomethyated similarly to male ESCs in 2i, and DNA methylation is further reduced in 2i. Regions with elevated DNA methylation in 2i strongly correlate with the presence of H3K9me3 on endogenous retroviruses (ERVs) and imprinted loci. The methylome of male ESCs in serum parallels postimplantation blastocyst cells, while 2i stalls ESCs in a hypomethylated, ICM-like state. WGBS analysis during adaptation of 2i ESCs to serum suggests that deposition of DNA methylation is largely random, while loss of DNA methylation during reversion to 2i occurs passively, initiating at TET1 binding sites. Together, our analysis provides insight into DNA methylation dynamics in cultured ESCs paralleling early developmental processes.

The Lrp family of transcriptional regulators

Genome analysis has revealed that members of the Lrp family of transcriptional regulators are widely distributed among prokaryotes, both bacteria and archaea. The archetype Leucine‐responsive Regulatory Protein from Escherichia coli is a global regulator involved in modulating a variety of metabolic functions, including the catabolism and anabolism of amino acids as well as pili synthesis. Most Lrp homologues, however, appear to act as specific regulators of amino acid metabolism‐related genes. Like most prokaryotic transcriptional regulators, Lrp‐like regulators consist of a DNA‐binding domain and a ligand‐binding domain. The crystal structure of the Pyrococcus furiosus LrpA revealed an N‐terminal domain with a common helix–turn–helix fold, and a C‐terminal domain with a typical αβ‐sandwich fold. The latter regulatory domain constitutes a novel ligand‐binding site and has been designated RAM. Database analysis reveals that the RAM domain is present in many prokaryotic genomes, potentially encoding (1) Lrp‐homologues, when fused to a DNA‐binding domain (2) enzymes, when fused as a potential regulatory domain to a catalytic domain, and (3) stand‐alone RAM modules with unknown function. The architecture of Lrp regulators with two distinct domains that harbour the regulatory (effector‐binding) site and the active (DNA‐binding) site, and their separation by a flexible hinge region, suggests a general allosteric switch of Lrp‐like regulators.

PML-RARα/RXR Alters the Epigenetic Landscape in Acute Promyelocytic Leukemia

Many different molecular mechanisms have been associated with PML-RARalpha-dependent transformation of hematopoietic progenitors. Here, we identified high confidence PML-RARalpha binding sites in an acute promyelocytic leukemia (APL) cell line and in two APL primary blasts. We found colocalization of PML-RARalpha with RXR to the vast majority of these binding regions. Genome-wide epigenetic studies revealed that treatment with pharmacological doses of all-trans retinoic acid induces changes in H3 acetylation, but not H3K27me3, H3K9me3, or DNA methylation at the PML-RARalpha/RXR binding sites or at nearby target genes. Our results suggest that PML-RARalpha/RXR functions as a local chromatin modulator and that specific recruitment of histone deacetylase activities to genes important for hematopoietic differentiation, RAR signaling, and epigenetic control is crucial to its transforming potential.

Whole-genome DNA methylation profiling using MethylCap-seq

MethylCap-seq is a robust procedure for genome-wide profiling of DNA methylation. The approach consists of the capture of methylated DNA using the MBD domain of MeCP2, and subsequent next-generation sequencing of eluted DNA. Elution of the captured methylated DNA is done in steps using a salt gradient, which stratifies the genome into fractions with different CpG density. The enrichment reached within the individual eluates allows for cost-effective deep sequence coverage. The profiles together yield a detailed genome-wide map of methylated regions and readily allows detection of DNA methylation in known and novel regions. Here, we describe principles and details of the MethylCap-seq procedure using different sources of starting material.

Crystal structure of the Lrp‐like transcriptional regulator from the archaeon Pyrococcus furiosus

The LrpA protein from the hyperthermophilic archaeon Pyrococcus furiosus belongs to the Lrp/AsnC family of transcriptional regulatory proteins, of which the Escherichia coli leucine‐responsive regulatory protein is the archetype. Its crystal structure has been determined at 2.9 Å resolution and is the first for a member of the Lrp/AsnC family, as well as one of the first for a transcriptional regulator from a hyperthermophile. The structure consists of an N‐terminal domain containing a helix–turn–helix (HtH) DNA‐binding motif, and a C‐terminal domain of mixed α/β character reminiscent of a number of RNA‐ and DNA‐binding domains. Pyrococcus furiosus LrpA forms a homodimer mainly through interactions between the antiparallel β‐sheets of the C‐terminal domain, and further interactions lead to octamer formation. The LrpA structure suggests how the protein might bind and possibly distort its DNA substrate through use of its HtH motifs and control gene expression. A possible location for an effector binding site is proposed by using sequence comparisons with other members of the family coupled to mutational analysis.