Independent and complementary methods for large-scale structural analysis of mammalian chromatin

The chromatin remodeling protein Lsh alters nucleosome occupancy at putative enhancers and modulates binding of lineage specific transcription factors

Dynamic regulation of chromatin accessibility is a key feature of cellular differentiation during embryogenesis, but the precise factors that control access to chromatin remain largely unknown. Lsh/HELLS is critical for normal development and mutations of Lsh in human cause the ICF (Immune deficiency, Centromeric instability, Facial anomalies) syndrome, a severe immune disorder with multiple organ deficiencies. We report here that Lsh, previously known to regulate DNA methylation level, has a genome wide chromatin remodeling function. Using micrococcal nuclease (MNase)-seq analysis, we demonstrate that Lsh protects MNase accessibility at transcriptional regulatory regions characterized by DNase I hypersensitivity and certain histone 3 (H3) tail modifications associated with enhancers. Using an auxin-inducible degron system, allowing proteolytical degradation of Lsh, we show that Lsh mediated changes in nucleosome occupancy are independent of DNA methylation level and are characterized by reduced H3 occupancy. While Lsh mediated nucleosome occupancy prevents binding sites for transcription factors in wild type cells, depletion of Lsh leads to an increase in binding of ectopically expressed tissue specific transcription factors to their respective binding sites. Our data suggests that Lsh mediated chromatin remodeling can modulate nucleosome positioning at a subset of putative enhancers contributing to the preservation of cellular identity through regulation of accessibility.

Comprehensive nucleosome mapping of the human genome in cancer progression

Altered chromatin structure is a hallmark of cancer, and inappropriate regulation of chromatin structure may represent the origin of transformation. Important studies have mapped human nucleosome distributions genome wide, but the role of chromatin structure in cancer progression has not been addressed. We developed a MNase-Transcription Start Site Sequence Capture method (mTSS-seq) to map the nucleosome distribution at human transcription start sites genome-wide in primary human lung and colon adenocarcinoma tissue. Here, we confirm that nucleosome redistribution is an early, widespread event in lung (LAC) and colon (CRC) adenocarcinoma. These altered nucleosome architectures are consistent between LAC and CRC patient samples indicating that they may serve as important early adenocarcinoma markers. We demonstrate that the nucleosome alterations are driven by the underlying DNA sequence and potentiate transcription factor binding. We conclude that DNA-directed nucleosome redistributions are widespread early in cancer progression. We have proposed an entirely new hierarchical model for chromatin-mediated genome regulation.

The role of Cockayne Syndrome Protein B in transcription regulation

We investigated the question if CSB (Cockayne Syndrome complementation B) protein actively regulates gene transcription and how mutations in CSB gene affect that regulatory role.

Here we describe how we processed and interpreted ChIP-seq data (deposited in Gene Expression Omnibus with accession number GSE50171) obtained during an investigation of that question, and how this analysis assisted in the generation of hypothesis that were subsequently validated using other types of experiment.

BAF250a Protein Regulates Nucleosome Occupancy and Histone Modifications in Priming Embryonic Stem 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.

DNA-Encoded Chromatin Structural Intron Boundary Signals Identify Conserved Genes with Common Function

The regulation of metazoan gene expression occurs in part by pre-mRNA splicing into mature RNAs. Signals affecting the efficiency and specificity with which introns are removed have not been completely elucidated. Splicing likely occurs cotranscriptionally, with chromatin structure playing a key regulatory role. We calculated DNA encoded nucleosome occupancy likelihood (NOL) scores at the boundaries between introns and exons across five metazoan species. We found that (i) NOL scores reveal a sequence-based feature at the introns on both sides of the intron-exon boundary; (ii) this feature is not part of any recognizable consensus sequence; (iii) this feature is conserved throughout metazoa; (iv) this feature is enriched in genes sharing similar functions: ATPase activity, ATP binding, helicase activity, and motor activity; (v) genes with these functions exhibit different genomic characteristics; (vi) in vivo nucleosome positioning data confirm ontological enrichment at this feature; and (vii) genes with this feature exhibit unique dinucleotide distributions at the intron-exon boundary. The NOL scores point toward a physical property of DNA that may play a role in the mechanism of pre-mRNA splicing. These results provide a foundation for identification of a new set of regulatory DNA elements involved in splicing regulation.

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

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.

Nucleosomal occupancy changes locally over key regulatory regions during cell differentiation and reprogramming

Chromatin structure determines DNA accessibility. We compare nucleosome occupancy in mouse and human embryonic stem cells (ESCs), induced-pluripotent stem cells (iPSCs) and differentiated cell types using MNase-seq. To address variability inherent in this technique, we developed a bioinformatic approach to identify regions of difference (RoD) in nucleosome occupancy between pluripotent and somatic cells. Surprisingly, most chromatin remains unchanged; a majority of rearrangements appear to affect a single nucleosome. RoDs are enriched at genes and regulatory elements, including enhancers associated with pluripotency and differentiation. RoDs co-localize with binding sites of key developmental regulators, including the reprogramming factors Klf4, Oct4/Sox2 and c-Myc. Nucleosomal landscapes in ESC enhancers are extensively altered, exhibiting lower nucleosome occupancy in pluripotent cells than in somatic cells. Most changes are reset during reprogramming. We conclude that changes in nucleosome occupancy are a hallmark of cell differentiation and reprogramming and likely identify regulatory regions essential for these processes.

Changes in chromatin structure impact gene expression programs by modulating accessibility to the transcription machinery. Here, West et al. explore differences in nucleosome occupancy between mammalian pluripotent and somatic cells and uncover regulatory regions likely to play key roles in determining cell identity.

The sequence-specific transcription factor c-Jun targets Cockayne syndrome protein B to regulate transcription and chromatin structure

Cockayne syndrome is an inherited premature aging disease associated with numerous developmental and neurological defects, and mutations in the gene encoding the CSB protein account for the majority of Cockayne syndrome cases. Accumulating evidence suggests that CSB functions in transcription regulation, in addition to its roles in DNA repair, and those defects in this transcriptional activity might contribute to the clinical features of Cockayne syndrome. Transcription profiling studies have so far uncovered CSB-dependent effects on gene expression; however, the direct targets of CSB's transcriptional activity remain largely unknown. In this paper, we report the first comprehensive analysis of CSB genomic occupancy during replicative cell growth. We found that CSB occupancy sites display a high correlation to regions with epigenetic features of promoters and enhancers. Furthermore, we found that CSB occupancy is enriched at sites containing the TPA-response element. Consistent with this binding site preference, we show that CSB and the transcription factor c-Jun can be found in the same protein-DNA complex, suggesting that c-Jun can target CSB to specific genomic regions. In support of this notion, we observed decreased CSB occupancy of TPA-response elements when c-Jun levels were diminished. By modulating CSB abundance, we found that CSB can influence the expression of nearby genes and impact nucleosome positioning in the vicinity of its binding site. These results indicate that CSB can be targeted to specific genomic loci by sequence-specific transcription factors to regulate transcription and local chromatin structure. Additionally, comparison of CSB occupancy sites with the MSigDB Pathways database suggests that CSB might function in peroxisome proliferation, EGF receptor transactivation, G protein signaling and NF-κB activation, shedding new light on the possible causes and mechanisms of Cockayne syndrome.

Cockayne syndrome is a devastating inherited disease, in which patients appear to age prematurely, have sun sensitivity and suffer from profound neurological and developmental defects. Mutations in the CSB gene account for the majority of Cockayne syndrome cases. CSB is an ATP-dependent chromatin remodeler, and these proteins can use energy from ATP-hydrolysis to alter contacts between DNA and histones of a nucleosome, the basic units of chromatin structure. CSB functions in DNA repair, but accumulating evidence reveals that CSB also functions in transcription regulation. Here, we determined the genomic localization of CSB to identify its gene targets and found that CSB occupancy displays high correlation to regions with epigenetic features of promoters and enhancers. Furthermore, CSB is enriched at genomic regions containing the binding site for the c-Jun transcription factor, and we found that these two proteins interact, uncovering a new targeting mechanism for CSB. We also demonstrate that CSB can influence gene expression in the vicinity of its binding sites and alter local chromatin structure. Together, this study supports the hypothesis that defects in the regulation of gene expression and chromatin structure by CSB might contribute to the diverse clinical features of Cockayne syndrome.

Genome-wide nucleosome map and cytosine methylation levels of an ancient human genome

Epigenetic information is available from contemporary organisms, but is difficult to track back in evolutionary time. Here, we show that genome-wide epigenetic information can be gathered directly from next-generation sequence reads of DNA isolated from ancient remains. Using the genome sequence data generated from hair shafts of a 4000-yr-old Paleo-Eskimo belonging to the Saqqaq culture, we generate the first ancient nucleosome map coupled with a genome-wide survey of cytosine methylation levels. The validity of both nucleosome map and methylation levels were confirmed by the recovery of the expected signals at promoter regions, exon/intron boundaries, and CTCF sites. The top-scoring nucleosome calls revealed distinct DNA positioning biases, attesting to nucleotide-level accuracy. The ancient methylation levels exhibited high conservation over time, clustering closely with modern hair tissues. Using ancient methylation information, we estimated the age at death of the Saqqaq individual and illustrate how epigenetic information can be used to infer ancient gene expression. Similar epigenetic signatures were found in other fossil material, such as 110,000- to 130,000-yr-old bones, supporting the contention that ancient epigenomic information can be reconstructed from a deep past. Our findings lay the foundation for extracting epigenomic information from ancient samples, allowing shifts in epialleles to be tracked through evolutionary time, as well as providing an original window into modern epigenomics.

Genome-wide prediction of nucleosome occupancy in maize reveals plant chromatin structural features at genes and other elements at multiple scales

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.

Dynamic nucleosome organization at hox promoters during zebrafish embryogenesis

Nucleosome organization at promoter regions plays an important role in regulating gene activity. Genome-wide studies in yeast, flies, worms, mammalian embryonic stem cells and transformed cell lines have found well-positioned nucleosomes flanking a nucleosome depleted region (NDR) at transcription start sites. This nucleosome arrangement depends on DNA sequence (cis-elements) as well as DNA binding factors and ATP-dependent chromatin modifiers (trans-factors). However, little is understood about how the nascent embryonic genome positions nucleosomes during development. This is particularly intriguing since the embryonic genome must undergo a broad reprogramming event upon fusion of sperm and oocyte. Using four stages of early embryonic zebrafish development, we map nucleosome positions at the promoter region of 37 zebrafish hox genes. We find that nucleosome arrangement at the hox promoters is a progressive process that takes place over several stages. At stages immediately after fertilization, nucleosomes appear to be largely disordered at hox promoter regions. At stages after activation of the embryonic genome, nucleosomes are detectable at hox promoters, with positions becoming more uniform and more highly occupied. Since the genomic sequence is invariant during embryogenesis, this progressive change in nucleosome arrangement suggests that trans-factors play an important role in organizing nucleosomes during embryogenesis. Separating hox genes into expressed and non-expressed groups shows that expressed promoters have better positioned and occupied nucleosomes, as well as distinct NDRs, than non-expressed promoters. Finally, by blocking the retinoic acid-signaling pathway, we disrupt early hox gene transcription, but observe no effect on nucleosome positions, suggesting that active hox transcription is not a driving force behind the arrangement of nucleosomes at the promoters of hox genes during early development.

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.

The roles of the monomer length and nucleotide context of plant tandem repeats in nucleosome positioning

Similar to regularly spaced nucleosomes in chromatin, long tandem DNA arrays are composed of regularly alternating monomers that have almost identical primary DNA structures. Such a similarity in the structural organization makes these arrays especially interesting for studying the role of intrinsic DNA preferences in nucleosome positioning. We have studied the nucleosome formation potential of DNA tandem repeat families with different monomer lengths (ML). In total, 165 plant tandem repeat families from the PlantSat database (http://w3lamc.umbr.cas.cz/PlantSat/) were divided into two classes based on the number of nucleosome repeats in one DNA monomer. For predicting nucleosome formation potential, we developed the Phase method, which combines the advantages of multiple bioinformatics models. The Phase method was able to distinguish interfamily differences and intrafamily monomer variation and identify the influence of nucleotide context on nucleosome formation potential. Three main types of nucleosome arrangement in DNA tandem repeat arrays--regular, partially regular (partial), and flexible--were distinguished among a great variety of Phase profiles. The regular type, in which all nucleosomes of the monomer array are positioned in a context-dependent manner, is the most representative type of the class 1 families, with ML equal to or a multiple of the nucleosome repeat length (NRL). In the partially regular type, nucleotide context influences the positioning of only a subset of nucleosomes. The influence of the nucleotide context on nucleosome positioning has the least effect in the flexible type, which contains the greatest number of families (65). The majority of these families belong to class 2 and have nonmultiple ML to NRL ratios.

Differences in genomic architecture between two distinct geographical strains of the blood fluke Schistosoma japonicum reveal potential phenotype basis

The Chinese (SjC) and Philippine (SjP) strains of the blood fluke Schistosoma japonicum have been shown to present clearly different phenotypes in fecundity, pathology, drug sensitivity and immunology. We used microarray based comparative genomic hybridisation (aCGH) to investigate structural differences in the genomes of the two strains and identified seven distinct regions of the S. japonicum genome that present differential aCGH representing either deletion or duplication regions in SjP. Within these regions, genes predicted to be associated with the recognised phenotypic differences were identified and that may provide new insights into the biology and evolution of the two strains, with implications for the epidemiology and control of schistosomiasis japonica in China and the Philippines.

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.

Characterization of enhancer function from genome-wide analyses

There has been a recent surge in the use of genome-wide methodologies to identify and annotate the transcriptional regulatory elements in the human genome. Here we review some of these methodologies and the conceptual insights about transcription regulation that have been gained from the use of genome-wide studies. It has become clear that the binding of transcription factors is itself a highly regulated process, and binding does not always appear to have functional consequences. Numerous properties have now been associated with regulatory elements that may be useful in their identification. Several aspects of enhancer function have been shown to be more widespread than was previously appreciated, including the highly combinatorial nature of transcription factor binding, the postinitiation regulation of many target genes, and the binding of enhancers at early stages to maintain their competence during development. Going forward, the integration of multiple genome-wide data sets should become a standard approach to elucidate higher-order regulatory interactions.

Epigenetic changes in fetal hypothalamic energy regulating pathways are associated with maternal undernutrition and twinning

Undernutrition during pregnancy is implicated in the programming of offspring for the development of obesity and diabetes. We hypothesized that maternal programming causes epigenetic changes in fetal hypothalamic pathways regulating metabolism. This study used sheep to examine the effect of moderate maternal undernutrition (60 d before to 30 d after mating) and twinning to investigate changes in the key metabolic regulators proopiomelanocortin (POMC) and the glucocorticoid receptor (GR) in fetal hypothalami. Methylation of the fetal hypothalamic POMC promoter was reduced in underfed singleton, fed twin, and underfed twin groups (60, 73, and 63% decrease, respectively). This was associated with reduced DNA methyltransferase activity and altered histone methylation and acetylation. Methylation of the hypothalamic GR promoter was decreased in both twin groups and in maternally underfed singleton fetuses (52, 65, and 55% decrease, respectively). This correlated with changes in histone methylation and acetylation and increased GR mRNA expression in the maternally underfed singleton group. Alterations in GR were hypothalamic specific, with no changes in hippocampi. Unaltered levels of OCT4 promoter methylation indicated gene-specific effects. In conclusion, twinning and periconceptional undernutrition are associated with epigenetic changes in fetal hypothalamic POMC and GR genes, potentially resulting in altered energy balance regulation in the offspring.

Basic procedures for epigenetic analysis in plant cell and tissue culture

In vitro culture is one of the most studied techniques, and it is used to study many developmental processes, especially in forestry species, because of growth timing and easy manipulation. Epigenetics has been shown as an important influence on many research analyses such as cancer in mammals and developmental processes in plants such as flowering, but regarding in vitro culture, techniques to study DNA methylation or chromatin modifications were mainly limited to identify somaclonal variation of the micropropagated material. Because in vitro culture is not only a way to generate plant material but also a bunch of differentially induced developmental processes, an approach of techniques and some research carried out to study the different changes regarding DNA methylation and chromatin and translational modifications that take place during these processes is reviewed.

Sequence-based classification using discriminatory motif feature selection

Most existing methods for sequence-based classification use exhaustive feature generation, employing, for example, all -mer patterns. The motivation behind such (enumerative) approaches is to minimize the potential for overlooking important features. However, there are shortcomings to this strategy. First, practical constraints limit the scope of exhaustive feature generation to patterns of length , such that potentially important, longer () predictors are not considered. Second, features so generated exhibit strong dependencies, which can complicate understanding of derived classification rules. Third, and most importantly, numerous irrelevant features are created. These concerns can compromise prediction and interpretation. While remedies have been proposed, they tend to be problem-specific and not broadly applicable. Here, we develop a generally applicable methodology, and an attendant software pipeline, that is predicated on discriminatory motif finding. In addition to the traditional training and validation partitions, our framework entails a third level of data partitioning, a discovery partition. A discriminatory motif finder is used on sequences and associated class labels in the discovery partition to yield a (small) set of features. These features are then used as inputs to a classifier in the training partition. Finally, performance assessment occurs on the validation partition. Important attributes of our approach are its modularity (any discriminatory motif finder and any classifier can be deployed) and its universality (all data, including sequences that are unaligned and/or of unequal length, can be accommodated). We illustrate our approach on two nucleosome occupancy datasets and a protein solubility dataset, previously analyzed using enumerative feature generation. Our method achieves excellent performance results, with and without optimization of classifier tuning parameters. A Python pipeline implementing the approach is available at http://www.epibiostat.ucsf.edu/biostat/sen/dmfs/.

Nucleosome positioning in Saccharomyces cerevisiae

The DNA of eukaryotic cells is spooled around large histone protein complexes, forming nucleosomes that make up the basis for a high-order packaging structure called chromatin. Compared to naked DNA, nucleosomal DNA is less accessible to regulatory proteins and regulatory processes. The exact positions of nucleosomes therefore influence several cellular processes, including gene expression, chromosome segregation, recombination, replication, and DNA repair. Here, we review recent technological advances enabling the genome-wide mapping of nucleosome positions in the model eukaryote Saccharomyces cerevisiae. We discuss the various parameters that determine nucleosome positioning in vivo, including cis factors like AT content, variable tandem repeats, and poly(dA:dT) tracts that function as chromatin barriers and trans factors such as chromatin remodeling complexes, transcription factors, histone-modifying enzymes, and RNA polymerases. In the last section, we review the biological role of chromatin in gene transcription, the evolution of gene regulation, and epigenetic phenomena.

Multiple distinct stimuli increase measured nucleosome occupancy around human promoters

Nucleosomes can block access to transcription factors. Thus the precise localization of nucleosomes relative to transcription start sites and other factor binding sites is expected to be a critical component of transcriptional regulation. Recently developed microarray approaches have allowed the rapid mapping of nucleosome positions over hundreds of kilobases (kb) of human genomic DNA, although these approaches have not yet been widely used to measure chromatin changes associated with changes in transcription. Here, we use custom tiling microarrays to reveal changes in nucleosome positions and abundance that occur when hormone-bound glucocorticoid receptor (GR) binds to sites near target gene promoters in human osteosarcoma cells. The most striking change is an increase in measured nucleosome occupancy at sites spanning ∼1 kb upstream and downstream of transcription start sites, which occurs one hour after addition of hormone, but is lost at 4 hours. Unexpectedly, this increase was seen both on GR-regulated and GR-non-regulated genes. In addition, the human SWI/SNF chromatin remodeling factor (a GR co-activator) was found to be important for increased occupancy upon hormone treatment and also for low nucleosome occupancy without hormone. Most surprisingly, similar increases in nucleosome occupancy were also seen on both regulated and non-regulated promoters during differentiation of human myeloid leukemia cells and upon activation of human CD4+ T-cells. These results indicate that dramatic changes in chromatin structure over ∼2 kb of human promoters may occur genomewide and in response to a variety of stimuli, and suggest novel models for transcriptional regulation.

Characterization of DXZ4 conservation in primates implies important functional roles for CTCF binding, array expression and tandem repeat organization on the X chromosome

BACKGROUND

Comparative sequence analysis is a powerful means with which to identify functionally relevant non-coding DNA elements through conserved nucleotide sequence. The macrosatellite DXZ4 is a polymorphic, uninterrupted, tandem array of 3-kb repeat units located exclusively on the human X chromosome. While not obviously protein coding, its chromatin organization suggests differing roles for the array on the active and inactive X chromosomes.

RESULTS

In order to identify important elements within DXZ4, we explored preservation of DNA sequence and chromatin conformation of the macrosatellite in primates. We found that DXZ4 DNA sequence conservation beyond New World monkeys is limited to the promoter and CTCF binding site, although DXZ4 remains a GC-rich tandem array. Investigation of chromatin organization in macaques revealed that DXZ4 in males and on the active X chromosome is packaged into heterochromatin, whereas on the inactive X, DXZ4 was euchromatic and bound by CTCF.

CONCLUSIONS

Collectively, these data suggest an important conserved role for DXZ4 on the X chromosome involving expression, CTCF binding and tandem organization.

Genome-wide maps of mononucleosomes and dinucleosomes containing hyperacetylated histones of Aspergillus fumigatus

It is suggested that histone modifications and/or histone variants influence the nucleosomal DNA length. We sequenced both ends of mononucleosomal and dinucleosomal DNA fragments of the filamentous fungus Aspergillus fumigatus, after treatment with the histone deacetylase inhibitor trichostatin A (TSA). After mapping the DNA fragments to the genome, we identified >7 million mononucleosome positions and >7 million dinucleosome positions. We showed that the distributions of the lengths of the mononucleosomal DNA fragments after 15-min and 30-min treatments with micrococcal nuclease (MNase) showed a single peak at 168 nt and 160 nt, respectively. The distributions of the lengths of the dinucleosomal DNA fragments after 15-min- and 30-min-treatment with MNase showed a single peak at 321 nt and 306 nt, respectively. The nucleosomal DNA fragments obtained from the TSA-treated cells were significantly longer than those obtained from the untreated cells. On the other hand, most of the genes did not undergo any change after treatment. Between the TSA-treated and untreated cells, only 77 genes had >or=2-fold change in expression levels. In addition, our results showed that the locations where mononucleosomes were frequently detected were conserved between the TSA-treated cells and untreated cells in the gene promoters (lower density of the nucleosomes). However, these locations were less conserved in the bodies (higher density of the nucleosomes) of genes with >or=2-fold changes. Our findings indicate that TSA influences the nucleosome positions, especially of the regions with high density of the nucleosomes by elongation of the nucleosomal DNA. However, most of the nucleosome positions are conserved in the gene promoters, even after treatment with TSA, because of the low density of nucleosomes in the gene promoters.

A region of the human HOXD cluster that confers polycomb-group responsiveness

Polycomb group (PcG) proteins are essential for accurate axial body patterning during embryonic development. PcG-mediated repression is conserved in metazoans and is targeted in Drosophila by Polycomb response elements (PREs). However, targeting sequences in humans have not been described. While analyzing chromatin architecture in the context of human embryonic stem cell (hESC) differentiation, we discovered a 1.8kb region between HOXD11 and HOXD12 (D11.12) that is associated with PcG proteins, becomes nuclease hypersensitive, and then shows alteration in nuclease sensitivity as hESCs differentiate. The D11.12 element repressed luciferase expression from a reporter construct and full repression required a highly conserved region and YY1 binding sites. Furthermore, repression was dependent on the PcG proteins BMI1 and EED and a YY1-interacting partner, RYBP. We conclude that D11.12 is a Polycomb-dependent regulatory region with similarities to Drosophila PREs, indicating conservation in the mechanisms that target PcG function in mammals and flies.

The organization of nucleosomes around splice sites

The occupancy of nucleosomes along chromosome is a key factor for gene regulation. However, except promoter regions, genome-wide properties and functions of nucleosome organization remain unclear in mammalian genomes. Using the computational model of Increment of Diversity with Quadratic Discriminant (IDQD) trained from the microarray data, the nucleosome occupancy score (NOScore) was defined and applied to splice junction regions of constitutive, cassette exon, alternative 3′ and 5′ splicing events in the human genome. We found an interesting relation between NOScore and RNA splicing: exon regions have higher NOScores compared with their flanking intron sequences in both constitutive and alternative splicing events, indicating the stronger nucleosome occupation potential of exon regions. In addition, NOScore valleys present at ∼25 bp upstream of the acceptor site in all splicing events. By defining folding diversity-to-energy ratio to describe RNA structural flexibility, we demonstrated that primary RNA transcripts from nucleosome occupancy regions are relatively rigid and those from nucleosome depleted regions are relatively flexible. The negative correlation between nucleosome occupation/depletion of DNA sequence and structural flexibility/rigidity of its primary transcript around splice junctions may provide clues to the deeper understanding of the unexpected role for nucleosome organization in the regulation of RNA splicing.

Epigenetic silencing in Friedreich ataxia is associated with depletion of CTCF (CCCTC-binding factor) and antisense transcription

BACKGROUND

Over 15 inherited diseases are caused by expansion of triplet-repeats. Friedreich ataxia (FRDA) patients are homozygous for an expanded GAA triplet-repeat sequence in intron 1 of the FXN gene. The expanded GAA triplet-repeat results in deficiency of FXN gene transcription, which is reversed via administration of histone deacetylase inhibitors indicating that transcriptional silencing is at least partially due to an epigenetic abnormality.

METHODOLOGY/PRINCIPAL FINDINGS

We found a severe depletion of the chromatin insulator protein CTCF (CCCTC-binding factor) in the 5'UTR of the FXN gene in FRDA, and coincident heterochromatin formation involving the +1 nucleosome via enrichment of H3K9me3 and recruitment of heterochromatin protein 1. We identified FAST-1 (FXNAntisense Transcript - 1), a novel antisense transcript that overlaps the CTCF binding site in the 5'UTR, which was expressed at higher levels in FRDA. The reciprocal relationship of deficient FXN transcript and higher levels of FAST-1 seen in FRDA was reproduced in normal cells via knockdown of CTCF.

CONCLUSIONS/SIGNIFICANCE

CTCF depletion constitutes an epigenetic switch that results in increased antisense transcription, heterochromatin formation and transcriptional deficiency in FRDA. These findings provide a mechanistic basis for the transcriptional silencing of the FXN gene in FRDA, and broaden our understanding of disease pathogenesis in triplet-repeat diseases.

Genome-wide maps of mono- and di-nucleosomes of Aspergillus fumigatus

We identified 5 975 567 mono- and 6 995 122 di-nucleosome positions of the fungus Aspergillus fumigatus, which was detected at high resolution based on the DNA sequence data obtained from both mono- and di-nucleosomal DNA fragments. We show that the distribution of lengths of the mono-nucleosomal DNA fragments has two peaks at 135 and 150 nt, whereas the distribution of di-nucleosomal DNA fragment lengths has a single peak at 285 nt. Although the gene bodies of the active and inactive genes and the inactive gene promoters had the two peaks of the mono-nucleosomal DNA fragment lengths, the active gene promoter lost the longer peak at 150 nt. Our findings strongly suggest that the nucleosomes protecting longer DNA fragments against MNase at the promoters, thereby inhibiting high gene expression.

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.

Inside the nucleus, DNA is wrapped into a complex molecular structure called chromatin, whose fundamental unit is ∼150 bp of DNA organized around the eight-histone protein complex known as the nucleosome. Understanding the local organization of nucleosomes is critical for understanding how chromatin impacts gene regulation. Here, we describe a computational model that predicts nucleosome placement from DNA sequence. We train the model using data derived from human cell lines, and we apply the model systematically to 1% of the human genome. We show that previously described models trained from yeast data correlate strongly with the human-trained model, suggesting a common mechanism for sequence-based determination of nucleosome occupancy. In addition, we observe a striking complementarity between models trained using data from weakly and strongly digested samples: one type of model recognizes nucleosome-free regions, whereas the other identifies well-positioned nucleosomes. Finally, our analysis of predicted nucleosome positions in the human genome allows us to identify common features of nucleosome-forming and inhibitory sequences. Overall, our results are consistent with the classical statistical positioning theory of nucleosome organization.

Dynamic regulation of nucleosome positioning in the human genome

The positioning of nucleosomes with respect to DNA plays an important role in regulating transcription. However, nucleosome mapping has been performed for only limited genomic regions in humans. We have generated genome-wide maps of nucleosome positions in both resting and activated human CD4+ T cells by direct sequencing of nucleosome ends using the Solexa high-throughput sequencing technique. We find that nucleosome phasing relative to the transcription start sites is directly correlated to RNA polymerase II (Pol II) binding. Furthermore, the first nucleosome downstream of a start site exhibits differential positioning in active and silent genes. TCR signaling induces extensive nucleosome reorganization in promoters and enhancers to allow transcriptional activation or repression. Our results suggest that H2A.Z-containing and modified nucleosomes are preferentially lost from the -1 nucleosome position. Our data provide a comprehensive view of the nucleosome landscape and its dynamic regulation in the human genome.