Regulation of Nuclear Epigenome by Mitochondrial DNA Heteroplasmy

Total Page:16

File Type:pdf, Size:1020Kb

Regulation of Nuclear Epigenome by Mitochondrial DNA Heteroplasmy Regulation of nuclear epigenome by mitochondrial DNA heteroplasmy Piotr K. Kopinskia,b,c, Kevin A. Janssend,e, Patrick M. Schaeferc, Sophie Trefelyf,g, Caroline E. Perryc, Prasanth Potluric, Jesus A. Tintos-Hernandezc, Larry N. Singhc, Kelly R. Karchd,e, Sydney L. Campbellf, Mary T. Doang, Helen Jiangg, Itzhak Nissime, Eiko Nakamaru-Ogisoe, Kathryn E. Wellenf, Nathaniel W. Snyderg, Benjamin A. Garciad,e, and Douglas C. Wallacec,h,1 aHoward Hughes Medical Institute, University of Pennsylvania, Philadelphia, PA 19104; bPerelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104; cCenter for Mitochondrial and Epigenomic Medicine, Children’s Hospital of Philadelphia, Philadelphia, PA 19104; dEpigenetics Institute, University of Pennsylvania, Philadelphia, PA 19104; eDepartment of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104; fDepartment of Cancer Biology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104; gA.J. Drexel Autism Institute, Drexel University, Philadelphia, PA 19104; and hDepartment of Pediatrics, Division of Human Genetics, The Children’s Hospital of Philadelphia, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104 Contributed by Douglas C. Wallace, May 27, 2019 (sent for review April 23, 2019; reviewed by Yidong Bai, Scott W. Ballinger, and Victor Faundez) Diseases associated with mitochondrial DNA (mtDNA) mutations are cardiomyopathy, lactic acidosis, and stroke-like episodes, pre- highly variable in phenotype, in large part because of differences in the senting separately or when in aggregate referred to as the percentage of normal and mutant mtDNAs (heteroplasmy) present MELAS syndrome. At 90% to 100% mutant, individuals manifest within the cell. For example, increasing heteroplasmy levels of the perinatal lethal diseases such as Leigh syndrome. Leu(UUR) mtDNA tRNA nucleotide (nt) 3243A > G mutation result succes- To investigate the molecular basis of the relationship between sively in diabetes, neuromuscular degenerative disease, and perinatal the mtDNA heteroplasmy levels and shifts in patient phenotype, lethality. These phenotypes are associated with differences in mitochon- we prepared transmitochondrial cybrids with the same nucleus, drial function and nuclear DNA (nDNA) gene expression, which are re- but with different heteroplasmic percentages of the tRNALeu(UUR) capitulated in cybrid cell lines with different percentages of m.3243G mutant mtDNAs. Using metabolic tracing, histone mass spectrometry, m.3243G mutation. Cell biological and biochemical character- and NADH fluorescence lifetime imaging microscopy in these cells, we ization of these cybrid cell lines revealed that the physiological now show that increasing levels of this single mtDNA mutation cause parameters of the cybrids mirrored those previously observed in profound changes in the nuclear epigenome. At high heteroplasmy, the patients with the same ranges of m.3243G heteroplasmy. mitochondrially derived acetyl-CoA levels decrease causing decreased Transcriptional analysis by RNA sequencing revealed that cybrids histone H4 acetylation, with glutamine-derived acetyl-CoA compensat- ing when glucose-derived acetyl-CoA is limiting. In contrast, Significance α-ketoglutarate levels increase at midlevel heteroplasmy and are in- versely correlated with histone H3 methylation. Inhibition of mitochon- Mitochondrial DNA (mtDNA) mutations have been associated drial protein synthesis induces acetylation and methylation changes, and with common metabolic and degenerative disease phenotypes, restoration of mitochondrial function reverses these effects. mtDNA het- + implying a bioenergetic etiology for these diseases. For example, eroplasmy also affects mitochondrial NAD /NADH ratio, which corre- tRNALeu(UUR) > + the mtDNA m.3243A G mutation manifests as di- lates with nuclear histone acetylation, whereas nuclear NAD /NADH abetes, neurodegenerative disease, or lethal pediatric disease, ratio correlates with changes in nDNA and mtDNA transcription. Thus, depending on the percentage of mutant mtDNAs within the cell mutations in the mtDNA cause distinct metabolic and epigenomic (heteroplasmy). Cultured cybrid cell lines harboring 3243G heter- changes at different heteroplasmy levels, potentially explaining tran- oplasmy levels corresponding to the different clinical phenotypes scriptional and phenotypic variability of mitochondrial disease. have distinct transcriptional profiles. Exhaustive metabolomic and histone posttranscriptional modification analysis of these 3243G mitochondria | epigenetics | metabolism | transcription | common diseases cybrids revealed that changes in mtDNA heteroplasmy cause changes in mitochondrial intermediates and redox state, which he mtDNA is maternally inherited and codes for 13 essential result in distinctive histone modification changes. Thus, changes in Tsubunits of the mitochondrial energy generating system, ox- the mitochondrial genotype change mitochondrial metabolism, idative phosphorylation (OXPHOS), plus the 2 rRNAs for the which change the epigenome and transcriptome, which induce mitochondrial, chloramphenicol (CP)-sensitive, ribosome and distinct clinical phenotypes. the 22 tRNAs required for mitochondrial mRNA translation. Mutations in the mtDNA tRNAs impair the biogenesis of all Author contributions: P.K.K. and D.C.W. designed research; P.K.K., K.A.J., P.M.S., S.T., mtDNA-encoded OXPHOS enzyme complexes, but particularly C.E.P., P.P., J.A.T.-H., K.R.K., S.L.C., M.T.D., H.J., I.N., E.N.-O., N.W.S., and B.A.G. performed research; P.K.K., K.A.J., P.M.S., L.N.S., K.E.W., N.W.S., B.A.G., and D.C.W. analyzed data; compromise complex I (NADH:CoQ oxidoreductase), which and P.K.K. and D.C.W. wrote the paper. incorporates 7 of the 13 mtDNA encoded polypeptides (1). Reviewers: Y.B., University of Texas Health Science Center at San Antonio; S.W.B., Univer- The overall prevalence of mtDNA mutations is about 1 in sity of Alabama at Birmingham; and V.F., Emory University. 5,000 (2, 3). This high prevalence of mtDNA diseases results The authors declare no conflict of interest. from the high mtDNA mutation rate, which is reflected in 1 in This open access article is distributed under Creative Commons Attribution-NonCommercial- 200 cord bloods harboring 1 of the 10 most common pathogenic NoDerivatives License 4.0 (CC BY-NC-ND). mtDNA mutations (4, 5) and 1 in 8 adults harboring a disease- Data deposition: Histone Post-Translational Modification Data: The current Chorus mass associated mutation at ≥1% heteroplasmy (6). spectrometry database is labeled project 1569. It is publicly shared at the URL https:// One of the most common pathogenic mtDNA mutations is the chorusproject.org/pages/dashboard.html#/projects/public/1569/experiments. Leu(UUR) tRNA m.3243 A > G mutation, which causes starkly See Commentary on page 15763. different clinical phenotypes depending on the heteroplasmy 1To whom correspondence may be addressed. Email: [email protected]. levels (1). At 20% to 30% m.3243G mutant, patients com- This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. monly present with type 1 or type 2 diabetes or manifest au- 1073/pnas.1906896116/-/DCSupplemental. tism. At 50% to 80% mutant, patients present with myopathy, Published online June 28, 2019. 16028–16035 | PNAS | August 6, 2019 | vol. 116 | no. 32 www.pnas.org/cgi/doi/10.1073/pnas.1906896116 Downloaded by guest on September 30, 2021 A R−value 1 + SEE COMMENTARY 0 Lactate (L) P/L Citrate aKG Succinyl-CoA Succinate Fumarate Malate Glutamate (E) Acetyl-CoA Pyruvate (P) - 1 Isocitrate aKG/suc Glutamine (Q) Q/E ac aKG ~ H3K9me3 H3 K9 me2 2. 0.2 R = -0.45 me3 p < 0.05 K14 ac 2 aKG/suc (ratio) ac K18 me1 0.1 K23 ac 1 ac me1 K27 me2 H3K9me3 (fraction) me3 0.0 0 me1 K36 0255075100 me3 A3243G (%) K27 me2 H3.3 K36 me3 K5 ac Acetyl-CoA ~ H4K16ac H4 K8 ac K12 ac1. R = 0.57 100 K16 ac 0.3 p < 0.01 [pmol] Acetyl-CoA ac 75 K9 H2A.1 me1 GENETICS ac 0.2 K13 50 me1 K15 ac 0.1 T16 ac 25 H4K16ac (fraction) K7 ac H2A.V K11 ac 0.0 0 K7 ac 0255075100 H2A.Z K11 ac A3243G (%) B 2 Relative + 1 Fraction C p-value p-value Histone Modification 0 (adjusted) 1 H3_27_40 K27ac 0.0010 0.0418 - 2 H4_4_17 K5acK16ac 0.0004 0.0418 H4_4_17 K8acK16ac 0.0003 0.0418 H4 K5acK8acK12ac H4_4_17 K8acK12acK16ac 0.0011 0.0418 H4 K5ac H2AZ_1_19 K4ac 0.0008 0.0418 H4 K5acK8acK16ac H4_4_17 K12acK16ac 0.0013 0.0434 H4 K5acK12acK16ac H4 K8acK12acK16ac H3_27_40 unmod 0.0018 0.0495 H4 K8acK16ac H4 K5acK16ac H4_4_17 K5acK12ac 0.0020 0.0495 H4 K12acK16ac H4_4_17 K5acK8acK16ac 0.0029 0.0570 H4 K5acK12ac H4_4_17 K5acK12acK16ac 0.0028 0.0570 H4 K5acK8ac H4 K16ac H33_27_40 K27ac 0.0033 0.0599 H4_4_17 unmod 0.0040 0.0666 H2AZ_1_19 K7ac 0.0049 0.0746 H3_18_26 K23me1 0.0093 0.1307 H4_4_17 K5ac 0.0105 0.1307 H4 K5acK8acK12acK16ac H4_4_17 K16ac 0.0103 0.1307 H3_9_17 K14ac 0.0138 0.1618 H2AZ_1_19 unmod 0.0196 0.2167 H4 unmodified H33_27_40 unmod 0.0307 0.2671 A3243G (%) 0 0 0 100 100 100 H14_25_32 unmod 0.0296 0.2671 Fig. 1. MtDNA heteroplasmy modulates metabolites and histone modifications. (A) Seven mtDNA m.3243A > G heteroplasmic cell lines plus the ρ0 parent cells were cultured in parallel for metabolite (j = 15) and histone modification (i = 107) quantification by liquid chromatography tandem mass spectrometry. Spearman correlation coefficient was calculated for each histone modification (unmod = unmodified peptide) against each metabolite (blue = positive, red = negative, and blank = no cor- relation), with the color intensity corresponding to correlation strength (significance = P < 0.05). (Example #1, blue box) Positive association between acetyl-CoA levels and histone H4 acetylation of lysines 8 (K8) and 16 (K16) across heteroplasmy levels. (Example #2, red box) Negative association between αKG/succinate ratio and histone H3 lysine 9 di- (me2) and trimethylation (me3) across heteroplasmy levels. Each point represents a measurement and lines represent local regression of the mean (R-software loess), n = 3. (B) Heat map of histone acetylation modifications of 0% and 100% m.3243G lines (n = 3), showing a dramatic effect of mtDNA genotype on the histone epigenome.
Recommended publications
  • The C-Terminal Region of Histone H4 Recognized by Antinucleosome Th Cells in Identification of a Minimal T Cell Epitope
    Identification of a Minimal T Cell Epitope Recognized by Antinucleosome Th Cells in the C-Terminal Region of Histone H4 This information is current as Patrice Decker, Anne Le Moal, Jean-Paul Briand and of October 2, 2021. Sylviane Muller J Immunol 2000; 165:654-662; ; doi: 10.4049/jimmunol.165.2.654 http://www.jimmunol.org/content/165/2/654 Downloaded from References This article cites 34 articles, 13 of which you can access for free at: http://www.jimmunol.org/content/165/2/654.full#ref-list-1 http://www.jimmunol.org/ Why The JI? Submit online. • Rapid Reviews! 30 days* from submission to initial decision • No Triage! Every submission reviewed by practicing scientists • Fast Publication! 4 weeks from acceptance to publication *average by guest on October 2, 2021 Subscription Information about subscribing to The Journal of Immunology is online at: http://jimmunol.org/subscription Permissions Submit copyright permission requests at: http://www.aai.org/About/Publications/JI/copyright.html Email Alerts Receive free email-alerts when new articles cite this article. Sign up at: http://jimmunol.org/alerts The Journal of Immunology is published twice each month by The American Association of Immunologists, Inc., 1451 Rockville Pike, Suite 650, Rockville, MD 20852 Copyright © 2000 by The American Association of Immunologists All rights reserved. Print ISSN: 0022-1767 Online ISSN: 1550-6606. Identification of a Minimal T Cell Epitope Recognized by Antinucleosome Th Cells in the C-Terminal Region of Histone H41 Patrice Decker, Anne Le Moal, Jean-Paul Briand, and Sylviane Muller2 Autoreactive T cells responding to systemic autoantigens have been characterized in patients and mice with autoimmune diseases and in healthy individuals.
    [Show full text]
  • Remodeling Adipose Tissue Through in Silico Modulation of Fat Storage For
    Chénard et al. BMC Systems Biology (2017) 11:60 DOI 10.1186/s12918-017-0438-9 RESEARCHARTICLE Open Access Remodeling adipose tissue through in silico modulation of fat storage for the prevention of type 2 diabetes Thierry Chénard2, Frédéric Guénard3, Marie-Claude Vohl3,4, André Carpentier5, André Tchernof4,6 and Rafael J. Najmanovich1* Abstract Background: Type 2 diabetes is one of the leading non-infectious diseases worldwide and closely relates to excess adipose tissue accumulation as seen in obesity. Specifically, hypertrophic expansion of adipose tissues is related to increased cardiometabolic risk leading to type 2 diabetes. Studying mechanisms underlying adipocyte hypertrophy could lead to the identification of potential targets for the treatment of these conditions. Results: We present iTC1390adip, a highly curated metabolic network of the human adipocyte presenting various improvements over the previously published iAdipocytes1809. iTC1390adip contains 1390 genes, 4519 reactions and 3664 metabolites. We validated the network obtaining 92.6% accuracy by comparing experimental gene essentiality in various cell lines to our predictions of biomass production. Using flux balance analysis under various test conditions, we predict the effect of gene deletion on both lipid droplet and biomass production, resulting in the identification of 27 genes that could reduce adipocyte hypertrophy. We also used expression data from visceral and subcutaneous adipose tissues to compare the effect of single gene deletions between adipocytes from each
    [Show full text]
  • METABOLIC EVOLUTION in GALDIERIA SULPHURARIA By
    METABOLIC EVOLUTION IN GALDIERIA SULPHURARIA By CHAD M. TERNES Bachelor of Science in Botany Oklahoma State University Stillwater, Oklahoma 2009 Submitted to the Faculty of the Graduate College of the Oklahoma State University in partial fulfillment of the requirements for the Degree of DOCTOR OF PHILOSOPHY May, 2015 METABOLIC EVOLUTION IN GALDIERIA SUPHURARIA Dissertation Approved: Dr. Gerald Schoenknecht Dissertation Adviser Dr. David Meinke Dr. Andrew Doust Dr. Patricia Canaan ii Name: CHAD M. TERNES Date of Degree: MAY, 2015 Title of Study: METABOLIC EVOLUTION IN GALDIERIA SULPHURARIA Major Field: PLANT SCIENCE Abstract: The thermoacidophilic, unicellular, red alga Galdieria sulphuraria possesses characteristics, including salt and heavy metal tolerance, unsurpassed by any other alga. Like most plastid bearing eukaryotes, G. sulphuraria can grow photoautotrophically. Additionally, it can also grow solely as a heterotroph, which results in the cessation of photosynthetic pigment biosynthesis. The ability to grow heterotrophically is likely correlated with G. sulphuraria ’s broad capacity for carbon metabolism, which rivals that of fungi. Annotation of the metabolic pathways encoded by the genome of G. sulphuraria revealed several pathways that are uncharacteristic for plants and algae, even red algae. Phylogenetic analyses of the enzymes underlying the metabolic pathways suggest multiple instances of horizontal gene transfer, in addition to endosymbiotic gene transfer and conservation through ancestry. Although some metabolic pathways as a whole appear to be retained through ancestry, genes encoding individual enzymes within a pathway were substituted by genes that were acquired horizontally from other domains of life. Thus, metabolic pathways in G. sulphuraria appear to be composed of a ‘metabolic patchwork’, underscored by a mosaic of genes resulting from multiple evolutionary processes.
    [Show full text]
  • The Role of Histone H2av Variant Replacement and Histone H4 Acetylation in the Establishment of Drosophila Heterochromatin
    The role of histone H2Av variant replacement and histone H4 acetylation in the establishment of Drosophila heterochromatin Jyothishmathi Swaminathan, Ellen M. Baxter, and Victor G. Corces1 Department of Biology, Johns Hopkins University, Baltimore, Maryland 21218, USA Activation and repression of transcription in eukaryotes involve changes in the chromatin fiber that can be accomplished by covalent modification of the histone tails or the replacement of the canonical histones with other variants. Here we show that the histone H2A variant of Drosophila melanogaster, H2Av, localizes to the centromeric heterochromatin, and it is recruited to an ectopic heterochromatin site formed by a transgene array. His2Av behaves genetically as a PcG gene and mutations in His2Av suppress position effect variegation (PEV), suggesting that this histone variant is required for euchromatic silencing and heterochromatin formation. His2Av mutants show reduced acetylation of histone H4 at Lys 12, decreased methylation of histone H3 at Lys 9, and a reduction in HP1 recruitment to the centromeric region. H2Av accumulation or histone H4 Lys 12 acetylation is not affected by mutations in Su(var)3-9 or Su(var)2-5. The results suggest an ordered cascade of events leading to the establishment of heterochromatin and requiring the recruitment of the histone H2Av variant followed by H4 Lys 12 acetylation as necessary steps before H3 Lys 9 methylation and HP1 recruitment can take place. [Keywords: Chromatin; silencing; transcription; histone; nucleus] Received September 8, 2004; revised version accepted November 4, 2004. The basic unit of chromatin is the nucleosome, which is guchi et al. 2004). The role of histone variants, and spe- made up of 146 bp of DNA wrapped around a histone cially those of H3 and H2A, in various nuclear processes octamer composed of two molecules each of the histones has been long appreciated (Wolffe and Pruss 1996; Ah- H2A, H2B, H3, and H4.
    [Show full text]
  • Supplementary Table S4. FGA Co-Expressed Gene List in LUAD
    Supplementary Table S4. FGA co-expressed gene list in LUAD tumors Symbol R Locus Description FGG 0.919 4q28 fibrinogen gamma chain FGL1 0.635 8p22 fibrinogen-like 1 SLC7A2 0.536 8p22 solute carrier family 7 (cationic amino acid transporter, y+ system), member 2 DUSP4 0.521 8p12-p11 dual specificity phosphatase 4 HAL 0.51 12q22-q24.1histidine ammonia-lyase PDE4D 0.499 5q12 phosphodiesterase 4D, cAMP-specific FURIN 0.497 15q26.1 furin (paired basic amino acid cleaving enzyme) CPS1 0.49 2q35 carbamoyl-phosphate synthase 1, mitochondrial TESC 0.478 12q24.22 tescalcin INHA 0.465 2q35 inhibin, alpha S100P 0.461 4p16 S100 calcium binding protein P VPS37A 0.447 8p22 vacuolar protein sorting 37 homolog A (S. cerevisiae) SLC16A14 0.447 2q36.3 solute carrier family 16, member 14 PPARGC1A 0.443 4p15.1 peroxisome proliferator-activated receptor gamma, coactivator 1 alpha SIK1 0.435 21q22.3 salt-inducible kinase 1 IRS2 0.434 13q34 insulin receptor substrate 2 RND1 0.433 12q12 Rho family GTPase 1 HGD 0.433 3q13.33 homogentisate 1,2-dioxygenase PTP4A1 0.432 6q12 protein tyrosine phosphatase type IVA, member 1 C8orf4 0.428 8p11.2 chromosome 8 open reading frame 4 DDC 0.427 7p12.2 dopa decarboxylase (aromatic L-amino acid decarboxylase) TACC2 0.427 10q26 transforming, acidic coiled-coil containing protein 2 MUC13 0.422 3q21.2 mucin 13, cell surface associated C5 0.412 9q33-q34 complement component 5 NR4A2 0.412 2q22-q23 nuclear receptor subfamily 4, group A, member 2 EYS 0.411 6q12 eyes shut homolog (Drosophila) GPX2 0.406 14q24.1 glutathione peroxidase
    [Show full text]
  • Generate Metabolic Map Poster
    Authors: Pallavi Subhraveti Ron Caspi Peter Midford Peter D Karp An online version of this diagram is available at BioCyc.org. Biosynthetic pathways are positioned in the left of the cytoplasm, degradative pathways on the right, and reactions not assigned to any pathway are in the far right of the cytoplasm. Transporters and membrane proteins are shown on the membrane. Ingrid Keseler Periplasmic (where appropriate) and extracellular reactions and proteins may also be shown. Pathways are colored according to their cellular function. Gcf_001591825Cyc: Bacillus vietnamensis NBRC 101237 Cellular Overview Connections between pathways are omitted for legibility. Anamika Kothari sn-glycerol phosphate phosphate pro phosphate phosphate phosphate thiamine molybdate D-xylose D-ribose glutathione 3-phosphate D-mannitol L-cystine L-djenkolate lanthionine α,β-trehalose phosphate phosphate [+ 3 more] α,α-trehalose predicted predicted ABC ABC FliY ThiT XylF RbsB RS10935 UgpC TreP PutP RS10200 PstB PstB RS10385 RS03335 RS20030 RS19075 transporter transporter of molybdate of phosphate α,β-trehalose 6-phosphate L-cystine D-xylose D-ribose sn-glycerol D-mannitol phosphate phosphate thiamine glutathione α α phosphate phosphate phosphate phosphate L-djenkolate 3-phosphate , -trehalose 6-phosphate pro 1-phosphate lanthionine molybdate phosphate [+ 3 more] Metabolic Regulator Amino Acid Degradation Amine and Polyamine Biosynthesis Macromolecule Modification tRNA-uridine 2-thiolation Degradation ATP biosynthesis a mature peptidoglycan a nascent β an N-terminal-
    [Show full text]
  • The Neurochemical Consequences of Aromatic L-Amino Acid Decarboxylase Deficiency
    The neurochemical consequences of aromatic L-amino acid decarboxylase deficiency Submitted By: George Francis Gray Allen Department of Molecular Neuroscience UCL Institute of Neurology Queen Square, London Submitted November 2010 Funded by the AADC Research Trust, UK Thesis submitted for the degree of Doctor of Philosophy, University College London (UCL) 1 I, George Allen confirm that the work presented in this thesis is my own. Where information has been derived from other sources, I confirm that this has been indicated in the thesis. Signed………………………………………………….Date…………………………… 2 Abstract Aromatic L-amino acid decarboxylase (AADC) catalyses the conversion of 5- hydroxytryptophan (5-HTP) and L-3,4-dihydroxyphenylalanine (L-dopa) to the neurotransmitters serotonin and dopamine respectively. The inherited disorder AADC deficiency leads to a severe deficit of serotonin and dopamine as well as an accumulation of 5-HTP and L-dopa. This thesis investigated the potential role of 5- HTP/L-dopa accumulation in the pathogenesis of AADC deficiency. Treatment of human neuroblastoma cells with L-dopa or dopamine was found to increase intracellular levels of the antioxidant reduced glutathione (GSH). However inhibiting AADC prevented the GSH increase induced by L-dopa. Furthermore dopamine but not L-dopa, increased GSH release from human astrocytoma cells, which do not express AADC activity. GSH release is the first stage of GSH trafficking from astrocytes to neurons. This data indicates dopamine may play a role in controlling brain GSH levels and consequently antioxidant status. The inability of L-dopa to influence GSH concentrations in the absence of AADC or with AADC inhibited indicates GSH trafficking/metabolism may be compromised in AADC deficiency.
    [Show full text]
  • Histone H4 Lysine 20 Mono-Methylation Directly Facilitates Chromatin Openness and Promotes Transcription of Housekeeping Genes
    ARTICLE https://doi.org/10.1038/s41467-021-25051-2 OPEN Histone H4 lysine 20 mono-methylation directly facilitates chromatin openness and promotes transcription of housekeeping genes Muhammad Shoaib 1,8,9, Qinming Chen2,9, Xiangyan Shi 3, Nidhi Nair1, Chinmayi Prasanna 2, Renliang Yang2,4, David Walter1, Klaus S. Frederiksen 5, Hjorleifur Einarsson1, J. Peter Svensson 6, ✉ ✉ ✉ Chuan Fa Liu 2, Karl Ekwall6, Mads Lerdrup 7 , Lars Nordenskiöld 2 & Claus S. Sørensen 1 1234567890():,; Histone lysine methylations have primarily been linked to selective recruitment of reader or effector proteins that subsequently modify chromatin regions and mediate genome functions. Here, we describe a divergent role for histone H4 lysine 20 mono-methylation (H4K20me1) and demonstrate that it directly facilitates chromatin openness and accessibility by disrupting chromatin folding. Thus, accumulation of H4K20me1 demarcates highly accessible chromatin at genes, and this is maintained throughout the cell cycle. In vitro, H4K20me1-containing nucleosomal arrays with nucleosome repeat lengths (NRL) of 187 and 197 are less compact than unmethylated (H4K20me0) or trimethylated (H4K20me3) arrays. Concordantly, and in contrast to trimethylated and unmethylated tails, solid-state NMR data shows that H4K20 mono-methylation changes the H4 conformational state and leads to more dynamic histone H4-tails. Notably, the increased chromatin accessibility mediated by H4K20me1 facilitates gene expression, particularly of housekeeping genes. Altogether, we show how the methy- lation state of a single histone H4 residue operates as a focal point in chromatin structure control. While H4K20me1 directly promotes chromatin openness at highly transcribed genes, it also serves as a stepping-stone for H4K20me3-dependent chromatin compaction.
    [Show full text]
  • (P -Value<0.05, Fold Change≥1.4), 4 Vs. 0 Gy Irradiation
    Table S1: Significant differentially expressed genes (P -Value<0.05, Fold Change≥1.4), 4 vs. 0 Gy irradiation Genbank Fold Change P -Value Gene Symbol Description Accession Q9F8M7_CARHY (Q9F8M7) DTDP-glucose 4,6-dehydratase (Fragment), partial (9%) 6.70 0.017399678 THC2699065 [THC2719287] 5.53 0.003379195 BC013657 BC013657 Homo sapiens cDNA clone IMAGE:4152983, partial cds. [BC013657] 5.10 0.024641735 THC2750781 Ciliary dynein heavy chain 5 (Axonemal beta dynein heavy chain 5) (HL1). 4.07 0.04353262 DNAH5 [Source:Uniprot/SWISSPROT;Acc:Q8TE73] [ENST00000382416] 3.81 0.002855909 NM_145263 SPATA18 Homo sapiens spermatogenesis associated 18 homolog (rat) (SPATA18), mRNA [NM_145263] AA418814 zw01a02.s1 Soares_NhHMPu_S1 Homo sapiens cDNA clone IMAGE:767978 3', 3.69 0.03203913 AA418814 AA418814 mRNA sequence [AA418814] AL356953 leucine-rich repeat-containing G protein-coupled receptor 6 {Homo sapiens} (exp=0; 3.63 0.0277936 THC2705989 wgp=1; cg=0), partial (4%) [THC2752981] AA484677 ne64a07.s1 NCI_CGAP_Alv1 Homo sapiens cDNA clone IMAGE:909012, mRNA 3.63 0.027098073 AA484677 AA484677 sequence [AA484677] oe06h09.s1 NCI_CGAP_Ov2 Homo sapiens cDNA clone IMAGE:1385153, mRNA sequence 3.48 0.04468495 AA837799 AA837799 [AA837799] Homo sapiens hypothetical protein LOC340109, mRNA (cDNA clone IMAGE:5578073), partial 3.27 0.031178378 BC039509 LOC643401 cds. [BC039509] Homo sapiens Fas (TNF receptor superfamily, member 6) (FAS), transcript variant 1, mRNA 3.24 0.022156298 NM_000043 FAS [NM_000043] 3.20 0.021043295 A_32_P125056 BF803942 CM2-CI0135-021100-477-g08 CI0135 Homo sapiens cDNA, mRNA sequence 3.04 0.043389246 BF803942 BF803942 [BF803942] 3.03 0.002430239 NM_015920 RPS27L Homo sapiens ribosomal protein S27-like (RPS27L), mRNA [NM_015920] Homo sapiens tumor necrosis factor receptor superfamily, member 10c, decoy without an 2.98 0.021202829 NM_003841 TNFRSF10C intracellular domain (TNFRSF10C), mRNA [NM_003841] 2.97 0.03243901 AB002384 C6orf32 Homo sapiens mRNA for KIAA0386 gene, partial cds.
    [Show full text]
  • Histones H3 and H4 Are Components of Upstream Activation Factor Required for the High-Level Transcription of Yeast Rdna by RNA Polymerase I
    Proc. Natl. Acad. Sci. USA Vol. 94, pp. 13458–13462, December 1997 Biochemistry Histones H3 and H4 are components of upstream activation factor required for the high-level transcription of yeast rDNA by RNA polymerase I JOHN KEENER*, JONATHAN A. DODD*, DOMINIQUE LALO, AND MASAYASU NOMURA† Department of Biological Chemistry, University of California, Irvine, CA 92697-1700 Contributed by Masayasu Nomura, October 16, 1997 ABSTRACT RNA polymerase I (Pol I) transcription in consisting of Rrn6p, Rrn7p, and Rrn11p; refs. 3, 9, and 10), the yeast Saccharomyces cerevisiae is greatly stimulated in vivo Rrn3p (5), and Pol I are required. In addition to these factors, and in vitro by the multiprotein complex, upstream activation upstream activation factor (UAF) and TATA box-binding factor (UAF). UAF binds tightly to the upstream element of the protein, as well as the upstream element, are required for a rDNA promoter, such that once bound (in vitro), UAF does not high level of transcription from the yeast rDNA promoter (4, readily exchange onto a competing template. Of the polypep- 11). Purified UAF previously was shown to contain three tides previously identified in purified UAF, three are encoded genetically defined subunits, Rrn5p, Rrn9p, and Rrn10p, and by genes required for Pol I transcription in vivo: RRN5, RRN9, two uncharacterized subunits, p30 and p18 (4). and RRN10. Two others, p30 and p18, have remained unchar- DNA in the eukaryotic nucleus is organized as chromatin, acterized. We report here that the N-terminal amino acid consisting mostly of regularly repeating nucleosomes in which sequence, its mobility in gel electrophoresis, and the immu- DNA is wrapped around an octameric structure of core noreactivity of p18 shows that it is histone H3.
    [Show full text]
  • The Histone Variant H3.3 Marks Active Chromatin by Replication-Independent Nucleosome Assembly
    Molecular Cell, Vol. 9, 1191–1200, June, 2002, Copyright 2002 by Cell Press The Histone Variant H3.3 Marks Active Chromatin by Replication-Independent Nucleosome Assembly Kami Ahmad and Steven Henikoff1 ment are not clear. A study in Tetrahymena concluded Fred Hutchinson Cancer Research Center that no protein difference between histone H3 variants 1100 Fairview Avenue North was required for replacement histone deposition and Seattle, Washington 98109 that expression of either variant outside of S phase ap- peared to be sufficient (Yu and Gorovsky, 1997). In con- trast, by examining the dynamics of histone proteins in Summary Drosophila nuclei we show that the major histone H3 and the replacement histone H3.3 proteins have distinct Two very similar H3 histones—differing at only four properties during in vivo chromatin assembly. Histone amino acid positions—are produced in Drosophila H3.3 participates in replication-independent (RI) nucleo- cells. Here we describe a mechanism of chromatin some assembly and is targeted to transcriptionally ac- regulation whereby the variant H3.3 is deposited at tive loci throughout the cell cycle. Transcription-coupled particular loci, including active rDNA arrays. While the deposition of H3.3-containing nucleosomes may be a major H3 is incorporated strictly during DNA replica- general mechanism for rapidly replacing permanently tion, amino acid changes toward H3.3 allow replica- modified nucleosomes and for heritably activating tion-independent (RI) deposition. In contrast to repli- genes. cation-coupled (RC) deposition, RI deposition does not require the N-terminal tail. H3.3 is the exclusive Results substrate for RI deposition, and its counterpart is the only substrate retained in yeast.
    [Show full text]
  • Loss of Arylformamidase with Reduced Thymidine Kinase Expression Leads to Impaired Glucose Tolerance Alison J
    © 2015. Published by The Company of Biologists Ltd | Biology Open (2015) 4, 1367-1375 doi:10.1242/bio.013342 RESEARCH ARTICLE Loss of arylformamidase with reduced thymidine kinase expression leads to impaired glucose tolerance Alison J. Hugill1,*, Michelle E. Stewart2,*, Marianne A. Yon1,*, Fay Probert2,3, I. Jane Cox3, Tertius A. Hough2, Cheryl L. Scudamore2, Liz Bentley1, Gary Wall1, Sara E. Wells2,* and Roger D. Cox1,*,‡ ABSTRACT steps, generating metabolites including kynurenic acid (KYNA), Tryptophan metabolites have been linked in observational studies quinolonic acid, xanthurenic acid (XA), and nicotinamide with type 2 diabetes, cognitive disorders, inflammation and immune (reviewed Stone and Darlington, 2002). system regulation. A rate-limiting enzyme in tryptophan conversion is It has been suggested that dysregulation of XA and kynurenine arylformamidase (Afmid), and a double knockout of this gene and (KYN) derivatives are involved in insulin resistance and insulin thymidine kinase (Tk) has been reported to cause renal failure and secretion (reviewed Oxenkrug, 2013). XA excretion has been abnormal immune system regulation. In order to further investigate reported to be elevated in the urine of patients with diabetes and in possible links between abnormal tryptophan catabolism and diabetes alloxan- and streptozotocin- induced diabetic rats (Hattori et al., and to examine the effect of single Afmid knockout, we have carried 1984; Ikeda and Kotake, 1986). Similarly, a metabolomic analysis out metabolic phenotyping of an exon 2 Afmid gene knockout. These in rhesus macaques, comparing normal and spontaneously diabetic mice exhibit impaired glucose tolerance, although their insulin animals, identified elevated KYNA in urine as a biomarker of sensitivity is unchanged in comparison to wild-type animals.
    [Show full text]