A genome-wide analysis in Saccharomyces cerevisiae demonstrates the influence of chromatin modifiers on transcription

Transcriptional repression by the histone tails in budding yeast is mediated by Rpd3, Tup1-Ssn6, and Bur6/NC2

Chromatin-mediated transcriptional regulation is modulated by post-translational modifications of the core histones, particularly the H3 and H4 unstructured amino termini, or "tails". In budding yeast, the H3 and H4 tails can be deacetylated by Rpd3 to repress specific target genes, and hypoacetylated histones can facilitate recruitment of the Tup1-Ssn6 complex to effect gene repression. However, the extent to which these mechanisms are used to effect repression by the histone tails, and whether other factors similarly collaborate with the tails to facilitate gene repression, has not been determined. Here, a chromatin modifier compendium of 170 gene expression profiles from yeast strains mutated for chromatin-related genes was used to query the effect of the corresponding mutations on gene cohorts repressed by the histone H3 and H4 tails and/or by Rpd3. The resulting analysis reveals that repression of nearly all of the genes repressed by the histone tails requires Rpd3 and/or the Tup1-Ssn6 complex. Repression by Rpd3 occurs via the Rpd3-L complex, and TFIID-dominated genes are underrepresented among genes repressed by mutations or deletions of the H3 or H4 tails, in accord with previous work. In addition, Bur6, the yeast homolog of human NC2α, is required for repression at ∼50% of genes repressed by the H3 or H4 tail. These results illuminate genome-wide repression mechanisms utilized by the histone tails in yeast and raise new questions regarding the role of Bur6 in histone tail-mediated repression and whether parallels exist in metazoan cells.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

CRISPR/dCas9-mediated epigenetic modification reveals differential regulation of histone acetylation on Aspergillus niger secondary metabolite

Epigenetic studies on secondary metabolites (SMs) mainly relied so far on non-selective epigenetic factors deletion or feeding epigenetic inhibitors in Aspergillus niger. Although technologies developed for epigenome editing at specific loci now enable the direct study of the functional relevance of precise gene regulation and epigenetic modification, relevant assays are limited in filamentous fungi. Herein, we show that CRISPR/dCas9-mediated histone epigenetic modification systems efficiently reprogramed the expression of target genes in A. niger. First, we constructed a p300-dCas9 system and demonstrated the activation of a EGFP fluorescent reporter. Second, by precisely locating histone acetylase p300 on ATG adjacent region of secondary metabolic gene breF, the transcription of breF was activated. Third, p300-dCas9 was guided to the native polyketide synthase (PKS) gene fuml, which increased production of the compound fumonisin B2 detected by HPLC and LC-MS. Then, endogenous histone acetylase GcnE-dCas9 and histone deacetylases HosA-dCas9 and RpdA-dCas9 repressed the transcription of breF. Finally, by targeting HosA-dCa9 fusion to pigment gene fwnA, we confirmed that histone deacetylase HosA activated the expression of fwnA, accelerated the synthesis of melanin. Targeted epigenome editing is a promising technology and this study is the first time to apply the epigenetic CRISPR/dCas9 system on regulating the expression of the secondary metabolic genes in A. niger.

Optogenetic Control Reveals Differential Promoter Interpretation of Transcription Factor Nuclear Translocation Dynamics

Gene expression is thought to be affected not only by the concentration of transcription factors (TFs) but also the dynamics of their nuclear translocation. Testing this hypothesis requires direct control of TF dynamics. Here, we engineer CLASP, an optogenetic tool for rapid and tunable translocation of a TF of interest. Using CLASP fused to Crz1, we observe that, for the same integrated concentration of nuclear TF over time, changing input dynamics changes target gene expression: pulsatile inputs yield higher expression than continuous inputs, or vice versa, depending on the target gene. Computational modeling reveals that a dose-response saturating at low TF input can yield higher gene expression for pulsatile versus continuous input, and that multi-state promoter activation can yield the opposite behavior. Our integrated tool development and modeling approach characterize promoter responses to Crz1 nuclear translocation dynamics, extracting quantitative features that may help explain the differential expression of target genes.

CLASP is a modular optogenetic strategy to control the nuclear localization of transcription factors (TFs) and elicit gene expression from their cognate promoters. CLASP control of Crz1 nuclear localization, coupled with computational modeling, revealed how promoters can differentially decode dynamic transcription factor signals. The integrated strategy of CLASP development and modeling presents a generalized approach to causally investigate the transcriptional consequences of dynamic TF nuclear shuttling.

The cis-regulatory codes of response to combined heat and drought stress in Arabidopsis thaliana

Plants respond to their environment by dynamically modulating gene expression. A powerful approach for understanding how these responses are regulated is to integrate information about cis-regulatory elements (CREs) into models called cis-regulatory codes. Transcriptional response to combined stress is typically not the sum of the responses to the individual stresses. However, cis-regulatory codes underlying combined stress response have not been established. Here we modeled transcriptional response to single and combined heat and drought stress in Arabidopsis thaliana. We grouped genes by their pattern of response (independent, antagonistic and synergistic) and trained machine learning models to predict their response using putative CREs (pCREs) as features (median F-measure = 0.64). We then developed a deep learning approach to integrate additional omics information (sequence conservation, chromatin accessibility and histone modification) into our models, improving performance by 6.2%. While pCREs important for predicting independent and antagonistic responses tended to resemble binding motifs of transcription factors associated with heat and/or drought stress, important synergistic pCREs resembled binding motifs of transcription factors not known to be associated with stress. These findings demonstrate how in silico approaches can improve our understanding of the complex codes regulating response to combined stress and help us identify prime targets for future characterization.

Complex genetic and epigenetic regulation deviates gene expression from a unifying global transcriptional program

Environmental or genetic perturbations lead to gene expression changes. While most analyses of these changes emphasize the presence of qualitative differences on just a few genes, we now know that changes are widespread. This large-scale variation has been linked to the exclusive influence of a global transcriptional program determined by the new physiological state of the cell. However, given the sophistication of eukaryotic regulation, we expect to have a complex architecture of specific control affecting this program. Here, we examine this architecture. Using data of Saccharomyces cerevisiae expression in different nutrient conditions, we first propose a five-sector genome partition, which integrates earlier models of resource allocation, as a framework to examine the deviations from the global control. In this scheme, we recognize invariant genes, whose regulation is dominated by physiology, specific genes, which substantially depart from it, and two additional classes that contain the frequently assumed growth-dependent genes. Whereas the invariant class shows a considerable absence of specific regulation, the rest is enriched by regulation at the level of transcription factors (TFs) and epigenetic modulators. We nevertheless find markedly different strategies in how these classes deviate. On the one hand, there are TFs that act in a unique way between partition constituents, and on the other, the action of chromatin modifiers is significantly diverse. The balance between regulatory strategies ultimately modulates the action of the general transcription machinery and therefore limits the possibility of establishing a unifying program of expression change at a genomic scale.

How can we understand expression changes observed as a result of environmental or genetic perturbations? This issue has been conventionally answered by examining small groups of genes whose expression becomes qualitatively altered after these perturbations. But this approach is too simplistic, as we now know that extensive variation is typically observed. To explain this variation, recent works proposed a model in which genome-wide changes are explained by the action of a general program of transcription. Our manuscript reasons that given the complexity of eukaryotic transcriptional control, a unifying program of regulation cannot be achievable. Instead, we propose within an integrated framework of resource allocation that a rich structure of deviations from it exists and that by characterizing these deviations we can fully appreciate large-scale expression change.

The Histone Deacetylases HosA and HdaA Affect the Phenotype and Transcriptomic and Metabolic Profiles of Aspergillus niger

Histone acetylation is an important modification for the regulation of chromatin accessibility and is controlled by two kinds of histone-modifying enzymes: histone acetyltransferases (HATs) and histone deacetylases (HDACs). In filamentous fungi, there is increasing evidence that HATs and HDACs are critical factors related to mycelial growth, stress response, pathogenicity and production of secondary metabolites (SMs). In this study, seven A. niger histone deacetylase-deficient strains were constructed to investigate their effects on the strain growth phenotype as well as the transcriptomic and metabolic profiles of secondary metabolic pathways. Phenotypic analysis showed that deletion of hosA in A. niger FGSC A1279 leads to a significant reduction in growth, pigment production, sporulation and stress resistance, and deletion of hdaA leads to an increase in pigment production in liquid CD medium. According to the metabolomic analysis, the production of the well-known secondary metabolite fumonisin was reduced in both the hosA and hdaA mutants, and the production of kojic acid was reduced in the hdaA mutant and slightly increased in the hosA mutant. Results suggested that the histone deacetylases HosA and HdaA play a role in development and SM biosynthesis in A. niger FGSC A1279. Histone deacetylases offer new strategies for regulation of SM synthesis.

Increased metabolite production by deletion of an HDA1-type histone deacetylase in the phytopathogenic fungi, Magnaporthe oryzae (Pyricularia oryzae) and Fusarium asiaticum

Histone deacetylases (HDACs) play an important role in the regulation of chromatin structure and gene expression. We found that dark pigmentation of Magnaporthe oryzae (anamorph Pyricularia oryzae) ΔMohda1, a mutant strain in which an orthologue of the yeast HDA1 was disrupted by double cross-over homologous recombination, was significantly stimulated in liquid culture. Analysis of metabolites in a ΔMohda1 mutant culture revealed that the accumulation of shunt products of the 1,8-dihydroxynaphthalene melanin and ergosterol pathways were significantly enhanced compared to the wild-type strain. Northern blot analysis of the ΔMohda1 mutant revealed transcriptional activation of three melanin genes that are dispersed throughout the genome of M. oryzae. The effect of deletion of the yeast HDA1 orthologue was also observed in Fusarium asiaticum from the Fusarium graminearum species complex; the HDF2 deletion mutant produced increased levels of nivalenol-type trichothecenes. These results suggest that histone modification via HDA1-type HDAC regulates the production of natural products in filamentous fungi.

SIGNIFICANCE AND IMPACT OF THE STUDY

Natural products of fungi have significant impacts on human welfare, in both detrimental and beneficial ways. Although HDA1-type histone deacetylase is not essential for vegetative growth, deletion of the gene affects the expression of clustered secondary metabolite genes in some fungi. Here, we report that such phenomena are also observed in physically unlinked genes required for melanin biosynthesis in the rice blast fungus. In addition, production of Fusarium trichothecenes, previously reported to be unaffected by HDA1 deletion, was significantly upregulated in another Fusarium species. Thus, the HDA1-inactivation strategy may be regarded as a general approach for overproduction and/or discovery of fungal metabolites.

FvSet2 regulates fungal growth, pathogenicity, and secondary metabolism in Fusarium verticillioides

Histone H3 lysine 36 methylation (H3K36me) is generally associated with activation of gene expression in most eukaryotic cells. However, the function of H3K36me in filamentous fungi is largely unknown. Set2 is the sole lysine histone methyltransferase (KHMTase) enzyme responsible for the methylation of H3K36 in Saccharomyces cerevisiae. In the current study, we identified a single ortholog of S. cerevisiae Set2 in Fusarium verticillioides. We report that FvSet2 is responsible for the trimethylation of H3K36 (H3K36me3). The FvSET2 deletion mutant (ΔFvSet2) showed significant defects in vegetative growth, FB₁ biosynthesis, pigmentation, and fungal virulence. Furthermore, trimethylation of H3K36 was found to be important for active transcription of genes involved in FB₁ and bikaverin biosyntheses. These data indicate that FvSet2 plays an important role in the regulation of secondary metabolism, vegetative growth and fungal virulence in F. verticillioides.

Global mapping of the regulatory interactions of histone residues

Histone residues can serve as platforms for specific regulatory function. Here we constructed a map of regulatory associations between histone residues and a wide spectrum of chromatin regulation factors based on gene expression changes by histone point mutations in Saccharomyces cerevisiae. Detailed analyses of this map revealed novel associations. Regarding the modulation of H3K4 and K36 methylation by Set1, Set2, or Jhd2, we proposed a role for H4K91 acetylation in early Pol II elongation, and for H4K16 deacetylation in late elongation and crosstalk with H3K4 demethylation for gene silencing. The association of H3K56 with nucleosome positioning suggested that this lysine residue and its acetylation might contribute to nucleosome mobility for transcription activation. Further insights into chromatin regulation are expected from this approach.

Contribution of Sequence Motif, Chromatin State, and DNA Structure Features to Predictive Models of Transcription Factor Binding in Yeast

Transcription factor (TF) binding is determined by the presence of specific sequence motifs (SM) and chromatin accessibility, where the latter is influenced by both chromatin state (CS) and DNA structure (DS) properties. Although SM, CS, and DS have been used to predict TF binding sites, a predictive model that jointly considers CS and DS has not been developed to predict either TF-specific binding or general binding properties of TFs. Using budding yeast as model, we found that machine learning classifiers trained with either CS or DS features alone perform better in predicting TF-specific binding compared to SM-based classifiers. In addition, simultaneously considering CS and DS further improves the accuracy of the TF binding predictions, indicating the highly complementary nature of these two properties. The contributions of SM, CS, and DS features to binding site predictions differ greatly between TFs, allowing TF-specific predictions and potentially reflecting different TF binding mechanisms. In addition, a "TF-agnostic" predictive model based on three DNA “intrinsic properties” (in silico predicted nucleosome occupancy, major groove geometry, and dinucleotide free energy) that can be calculated from genomic sequences alone has performance that rivals the model incorporating experiment-derived data. This intrinsic property model allows prediction of binding regions not only across TFs, but also across DNA-binding domain families with distinct structural folds. Furthermore, these predicted binding regions can help identify TF binding sites that have a significant impact on target gene expression. Because the intrinsic property model allows prediction of binding regions across DNA-binding domain families, it is TF agnostic and likely describes general binding potential of TFs. Thus, our findings suggest that it is feasible to establish a TF agnostic model for identifying functional regulatory regions in potentially any sequenced genome.

Identification of transcription factor binding sites based on sequence motifs is typically accompanied by a high false positive rate. Increasing evidence suggests that there are many other factors besides DNA sequence that may affect the binding and interaction of TFs with DNA. Through the integration of sequence motif, chromatin state, and DNA structure properties, we show that TF binding can be better predicted. Moreover, considering chromatin state and DNA structure properties simultaneously yields a significant improvement. While the binding of some TFs can be readily predicted using either chromatin state information or DNA structure, other TFs need both. Thus, our findings provide insights on how different histone modifications and DNA structure properties may influence the binding of a particular TF and thus how TFs regulate gene expression. These features are referred to as sequence “intrinsic properties” because they can be predicted from sequences alone. These intrinsic properties can be used to build a TF binding prediction model that has a similar performance to considering all features. Moreover, the intrinsic property model allows TFBS predictions not only across TFs, but also across DNA-binding domain families that are present in most eukaryotes, suggesting that the model likely can be used across species.

The conserved histone deacetylase Rpd3 and its DNA binding subunit Ume6 control dynamic transcript architecture during mitotic growth and meiotic development

It was recently reported that the sizes of many mRNAs change when budding yeast cells exit mitosis and enter the meiotic differentiation pathway. These differences were attributed to length variations of their untranslated regions. The function of UTRs in protein translation is well established. However, the mechanism controlling the expression of distinct transcript isoforms during mitotic growth and meiotic development is unknown. In this study, we order developmentally regulated transcript isoforms according to their expression at specific stages during meiosis and gametogenesis, as compared to vegetative growth and starvation. We employ regulatory motif prediction, in vivo protein-DNA binding assays, genetic analyses and monitoring of epigenetic amino acid modification patterns to identify a novel role for Rpd3 and Ume6, two components of a histone deacetylase complex already known to repress early meiosis-specific genes in dividing cells, in mitotic repression of meiosis-specific transcript isoforms. Our findings classify developmental stage-specific early, middle and late meiotic transcript isoforms, and they point to a novel HDAC-dependent control mechanism for flexible transcript architecture during cell growth and differentiation. Since Rpd3 is highly conserved and ubiquitously expressed in many tissues, our results are likely relevant for development and disease in higher eukaryotes.

Spatial localization of co-regulated genes exceeds genomic gene clustering in the Saccharomyces cerevisiae genome

While it has been long recognized that genes are not randomly positioned along the genome, the degree to which its 3D structure influences the arrangement of genes has remained elusive. In particular, several lines of evidence suggest that actively transcribed genes are spatially co-localized, forming transcription factories; however, a generalized systematic test has hitherto not been described. Here we reveal transcription factories using a rigorous definition of genomic structure based on Saccharomyces cerevisiae chromosome conformation capture data, coupled with an experimental design controlling for the primary gene order. We develop a data-driven method for the interpolation and the embedding of such datasets and introduce statistics that enable the comparison of the spatial and genomic densities of genes. Combining these, we report evidence that co-regulated genes are clustered in space, beyond their observed clustering in the context of gene order along the genome and show this phenomenon is significant for 64 out of 117 transcription factors. Furthermore, we show that those transcription factors with high spatially co-localized targets are expressed higher than those whose targets are not spatially clustered. Collectively, our results support the notion that, at a given time, the physical density of genes is intimately related to regulatory activity.

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

Transcription factors (TFs) and histone octamers are two abundant classes of DNA binding proteins that coordinate the transcriptional program in cells. Detailed studies of individual TFs have shown that TFs bind to nucleosome-occluded DNA sequences and induce nucleosome disruption/repositioning, while recent global studies suggest this is not the only mechanism used by all TFs. We have analyzed to what extent the intrinsic DNA binding preferences of TFs and histones play a role in determining nucleosome occupancy, in addition to nonintrinsic factors such as the enzymatic activity of chromatin remodelers. The majority of TFs in budding yeast have an intrinsic sequence preference overlapping with nucleosomal histones. TFs with intrinsic DNA binding properties highly correlated with those of histones tend to be associated with gene activation and might compete with histones to bind to genomic DNA. Consistent with this, we show that activators induce more nucleosome disruption upon transcriptional activation than repressors.

Characterization of Chromatin Structure-associated Histone Modifications in Breast Cancer Cells

Chromatin structure and dynamics that are influenced by epigenetic marks, such as histone modification and DNA methylation, play a crucial role in modulating gene transcription. To understand the relationship between histone modifications and regulatory elements in breast cancer cells, we compared our chromatin immunoprecipitation sequencing (ChIP-Seq) histone modification patterns for histone H3K4me1, H3K4me3, H3K9/16ac, and H3K27me3 in MCF-7 cells with publicly available formaldehyde-assisted isolation of regulatory elements (FAIRE)-chip signals in human chromosomes 8, 11, and 12, identified by a method called FAIRE. Active regulatory elements defined by FAIRE were highly associated with active histone modifications, like H3K4me3 and H3K9/16ac, especially near transcription start sites. The H3K9/16ac-enriched genes that overlapped with FAIRE signals (FAIRE-H3K9/14ac) were moderately correlated with gene expression levels. We also identified functional sequence motifs at H3K4me1-enriched FAIRE sites upstream of putative promoters, suggesting that regulatory elements could be associated with H3K4me1 to be regarded as distal regulatory elements. Our results might provide an insight into epigenetic regulatory mechanisms explaining the association of histone modifications with open chromatin structure in breast cancer cells.

Balancing noise and plasticity in eukaryotic gene expression

BACKGROUND

Coupling the control of expression stochasticity (noise) to the ability of expression change (plasticity) can alter gene function and influence adaptation. A number of factors, such as transcription re-initiation, strong chromatin regulation or genome neighboring organization, underlie this coupling. However, these factors do not necessarily combine in equivalent ways and strengths in all genes. Can we identify then alternative architectures that modulate in distinct ways the linkage of noise and plasticity?

RESULTS

Here we first show that strong chromatin regulation, commonly viewed as a source of coupling, can lead to plasticity without noise. The nature of this regulation is relevant too, with plastic but noiseless genes being subjected to general activators whereas plastic and noisy genes experience more specific repression. Contrarily, in genes exhibiting poor transcriptional control, it is translational efficiency what separates noise from plasticity, a pattern related to transcript length. This additionally implies that genome neighboring organization -as modifier- appears only effective in highly plastic genes. In this class, we confirm bidirectional promoters (bipromoters) as a configuration capable to reduce coupling by abating noise but also reveal an important trade-off, since bipromoters also decrease plasticity. This presents ultimately a paradox between intergenic distances and modulation, with short intergenic distances both associated and disassociated to noise at different plasticity levels.

CONCLUSIONS

Balancing the coupling among different types of expression variability appears as a potential shaping force of genome regulation and organization. This is reflected in the use of different control strategies at genes with different sets of functional constraints.

Histone modification pattern evolution after yeast gene duplication

Background

Gene duplication and subsequent functional divergence especially expression divergence have been widely considered as main sources for evolutionary innovations. Many studies evidenced that genetic regulatory network evolved rapidly shortly after gene duplication, thus leading to accelerated expression divergence and diversification. However, little is known whether epigenetic factors have mediated the evolution of expression regulation since gene duplication. In this study, we conducted detailed analyses on yeast histone modification (HM), the major epigenetics type in this organism, as well as other available functional genomics data to address this issue.

Results

Duplicate genes, on average, share more common HM-code patterns than random singleton pairs in their promoters and open reading frames (ORF). Though HM-code divergence between duplicates in both promoter and ORF regions increase with their sequence divergence, the HM-code in ORF region evolves slower than that in promoter region, probably owing to the functional constraints imposed on protein sequences. After excluding the confounding effect of sequence divergence (or evolutionary time), we found the evidence supporting the notion that in yeast, the HM-code may co-evolve with cis- and trans-regulatory factors. Moreover, we observed that deletion of some yeast HM-related enzymes increases the expression divergence between duplicate genes, yet the effect is lower than the case of transcription factor (TF) deletion or environmental stresses.

Conclusions

Our analyses demonstrate that after gene duplication, yeast histone modification profile between duplicates diverged with evolutionary time, similar to genetic regulatory elements. Moreover, we found the evidence of the co-evolution between genetic and epigenetic elements since gene duplication, together contributing to the expression divergence between duplicate genes.

Mcm2 phosphorylation and the response to replicative stress

Background

The replicative helicase in eukaryotic cells is comprised of minichromosome maintenance (Mcm) proteins 2 through 7 (Mcm2-7) and is a key target for regulation of cell proliferation. In addition, it is regulated in response to replicative stress. One of the protein kinases that targets Mcm2-7 is the Dbf4-dependent kinase Cdc7 (DDK). In a previous study, we showed that alanine mutations of the DDK phosphorylation sites at S164 and S170 in Saccharomyces cerevisiae Mcm2 result in sensitivity to caffeine and methyl methanesulfonate (MMS) leading us to suggest that DDK phosphorylation of Mcm2 is required in response to replicative stress.

Results

We show here that a strain with the mcm2 allele lacking DDK phosphorylation sites (mcm2_AA) is also sensitive to the ribonucleotide reductase inhibitor, hydroxyurea (HU) and to the base analogue 5-fluorouracil (5-FU) but not the radiomimetic drug, phleomycin. We screened the budding yeast non-essential deletion collection for synthetic lethal interactions with mcm2_AA and isolated deletions that include genes involved in the control of genome integrity and oxidative stress. In addition, the spontaneous mutation rate, as measured by mutations in CAN1, was increased in the mcm2_AA strain compared to wild type, whereas with a phosphomimetic allele (mcm2_EE) the mutation rate was decreased. These results led to the idea that the mcm2_AA strain is unable to respond properly to DNA damage. We examined this by screening the deletion collection for suppressors of the caffeine sensitivity of mcm2_AA. Deletions that decrease spontaneous DNA damage, increase homologous recombination or slow replication forks were isolated. Many of the suppressors of caffeine sensitivity suppressed other phenotypes of mcm2_AA including sensitivity to genotoxic drugs, the increased frequency of cells with RPA foci and the increased mutation rate.

Conclusions

Together these observations point to a role for DDK-mediated phosphorylation of Mcm2 in the response to replicative stress, including some forms of DNA damage. We suggest that phosphorylation of Mcm2 modulates Mcm2-7 activity resulting in the stabilization of replication forks in response to replicative stress.

A genomic model of condition-specific nucleosome behavior explains transcriptional activity in yeast

Nucleosomes play an important role in gene regulation. Molecular studies observed that nucleosome binding in promoters tends to be repressive. In contrast, genomic studies have delivered conflicting results: An analysis of yeast grown on diverse carbon sources reported that nucleosome occupancies remain largely unchanged between conditions, whereas a study of the heat-shock response suggested that nucleosomes get evicted at promoters of genes with increased expression. Consequently, there are few general principles that capture the relationship between chromatin organization and transcriptional regulation. Here, we present a qualitative model for nucleosome positioning in Saccharomyces cerevisiae that helps explain important properties of gene expression. By integrating publicly available data sets, we observe that promoter-bound nucleosomes assume one of four discrete configurations that determine the active and silent transcriptional states of a gene, but not its expression level. In TATA-box-containing promoters, nucleosome architecture indicates the amount of transcriptional noise. We show that >20% of genes switch promoter states upon changes in cellular conditions. The data suggest that DNA-binding transcription factors together with chromatin-remodeling enzymes are primarily responsible for the nucleosome architecture. Our model for promoter nucleosome architecture reconciles genome-scale findings with molecular studies; in doing so, we establish principles for nucleosome positioning and gene expression that apply not only to individual genes, but across the entire genome. The study provides a stepping stone for future models of transcriptional regulation that encompass the intricate interplay between cis- and trans-acting factors, chromatin, and the core transcriptional machinery.

Genome-wide analysis of chromatin features identifies histone modification sensitive and insensitive yeast transcription factors

We propose a method to predict yeast transcription factor targets by integrating histone modification profiles with transcription factor binding motif information. It shows improved predictive power compared to a binding motif-only method. We find that transcription factors cluster into histone-sensitive and -insensitive classes. The target genes of histone-sensitive transcription factors have stronger histone modification signals than those of histone-insensitive ones. The two classes also differ in tendency to interact with histone modifiers, degree of connectivity in protein-protein interaction networks, position in the transcriptional regulation hierarchy, and in a number of additional features, indicating possible differences in their transcriptional regulation mechanisms.

Distinct role of Mediator tail module in regulation of SAGA-dependent, TATA-containing genes in yeast

The evolutionarily conserved Mediator complex is required for transcription of nearly all RNA Pol II-dependent promoters, with the tail module serving to recruit Mediator to active promoters in current models. However, transcriptional dependence on tail module subunits varies in a gene-specific manner, and the generality of the tail module requirement for transcriptional activation has not been explored. Here, we show that tail module subunits function redundantly to recruit Mediator to promoters in yeast, and transcriptome analysis shows stronger effects on genome-wide expression in a double-tail subunit deletion mutant than in single-subunit deletion mutants. Unexpectedly, TATA-containing and SAGA-dependent genes were much more affected by impairment of tail module function than were TFIID-dependent genes. Consistent with this finding, Mediator and preinitiation complex association with SAGA-dependent promoters is substantially reduced in gal11/med15Δ med3Δ yeast, whereas association of TBP, Pol II, and other Mediator modules with TFIID-dependent genes is largely independent of the tail module. Thus, we have identified a connection between the Mediator tail module and the division of promoter dependence between TFIID and SAGA.

Nuclear colocalization of transcription factor target genes strengthens coregulation in yeast

Eukaryotic chromosomes are not randomly distributed in the interphase nucleus, but instead occupy distinct territories. Nonetheless, the genome-wide relationships of gene regulation to gene nuclear location remain poorly understood in yeast. In the three-dimensional view of gene regulation, we found that a considerable number of transcription factors (TFs) regulate genes that are colocalized in the nucleus. Colocalized TF target genes are more strongly coregulated compared with the other TF target genes. Target genes of chromatin regulators are also colocalized. These results demonstrate that colocalization of coregulated genes is a common process, and three-dimensional gene positioning is an important part of gene regulation. Our findings will have implications in understanding nuclear architecture and function.

The stress response factors Yap6, Cin5, Phd1, and Skn7 direct targeting of the conserved co-repressor Tup1-Ssn6 in S. cerevisiae

Maintaining the proper expression of the transcriptome during development or in response to a changing environment requires a delicate balance between transcriptional regulators with activating and repressing functions. The budding yeast transcriptional co-repressor Tup1-Ssn6 is a model for studying similar repressor complexes in multicellular eukaryotes. Tup1-Ssn6 does not bind DNA directly, but is directed to individual promoters by one or more DNA-binding proteins, referred to as Tup1 recruiters. This functional architecture allows the Tup1-Ssn6 to modulate the expression of genes required for the response to a variety of cellular stresses. To understand the targeting or the Tup1-Ssn6 complex, we determined the genomic distribution of Tup1 and Ssn6 by ChIP-chip. We found that most loci bound by Tup1-Ssn6 could not be explained by co-occupancy with a known recruiting cofactor and that deletion of individual known Tup1 recruiters did not significantly alter the Tup1 binding profile. These observations suggest that new Tup1 recruiting proteins remain to be discovered and that Tup1 recruitment typically depends on multiple recruiting cofactors. To identify new recruiting proteins, we computationally screened for factors with binding patterns similar to the observed Tup1-Ssn6 genomic distribution. Four top candidates, Cin5, Skn7, Phd1, and Yap6, all known to be associated with stress response gene regulation, were experimentally confirmed to physically interact with Tup1 and/or Ssn6. Incorporating these new recruitment cofactors with previously characterized cofactors now explains the majority of Tup1 targeting across the genome, and expands our understanding of the mechanism by which Tup1-Ssn6 is directed to its targets.

Transcriptional regulation via TF-modifying enzymes: an integrative model-based analysis

Transcription factor activity is largely regulated through post-translational modification. Here, we report the first integrative model of transcription that includes both interactions between transcription factors and promoters, and between transcription factors and modifying enzymes. Simulations indicate that our method is robust against noise. We validated our tool on a well-studied stress response network in yeast and on a STAT1-mediated regulatory network in human B cells. Our work represents a significant step toward a comprehensive model of gene transcription.

Gene clustering pattern, promoter architecture, and gene expression stability in eukaryotic genomes

A balance between gene expression stability and evolvability is essential for the long-term maintenance of a living system. In this paper, we studied whether the genetic and epigenetic properties of the promoter affect gene expression variability. We hypothesized that upstream distance and orientation (head-to-head or head-to-tail) are important for the promoter architecture and gene expression variability. We found that in budding yeast genes with a short upstream distance tend to have low gene expression variability, and their promoter is flanked by strongly positioned nucleosomes and tends to have low nucleosome occupancy. These observations suggest that in vivo positioning of the flanking nucleosomes facilitates stable nucleosome depletion at the core promoter region and enhances gene expression stability. Head-to-head genes have, on average, lower gene expression variability, greater nucleosome depletion at the core promoter region, and more strongly positioned nucleosomes that flank the core promoter than do head-to-tail genes. These observations hold for diverse eukaryotes. In complex organisms such as mammals, only a small fraction of head-to-tail genes have retained a short upstream distance, probably because the promoter may not be flanked by a strongly positioned nucleosome on the upstream side.

Coordinated histone modifications are associated with gene expression variation within and between species

Histone modifications regulate gene expression in eukaryotes, but their effects on transcriptomes of a multicellular organism and on transcriptomic divergence between species are poorly understood. Here we present the first nucleotide-resolution maps of histone acetylation, methylation, and core histone in Arabidopsis thaliana and a comprehensive analysis of these and all other available maps with gene expression data in A. thaliana, Arabidopsis arenosa, and allotetraploids. H3K9 acetylation (H3K9ac) and H3K4 trimethylation (H3K4me3) are correlated, and their distribution patterns are associated with Gene Ontology (GO) functional classifications. Highly dense and narrow distributions of these modifications near transcriptional start sites are associated with constitutive expression of genes involved in translation, whereas broad distributions toward coding regions correlate with expression variation of the genes involved in photosynthesis, carbohydrate metabolism, and defense responses. Compared to animal stem cells, dispersed distributions of H3K27me3 without bivalent H3K4me3 and H3K9ac marks correlate with developmentally repressed genes in Arabidopsis. Finally, genes affected by A. thaliana histone deacetylase 1 mutation tend to show high levels of expression variation within and between species. The data suggest that genome-wide coordinated modifications of histone acetylation and methylation provide a general mechanism for gene expression changes within and between species and in allopolyploids.

Differential effects of chromatin regulators and transcription factors on gene regulation: a nucleosomal perspective

MOTIVATION

Chromatin regulators (CR) and transcription factors (TF) are important trans-acting factors regulating transcription process, and many efforts have been devoted to understand their underlying mechanisms in gene regulation. However, the influences of CR and TF regulation effects on nucleosomes during transcription are still minimally understood, and it remains to be determined the extent to which CR and TF regulatory effect shape the organization of nucleosomes in the genome. In this article we attempted to address this problem and examine the patterns of CR and TF regulation effects from the nucleosome perspective.

RESULTS

Our results show that the CR and TF regulatory effects exhibit different paradigms of transcriptional control in Saccharomyces cerevisiae. We grouped yeast genes into two categories, 'CR-sensitive' genes and 'TF-sensitive' genes, based on how their expression profiles change upon deletion of CRs or TFs. We found that genes in these two groups have very different patterns of nucleosome organization. The promoters of CR-sensitive genes tend to have higher nucleosome occupancy, whereas the promoters of TF-sensitive genes are depleted of nucleosomes. Furthermore, the nucleosome profiles of CR-sensitive genes tend to show more dynamic characteristics than TF-sensitive genes. These results reveal that the nucleosome organizations of yeast genes have a strong impact on their mode of regulation, and there are differential regulation effects on nucleosomes between CRs and TFs.

AVAILABILITY

http://www.utoronto.ca/zhanglab/papers/bioinfo_2010/.

Identifying the genetic determinants of transcription factor activity

Genome-wide messenger RNA expression levels are highly heritable. However, the molecular mechanisms underlying this heritability are poorly understood.

The influence of trans-acting polymorphisms is often mediated by changes in the regulatory activity of one or more sequence-specific transcription factors (TFs). We use a method that exploits prior information about the DNA-binding specificity of each TF to estimate its genotype-specific regulatory activity. To this end, we perform linear regression of genotype-specific differential mRNA expression on TF-specific promoter-binding affinity.

Treating inferred TF activity as a quantitative trait and mapping it across a panel of segregants from an experimental genetic cross allows us to identify trans-acting loci (‘aQTLs') whose allelic variation modulates the TF. A few of these aQTL regions contain the gene encoding the TF itself; several others contain a gene whose protein product is known to interact with the TF.

Our method is strictly causal, as it only uses sequence-based features as predictors. Application to budding yeast demonstrates a dramatic increase in statistical power, compared with existing methods, to detect locus-TF associations and trans-acting loci. Our aQTL mapping strategy also succeeds in mouse.

Genetic sequence variation naturally perturbs mRNA expression levels in the cell. In recent years, analysis of parallel genotyping and expression profiling data for segregants from genetic crosses between parental strains has revealed that mRNA expression levels are highly heritable. Expression quantitative trait loci (eQTLs), whose allelic variation regulates the expression level of individual genes, have successfully been identified (Brem et al, 2002; Schadt et al, 2003). The molecular mechanisms underlying the heritability of mRNA expression are poorly understood. However, they are likely to involve mediation by transcription factors (TFs). We present a new transcription-factor-centric method that greatly increases our ability to understand what drives the genetic variation in mRNA expression (Figure 1). Our method identifies genomic loci (‘aQTLs') whose allelic variation modulates the protein-level activity of specific TFs. To map aQTLs, we integrate genotyping and expression profiling data with quantitative prior information about DNA-binding specificity of transcription factors in the form of position-specific affinity matrices (Bussemaker et al, 2007). We applied our method in two different organisms: budding yeast and mouse.

In our approach, the inferred TF activity is explicitly treated as a quantitative trait, and genetically mapped. The decrease of ‘phenotype space' from that of all genes (in the eQTL approach) to that of all TFs (in our aQTL approach) increases the statistical power to detect trans-acting loci in two distinct ways. First, as each inferred TF activity is derived from a large number of genes, it is far less noisy than mRNA levels of individual genes. Second, the number of trait/marker combinations that needs to be tested for statistical significance in parallel is roughly two orders of magnitude smaller than for eQTLs. We identified a total of 103 locus-TF associations, a more than six-fold improvement over the 17 locus-TF associations identified by several existing methods (Brem et al, 2002; Yvert et al, 2003; Lee et al, 2006; Smith and Kruglyak, 2008; Zhu et al, 2008). The total number of distinct genomic loci identified as an aQTL equals 31, which includes 11 of the 13 previously identified eQTL hotspots (Smith and Kruglyak, 2008).

To better understand the mechanisms underlying the identified genetic linkages, we examined the genes within each aQTL region. First, we found four ‘local' aQTLs, which encompass the gene encoding the TF itself. This includes the known polymorphism in the HAP1 gene (Brem et al, 2002), but also novel predictions of trans-acting polymorphisms in RFX1, STB5, and HAP4. Second, using high-throughput protein–protein interaction data, we identified putative causal genes for several aQTLs. For example, we predict that a polymorphism in the cyclin-dependent kinase CDC28 antagonistically modulates the functionally distinct cell cycle regulators Fkh1 and Fkh2. In this and other cases, our approach naturally accounts for post-translational modulation of TF activity at the protein level.

We validated our ability to predict locus-TF associations in yeast using gene expression profiles of allele replacement strains from a previous study (Smith and Kruglyak, 2008). Chromosome 15 contains an aQTL whose allelic status influences the activity of no fewer than 30 distinct TFs. This locus includes IRA2, which controls intracellular cAMP levels. We used the gene expression profile of IRA2 replacement strains to confirm that the polymorphism within IRA2 indeed modulates a subset of the TFs whose activity was predicted to link to this locus, and no other TFs.

Application of our approach to mouse data identified an aQTL modulating the activity of a specific TF in liver cells. We identified an aQTL on mouse chromosome 7 for Zscan4, a transcription factor containing four zinc finger domains and a SCAN domain. Even though we could not detect a candidate causal gene for Zscan4p because of lack of information about the mouse genome, our result demonstrates that our method also works in higher eukaryotes.

In summary, aQTL mapping has a greatly improved sensitivity to detect molecular mechanisms underlying the heritability of gene expression. The successful application of our approach to yeast and mouse data underscores the value of explicitly treating the inferred TF activity as a quantitative trait for increasing statistical power of detecting trans-acting loci. Furthermore, our method is computationally efficient, and easily applicable to any other organism whenever prior information about the DNA-binding specificity of TFs is available.

Analysis of parallel genotyping and expression profiling data has shown that mRNA expression levels are highly heritable. Currently, only a tiny fraction of this genetic variance can be mechanistically accounted for. The influence of trans-acting polymorphisms on gene expression traits is often mediated by transcription factors (TFs). We present a method that exploits prior knowledge about the in vitro DNA-binding specificity of a TF in order to map the loci (‘aQTLs') whose inheritance modulates its protein-level regulatory activity. Genome-wide regression of differential mRNA expression on predicted promoter affinity is used to estimate segregant-specific TF activity, which is subsequently mapped as a quantitative phenotype. In budding yeast, our method identifies six times as many locus-TF associations and more than twice as many trans-acting loci as all existing methods combined. Application to mouse data from an F2 intercross identified an aQTL on chromosome VII modulating the activity of Zscan4 in liver cells. Our method has greatly improved statistical power over existing methods, is mechanism based, strictly causal, computationally efficient, and generally applicable.

Evidences for increased expression variation of duplicate genes in budding yeast: from cis- to trans-regulation effects

Duplicate genes tend to have a more variable expression program than singleton genes, which was thought to be an important way for the organism to respond and adapt to fluctuating environment. However, the underlying molecular mechanisms driving such expression variation remain largely unexplored. In this work, we first rigorously confirmed that duplicate genes indeed have higher gene expression variation than singleton genes in several aspects, i.e. responses to environmental perturbation, between-strain divergence, and expression noise. To investigate the underlying mechanism, we further analyzed a previously published expression dataset of yeast segregants produced from genetic crosses. We dissected the observed expression divergence between segregant strains into cis- and trans-variabilities, and demonstrated that trans-regulation effect can explain larger fraction of the expression variation than cis-regulation effect. This is true for both duplicate genes and singleton genes. In contrast, we found, between a pair of sister paralogs, cis-variability explains more of the expression divergence between the paralogs than trans-variability. We next investigated the presence of cis- and trans-features that are associated with elevated expression variations. For cis-acting regulation, duplicate genes have higher genetic diversity in their promoters and coding regions than singleton genes. For trans-acting regulation, duplicate and singleton genes are differentially regulated by chromatin regulators and transcription factors, and duplicate genes are more severely affected by the deletion of histone tails. These results showed that both cis-and trans-factors have great effect in causing the increased expression variation of duplicate genes, and explained the previously observed differences in transcription regulation between duplicate genes and singleton genes.

Gene expression variations are predictive for stochastic noise

Fluctuations in protein abundance among single cells are primarily due to the inherent stochasticity in transcription and translation processes, such stochasticity can often confer phenotypic heterogeneity among isogenic cells. It has been proposed that expression noise can be triggered as an adaptation to environmental stresses and genetic perturbations, and as a mechanism to facilitate gene expression evolution. Thus, elucidating the relationship between expression noise, measured at the single-cell level, and expression variation, measured on population of cells, can improve our understanding on the variability and evolvability of gene expression. Here, we showed that noise levels are significantly correlated with conditional expression variations. We further demonstrated that expression variations are highly predictive for noise level, especially in TATA-box containing genes. Our results suggest that expression variabilities can serve as a proxy for noise level, suggesting that these two properties share the same underlining mechanism, e.g. chromatin regulation. Our work paves the way for the study of stochastic noise in other single-cell organisms.

Identifying the combinatorial effects of histone modifications by association rule mining in yeast

Eukaryotic genomes are packaged into chromatin by histone proteins whose chemical modification can profoundly influence gene expression. The histone modifications often act in combinations, which exert different effects on gene expression. Although a number of experimental techniques and data analysis methods have been developed to study histone modifications, it is still very difficult to identify the relationships among histone modifications on a genome-wide scale.We proposed a method to identify the combinatorial effects of histone modifications by association rule mining. The method first identified Functional Modification Transactions (FMTs) and then employed association rule mining algorithm and statistics methods to identify histone modification patterns. We applied the proposed methodology to Pokholok et al's data with eight sets of histone modifications and Kurdistani et al's data with eleven histone acetylation sites. Our method succeeds in revealing two different global views of histone modification landscapes on two datasets and identifying a number of modification patterns some of which are supported by previous studies.We concentrate on combinatorial effects of histone modifications which significantly affect gene expression. Our method succeeds in identifying known interactions among histone modifications and uncovering many previously unknown patterns. After in-depth analysis of possible mechanism by which histone modification patterns can alter transcriptional states, we infer three possible modification pattern reading mechanism ('redundant', 'trivial', 'dominative'). Our results demonstrate several histone modification patterns which show significant correspondence between yeast and human cells.

New insights into two distinct nucleosome distributions: comparison of cross-platform positioning datasets in the yeast genome

Background

Recently, a number of high-resolution genome-wide maps of nucleosome locations in S. cerevisiae have been derived experimentally. However, nucleosome positions are determined in vivo by the combined effects of numerous factors. Consequently, nucleosomes are not simple static units, which may explain the discrepancies in reported nucleosome positions as measured by different experiments. In order to more accurately depict the genome-wide nucleosome distribution, we integrated multiple nucleosomal positioning datasets using a multi-angle analysis strategy.

Results

To evaluate the contribution of chromatin structure to transcription, we used the vast amount of available nucleosome analyzed data. Analysis of this data allowed for the comprehensive identification of the connections between promoter nucleosome positioning patterns and various transcription-dependent properties. Further, we characterised the function of nucleosome destabilisation in the context of transcription regulation. Our results indicate that genes with similar nucleosome occupancy patterns share general transcription attributes. We identified the local regulatory correlation (LRC) regions for two distinct types of nucleosomes and we assessed their regulatory properties. We also estimated the nucleosome reproducibility and measurement accuracy for high-confidence transcripts. We found that by maintaining a distance of ~13 bp between the upstream border of the +1 nucleosome and the transcription start sites (TSSs), the stable +1 nucleosome may form a barrier against the accessibility of the TSS and shape an optimum chromatin conformation for gene regulation. An in-depth analysis of nucleosome positioning in normally growing and heat shock cells suggested that the extent and patterns of nucleosome sliding are associated with gene activation.

Conclusions

Our results, which combine different types of data, suggest that cross-platform information, including discrepancy and consistency, reflects the mechanisms of nucleosome packaging in vivo more faithfully than individual studies. Furthermore, nucleosomes can be divided into two classes according to their stable and dynamic characteristics. We found that two different nucleosome-positioning characteristics may significantly impact transcription programs. Besides, some positioned-nucleosomes are involved in the transition from stable state to dynamic state in response to abrupt environmental changes.

Mimosa: mixture model of co-expression to detect modulators of regulatory interaction

BACKGROUND

Functionally related genes tend to be correlated in their expression patterns across multiple conditions and/or tissue-types. Thus co-expression networks are often used to investigate functional groups of genes. In particular, when one of the genes is a transcription factor (TF), the co-expression-based interaction is interpreted, with caution, as a direct regulatory interaction. However, any particular TF, and more importantly, any particular regulatory interaction, is likely to be active only in a subset of experimental conditions. Moreover, the subset of expression samples where the regulatory interaction holds may be marked by presence or absence of a modifier gene, such as an enzyme that post-translationally modifies the TF. Such subtlety of regulatory interactions is overlooked when one computes an overall expression correlation.

RESULTS

Here we present a novel mixture modeling approach where a TF-Gene pair is presumed to be significantly correlated (with unknown coefficient) in an (unknown) subset of expression samples. The parameters of the model are estimated using a Maximum Likelihood approach. The estimated mixture of expression samples is then mined to identify genes potentially modulating the TF-Gene interaction. We have validated our approach using synthetic data and on four biological cases in cow, yeast, and humans.

CONCLUSIONS

While limited in some ways, as discussed, the work represents a novel approach to mine expression data and detect potential modulators of regulatory interactions.

Regulating the regulators: modulators of transcription factor activity

Gene transcription is largely regulated by DNA-binding transcription factors (TFs). However, the TF activity itself is modulated via, among other things, post-translational modifications (PTMs) by specific modification enzymes in response to cellular stimuli. TF-PTMs thus serve as "molecular switchboards" that map upstream signaling events to the downstream transcriptional events. An important long-term goal is to obtain a genome-wide map of "regulatory triplets" consisting of a TF, target gene, and a modulator gene that specifically modulates the regulation of the target gene by the TF. A variety of genome-wide data sets can be exploited by computational methods to obtain a rough map of regulatory triplets, which can guide directed experiments. However, a prerequisite to developing such computational tools is a systematic catalog of known instances of regulatory triplets. We first describe PTM-Switchboard, a recent database that stores triplets of genes such that the ability of one gene (the TF) to regulate a target gene is dependent on one or more PTMs catalyzed by a third gene, the modifying enzyme. We also review current computational approaches to infer regulatory triplets from genome-wide data sets and conclude with a discussion of potential future research. PTM-Switchboard is accessible at http://cagr.pcbi.upenn.edu/PTMswitchboard /

Nucleosome positioning by genomic excluding-energy barriers

Recent genome-wide nucleosome mappings along with bioinformatics studies have confirmed that the DNA sequence plays a more important role in the collective organization of nucleosomes in vivo than previously thought. Yet in living cells, this organization also results from the action of various external factors like DNA-binding proteins and chromatin remodelers. To decipher the code for intrinsic chromatin organization, there is thus a need for in vitro experiments to bridge the gap between computational models of nucleosome sequence preferences and in vivo nucleosome occupancy data. Here we combine atomic force microscopy in liquid and theoretical modeling to demonstrate that a major sequence signaling in vivo are high-energy barriers that locally inhibit nucleosome formation rather than favorable positioning motifs. We show that these genomic excluding-energy barriers condition the collective assembly of neighboring nucleosomes consistently with equilibrium statistical ordering principles. The analysis of two gene promoter regions in Saccharomyces cerevisiae and the human genome indicates that these genomic barriers direct the intrinsic nucleosome occupancy of regulatory sites, thereby contributing to gene expression regulation.

Bayesian network analysis of targeting interactions in chromatin

In eukaryotes, many chromatin proteins together regulate gene expression. Chromatin proteins often direct the genomic binding pattern of other chromatin proteins, for example, by recruitment or competition mechanisms. The network of such targeting interactions in chromatin is complex and still poorly understood. Based on genome-wide binding maps, we constructed a Bayesian network model of the targeting interactions among a broad set of 43 chromatin components in Drosophila cells. This model predicts many novel functional relationships. For example, we found that the homologous proteins HP1 and HP1C each target the heterochromatin protein HP3 to distinct sets of genes in a competitive manner. We also discovered a central role for the remodeling factor Brahma in the targeting of several DNA-binding factors, including GAGA factor, JRA, and SU(VAR)3-7. Our network model provides a global view of the targeting interplay among dozens of chromatin components.

Chromatin regulation and gene centrality are essential for controlling fitness pleiotropy in yeast

Background

There are a wide range of phenotypes that are due to loss-of-function or null mutations. Previously, the functions of gene products that distinguish essential from nonessential genes were characterized. However, the functions of products of non-essential genes that contribute to fitness remain minimally understood.

Principal Findings

Using data from Saccharomyces cerevisiae, we investigated several gene characteristics, which we are able to measure, that are significantly associated with a gene's fitness pleiotropy. Fitness pleiotropy is a measurement of the gene's importance to fitness. These characteristics include: 1) whether the gene's product functions in chromatin regulation, 2) whether the regulation of the gene is influenced by chromatin state, measured by chromatin regulation effect (CRE), 3) whether the gene's product functions as a transcription factor (TF) and the number of genes a TF regulates, 4) whether the gene contains TATA-box, and 5) whether the gene's product is central in a protein interaction network. Partial correlation analysis was used to study how these characteristics interact to influence fitness pleiotropy. We show that all five characteristics that were measured are statistically significantly associated with fitness pleiotropy. However, fitness pleiotropy is not associated with the presence of TATA-box when CRE is controlled. In particular, two characteristics: 1) whether the regulation of a gene is more likely to be influenced by chromatin state, and 2) whether the gene product is central in a protein interaction network measured by the number of protein interactions were found to play the most important roles affecting a gene's fitness pleiotropy.

Conclusions

These findings highlight the significance of both epigenetic gene regulation and protein interaction networks in influencing the fitness pleiotropy.

A novel strategy of transcription regulation by intragenic nucleosome ordering

Numerous studies of chromatin structure showed that nucleosome free regions (NFRs) located at 5' gene ends contribute to transcription initiation regulation. Here, we determine the role of intragenic chromatin structure on gene expression regulation. We show that, along Saccharomyces cerevisiae genes, nucleosomes are highly organized following two types of architecture that depend only on the distance between the NFRs located at the 5' and 3' gene ends. In the first type, this distance constrains in vivo the positioning of n nucleosomes regularly organized in a "crystal-like" array. In the second type, this distance is such that the corresponding genes can accommodate either n or (n + 1) nucleosomes, thereby displaying two possible crystal-like arrays of n weakly compacted or n + 1 highly compacted nucleosomes. This adaptability confers "bi-stable" properties to chromatin and is a key to its dynamics. Compared to crystal-like genes, bi-stable genes present higher transcriptional plasticity, higher sensitivity to chromatin regulators, higher H3 turnover rate, and lower H2A.Z enrichment. The results strongly suggest that transcription elongation is facilitated by higher chromatin compaction. The data allow us to propose a new paradigm of transcriptional control mediated by the stability and the level of compaction of the intragenic chromatin architecture and open new ways for investigating eukaryotic gene expression regulation.

Properties of untranslated regions of the S. cerevisiae genome

Background

During evolution selection forces such as changing environments shape the architecture of genomes. The distribution of genes along chromosomes and the length of intragenic regions are basic genomic features known to play a major role in the regulation of gene transcription and translation.

Results

In this work we perform the first large scale analysis of the length distribution of untranslated regions (promoters, 5' and 3' untranslated regions, terminators) in the genome of the yeast Saccharomyces cerevisiae. Our analysis shows that the length of each open reading frame (ORF) and that of its associated regulatory and untranslated regions significantly correlate with each other. Moreover, significant correlations with other features related to gene expression and evolution (number of regulating transcription factors, mRNA and protein abundance, evolutionary rate, etc) were observed. Furthermore, the function of genes seems to have an important role in the evolution of these lengths. Notably, genes that are related to RNA metabolism tend to have shorter untranslated regions and thus tend to be closer to their neighbouring genes while genes coding for cell wall proteins tend to be isolated in the genome.

Conclusion

These results indicate that genome architecture has a significant role in regulating gene expression, and in shaping the characteristics and functionality of proteins.

Genome-wide analysis of interactions between ATP-dependent chromatin remodeling and histone modifications

BACKGROUND

ATP-dependent chromatin remodeling and the covalent modification of histones play central roles in determining chromatin structure and function. Although several specific interactions between these two activities have been elaborated, the global landscape remains to be elucidated.

RESULTS

In this paper, we have developed a computational method to generate the first genome-wide landscape of interactions between ATP-dependent chromatin remodeling and the covalent modification of histones in Saccharomyces cerevisiae. Our method succeeds in identifying known interactions and uncovers many previously unknown interactions between these two activities. Analysis of the genome-wide picture revealed that transcription-related modifications tend to interact with more chromatin remodelers. Our results also demonstrate that most chromatin remodeling-modification interactions act via interactions of remodelers with both histone-modifying enzymes and histone residues. We also found that the co-occurrence of both modification and remodeling has significantly different influences on multiple gene features (e.g. nucleosome occupancy) compared with the presence of either one.

CONCLUSION

We gave the first genome-wide picture of ATP-dependent chromatin remodeling-histone modification interactions. We also revealed how these two activities work together to regulate chromatin structure and function. Our results suggest that distinct strategies for regulating chromatin activity are selectively employed by genes with different properties.

HdaA, a class 2 histone deacetylase of Aspergillus fumigatus, affects germination and secondary metabolite production

Histone deacetylases (HDACs) play an important role in regulation of gene expression through histone modifications. Here we show that the Aspergillus fumigatus HDAC HdaA is involved in regulation of secondary metabolite production and is required for normal germination and vegetative growth. Deletion of the hdaA gene increased the production of several secondary metabolites but decreased production of gliotoxin whereas over-expression hdaA increased production of gliotoxin. RT-PCR analysis of 14 nonribosomal peptide synthases indicated HdaA regulation of up to nine of them. A mammalian cell toxicity assay indicated increased activity in the over-expression strain. Neither mutant affected virulence of the fungus as measured by macrophage engulfment of conidia or virulence in a neutropenic mouse model.

Unstable tandem repeats in promoters confer transcriptional evolvability

Relative to most regions of the genome, tandemly repeated DNA sequences display a greater propensity to mutate. A search for tandem repeats in the Saccharomyces cerevisiae genome revealed that the nucleosome-free region directly upstream of genes (the promoter region) is enriched in repeats. As many as 25% of all gene promoters contain tandem repeat sequences. Genes driven by these repeat-containing promoters show significantly higher rates of transcriptional divergence. Variations in repeat length result in changes in expression and local nucleosome positioning. Tandem repeats are variable elements in promoters that may facilitate evolutionary tuning of gene expression by affecting local chromatin structure.

Dynamic reprogramming of transcription factors to and from the subtelomere

Transcription factors are most commonly thought of as proteins that regulate expression of specific genes, independently of the order of those genes along the chromosome. By screening genome-wide chromatin immunoprecipitation (ChIP) profiles in yeast, we find that more than 10% of DNA-binding transcription factors concentrate at the subtelomeric regions near to chromosome ends. None of the proteins identified were previously implicated in regulation at telomeres, yet genomic and proteomic studies reveal that a subset of factors show many interactions with established telomere binding complexes. For many factors, the subtelomeric binding pattern is dynamic and undergoes flux toward or away from the telomere as physiological conditions shift. We find that subtelomeric binding is dependent on environmental conditions and correlates with the induction of gene expression in response to stress. Taken together, these results underscore the importance of genome structure in understanding the regulatory dynamics of transcriptional networks.

Two distinct modes of nucleosome modulation associated with different degrees of dependence of nucleosome positioning on the underlying DNA sequence

BACKGROUND

The nucleosome is the fundamental unit of eukaryotic genomes. Its positioning plays a central role in diverse cellular processes that rely on access to genomic DNA. Experimental evidence suggests that the genomic DNA sequence is one important determinant of nucleosome positioning. Yet it is less clear whether the role of the underlying DNA sequence in nucleosome positioning varies across different promoters. Whether different determinants of nucleosome positioning have characteristic influences on nucleosome modulation also remains to be elucidated.

RESULTS

We identified two typical promoter classes in yeast associated with high or low dependence of nucleosome positioning on the underlying DNA sequence, respectively. Importantly, the two classes have low or high intrinsic sequence preferences for nucleosomes, respectively. The two classes are further distinguished by multiple promoter features, including nucleosome occupancy, nucleosome fuzziness, H2A.Z occupancy, changes in nucleosome positions before and after transcriptional perturbation, and gene activity. Both classes have significantly high turnover rates of histone H3, but employ distinct modes of nucleosome modulation: The first class is characterized by hyperacetylation, whereas the second class is highly regulated by ATP-dependent chromatin remodelling.

CONCLUSION

Our results, coupled with the known features of nucleosome modulation, suggest that the two distinct modes of nucleosome modulation selectively employed by different genes are linked with the intrinsic sequence preferences for nucleosomes. The difference in modes of nucleosome modulation can account for the difference in the contribution of DNA sequence to nucleosome positioning between both promoter classes.

Cell-type selective chromatin remodeling defines the active subset of FOXA1-bound enhancers

Selective activity of a specific set of enhancers defines tissue-specific gene transcription. The pioneer factor FOXA1 has been shown to induce functional enhancer competency through chromatin openings. We have previously found that FOXA1 is recruited to thousands of regions across the genome of a given cell type. Here, we monitored the chromatin structure at FOXA1 binding sites on a chromosome-wide scale using formaldehyde assisted isolation of regulatory elements (FAIRE). Surprisingly, we find that a significant fraction of FOXA1-bound sites have a relatively closed chromatin conformation linked to a shift of the epigenetic signature toward repressive histone marks. Importantly, these sites are not correlated with gene expression in a given cell type suggesting that FOXA1 is required, but not sufficient, for the functional activity of bound enhancers. Interestingly, we find that a significant proportion of the inactive FOXA1-bound regulatory sites in one cell type are actually functional in another cellular context. We found that at least half of the FOXA1 binding sites from a given cell type are shared with another cell lineage. Mechanisms that restrict the activity of shared FOXA1-bound enhancers likely play a significant role in defining the cell-type-specific functions of FOXA1.

Transcription factor control of growth rate dependent genes in Saccharomyces cerevisiae: a three factor design

Background

Characterization of cellular growth is central to understanding living systems. Here, we applied a three-factor design to study the relationship between specific growth rate and genome-wide gene expression in 36 steady-state chemostat cultures of Saccharomyces cerevisiae. The three factors we considered were specific growth rate, nutrient limitation, and oxygen availability.

Results

We identified 268 growth rate dependent genes, independent of nutrient limitation and oxygen availability. The transcriptional response was used to identify key areas in metabolism around which mRNA expression changes are significantly associated. Among key metabolic pathways, this analysis revealed de novo synthesis of pyrimidine ribonucleotides and ATP producing and consuming reactions at fast cellular growth. By scoring the significance of overlap between growth rate dependent genes and known transcription factor target sets, transcription factors that coordinate balanced growth were also identified. Our analysis shows that Fhl1, Rap1, and Sfp1, regulating protein biosynthesis, have significantly enriched target sets for genes up-regulated with increasing growth rate. Cell cycle regulators, such as Ace2 and Swi6, and stress response regulators, such as Yap1, were also shown to have significantly enriched target sets.

Conclusion

Our work, which is the first genome-wide gene expression study to investigate specific growth rate and consider the impact of oxygen availability, provides a more conservative estimate of growth rate dependent genes than previously reported. We also provide a global view of how a small set of transcription factors, 13 in total, contribute to control of cellular growth rate. We anticipate that multi-factorial designs will play an increasing role in elucidating cellular regulation.

The nuclear actin-related protein Act3p/Arp4 influences yeast cell shape and bulk chromatin organization

ACT3/ARP4 is an essential gene, coding for the actin-related protein Act3p/Arp4 of Saccharomyces cerevisiae located within the nucleus. Act3p/Arp4 is a stoichiometric component of the NuA4, INO80, and SWR1 chromatin modulating complexes, and recruits these complexes onto chromatin for their proper chromatin functions. Mutated Act3p/Arp4 leads to impairment of the functions of these complexes and affects transcription of specific genes. Our results revealed significant disorder in the cell size and shape of act3/arp4 mutant cells, when grown at permissive temperature. act3/arp4 mutants have also demonstrated an increase in their nuclear diameters, thus suggesting that Act3p/Arp4 is a key regulator in the maintenance of cellular shape and nuclear organization. Furthermore, the use of Chromatin Yeast Comet Assay (ChYCA) for assessment of single-cell bulk chromatin organization in act3/arp4 mutant cells allowed us to detect an elevated sensitivity toward nuclease action, denoting differences in higher-order chromatin structure of the mutants.

Histone modifications and chromatin dynamics: a focus on filamentous fungi

The readout of the genetic information of eukaryotic organisms is significantly regulated by modifications of DNA and chromatin proteins. Chromatin alterations induce genome-wide and local changes in gene expression and affect a variety of processes in response to internal and external signals during growth, differentiation, development, in metabolic processes, diseases, and abiotic and biotic stresses. This review aims at summarizing the roles of histone H1 and the acetylation and methylation of histones in filamentous fungi and links this knowledge to the huge body of data from other systems. Filamentous fungi show a wide range of morphologies and have developed a complex network of genes that enables them to use a great variety of substrates. This fact, together with the possibility of simple and quick genetic manipulation, highlights these organisms as model systems for the investigation of gene regulation. However, little is still known about regulation at the chromatin level in filamentous fungi. Understanding the role of chromatin in transcriptional regulation would be of utmost importance with respect to the impact of filamentous fungi in human diseases and agriculture. The synthesis of compounds (antibiotics, immunosuppressants, toxins, and compounds with adverse effects) is also likely to be regulated at the chromatin level.