Distinct silencer states generate epigenetic states of heterochromatin.

Heterochromatic loci can exhibit different transcriptional states in genetically identical cells. A popular model posits that the inheritance of modified histones is sufficient for inheritance of the silenced state. However, silencing inheritance requires silencers and therefore cannot be driven by the inheritance of modified histones alone. To address these observations, we determined the chromatin architectures produced by strong and weak silencers in Saccharomyces. Strong silencers recruited Sir proteins and silenced the locus in all cells. Strikingly, weakening these silencers reduced Sir protein recruitment and stably silenced the locus in some cells; however, this silenced state could probabilistically convert to an expressed state that lacked Sir protein recruitment. Additionally, changes in the constellation of silencer-bound proteins or the concentration of a structural Sir protein modulated the probability that a locus exhibited the silenced or expressed state. These findings argued that distinct silencer states generate epigenetic states and regulate their dynamics.

Metabolism and epigenetics.

Epigenetic mechanisms by which cells inherit information are, to a large extent, enabled by DNA methylation and posttranslational modifications of histone proteins. These modifications operate both to influence the structure of chromatin per se and to serve as recognition elements for proteins with motifs dedicated to binding particular modifications. Each of these modifications results from an enzyme that consumes one of several important metabolites during catalysis. Likewise, the removal of these marks often results in the consumption of a different metabolite. Therefore, these so-called epigenetic marks have the capacity to integrate the expression state of chromatin with the metabolic state of the cell. This review focuses on the central roles played by acetyl-CoA, S-adenosyl methionine, NAD(+), and a growing list of other acyl-CoA derivatives in epigenetic processes. We also review how metabolites that accumulate as a result of oncogenic mutations are thought to subvert the epigenetic program.

Two-way feedback between chromatin compaction and histone modification state explains Saccharomyces cerevisiae heterochromatin bistability.

Compact chromatin is closely linked with gene silencing in part by sterically masking access to promoters, inhibiting transcription factor binding and preventing polymerase from efficiently transcribing a gene. However, a broader hypothesis suggests that chromatin compaction can be both a cause and a consequence of the locus histone modification state, with a tight bidirectional interaction underpinning bistable transcriptional states. To rigorously test this hypothesis, we developed a mathematical model for the dynamics of the HMR locus in Saccharomyces cerevisiae, that incorporates activating histone modifications, silencing proteins, and a dynamic, acetylation-dependent, three-dimensional locus size. Chromatin compaction enhances silencer protein binding, which in turn feeds back to remove activating histone modifications, leading to further compaction. The bistable output of the model was in good agreement with prior quantitative data, including switching rates from expressed to silent states (and vice versa), and protein binding/histone modification levels within the locus. We then tested the model by predicting changes in switching rates as the genetic length of the locus was increased, which were then experimentally verified. Such bidirectional feedback between chromatin compaction and the histone modification state may be a widespread and important regulatory mechanism given the hallmarks of many heterochromatic regions: physical chromatin compaction and dimerizing (or multivalent) silencing proteins.

Context-dependent function of the transcriptional regulator Rap1 in gene silencing and activation in Saccharomyces cerevisiae.

In Saccharomyces cerevisiae, heterochromatin is formed through interactions between site-specific DNA-binding factors, including the transcriptional activator Repressor Activator Protein (Rap1), and Sir proteins. Despite an understanding of the establishment and maintenance of Sir-silenced chromatin, the mechanism of gene silencing by Sir proteins has remained a mystery. Utilizing high-resolution chromatin immunoprecipitation, we found that Rap1, the native activator of the bidirectional HMLα promoter, bound its recognition sequence in silenced chromatin, and its binding was enhanced by the presence of Sir proteins. In contrast to prior results, various components of transcription machinery were not able to access HMLα in the silenced state. These findings disproved the long-standing model of indiscriminate steric occlusion by Sir proteins and led to investigation of the role of the transcriptional activator Rap1 in Sir-silenced chromatin. Using a highly sensitive assay that monitors loss-of-silencing events, we identified a role for promoter-bound Rap1 in the maintenance of silent chromatin through interactions with the Sir complex. We also found that promoter-bound Rap1 activated HMLα when in an expressed state, and aided in the transition from transcription initiation to elongation. Highlighting the importance of epigenetic context in transcription factor function, these results point toward a model in which the duality of Rap1 function was mediated by local chromatin environment rather than binding-site availability.

Nucleosome Positioning Regulates the Establishment, Stability, and Inheritance of Heterochromatin in Saccharomyces cerevisiae.

Heterochromatic domains are complex structures composed of nucleosome arrays that are bound by silencing factors. This composition raises the possibility that certain configurations of nucleosome arrays facilitate heterochromatic silencing. We tested this possibility in Saccharomyces cerevisiae by systematically altering the distance between heterochromatic nucleosome-depleted regions (NDRs), which is predicted to affect local nucleosome positioning by limiting how nucleosomes can be packed between NDRs. Consistent with this prediction, serial deletions that altered the distance between heterochromatic NDRs revealed a striking oscillatory relationship between inter-NDR distance and defects in nucleosome positioning. Furthermore, conditions that caused poor nucleosome positioning also led to defects in both heterochromatin stability and the ability of cells to generate and inherit epigenetic transcriptional states. These findings strongly suggest that nucleosome positioning can contribute to formation and maintenance of functional heterochromatin and point to previously unappreciated roles of NDR positioning within heterochromatic domains.

The nucleosome core particle remembers its position through DNA replication and RNA transcription.

Nucleosomes are the fundamental structural unit of chromatin. In addition to stabilizing the DNA polymer, nucleosomes are modified in ways that reflect and affect gene expression in their vicinity. It has long been assumed that nucleosomes can transmit memory of gene expression through their covalent posttranslational modifications. An unproven assumption of this model, which is essential to most models of epigenetic inheritance, is that a nucleosome present at a locus reoccupies the same locus after DNA replication. We tested this assumption by nucleating a synthetic chromatin domain in vivo, in which ∼4 nucleosomes at an arbitrary locus were covalently labeled with biotin. We tracked the fate of labeled nucleosomes through DNA replication, and established that nucleosomes present at a locus remembered their position during DNA replication. The replication-associated histone chaperones Dpb3 and Mcm2 were essential for nucleosome position memory, and in the absence of both Dpb3 and Mcm2 histone chaperone activity, nucleosomes did not remember their position. Using the same approach, we tested the model that transcription results in retrograde transposition of nucleosomes along a transcription unit. We found no evidence of retrograde transposition. Our results suggest that nucleosomes have the capacity to transmit epigenetic memory across mitotic generations with exquisite spatial fidelity.

The molecular topography of silenced chromatin in Saccharomyces cerevisiae.

Heterochromatin imparts regional, promoter-independent repression of genes and is epigenetically heritable. Understanding how silencing achieves this regional repression is a fundamental problem in genetics and development. Current models of yeast silencing posit that Sir proteins, recruited by transcription factors bound to the silencers, spread throughout the silenced region. To test this model directly at high resolution, we probed the silenced chromatin architecture by chromatin immunoprecipitation (ChIP) followed by next-generation sequencing (ChIP-seq) of Sir proteins, histones, and a key histone modification, H4K16-acetyl. These analyses revealed that Sir proteins are strikingly concentrated at and immediately adjacent to the silencers, with lower levels of enrichment over the promoters at HML and HMR, the critical targets for transcriptional repression. The telomeres also showed discrete peaks of Sir enrichment yet a continuous domain of hypoacetylated histone H4K16. Surprisingly, ChIP-seq of cross-linked chromatin revealed a distribution of nucleosomes at silenced loci that was similar to Sir proteins, whereas native nucleosome maps showed a regular distribution throughout silenced loci, indicating that cross-linking captured a specialized chromatin organization imposed by Sir proteins. This specialized chromatin architecture observed in yeast informs the importance of a steric contribution to regional repression in other organisms.

Pivotal roles of PCNA loading and unloading in heterochromatin function.

In Saccharomyces cerevisiae, heterochromatin structures required for transcriptional silencing of the HML and HMR loci are duplicated in coordination with passing DNA replication forks. Despite major reorganization of chromatin structure, the heterochromatic, transcriptionally silent states of HML and HMR are successfully maintained throughout S-phase. Mutations of specific components of the replisome diminish the capacity to maintain silencing of HML and HMR through replication. Similarly, mutations in histone chaperones involved in replication-coupled nucleosome assembly reduce gene silencing. Bridging these observations, we determined that the proliferating cell nuclear antigen (PCNA) unloading activity of Elg1 was important for coordinating DNA replication forks with the process of replication-coupled nucleosome assembly to maintain silencing of HML and HMR through S-phase. Collectively, these data identified a mechanism by which chromatin reassembly is coordinated with DNA replication to maintain silencing through S-phase.

Assessing computational predictions of the phenotypic effect of cystathionine-beta-synthase variants.

Accurate prediction of the impact of genomic variation on phenotype is a major goal of computational biology and an important contributor to personalized medicine. Computational predictions can lead to a better understanding of the mechanisms underlying genetic diseases, including cancer, but their adoption requires thorough and unbiased assessment. Cystathionine-beta-synthase (CBS) is an enzyme that catalyzes the first step of the transsulfuration pathway, from homocysteine to cystathionine, and in which variations are associated with human hyperhomocysteinemia and homocystinuria. We have created a computational challenge under the CAGI framework to evaluate how well different methods can predict the phenotypic effect(s) of CBS single amino acid substitutions using a blinded experimental data set. CAGI participants were asked to predict yeast growth based on the identity of the mutations. The performance of the methods was evaluated using several metrics. The CBS challenge highlighted the difficulty of predicting the phenotype of an ex vivo system in a model organism when classification models were trained on human disease data. We also discuss the variations in difficulty of prediction for known benign and deleterious variants, as well as identify methodological and experimental constraints with lessons to be learned for future challenges.

Accumulation of rare coding variants in genes implicated in risk of human cleft lip with or without cleft palate.

Cleft lip with/without cleft palate (CLP) is a common craniofacial malformation with complex etiologies, reflecting both genetic and environmental factors. Most of the suspected genetic risk for CLP has yet to be identified. To further classify risk loci and estimate the contribution of rare variants, we sequenced the exons in 49 candidate genes in 323 CLP cases and 211 nonmalformed controls. Our findings indicated that rare, protein-altering variants displayed markedly higher burdens in CLP cases at relevant loci. First, putative loss-of-function mutations (nonsense, frameshift) were significantly enriched among cases: 13 of 323 cases (~4%) harbored such alleles within these 49 genes, versus one such change in controls (p = 0.01). Second, in gene-level analyses, the burden of rare alleles showed greater case-association for several genes previously implicated in cleft risk. For example, BHMT displayed a 10-fold increase in protein-altering variants in CLP cases (p = .03), including multiple case occurrences of a rare frameshift mutation (K400 fs). Other loci with greater rare, coding allele burdens in cases were in signaling pathways relevant to craniofacial development (WNT9B, BMP4, BMPR1B) as well as the methionine cycle (MTRR). We conclude that rare coding variants may confer risk for isolated CLP.

Riches of phenotype computationally extracted from microbial colonies.

The genetic, epigenetic, and physiological differences among cells in clonal microbial colonies are underexplored opportunities for discovery. A recently developed genetic assay reveals that transient losses of heterochromatic repression, a heritable form of gene silencing, occur throughout the growth of Saccharomyces colonies. This assay requires analyzing two-color fluorescence patterns in yeast colonies, which is qualitatively appealing but quantitatively challenging. In this paper, we developed a suite of automated image processing, visualization, and classification algorithms (MORPHE) that facilitated the analysis of heterochromatin dynamics in the context of colonial growth and that can be broadly adapted to many colony-based assays in Saccharomyces and other microbes. Using the features that were automatically extracted from fluorescence images, our classification method distinguished loss-of-silencing patterns between mutants and wild type with unprecedented precision. Application of MORPHE revealed subtle but significant differences in the stability of heterochromatic repression between various environmental conditions, revealed that haploid cells experienced higher rates of silencing loss than diploids, and uncovered the unexpected contribution of a sirtuin to heterochromatin dynamics.

Switching the mechanism of mating type switching: a domesticated transposase supplants a domesticated homing endonuclease.

Programmed DNA rearrangements are critical for the development of many organisms and, intriguingly, can be catalyzed by domesticated mobile genetic elements. In this issue of Genes & Development, Barsoum and colleagues (pp. 33-44) demonstrate that, in the budding yeast Kluyveromyces lactis, a DNA rearrangement associated with mating type switching requires a domesticated transposase and occurs through a mechanism distinct from that in the related yeast, Saccharomyces cerevisiae. Thus, mechanisms for mating type switching have evolved multiple times, indicating the relative ease with which mobile genetic elements can be captured.

Sequence variation in folate pathway genes and risks of human cleft lip with or without cleft palate.

In an effort to comprehensively interrogate genetic variation in the folate pathway for risk of cleft lip with or without cleft palate (CLP), we evaluated 504 common and rare variants in 35 folate-related genes in a panel of 330 infants with CLP and 367 non-malformed controls. Odds ratios (OR) with 95% confidence intervals were computed for common genotypes. A Case-Control Difference metric was calculated for rare variants to highlight differentially occurring alleles. Interactions between variants and a maternal folate intake variable were also evaluated. In gene-only results, significant odds ratios were observed for multiple variants in the BHMT/BHMT2/DMGDH gene cluster, particularly in Hispanic infants. Also in this cluster, rare variant analysis highlighted a substantial case-control difference in BHMT rs60340837 (synonymous Y284Y). In Hispanics, the ALDH1L1 I812V variant (rs4646750) was the most significant risk allele: OR = 3.8 (95%CI = 1.6-9.2) when heterozygous. In non-Hispanic white infants, we observed significant risk for AHCYL2 rs1095423 (homozygous OR = 3.0, 95%CI 1.1-7.8) and the 68 bp CBS insertion (c.844ins68; heterozygous OR = 2.4, 95%CI = 1.1-5.3). Rare variant analysis in this group revealed case-control differences in MTRR and several other methionine cycle genes, a process implicated previously in clefting risk. In women with low folate intake specifically, increased risks were observed for CBS rs2851391 (OR = 3.6, 95%CI = 1.3-9.6) and the R259P nonsynonymous variant of TCN2 (rs1801198; OR = 2.8, 95%CI = 1.2-6.3). This comprehensive study provides further direction on candidate loci to help disentangle the folate-related developmental phenomena in human clefting risk. © 2016 Wiley Periodicals, Inc.

Highly expressed loci are vulnerable to misleading ChIP localization of multiple unrelated proteins.

Chromatin immunoprecipitation (ChIP) is the gold-standard technique for localizing nuclear proteins in the genome. We used ChIP, in combination with deep sequencing (Seq), to study the genome-wide distribution of the Silent information regulator (Sir) complex in Saccharomyces cerevisiae. We analyzed ChIP-Seq peaks of the Sir2, Sir3, and Sir4 silencing proteins and discovered 238 unexpected euchromatic loci that exhibited enrichment of all three. Surprisingly, published ChIP-Seq datasets for the Ste12 transcription factor and the centromeric Cse4 protein indicated that these proteins were also enriched in the same euchromatic regions with the high Sir protein levels. The 238 loci, termed "hyper-ChIPable", were in highly expressed regions with strong polymerase II and polymerase III enrichment signals, and the correlation between transcription level and ChIP enrichment was not limited to these 238 loci but extended genome-wide. The apparent enrichment of various proteins at hyper-ChIPable loci was not a consequence of artifacts associated with deep sequencing methods, as confirmed by ChIP-quantitative PCR. The localization of unrelated proteins, including the entire silencing complex, to the most highly transcribed genes was highly suggestive of a technical issue with the immunoprecipitations. ChIP-Seq on chromatin immunoprecipitated with a nuclear-localized GFP reproduced the above enrichment in an expression-dependent manner: induction of the GAL genes resulted in an increased ChIP signal of the GFP protein at these loci, with presumably no biological relevance. Whereas ChIP is a broadly valuable technique, some published conclusions based upon ChIP procedures may merit reevaluation in light of these findings.

Two surfaces on the histone chaperone Rtt106 mediate histone binding, replication, and silencing.

The histone chaperone Rtt106 binds histone H3 acetylated at lysine 56 (H3K56ac) and facilitates nucleosome assembly during several molecular processes. Both the structural basis of this modification-specific recognition and how this recognition informs Rtt106 function are presently unclear. Guided by our crystal structure of Rtt106, we identified two regions on its double-pleckstrin homology domain architecture that mediated histone binding. When histone binding was compromised, Rtt106 localized properly to chromatin but failed to deliver H3K56ac, leading to replication and silencing defects. By mutating analogous regions in the structurally homologous chromatin-reorganizer Pob3, we revealed a conserved histone-binding function for a basic patch found on both proteins. In contrast, a loop connecting two ß-strands was required for histone binding by Rtt106 but was dispensable for Pob3 function. Unlike Rtt106, Pob3 histone binding was modification-independent, implicating the loop of Rtt106 in H3K56ac-specific recognition in vivo. Our studies described the structural origins of Rtt106 function, identified a conserved histone-binding surface, and defined a critical role for Rtt106:H3K56ac-binding specificity in silencing and replication-coupled nucleosome turnover.

Direct regulation of nucleosome density by the conserved AAA-ATPase Yta7.

Yta7 is a highly conserved bromodomain-containing protein with AAA-ATPase homology originally implicated in heterochromatin boundary function in Saccharomyces cerevisiae. Although increased activity of the human ortholog has been implicated in malignant breast tumors, Yta7's precise mode of action is unknown. Transcriptional analysis in yeast cells revealed a role for Yta7 and its ATPase function in gene induction, including galactose- and sporulation-induced transcription. This requirement was direct and activating, because Yta7 associated with the GAL gene cluster only upon transcriptional induction. Suggestive of a role in transcriptional elongation, Yta7 localized to the ORFs of highly transcribed genes. Intriguingly, the yta7Δ mutant's transcriptional defects were partially suppressed by decreased dosage of histones H3 and H4. Consistent with this suppression, cells lacking Yta7 exhibited both increased levels of chromatin-incorporated histone H3 and decreased nucleosome spacing. Importantly, this modulation of H3 levels occurred independently of changes in H3 transcript level. Because Yta7 binds histone H3 in vitro, these results suggested a direct role for Yta7 in H3 eviction or degradation. Further, local loss of Yta7 activity at a long inducible gene resulted in accumulation of H3 at the 3' end upon transcriptional activation, implying Yta7 may regulate H3 cotranscriptionally.

Symmetry, asymmetry, and kinetics of silencing establishment in Saccharomyces cerevisiae revealed by single-cell optical assays.

In Saccharomyces cerevisiae, silent chromatin inhibits the expression of genes at the HML, HMR, and telomeric loci. When silent chromatin forms de novo, the rate of its establishment is influenced by different chromatin states. In particular, loss of the enzyme Dot1, an H3 K79 methyltransferase, leads to rapid silencing establishment. We tested whether silencing establishment was antagonized by H3 K79 methylation or by the Dot1 protein itself competing with Sir3 for binding sites on nucleosomes. To do so, we monitored fluorescence activity in cells containing a GFP gene within the HML locus during silencing establishment in a series of dot1 and histone mutant backgrounds. Silencing establishment rate was correlated with Dot1's enzymatic function rather than with the Dot1 protein itself. In addition, histone mutants that mimicked the conformation of unmethylated H3 K79 increased the rate of silencing establishment, indicating that the H3 K79 residue affected silencing independently of Dot1 abundance. Using fluorophore-based reporters, we confirmed that mother and daughter cells often silence in concert, but in instances where asymmetric silencing occurs, daughter cells established silencing earlier than their mothers. This noninvasive technique enabled us to demonstrate an asymmetry in silencing establishment of a key regulatory locus controlling cell fate.

Co-evolution of transcriptional silencing proteins and the DNA elements specifying their assembly.

Co-evolution of transcriptional regulatory proteins and their sites of action has been often hypothesized but rarely demonstrated. Here we provide experimental evidence of such co-evolution in yeast silent chromatin, a finding that emerged from studies of hybrids formed between two closely related Saccharomyces species. A unidirectional silencing incompatibility between S. cerevisiae and S. bayanus led to a key discovery: asymmetrical complementation of divergent orthologs of the silent chromatin component Sir4. In S. cerevisiae/S. bayanus interspecies hybrids, ChIP-Seq analysis revealed a restriction against S. cerevisiae Sir4 associating with most S. bayanus silenced regions; in contrast, S. bayanus Sir4 associated with S. cerevisiae silenced loci to an even greater degree than did S. cerevisiae's own Sir4. Functional changes in silencer sequences paralleled changes in Sir4 sequence and a reduction in Sir1 family members in S. cerevisiae. Critically, species-specific silencing of the S. bayanus HMR locus could be reconstituted in S. cerevisiae by co-transfer of the S. bayanus Sir4 and Kos3 (the ancestral relative of Sir1) proteins. As Sir1/Kos3 and Sir4 bind conserved silencer-binding proteins, but not specific DNA sequences, these rapidly evolving proteins served to interpret differences in the two species' silencers presumably involving emergent features created by the regulatory proteins that bind sequences within silencers. The results presented here, and in particular the high resolution ChIP-Seq localization of the Sir4 protein, provided unanticipated insights into the mechanism of silent chromatin assembly in yeast.

Roles for H2A.Z and its acetylation in GAL1 transcription and gene induction, but not GAL1-transcriptional memory.

H2A.Z is a histone H2A variant conserved from yeast to humans, and is found at 63% of promoters in Saccharomyces cerevisiae. This pattern of localization suggests that H2A.Z is somehow important for gene expression or regulation. H2A.Z can be acetylated at up to four lysine residues on its amino-terminal tail, and acetylated-H2A.Z is enriched in chromatin containing promoters of active genes. We investigated whether H2A.Z's role in GAL1 gene regulation and gene expression depends on H2A.Z acetylation. Our findings suggested that H2A.Z functioned both in gene regulation and in gene expression and that only its role in gene regulation depended upon its acetylation. Our findings provided an alternate explanation for results that were previously interpreted as evidence that H2A.Z plays a role in GAL1 transcriptional memory. Additionally, our findings provided new insights into the phenotypes of htz1Delta mutants: in the absence of H2A.Z, the SWR1 complex, which deposits H2A.Z into chromatin, was deleterious to the cell, and many of the phenotypes of cells lacking H2A.Z were due to the SWR1 complex's activity rather than to the absence of H2A.Z per se. These results highlight the need to reevaluate all studies on the phenotypes of cells lacking H2A.Z.

Genetic control of immune cell types in fungal disease.

Millions of people harbor latent infections of the fungus Histoplasma capsulatum. Such persistent infections represent a stalemate between mechanisms of virulence and the immune response. The differing responses of inbred mouse strains to the same pathogen reflect variation in the genes that control the outcome of infection. Here we show that a 250-fold difference in H. capsulatum susceptibility between inbred mouse strains is attributable to the genotype at the MHC H2 locus. Gene expression analysis of strains varying only at the H2 locus identified genotype-specific and genotype-independent expression signatures, including infection-induced genes such as the fungal pattern recognition receptor Clec7a. Surprisingly, B-cell-specific gene expression was negatively correlated with fungal burden, whereas neutrophil-specific genes were correlated with superior disease outcome. Indeed, disease outcome improved when B cells were eliminated and neutrophils were more active, a previously unknown aspect of the host response. These data refine the understanding of genetic influences on histoplasmosis, reveal how shifts in the composition of immune cell populations compel different disease outcomes, and uncover how innate immunity modulation alters histoplasmosis.