A Polycomb repressive complex is required for RNAi-mediated heterochromatin formation and dynamic distribution of nuclear bodies.

Polycomb group (PcG) proteins are widely utilized for transcriptional repression in eukaryotes. Here, we characterize, in the protist Tetrahymena thermophila, the EZL1 (E(z)-like 1) complex, with components conserved in metazoan Polycomb Repressive Complexes 1 and 2 (PRC1 and PRC2). The EZL1 complex is required for histone H3 K27 and K9 methylation, heterochromatin formation, transposable element control, and programmed genome rearrangement. The EZL1 complex interacts with EMA1, a helicase required for RNA interference (RNAi). This interaction is implicated in co-transcriptional recruitment of the EZL1 complex. Binding of H3K27 and H3K9 methylation by PDD1-another PcG protein interacting with the EZL1 complex-reinforces its chromatin association. The EZL1 complex is an integral part of Polycomb bodies, which exhibit dynamic distribution in Tetrahymena development: Their dispersion is driven by chromatin association, while their coalescence by PDD1, likely via phase separation. Our results provide a molecular mechanism connecting RNAi and Polycomb repression, which coordinately regulate nuclear bodies and reorganize the genome.

Global gene expression and chromatin accessibility of the peripheral nervous system in animal models of persistent pain.

BACKGROUND: Efforts to understand genetic variability involved in an individual's susceptibility to chronic pain support a role for upstream regulation by epigenetic mechanisms. METHODS: To examine the transcriptomic and epigenetic basis of chronic pain that resides in the peripheral nervous system, we used RNA-seq and ATAC-seq of the rat dorsal root ganglion (DRG) to identify novel molecular pathways associated with pain hypersensitivity in two well-studied persistent pain models induced by chronic constriction injury (CCI) of the sciatic nerve and intra-plantar injection of complete Freund's adjuvant (CFA) in rats. RESULTS: Our RNA-seq studies identify a variety of biological process related to synapse organization, membrane potential, transmembrane transport, and ion binding. Interestingly, genes that encode transcriptional regulators were disproportionately downregulated in both models. Our ATAC-seq data provide a comprehensive map of chromatin accessibility changes in the DRG. A total of 1123 regions showed changes in chromatin accessibility in one or both models when compared to the naïve and 31 shared differentially accessible regions (DAR)s. Functional annotation of the DARs identified disparate molecular functions enriched for each pain model which suggests that chromatin structure may be altered differently following sciatic nerve injury and hind paw inflammation. Motif analysis identified 17 DNA sequences known to bind transcription factors in the CCI DARs and 33 in the CFA DARs. Two motifs were significantly enriched in both models. CONCLUSIONS: Our improved understanding of the changes in chromatin accessibility that occur in chronic pain states may identify regulatory genomic elements that play essential roles in modulating gene expression in the DRG.

Sex differences in gene regulation in the dorsal root ganglion after nerve injury.

BACKGROUND: Pain is a subjective experience derived from complex interactions among biological, environmental, and psychosocial pathways. Sex differences in pain sensitivity and chronic pain prevalence are well established. However, the molecular basis underlying these sex dimorphisms are poorly understood particularly with regard to the role of the peripheral nervous system. Here we sought to identify shared and distinct gene networks functioning in the peripheral nervous systems that may contribute to sex differences of pain in rats after nerve injury. RESULTS: We performed RNA-seq on dorsal root ganglia following chronic constriction injury of the sciatic nerve in male and female rats. Analysis from paired naive and injured tissues showed that 1513 genes were differentially expressed between sexes. Genes which facilitated synaptic transmission in naïve and injured females did not show increased expression in males. CONCLUSIONS: Appreciating sex-related gene expression differences and similarities in neuropathic pain models may help to improve the translational relevance to clinical populations and efficacy of clinical trials of this major health issue.

Hydrazide Mimics for Protein Lysine Acylation To Assess Nucleosome Dynamics and Deubiquitinase Action.

A range of acyl-lysine (acyl-Lys) modifications on histones and other proteins have been mapped over the past decade but for most, their functional and structural significance remains poorly characterized. One limitation in the study of acyl-Lys containing proteins is the challenge of producing them or their mimics in site-specifically modified forms. We describe a cysteine alkylation-based method to install hydrazide mimics of acyl-Lys post-translational modifications (PTMs) on proteins. We have applied this method to install mimics of acetyl-Lys, 2-hydroxyisobutyryl-Lys, and ubiquityl-Lys that could be recognized selectively by relevant acyl-Lys modification antibodies. The acyl-Lys modified histone H3 proteins were reconstituted into nucleosomes to study nucleosome dynamics and stability as a function of modification type and site. We also installed a ubiquityl-Lys mimic in histone H2B and generated a diubiquitin analog, both of which could be cleaved by deubiquitinating enzymes. Nucleosomes containing the H2B ubiquityl-Lys mimic were used to study the SAGA deubiquitinating module's molecular recognition. These results suggest that acyl-Lys mimics offer a relatively simple and promising strategy to study the role of acyl-Lys modifications in the function, structure, and regulation of proteins and protein complexes.

RNA-seq of spinal cord from nerve-injured rats after spinal cord stimulation.

Spinal cord stimulation has become an important modality in pain treatment especially for neuropathic pain conditions refractory to pharmacotherapy. However, the molecular control of inhibitory and excitatory mechanisms observed after spinal cord stimulation are poorly understood. Here, we used RNA-seq to identify differences in the expression of genes and gene networks in spinal cord tissue from nerve-injured rats with and without repetitive conventional spinal cord stimulation treatment. Five weeks after chronic constrictive injury to the left sciatic nerve, male and female rats were randomized to receive repetitive spinal cord stimulation or no treatment. Rats receiving spinal cord stimulation underwent epidural placement of a miniature stimulating electrode and received seven sessions of spinal cord stimulation (50 Hz, 80% motor threshold, 0.2 ms, constant current bipolar stimulation, 120 min/session) over four consecutive days. Within 2 h after the last spinal cord stimulation treatment, the L4-L6 spinal segments ipsilateral to the side of nerve injury were harvested and used to generate libraries for RNA-seq. Our RNA-seq data suggest further increases of many existing upregulated immune responses in chronic constrictive injury rats after repetitive spinal cord stimulation, including transcription of cell surface receptors and activation of non-neuronal cells. We also demonstrate that repetitive spinal cord stimulation represses transcription of several key synaptic signaling genes that encode scaffold proteins in the post-synaptic density. Our transcriptional studies suggest a potential relationship between specific genes and the therapeutic effects observed in patients undergoing conventional spinal cord stimulation after nerve injury. Furthermore, our results may help identify new therapeutic targets for improving the efficacy of conventional spinal cord stimulation and other chronic pain treatments.

Dissecting relative contributions of cis- and trans-determinants to nucleosome distribution by comparing Tetrahymena macronuclear and micronuclear chromatin.

The ciliate protozoan Tetrahymena thermophila contains two types of structurally and functionally differentiated nuclei: the transcriptionally active somatic macronucleus (MAC) and the transcriptionally silent germ-line micronucleus (MIC). Here, we demonstrate that MAC features well-positioned nucleosomes downstream of transcription start sites and flanking splice sites. Transcription-associated trans-determinants promote nucleosome positioning in MAC. By contrast, nucleosomes in MIC are dramatically delocalized. Nucleosome occupancy in MAC and MIC are nonetheless highly correlated with each other, as well as with in vitro reconstitution and predictions based upon DNA sequence features, revealing unexpectedly strong contributions from cis-determinants. In particular, well-positioned nucleosomes are often matched with GC content oscillations. As many nucleosomes are coordinately accommodated by both cis- and trans-determinants, we propose that their distribution is shaped by the impact of these nucleosomes on the mutational and transcriptional landscape, and driven by evolutionary selection.

Nucleosome competition reveals processive acetylation by the SAGA HAT module.

The Spt-Ada-Gcn5 acetyltransferase (SAGA) coactivator complex hyperacetylates histone tails in vivo in a manner that depends upon histone 3 lysine 4 trimethylation (H3K4me3), a histone mark enriched at promoters of actively transcribed genes. SAGA contains a separable subcomplex known as the histone acetyltransferase (HAT) module that contains the HAT, Gcn5, bound to Sgf29, Ada2, and Ada3. Sgf29 contains a tandem Tudor domain that recognizes H3K4me3-containing peptides and is required for histone hyperacetylation in vivo. However, the mechanism by which H3K4me3 recognition leads to lysine hyperacetylation is unknown, as in vitro studies show no effect of the H3K4me3 modification on histone peptide acetylation by Gcn5. To determine how H3K4me3 binding by Sgf29 leads to histone hyperacetylation by Gcn5, we used differential fluorescent labeling of histones to monitor acetylation of individual subpopulations of methylated and unmodified nucleosomes in a mixture. We find that the SAGA HAT module preferentially acetylates H3K4me3 nucleosomes in a mixture containing excess unmodified nucleosomes and that this effect requires the Tudor domain of Sgf29. The H3K4me3 mark promotes processive, multisite acetylation of histone H3 by Gcn5 that can account for the different acetylation patterns established by SAGA at promoters versus coding regions. Our results establish a model for Sgf29 function at gene promoters and define a mechanism governing crosstalk between histone modifications.

The bromodomain of Gcn5 regulates site specificity of lysine acetylation on histone H3.

In yeast, the conserved histone acetyltransferase (HAT) Gcn5 associates with Ada2 and Ada3 to form the catalytic module of the ADA and SAGA transcriptional coactivator complexes. Gcn5 also contains an acetyl-lysine binding bromodomain that has been implicated in regulating nucleosomal acetylation in vitro, as well as at gene promoters in cells. However, the contribution of the Gcn5 bromodomain in regulating site specificity of HAT activity remains unclear. Here, we used a combined acid-urea gel and quantitative mass spectrometry approach to compare the HAT activity of wild-type and Gcn5 bromodomain-mutant ADA subcomplexes (Gcn5-Ada2-Ada3). Wild-type ADA subcomplex acetylated H3 lysines with the following specificity; H3K14 > H3K23 > H3K9 ≈ H3K18 > H3K27 > H3K36. However, when the Gcn5 bromodomain was defective in acetyl-lysine binding, the ADA subcomplex demonstrated altered site-specific acetylation on free and nucleosomal H3, with H3K18ac being the most severely diminished. H3K18ac was also severely diminished on H3K14R, but not H3K23R, substrates in wild-type HAT reactions, further suggesting that Gcn5-catalyzed acetylation of H3K14 and bromodomain binding to H3K14ac are important steps preceding H3K18ac. In sum, this work details a previously uncharacterized cross-talk between the Gcn5 bromodomain "reader" function and enzymatic HAT activity that might ultimately affect gene expression. Future studies of how mutations in bromodomains or other histone post-translational modification readers can affect chromatin-templated enzymatic activities will yield unprecedented insight into a potential "histone/epigenetic code." MS data are available via ProteomeXchange with identifier PXD001167.

A PWWP domain-containing protein targets the NuA3 acetyltransferase complex via histone H3 lysine 36 trimethylation to coordinate transcriptional elongation at coding regions.

Post-translational modifications of histones, such as acetylation and methylation, are differentially positioned in chromatin with respect to gene organization. For example, although histone H3 is often trimethylated on lysine 4 (H3K4me3) and acetylated on lysine 14 (H3K14ac) at active promoter regions, histone H3 lysine 36 trimethylation (H3K36me3) occurs throughout the open reading frames of transcriptionally active genes. The conserved yeast histone acetyltransferase complex, NuA3, specifically binds H3K4me3 through a plant homeodomain (PHD) finger in the Yng1 subunit, and subsequently catalyzes the acetylation of H3K14 through the histone acetyltransferase domain of Sas3, leading to transcription initiation at a subset of genes. We previously found that Ylr455w (Pdp3), an uncharacterized proline-tryptophan-tryptophan-proline (PWWP) domain-containing protein, copurifies with stable members of NuA3. Here, we employ mass-spectrometric analysis of affinity purified Pdp3, biophysical binding assays, and genetic analyses to classify NuA3 into two functionally distinct forms: NuA3a and NuA3b. Although NuA3a uses the PHD finger of Yng1 to interact with H3K4me3 at the 5'-end of open reading frames, NuA3b contains the unique member, Pdp3, which regulates an interaction between NuA3b and H3K36me3 at the transcribed regions of genes through its PWWP domain. We find that deletion of PDP3 decreases NuA3-directed transcription and results in growth defects when combined with transcription elongation mutants, suggesting NuA3b acts as a positive elongation factor. Finally, we determine that NuA3a, but not NuA3b, is synthetically lethal in combination with a deletion of the histone acetyltransferase GCN5, indicating NuA3b has a specialized role at coding regions that is independent of Gcn5 activity. Collectively, these studies define a new form of the NuA3 complex that associates with H3K36me3 to effect transcriptional elongation. MS data are available via ProteomeXchange with identifier PXD001156.

Purification of a specific native genomic locus for proteomic analysis.

Here, we describe an approach to isolate native chromatin sections without genomic engineering for label-free proteomic identification of associated proteins and histone post-translational modifications. A transcription activator-like (TAL) protein A fusion protein was designed to recognize a unique site in the yeast GAL1 promoter. The TAL-PrA fusion enabled chromatin affinity purification (ChAP) of a small section of native chromatin upstream from the GAL1 locus, permitting mass spectrometric (MS) identification of proteins and histone post-translational modifications regulating galactose-induced transcription. This TAL-ChAP-MS approach allows the biochemical isolation of a specific native genomic locus for proteomic studies and will provide for unprecedented objective insight into protein and epigenetic mechanisms regulating site-specific chromosome metabolism.

Using S-adenosyl-L-homocysteine capture compounds to characterize S-adenosyl-L-methionine and S-adenosyl-L-homocysteine binding proteins.

S-Adenosyl-l-methionine (SAM) is recognized as an important cofactor in a variety of biochemical reactions. As more proteins and pathways that require SAM are discovered, it is important to establish a method to quickly identify and characterize SAM binding proteins. The affinity of S-adenosyl-l-homocysteine (SAH) for SAM binding proteins was used to design two SAH-derived capture compounds (CCs). We demonstrate interactions of the proteins COMT and SAHH with SAH-CC with biotin used in conjunction with streptavidin-horseradish peroxidase. After demonstrating SAH-dependent photo-crosslinking of the CC to these proteins, we used a CC labeled with a fluorescein tag to measure binding affinity via fluorescence anisotropy. We then used this approach to show and characterize binding of SAM to the PR domain of PRDM2, a lysine methyltransferase with putative tumor suppressor activity. We calculated the Kd values for COMT, SAHH, and PRDM2 (24.1 ± 2.2 µM, 6.0 ± 2.9 µM, and 10.06 ± 2.87 µM, respectively) and found them to be close to previously established Kd values of other SAM binding proteins. Here, we present new methods to discover and characterize SAM and SAH binding proteins using fluorescent CCs.

Identification of Regulatory Elements in Primary Sensory Neurons Involved in Trauma-Induced Neuropathic Pain.

Chronic pain is a significant public health issue that is often refractory to existing therapies. Here we use a multiomic approach to identify cis-regulatory elements that show differential chromatin accessibility and reveal transcription factor (TF) binding motifs with functional regulation in the rat dorsal root ganglion (DRG), which contain cell bodies of primary sensory neurons, after nerve injury. We integrated RNA-seq to understand how differential chromatin accessibility after nerve injury may influence gene expression. Using TF protein arrays and chromatin immunoprecipitation-qPCR, we confirmed C/EBPγ binding to a differentially accessible sequence and used RNA-seq to identify processes in which C/EBPγ plays an important role. Our findings offer insights into TF motifs that are associated with chronic pain. These data show how interactions between chromatin landscapes and TF expression patterns may work together to determine gene expression programs in rat DRG neurons after nerve injury.

Molecular insight into interactions between the Taf14, Yng1 and Sas3 subunits of the NuA3 complex.

The NuA3 complex is a major regulator of gene transcription and the cell cycle in yeast. Five core subunits are required for complex assembly and function, but it remains unclear how these subunits interact to form the complex. Here, we report that the Taf14 subunit of the NuA3 complex binds to two other subunits of the complex, Yng1 and Sas3, and describe the molecular mechanism by which the extra-terminal domain of Taf14 recognizes the conserved motif present in Yng1 and Sas3. Structural, biochemical, and mutational analyses show that two motifs are sandwiched between the two extra-terminal domains of Taf14. The head-to-toe dimeric complex enhances the DNA binding activity of Taf14, and the formation of the hetero-dimer involving the motifs of Yng1 and Sas3 is driven by sequence complementarity. In vivo assays in yeast demonstrate that the interactions of Taf14 with both Sas3 and Yng1 are required for proper function of the NuA3 complex in gene transcription and DNA repair. Our findings suggest a potential basis for the assembly of three core subunits of the NuA3 complex, Taf14, Yng1 and Sas3.

Azalysine analogues as probes for protein lysine deacetylation and demethylation.

Reversible lysine acetylation and methylation regulate the function of a wide variety of proteins, including histones. Here, we have synthesized azalysine-containing peptides in acetylated and unacetylated forms as chemical probes of the histone deacetylases (HDAC8, Sir2Tm, and SIRT1) and the histone demethylase, LSD1. We have shown that the acetyl-azalysine modification is a fairly efficient substrate for the sirtuins, but a weaker substrate for HDAC8, a classical HDAC. In addition to deacetylation by sirtuins, the acetyl-azalysine analogue generates a novel ADP-ribose adduct that was characterized by mass spectrometry, Western blot analysis, and nuclear magnetic resonance spectroscopy. This peptide-ADP-ribose adduct is proposed to correspond to a derailed reaction intermediate, providing unique evidence for the direct 2'-hydroxyl attack on the O-alkylimidate intermediate that is formed in the course of sirtuin catalyzed deacetylation. An unacetylated azalysine-containing H3 peptide proved to be a potent inhibitor of the LSD1 demethylase, forming an FAD adduct characteristic of previously reported related structures, providing a new chemical probe for mechanistic analysis.

De novo methylation of histone H3K23 by the methyltransferases EHMT1/GLP and EHMT2/G9a.

Epigenetic modifications to histone proteins serve an important role in regulating permissive and repressive chromatin states, but despite the identification of many histone PTMs and their perceived role, the epigenetic writers responsible for generating these chromatin signatures are not fully characterized. Here, we report that the canonical histone H3K9 methyltransferases EHMT1/GLP and EHMT2/G9a are capable of catalyzing methylation of histone H3 lysine 23 (H3K23). Our data show that while both enzymes can mono- and di-methylate H3K23, only EHMT1/GLP can tri-methylate H3K23. We also show that pharmacologic inhibition or genetic ablation of EHMT1/GLP and/or EHMT2/G9a leads to decreased H3K23 methylation in mammalian cells. Taken together, this work identifies H3K23 as a new direct methylation target of EHMT1/GLP and EHMT2/G9a, and highlights the differential activity of these enzymes on H3K23 as a substrate.

Human brain region-specific variably methylated regions are enriched for heritability of distinct neuropsychiatric traits.

BACKGROUND: DNA methylation dynamics in the brain are associated with normal development and neuropsychiatric disease and differ across functionally distinct brain regions. Previous studies of genome-wide methylation differences among human brain regions focus on limited numbers of individuals and one to two brain regions. RESULTS: Using GTEx samples, we generate a resource of DNA methylation in purified neuronal nuclei from 8 brain regions as well as lung and thyroid tissues from 12 to 23 donors. We identify differentially methylated regions between brain regions among neuronal nuclei in both CpG (181,146) and non-CpG (264,868) contexts, few of which were unique to a single pairwise comparison. This significantly expands the knowledge of differential methylation across the brain by 10-fold. In addition, we present the first differential methylation analysis among neuronal nuclei from basal ganglia tissues and identify unique CpG differentially methylated regions, many associated with ion transport. We also identify 81,130 regions of variably CpG methylated regions, i.e., variable methylation among individuals in the same brain region, which are enriched in regulatory regions and in CpG differentially methylated regions. Many variably methylated regions are unique to a specific brain region, with only 202 common across all brain regions, as well as lung and thyroid. Variably methylated regions identified in the amygdala, anterior cingulate cortex, and hippocampus are enriched for heritability of schizophrenia. CONCLUSIONS: These data suggest that epigenetic variation in these particular human brain regions could be associated with the risk for this neuropsychiatric disorder.

Phosphorylation of the SQ H2A.X motif is required for proper meiosis and mitosis in Tetrahymena thermophila.

Phosphorylation of the C terminus SQ motif that defines H2A.X variants is required for efficient DNA double-strand break (DSB) repair in diverse organisms but has not been studied in ciliated protozoa. Tetrahymena H2A.X is one of two similarly expressed major H2As, thereby differing both from mammals, where H2A.X is a quantitatively minor component, and from Saccharomyces cerevisiae where it is the only type of major H2A. Tetrahymena H2A.X is phosphorylated in the SQ motif in both the mitotic micronucleus and the amitotic macronucleus in response to DSBs induced by chemical agents and in the micronucleus during prophase of meiosis, which occurs in the absence of a synaptonemal complex. H2A.X is phosphorylated when programmed DNA rearrangements occur in developing macronuclei, as for immunoglobulin gene rearrangements in mammals, but not during the DNA fragmentation that accompanies breakdown of the parental macronucleus during conjugation, correcting the previous interpretation that this process is apoptosis-like. Using strains containing a mutated (S134A) SQ motif, we demonstrate that phosphorylation of this motif is important for Tetrahymena cells to recover from exogenous DNA damage and is required for normal micronuclear meiosis and mitosis and, to a lesser extent, for normal amitotic macronuclear division; its absence, while not lethal, leads to the accumulation of DSBs in both micro- and macronuclei. These results demonstrate multiple roles of H2A.X phosphorylation in maintaining genomic integrity in different phases of the Tetrahymena life cycle.

Genome-wide Cas9 binding specificity in Saccharomyces cerevisiae.

The CRISPR system has become heavily utilized in biomedical research as a tool for genomic editing as well as for site-specific chromosomal localization of specific proteins. For example, we developed a CRISPR-based methodology for enriching a specific genomic locus of interest for proteomic analysis in Saccharomyces cerevisiae, which utilized a guide RNA-targeted, catalytically dead Cas9 (dCas9) as an affinity reagent. To more comprehensively evaluate the genomic specificity of using dCas9 as a site-specific tool for chromosomal studies, we performed dCas9-mediated locus enrichment followed by next-generation sequencing on a genome-wide scale. As a test locus, we used the ARS305 origin of replication on chromosome III in S. cerevisiae. We found that enrichment of this site is highly specific, with virtually no off-target enrichment of unique genomic sequences. The high specificity of genomic localization and enrichment suggests that dCas9-mediated technologies have promising potential for site-specific chromosomal studies in organisms with relatively small genomes such as yeasts.

Saccharomyces cerevisiae Yta7 regulates histone gene expression.

The Saccharomyces cerevisiae Yta7 protein is a component of a nucleosome bound protein complex that maintains distinct transcriptional zones of chromatin. We previously found that one protein copurifying with Yta7 is the yFACT member Spt16. Epistasis analyses revealed a link between Yta7, Spt16, and other previously identified members of the histone regulatory pathway. In concurrence, Yta7 was found to regulate histone gene transcription in a cell-cycle-dependent manner. Association at the histone gene loci appeared to occur through binding of the bromodomain-like region of Yta7 with the N-terminal tail of histone H3. Our work suggests a mechanism in which Yta7 is localized to chromatin to establish regions of transcriptional silencing, and that one facet of this cellular mechanism is to modulate transcription of histone genes.