Multifactorial profiling of epigenetic landscapes at single-cell resolution using MulTI-Tag

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.

Distinct roles for canonical and variant histone H3 lysine-36 in Polycomb silencing

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.2K36R) or replication-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 the H3.2K36R mutation disrupts H3K27me3 levels broadly throughout silenced domains, whereas these regions are mostly unaffected in H3.3K36R 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.

Nucleosome Patterns in Circulating Tumor DNA Reveal Transcriptional Regulation of Advanced Prostate Cancer Phenotypes

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.

CUT&Tag2for1: a modified method for simultaneous profiling of the accessible and silenced regulome in single cells

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.

Histone deposition pathways determine the chromatin landscapes of H3.1 and H3.3 K27M oncohistones

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 DNAreplication-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.

MYCN amplification and ATRX mutations are incompatible in neuroblastoma

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.

Pioneer Factor-Nucleosome Binding Events during Differentiation Are Motif Encoded

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.

Old cogs, new tricks: the evolution of gene expression in a chromatin context

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.

Peak calling by Sparse Enrichment Analysis for CUT&RUN chromatin profiling

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.

Improved CUT&RUN chromatin profiling tools

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.

Automated in situ chromatin profiling efficiently resolves cell types and gene regulatory programs

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.

Transcription start site profiling uncovers divergent transcription and enhancer-associated RNAs in Drosophila melanogaster

Background

High-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.

Results

The 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 developmental transcription programs.

Conclusions

High-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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4510-7) contains supplementary material, which is available to authorized users.

An Animal Model for Genetic Analysis of Multi-Gene Families: Cloning and Transgenesis of Large Tandemly Repeated Histone Gene Clusters

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 gene replacement reveals a post-transcriptional role for H3K36 in maintaining metazoan transcriptome fidelity

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

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

Interrogating the function of metazoan histones using engineered gene clusters

Histones 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.

Poly(A) RNA-binding proteins and polyadenosine RNA: new members and novel functions

Poly(A) RNA-binding proteins (Pabs) bind with high affinity and specificity to polyadenosine RNA. Textbook models show a nuclear Pab, PABPN1, and a cytoplasmic Pab, PABPC, where the nuclear PABPN1 modulates poly(A) tail length and the cytoplasmic PABPC stabilizes poly(A) RNA in the cytoplasm and also enhances translation. While these conventional roles are critically important, the Pab family has expanded recently both in number and in function. A number of novel roles have emerged for both PAPBPN1 and PABPC that contribute to the fine-tuning of gene expression. Furthermore, as the characterization of the nucleic acid binding properties of RNA-binding proteins advances, additional proteins that show high affinity and specificity for polyadenosine RNA are being discovered. With this expansion of the Pab family comes a concomitant increase in the potential for Pabs to modulate gene expression. Further complication comes from an expansion of the potential binding sites for Pab proteins as revealed by an analysis of templated polyadenosine stretches present within the transcriptome. Thus, Pabs could influence mRNA fate and function not only by binding to the nontemplated poly(A) tail but also to internal stretches of adenosine. Understanding the diverse functions of Pab proteins is not only critical to understand how gene expression is regulated but also to understand the molecular basis for tissue-specific diseases that occur when Pab proteins are altered. Here we describe both conventional and recently emerged functions for PABPN1 and PABPC and then introduce and discuss three new Pab family members, ZC3H14, hnRNP-Q1, and LARP4.

Developmental arrest of Drosophila survival motor neuron (Smn) mutants accounts for differences in expression of minor intron-containing genes

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.