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.

Two-way feedback between chromatin compaction and histone modification state explains S. 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. Here, we propose a broader view: chromatin compaction can be both a cause and a consequence of the histone modification state, and this tight bidirectional interaction can underpin bistable transcriptional states. To test this theory, we developed a mathematical model for the dynamics of the HMR locus in S. 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 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. This bidirectional feedback between chromatin compaction and the histone modification state may be an important regulatory mechanism at many loci.

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 Rap1, and Sir proteins. Despite a vast 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 bi-directional 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 transcriptional activator Rap1 in Sir-silenced chromatin. Using a highly sensitive assay that monitors loss-of-silencing events, we identified a novel 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. Significance Statement: The coarse partitioning of the genome into regions of active euchromatin and repressed heterochromatin is an important, and conserved, level gene expression regulation in eukaryotes. Repressor Activator Protein (Rap1) is a transcription factor that promotes the activation of genes when recruited to promoters, and aids in the establishment of heterochromatin through interactions with silencer elements. Here, we investigate the role of Rap1 when bound to a promoter in silent chromatin and dissect the context-specific epigenetic cues that regulate the dual properties of this transcription factor. Together, our data highlight the importance of protein-protein interactions and local chromatin state on transcription factor function.

Limits to transcriptional silencing in Saccharomyces cerevisiae.

Mating-type switching in the budding yeast Saccharomyces cerevisiae relies on the Sir protein complex to silence HML and HMR, the two loci containing copies of the alleles of the mating type locus, MAT. Sir-based transcriptional silencing has been considered locus-specific, but the recent discovery of rare and transient escapes from silencing at HMLα2 with a sensitive assay called to question if these events extend to the whole locus. Adapting the same assay, we measured that transient silencing failures at HML were more frequent for the α2 gene than α1, similarly to their expression level in unsilenced cells. By coupling a mating assay, at HML we found that one of the two genes at that locus can be transiently expressed while the other gene is maintained silent. Thus, transient silencing loss can be a property of the gene rather than the locus. Cells lacking the SIR1 gene experience epigenetic bistability at HML and HMR. Our previous result led us to ask if HML could allow for two independent epigenetic states within the locus in a sir1Δ mutant. A simple construct using a double fluorescent reporter at HMLα1 and HMLα2 ruled out this possibility. Each HML locus displayed a single epigenetic state. We revisited the question of the correlation between the states of two HML loci in diploid cells, and showed they were independent. Finally, we determined the relative strength of gene repression achieved by Sir-based silencing with that achieved by the a1-α2 repressor.

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.

Distinguishing between recruitment and spread of silent chromatin structures in Saccharomyces cerevisiae.

The formation of heterochromatin at HML, HMR, and telomeres in Saccharomyces cerevisiae involves two main steps: the recruitment of Sir proteins to silencers and their spread throughout the silenced domain. We developed a method to study these two processes at single basepair resolution. Using a fusion protein between the heterochromatin protein Sir3 and the nonsite-specific bacterial adenine methyltransferase M.EcoGII, we mapped sites of Sir3-chromatin interactions genome-wide using long-read Nanopore sequencing to detect adenines methylated by the fusion protein and by ChIP-seq to map the distribution of Sir3-M.EcoGII. A silencing-deficient mutant of Sir3 lacking its Bromo-Adjacent Homology (BAH) domain, sir3-bah∆, was still recruited to HML, HMR, and telomeres. However, in the absence of the BAH domain, it was unable to spread away from those recruitment sites. Overexpression of Sir3 did not lead to further spreading at HML, HMR, and most telomeres. A few exceptional telomeres, like 6R, exhibited a small amount of Sir3 spreading, suggesting that boundaries at telomeres responded variably to Sir3-M.EcoGII overexpression. Finally, by using a temperature-sensitive allele of SIR3 fused to M.ECOGII, we tracked the positions first methylated after induction and found that repression of genes at HML and HMR began before Sir3 occupied the entire locus.

A novel allele of SIR2 reveals a heritable intermediate state of gene silencing.

Genetic information acquires additional meaning through epigenetic regulation, the process by which genetically identical cells can exhibit heritable differences in gene expression and phenotype. Inheritance of epigenetic information is a critical step in maintaining cellular identity and organismal health. In Saccharomyces cerevisiae, one form of epigenetic regulation is the transcriptional silencing of two mating-type loci, HML and HMR, by the SIR-protein complex. To focus on the epigenetic dimension of this gene regulation, we conducted a forward mutagenesis screen to identify mutants exhibiting an epigenetic or metastable silencing defect. We utilized fluorescent reporters at HML and HMR, and screened yeast colonies for epigenetic silencing defects. We uncovered numerous independent sir1 alleles, a gene known to be required for stable epigenetic inheritance. More interestingly, we recovered a missense mutation within SIR2, which encodes a highly conserved histone deacetylase. In contrast to sir1Δ, which exhibits states that are either fully silenced or fully expressed, this sir2 allele exhibited heritable states that were either fully silenced or expressed at an intermediate level. The heritable nature of this unique silencing defect was influenced by, but not completely dependent on, changes in rDNA copy number. Therefore, this study revealed a heritable state of intermediate silencing and linked this state to a central silencing factor, Sir2.

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.

S-phase-independent silencing establishment in Saccharomyces cerevisiae.

The establishment of silent chromatin, a heterochromatin-like structure at HML and HMR in Saccharomyces cerevisiae, depends on progression through S phase of the cell cycle, but the molecular nature of this requirement has remained elusive despite intensive study. Using high-resolution chromatin immunoprecipitation and single-molecule RNA analysis, we found that silencing establishment proceeded via gradual repression of transcription in individual cells over several cell cycles, and that the cell-cycle-regulated step was downstream of Sir protein recruitment. In contrast to prior results, HML and HMR had identical cell-cycle requirements for silencing establishment, with no apparent contribution from a tRNA gene adjacent to HMR. We identified the cause of the S-phase requirement for silencing establishment: removal of transcription-favoring histone modifications deposited by Dot1, Sas2, and Rtt109. These results revealed that silencing establishment was absolutely dependent on the cell-cycle-regulated interplay between euchromatic and heterochromatic histone modifications.

A proactive genotype-to-patient-phenotype map for cystathionine beta-synthase.

BACKGROUND: For the majority of rare clinical missense variants, pathogenicity status cannot currently be classified. Classical homocystinuria, characterized by elevated homocysteine in plasma and urine, is caused by variants in the cystathionine beta-synthase (CBS) gene, most of which are rare. With early detection, existing therapies are highly effective. METHODS: Damaging CBS variants can be detected based on their failure to restore growth in yeast cells lacking the yeast ortholog CYS4. This assay has only been applied reactively, after first observing a variant in patients. Using saturation codon-mutagenesis, en masse growth selection, and sequencing, we generated a comprehensive, proactive map of CBS missense variant function. RESULTS: Our CBS variant effect map far exceeds the performance of computational predictors of disease variants. Map scores correlated strongly with both disease severity (Spearman's ϱ = 0.9) and human clinical response to vitamin B6 (ϱ = 0.93). CONCLUSIONS: We demonstrate that highly multiplexed cell-based assays can yield proactive maps of variant function and patient response to therapy, even for rare variants not previously seen in the clinic.

Epigenetic memory independent of symmetric histone inheritance.

Heterochromatic gene silencing is an important form of gene regulation that usually requires specific histone modifications. A popular model posits that inheritance of modified histones, especially in the form of H3-H4 tetramers, underlies inheritance of heterochromatin. Because H3-H4 tetramers are randomly distributed between daughter chromatids during DNA replication, rare occurrences of asymmetric tetramer inheritance within a heterochromatic domain would have the potential to destabilize heterochromatin. This model makes a prediction that shorter heterochromatic domains would experience unbalanced tetramer inheritance more frequently, and thereby be less stable. In contrast to this prediction, we found that shortening a heterochromatic domain in Saccharomyces had no impact on the strength of silencing nor its heritability. Additionally, we found that replisome mutations that disrupt inheritance of H3-H4 tetramers had only minor effects on heterochromatin stability. These findings suggest that histones carry little or no memory of the heterochromatin state through DNA replication.

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.

Mutations in the PCNA DNA Polymerase Clamp of Saccharomyces cerevisiae Reveal Complexities of the Cell Cycle and Ploidy on Heterochromatin Assembly.

In Saccharomyces cerevisiae, transcriptional silencing at HML and HMR maintains mating-type identity. The repressive chromatin structure at these loci is replicated every cell cycle and must be re-established quickly to prevent transcription of the genes at these loci. Mutations in a component of the replisome, the proliferating cell nuclear antigen (PCNA), encoded by POL30, cause a loss of transcriptional silencing at HMR We used an assay that captures transient losses of silencing at HML and HMR to perform extended genetic analyses of the pol30-6, pol30-8, and pol30-79 alleles. All three alleles destabilized silencing only transiently and only in cycling cells. Whereas pol30-8 caused loss of silencing by disrupting the function of Chromatin Assembly Factor 1, pol30-6 and pol30-79 acted through a separate genetic pathway, but one still dependent on histone chaperones. Surprisingly, the silencing-loss phenotypes of pol30-6 and pol30-79 depended on ploidy, but not on POL30 dosage or mating-type identity. Separately from silencing loss, the pol30-6 and pol30-79 alleles also displayed high levels of mitotic recombination in diploids. These results established that histone trafficking involving PCNA at replication forks is crucial to the maintenance of chromatin state and genome stability during DNA replication. They also raised the possibility that increased ploidy may protect chromatin states when the replisome is perturbed.

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.

Epigenomic profiling of newborns with isolated orofacial clefts reveals widespread DNA methylation changes and implicates metastable epiallele regions in disease risk.

Cleft lip with or without cleft palate (CL/P) is a common human birth defect whose etiologies remain largely unknown. Several studies have demonstrated that periconceptional supplementation of folic acid can reduce risk of CL/P in offspring. In this study, we tested the hypothesis that the preventive effect of folic acid is manifested through epigenetic modifications by determining whether DNA methylation changes are associated with CL/P. To more readily observe the potential effects of maternal folate on the offspring epigenome, we focused on births prior to mandatory dietary folate fortification in the United States (i.e. birth year 1997 or earlier). Genomic DNA methylation levels were assessed from archived newborn bloodspots in a 182-member case-control study using the Illumina® Human Beadchip 450K array. CL/P cases displayed striking epigenome-wide hypomethylation relative to controls: 63% of CpGs interrogated had lower methylation levels in case newborns, a trend which held up in racially stratified sub-groups. 28 CpG sites reached epigenome-wide significance and all were case-hypomethylated. The most significant CL/P-associated differentially methylated region encompassed the VTRNA2-1 gene, which was also hypomethylated in cases (FWER p = 0.014). This region has been previously characterized as a nutritionally-responsive, metastable epiallele and CL/P-associated methylation changes, in general, were greater at or near putative metastable epiallelic regions. Gene Set Enrichment Analysis of CL/P-associated DMRs showed an over-representation of genes involved in palate development such as WNT9B, MIR140 and LHX8. CL/P-associated DNA methylation changes may partly explain the mechanism by which orofacial clefts are responsive to maternal folate levels.

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.

Impact of Homologous Recombination on Silent Chromatin in Saccharomyces cerevisiae.

Specialized chromatin domains repress transcription of genes within them and present a barrier to many DNA-protein interactions. Silent chromatin in the budding yeast Saccharomyces cerevisiae, akin to heterochromatin of metazoans and plants, inhibits transcription of PolII- and PolIII-transcribed genes, yet somehow grants access to proteins necessary for DNA transactions like replication and homologous recombination. In this study, we adapted a novel assay to detect even transient changes in the dynamics of transcriptional silencing at HML after it served as a template for homologous recombination. Homologous recombination specifically targeted to HML via double-strand-break formation at a homologous locus often led to transient loss of transcriptional silencing at HML Interestingly, many cells could template homology-directed repair at HML without an obligate loss of silencing, even in recombination events with extensive gene conversion tracts. In a population of cells that experienced silencing loss following recombination, transcription persisted for 2-3 hr after all double-strand breaks were repaired. mRNA levels from cells that experienced recombination-induced silencing loss did not approach the amount of mRNA seen in cells lacking transcriptional silencing. Thus, silencing loss at HML after homologous recombination was short-lived and limited.

DNA Methylome Profiling on the Infinium HumanMethylation450 Array from Limiting Quantities of Genomic DNA from a Single, Small Archived Bloodspot.

AIMS: Archived newborn bloodspots are valuable sample collections for genetic and epigenetic disease research. However, they have often been stored for long periods of time, under less than ideal circumstances, and nucleic acid yields can be low, particularly when samples become limiting. We wished to determine whether the quantity and quality of genomic DNA (gDNA) isolated from a single, surgical bloodspot punch (2 mm dia.) was adequate for accurate and reliable DNA methylome profiling on the Illumina HumanMethylation450 array. METHODS: A total of 25-750 ng of archived bloodspot or Jurkat cell gDNA were bisulfite converted and analyzed on the array without any additional DNA amplification steps. Methylation profiles were assessed for call rate, call confidence (detection p-value), and reproducibility. RESULTS: Using 25 ng gDNA from either Jurkat cells or dried bloodspots, array-wide call rates (∼99.9%) and detection p-values (99.9% with p < 5 × 10-6) were excellent. There was good agreement between methylation profiles generated from 25 ng gDNA and those generated from 750 ng (ρ > 0.98), although a fraction of CpG sites (2-8% depending on experiment) exhibited quantitative differences. Genome-wide methylation levels were strikingly reproducible from 25 ng DNA in both replicate and interindividual samples (ρ > 0.98). CONCLUSIONS: Twenty-five nanograms of gDNA, isolated from a single, surgical punch (2 mm dia.) of an archived newborn bloodspot, generate a genome-wide methylation profile on the Illumina HumanMethylation450 array that is robust, reproducible, and suitable for differential methylation studies.