Nucleosome positioning signals in the DNA sequence of the human and mouse H19 imprinting control regions

Genome wide nucleosome mapping for HSV-1 shows nucleosomes are deposited at preferred positions during lytic infection

HSV is a large double stranded DNA virus, capable of causing a variety of diseases from the common cold sore to devastating encephalitis. Although DNA within the HSV virion does not contain any histone protein, within 1 h of infecting a cell and entering its nucleus the viral genome acquires some histone protein (nucleosomes). During lytic infection, partial micrococcal nuclease (MNase) digestion does not give the classic ladder band pattern, seen on digestion of cell DNA or latent viral DNA. However, complete digestion does give a mono-nucleosome band, strongly suggesting that there are some nucleosomes present on the viral genome during the lytic infection, but that they are not evenly positioned, with a 200 bp repeat pattern, like cell DNA. Where then are the nucleosomes positioned? Here we perform HSV-1 genome wide nucleosome mapping, at a time when viral replication is in full swing (6 hr PI), using a microarray consisting of 50mer oligonucleotides, covering the whole viral genome (152 kb). Arrays were probed with MNase-protected fragments of DNA from infected cells. Cells were not treated with crosslinking agents, thus we are only mapping tightly bound nucleosomes. The data show that nucleosome deposition is not random. The distribution of signal on the arrays suggest that nucleosomes are located at preferred positions on the genome, and that there are some positions that are not occupied (nucleosome free regions -NFR or Nucleosome depleted regions -NDR), or occupied at frequency below our limit of detection in the population of genomes. Occupancy of only a fraction of the possible sites may explain the lack of a typical MNase partial digestion band ladder pattern for HSV DNA during lytic infection. On average, DNA encoding Immediate Early (IE), Early (E) and Late (L) genes appear to have a similar density of nucleosomes.

Epigenetic control of SPI1 gene by CTCF and ISWI ATPase SMARCA5

CCCTC-binding factor (CTCF) can both activate as well as inhibit transcription by forming chromatin loops between regulatory regions and promoters. In this regard, Ctcf binding on non-methylated DNA and its interaction with the Cohesin complex results in differential regulation of the H19/Igf2 locus. Similarly, a role for CTCF has been established in normal hematopoietic development; however its involvement in leukemia remains elusive. Here, we show that Ctcf binds to the imprinting control region of H19/Igf2 in AML blasts. We also demonstrate that Smarca5, which also associates with the Cohesin complex, facilitates Ctcf binding to its target sites on DNA. Furthermore, Smarca5 supports Ctcf functionally and is needed for enhancer-blocking effect at ICR. We next asked whether CTCF and SMARCA5 control the expression of key hematopoiesis regulators. In normally differentiating myeloid cells both CTCF and SMARCA5 together with members of the Cohesin complex are recruited to the SPI1 gene, a key hematopoiesis regulator and leukemia suppressor. Due to DNA methylation, CTCF binding to the SPI1 gene is blocked in AML blasts. Upon AZA-mediated DNA demethylation of human AML blasts, CTCF and SMARCA5 are recruited to the −14.4 Enhancer of SPI1 gene and block its expression. Our data provide new insight into complex SPI1 gene regulation now involving additional key epigenetic factors, CTCF and SMARCA5 that control PU.1 expression at the −14.4 Enhancer.

A co-ordinated interaction between CTCF and ER in breast cancer cells

Background

CCCTC-binding factor (CTCF) is a conserved zinc finger transcription factor that is involved in both intra- and interchromasomal looping. Recent research has shown a role for CTCF in estrogen receptor (ER) biology, at some individual loci, but a multi-context global analysis of CTCF binding and transcription activity is lacking.

Results

We now map CTCF binding genome wide in breast cancer cells and find that CTCF binding is unchanged in response to estrogen or tamoxifen treatment. We find a small but reproducible set of CTCF binding events that overlap with both the nuclear receptor, estrogen receptor, and the forkhead protein FOXA1. These overlapping binding events are likely functional as they are biased towards estrogen-regulated genes, compared to regions lacking either CTCF or ER binding. In addition we identify cell-line specific CTCF binding events. These binding events are more likely to be associated with cell-line specific ER binding events and are also more likely to be adjacent to genes that are expressed in that particular cell line.

Conclusion

The evolving role for CTCF in ER biology is complex, but is likely to be multifunctional and possibly influenced by the specific genomic locus. Our data suggest a positive, pro-transcriptional role for CTCF in ER-mediated gene expression in breast cancer cells. CTCF not only provides boundaries for accessible and 'protected' transcriptional blocks, but may also influence the actual binding of ER to the chromatin, thereby modulating the estrogen-mediated gene expression changes observed in breast cancer cells.

ArchAlign: coordinate-free chromatin alignment reveals novel architectures

To facilitate identification and characterization of genomic functional elements, we have developed a chromatin architecture alignment algorithm (ArchAlign). ArchAlign identifies shared chromatin structural patterns from high-resolution chromatin structural datasets derived from next-generation sequencing or tiled microarray approaches for user defined regions of interest. We validated ArchAlign using well characterized functional elements, and used it to explore the chromatin structural architecture at CTCF binding sites in the human genome. ArchAlign is freely available at http://www.acsu.buffalo.edu/~mjbuck/ArchAlign.html.

Does CTCF mediate between nuclear organization and gene expression?

The multifunctional zinc-finger protein CCCTC-binding factor (CTCF) is a very strong candidate for the role of coordinating the expression level of coding sequences with their three-dimensional position in the nucleus, apparently responding to a "code" in the DNA itself. Dynamic interactions between chromatin fibers in the context of nuclear architecture have been implicated in various aspects of genome functions. However, the molecular basis of these interactions still remains elusive and is a subject of intense debate. Here we discuss the nature of CTCF-DNA interactions, the CTCF-binding specificity to its binding sites and the relationship between CTCF and chromatin, and we examine data linking CTCF with gene regulation in the three-dimensional nuclear space. We discuss why these features render CTCF a very strong candidate for the role and propose a unifying model, the "CTCF code," explaining the mechanistic basis of how the information encrypted in DNA may be interpreted by CTCF into diverse nuclear functions.

Multiple sequence-directed possibilities provide a pool of nucleosome position choices in different states of activity of a gene

BACKGROUND

Genome-wide mappings of nucleosome occupancy in different species have shown presence of well-positioned nucleosomes. While the DNA sequences may help decide their locations, the observed positions in vivo are end-results of chromatin remodeling, the state of gene activity and binding of the sequence-specific factors to the DNA, all of which influence nucleosome positions. Thus, the observed nucleosome locations in vivo do not reflect the true contribution of DNA sequence to the mapped position. Moreover, the naturally occurring nucleosome-positioning sequences are known to guide multiple translational positionings.

RESULTS

We show that yeast SNR6, a gene transcribed by RNA polymerase III, constitutes nucleosome-positioning sequence. In the absence of a chromatin remodeler or any factor binding, the gene sequence confers a unique rotational phase to nucleosomes in the gene region, and directs assembly of several translationally positioned nucleosomes on approximately 1.2 kb DNA from the gene locus, including the short approximately 250 bp gene region. Mapping of all these gene sequence-directed nucleosome positions revealed that the array of nucleosomes in the gene upstream region occupy the same positions as those observed in vivo but the nucleosomes on the gene region can be arranged in three distinct registers. Two of these arrangements differ from each other in the position of only one nucleosome, and match with the nucleosome positions on the gene in repressed and active states in vivo, where the gene-specific factor is known to occupy the gene in both the states. The two positions are interchanged by an ATP-dependent chromatin remodeler in vivo. The third register represents the positions which block the access of the factor to the gene promoter elements.

CONCLUSION

On a gene locus, multiple nucleosome positions are directed by a gene sequence to provide a pool of possibilities, out of which the preferred ones are selected by the chromatin remodeler and transcription factor of the gene under different states of activity of the gene.

Role of hMOF-dependent histone H4 lysine 16 acetylation in the maintenance of TMS1/ASC gene activity

Epigenetic silencing of tumor suppressor genes in human cancers is associated with aberrant methylation of promoter region CpG islands and local alterations in histone modifications. However, the mechanisms that drive these events remain unclear. Here, we establish an important role for histone H4 lysine 16 acetylation (H4K16Ac) and the histone acetyltransferase hMOF in the regulation of TMS1/ASC, a proapoptotic gene that undergoes epigenetic silencing in human cancers. In the unmethylated and active state, the TMS1 CpG island is spanned by positioned nucleosomes and marked by histone H3K4 methylation. H4K16Ac was uniquely localized to two sharp peaks that flanked the unmethylated CpG island and corresponded to strongly positioned nucleosomes. Aberrant methylation and silencing of TMS1 was accompanied by loss of the H4K16Ac peaks, loss of nucleosome positioning, hypomethylation of H3K4, and hypermethylation of H3K9. In addition, a single peak of histone H4 lysine 20 trimethylation was observed near the transcription start site. Down-regulation of hMOF or another component of the MSL complex resulted in a gene-specific decrease in H4K16Ac, loss of nucleosome positioning, and silencing of TMS1. Gene silencing induced by H4K16 deacetylation occurred independently of changes in histone methylation and DNA methylation and was reversed on hMOF reexpression. These results indicate that the selective marking of nucleosomes flanking the CpG island by hMOF is required to maintain TMS1 gene activity and suggest that the loss of H4K16Ac, mobilization of nucleosomes, and transcriptional down-regulation may be important events in the epigenetic silencing of certain tumor suppressor genes in cancer.

The testis-specific factor CTCFL cooperates with the protein methyltransferase PRMT7 in H19 imprinting control region methylation

Expression of imprinted genes is restricted to a single parental allele as a result of epigenetic regulation—DNA methylation and histone modifications. Igf2/H19 is a reciprocally imprinted locus exhibiting paternal Igf2 and maternal H19 expression. Their expression is regulated by a paternally methylated imprinting control region (ICR) located between the two genes. Although the de novo DNA methyltransferases have been shown to be necessary for the establishment of ICR methylation, the mechanism by which they are targeted to the region remains unknown. We demonstrate that CTCFL/BORIS, a paralog of CTCF, is an ICR-binding protein expressed during embryonic male germ cell development, coinciding with the timing of ICR methylation. PRMT7, a protein arginine methyltransferase with which CTCFL interacts, is also expressed during embryonic testis development. Symmetrical dimethyl arginine 3 of histone H4, a modification catalyzed by PRMT7, accumulates in germ cells during this developmental period. This modified histone is also found enriched in both H19 ICR and Gtl2 differentially methylated region (DMR) chromatin of testis by chromatin immunoprecipitation (ChIP) analysis. In vitro studies demonstrate that CTCFL stimulates the histone-methyltransferase activity of PRMT7 via interactions with both histones and PRMT7. Finally, H19 ICR methylation is demonstrated by nuclear co-injection of expression vectors encoding CTCFL, PRMT7, and the de novo DNA methyltransferases, Dnmt3a, -b and -L, in Xenopus oocytes. These results suggest that CTCFL and PRMT7 may play a role in male germline imprinted gene methylation.

The testes-specific factor CTCFL (also called BORIS) can cooperate with the arginine histone methyltransferase (PRMT7) to induce de novo DNA methylation at an imprinting control region.

In silico approaches reveal the potential for DNA sequence-dependent histone octamer affinity to influence chromatin structure in vivo

Nucleosome positioning signals embedded within the DNA sequence have the potential to influence the detailed structure of the higher-order chromatin fibre. In two previous studies of long stretches of DNA, encompassing the chicken beta-globin and ovine beta-lactoglobulin genes, respectively, we mapped the relative affinity of every site for the core histone octamer. In both cases a periodic arrangement of the in vitro positioning sites suggests that they might influence the folding of a nucleosome chain into higher-order structure; this hypothesis was borne out in the case of the beta-lactoglobulin gene, where the distribution of the in vitro positioning sites is related to the positions nucleosomes actually occupy in sheep liver cells. Here, we have exploited the in vitro nucleosome positioning datasets to simulate nucleosomal organisation using in silico approaches. We use the high-resolution, quantitative positioning maps to define a one-dimensional positioning energy lattice, which can be populated with a defined number of nucleosomes. Monte Carlo techniques are employed to simulate the behaviour of the model at equilibrium to produce a set of configurations, which provide a probability-based occupancy map. Employing a variety of techniques we show that the occupancy maps are a sensitive function of the histone octamer density (nucleosome repeat length) and find that a minimal change in this property can produce dramatic localised changes in structure. Although simulations generally give rise to regular periodic nucleosomal arrangements, they often show octamer density-dependent discontinuities, which tend to co-localise with sequences that adopt distinctive chromatin structure in vivo. Furthermore, the overall organisation of simulated chromatin structures are more closely related to the situation in vivo than is the original in vitro positioning data, particularly at a nucleosome density corresponding to the in vivo state. Although our model is simplified, we argue that it provides a unique insight into the influence that DNA sequence can have in determining chromatin structure and could serve as a useful basis for the incorporation of other parameters.

NXSensor web tool for evaluating DNA for nucleosome exclusion sequences and accessibility to binding factors

Nucleosomes, a basic structural unit of eukaryotic chromatin, play a significant role in regulating gene expression. We have developed a web tool based on DNA sequences known from empirical and theoretical studies to influence DNA bending and flexibility, and to exclude nucleosomes. NXSensor (available at ) finds nucleosome exclusion sequences, evaluates their length and spacing, and computes an ‘accessibility score’ giving the proportion of base pairs likely to be nucleosome-free. Application of NXSensor to the promoter regions of housekeeping (HK) genes and those of tissue-specific (TS) genes revealed a significant difference between the two classes of gene, the former being significantly more open, on average, particularly near transcription start sites (TSSs). NXSensor should be a useful tool in assessing the likelihood of nucleosome formation in regions involved in gene regulation and other aspects of chromatin function.

Biochemical screening of stable dinucleosomes using DNA fragments from a dinucleosome DNA library

The dinucleosome is an informative unit for analysis of the higher-order chromatin structure. DNA fragments forming stable dinucleosomes were screened from a dinucleosome DNA library after the reconstitution of nucleosomes in vitro and digestion with micrococcal nuclease. Reconstituted dinucleosomes showed a diversity of sensitivity to micrococcal nuclease, suggesting that the biochemical stability of a dinucleosome depends, in part, on the DNA fragments. The DNA fragments after the screening were classified into three groups represented by clones bf10, af14 and af32 according to the sensitivity to micrococcal nuclease. Mapping of the nucleosome boundaries by Southern blotting of the DNA after restriction digestion and by primer extension analysis showed that each nucleosome position of clone af32 was fixed. Analysis of reconstituted dinucleosomes using mutant DNA fragments of clone af32 revealed a unique property characteristic of a key nucleosome, given that the replacement of a DNA fragment corresponding to the right nucleosome position resulted in marked sensitivity to micrococcal nuclease, whereas the replacement of the other nucleosome fragment had almost no effect on sensitivity as compared to the original af32 construct. The mutant construct in which the right nucleosome was removed showed multiple nucleosome phases, suggesting that the right nucleosome stabilized first each mononucleosome and then the dinucleosome. An oligonucleotide bending assay revealed that the DNA fragment in the right nucleosome included curved DNA, suggesting that the positioning activity of the nucleosome was attributed to its DNA structure. These results suggest that information for forming stable dinucleosome is embedded in the genomic DNA and that a further characterization of the key nucleosome is useful for understanding the building up of the chromatin structure.

Microdeletion of target sites for insulator protein CTCF in a chromosome 11p15 imprinting center in Beckwith-Wiedemann syndrome and Wilms' tumor

We have analyzed several cases of Beckwith-Wiedemann syndrome (BWS) with Wilms' tumor in a familial setting, which give insight into the complex controls of imprinting and gene expression in the chromosome 11p15 region. We describe a 2.2-kbp microdeletion in the H19/insulin-like growth factor 2 (IGF2)-imprinting center eliminating three target sites of the chromatin insulator protein CTCF that we believe here is necessary, but not sufficient, to cause BWS and Wilms' tumor. Maternal inheritance of the deletion is associated with IGF2 loss of imprinting and up-regulation of IGF2 mRNA. However, in at least one affected family member a second genetic lesion (a duplication of maternal 11p15) was identified and accompanied by a further increase in IGF2 mRNA levels 35-fold higher than control values. Our results suggest that the combined effects of the H19/IGF2-imprinting center microdeletion and 11p15 chromosome duplication were necessary for manifestation of BWS.

A determining influence for CpG dinucleotides on nucleosome positioning in vitro

DNA sequence information that directs the translational positioning of nucleosomes can be attenuated by cytosine methylation when a short run of CpG dinucleotides is located close to the dyad axis of the nucleosome. Here, we show that point mutations introduced to re-pattern methylation at the (CpG)3 element in the chicken betaA-globin promoter sequence themselves strongly influenced nucleosome formation in reconstituted chromatin. The disruptive effect of cytosine methylation on nucleosome formation was found to be determined by the sequence context of CpG dinucleotides, not just their location in the positioning sequence. Additional mutations indicated that methylation can also promote the occupation of certain nucleosome positions. DNase I analysis demonstrated that these genetic and epigenetic modifications altered the structural characteristics of the (CpG)3 element. Our findings support a proposal that the intrinsic structural properties of the DNA at the -1.5 site, as occupied by (CpG)3 in the nucleosome studied, can be decisive for nucleosome formation and stability, and that changes in anisotropic DNA bending or flexibility at this site explain why nucleosome positioning can be exquisitely sensitive to genetic and epigenetic modification of the DNA sequence.

Suv39h histone methyltransferases interact with Smads and cooperate in BMP-induced repression

Smad proteins transduce signals from transforming growth factor-beta (TGF-beta) superfamily ligands to regulate the expression of target genes. In order to identify novel partners of Smad proteins in transcriptional regulation, we performed a two-hybrid screen using Smad5, a protein that is activated predominantly by bone morphogenetic protein (BMP) signaling. We identified an interaction between Smad5 and suppressor of variegation 3-9 homolog 2 (Suv39h2), a chromatin modifier enzyme. Suv39h proteins are histone methyltransferases that methylate histone H3 on lysine 9, resulting in transcriptional repression or silencing of target genes. Biochemical studies in mammalian cells demonstrated that Smad5 binds to both known mammalian isoforms of Suv39h proteins, and that Smad proteins activated by the TGF-beta signaling pathway, Smad2 and Smad3, do not bind with significant affinity. Functional studies using the muscle creatine kinase (MCK) promoter, which is suppressed by BMP signaling, demonstrate that Suv39h proteins and Smads cooperate to repress promoter activity. These data suggest a model where association of Smad proteins with Suv39h methyltransferases can repress or silence genes involved in developmental processes, and argues that inefficient gene repression may result in the alteration of the differentiated phenotype. Thus, examination of the Smad-Suv interaction may provide insight into the mechanism of phenotypic determination mediated by BMP signaling.

Nucleosome positioning signals and potential H-DNA within the DNA sequence of the imprinting control region of the mouse Igf2r gene

The imprinting control region within the second intron of the mouse Igf2r gene contains a CpG island comprising direct repeats, an imprinting box and the Air antisense promoter which is blocked by the methylation imprint on the active maternal allele. We have investigated the structural features of this DNA, including a mapping of all nucleosome positioning signals within the nucleotide sequence. A discrete series of strong positioning signals distinguished the direct repeat region from the much more diverse positioning capacity of the sequence encompassing the known regulatory elements. At only a few locations did CpG methylation modulate the use of this positioning information. Direct effects upon histone-DNA interactions are therefore unlikely to contribute significantly to the means by which the imprint may establish allele-specific chromatin architecture and determine Air expression. A strand-specific obstruction to DNA polymerase was observed between the repeat and regulatory regions. The same region adopts triple-stranded H-DNA structures in supercoiled DNA, according to pH and divalent cation exposure. Methylation did not modulate the occurrence or form of this structure under the conditions tested. This finding nevertheless adds to the repertoire of potential H-DNA structures found in the vicinity of regulatory sequences-here, in an imprinting context.