Christoph Bock

Somatic Mutations of Calreticulin in Myeloproliferative Neoplasms

BACKGROUND Approximately 50 to 60% of patients with essential thrombocythemia or primary myelofibrosis carry a mutation in the Janus kinase 2 gene (JAK2), and an additional 5 to 10% have activating mutations in the thrombopoietin receptor gene (MPL). So far, no specific molecular marker has been identified in the remaining 30 to 45% of patients. METHODS We performed whole-exome sequencing to identify somatically acquired mutations in six patients who had primary myelofibrosis without mutations in JAK2 or MPL. Resequencing of CALR, encoding calreticulin, was then performed in cohorts of patients with myeloid neoplasms. RESULTS Somatic insertions or deletions in exon 9 of CALR were detected in all patients who underwent whole-exome sequencing. Resequencing in 1107 samples from patients with myeloproliferative neoplasms showed that CALR mutations were absent in polycythemia vera. In essential thrombocythemia and primary myelofibrosis, CALR mutations and JAK2 and MPL mutations were mutually exclusive. Among patients with essential thrombocythemia or primary myelofibrosis with nonmutated JAK2 or MPL, CALR mutations were detected in 67% of those with essential thrombocythemia and 88% of those with primary myelofibrosis. A total of 36 types of insertions or deletions were identified that all cause a frameshift to the same alternative reading frame and generate a novel C-terminal peptide in the mutant calreticulin. Overexpression of the most frequent CALR deletion caused cytokine-independent growth in vitro owing to the activation of signal transducer and activator of transcription 5 (STAT5) by means of an unknown mechanism. Patients with mutated CALR had a lower risk of thrombosis and longer overall survival than patients with mutated JAK2. CONCLUSIONS Most patients with essential thrombocythemia or primary myelofibrosis that was not associated with a JAK2 or MPL alteration carried a somatic mutation in CALR. The clinical course in these patients was more indolent than that in patients with the JAK2 V617F mutation. (Funded by the MPN Research Foundation and Associazione Italiana per la Ricerca sul Cancro.).

Reference Maps of Human ES and iPS Cell Variation Enable High-Throughput Characterization of Pluripotent Cell Lines

The developmental potential of human pluripotent stem cells suggests that they can produce disease-relevant cell types for biomedical research. However, substantial variation has been reported among pluripotent cell lines, which could affect their utility and clinical safety. Such cell-line-specific differences must be better understood before one can confidently use embryonic stem (ES) or induced pluripotent stem (iPS) cells in translational research. Toward this goal we have established genome-wide reference maps of DNA methylation and gene expression for 20 previously derived human ES lines and 12 human iPS cell lines, and we have measured the in vitro differentiation propensity of these cell lines. This resource enabled us to assess the epigenetic and transcriptional similarity of ES and iPS cells and to predict the differentiation efficiency of individual cell lines. The combination of assays yields a scorecard for quick and comprehensive characterization of pluripotent cell lines.

Dnmt3a is essential for hematopoietic stem cell differentiation

Loss of the de novo DNA methyltransferases Dnmt3a and Dnmt3b in embryonic stem cells obstructs differentiation; however, the role of these enzymes in somatic stem cells is largely unknown. Using conditional ablation, we show that Dnmt3a loss progressively impairs hematopoietic stem cell (HSC) differentiation over serial transplantation, while simultaneously expanding HSC numbers in the bone marrow. Dnmt3a-null HSCs show both increased and decreased methylation at distinct loci, including substantial CpG island hypermethylation. Dnmt3a-null HSCs upregulate HSC multipotency genes and downregulate differentiation factors, and their progeny exhibit global hypomethylation and incomplete repression of HSC-specific genes. These data establish Dnmt3a as a critical participant in the epigenetic silencing of HSC regulatory genes, thereby enabling efficient differentiation.

Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.

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.

Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling

Genome-wide mapping of 5-methylcytosine is of broad interest to many fields of biology and medicine. A variety of methods have been developed, and several have recently been advanced to genome-wide scale using arrays and next-generation sequencing approaches. We have previously reported reduced representation bisulfite sequencing (RRBS), a bisulfite-based protocol that enriches CG-rich parts of the genome, thereby reducing the amount of sequencing required while capturing the majority of promoters and other relevant genomic regions. The approach provides single-nucleotide resolution, is highly sensitive and provides quantitative DNA methylation measurements. This protocol should enable any standard molecular biology laboratory to generate RRBS libraries of high quality. Briefly, purified genomic DNA is digested by the methylation-insensitive restriction enzyme MspI to generate short fragments that contain CpG dinucleotides at the ends. After end-repair, A-tailing and ligation to methylated Illumina adapters, the CpG-rich DNA fragments (40–220 bp) are size selected, subjected to bisulfite conversion, PCR amplified and end sequenced on an Illumina Genome Analyzer. Note that alignment and analysis of RRBS sequencing reads are not covered in this protocol. The extremely low input requirements (10–300 ng), the applicability of the protocol to formalin-fixed and paraffin-embedded samples, and the techniques single-nucleotide resolution extends RRBS to a wide range of biological and clinical samples and research applications. The entire process of RRBS library construction takes ∼9 d.

Analysing and interpreting DNA methylation data

DNA methylation is an epigenetic mark that has suspected regulatory roles in a broad range of biological processes and diseases. The technology is now available for studying DNA methylation genome-wide, at a high resolution and in a large number of samples. This Review discusses relevant concepts, computational methods and software tools for analysing and interpreting DNA methylation data. It focuses not only on the bioinformatic challenges of large epigenome-mapping projects and epigenome-wide association studies but also highlights software tools that make genome-wide DNA methylation mapping more accessible for laboratories with limited bioinformatics experience.

Genomic Distribution and Inter-Sample Variation of Non- CpG Methylation across Human Cell Types

DNA methylation plays an important role in development and disease. The primary sites of DNA methylation in vertebrates are cytosines in the CpG dinucleotide context, which account for roughly three quarters of the total DNA methylation content in human and mouse cells. While the genomic distribution, inter-individual stability, and functional role of CpG methylation are reasonably well understood, little is known about DNA methylation targeting CpA, CpT, and CpC (non-CpG) dinucleotides. Here we report a comprehensive analysis of non-CpG methylation in 76 genome-scale DNA methylation maps across pluripotent and differentiated human cell types. We confirm non-CpG methylation to be predominantly present in pluripotent cell types and observe a decrease upon differentiation and near complete absence in various somatic cell types. Although no function has been assigned to it in pluripotency, our data highlight that non-CpG methylation patterns reappear upon iPS cell reprogramming. Intriguingly, the patterns are highly variable and show little conservation between different pluripotent cell lines. We find a strong correlation of non-CpG methylation and DNMT3 expression levels while showing statistical independence of non-CpG methylation from pluripotency associated gene expression. In line with these findings, we show that knockdown of DNMTA and DNMT3B in hESCs results in a global reduction of non-CpG methylation. Finally, non-CpG methylation appears to be spatially correlated with CpG methylation. In summary these results contribute further to our understanding of cytosine methylation patterns in human cells using a large representative sample set.

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.

Proliferation-Dependent Alterations of the DNA Methylation Landscape Underlie Hematopoietic Stem Cell Aging

The functional potential of hematopoietic stem cells (HSCs) declines during aging, and in doing so, significantly contributes to hematopoietic pathophysiology in the elderly. To explore the relationship between age-associated HSC decline and the epigenome, we examined global DNA methylation of HSCs during ontogeny in combination with functional analysis. Although the DNA methylome is generally stable during aging, site-specific alterations of DNA methylation occur at genomic regions associated with hematopoietic lineage potential and selectively target genes expressed in downstream progenitor and effector cells. We found that age-associated HSC decline, replicative limits, and DNA methylation are largely dependent on the proliferative history of HSCs, yet appear to be telomere-length independent. Physiological aging and experimentally enforced proliferation of HSCs both led to DNA hypermethylation of genes regulated by Polycomb Repressive Complex 2. Our results provide evidence that epigenomic alterations of the DNA methylation landscape contribute to the functional decline of HSCs during aging.