DNA sequence preferences of transcriptional activators correlate more strongly than repressors with nucleosomes

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

In vitro versus in vivo compositional landscapes of histone sequence preferences in eucaryotic genomes

Motivation

Although the nucleosome occupancy along a genome can be in part predicted by in vitro experiments, it has been recently observed that the chromatin organization presents important differences in vitro with respect to in vivo. Such differences mainly regard the hierarchical and regular structures of the nucleosome fiber, whose existence has long been assumed, and in part also observed in vitro, but that does not apparently occur in vivo. It is also well known that the DNA sequence has a role in determining the nucleosome occupancy. Therefore, an important issue is to understand if, and to what extent, the structural differences in the chromatin organization between in vitro and in vivo have a counterpart in terms of the underlying genomic sequences.

Results

We present the first quantitative comparison between the in vitro and in vivo nucleosome maps of two model organisms (S. cerevisiae and C. elegans). The comparison is based on the construction of weighted k-mer dictionaries. Our findings show that there is a good level of sequence conservation between in vitro and in vivo in both the two organisms, in contrast to the abovementioned important differences in chromatin structural organization. Moreover, our results provide evidence that the two organisms predispose themselves differently, in terms of sequence composition and both in vitro and in vivo, for the nucleosome occupancy. This leads to the conclusion that, although the notion of a genome encoding for its own nucleosome occupancy is general, the intrinsic histone k-mer sequence preferences tend to be species-specific.

Availability and implementation

The files containing the dictionaries and the main results of the analysis are available at http://math.unipa.it/rombo/material.

Supplementary information

Supplementary data are available at Bioinformatics online.

Application of MNase-Seq in the Global Mapping of Nucleosome Positioning in Plants

The precise positioning of nucleosomes along the underlying DNA is critical for a variety of biological processes, especially in regulating transcription. The interplay between nucleosomes and transcription factors for accessing the underlying DNA sequences is one of the key determinants that affect transcriptional regulation. Moreover, nucleosomes with various packing statuses confer distinct functions in regulating gene expressions in response to various internal or external signals. Therefore, global mapping of nucleosome positions is one informative way to elucidate the relationship between patterns of nucleosome positioning/occupancy and transcriptional regulations. MNase digestion coupled with high-throughput sequencing (MNase-seq) has been utilized widely for global mapping of nucleosome positioning in eukaryotes that have a sequenced genome. We have developed a robust MNase-seq procedure in plants. It mainly includes plant nuclei isolation, treatment of purified nuclei with MNase, gel recovery of MNase-trimmed mononucleosomal DNA with an approximate size of 150 bp, MNase-seq library preparation followed by Illumina sequencing, and data analysis. MNase-seq has already been successfully applied to identify genome-wide nucleosome positioning in model plants, rice, and Arabidopsis thaliana.

Gene expression and promoter characterization of heat-shock protein 90B gene (HSP90B) in the model unicellular green alga Chlamydomonas reinhardtii

Molecular chaperones or heat shock proteins are a large protein family with important functions in every cellular organism. Among all types of the heat shock proteins, information on the ER-localized HSP90 protein (HSP90B) and its encoding gene is relatively scarce in the literature, especially in photosynthetic organisms. In this study, expression profiles as well as promoter sequence of the HSP90B gene were investigated in the model green alga Chlamydomonas reinhardtii. We have found that HSP90B is strongly induced by heat and ER stresses, while other short-term exposure to abiotic stresses, such as salinity, dark-to-light transition or light stress does not appear to affect the expression. Promoter truncation analysis as well as chromatin immunoprecipitation using the antibodies recognizing histone H3 and acetylated histone H3, revealed a putative core constitutive promoter sequence between -1 to -253 bp from the transcription start site. Our results also suggested that the nucleotides upstream of the core promoter may contain repressive elements such as putative repressor binding site(s).

RNA Nuclear Export: From Neurological Disorders to Cancer

The presence of a nuclear envelope, also known as nuclear membrane, defines the structural framework of all eukaryotic cells by separating the nucleus, which contains the genetic material, from the cytoplasm where the synthesis of proteins takes place. Translation of proteins in Eukaryotes is thus dependent on the active transport of DNA-encoded RNA molecules through pores embedded within the nuclear membrane. Several mechanisms are involved in this process generally referred to as RNA nuclear export or nucleocytoplasmic transport of RNA. The regulated expression of genes requires the nuclear export of protein-coding messenger RNA molecules (mRNAs) as well as non-coding RNAs (ncRNAs) together with proteins and pre-assembled ribosomal subunits. The nuclear export of mRNAs is intrinsically linked to the co-transcriptional processing of nascent transcripts synthesized by the RNA polymerase II. This functional coupling is essential for the survival of cells allowing for timely nuclear export of fully processed transcripts, which could otherwise cause the translation of abnormal proteins such as the polymeric repeat proteins produced in some neurodegenerative diseases. Alterations of the mRNA nuclear export pathways can also lead to genome instability and to various forms of cancer. This chapter will describe the molecular mechanisms driving the nuclear export of RNAs with a particular emphasis on mRNAs. It will also review their known alterations in neurological disorders and cancer, and the recent opportunities they offer for the potential development of novel therapeutic strategies.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Transcriptional Regulation of the Ambient Temperature Response by H2A.Z Nucleosomes and HSF1 Transcription Factors in Arabidopsis

Temperature influences the distribution, range, and phenology of plants. The key transcriptional activators of heat shock response in eukaryotes, the heat shock factors (HSFs), have undergone large-scale gene amplification in plants. While HSFs are central in heat stress responses, their role in the response to ambient temperature changes is less well understood. We show here that the warm ambient temperature transcriptome is dependent upon the HSFA1 clade of Arabidopsis HSFs, which cause a rapid and dynamic eviction of H2A.Z nucleosomes at target genes. A transcriptional cascade results in the activation of multiple downstream stress-responsive transcription factors, triggering large-scale changes to the transcriptome in response to elevated temperature. H2A.Z nucleosomes are enriched at temperature-responsive genes at non-inducible temperature, and thus likely confer inducibility of gene expression and higher responsive dynamics. We propose that the antagonistic effects of H2A.Z and HSF1 provide a mechanism to activate gene expression rapidly and precisely in response to temperature, while preventing leaky transcription in the absence of an activation signal.

Negative selection maintains transcription factor binding motifs in human cancer

Background

Somatic mutations in cancer cells affect various genomic elements disrupting important cell functions. In particular, mutations in DNA binding sites recognized by transcription factors can alter regulator binding affinities and, consequently, expression of target genes. A number of promoter mutations have been linked with an increased risk of cancer. Cancer somatic mutations in binding sites of selected transcription factors have been found under positive selection. However, action and significance of negative selection in non-coding regions remain controversial.

Results

Here we present analysis of transcription factor binding motifs co-localized with non-coding variants. To avoid statistical bias we account for mutation signatures of different cancer types. For many transcription factors, including multiple members of FOX, HOX, and NR families, we show that human cancers accumulate fewer mutations than expected by chance that increase or decrease affinity of predicted binding sites. Such stability of binding motifs is even more exhibited in DNase accessible regions.

Conclusions

Our data demonstrate negative selection against binding sites alterations and suggest that such selection pressure protects cancer cells from rewiring of regulatory circuits. Further analysis of transcription factors with conserved binding motifs can reveal cell regulatory pathways crucial for the survivability of various human cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2728-9) contains supplementary material, which is available to authorized users.

How motif environment influences transcription factor search dynamics: Finding a needle in a haystack

Transcription factors (TFs) have to find their binding sites, which are distributed throughout the genome. Facilitated diffusion is currently the most widely accepted model for this search process. Based on this model the TF alternates between one-dimensional sliding along the DNA, and three-dimensional bulk diffusion. In this view, the non-specific associations between the proteins and the DNA play a major role in the search dynamics. However, little is known about how the DNA properties around the motif contribute to the search. Accumulating evidence showing that TF binding sites are embedded within a unique environment, specific to each TF, leads to the hypothesis that the search process is facilitated by favorable DNA features that help to improve the search efficiency. Here, we review the field and present the hypothesis that TF-DNA recognition is dictated not only by the motif, but is also influenced by the environment in which the motif resides.

Promoter architectures and developmental gene regulation

Core promoters are minimal regions sufficient to direct accurate initiation of transcription and are crucial for regulation of gene expression. They are highly diverse in terms of associated core promoter motifs, underlying sequence composition and patterns of transcription initiation. Distinctive features of promoters are also seen at the chromatin level, including nucleosome positioning patterns and presence of specific histone modifications. Recent advances in identifying and characterizing promoters using next-generation sequencing-based technologies have provided the basis for their classification into functional groups and have shed light on their modes of regulation, with important implications for transcriptional regulation in development. This review discusses the methodology and the results of genome-wide studies that provided insight into the diversity of RNA polymerase II promoter architectures in vertebrates and other Metazoa, and the association of these architectures with distinct modes of regulation in embryonic development and differentiation.

"Hit-and-run": Transcription factors get caught in the act

A key challenge for understanding transcriptional regulation is being able to measure transcription factor (TF)-DNA binding events with sufficient spatial and temporal resolution; that is, when and where TFs occupy their cognate sites. A recent study by Para et al. has highlighted the dynamics underlying the activation of gene expression by a master regulator TF. This study provides concrete evidence for a long-standing hypothesis in biology, the "hit-and-run" mechanism, which was first proposed decades ago. That is, gene expression is dynamically controlled by a TF that transiently binds and activates a target gene, which might stay in a transcriptionally active state after the initial binding event has ended. Importantly, the experimental procedure introduced, TARGET, provides a useful way for identifying multiple target genes transiently bound by their regulators, which can be used in conjunction with other well-established methods to improve our understanding of transcriptional regulatory dynamics.

Differential Nucleosome Occupancies across Oct4-Sox2 Binding Sites in Murine Embryonic Stem Cells

The binding sequence for any transcription factor can be found millions of times within a genome, yet only a small fraction of these sequences encode functional transcription factor binding sites. One of the reasons for this dichotomy is that many other factors, such as nucleosomes, compete for binding. To study how the competition between nucleosomes and transcription factors helps determine a functional transcription factor site from a predicted transcription factor site, we compared experimentally-generated in vitro nucleosome occupancy with in vivo nucleosome occupancy and transcription factor binding in murine embryonic stem cells. Using a solution hybridization enrichment technique, we generated a high-resolution nucleosome map from targeted regions of the genome containing predicted sites and functional sites of Oct4/Sox2 regulation. We found that at Pax6 and Nes, which are bivalently poised in stem cells, functional Oct4 and Sox2 sites show high amounts of in vivo nucleosome displacement compared to in vitro. Oct4 and Sox2, which are active, show no significant displacement of in vivo nucleosomes at functional sites, similar to nonfunctional Oct4/Sox2 binding. This study highlights a complex interplay between Oct4 and Sox2 transcription factors and nucleosomes among different target genes, which may result in distinct patterns of stem cell gene regulation.

Promoter-proximal transcription factor binding is transcriptionally active when coupled with nucleosome repositioning in immediate vicinity

Previous studies have analyzed patterns of transcription, transcription factor (TF) binding or mapped nucleosome occupancy across the genome. These suggest that the three aspects are genetically connected but the cause and effect relationships are still unknown. For example, physiologic TF binding studies involve many TFs, consequently, it is difficult to assign nucleosome reorganization to the binding site occupancy of any particular TF. Therefore, several aspects remain unclear: does TF binding influence nucleosome (re)organizations locally or impact the chromatin landscape at a more global level; are all or only a fraction of TF binding a result of reorganization in nucleosome occupancy and do all TF binding and associated changes in nucleosome occupancy result in altered gene expression? With these in mind, following characterization of two states (before and after induction of a single TF of choice) we determined: (i) genomic binding sites of the TF, (ii) promoter nucleosome occupancy and (iii) transcriptome profiles. Results demonstrated that promoter-proximal TF binding influenced expression of the target gene when it was coupled to nucleosome repositioning at or close to its binding site in most cases. In contrast, only in few cases change in target gene expression was found when TF binding occurred without local nucleosome reorganization.

Nucleosome organization in the vicinity of transcription factor binding sites in the human genome

Background

The binding of transcription factors (TFs) to specific DNA sequences is an initial and crucial step of transcription. In eukaryotes, this process is highly dependent on the local chromatin state, which can be modified by recruiting chromatin remodelers. However, previous studies have focused mainly on nucleosome occupancy around the TF binding sites (TFBSs) of a few specific TFs. Here, we investigated the nucleosome occupancy profiles around computationally inferred binding sites, based on 519 TF binding motifs, in human GM12878 and K562 cells.

Results

Although high nucleosome occupancy is intrinsically encoded at TFBSs in vitro, nucleosomes are generally depleted at TFBSs in vivo, and approximately a quarter of TFBSs showed well-positioned in vivo nucleosomes on both sides. RNA polymerase near the transcription start site (TSS) has a large effect on the nucleosome occupancy distribution around the binding sites located within one kilobase to the nearest TSS; fuzzier nucleosome positioning was thus observed around these sites. In addition, in contrast to yeast, repressors, rather than activators, were more likely to bind to nucleosomal DNA in the human cells, and nucleosomes around repressor sites were better positioned in vivo. Genes with repressor sites exhibiting well-positioned nucleosomes on both sides, and genes with activator sites occupied by nucleosomes had significantly lower expression, suggesting that actions of activators and repressors are associated with the nucleosome occupancy around their binding sites. It was also interesting to note that most of the binding sites, which were not in the DNase I-hypersensitive regions, were cell-type specific, and higher in vivo nucleosome occupancy were observed at these binding sites.

Conclusions

This study demonstrated that RNA polymerase and the functions of bound TFs affected the local nucleosome occupancy around TFBSs, and nucleosome occupancy patterns around TFBSs were associated with the expression levels of target genes.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-493) contains supplementary material, which is available to authorized users.

Coregulation of transcription factor binding and nucleosome occupancy through DNA features of mammalian enhancers

Transcription factors (TFs) preferentially bind sites contained in regions of computationally predicted high nucleosomal occupancy, suggesting that nucleosomes are gatekeepers of TF binding sites. However, because of their complexity mammalian genomes contain millions of randomly occurring, unbound TF consensus binding sites. We hypothesized that the information controlling nucleosome assembly may coincide with the information that enables TFs to bind cis-regulatory elements while ignoring randomly occurring sites. Hence, nucleosomes would selectively mask genomic sites that can be contacted by TFs and thus be potentially functional. The hematopoietic pioneer TF Pu.1 maintained nucleosome depletion at macrophage-specific enhancers that displayed a broad range of nucleosome occupancy in other cell types and in reconstituted chromatin. We identified a minimal set of DNA sequence and shape features that accurately predicted both Pu.1 binding and nucleosome occupancy genome-wide. These data reveal a basic organizational principle of mammalian cis-regulatory elements whereby TF recruitment and nucleosome deposition are controlled by overlapping DNA sequence features.

Positive and negative design for nonconsensus protein-DNA binding affinity in the vicinity of functional binding sites

Recent experiments provide an unprecedented view of protein-DNA binding in yeast and human genomes at single-nucleotide resolution. These measurements, performed over large cell populations, show quite generally that sequence-specific transcription regulators with well-defined protein-DNA consensus motifs bind only a fraction among all consensus motifs present in the genome. Alternatively, proteins in vivo often bind DNA regions lacking known consensus sequences. The rules determining whether a consensus motif is functional remain incompletely understood. Here we predict that genomic background surrounding specific protein-DNA binding motifs statistically modulates the binding of sequence-specific transcription regulators to these motifs. In particular, we show that nonconsensus protein-DNA binding in yeast is statistically enhanced, on average, around functional Reb1 motifs that are bound as compared to nonfunctional Reb1 motifs that are unbound. The landscape of nonconsensus protein-DNA binding around functional CTCF motifs in human demonstrates a more complex behavior. In particular, human genomic regions characterized by the highest CTCF occupancy, show statistically reduced level of nonconsensus protein-DNA binding. Our findings suggest that nonconsensus protein-DNA binding is fine-tuned around functional binding sites using a variety of design strategies.

Clustered regulatory elements at nucleosome-depleted regions punctuate a constant nucleosomal landscape in Schizosaccharomyces pombe

BACKGROUND

Nucleosomes facilitate the packaging of the eukaryotic genome and modulate the access of regulators to DNA. A detailed description of the nucleosomal organization under different transcriptional programmes is essential to understand their contribution to genomic regulation.

RESULTS

To visualize the dynamics of individual nucleosomes under different transcriptional programmes we have generated high-resolution nucleosomal maps in Schizosaccharomyces pombe. We show that 98.5% of the genome remains almost invariable during mitosis and meiosis while remodelling is limited to approximately 1100 nucleosomes in the promoters of a subset of meiotic genes. These inducible nucleosome-depleted regions (NDR) and also those constitutively present in the genome overlap precisely with clusters of binding sites for transcription factors (TF) specific for meiosis and for different functional classes of genes, respectively. Deletion of two TFs affects only a small fraction of all the NDRs to which they bind in vivo, indicating that TFs collectively contribute to NDR maintenance.

CONCLUSIONS

Our results show that the nucleosomal profile in S. pombe is largely maintained under different physiological conditions and patterns of gene expression. This relatively constant landscape favours the concentration of regulators in constitutive and inducible NDRs. The combinatorial analysis of binding motifs in this discrete fraction of the genome will facilitate the definition of the transcriptional regulatory networks.

Reprogramming cells for brain repair

At present there are no clinical therapies that can repair traumatic brain injury, spinal cord injury or degenerative brain disease. While redundancy and rewiring of surviving circuits can recover some lost function, the brain and spinal column lack sufficient endogenous stem cells to replace lost neurons or their supporting glia. In contrast, pre-clinical studies have demonstrated that exogenous transplants can have remarkable efficacy for brain repair in animal models. Mesenchymal stromal cells (MSCs) can provide paracrine factors that repair damage caused by ischemic injury, and oligodendrocyte progenitor cell (OPC) grafts give dramatic functional recovery from spinal cord injury. These studies have progressed to clinical trials, including human embryonic stem cell (hESC)-derived OPCs for spinal cord repair. However, ESC-derived allografts are less than optimal, and we need to identify a more appropriate donor graft population. The cell reprogramming field has developed the ability to trans-differentiate somatic cells into distinct cell types, a technology that has the potential to generate autologous neurons and glia which address the histocompatibility concerns of allografts and the tumorigenicity concerns of ESC-derived grafts. Further clarifying how cell reprogramming works may lead to more efficient direct reprogram approaches, and possibly in vivo reprogramming, in order to promote brain and spinal cord repair.

Contribution of nucleosome binding preferences and co-occurring DNA sequences to transcription factor binding

BACKGROUND

Chromatin plays a critical role in regulating transcription factors (TFs) binding to their canonical transcription factor binding sites (TFBS). Recent studies in vertebrates show that many TFs preferentially bind to genomic regions that are well bound by nucleosomes in vitro. Co-occurring secondary motifs sometimes correlated with functional TFBS.

RESULTS

We used a logistic regression to evaluate how well the propensity for nucleosome binding and co-occurrence of a secondary motif identify which canonical motifs are bound in vivo. We used ChIP-seq data for three transcription factors binding to their canonical motifs: c-Jun binding the AP-1 motif (TGA(C)/(G)TCA), GR (glucocorticoid receptor) binding the GR motif (G-ACA---(T)/(C)GT-C), and Hoxa2 (homeobox a2) binding the Pbx (Pre-B-cell leukemia homeobox) motif (TGATTGAT). For all canonical TFBS in the mouse genome, we calculated intrinsic nucleosome occupancy scores (INOS) for its surrounding 150-bps DNA and examined the relationship with in vivo TF binding. In mouse mammary 3134 cells, c-Jun and GR proteins preferentially bound regions calculated to be well-bound by nucleosomes in vitro with the canonical AP-1 and GR motifs themselves contributing to the high INOS. Functional GR motifs are enriched for AP-1 motifs if they are within a nucleosome-sized 150-bps region. GR and Hoxa2 also bind motifs with low INOS, perhaps indicating a different mechanism of action.

CONCLUSION

Our analysis quantified the contribution of INOS and co-occurring sequence to the identification of functional canonical motifs in the genome. This analysis revealed an inherent competition between some TFs and nucleosomes for binding canonical TFBS. GR and c-Jun cooperate if they are within 150-bps. Binding of Hoxa2 and a fraction of GR to motifs with low INOS values suggesting they are not in competition with nucleosomes and may function using different mechanisms.

Nucleosome assembly factors CAF-1 and HIR modulate epigenetic switching frequencies in an H3K56 acetylation-associated manner in Candida albicans

CAF-1 and HIR are highly conserved histone chaperone protein complexes that function in the assembly of nucleosomes onto chromatin. CAF-1 is characterized as having replication-coupled nucleosome activity, whereas the HIR complex can assemble nucleosomes independent of replication. Histone H3K56 acetylation, controlled by the acetyltransferase Rtt109 and deacetylase Hst3, also plays a significant role in nucleosome assembly. In this study, we generated a set of deletion mutants to genetically characterize pathway-specific and overlapping functions of CAF-1 and HIR in C. albicans. Their roles in epigenetic maintenance of cell type were examined by using the white-opaque switching system in C. albicans. We show that CAF-1 and HIR play conserved roles in UV radiation recovery, repression of histone gene expression, correct chromosome segregation, and stress responses. Unique to C. albicans, the cac2Δ/Δ mutant shows increased sensitivity to the Hst3 inhibitor nicotinamide, while the rtt109Δ/Δ cac2Δ/Δ and hir1Δ/Δ cac2Δ/Δ mutants are resistant to nicotinamide. CAF-1 plays a major role in maintaining cell types, as the cac2Δ/Δ mutant exhibited increased switching frequencies in both directions and switched at a high frequency to opaque in response to nicotinamide. Like the rtt109Δ/Δ mutant, the hir1Δ/Δ cac2Δ/Δ double mutant is defective in maintaining the opaque cell fate and blocks nicotinamide-induced opaque formation, and the defects are suppressed by ectopic expression of the master white-opaque regulator Wor1. Our data suggest an overlapping function of CAF-1 and HIR in epigenetic regulation of cell fate determination in an H3K56 acetylation-associated manner.