Search

Abstract Chromosomal translocations involving the mixed-lineage leukemia (MLL) locus generate potent oncogenic fusion proteins (oncoproteins) that disrupt regulation of developmental gene expression. By profiling the oncoprotein-target sites of 36 broadly representative MLL-rearranged leukemia samples, including three samples that underwent a lymphoid-to-myeloid lineage-switching event in response to therapy, we find the genomic enrichment of the oncoprotein is highly variable between samples and subject to dynamic regulation. At high levels of expression, the oncoproteins preferentially activate either an acute lymphoblastic leukemia (ALL) program, enriched for pro-B-cell genes, or an acute myeloid leukemia (AML) program, enriched for hematopoietic-stem-cell genes. The fusion-partner-specific-binding patterns over these gene sets are highly correlated with the prevalence of each mutation in ALL versus AML. In lineage-switching samples the oncoprotein levels are reduced and the oncoproteins preferentially activate granulocyte-monocyte progenitor (GMP) genes. In a sample that lineage switched during treatment with the menin inhibitor revumenib, the oncoprotein and menin are reduced to undetectable levels, but ENL, a transcriptional cofactor of the oncoprotein, persists on numerous oncoprotein-target loci, including genes in the GMP-like lineage-switching program. We propose MLL oncoproteins promote lineage-switching events through dynamic chromatin binding at lineage-specific target genes, and may support resistance to menin inhibitors through similar changes in chromatin occupancy.

Summary: Centromere length changes occurring during somatic cell divisions can be estimated by quantifying the copy numbers (CNs) of higher-order repeats (HORs), which are nested repeats of monomers that comprise centromeric arrays. Here, we present a protocol for single-cell isolation for clonal evolution followed by droplet digital PCR-based quantification. The assay measures HOR CNs across subclones to determine the frequency and degree of changes in HOR CNs. This protocol tests the underlying molecular mechanisms responsible for rapid centromere sequence evolution.For complete details on the use and execution of this protocol, please refer to Showman et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Summary: Human centromeres are located within α-satellite arrays and evolve rapidly, which can lead to individual variation in array length. Proposed mechanisms for such alterations in length are unequal crossover between sister chromatids, gene conversion, and break-induced replication. However, the underlying molecular mechanisms responsible for the massive, complex, and homogeneous organization of centromeric arrays have not been experimentally validated. Here, we use droplet digital PCR assays to demonstrate that centromeric arrays can expand and contract within ∼20 somatic cell divisions of an alternative lengthening of telomere (ALT)-positive cell line. We find that the frequency of array variation among single-cell-derived subclones ranges from a minimum of ∼7% to a maximum of ∼100%. Further clonal evolution revealed that centromere expansion is favored over contraction. We find that the homologous recombination protein RAD52 and the helicase PIF1 are required for extensive array change, suggesting that centromere sequence evolution can occur via break-induced replication.

Abstract For more than a century, formalin-fixed paraffin-embedded (FFPE) sample preparation has been the preferred method for long-term preservation of biological material. However, the use of FFPE samples for epigenomic studies has been difficult because of chromatin damage from long exposure to high concentrations of formaldehyde. Previously, we introduced Cleavage Under Targeted Accessible Chromatin (CUTAC), an antibody-targeted chromatin accessibility mapping protocol based on CUT&Tag. Here we show that simple modifications of our CUTAC protocol either in single tubes or directly on slides produce high-resolution maps of paused RNA Polymerase II at enhancers and promoters using FFPE samples. We find that transcriptional regulatory element differences produced by FFPE-CUTAC distinguish between mouse brain tumors and identify and map regulatory element markers with high confidence and precision, including microRNAs not detectable by RNA-seq. Our simple workflows make possible affordable epigenomic profiling of archived biological samples for biomarker identification, clinical applications and retrospective studies.

Spermatogenesis in the Drosophila male germline proceeds through a unique transcriptional program controlled both by germline-specific transcription factors and by testis-specific versions of core transcriptional machinery. This program includes the activation of genes on the heterochromatic Y chromosome, and reduced transcription from the X chromosome, but how expression from these sex chromosomes is regulated has not been defined. To resolve this, we profiled active chromatin features in the testes from wildtype and meiotic arrest mutants and integrate this with single-cell gene expression data from the Fly Cell Atlas. These data assign the timing of promoter activation for genes with germline-enriched expression throughout spermatogenesis, and general alterations of promoter regulation in germline cells. By profiling both active RNA polymerase II and histone modifications in isolated spermatocytes, we detail widespread patterns associated with regulation of the sex chromosomes. Our results demonstrate that the X chromosome is not enriched for silencing histone modifications, implying that sex chromosome inactivation does not occur in the Drosophila male germline. Instead, a lack of dosage compensation in spermatocytes accounts for the reduced expression from this chromosome. Finally, profiling uncovers dramatic ubiquitinylation of histone H2A and lysine-16 acetylation of histone H4 across the Y chromosome in spermatocytes that may contribute to the activation of this heterochromatic chromosome.

Abstract The common histones H2A, H2B, H3, and H4 are the characteristic components of eukaryotic nucleosomes, which function to wrap DNA and compact the genome as well as to regulate access to DNA for transcription and replication in all eukaryotes. In the past two decades, histones have also been found to be encoded in some DNA viruses, where their functions and properties are largely unknown, though recently histones from two related viruses have been shown to form nucleosome-like structures in vitro. Viral histones can be highly similar to eukaryotic histones in primary sequence, suggesting they have been recently picked up from eukaryotic hosts, or they can be radically divergent in primary sequence and may occur as conjoined histone doublets, triplets, or quadruplets, suggesting ancient origins prior to the divergence of modern eukaryotes. Here, we review what is known of viral histones and discuss their possible origins and functions. We consider how the viral life cycle may affect their properties and histories, and reflect on the possible roles of viruses in the origin of the nucleus of modern eukaryotic cells.

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

Short H2A variants are testis-specific histones that destabilize nucleosomes during spermatogenesis. In this study, the authors show that these variants are expressed in an array of different cancers and identify splicing changes associated with nucleosome instability in these malignancies.

Firre encodes a lncRNA involved in nuclear organization in mammals. Here, the authors find that allelic deletion of Firre on the active X chromosome (Xa) results in dose-dependent loss of histone H3K27me3 on the inactive X chromosome (Xi), along with other trans-acting effects, including disruption of the perinuclear location and minor dysregulation of gene expression.

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.

The Fbw7 ubiquitin ligase targets many proteins for proteasomal degradation, which include oncogenic transcription factors (TFs) (e.g., c-Myc, c-Jun, and Notch). Fbw7 is a tumor suppressor and tumors often contain mutations in FBXW7, the gene that encodes Fbw7. The complexity of its substrate network has obscured the mechanisms of Fbw7-associated tumorigenesis, yet this understanding is needed for developing therapies. We used an integrated approach employing RNA-Seq and high-resolution mapping (cleavage under target and release using nuclease) of histone modifications and TF occupancy (c-Jun and c-Myc) to examine the combinatorial effects of misregulated Fbw7 substrates in colorectal cancer (CRC) cells with engineered tumor-associated FBXW7 null or missense mutations. Both Fbw7 mutations caused widespread transcriptional changes associated with active chromatin and altered TF occupancy: some were common to both Fbw7 mutant cell lines, whereas others were mutation specific. We identified loci where both Jun and Myc were coregulated by Fbw7, suggesting that substrates may have synergistic effects. One coregulated gene was CIITA, the master regulator of MHC Class II gene expression. Fbw7 loss increased MHC Class II expression and Fbw7 mutations were correlated with increased CIITA expression in TCGA colorectal tumors and cell lines, which may have immunotherapeutic implications for Fbw7-associated cancers. Analogous studies in neural stem cells in which FBXW7 had been acutely deleted closely mirrored the results in CRC cells. Gene set enrichment analyses revealed Fbw7-associated pathways that were conserved across both cell types that may reflect fundamental Fbw7 functions. These analyses provide a framework for understanding normal and neoplastic context-specific Fbw7 functions.

Development proceeds by the activation of genes by transcription factors and the inactivation of others by chromatin-mediated gene silencing. In certain cases development can be reversed or redirected by mis-expression of master regulator transcription factors. This must involve the activation of previously silenced genes, and such developmental aberrations are thought to underlie a variety of cancers. Here, we express the wing-specific Vestigial master regulator to reprogram the developing eye, and test the role of silencing in reprogramming using an H3.3K27M oncohistone mutation that dominantly inhibits histone H3K27 trimethylation. We find that production of the oncohistone blocks eye-to-wing reprogramming. CUT&Tag chromatin profiling of mutant tissues shows that H3K27me3 of domains is generally reduced upon oncohistone production, suggesting that a previous developmental program must be silenced for effective transformation. Strikingly, Vg and H3.3K27M synergize to stimulate overgrowth of eye tissue, a phenotype that resembles that of mutations in Polycomb silencing components. Transcriptome profiling of elongating RNA Polymerase II implicates the mis-regulation of signaling factors in overgrowth. Our results demonstrate that growth dysregulation can result from the simple combination of crippled silencing and transcription factor mis-expression, an effect that may explain the origins of oncohistone-bearing cancers.

We previously introduced Cleavage Under Targets & Tagmentation (CUT&Tag), an epigenomic profiling method in which antibody tethering of the Tn5 transposase to a chromatin epitope of interest maps specific chromatin features in small samples and single cells. With CUT&Tag, intact cells or nuclei are permeabilized, followed by successive addition of a primary antibody, a secondary antibody, and a chimeric Protein A-Transposase fusion protein that binds to the antibody. Addition of Mg++ activates the transposase and inserts sequencing adapters into adjacent DNA in situ. We have since adapted CUT&Tag to also map chromatin accessibility by simply modifying the transposase activation conditions when using histone H3K4me2, H3K4me3, or Serine-5-phosphorylated RNA Polymerase II antibodies. Using these antibodies, we redirect the tagmentation of accessible DNA sites to produce chromatin accessibility maps with exceptionally high signal-to-noise and resolution. All steps from nuclei to amplified sequencing-ready libraries are performed in single PCR tubes using non-toxic reagents and inexpensive equipment, making our simplified strategy for simultaneous chromatin profiling and accessibility mapping suitable for the lab, home workbench, or classroom.

Abstract Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, its aggressive endo-/exo-nuclease activities make MNase-seq unreliable for determining nucleosome occupancies, because cleavages within linker regions produce oligo- and mono-nucleosomes, whereas cleavages within nucleosomes destroy them. Here, we introduce a theoretical framework for predicting nucleosome occupancies and an experimental protocol with appropriate spike-in normalization that confirms our theory and provides accurate occupancy levels over an MNase digestion time course. As with human cells, we observe no overall differences in nucleosome occupancies between Drosophila euchromatin and heterochromatin, which implies that heterochromatic compaction does not reduce MNase accessibility of linker DNA.

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

Understanding gene regulation will require mapping specific chromain features in a small number of cells at high resolution. Here the authors describe CUT&Tag, which uses antibody-mediated tethering of Tn5 transposase to a chromatin protein to generate high resolution libraries.

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

Histone variants expand chromatin functions in eukaryote genomes. H2A.B genes are testis-expressed short histone H2A variants that arose in placental mammals. Their biological functions remain largely unknown. To investigate their function, we generated a knockout (KO) model that disrupts all 3 H2A.B genes in mice. We show that H2A.B KO males have globally altered chromatin structure in postmeiotic germ cells. Yet, they do not show impaired spermatogenesis or testis function. Instead, we find that H2A.B plays a crucial role postfertilization. Crosses between H2A.B KO males and females yield embryos with lower viability and reduced size. Using a series of genetic crosses that separate parental and zygotic contributions, we show that the H2A.B status of both the father and mother, but not of the zygote, affects embryonic viability and growth during gestation. We conclude that H2A.B is a novel parental-effect gene, establishing a role for short H2A histone variants in mammalian development. We posit that parental antagonism over embryonic growth drove the origin and ongoing diversification of short histone H2A variants in placental mammals.

Chromatin accessibility mapping is a powerful approach to identify potential regulatory elements. A popular example is ATAC-seq, whereby Tn5 transposase inserts sequencing adapters into accessible DNA (‘tagmentation’). CUT&Tag is a tagmentation-based epigenomic profiling method in which antibody tethering of Tn5 to a chromatin epitope of interest profiles specific chromatin features in small samples and single cells. Here, we show that by simply modifying the tagmentation conditions for histone H3K4me2 or H3K4me3 CUT&Tag, antibody-tethered tagmentation of accessible DNA sites is redirected to produce chromatin accessibility maps that are indistinguishable from the best ATAC-seq maps. Thus, chromatin accessibility maps can be produced in parallel with CUT&Tag maps of other epitopes with all steps from nuclei to amplified sequencing-ready libraries performed in single PCR tubes in the laboratory or on a home workbench. As H3K4 methylation is produced by transcription at promoters and enhancers, our method identifies transcription-coupled accessible regulatory sites.

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.

Abstract We developed a chemical cleavage method that releases single nucleosome dyad-containing fragments, allowing us to precisely map both single nucleosomes and linkers with high accuracy genome-wide in yeast. Our single nucleosome positioning data reveal that nucleosomes occupy preferred positions that differ by integral multiples of the DNA helical repeat. By comparing nucleosome dyad positioning maps to existing genomic and transcriptomic data, we evaluated the contributions of sequence, transcription, and histones H1 and H2A.Z in defining the chromatin landscape. We present a biophysical model that neglects DNA sequence and shows that steric occlusion suffices to explain the salient features of nucleosome positioning.

Summary: The DUX4 transcription factor is briefly expressed in the early cleavage-stage embryo, where it induces an early wave of zygotic gene transcription, whereas its mis-expression in skeletal muscle causes the muscular dystrophy facioscapulohumeral dystrophy (FSHD). Here, we show that DUX4 induces the expression of the histone variants H3.X and H3.Y. We have used a myoblast cell line with doxycycline-inducible DUX4 to show that these histone variants are incorporated throughout the body of DUX4-induced genes. Following a brief pulse of DUX4, these histones contribute to greater perdurance and to enhanced re-activation of DUX4 target gene expression. These findings provide a model for H3.X/Y as a chromatin mechanism that facilitates the expression of DUX4 target genes subsequent to a brief pulse of DUX4 expression. : Resnick et al. show that the DUX4-induced histone variants, H3.X and H3.Y, incorporate into the chromatin of DUX4-induced genes, making them more sensitive to subsequent expression. This suggests a mechanism for how the brief expression of DUX4 can establish a memory of its transcriptional network. Keywords: DUX4, H3.X, H3.Y, facioscapulohumeral dystrophy, transcription, histone variants

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.

Satellite DNA consists of abundant tandem repeats that play important roles in cellular processes, including chromosome segregation, genome organization and chromosome end protection. Most satellite DNA repeat units are either of nucleosomal length or 5–10 bp long and occupy centromeric, pericentromeric or telomeric regions. Due to high repetitiveness, satellite DNA sequences have largely been absent from genome assemblies. Although few conserved satellite-specific sequence motifs have been identified, DNA curvature, dyad symmetries and inverted repeats are features of various satellite DNAs in several organisms. Satellite DNA sequences are either embedded in highly compact gene-poor heterochromatin or specialized chromatin that is distinct from euchromatin. Nevertheless, some satellite DNAs are transcribed into non-coding RNAs that may play important roles in satellite DNA function. Intriguingly, satellite DNAs are among the most rapidly evolving genomic elements, such that a large fraction is species-specific in most organisms. Here we describe the different classes of satellite DNA sequences, their satellite-specific chromatin features, and how these features may contribute to satellite DNA biology and evolution. We also discuss how the evolution of functional satellite DNA classes may contribute to speciation in plants and animals.

We describe Cleavage Under Targets and Release Using Nuclease (CUT&RUN), a chromatin profiling strategy in which antibody-targeted controlled cleavage by micrococcal nuclease releases specific protein-DNA complexes into the supernatant for paired-end DNA sequencing. Unlike Chromatin Immunoprecipitation (ChIP), which fragments and solubilizes total chromatin, CUT&RUN is performed in situ, allowing for both quantitative high-resolution chromatin mapping and probing of the local chromatin environment. When applied to yeast and human nuclei, CUT&RUN yielded precise transcription factor profiles while avoiding crosslinking and solubilization issues. CUT&RUN is simple to perform and is inherently robust, with extremely low backgrounds requiring only ~1/10th the sequencing depth as ChIP, making CUT&RUN especially cost-effective for transcription factor and chromatin profiling. When used in conjunction with native ChIP-seq and applied to human CTCF, CUT&RUN mapped directional long range contact sites at high resolution. We conclude that in situ mapping of protein-DNA interactions by CUT&RUN is an attractive alternative to ChIP-seq.

Chromatin immunoprecipitation (ChIP) and its derivatives are the main techniques used to determine transcription factor binding sites. However, conventional ChIP with sequencing (ChIP-seq) has problems with poor resolution, and newer techniques require significant experimental alterations and complex bioinformatics. Previously, we have used a new crosslinking ChIP-seq protocol (X-ChIP-seq) to perform high-resolution mapping of RNA Polymerase II (Skene et al., 2014). Here, we build upon this work and compare X-ChIP-seq to existing methodologies. By using micrococcal nuclease, which has both endo- and exo-nuclease activity, to fragment the chromatin and thereby generate precise protein–DNA footprints, high-resolution X-ChIP-seq achieves single base-pair resolution of transcription factor binding. A significant advantage of this protocol is the minimal alteration to the conventional ChIP-seq workflow and simple bioinformatic processing.

Genomes must balance active suppression of transposable elements (TEs) with the need to maintain gene expression. In Arabidopsis, euchromatic TEs are targeted by RNA-directed DNA methylation (RdDM). Conversely, active DNA demethylation prevents accumulation of methylation at genes proximal to these TEs. It is unknown how a cellular balance between methylation and demethylation activities is achieved. Here we show that both RdDM and DNA demethylation are highly active at a TE proximal to the major DNA demethylase gene ROS1. Unexpectedly, and in contrast to most other genomic targets, expression of ROS1 is promoted by DNA methylation and antagonized by DNA demethylation. We demonstrate that inducing methylation in the ROS1 proximal region is sufficient to restore ROS1 expression in an RdDM mutant. Additionally, methylation-sensitive expression of ROS1 is conserved in other species, suggesting it is adaptive. We propose that the ROS1 locus functions as an epigenetic rheostat, tuning the level of demethylase activity in response to methylation alterations, thus ensuring epigenomic stability.

Heterochromatin comprises a significant component of many eukaryotic genomes. In comparison to euchromatin, heterochromatin is gene poor, transposon rich, and late replicating. It serves many important biological roles, from gene silencing to accurate chromosome segregation, yet little is known about the evolutionary constraints that shape heterochromatin. A complementary approach to the traditional one of directly studying heterochromatic DNA sequence is to study the evolution of proteins that bind and define heterochromatin. One of the best markers for heterochromatin is the heterochromatin protein 1 (HP1), which is an essential, nonhistone chromosomal protein. Here we investigate the molecular evolution of five HP1 paralogs present in Drosophila melanogaster. Three of these paralogs have ubiquitous expression patterns in adult Drosophila tissues, whereas HP1D/rhino and HP1E are expressed predominantly in ovaries and testes respectively. The HP1 paralogs also have distinct localization preferences in Drosophila cells. Thus, Rhino localizes to the heterochromatic compartment in Drosophila tissue culture cells, but in a pattern distinct from HP1A and lysine-9 dimethylated H3. Using molecular evolution and population genetic analyses, we find that rhino has been subject to positive selection in all three domains of the protein: the N-terminal chromo domain, the C-terminal chromo-shadow domain, and the hinge region that connects these two modules. Maximum likelihood analysis of rhino sequences from 20 species of Drosophila reveals that a small number of residues of the chromo and shadow domains have been subject to repeated positive selection. The rapid and positive selection of rhino is highly unusual for a gene encoding a chromosomal protein and suggests that rhino is involved in a genetic conflict that affects the germline, belying the notion that heterochromatin is simply a passive recipient of "junk DNA" in eukaryotic genomes.

Genes occupy only a minor fraction of genomes such as ours; however, histone and nonhistone chromosomal proteins and methylated DNA bases are distributed over both genic and nongenic regions. These widespread “epigenomic” features can be mapped and characterized by alternative applications of the same microarray technologies that have been used for conventional transcriptional profiling. Here we describe diverse microarray-based strategies for profiling patterns of DNA methylation, DNA replication, DNA binding, and chromatin-associated proteins and histone modifications. The rapid progress that is being made in developing and applying epigenomic profiling methods and the increasing availability of microarrays mean that epigenomic profiling is likely to become a standard research tool for understanding chromatin structure and gene expression during development.

Nature Communications6: Article number:8733 (2015); Published: 22 October 2015; Updated 5 June 2017 We have been alerted to an inconsistency regarding the control used as input in the DNA shape analysis shown in Fig. 7 of this paper. As described in the Methods section of the Article, ‘FASTA sequences generated from the previously generated 100 bp windows centered on motif match midpoints were used as input for our DNA shape method for high-throughput prediction of DNA structural features.

Faithful chromosome segregation in all eukaryotes relies on centromeres, the chromosomal sites that recruit kinetochore proteins and mediate spindle attachment during cell division. The centromeric histone H3 variant, CenH3, is the defining chromatin component of centromeres in most eukaryotes, including animals, fungi, plants, and protists. In this study, using detailed genomic and transcriptome analyses, we show that CenH3 was lost independently in at least four lineages of insects. Each of these lineages represents an independent transition from monocentricity (centromeric determinants localized to a single chromosomal region) to holocentricity (centromeric determinants extended over the entire chromosomal length) as ancient as 300 million years ago. Holocentric insects therefore contain a CenH3-independent centromere, different from almost all the other eukaryotes. We propose that ancient transitions to holocentricity in insects obviated the need to maintain CenH3, which is otherwise essential in most eukaryotes, including other holocentrics.

In budding yeast, a single cenH3 (Cse4) nucleosome occupies the ∼120-bp functional centromere, however conflicting structural models for the particle have been proposed. To resolve this controversy, we have applied H4S47C-anchored cleavage mapping, which reveals the precise position of histone H4 in every nucleosome in the genome. We find that cleavage patterns at centromeres are unique within the genome and are incompatible with symmetrical structures, including octameric nucleosomes and (Cse4/H4)2 tetrasomes. Centromere cleavage patterns are compatible with a precisely positioned core structure, one in which each of the 16 yeast centromeres is occupied by oppositely oriented Cse4/H4/H2A/H2B hemisomes in two rotational phases within the population. Centromere-specific hemisomes are also inferred from distances observed between closely-spaced H4 cleavages, as predicted from structural modeling. Our results indicate that the orientation and rotational position of the stable hemisome at each yeast centromere is not specified by the functional centromere sequence.

Centromeres vary greatly in size and sequence composition, ranging from ‘point’ centromeres with a single cenH3-containing nucleosome to ‘regional’ centromeres embedded in tandemly repeated sequences to holocentromeres that extend along the length of entire chromosomes. Point centromeres are defined by sequence, whereas regional and holocentromeres are epigenetically defined by the location of cenH3-containing nucleosomes. In this study, we show that Caenorhabditis elegans holocentromeres are organized as dispersed but discretely localized point centromeres, each forming a single cenH3-containing nucleosome. These centromeric sites co-localize with kinetochore components, and their occupancy is dependent on the cenH3 loading machinery. These sites coincide with non-specific binding sites for multiple transcription factors (‘HOT’ sites), which become occupied when cenH3 is lost. Our results show that the point centromere is the basic unit of holocentric organization in support of the classical polycentric model for holocentromeres, and provide a mechanistic basis for understanding how centromeric chromatin might be maintained.

RNA polymerase II (PolII) transcribes RNA within a chromatin context, with nucleosomes acting as barriers to transcription. Despite these barriers, transcription through chromatin in vivo is highly efficient, suggesting the existence of factors that overcome this obstacle. To increase the resolution obtained by standard chromatin immunoprecipitation, we developed a novel strategy using micrococcal nuclease digestion of cross-linked chromatin. We find that the chromatin remodeler Chd1 is recruited to promoter proximal nucleosomes of genes undergoing active transcription, where Chd1 is responsible for the vast majority of PolII-directed nucleosome turnover. The expression of a dominant negative form of Chd1 results in increased stalling of PolII past the entry site of the promoter proximal nucleosomes. We find that Chd1 evicts nucleosomes downstream of the promoter in order to overcome the nucleosomal barrier and enable PolII promoter escape, thus providing mechanistic insight into the role of Chd1 in transcription and pluripotency.

ATP-dependent nucleosome remodelers influence genetic processes by altering nucleosome occupancy, positioning, and composition. In vitro, Saccharomyces cerevisiae ISWI and CHD remodelers require ∼30-85 bp of extranucleosomal DNA to reposition nucleosomes, but linker DNA in S. cerevisiae averages

Differential expression of maternally and paternally inherited alleles of a gene is referred to as gene imprinting, a form of epigenetic gene regulation common to flowering plants and mammals. In plants, imprinting primarily occurs in the endosperm, a seed tissue that supports the embryo during its growth and development. Previously, we demonstrated that widespread DNA demethylation at remnants of transposable elements accompanies endosperm development and that a subset of these methylation changes are associated with gene imprinting. Here we assay imprinted gene expression genome-wide by performing high-throughput sequencing of RNA derived from seeds of reciprocal intraspecific crosses. We identify more than 200 loci that exhibit parent-of-origin effects on gene expression in the endosperm, including a large number of transcription factors, hormone biosynthesis and response genes, and genes that encode regulators of epigenetic information, such as methylcytosine binding proteins, histone methyltransferases, and chromatin remodelers. The majority of these genes are partially, rather than completely, imprinted, suggesting that gene dosage regulation is an important aspect of imprinted gene expression.

Insulators are DNA sequences that control the interactions among genomic regulatory elements and act as chromatin boundaries. A thorough understanding of their location and function is necessary to address the complexities of metazoan gene regulation. We studied by ChIP-chip the genome-wide binding sites of 6 insulator-associated proteins-dCTCF, CP190, BEAF-32, Su(Hw), Mod(mdg4), and GAF-to obtain the first comprehensive map of insulator elements in Drosophila embryos. We identify over 14,000 putative insulators, including all classically defined insulators. We find two major classes of insulators defined by dCTCF/CP190/BEAF-32 and Su(Hw), respectively. Distributional analyses of insulators revealed that particular sub-classes of insulator elements are excluded between cis-regulatory elements and their target promoters; divide differentially expressed, alternative, and divergent promoters; act as chromatin boundaries; are associated with chromosomal breakpoints among species; and are embedded within active chromatin domains. Together, these results provide a map demarcating the boundaries of gene regulatory units and a framework for understanding insulator function during the development and evolution of Drosophila.

Centromeres are sites for assembly of the chromosomal structures that mediate faithful segregation at mitosis and meiosis. Plant and animal centromeres are typically located in megabase-sized arrays of tandem satellite repeats, making their precise mapping difficult. However, some rice centromeres are largely embedded in nonsatellite DNA, providing an excellent model to study centromere structure and evolution. We used chromatin immunoprecipitation and 454 sequencing to define the boundaries of nine of the 12 centromeres of rice. Centromere regions from chromosomes 8 and 9 were found to share synteny, most likely reflecting an ancient genome duplication. For four centromeres, we mapped discrete subdomains of binding by the centromeric histone variant CENH3. These subdomains were depleted in both intact and nonfunctional genes relative to interspersed subdomains lacking CENH3. The intergenic location of rice centromeric chromatin resembles the situation for human neocentromeres and supports a model of the evolution of centromeres from gene-poor regions.

Centromeres, the specialized chromatin structures that are responsible for equal segregation of chromosomes at mitosis, are epigenetically maintained by a centromere-specific histone H3 variant (CenH3). However, the mechanistic basis for centromere maintenance is unknown. We investigated biochemical properties of CenH3 nucleosomes from Drosophila melanogaster cells. Cross-linking of CenH3 nucleosomes identifies heterotypic tetramers containing one copy of CenH3, H2A, H2B, and H4 each. Interphase CenH3 particles display a stable association of approximately 120 DNA base pairs. Purified centromeric nucleosomal arrays have typical "beads-on-a-string" appearance by electron microscopy but appear to resist condensation under physiological conditions. Atomic force microscopy reveals that native CenH3-containing nucleosomes are only half as high as canonical octameric nucleosomes are, confirming that the tetrameric structure detected by cross-linking comprises the entire interphase nucleosome particle. This demonstration of stable half-nucleosomes in vivo provides a possible basis for the instability of centromeric nucleosomes that are deposited in euchromatic regions, which might help maintain centromere identity.

Germline chromatin undergoes dramatic remodeling events involving histone variants during the life cycle of an organism. A universal histone variant, H3.3, is incorporated at sites of active transcription throughout the cell cycle. The presence of H3.3 in chromatin indicates histone turnover, which is the energy-dependent removal of preexisting histones and replacement with new histones. H3.3 is also incorporated during decondensation of the Drosophila sperm pronucleus, indicating a direct role in chromatin remodeling upon fertilization. Here we present a system to monitor histone turnover and chromatin remodeling during Caenorhabditis elegans development by following the developmental dynamics of H3.3. We generated worm strains expressing green fluorescent protein- or yellow fluorescent protein-fused histone H3.3 proteins, HIS-71 and HIS-72. We found that H3.3 is retained in mature sperm chromatin, raising the possibility that it transmits epigenetic information via the male germline. Upon fertilization, maternal H3.3 enters both male and female pronuclei and is incorporated into paternal chromatin, apparently before the onset of embryonic transcription, suggesting that H3.3 can be incorporated independent of transcription. In early embryos, H3.3 becomes specifically depleted from primordial germ cells. Strikingly, the X chromosome becomes deficient in H3.3 during gametogenesis, indicating a low level of histone turnover. These results raise the possibility that the asymmetry in histone turnover between the X chromosome and autosomes is established during gametogenesis. H3.3 patterns are similar to patterns of H3K4 methylation in the primordial germ cells and on the X chromosome during gametogenesis, suggesting that histone turnover and modification are coupled processes. Our demonstration of dynamic H3.3 incorporation in nondividing cells provides a mechanistic basis for chromatin changes during germ cell development.

Eukaryotic genomes can usurp enzymatic functions encoded by mobile elements for their own use. A particularly interesting kind of acquisition involves the domestication of retroviral envelope genes, which confer infectious membrane-fusion ability to retroviruses. So far, these examples have been limited to vertebrate genomes, including primates where the domesticated envelope is under purifying selection to assist placental function. Here, we show that in Drosophila genomes, a previously unannotated gene (CG4715, renamed Iris) was domesticated from a novel, active Kanga lineage of insect retroviruses at least 25 million years ago, and has since been maintained as a host gene that is expressed in all adult tissues. Iris and the envelope genes from Kanga retroviruses are homologous to those found in insect baculoviruses and gypsy and roo insect retroviruses. Two separate envelope domestications from the Kanga and roo retroviruses have taken place, in fruit fly and mosquito genomes, respectively. Whereas retroviral envelopes are proteolytically cleaved into the ligand-interaction and membrane-fusion domains, Iris appears to lack this cleavage site. In the takahashii/suzukii species groups of Drosophila, we find that Iris has tandemly duplicated to give rise to two genes (Iris-A and Iris-B). Iris-B has significantly diverged from the Iris-A lineage, primarily because of the "invention" of an intron de novo in what was previously exonic sequence. Unlike domesticated retroviral envelope genes in mammals, we find that Iris has been subject to strong positive selection between Drosophila species. The rapid, adaptive evolution of Iris is sufficient to unambiguously distinguish the phylogenies of three closely related sibling species of Drosophila (D. simulans, D. sechellia, and D. mauritiana), a discriminative power previously described only for a putative "speciation gene." Iris represents the first instance of a retroviral envelope-derived host gene outside vertebrates. It is also the first example of a retroviral envelope gene that has been found to be subject to positive selection following its domestication. The unusual selective pressures acting on Iris suggest that it is an active participant in an ongoing genetic conflict. We propose a model in which Iris has "switched sides," having been recruited by host genomes to combat baculoviruses and retroviruses, which employ homologous envelope genes to mediate infection.