Jonathan H. Dennis

Topologically associating domains are stable units of replication-timing regulation

Eukaryotic chromosomes replicate in a temporal order known as the replication-timing program. In mammals, replication timing is cell-type-specific with at least half the genome switching replication timing during development, primarily in units of 400–800 kilobases (‘replication domains’), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements. Early and late replication correlate, respectively, with open and closed three-dimensional chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, late replication correlates with lamina-associated domains (LADs). Recent Hi-C mapping has unveiled substructure within chromatin compartments called topologically associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to replication domains. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure. Here we localize boundaries of replication domains to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, replication domain boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure replication domain boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type-specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell-type-specific sub-nuclear compartmentalization and replication timing with developmentally stable structural domains and offer a unified model for large-scale chromosome structure and function.

Predicting Human Nucleosome Occupancy from Primary Sequence

Nucleosomes are the fundamental repeating unit of chromatin and comprise the structural building blocks of the living eukaryotic genome. Micrococcal nuclease (MNase) has long been used to delineate nucleosomal organization. Microarray-based nucleosome mapping experiments in yeast chromatin have revealed regularly-spaced translational phasing of nucleosomes. These data have been used to train computational models of sequence-directed nuclesosome positioning, which have identified ubiquitous strong intrinsic nucleosome positioning signals. Here, we successfully apply this approach to nucleosome positioning experiments from human chromatin. The predictions made by the human-trained and yeast-trained models are strongly correlated, suggesting a shared mechanism for sequence-based determination of nucleosome occupancy. In addition, we observed striking complementarity between classifiers trained on experimental data from weakly versus heavily digested MNase samples. In the former case, the resulting model accurately identifies nucleosome-forming sequences; in the latter, the classifier excels at identifying nucleosome-free regions. Using this model we are able to identify several characteristics of nucleosome-forming and nucleosome-disfavoring sequences. First, by combining results from each classifier applied de novo across the human ENCODE regions, the classifier reveals distinct sequence composition and periodicity features of nucleosome-forming and nucleosome-disfavoring sequences. Short runs of dinucleotide repeat appear as a hallmark of nucleosome-disfavoring sequences, while nucleosome-forming sequences contain short periodic runs of GC base pairs. Second, we show that nucleosome phasing is most frequently predicted flanking nucleosome-free regions. The results suggest that the major mechanism of nucleosome positioning in vivo is boundary-event-driven and affirm the classical statistical positioning theory of nucleosome organization.

Pre-replication complex proteins assemble at regions of low nucleosome occupancy within the Chinese hamster dihydrofolate reductase initiation zone

Genome-scale mapping of pre-replication complex proteins has not been reported in mammalian cells. Poor enrichment of these proteins at specific sites may be due to dispersed binding, poor epitope availability or cell cycle stage-specific binding. Here, we have mapped sites of biotin-tagged ORC and MCM protein binding in G1-synchronized populations of Chinese hamster cells harboring amplified copies of the dihydrofolate reductase (DHFR) locus, using avidin-affinity purification of biotinylated chromatin followed by high-density microarray analysis across the DHFR locus. We have identified several sites of significant enrichment for both complexes distributed throughout the previously identified initiation zone. Analysis of the frequency of initiations across stretched DNA fibers from the DHFR locus confirmed a broad zone of de-localized initiation activity surrounding the sites of ORC and MCM enrichment. Mapping positions of mononucleosomal DNA empirically and computing nucleosome-positioning information in silico revealed that ORC and MCM map to regions of low measured and predicted nucleosome occupancy. Our results demonstrate that specific sites of ORC and MCM enrichment can be detected within a mammalian intitiation zone, and suggest that initiation zones may be regions of generally low nucleosome occupancy where flexible nucleosome positioning permits flexible pre-RC assembly sites.

Chromatin-interaction compartment switch at developmentally regulated chromosomal domains reveals an unusual principle of chromatin folding

Several 400- to 800-kb murine chromosome domains switch from early to late replication during loss of pluripotency, accompanied by a stable form of gene silencing that is resistant to reprogramming. We found that, whereas enhanced nuclease accessibility correlated with early replication genome-wide, domains that switch replication timing during differentiation were exceptionally inaccessible even when early-replicating. Nonetheless, two domains studied in detail exhibited substantial changes in transcriptional activity and higher-order chromatin unfolding confined to the region of replication timing change. Chromosome conformation capture (4C) data revealed that in the unfolded state in embryonic stem cells, these domains interacted preferentially with the early-replicating chromatin compartment, rarely interacting even with flanking late-replicating domains, whereas after differentiation, these same domains preferentially associated with late-replicating chromatin, including flanking domains. In both configurations they retained local boundaries of self-interaction, supporting the replication domain model of replication-timing regulation. Our results reveal a principle of developmentally regulated, large-scale chromosome folding involving a subnuclear compartment switch of inaccessible chromatin. This unusual level of regulation may underlie resistance to reprogramming in replication-timing switch regions.

Expression of the Brn-3b transcription factor correlates with expression of HSP-27 in breast cancer biopsies and is required for maximal activation of the HSP-27 promoter

In breast cancer, overexpression of the small heat shock protein, HSP-27, is associated with increased anchorage-independent growth, increased invasiveness, and resistance to chemotherapeutic drugs and is associated with poor prognosis and reduced disease-free survival. Therefore, factors that increase the expression of HSP-27 in breast cancer are likely to affect the prognosis and outcome of treatment. In this study, we show a strong correlation between elevated levels of the Brn-3b POU transcription factor and high levels of HSP-27 protein in manipulated MCF-7 breast cancer cells as well as in human breast biopsies. Conversely, HSP-27 is decreased on loss of Brn-3b. In cotransfection assays, Brn-3b can strongly transactivate the HSP-27 promoter, supporting a role for direct regulation of HSP-27 expression. Brn-3b also cooperates with the estrogen receptor (ER) to facilitate maximal stimulation of the HSP-27 promoter, with significantly enhanced activity of this promoter observed on coexpression of Brn-3b and ER compared with either alone. RNA interference and site-directed mutagenesis support the requirement for the Brn-3b binding site on the HSP-27 promoter, which facilitates maximal transactivation either alone or on interaction with the ER. Chromatin immunoprecipitation provides evidence for association of Brn-3b with the HSP-27 promoter in the intact cell. Thus, Brn-3b can, directly and indirectly (via interaction with the ER), activate HSP-27 expression, and this may represent one mechanism by which Brn-3b mediates its effects in breast cancer cells.

Differential nuclease sensitivity profiling of chromatin reveals biochemical footprints coupled to gene expression and functional DNA elements in maize.

This work presents a technique to identify specific regions of maize chromatin that are hypersensitive to digestion by micrococcal nuclease, which preferentially cleaves DNA not bound into nucleosomes. This chromatin profiling technique helps elucidate how global changes in gene expression relate to changes in nucleosome position and should be broadly applicable to other species. The eukaryotic genome is organized into nucleosomes, the fundamental units of chromatin. The positions of nucleosomes on DNA regulate protein-DNA interactions and in turn influence DNA-templated events. Despite the increasing number of genome-wide maps of nucleosome position, how global changes in gene expression relate to changes in nucleosome position is poorly understood. We show that in nucleosome occupancy mapping experiments in maize (Zea mays), particular genomic regions are highly susceptible to variation introduced by differences in the extent to which chromatin is digested with micrococcal nuclease (MNase). We exploited this digestion-linked variation to identify protein footprints that are hypersensitive to MNase digestion, an approach we term differential nuclease sensitivity profiling (DNS-chip). Hypersensitive footprints were enriched at the 5′ and 3′ ends of genes, associated with gene expression levels, and significantly overlapped with conserved noncoding sequences and the binding sites of the transcription factor KNOTTED1. We also found that the tissue-specific regulation of gene expression was linked to tissue-specific hypersensitive footprints. These results reveal biochemical features of nucleosome organization that correlate with gene expression levels and colocalize with functional DNA elements. This approach to chromatin profiling should be broadly applicable to other species and should shed light on the relationships among chromatin organization, protein-DNA interactions, and genome regulation.

Genome-Wide Prediction of Nucleosome Occupancy in Maize Reveals Plant Chromatin Structural Features at Genes and Other Elements at Multiple Scales

Genome-wide prediction of nucleosome occupancy in maize was performed using a computational approach, and the predictions were validated using biochemical data and reveal novel features of maize promoters. The nucleosome is a fundamental structural and functional chromatin unit that affects nearly all DNA-templated events in eukaryotic genomes. It is also a biochemical substrate for higher order, cis-acting gene expression codes and the monomeric structural unit for chromatin packaging at multiple scales. To predict the nucleosome landscape of a model plant genome, we used a support vector machine computational algorithm trained on human chromatin to predict the nucleosome occupancy likelihood (NOL) across the maize (Zea mays) genome. Experimentally validated NOL plots provide a novel genomic annotation that highlights gene structures, repetitive elements, and chromosome-scale domains likely to reflect regional gene density. We established a new genome browser (http://www.genomaize.org) for viewing support vector machine-based NOL scores. This annotation provides sequence-based comprehensive coverage across the entire genome, including repetitive genomic regions typically excluded from experimental genomics data. We find that transposable elements often displayed family-specific NOL profiles that included distinct regions, especially near their termini, predicted to have strong affinities for nucleosomes. We examined transcription start site consensus NOL plots for maize gene sets and discovered that most maize genes display a typical +1 nucleosome positioning signal just downstream of the start site but not upstream. This overall lack of a –1 nucleosome positioning signal was also predicted by our method for Arabidopsis (Arabidopsis thaliana) genes and verified by additional analysis of previously published Arabidopsis MNase-Seq data, revealing a general feature of plant promoters. Our study advances plant chromatin research by defining the potential contribution of the DNA sequence to observed nucleosome positioning and provides an invariant baseline annotation against which other genomic data can be compared.

BAF250a Protein Regulates Nucleosome Occupancy and Histone Modifications in Priming Embryonic Stem Cell Differentiation

Background: How BAF250a regulates nucleosome configuration in ES cells is not clear. Results: BAF250a regulates nucleosome occupancy and H3K27me3 to control gene expression during ES cell differentiation. Conclusion: BAF250a plays a key role in poised chromatin regulation. Significance: Understanding the mechanisms of chromatin remodeling in poised chromatin regulation provides epigenetic insights into ES cell differentiation. The unique chromatin signature of ES cells is fundamental to the pluripotency and differentiation of ES cells. One key feature is the poised chromatin state of master developmental genes that are transcriptionally repressed in ES cells but ready to be activated in response to differentiation signals. Poised chromatin in ES cells contains both H3 Lys-4 trimethylation (H3K4me3) and H3 Lys-27 trimethylation (H3K27me3) methylation, indicating activating and repressing potential. However, the contribution of non-covalent chromatin structure to the poised state is not well understood. To address whether remodeling of nucleosomes is important to the poised state, we characterized the function of BAF250a, a key regulatory subunit of the ES cell ATP-dependent Brahma-associated factor (BAF) chromatin remodeling complex (esBAF). Acute deletion of BAF250a disrupted the differentiation potential of ES cells by altering the expression timing of key developmental genes and pluripotent genes. Our genome-wide nucleosome and histone modification analyses indicated that the disruption of gene expression timing was largely due to changes of chromatin structures at poised genes, particularly those key developmental genes mediated by BAF250a. Specifically, BAF250a deletion caused a nucleosome occupancy increase at H3K4me3- and/or H3K27me3-associated promoters. Moreover, H3K27me3 levels and the number of bivalent promoter genes were reduced in BAF250a KO ES cells. We revealed that BAF250a ablation led to elevated Brg1 but reduced Suz12 recruitment at nucleosome occupancy-increased regions, indicating an unexpected and complicated role of BAF250a in regulating esBAF and Polycomb repressive complex (PRC) activities. Together, our studies identified that BAF250a mediates esBAF and PRC functions to establish the poised chromatin configuration in ES cells, which is essential for the proper differentiation of ES cells.

Chromatin patterns associated with lung adenocarcinoma progression

The development and progression of lung adenocarcinoma, one of the most common cancers, is driven by the interplay of genetic and epigenetic changes and the role of chromatin structure in malignant transformation remains poorly understood. We used systematic nucleosome distribution and chromatin accessibility microarray mapping platforms to analyze the genome-wide chromatin structure from normal tissues and from primary lung adenocarcinoma of different grades and stages. We identified chromatin-based patterns across different patients with lung adenocarcinoma of different cancer grade and stage. Low-grade cancers had nucleosome distributions very different compared with the corresponding normal tissue but had nearly identical chromatin accessibility. Conversely, nucleosome distributions of high-grade cancers showed few differences. Substantial disruptions in chromosomal accessibility were seen in a patient with a high-grade and high-stage tumor. These data imply that chromatin structure changes during the progression of lung adenocarcinoma. We have therefore developed a model in which low-grade lung adenocarcinomas are linked to changes in nucleosome distributions, whereas higher-grade tumors are linked to large-scale chromosomal changes. These results provide a foundation for the development of a comprehensive framework linking the general and locus-specific roles of chromatin structure to lung cancer progression. We propose that this strategy has the potential to identify a new class of chromatin-based diagnostic, prognostic and therapeutic markers in cancer progression.

Regulated large-scale nucleosome density patterns and precise nucleosome positioning correlate with V(D)J recombination

Significance Immunoglobulin and T-cell receptor genes are assembled in lymphoid cells from gene fragments by the process known as V(D)J recombination, which is initiated by the recombination activating gene (RAG)1/RAG2 recombinase. To ensure that recombination occurs only in the correct cell type and at the right developmental stage, multiple layers of regulation are necessary, including specific modifications of chromatin. We show that nucleosome positioning is another important factor in this regulation. Developmentally regulated changes in nucleosome positioning help to guide RAG1/RAG2 to the correct sites in recombinationally active cells. These changes occur on the scale of hundreds of kilobases, a form of regulation not typically seen in the rest of the mammalian genome. We show that the physical distribution of nucleosomes at antigen receptor loci is subject to regulated cell type-specific and lineage-specific positioning and correlates with the accessibility of these gene segments to recombination. At the Ig heavy chain locus (IgH), a nucleosome in pro-B cells is generally positioned over each IgH variable (VH) coding segment, directly adjacent to the recombination signal sequence (RSS), placing the RSS in a position accessible to the recombination activating gene (RAG) recombinase. These changes result in establishment of a specific chromatin organization at the RSS that facilitates accessibility of the genomic DNA for the RAG recombinase. In contrast, in mouse embryonic fibroblasts the coding segment is depleted of nucleosomes, which instead cover the RSS, thereby rendering it inaccessible. Pro-T cells exhibit a pattern intermediate between pro-B cells and mouse embryonic fibroblasts. We also find large-scale variations of nucleosome density over hundreds of kilobases, delineating chromosomal domains within IgH, in a cell type-dependent manner. These findings suggest that developmentally regulated changes in nucleosome location and occupancy, in addition to the known chromatin modifications, play a fundamental role in regulating V(D)J recombination. Nucleosome positioning—which has previously been observed to vary locally at individual enhancers and promoters—may be a more general mechanism by which cells can regulate the accessibility of the genome during development, at scales ranging from several hundred base pairs to many kilobases.