Search

Chromatin profiling at locus resolution uncovers gene regulatory features that define cell types and developmental trajectories, but it remains challenging to map and compare different chromatin-associated proteins in the same sample. Here we describe Multiple Target Identification by Tagmentation (MulTI-Tag), an antibody barcoding approach for profiling multiple chromatin features simultaneously in single cells. We optimized MulTI-Tag to retain high sensitivity and specificity, and we demonstrate detection of up to three histone modifications in the same cell: H3K27me3, H3K4me1/2 and H3K36me3. We apply MulTI-Tag to resolve distinct cell types and developmental trajectories; to distinguish unique, coordinated patterns of active and repressive element regulatory usage associated with differentiation outcomes; and to uncover associations between histone marks. Multifactorial epigenetic profiling holds promise for comprehensively characterizing cell-specific gene regulatory landscapes in development and disease.

Feature-region combination AUCs and benchmarking Sheet 1: Log2 fold-change, p-value, and q-value between ARPC and NEPC lines for central mean coverage in all queried (338) TFs (two tailed Mann-Whitney U test, Benjamini-Hochberg adjusted). Sheet 2: 100-fold cross-validation AUCs for all region and feature combinations subset by the ‘AR10’ overlapping features (see Methods). Sheet 3: 100-fold cross-validation AUCs for all region and feature combinations subset by the ‘Phenotype 47-defining’ overlapping features. Sheet 4: 100-fold cross-validation AUCs for all region and feature combinations (global). Sheet 5: Predictions scores, tumor fraction, depth of coverage, and subtype for benchmarking admixtures. Sheet 6: AUCs for unsupervised prediction of admixture subtypes grouped by tumor fraction and depth. Sheet 7: NEPC:ARPC ratio, tumor fraction, and ARPC and NEPC fraction predictions for mixed phenotype admixtures using Keraon (see Methods).

Supplementary Figures S1-S18 and Supplementary Table legends.

Polycomb complexes regulate cell type–specific gene expression programs through heritable silencing of target genes. Trimethylation of histone H3 lysine 27 (H3K27me3) is essential for this process. Perturbation of H3K36 is thought to interfere with H3K27me3. We show that mutants of Drosophila replication-dependent ( H3.2 K36R ) or replication-independent ( H3.3 K36R ) histone H3 genes generally maintain Polycomb silencing and reach later stages of development. In contrast, combined ( H3.3 K36R H3.2 K36R ) mutants display widespread Hox gene misexpression and fail to develop past the first larval stage. Chromatin profiling revealed that the H3.2 K36R mutation disrupts H3K27me3 levels broadly throughout silenced domains, whereas these regions are mostly unaffected in H3.3 K36R animals. Analysis of H3.3 distributions showed that this histone is enriched at presumptive Polycomb response elements located outside of silenced domains but relatively depleted from those inside. We conclude that H3.2 and H3.3 K36 residues collaborate to repress Hox genes using different mechanisms.

Advanced prostate cancers comprise distinct phenotypes, but tumor classification remains clinically challenging. Here, we harnessed circulating tumor DNA (ctDNA) to study tumor phenotypes by ascertaining nucleosome positioning patterns associated with transcription regulation. We sequenced plasma ctDNA whole genomes from patient-derived xenografts representing a spectrum of androgen receptor active (ARPC) and neuroendocrine (NEPC) prostate cancers. Nucleosome patterns associated with transcriptional activity were reflected in ctDNA at regions of genes, promoters, histone modifications, transcription factor binding, and accessible chromatin. We identified the activity of key phenotype-defining transcriptional regulators from ctDNA, including AR, ASCL1, HOXB13, HNF4G, and GATA2. To distinguish NEPC and ARPC in patient plasma samples, we developed prediction models that achieved accuracies of 97% for dominant phenotypes and 87% for mixed clinical phenotypes. Although phenotype classification is typically assessed by IHC or transcriptome profiling from tumor biopsies, we demonstrate that ctDNA provides comparable results with diagnostic advantages for precision oncology. Significance: This study provides insights into the dynamics of nucleosome positioning and gene regulation associated with cancer phenotypes that can be ascertained from ctDNA. New methods for classification in phenotype mixtures extend the utility of ctDNA beyond assessments of somatic DNA alterations with important implications for molecular classification and precision oncology. This article is highlighted in the In This Issue feature, p. 517

Polycomb complexes regulate cell-type specific gene expression programs through heritable silencing of target genes. Trimethylation of histone H3 lysine 27 (H3K27me3) is essential for this process. Perturbation of H3K36 is thought to interfere with H3K27me3. We show that mutants ofDrosophilareplication-dependent(H3.2K36R)or -independent(H3.3K36R)histone H3 genes generally maintain Polycomb silencing and reach later stages of development. In contrast, combined(H3.3K36RH3.2K36R)mutants display widespread Hox gene misexpression and fail to develop past the first larval stage. Chromatin profiling revealed that theH3.2K36Rmutation disrupts H3K27me3 levels broadly throughout silenced domains, whereas these regions are mostly unaffected inH3.3K36Ranimals. Analysis of H3.3 distributions showed that this histone is enriched at presumptive PREs (Polycomb Response Elements) located outside of silenced domains but relatively depleted from those inside. We conclude that H3.2 and H3.3 K36 residues collaborate to repress Hox genes using different mechanisms.Short summaryHistone H3.2 and H3.3 K36 residues ensure Hox gene silencing and enable development by different, but synergistic mechanisms.

SummaryCleavage Under Targets & Tagmentation (CUT&Tag) is an antibody-directed transposase tethering strategy for in situ chromatin profiling in small samples and single cells. We describe a modified CUT&Tag protocol using a mixture of an antibody to the initiation form of RNA Polymerase II (Pol2 Serine-5 phosphate) and an antibody to repressive Polycomb domains (H3K27me3) followed by computational signal deconvolution to produce high-resolution maps of both the active and repressive regulomes in single cells. The ability to seamlessly map active promoters, enhancers and repressive regulatory elements using a single workflow provides a complete regulome profiling strategy suitable for high-throughput single-cell platforms.

Cleavage Under Targets and Tagmentation (CUT&Tag) is an antibody-directed transposase tethering strategy for in situ chromatin profiling in small samples and single cells. We describe a modified CUT&Tag protocol using a mixture of an antibody to the initiation form of RNA polymerase II (Pol2 Serine-5 phosphate) and an antibody to repressive Polycomb domains (H3K27me3) followed by computational signal deconvolution to produce high-resolution maps of both the active and repressive regulomes in single cells. The ability to seamlessly map active promoters, enhancers, and repressive regulatory elements using a single workflow provides a complete regulome profiling strategy suitable for high-throughput single-cell platforms.

Chromatin profiling at locus resolution uncovers gene regulatory features that define cell types and developmental trajectories, but it remains challenging to map and compare different chromatin-associated proteins in the same sample. Here we describe Multiple Target Identification by Tagmentation (MulTI-Tag), an antibody barcoding approach for profiling multiple chromatin features simultaneously in single cells. We optimized MulTI-Tag to retain high sensitivity and specificity, and we demonstrate detection of up to three histone modifications in the same cell: H3K27me3, H3K4me1/2 and H3K36me3. We apply MulTI-Tag to resolve distinct cell types and developmental trajectories; to distinguish unique, coordinated patterns of active and repressive element regulatory usage associated with differentiation outcomes; and to uncover associations between histone marks. Multifactorial epigenetic profiling holds promise for comprehensively characterizing cell-specific gene regulatory landscapes in development and disease.

Background Our understanding of eukaryotic gene regulation is limited by the complexity of protein–DNA interactions that comprise the chromatin landscape and by inefficient methods for characterizing these interactions. We recently introduced CUT&RUN, an antibody-targeted nuclease cleavage method that profiles DNA-binding proteins, histones and chromatin-modifying proteins in situ with exceptional sensitivity and resolution. Results Here, we describe an automated CUT&RUN platform and apply it to characterize the chromatin landscapes of human cells. We find that automated CUT&RUN profiles of histone modifications crisply demarcate active and repressed chromatin regions, and we develop a continuous metric to identify cell-type-specific promoter and enhancer activities. We test the ability of automated CUT&RUN to profile frozen tumor samples and find that our method readily distinguishes two pediatric glioma xenografts by their subtype-specific gene expression programs. Conclusions The easy, cost-effective workflow makes automated CUT&RUN an attractive tool for high-throughput characterization of cell types and patient samples. Electronic supplementary material The online version of this article (10.1186/s13072-018-0243-8) contains supplementary material, which is available to authorized users.

Prostate cancer (PC) is driven by androgen receptor (AR) activity, a master regulator of prostate development and homeostasis. Frontline therapies for metastatic PC deprive the AR of the activating ligands testosterone (T) and dihydrotestosterone (DHT) by limiting their biosynthesis or blocking AR binding. Notably, AR signaling is dichotomous, inducing growth at lower activity levels, while suppressing growth at higher levels. Recent clinical studies have exploited this effect by administration of supraphysiological concentrations of T, resulting in clinical responses and improvements in quality of life. However, the use of T as a therapeutic agent in oncology is limited by poor drug-like properties as well as rapid and variable metabolism. Here, we investigated the antitumor effects of selective AR modulators (SARMs), which are small-molecule nonsteroidal AR agonists developed to treat muscle wasting and cachexia. Several orally administered SARMs activated the AR program in PC models. AR cistromes regulated by steroidal androgens and SARMs were superimposable. Coregulatory proteins including HOXB13 and GRHL2 comprised AR complexes assembled by both androgens and SARMs. At bioavailable concentrations, SARMs repressed MYC oncoprotein expression and inhibited the growth of castration-sensitive and castration-resistant PC in vitro and in vivo. These results support further clinical investigation of SARMs for treating advanced PC.

Acute myeloid and lymphoid leukemias often harbor chromosomal translocations involving the Mixed Lineage Leukemia-1 gene, encoding the KMT2A lysine methyltransferase. The most common translocations produce in-frame fusions of KMT2A to other chromatin regulatory proteins. Here we develop a strategy to map the genome-wide occupancy of oncogenic KMT2A fusion proteins in primary patient samples regardless of fusion partner. By modifying the versatile CUT&Tag method for full automation we identify common and tumor-specific patterns of aberrant chromatin regulation induced by different KMT2A fusion proteins. Integration of automated and single-cell CUT&Tag uncovers epigenomic heterogeneity within patient samples and predicts sensitivity to therapeutic agents.

Lysine 27-to-methionine (K27M) mutations in the H3.1 or H3.3 histone genes are characteristic of pediatric diffuse midline gliomas (DMGs). These oncohistone mutations dominantly inhibit histone H3K27 trimethylation and silencing, but it is unknown how oncohistone type affects gliomagenesis. We show that the genomic distributions of H3.1 and H3.3 oncohistones in human patient-derived DMG cells are consistent with the DNA replication-coupled deposition of histone H3.1 and the predominant replication-independent deposition of histone H3.3. Although H3K27 trimethylation is reduced for both oncohistone types, H3.3K27M-bearing cells retain some domains, and only H3.1K27M-bearing cells lack H3K27 trimethylation. Neither oncohistone interferes with PRC2 binding. Using Drosophila as a model, we demonstrate that inhibition of H3K27 trimethylation occurs only when H3K27M oncohistones are deposited into chromatin and only when expressed in cycling cells. We propose that oncohistones inhibit the H3K27 methyltransferase as chromatin patterns are being duplicated in proliferating cells, predisposing them to tumorigenesis.

Aggressive cancers often have activating mutations in growth-controlling oncogenes and inactivating mutations in tumor-suppressor genes. In neuroblastoma, amplification of the MYCN oncogene and inactivation of the ATRX tumor-suppressor gene correlate with high-risk disease and poor prognosis. Here we show that ATRX mutations and MYCN amplification are mutually exclusive across all ages and stages in neuroblastoma. Using human cell lines and mouse models, we found that elevated MYCN expression and ATRX mutations are incompatible. Elevated MYCN levels promote metabolic reprogramming, mitochondrial dysfunction, reactive-oxygen species generation, and DNA-replicative stress. The combination of replicative stress caused by defects in the ATRX–histone chaperone complex, and that induced by MYCN-mediated metabolic reprogramming, leads to synthetic lethality. Therefore, ATRX and MYCN represent an unusual example, where inactivation of a tumor-suppressor gene and activation of an oncogene are incompatible. This synthetic lethality may eventually be exploited to improve outcomes for patients with high-risk neuroblastoma., In most cancers, mutations that lead to oncogene activation and tumor suppressor inactivation synergize to promote tumorigenesis. However, in neuroblastomas, MYCN amplification and ATRX mutations are mutually exclusive and incompatible.

Background CUT&RUN is an efficient epigenome profiling method that identifies sites of DNA binding protein enrichment genome-wide with high signal to noise and low sequencing requirements. Currently, the analysis of CUT&RUN data is complicated by its exceptionally low background, which renders programs designed for analysis of ChIP-seq data vulnerable to oversensitivity in identifying sites of protein binding. Results Here we introduce Sparse Enrichment Analysis for CUT&RUN (SEACR), an analysis strategy that uses the global distribution of background signal to calibrate a simple threshold for peak calling. SEACR discriminates between true and false-positive peaks with near-perfect specificity from “gold standard” CUT&RUN datasets and efficiently identifies enriched regions for several different protein targets. We also introduce a web server (http://seacr.fredhutch.org) for plug-and-play analysis with SEACR that facilitates maximum accessibility across users of all skill levels. Conclusions SEACR is a highly selective peak caller that definitively validates the accuracy of CUT&RUN for datasets with known true negatives. Its ease of use and performance in comparison with existing peak calling strategies make it an ideal choice for analyzing CUT&RUN data. Electronic supplementary material The online version of this article (10.1186/s13072-019-0287-4) contains supplementary material, which is available to authorized users.

Summary Although the in vitro structural and in vivo spatial characteristics of transcription factor (TF) binding are well defined, TF interactions with chromatin and other companion TFs during development are poorly understood. To analyze such interactions in vivo, we profiled several TFs across a time course of human embryonic stem cell differentiation and studied their interactions with nucleosomes and co-occurring TFs by enhanced chromatin occupancy (EChO), a computational strategy for classifying TF interactions with chromatin. EChO shows that multiple individual TFs can employ either direct DNA binding or “pioneer” nucleosome binding at different enhancer targets. Nucleosome binding is not exclusively confined to inaccessible chromatin but rather correlated with local binding of other TFs and degeneracy at key bases in the pioneer factor target motif responsible for direct DNA binding. Our strategy reveals a dynamic exchange of TFs at enhancers across developmental time that is aided by pioneer nucleosome binding.

Sophisticated gene-regulatory mechanisms probably evolved in prokaryotes billions of years before the emergence of modern eukaryotes, which inherited the same basic enzymatic machineries. However, the epigenomic landscapes of eukaryotes are dominated by nucleosomes, which have acquired roles in genome packaging, mitotic condensation and silencing parasitic genomic elements. Although the molecular mechanisms by which nucleosomes are displaced and modified have been described, just how transcription factors, histone variants and modifications and chromatin regulators act on nucleosomes to regulate transcription is the subject of considerable ongoing study. We explore the extent to which these transcriptional regulatory components function in the context of the evolutionarily ancient role of chromatin as a barrier to processes acting on DNA and how chromatin proteins have diversified to carry out evolutionarily recent functions that accompanied the emergence of differentiation and development in multicellular eukaryotes.

We previously described a novel alternative to Chromatin Immunoprecipitation, Cleavage Under Targets & Release Using Nuclease (CUT&RUN), in which unfixed permeabilized cells are incubated with antibody, followed by binding of a Protein A-Micrococcal Nuclease (pA/MNase) fusion protein (1). Upon activation of tethered MNase, the bound complex is excised and released into the supernatant for DNA extraction and sequencing. Here we introduce four enhancements to CUT&RUN: 1) a hybrid Protein A-Protein G-MNase construct that expands antibody compatibility and simplifies purification; 2) a modified digestion protocol that inhibits premature release of the nuclease-bound complex; 3) a calibration strategy based on carry-over ofE. coliDNA introduced with the fusion protein; and 4) a novel peak-calling strategy customized for the low-background profiles obtained using CUT&RUN. These new features, coupled with the previously described low-cost, high efficiency, high reproducibility and high-throughput capability of CUT&RUN make it the method of choice for routine epigenomic profiling.

Previously, we described a novel alternative to chromatin immunoprecipitation, CUT&RUN, in which unfixed permeabilized cells are incubated with antibody, followed by binding of a protein A-Micrococcal Nuclease (pA/MNase) fusion protein (Skene and Henikoff, 2017). Here we introduce three enhancements to CUT&RUN: A hybrid protein A-Protein G-MNase construct that expands antibody compatibility and simplifies purification, a modified digestion protocol that inhibits premature release of the nuclease-bound complex, and a calibration strategy based on carry-over of E. coli DNA introduced with the fusion protein. These new features, coupled with the previously described low-cost, high efficiency, high reproducibility and high-throughput capability of CUT&RUN make it the method of choice for routine epigenomic profiling.

Histone post-translational modifications (PTMs) are thought to participate in a range of essential molecular and cellular processes, including gene expression, replication, and nuclear organization. Importantly, histone PTMs are also thought to be prime candidates for carriers of epigenetic information across cell cycles and generations. However, directly testing the necessity of histone PTMs themselves in these processes by mutagenesis has been extremely difficult to carry out because of the highly repetitive nature of histone genes in animal genomes. We developed a transgenic system to generate Drosophila melanogaster genotypes in which the entire complement of replication-dependent histone genes is mutant at a residue of interest. We built a BAC vector containing a visible marker for lineage tracking along with the capacity to clone large (60-100 kb) inserts that subsequently can be site-specifically integrated into the D. melanogaster genome. We demonstrate that artificial tandem arrays of the core 5 kb replication-dependent histone repeat can be generated with relative ease. This genetic platform represents the first histone replacement system to leverage a single tandem transgenic insertion for facile genetics and analysis of molecular and cellular phenotypes. We demonstrate the utility of our system for directly preventing histone residues from being modified, and studying the consequent phenotypes. This system can be generalized to the cloning and transgenic insertion of any tandemly repeated sequence of biological interest.

Histone post-translational modifications (PTMs) are thought to participate in a range of essential molecular and cellular processes, including gene expression, replication, and nuclear organization. Importantly, histone PTMs are also thought to be prime candidates for carriers of epigenetic information across cell cycles and generations. However, directly testing the necessity of histone PTMs themselves in these processes by mutagenesis has been extremely difficult to carry out because of the highly repetitive nature of histone genes in animal genomes. We developed a transgenic system to generate Drosophila melanogaster genotypes in which the entire complement of replication-dependent histone genes is mutant at a residue of interest. We built a BAC vector containing a visible marker for lineage tracking along with the capacity to clone large (60-100 kb) inserts that subsequently can be site-specifically integrated into the D. melanogaster genome. We demonstrate that artificial tandem arrays of the core 5 kb replication-dependent histone repeat can be generated with relative ease. This genetic platform represents the first histone replacement system to leverage a single tandem transgenic insertion for facile genetics and analysis of molecular and cellular phenotypes. We demonstrate the utility of our system for directly preventing histone residues from being modified, and studying the consequent phenotypes. This system can be generalized to the cloning and transgenic insertion of any tandemly repeated sequence of biological interest.

BackgroundHigh-resolution transcription start site (TSS) mapping in D. melanogaster embryos and cell lines has revealed a rich and detailed landscape of both cis- and trans-regulatory elements and factors. However, TSS profiling has not been investigated in an orthogonal in vivo setting. Here, we present a comprehensive dataset that links TSS dynamics with nucleosome occupancy and gene expression in the wandering third instar larva, a developmental stage characterized by large-scale shifts in transcriptional programs in preparation for metamorphosis.ResultsThe data recapitulate major regulatory classes of TSSs, based on peak width, promoter-proximal polymerase pausing, and cis-regulatory element density. We confirm the paucity of divergent transcription units in D. melanogaster, but also identify notable exceptions. Furthermore, we identify thousands of novel initiation events occurring at unannotated TSSs that can be classified into functional categories by their local density of histone modifications. Interestingly, a sub-class of these unannotated TSSs overlaps with functionally validated enhancer elements, consistent with a regulatory role for “enhancer RNAs” (eRNAs) in defining transcriptional programs that are important for animal development.ConclusionsHigh-depth TSS mapping is a powerful strategy for identifying and characterizing low-abundance and/or low-stability RNAs. Global analysis of transcription initiation patterns in a developing organism reveals a vast number of novel initiation events that identify potential eRNAs as well as other non-coding transcripts critical for animal development.

Histone H3 lysine 36 methylation (H3K36me) is thought to participate in a host of co-transcriptional regulatory events. To study the function of this residue independent from the enzymes that modify it, we used a ‘histone replacement’ system in Drosophila to generate a non-modifiable H3K36 lysine-to-arginine (H3K36R) mutant. We observed global dysregulation of mRNA levels in H3K36R animals that correlates with the incidence of H3K36me3. Similar to previous studies, we found that mutation of H3K36 also resulted in H4 hyperacetylation. However, neither cryptic transcription initiation, nor alternative pre-mRNA splicing, contributed to the observed changes in expression, in contrast with previously reported roles for H3K36me. Interestingly, knockdown of the RNA surveillance nuclease, Xrn1, and members of the CCR4-Not deadenylase complex, restored mRNA levels for a class of downregulated, H3K36me3-rich genes. We propose a post-transcriptional role for modification of replication-dependent H3K36 in the control of metazoan gene expression. DOI: http://dx.doi.org/10.7554/eLife.23249.001, eLife digest In a single human cell there is enough DNA to stretch over a meter if laid end to end. To fit this DNA inside the cell – which is less than 20 micrometers in diameter – the DNA is tightly wrapped around millions of proteins known as histones, which look like “beads” along a “string” of DNA. These histones can prevent other proteins from binding to DNA and activating specific genes. Therefore, cells use enzymes to chemically modify histones to allow particular stretches of DNA to be unwrapped at specific times. Proteins are made up of building blocks called amino acids. A specific amino acid on histones known as H3K36 is modified in certain sections of DNA that suggest it affects the activities of many genes. However, the precise role of this amino acid remains unclear. Previous studies have tried to investigate this by removing the enzymes that modify it, but these enzymes can also modify many other proteins, making it difficult to know what exactly causes the changes in gene activity. Fruit flies are often used in experiments as models of how genetic processes work in humans and other animals. Like us, fruit flies also package their DNA using histones. To investigate the role of H3K36, Meers et al. generated a mutant fruit fly that has a version of the amino acid that cannot be chemically modified by the normal enzymes. Unexpectedly, the experiments suggest that some changes in gene activity that have been previously reported to be caused by modifying H3K36 might actually be due to other factors. Meers et al. found that H3K36 modifications may instead “mark” certain genes to be more active than they otherwise would be. These findings provide a starting point for understanding exactly how H3K36 regulates gene activity. The next challenge is to refine our understanding of how H3K36 modification affects genes in cancer and other diseases, which may aid the development of new therapies to treat these conditions. DOI: http://dx.doi.org/10.7554/eLife.23249.002

Reduced levels of survival motor neuron (SMN) protein lead to a neuromuscular disease called spinal muscular atrophy (SMA). Animal models of SMA recapitulate many aspects of the human disease, including locomotion and viability defects, but have thus far failed to uncover the causative link between a lack of SMN protein and neuromuscular dysfunction. While SMN is known to assemble small nuclear ribonucleoproteins (snRNPs) that catalyze pre-mRNA splicing, it remains unclear whether disruptions in splicing are etiologic for SMA. To investigate this issue, we carried out RNA deep-sequencing (RNA-seq) on age-matched Drosophila Smn-null and wild-type larvae. Comparison of genome-wide mRNA expression profiles with publicly available data sets revealed the timing of a developmental arrest in the Smn mutants. Furthermore, genome-wide differences in splicing between wild-type and Smn animals did not correlate with changes in mRNA levels. Specifically, we found that mRNA levels of genes that contain minor introns vary more over developmental time than they do between wild-type and Smn mutants. An analysis of reads mapping to minor-class intron–exon junctions revealed only small changes in the splicing of minor introns in Smn larvae, within the normal fluctuations that occur throughout development. In contrast, Smn mutants displayed a prominent increase in levels of stress-responsive transcripts, indicating a systemic response to the developmental arrest induced by loss of SMN protein. These findings not only provide important mechanistic insight into the developmental arrest displayed by Smn mutants, but also argue against a minor-intron-dependent etiology for SMA.

SummaryHistones and their posttranslational modifications influence the regulation of many DNA-dependent processes. Although an essential role for histone-modifying enzymes in these processes is well established, defining the specific contribution of individual histone residues remains a challenge because many histone-modifying enzymes have nonhistone targets. This challenge is exacerbated by the paucity of suitable approaches to genetically engineer histone genes in metazoans. Here, we describe a platform in Drosophila for generating and analyzing any desired histone genotype, and we use it to test the in vivo function of three histone residues. We demonstrate that H4K20 is neither essential for DNA replication nor for completion of development, unlike inferences drawn from analyses of H4K20 methyltransferases. We also show that H3K36 is required for viability and H3K27 is essential for maintenance of cellular identity but not for gene activation. These findings highlight the power of engineering histones to interrogate genome structure and function in animals.