Genome-wide mapping of histone modifications and mass spectrometry reveal H4 acetylation bias and H3K36 methylation at gene promoters in fission yeast

KDM8 epigenetically controls cardiac metabolism to prevent initiation of dilated cardiomyopathy

Cardiac metabolism is deranged in heart failure, but underlying mechanisms remain unclear. Here, we show that lysine demethylase 8 (Kdm8) maintains an active mitochondrial gene network by repressing Tbx15, thus preventing dilated cardiomyopathy leading to lethal heart failure. Deletion of Kdm8 in mouse cardiomyocytes increased H3K36me2 with activation of Tbx15 and repression of target genes in the NAD+ pathway before dilated cardiomyopathy initiated. NAD+ supplementation prevented dilated cardiomyopathy in Kdm8 mutant mice, and TBX15 overexpression blunted NAD+-activated cardiomyocyte respiration. Furthermore, KDM8 was downregulated in human hearts affected by dilated cardiomyopathy, and higher TBX15 expression defines a subgroup of affected hearts with the strongest downregulation of genes encoding mitochondrial proteins. Thus, KDM8 represses TBX15 to maintain cardiac metabolism. Our results suggest that epigenetic dysregulation of metabolic gene networks initiates myocardium deterioration toward heart failure and could underlie heterogeneity of dilated cardiomyopathy.

H3K36 methylation reprograms gene expression to drive early gametocyte development in Plasmodium falciparum

Background

The Plasmodium sexual gametocyte stages are the only transmissible form of the malaria parasite and are thus responsible for the continued transmission of the disease. Gametocytes undergo extensive functional and morphological changes from commitment to maturity, directed by an equally extensive control program. However, the processes that drive the differentiation and development of the gametocyte post-commitment, remain largely unexplored. A previous study reported enrichment of H3K36 di- and tri-methylated (H3K36me2&3) histones in early-stage gametocytes. Using chromatin immunoprecipitation followed by high-throughput sequencing, we identify a stage-specific association between these repressive histone modifications and transcriptional reprogramming that define a stage II gametocyte transition point.

Results

Here, we show that H3K36me2 and H3K36me3 from stage II gametocytes are associated with repression of genes involved in asexual proliferation and sexual commitment, indicating that H3K36me2&3-mediated repression of such genes is essential to the transition from early gametocyte differentiation to intermediate development. Importantly, we show that the gene encoding the transcription factor AP2-G as commitment master regulator is enriched with H3K36me2&3 and actively repressed in stage II gametocytes, providing the first evidence of ap2-g gene repression in post-commitment gametocytes. Lastly, we associate the enhanced potency of the pan-selective Jumonji inhibitor JIB-04 in gametocytes with the inhibition of histone demethylation including H3K36me2&3 and a disruption of normal transcriptional programs.

Conclusions

Taken together, our results provide the first description of an association between global gene expression reprogramming and histone post-translational modifications during P. falciparum early sexual development. The stage II gametocyte-specific abundance of H3K36me2&3 manifests predominantly as an independent regulatory mechanism targeted towards genes that are repressed post-commitment. H3K36me2&3-associated repression of genes is therefore involved in key transcriptional shifts that accompany the transition from early gametocyte differentiation to intermediate development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00393-9.

Surprising phenotypic diversity of cancer-associated mutations of Gly 34 in the histone H3 tail

Sequencing of cancer genomes has identified recurrent somatic mutations in histones, termed oncohistones, which are frequently poorly understood. Previously we showed that fission yeast expressing only the H3.3G34R mutant identified in aggressive pediatric glioma had reduced H3K36 trimethylation and acetylation, increased genomic instability and replicative stress, and defective homology-dependent DNA damage repair. Here we show that surprisingly distinct phenotypes result from G34V (also in glioma) and G34W (giant cell tumors of bone) mutations, differentially affecting H3K36 modifications, subtelomeric silencing, genomic stability; sensitivity to irradiation, alkylating agents, and hydroxyurea; and influencing DNA repair. In cancer, only 1 of 30 alleles encoding H3 is mutated. Whilst co-expression of wild-type H3 rescues most G34 mutant phenotypes, G34R causes dominant hydroxyurea sensitivity, homologous recombination defects, and dominant subtelomeric silencing. Together, these studies demonstrate the complexity associated with different substitutions at even a single residue in H3 and highlight the utility of genetically tractable systems for their analysis.

The inhibition of checkpoint activation by telomeres does not involve exclusion of dimethylation of histone H4 lysine 20 (H4K20me2)

DNA double-strand (DSBs) breaks activate the DNA damage checkpoint machinery to pause or halt the cell cycle.  Telomeres, the specific DNA-protein complexes at linear eukaryotic chromosome ends, are capped DSBs that do not activate DNA damage checkpoints.  This "checkpoint privileged" status of telomeres was previously investigated in the yeast Schizosaccharomyces pombe lacking the major double-stranded telomere DNA binding protein Taz1.  Telomeric DNA repeats in cells lacking Taz1 are 10 times longer than normal and contain single-stranded DNA regions.  DNA damage checkpoint proteins associate with these damaged telomeres, but the DNA damage checkpoint is not activated.  This severing of the DNA damage checkpoint signaling pathway was reported to stem from exclusion of histone H4 lysine 20 dimethylation (H4K20me2) from telomeric nucleosomes in both wild type cells and cells lacking Taz1.  However, experiments to identify the mechanism of this exclusion failed, prompting our re-evaluation of H4K20me2 levels at telomeric chromatin.  In this short report, we used an extensive series of controls to identify an antibody specific for the H4K20me2 modification and show that the level of this modification is the same at telomeres and internal loci in both wild type cells and those lacking Taz1.  Consequently, telomeres must block activation of the DNA Damage Response by another mechanism that remains to be determined.

The inhibition of checkpoint activation by telomeres does not involve exclusion of dimethylation of histone H4 lysine 20 (H4K20me2)

DNA double-strand breaks (DSBs) activate the DNA damage checkpoint machinery to pause or halt the cell cycle.  Telomeres, the specific DNA-protein complexes at linear eukaryotic chromosome ends, are capped DSBs that do not activate DNA damage checkpoints.  This "checkpoint privileged" status of telomeres was previously investigated in the yeast  Schizosaccharomyces pombelacking the major double-stranded telomere DNA binding protein Taz1. Telomeric DNA repeats in cells lacking Taz1 are 10 times longer than normal and contain single-stranded DNA regions. DNA damage checkpoint proteins associate with these damaged telomeres, but the DNA damage checkpoint is not activated. This severing of the DNA damage checkpoint signaling pathway was reported to stem from exclusion of histone H4 lysine 20 dimethylation (H4K20me2) from telomeric nucleosomes in both wild type cells and cells lacking Taz1.  However, experiments to identify the mechanism of this exclusion failed, prompting our re-evaluation of H4K20me2 levels at telomeric chromatin.  In this short report, we used an extensive series of controls to identify an antibody specific for the H4K20me2 modification and show that the level of this modification is the same at telomeres and internal loci in both wild type cells and those lacking Taz1.  Consequently, telomeres must block activation of the DNA Damage Response by another mechanism that remains to be determined.

Transcriptional reprogramming in cellular quiescence

Most cells in nature are not actively dividing, yet are able to return to the cell cycle given the appropriate environmental signals. There is now ample evidence that quiescent G0 cells are not shut-down but still metabolically and transcriptionally active. Quiescent cells must maintain a basal transcriptional capacity to maintain transcripts and proteins necessary for survival. This implies a tight control over RNA polymerases: RNA pol II for mRNA transcription during G0, but especially RNA pol I and RNA pol III to maintain an appropriate level of structural RNAs, raising the possibility that specific transcriptional control mechanisms evolved in quiescent cells. In accordance with this, we recently discovered that RNA interference is necessary to control RNA polymerase I transcription during G0. While this mini-review focuses on yeast model organisms (Saccharomyces cerevisiae and Schizosaccharomyces pombe), parallels are drawn to other eukaryotes and mammalian systems, in particular stem cells.

Histone H3G34R mutation causes replication stress, homologous recombination defects and genomic instability in S. pombe

Recurrent somatic mutations of H3F3A in aggressive pediatric high-grade gliomas generate K27M or G34R/V mutant histone H3.3. H3.3-G34R/V mutants are common in tumors with mutations in p53 and ATRX, an H3.3-specific chromatin remodeler. To gain insight into the role of H3-G34R, we generated fission yeast that express only the mutant histone H3. H3-G34R specifically reduces H3K36 tri-methylation and H3K36 acetylation, and mutants show partial transcriptional overlap with set2 deletions. H3-G34R mutants exhibit genomic instability and increased replication stress, including slowed replication fork restart, although DNA replication checkpoints are functional. H3-G34R mutants are defective for DNA damage repair by homologous recombination (HR), and have altered HR protein dynamics in both damaged and untreated cells. These data suggest H3-G34R slows resolution of HR-mediated repair and that unresolved replication intermediates impair chromosome segregation. This analysis of H3-G34R mutant fission yeast provides mechanistic insight into how G34R mutation may promote genomic instability in glioma.

Transcription-coupled changes to chromatin underpin gene silencing by transcriptional interference

Long non-coding RNA (lncRNA) transcription into a downstream promoter frequently results in transcriptional interference. However, the mechanism of this repression is not fully understood. We recently showed that drug tolerance in fission yeast Schizosaccharomyces pombe is controlled by lncRNA transcription upstream of the tgp1⁺ permease gene. Here we demonstrate that transcriptional interference of tgp1⁺ involves several transcription-coupled chromatin changes mediated by conserved elongation factors Set2, Clr6CII, Spt6 and FACT. These factors are known to travel with RNAPII and establish repressive chromatin in order to limit aberrant transcription initiation from cryptic promoters present in gene bodies. We therefore conclude that conserved RNAPII-associated mechanisms exist to both suppress intragenic cryptic promoters during genic transcription and to repress gene promoters by transcriptional interference. Our analyses also demonstrate that key mechanistic features of transcriptional interference are shared between S. pombe and the highly divergent budding yeast Saccharomyces cerevisiae. Thus, transcriptional interference is an ancient, conserved mechanism for tightly controlling gene expression. Our mechanistic insights allowed us to predict and validate a second example of transcriptional interference involving the S. pombe pho1⁺ gene. Given that eukaryotic genomes are pervasively transcribed, transcriptional interference likely represents a more general feature of gene regulation than is currently appreciated.

Unraveling Site-Specific and Combinatorial Histone Modifications Using High-Resolution Mass Spectrometry in Histone Deacetylase Mutants of Fission Yeast

Histone deacetylases (HDACs) catalyze the removal of acetylation marks from lysine residues on histone and nonhistone substrates. Their activity is generally associated with essential cellular processes such as transcriptional repression and heterochromatin formation. Interestingly, abnormal activity of HDACs has been reported in various types of cancers, which makes them a promising therapeutic target for cancer treatment. In the current study, we aim to understand the mechanisms underlying the function of HDACs using an in-depth quantitative analysis of changes in histone acetylation levels in Schizosaccharomyces pombe (S. pombe) lacking major HDAC activities. We employed a targeted quantitative mass spectrometry approach to profile changes of acetylation and methylation at multiple lysine residues on the N-terminal tail of histones H3 and H4. Our analyses identified a number of histone acetylation sites that are significantly affected by S. pombe HDAC mutations. We discovered that mutation of the Class I HDAC known as Clr6 causes a major increase in the abundance of triacetylated H4 molecules at K5, K8, and K12. A clr6-1 hypomorphic mutation also increased the abundance of multiple acetyl-lysines in histone H3. In addition, our study uncovered a few crosstalks between histone acetylation and methylation upon deletion of HDACs Hos2 and Clr3. We anticipate that the results from this study will greatly improve our current understanding of the mechanisms involved in HDAC-mediated gene regulation and heterochromatin assembly.

Nucleosome alterations caused by mutations at modifiable histone residues in Saccharomyces cerevisiae

Nucleosome organization exhibits dynamic properties depending on the cell state and environment. Histone proteins, fundamental components of nucleosomes, are subject to chemical modifications on particular residues. We examined the effect of substituting modifiable residues of four core histones with the non-modifiable residue alanine on nucleosome dynamics. We mapped the genome-wide nucleosomes in 22 histone mutants of Saccharomyces cerevisiae and compared the nucleosome alterations relative to the wild-type strain. Our results indicated that different types of histone mutation resulted in different phenotypes and a distinct reorganization of nucleosomes. Nucleosome occupancy was altered at telomeres, but not at centromeres. The first nucleosomes upstream (−1) and downstream (+1) of the transcription start site (TSS) were more dynamic than other nucleosomes. Mutations in histones affected the nucleosome array downstream of the TSS. Highly expressed genes, such as ribosome genes and genes involved in glycolysis, showed increased nucleosome occupancy in many types of histone mutant. In particular, the H3K56A mutant exhibited a high percentage of dynamic genomic regions, decreased nucleosome occupancy at telomeres, increased occupancy at the +1 and −1 nucleosomes, and a slow growth phenotype under stress conditions. Our findings provide insight into the influence of histone mutations on nucleosome dynamics.

The SET Domain Is Essential for Metnase Functions in Replication Restart and the 5' End of SS-Overhang Cleavage

Metnase (also known as SETMAR) is a chimeric SET-transposase protein that plays essential role(s) in non-homologous end joining (NHEJ) repair and replication fork restart. Although the SET domain possesses histone H3 lysine 36 dimethylation (H3K36me2) activity associated with an improved association of early repair components for NHEJ, its role in replication restart is less clear. Here we show that the SET domain is necessary for the recovery from DNA damage at the replication forks following hydroxyurea (HU) treatment. Cells overexpressing the SET deletion mutant caused a delay in fork restart after HU release. Our In vitro study revealed that the SET domain but not the H3K36me2 activity is required for the 5’ end of ss-overhang cleavage with fork and non-fork DNA without affecting the Metnase-DNA interaction. Together, our results suggest that the Metnase SET domain has a positive role in restart of replication fork and the 5’ end of ss-overhang cleavage, providing a new insight into the functional interaction of the SET and the transposase domains.

Global and specific responses of the histone acetylome to systematic perturbation

Regulation of histone acetylation is fundamental to the utilization of eukaryotic genomes in chromatin. Aberrant acetylation contributes to disease and can be clinically combated by inhibiting the responsible enzymes. Our knowledge of the histone acetylation system is patchy because we so far lacked the methodology to describe acetylation patterns and their genesis by integrated enzyme activities. We devised a generally applicable, mass spectrometry-based strategy to precisely and accurately quantify combinatorial modification motifs. This was applied to generate a comprehensive inventory of acetylation motifs on histones H3 and H4 in Drosophila cells. Systematic depletion of known or suspected acetyltransferases and deacetylases revealed specific alterations of histone acetylation signatures, established enzyme-substrate relationships, and unveiled an extensive crosstalk between neighboring modifications. Unexpectedly, overall histone acetylation levels remained remarkably constant upon depletion of individual acetyltransferases. Conceivably, the acetylation level is adjusted to maintain the global charge neutralization of chromatin and the stability of nuclei.

Quantitative dynamics of the link between cellular metabolism and histone acetylation

Acetylation on the tails of histones plays an important role in controlling transcription initiation. Although the steady-state abundances of histone acetyl groups have been reported, the rate at which histones are acetylated and deacetylated on a residue-specific basis has not been quantitatively established. We added [(13)C]glucose to human cells and monitored the dynamic incorporation of (13)C-labeled acetyl groups onto specific histone lysines with quantitative mass spectrometry. We determined the turnover of acetylation to be generally slower than phosphorylation, but fast relative to methylation, and that the rate varied depending on the histone, the residue modified, and also the neighboring modifications. Cells were also treated with a deacetylase inhibitor to determine the rate due to histone acetyltransferase activity alone and in the absence of deacetylase activity. Introduction of (13)C-labeled glucose also resulted in the incorporation of (13)C into alanine, which allowed us to partition histones into existing and newly synthesized protein categories. Newly synthesized histones were slower to accumulate histone modifications, especially modifications associated with silent chromatin. Finally, we applied our new approaches to find that quiescent fibroblasts exhibited lower levels of labeled acetyl accumulation compared with proliferating fibroblasts. This suggests that acetylation rates can be modulated in cells in different biological states and that these changes can be detected with the approach presented here. The methods we describe can be broadly applied to defining the turnover of histone acetylation in other cell states such as during cellular reprogramming and to quantify non-histone protein acetylation dynamics.

Mechanism for epigenetic variegation of gene expression at yeast telomeric heterochromatin

Yeast contains heterochromatin at telomeres and the silent mating-type loci (HML/HMR). Genes positioned within the telomeric heterochromatin of Saccharomyces cerevisiae switch stochastically between epigenetically bistable ON and OFF expression states. Important aspects of the mechanism of variegated gene expression, including the chromatin structure of the natural ON state and the mechanism by which it is maintained, are unknown. To address this issue, we developed approaches to select cells in the ON and OFF states. We found by chromatin immunoprecipitation (ChIP) that natural ON telomeres are associated with Rap1 binding and, surprisingly, also contain known characteristics of OFF telomeres, including significant amounts of Sir3 and H4K16 deacetylated nucleosomes. Moreover, we found that H3K79 methylation (H3K79me), H3K4me, and H3K36me, which are depleted from OFF telomeres, are enriched at ON telomeres. We demonstrate in vitro that H3K79me, but not H3K4me or H3K36me, disrupts transcriptional silencing. Importantly, H3K79me does not significantly reduce Sir complex binding in vivo or in vitro. Finally, we show that maintenance of H3K79me at ON telomeres is dependent on transcription. Therefore, although Sir proteins are required for silencing, we propose that epigenetic variegation of telomeric gene expression is due to the bistable enrichment/depletion of H3K79me and not the fluctuation in the amount of Sir protein binding to nucleosomes.

CHD1 remodelers regulate nucleosome spacing in vitro and align nucleosomal arrays over gene coding regions in S. pombe

Nucleosome positioning governs access to eukaryotic genomes. Many genes show a stereotypic organisation at their 5'end: a nucleosome free region just upstream of the transcription start site (TSS) followed by a regular nucleosomal array over the coding region. The determinants for this pattern are unclear, but nucleosome remodelers are likely critical. Here we study the role of remodelers in global nucleosome positioning in S. pombe and the corresponding changes in expression. We find a striking evolutionary shift in remodeler usage between budding and fission yeast. The S. pombe RSC complex does not seem to be involved in nucleosome positioning, despite its prominent role in S. cerevisiae. While S. pombe lacks ISWI-type remodelers, it has two CHD1-type ATPases, Hrp1 and Hrp3. We demonstrate nucleosome spacing activity for Hrp1 and Hrp3 in vitro, and that together they are essential for linking regular genic arrays to most TSSs in vivo. Impaired arrays in the absence of either or both remodelers may lead to increased cryptic antisense transcription, but overall gene expression levels are only mildly affected.

Mapping Post-translational Modifications of Histones H2A, H2B and H4 in Schizosaccharomyces pombe

Core histones are known to carry a variety of post-translational modifications (PTMs), including acetylation, phosphorylation, methylation and ubiquitination, which play important roles in the epigenetic control of gene expression. The nature and biological functions of these PTMs in histones from plants, animals and budding yeast have been extensively investigated. In contrast, the corresponding studies for fission yeast were mainly focused on histone H3. In the present study, we applied LC-nano-ESI-MS/MS, coupled with multiple protease digestion, to identify PTMs in histones H2A, H2B and H4 from Schizosaccharomyces pombe (S. pombe), the typical model organism of fission yeast. Various protease digestions provided high sequence coverage for PTM mapping, and accurate mass measurement of fragment ions allowed for unambiguous differentiation of acetylation from tri-methylation. Many modification sites conserved in other organisms were identified in S. pombe. In addition, some unique modification sites, including N-terminal acetylation in H2A and H2B as well as K123 acetylation in H2A.β, were observed. Our results provide a comprehensive picture of the PTMs of histones H2A, H2B and H4 in S. pombe, which serves as a foundation for future investigations on the regulation and functions of histone modifications in this important model organism.