1
|
Tremblay BJM, Santini CP, Cheng Y, Zhang X, Rosa S, Qüesta JI. Interplay between coding and non-coding regulation drives the Arabidopsis seed-to-seedling transition. Nat Commun 2024; 15:1724. [PMID: 38409232 PMCID: PMC10897432 DOI: 10.1038/s41467-024-46082-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 02/14/2024] [Indexed: 02/28/2024] Open
Abstract
Translation of seed stored mRNAs is essential to trigger germination. However, when RNAPII re-engages RNA synthesis during the seed-to-seedling transition has remained in question. Combining csRNA-seq, ATAC-seq and smFISH in Arabidopsis thaliana we demonstrate that active transcription initiation is detectable during the entire germination process. Features of non-coding regulation such as dynamic changes in chromatin accessible regions, antisense transcription, as well as bidirectional non-coding promoters are widespread throughout the Arabidopsis genome. We show that sensitivity to exogenous ABSCISIC ACID (ABA) during germination depends on proximal promoter accessibility at ABA-responsive genes. Moreover, we provide genetic validation of the existence of divergent transcription in plants. Our results reveal that active enhancer elements are transcribed producing non-coding enhancer RNAs (eRNAs) as widely documented in metazoans. In sum, this study defining the extent and role of coding and non-coding transcription during key stages of germination expands our understanding of transcriptional mechanisms underlying plant developmental transitions.
Collapse
Affiliation(s)
- Benjamin J M Tremblay
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus UAB, Bellaterra, Barcelona, Spain
| | - Cristina P Santini
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus UAB, Bellaterra, Barcelona, Spain
| | - Yajiao Cheng
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus UAB, Bellaterra, Barcelona, Spain
| | - Xue Zhang
- Plant Biology Department, Swedish University of Agricultural Sciences (SLU), Uppsala, Sweden
| | - Stefanie Rosa
- Plant Biology Department, Swedish University of Agricultural Sciences (SLU), Uppsala, Sweden
| | - Julia I Qüesta
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus UAB, Bellaterra, Barcelona, Spain.
| |
Collapse
|
2
|
Malfait J, Wan J, Spicuglia S. Epromoters are new players in the regulatory landscape with potential pleiotropic roles. Bioessays 2023; 45:e2300012. [PMID: 37246247 DOI: 10.1002/bies.202300012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Revised: 05/11/2023] [Accepted: 05/15/2023] [Indexed: 05/30/2023]
Abstract
Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.
Collapse
Affiliation(s)
- Juliette Malfait
- Aix-Marseille University, Inserm, TAGC, UMR1090, Marseille, France
- Equipe Labélisée Ligue Contre le Cancer, LIGUE, Marseille, France
| | - Jing Wan
- Aix-Marseille University, Inserm, TAGC, UMR1090, Marseille, France
- Equipe Labélisée Ligue Contre le Cancer, LIGUE, Marseille, France
| | - Salvatore Spicuglia
- Aix-Marseille University, Inserm, TAGC, UMR1090, Marseille, France
- Equipe Labélisée Ligue Contre le Cancer, LIGUE, Marseille, France
| |
Collapse
|
3
|
Kravchuk EV, Ashniev GA, Gladkova MG, Orlov AV, Vasileva AV, Boldyreva AV, Burenin AG, Skirda AM, Nikitin PI, Orlova NN. Experimental Validation and Prediction of Super-Enhancers: Advances and Challenges. Cells 2023; 12:cells12081191. [PMID: 37190100 DOI: 10.3390/cells12081191] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 04/07/2023] [Accepted: 04/14/2023] [Indexed: 05/17/2023] Open
Abstract
Super-enhancers (SEs) are cis-regulatory elements of the human genome that have been widely discussed since the discovery and origin of the term. Super-enhancers have been shown to be strongly associated with the expression of genes crucial for cell differentiation, cell stability maintenance, and tumorigenesis. Our goal was to systematize research studies dedicated to the investigation of structure and functions of super-enhancers as well as to define further perspectives of the field in various applications, such as drug development and clinical use. We overviewed the fundamental studies which provided experimental data on various pathologies and their associations with particular super-enhancers. The analysis of mainstream approaches for SE search and prediction allowed us to accumulate existing data and propose directions for further algorithmic improvements of SEs' reliability levels and efficiency. Thus, here we provide the description of the most robust algorithms such as ROSE, imPROSE, and DEEPSEN and suggest their further use for various research and development tasks. The most promising research direction, which is based on topic and number of published studies, are cancer-associated super-enhancers and prospective SE-targeted therapy strategies, most of which are discussed in this review.
Collapse
Affiliation(s)
- Ekaterina V Kravchuk
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
- Faculty of Biology, Lomonosov Moscow State University, Leninskiye Gory, MSU, 1-12, 119991 Moscow, Russia
| | - German A Ashniev
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
- Faculty of Biology, Lomonosov Moscow State University, Leninskiye Gory, MSU, 1-12, 119991 Moscow, Russia
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, GSP-1, Leninskiye Gory, MSU, 1-73, 119234 Moscow, Russia
| | - Marina G Gladkova
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, GSP-1, Leninskiye Gory, MSU, 1-73, 119234 Moscow, Russia
| | - Alexey V Orlov
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
| | - Anastasiia V Vasileva
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
| | - Anna V Boldyreva
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
| | - Alexandr G Burenin
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
| | - Artemiy M Skirda
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
| | - Petr I Nikitin
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
| | - Natalia N Orlova
- Prokhorov General Physics Institute of the Russian Academy of Sciences, 38 Vavilov St., 119991 Moscow, Russia
| |
Collapse
|
4
|
Anreiter I, Tian YW, Soller M. The cap epitranscriptome: Early directions to a complex life as mRNA. Bioessays 2023; 45:e2200198. [PMID: 36529693 DOI: 10.1002/bies.202200198] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Revised: 12/07/2022] [Accepted: 12/08/2022] [Indexed: 12/23/2022]
Abstract
Animal, protist and viral messenger RNAs (mRNAs) are most prominently modified at the beginning by methylation of cap-adjacent nucleotides at the 2'-O-position of the ribose (cOMe) by dedicated cap methyltransferases (CMTrs). If the first nucleotide of an mRNA is an adenosine, PCIF1 can methylate at the N6 -position (m6 A), while internally the Mettl3/14 writer complex can methylate. These modifications are introduced co-transcriptionally to affect many aspects of gene expression including localisation to synapses and local translation. Of particular interest, transcription start sites of many genes are heterogeneous leading to sequence diversity at the beginning of mRNAs, which together with cOMe and m6 Am could constitute an extensive novel layer of gene expression control. Given the role of cOMe and m6 A in local gene expression at synapses and higher brain functions including learning and memory, such code could be implemented at the transcriptional level for lasting memories through local gene expression at synapses.
Collapse
Affiliation(s)
- Ina Anreiter
- Department of Biological Sciences, University of Toronto Scarborough, Toronto, Canada
| | - Yuan W Tian
- Birmingham Centre for Genome Biology, University of Birmingham, Birmingham, UK.,School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Birmingham, UK
| | - Matthias Soller
- Birmingham Centre for Genome Biology, University of Birmingham, Birmingham, UK.,School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Birmingham, UK
| |
Collapse
|
5
|
Hamamoto K, Umemura Y, Makino S, Fukaya T. Dynamic interplay between non-coding enhancer transcription and gene activity in development. Nat Commun 2023; 14:826. [PMID: 36805453 PMCID: PMC9941499 DOI: 10.1038/s41467-023-36485-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2022] [Accepted: 02/03/2023] [Indexed: 02/22/2023] Open
Abstract
Non-coding transcription at the intergenic regulatory regions is a prevalent feature of metazoan genomes, but its biological function remains uncertain. Here, we devise a live-imaging system that permits simultaneous visualization of gene activity along with intergenic non-coding transcription at single-cell resolution in Drosophila. Quantitative image analysis reveals that elongation of RNA polymerase II across the internal core region of enhancers leads to suppression of transcriptional bursting from linked genes. Super-resolution imaging and genome-editing analysis further demonstrate that enhancer transcription antagonizes molecular crowding of transcription factors, thereby interrupting the formation of a transcription hub at the gene locus. We also show that a certain class of developmental enhancers are structurally optimized to co-activate gene transcription together with non-coding transcription effectively. We suggest that enhancer function is flexibly tunable through the modulation of hub formation via surrounding non-coding transcription during development.
Collapse
Affiliation(s)
- Kota Hamamoto
- Laboratory of Transcription Dynamics, Research Center for Biological Visualization, Institute for Quantitative Biosciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan.,Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Yusuke Umemura
- Laboratory of Transcription Dynamics, Research Center for Biological Visualization, Institute for Quantitative Biosciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan.,Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Shiho Makino
- Laboratory of Transcription Dynamics, Research Center for Biological Visualization, Institute for Quantitative Biosciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Takashi Fukaya
- Laboratory of Transcription Dynamics, Research Center for Biological Visualization, Institute for Quantitative Biosciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan. .,Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan.
| |
Collapse
|
6
|
Current challenges in understanding the role of enhancers in disease. Nat Struct Mol Biol 2022; 29:1148-1158. [PMID: 36482255 DOI: 10.1038/s41594-022-00896-3] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 11/04/2022] [Indexed: 12/13/2022]
Abstract
Enhancers play a central role in the spatiotemporal control of gene expression and tend to work in a cell-type-specific manner. In addition, they are suggested to be major contributors to phenotypic variation, evolution and disease. There is growing evidence that enhancer dysfunction due to genetic, structural or epigenetic mechanisms contributes to a broad range of human diseases referred to as enhanceropathies. Such mechanisms often underlie the susceptibility to common diseases, but can also play a direct causal role in cancer or Mendelian diseases. Despite the recent gain of insights into enhancer biology and function, we still have a limited ability to predict how enhancer dysfunction impacts gene expression. Here we discuss the major challenges that need to be overcome when studying the role of enhancers in disease etiology and highlight opportunities and directions for future studies, aiming to disentangle the molecular basis of enhanceropathies.
Collapse
|
7
|
van den Elzen AMG, Watson MJ, Thoreen CC. mRNA 5' terminal sequences drive 200-fold differences in expression through effects on synthesis, translation and decay. PLoS Genet 2022; 18:e1010532. [PMID: 36441824 PMCID: PMC9731452 DOI: 10.1371/journal.pgen.1010532] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Revised: 12/08/2022] [Accepted: 11/15/2022] [Indexed: 11/30/2022] Open
Abstract
mRNA regulatory sequences control gene expression at multiple levels including translation initiation and mRNA decay. The 5' terminal sequences of mRNAs have unique regulatory potential because of their proximity to key post-transcriptional regulators. Here we have systematically probed the function of 5' terminal sequences in gene expression in human cells. Using a library of reporter mRNAs initiating with all possible 7-mer sequences at their 5' ends, we find an unexpected impact on transcription that underlies 200-fold differences in mRNA expression. Library sequences that promote high levels of transcription mirrored those found in native mRNAs and define two basic classes with similarities to classic Initiator (Inr) and TCT core promoter motifs. By comparing transcription, translation and decay rates, we identify sequences that are optimized for both efficient transcription and growth-regulated translation and stability, including variants of terminal oligopyrimidine (TOP) motifs. We further show that 5' sequences of endogenous mRNAs are enriched for multi-functional TCT/TOP hybrid sequences. Together, our results reveal how 5' sequences define two general classes of mRNAs with distinct growth-responsive profiles of expression across synthesis, translation and decay.
Collapse
Affiliation(s)
- Antonia M. G. van den Elzen
- Department of Cellular and Molecular Physiology, Yale School of Medicine, New Haven, Connecticut, United States of America
| | - Maegan J. Watson
- Department of Cellular and Molecular Physiology, Yale School of Medicine, New Haven, Connecticut, United States of America
| | - Carson C. Thoreen
- Department of Cellular and Molecular Physiology, Yale School of Medicine, New Haven, Connecticut, United States of America
- * E-mail:
| |
Collapse
|
8
|
CMTr cap-adjacent 2'-O-ribose mRNA methyltransferases are required for reward learning and mRNA localization to synapses. Nat Commun 2022; 13:1209. [PMID: 35260552 PMCID: PMC8904806 DOI: 10.1038/s41467-022-28549-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Accepted: 01/18/2022] [Indexed: 11/25/2022] Open
Abstract
Cap-adjacent nucleotides of animal, protist and viral mRNAs can be O-methylated at the 2‘ position of the ribose (cOMe). The functions of cOMe in animals, however, remain largely unknown. Here we show that the two cap methyltransferases (CMTr1 and CMTr2) of Drosophila can methylate the ribose of the first nucleotide in mRNA. Double-mutant flies lack cOMe but are viable. Consistent with prominent neuronal expression, they have a reward learning defect that can be rescued by conditional expression in mushroom body neurons before training. Among CMTr targets are cell adhesion and signaling molecules. Many are relevant for learning, and are also targets of Fragile X Mental Retardation Protein (FMRP). Like FMRP, cOMe is required for localization of untranslated mRNAs to synapses and enhances binding of the cap binding complex in the nucleus. Hence, our study reveals a mechanism to co-transcriptionally prime mRNAs by cOMe for localized protein synthesis at synapses. The two cap methyltransferases (CMTrs) redundantly methylate riboses of first cap adjacent nucleotides in messenger RNAs in Drosophila. Here, CMTrs are required for reward learning and localization of untranslated messenger RNAs to synapses.
Collapse
|
9
|
Santiago-Algarra D, Souaid C, Singh H, Dao LTM, Hussain S, Medina-Rivera A, Ramirez-Navarro L, Castro-Mondragon JA, Sadouni N, Charbonnier G, Spicuglia S. Epromoters function as a hub to recruit key transcription factors required for the inflammatory response. Nat Commun 2021; 12:6660. [PMID: 34795220 PMCID: PMC8602369 DOI: 10.1038/s41467-021-26861-0] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Accepted: 10/14/2021] [Indexed: 12/14/2022] Open
Abstract
Gene expression is controlled by the involvement of gene-proximal (promoters) and distal (enhancers) regulatory elements. Our previous results demonstrated that a subset of gene promoters, termed Epromoters, work as bona fide enhancers and regulate distal gene expression. Here, we hypothesized that Epromoters play a key role in the coordination of rapid gene induction during the inflammatory response. Using a high-throughput reporter assay we explored the function of Epromoters in response to type I interferon. We find that clusters of IFNa-induced genes are frequently associated with Epromoters and that these regulatory elements preferentially recruit the STAT1/2 and IRF transcription factors and distally regulate the activation of interferon-response genes. Consistently, we identified and validated the involvement of Epromoter-containing clusters in the regulation of LPS-stimulated macrophages. Our findings suggest that Epromoters function as a local hub recruiting the key TFs required for coordinated regulation of gene clusters during the inflammatory response.
Collapse
Affiliation(s)
- David Santiago-Algarra
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Equipe Labellisée Ligue Contre le Cancer, Paris, France
| | - Charbel Souaid
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Equipe Labellisée Ligue Contre le Cancer, Paris, France
| | - Himanshu Singh
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Equipe Labellisée Ligue Contre le Cancer, Paris, France
| | - Lan T M Dao
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Equipe Labellisée Ligue Contre le Cancer, Paris, France
- Vinmec Research Institute of Stem cell and Gene technology, Vinmec Healthcare System, Hanoi, Vietnam
| | - Saadat Hussain
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Equipe Labellisée Ligue Contre le Cancer, Paris, France
| | - Alejandra Medina-Rivera
- Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Juriquilla, Mexico
| | - Lucia Ramirez-Navarro
- Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Juriquilla, Mexico
| | - Jaime A Castro-Mondragon
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, 0318, Oslo, Norway
| | - Nori Sadouni
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Equipe Labellisée Ligue Contre le Cancer, Paris, France
| | - Guillaume Charbonnier
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France
- Equipe Labellisée Ligue Contre le Cancer, Paris, France
| | - Salvatore Spicuglia
- Aix-Marseille University, INSERM, TAGC, UMR 1090, Marseille, France.
- Equipe Labellisée Ligue Contre le Cancer, Paris, France.
| |
Collapse
|
10
|
Wolfe JC, Mikheeva LA, Hagras H, Zabet NR. An explainable artificial intelligence approach for decoding the enhancer histone modifications code and identification of novel enhancers in Drosophila. Genome Biol 2021; 22:308. [PMID: 34749786 PMCID: PMC8574042 DOI: 10.1186/s13059-021-02532-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Accepted: 10/29/2021] [Indexed: 12/16/2022] Open
Abstract
BACKGROUND Enhancers are non-coding regions of the genome that control the activity of target genes. Recent efforts to identify active enhancers experimentally and in silico have proven effective. While these tools can predict the locations of enhancers with a high degree of accuracy, the mechanisms underpinning the activity of enhancers are often unclear. RESULTS Using machine learning (ML) and a rule-based explainable artificial intelligence (XAI) model, we demonstrate that we can predict the location of known enhancers in Drosophila with a high degree of accuracy. Most importantly, we use the rules of the XAI model to provide insight into the underlying combinatorial histone modifications code of enhancers. In addition, we identified a large set of putative enhancers that display the same epigenetic signature as enhancers identified experimentally. These putative enhancers are enriched in nascent transcription, divergent transcription and have 3D contacts with promoters of transcribed genes. However, they display only intermediary enrichment of mediator and cohesin complexes compared to previously characterised active enhancers. We also found that 10-15% of the predicted enhancers display similar characteristics to super enhancers observed in other species. CONCLUSIONS Here, we applied an explainable AI model to predict enhancers with high accuracy. Most importantly, we identified that different combinations of epigenetic marks characterise different groups of enhancers. Finally, we discovered a large set of putative enhancers which display similar characteristics with previously characterised active enhancers.
Collapse
Affiliation(s)
- Jareth C Wolfe
- School of Life Sciences, University of Essex, Colchester, CO4 3SQ, UK
- School of Computer Science and Electronic Engineering, University of Essex, Colchester, CO4 3SQ, UK
- Blizard Institute, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, E1 2AT, London, UK
| | - Liudmila A Mikheeva
- School of Life Sciences, University of Essex, Colchester, CO4 3SQ, UK
- Blizard Institute, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, E1 2AT, London, UK
- Department of Mathematical Sciences, University of Essex, Colchester, CO4 3SQ, UK
| | - Hani Hagras
- School of Computer Science and Electronic Engineering, University of Essex, Colchester, CO4 3SQ, UK.
| | - Nicolae Radu Zabet
- School of Life Sciences, University of Essex, Colchester, CO4 3SQ, UK.
- Blizard Institute, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, E1 2AT, London, UK.
| |
Collapse
|
11
|
Lozano R, Booth GT, Omar BY, Li B, Buckler ES, Lis JT, Del Carpio DP, Jannink JL. RNA polymerase mapping in plants identifies intergenic regulatory elements enriched in causal variants. G3-GENES GENOMES GENETICS 2021; 11:6364897. [PMID: 34499719 PMCID: PMC8527479 DOI: 10.1093/g3journal/jkab273] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Accepted: 07/04/2021] [Indexed: 12/14/2022]
Abstract
Control of gene expression is fundamental at every level of cell function. Promoter-proximal pausing and divergent transcription at promoters and enhancers, which are prominent features in animals, have only been studied in a handful of research experiments in plants. PRO-Seq analysis in cassava (Manihot esculenta) identified peaks of transcriptionally engaged RNA polymerase at both the 5' and 3' end of genes, consistent with paused or slowly moving Polymerase. In addition, we identified divergent transcription at intergenic sites. A full genome search for bi-directional transcription using an algorithm for enhancer detection developed in mammals (dREG) identified many intergenic regulatory element (IRE) candidates. These sites showed distinct patterns of methylation and nucleotide conservation based on genomic evolutionary rate profiling (GERP). SNPs within these IRE candidates explained significantly more variation in fitness and root composition than SNPs in chromosomal segments randomly ascertained from the same intergenic distribution, strongly suggesting a functional importance of these sites. Maize GRO-Seq data showed RNA polymerase occupancy at IREs consistent with patterns in cassava. Furthermore, these IREs in maize significantly overlapped with sites previously identified on the basis of open chromatin, histone marks, and methylation, and were enriched for reported eQTL. Our results suggest that bidirectional transcription can identify intergenic genomic regions in plants that play an important role in transcription regulation and whose identification has the potential to aid crop improvement.
Collapse
Affiliation(s)
- Roberto Lozano
- Plant Breeding and Genetics, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA
| | - Gregory T Booth
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY 14853, USA
| | | | - Bo Li
- State Key Laboratory of Plant Genomics and National Center for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Science, Beijing 100101, China
| | - Edward S Buckler
- Plant Breeding and Genetics, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA.,Institute for Genomic Diversity, Cornell University, Ithaca, NY 14853, USA.,United States Department of Agriculture, Agricultural Research Service (USDA-ARS) R.W. Holley Center for Agriculture and Health, Ithaca, NY 14853, USA
| | - John T Lis
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY 14853, USA
| | - Dunia Pino Del Carpio
- Plant Breeding and Genetics, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA
| | - Jean-Luc Jannink
- Plant Breeding and Genetics, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA.,United States Department of Agriculture, Agricultural Research Service (USDA-ARS) R.W. Holley Center for Agriculture and Health, Ithaca, NY 14853, USA
| |
Collapse
|
12
|
Neville MC, Eastwood A, Allen AM, de Haan A, Nojima T, Goodwin SF. Generation and characterization of fruitless P1 promoter mutant in Drosophila melanogaster. J Neurogenet 2021; 35:285-294. [PMID: 34338589 PMCID: PMC8477730 DOI: 10.1080/01677063.2021.1931179] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
The identification of mutations in the gene fruitless (fru) paved the way for understanding the genetic basis of male sexual behavior in the vinegar fly Drosophila melanogaster. D. melanogaster males perform an elaborate courtship display to the female, ultimately leading to copulation. Mutations in fru have been shown to disrupt most aspects of the male's behavioral display, rendering males behaviorally sterile. The fru genomic locus encodes for multiple transcription factor isoforms from several promoters; only those under the regulation of the most distal P1 promoter are under the control of the sex determination hierarchy and play a role in male-specific behaviors. In this study, we used CRISPR/Cas9-based targeted genome editing of the fru gene, to remove the P1 promoter region. We have shown that removal of the P1 promoter leads to a dramatic decrease in male courtship displays towards females and male-specific sterility. We have expanded the analysis of fru P1-dependent behaviors, examining male's response to courtship song and general activity levels during12-hour light: dark cycles. Our novel allele expands the mutant repertoire available for future studies of fru P1-derived function in D. melanogaster. Our fruΔP1 mutant will be useful for future studies of fru P1-derived function, as it can be homozygosed without disrupting additional downstream promoter function and can be utilized in heterozygous combinations with other extant fru alleles.
Collapse
Affiliation(s)
- Megan C. Neville
- University of Oxford, Centre for Neural Circuits & Behaviour, Oxford, UK,CONTACT Megan C. Neville University of Oxford, Centre for Neural Circuits & Behaviour, Oxford, UK
| | - Alexander Eastwood
- University of Oxford, Centre for Neural Circuits & Behaviour, Oxford, UK
| | - Aaron M. Allen
- University of Oxford, Centre for Neural Circuits & Behaviour, Oxford, UK
| | - Ammerins de Haan
- University of Oxford, Centre for Neural Circuits & Behaviour, Oxford, UK
| | - Tetsuya Nojima
- University of Oxford, Centre for Neural Circuits & Behaviour, Oxford, UK
| | - Stephen F. Goodwin
- University of Oxford, Centre for Neural Circuits & Behaviour, Oxford, UK
| |
Collapse
|
13
|
Wang R, Tang Q. Current Advances on the Important Roles of Enhancer RNAs in Molecular Pathways of Cancer. Int J Mol Sci 2021; 22:5640. [PMID: 34073237 PMCID: PMC8198447 DOI: 10.3390/ijms22115640] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2021] [Revised: 05/22/2021] [Accepted: 05/23/2021] [Indexed: 01/08/2023] Open
Abstract
Enhancers are critical genomic elements that can cooperate with promoters to regulate gene transcription in both normal and cancer cells. Recent studies reveal that enhancer regions are transcribed to produce a class of noncoding RNAs referred to as enhancer RNAs (eRNAs). Emerging evidence shows that eRNAs play important roles in enhancer activation and enhancer-driven gene regulation, and the expression of eRNAs may be a critical factor in tumorigenesis. The important roles of eRNAs in cancer signaling pathways are also gradually unveiled, providing a new insight into cancer therapy. Here, we review the roles of eRNAs in regulating cancer signaling pathways and discuss the potential of eRNA-targeted therapy for human cancers.
Collapse
Affiliation(s)
| | - Qianzi Tang
- Institute of Animal Genetics and Breeding, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China;
| |
Collapse
|
14
|
Steinhaus R, Gonzalez T, Seelow D, Robinson PN. Pervasive and CpG-dependent promoter-like characteristics of transcribed enhancers. Nucleic Acids Res 2020; 48:5306-5317. [PMID: 32338759 PMCID: PMC7261191 DOI: 10.1093/nar/gkaa223] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2019] [Revised: 03/23/2020] [Accepted: 03/25/2020] [Indexed: 12/17/2022] Open
Abstract
The temporal and spatial expression of genes is controlled by promoters and enhancers. Findings obtained over the last decade that not only promoters but also enhancers are characterized by bidirectional, divergent transcription have challenged the traditional notion that promoters and enhancers represent distinct classes of regulatory elements. Over half of human promoters are associated with CpG islands (CGIs), relatively CpG-rich stretches of generally several hundred nucleotides that are often associated with housekeeping genes. Only about 6% of transcribed enhancers defined by CAGE-tag analysis are associated with CGIs. Here, we present an analysis of enhancer and promoter characteristics and relate them to the presence or absence of CGIs. We show that transcribed enhancers share a number of CGI-dependent characteristics with promoters, including statistically significant local overrepresentation of core promoter elements. CGI-associated enhancers are longer, display higher directionality of transcription, greater expression, a lesser degree of tissue specificity, and a higher frequency of transcription-factor binding events than non-CGI-associated enhancers. Genes putatively regulated by CGI-associated enhancers are enriched for transcription regulator activity. Our findings show that CGI-associated transcribed enhancers display a series of characteristics related to sequence, expression and function that distinguish them from enhancers not associated with CGIs.
Collapse
Affiliation(s)
- Robin Steinhaus
- Berlin Institute of Health, Charitéplatz 1, 10117 Berlin, Germany.,Institute of Medical Genetics and Human Genetics, Charité - Universitätsmedizin Berlin, Augustenburger Platz 1, 13353 Berlin, Germany
| | - Tonatiuh Gonzalez
- The Jackson Laboratory for Genomic Medicine, 10 Discovery Drive, Farmington, CT 06032, USA.,Harvey Mudd College, 301 Platt Boulevard, Claremont, CA 91711, USA
| | - Dominik Seelow
- Berlin Institute of Health, Charitéplatz 1, 10117 Berlin, Germany.,Institute of Medical Genetics and Human Genetics, Charité - Universitätsmedizin Berlin, Augustenburger Platz 1, 13353 Berlin, Germany
| | - Peter N Robinson
- The Jackson Laboratory for Genomic Medicine, 10 Discovery Drive, Farmington, CT 06032, USA.,Institute for Systems Genomics, University of Connecticut, 263 Farmington Avenue, Farmington, CT 06030, USA
| |
Collapse
|
15
|
Ibragimov AN, Bylino OV, Shidlovskii YV. Molecular Basis of the Function of Transcriptional Enhancers. Cells 2020; 9:E1620. [PMID: 32635644 PMCID: PMC7407508 DOI: 10.3390/cells9071620] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2020] [Revised: 07/03/2020] [Accepted: 07/03/2020] [Indexed: 02/06/2023] Open
Abstract
Transcriptional enhancers are major genomic elements that control gene activity in eukaryotes. Recent studies provided deeper insight into the temporal and spatial organization of transcription in the nucleus, the role of non-coding RNAs in the process, and the epigenetic control of gene expression. Thus, multiple molecular details of enhancer functioning were revealed. Here, we describe the recent data and models of molecular organization of enhancer-driven transcription.
Collapse
Affiliation(s)
- Airat N. Ibragimov
- Laboratory of Gene Expression Regulation in Development, Institute of Gene Biology, Russian Academy of Sciences, 34/5 Vavilov St., 119334 Moscow, Russia; (A.N.I.); (O.V.B.)
- Center for Precision Genome Editing and Genetic Technologies for Biomedicine, Institute of Gene Biology, Russian Academy of Sciences, 34/5 Vavilov St., 119334 Moscow, Russia
| | - Oleg V. Bylino
- Laboratory of Gene Expression Regulation in Development, Institute of Gene Biology, Russian Academy of Sciences, 34/5 Vavilov St., 119334 Moscow, Russia; (A.N.I.); (O.V.B.)
| | - Yulii V. Shidlovskii
- Laboratory of Gene Expression Regulation in Development, Institute of Gene Biology, Russian Academy of Sciences, 34/5 Vavilov St., 119334 Moscow, Russia; (A.N.I.); (O.V.B.)
- I.M. Sechenov First Moscow State Medical University, 8, bldg. 2 Trubetskaya St., 119048 Moscow, Russia
| |
Collapse
|
16
|
Abstract
Cell-type- and condition-specific profiles of gene expression require coordination between protein-coding gene promoters and cis-regulatory sequences called enhancers. Enhancers can stimulate gene activity at great genomic distances from their targets, raising questions about how enhancers communicate with specific gene promoters and what molecular mechanisms underlie enhancer function. Characterization of enhancer loci has identified the molecular features of active enhancers that accompany the binding of transcription factors and local opening of chromatin. These characteristics include coactivator recruitment, histone modifications, and noncoding RNA transcription. However, it remains unclear which of these features functionally contribute to enhancer activity. Here, we discuss what is known about how enhancers regulate their target genes and how enhancers and promoters communicate. Further, we describe recent data demonstrating many similarities between enhancers and the gene promoters they control, and we highlight unanswered questions in the field, such as the potential roles of transcription at enhancers.
Collapse
Affiliation(s)
- Andrew Field
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Karen Adelman
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, Massachusetts 02115, USA
| |
Collapse
|
17
|
Thieffry A, Vigh ML, Bornholdt J, Ivanov M, Brodersen P, Sandelin A. Characterization of Arabidopsis thaliana Promoter Bidirectionality and Antisense RNAs by Inactivation of Nuclear RNA Decay Pathways. THE PLANT CELL 2020; 32:1845-1867. [PMID: 32213639 PMCID: PMC7268790 DOI: 10.1105/tpc.19.00815] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Revised: 02/03/2020] [Accepted: 03/20/2020] [Indexed: 05/20/2023]
Abstract
In animals, RNA polymerase II initiates transcription bidirectionally from gene promoters to produce pre-mRNAs on the forward strand and promoter upstream transcripts (PROMPTs) on the reverse strand. PROMPTs are degraded by the nuclear exosome. Previous studies based on nascent RNA approaches concluded that Arabidopsis (Arabidopsis thaliana) does not produce PROMPTs. Here, we used steady-state RNA sequencing in mutants defective in nuclear RNA decay including the exosome to reassess the existence of Arabidopsis PROMPTs. While they are rare, we identified ∼100 cases of exosome-sensitive PROMPTs in Arabidopsis. Such PROMPTs are sources of small interfering RNAs in exosome-deficient mutants, perhaps explaining why plants have evolved mechanisms to suppress PROMPTs. In addition, we found ∼200 long, unspliced and exosome-sensitive antisense RNAs that arise from transcription start sites within parts of the genome encoding 3'-untranslated regions on the sense strand. The previously characterized noncoding RNA that regulates expression of the key seed dormancy regulator, DELAY OF GERMINATION1, is a typical representative of this class of RNAs. Transcription factor genes are overrepresented among loci with exosome-sensitive antisense RNAs, suggesting a potential for widespread control of gene expression via this class of noncoding RNAs. Lastly, we assess the use of alternative promoters in Arabidopsis and compare the accuracy of existing TSS annotations.
Collapse
Affiliation(s)
- Axel Thieffry
- Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
- Biotech Research and Innovation Centre, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Maria Louisa Vigh
- Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Jette Bornholdt
- Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
- Biotech Research and Innovation Centre, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Maxim Ivanov
- Department of Plant and Environmental Sciences, University of Copenhagen, DK-1871 Frederiksberg C, Denmark
| | - Peter Brodersen
- Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Albin Sandelin
- Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
- Biotech Research and Innovation Centre, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
18
|
Fitz J, Neumann T, Steininger M, Wiedemann EM, Garcia AC, Athanasiadis A, Schoeberl UE, Pavri R. Spt5-mediated enhancer transcription directly couples enhancer activation with physical promoter interaction. Nat Genet 2020; 52:505-515. [PMID: 32251373 DOI: 10.1038/s41588-020-0605-6] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2019] [Accepted: 02/28/2020] [Indexed: 12/19/2022]
Abstract
Active enhancers are frequently transcribed, yet the regulatory role of enhancer transcription remains debated. Here, we depleted the RNA polymerase II pausing and elongation factor Spt5 in activated mouse B cells and found that approximately 50% of enhancer-gene pairs showed co-regulated transcription, consistent with a potential functional requirement for enhancer transcription. In particular, Spt5 depletion led to loss of super-enhancer-promoter physical interaction and gene expression at the immunoglobulin heavy-chain locus (Igh), abrogating antibody class switch recombination. This defect correlated strictly with loss of enhancer transcription but did not affect acetylation of histone H3 at lysine 27, chromatin accessibility and occupancy of Mediator and cohesin at the enhancer. Strikingly, CRISPRa-mediated rescue of enhancer transcription in Spt5-depleted cells restored Igh gene expression. Our work suggests that Spt5-mediated enhancer transcription underlies the physical and functional interaction between a subset of active enhancers and their target promoters.
Collapse
Affiliation(s)
- Johanna Fitz
- Research Institute of Molecular Pathology (IMP), Vienna, Austria
| | - Tobias Neumann
- Research Institute of Molecular Pathology (IMP), Vienna, Austria
| | | | | | | | | | | | - Rushad Pavri
- Research Institute of Molecular Pathology (IMP), Vienna, Austria.
| |
Collapse
|
19
|
Developmental regulation of cell type-specific transcription by novel promoter-proximal sequence elements. Genes Dev 2020; 34:663-677. [PMID: 32217666 PMCID: PMC7197356 DOI: 10.1101/gad.335331.119] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2019] [Accepted: 03/04/2020] [Indexed: 02/06/2023]
Abstract
Here, Lu et al. use Drosophila to show that opening of promoters from their closed state in precursor cells requires function of the spermatocyte-specific tMAC complex, localized at the promoters. Their findings provide novel insight into how promoter-proximal sequence elements that recruit and are acted on by cell type-specific chromatin-binding complexes help establish a robust, cell type-specific transcription program for terminal differentiation. Cell type-specific transcriptional programs that drive differentiation of specialized cell types are key players in development and tissue regeneration. One of the most dramatic changes in the transcription program in Drosophila occurs with the transition from proliferating spermatogonia to differentiating spermatocytes, with >3000 genes either newly expressed or expressed from new alternative promoters in spermatocytes. Here we show that opening of these promoters from their closed state in precursor cells requires function of the spermatocyte-specific tMAC complex, localized at the promoters. The spermatocyte-specific promoters lack the previously identified canonical core promoter elements except for the Inr. Instead, these promoters are enriched for the binding site for the TALE-class homeodomain transcription factors Achi/Vis and for a motif originally identified under tMAC ChIP-seq peaks. The tMAC motif resembles part of the previously identified 14-bp β2UE1 element critical for spermatocyte-specific expression. Analysis of downstream sequences relative to transcription start site usage suggested that ACA and CNAAATT motifs at specific positions can help promote efficient transcription initiation. Our results reveal how promoter-proximal sequence elements that recruit and are acted upon by cell type-specific chromatin binding complexes help establish a robust, cell type-specific transcription program for terminal differentiation.
Collapse
|
20
|
Arnold PR, Wells AD, Li XC. Diversity and Emerging Roles of Enhancer RNA in Regulation of Gene Expression and Cell Fate. Front Cell Dev Biol 2020; 7:377. [PMID: 31993419 PMCID: PMC6971116 DOI: 10.3389/fcell.2019.00377] [Citation(s) in RCA: 120] [Impact Index Per Article: 30.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Accepted: 12/17/2019] [Indexed: 12/27/2022] Open
Abstract
Enhancers are cis-regulatory elements in the genome that cooperate with promoters to control target gene transcription. Unlike promoters, enhancers are not necessarily adjacent to target genes and can exert their functions regardless of enhancer orientations, positions and spatial segregations from target genes. Thus, for a long time, the question as to how enhancers act in a temporal and spatial manner attracted considerable attention. The recent discovery that enhancers are also abundantly transcribed raises interesting questions about the exact roles of enhancer RNA (eRNA) in gene regulation. In this review, we highlight the process of enhancer transcription and the diverse features of eRNA. We review eRNA functions, which include enhancer-promoter looping, chromatin modifying, and transcription regulating. As eRNA are transcribed from active enhancers, they exhibit tissue and lineage specificity, and serve as markers of cell state and function. Finally, we discuss the unique relationship between eRNA and super enhancers in phase separation wherein eRNA may contribute significantly to cell fate decisions.
Collapse
Affiliation(s)
- Preston R Arnold
- Texas A&M Health Science Center, College of Medicine, Bryan, TX, United States.,Immunobiology and Transplant Sciences, Department of Surgery, Houston Methodist Hospital, Houston, TX, United States
| | - Andrew D Wells
- Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States
| | - Xian C Li
- Immunobiology and Transplant Sciences, Department of Surgery, Houston Methodist Hospital, Houston, TX, United States
| |
Collapse
|
21
|
Cardiello JF, Sanchez GJ, Allen MA, Dowell RD. Lessons from eRNAs: understanding transcriptional regulation through the lens of nascent RNAs. Transcription 2019; 11:3-18. [PMID: 31856658 DOI: 10.1080/21541264.2019.1704128] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Nascent transcription assays, such as global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), have uncovered a myriad of unstable RNAs being actively produced from numerous sites genome-wide. These transcripts provide a more complete and immediate picture of the impact of regulatory events. Transcription factors recruit RNA polymerase II, effectively initiating the process of transcription; repressors inhibit polymerase recruitment. Efficiency of recruitment is dictated by sequence elements in and around the RNA polymerase loading zone. A combination of sequence elements and RNA binding proteins subsequently influence the ultimate stability of the resulting transcript. Some of these transcripts are capable of providing feedback on the process, influencing subsequent transcription. By monitoring RNA polymerase activity, nascent assays provide insights into every step of the regulated process of transcription.
Collapse
Affiliation(s)
| | - Gilson J Sanchez
- BioFrontiers Institute, University of Colorado, Boulder, CO, USA
| | - Mary A Allen
- BioFrontiers Institute, University of Colorado, Boulder, CO, USA
| | - Robin D Dowell
- BioFrontiers Institute, University of Colorado, Boulder, CO, USA.,Molecular, Cellular, and Developmental Biology, University of Colorado, Boulder, CO, USA
| |
Collapse
|
22
|
Determinants of enhancer and promoter activities of regulatory elements. Nat Rev Genet 2019; 21:71-87. [DOI: 10.1038/s41576-019-0173-8] [Citation(s) in RCA: 284] [Impact Index Per Article: 56.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/04/2019] [Indexed: 12/13/2022]
|
23
|
Esse R, Gushchanskaia ES, Lord A, Grishok A. DOT1L complex suppresses transcription from enhancer elements and ectopic RNAi in Caenorhabditis elegans. RNA (NEW YORK, N.Y.) 2019; 25:1259-1273. [PMID: 31300558 PMCID: PMC6800474 DOI: 10.1261/rna.070292.119] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/07/2019] [Accepted: 07/10/2019] [Indexed: 05/14/2023]
Abstract
Methylation of histone H3 on lysine 79 (H3K79) by DOT1L is associated with actively transcribed genes. Earlier, we described that DOT-1.1, the Caenorhabditis elegans homolog of mammalian DOT1L, cooperates with the chromatin-binding protein ZFP-1 (AF10 homolog) to negatively modulate transcription of highly and widely expressed target genes. Also, the reduction of ZFP-1 levels has consistently been associated with lower efficiency of RNA interference (RNAi) triggered by exogenous double-stranded RNA (dsRNA), but the reason for this is not clear. Here, we demonstrate that the DOT1L complex suppresses transcription originating from enhancer elements and antisense transcription, thus potentiating the expression of enhancer-regulated genes. We also show that worms lacking H3K79 methylation do not survive, and this lethality is suppressed by a loss of caspase-3 or Dicer complex components that initiate gene silencing response to exogenous dsRNA. Our results suggest that ectopic elevation of endogenous dsRNA directly or indirectly resulting from global misregulation of transcription in DOT1L complex mutants may engage the Dicer complex and, therefore, limit the efficiency of exogenous RNAi.
Collapse
Affiliation(s)
- Ruben Esse
- Department of Biochemistry, Boston University School of Medicine, Boston, Massachusetts 02118, USA
| | | | - Avery Lord
- Department of Biochemistry, Boston University School of Medicine, Boston, Massachusetts 02118, USA
| | - Alla Grishok
- Department of Biochemistry, Boston University School of Medicine, Boston, Massachusetts 02118, USA
- Genome Science Institute, Boston University School of Medicine, Boston, Massachusetts 02118, USA
| |
Collapse
|
24
|
Chereji RV, Bryson TD, Henikoff S. Quantitative MNase-seq accurately maps nucleosome occupancy levels. Genome Biol 2019; 20:198. [PMID: 31519205 PMCID: PMC6743174 DOI: 10.1186/s13059-019-1815-z] [Citation(s) in RCA: 74] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2019] [Accepted: 09/04/2019] [Indexed: 12/15/2022] Open
Abstract
Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, its aggressive endo-/exo-nuclease activities make MNase-seq unreliable for determining nucleosome occupancies, because cleavages within linker regions produce oligo- and mono-nucleosomes, whereas cleavages within nucleosomes destroy them. Here, we introduce a theoretical framework for predicting nucleosome occupancies and an experimental protocol with appropriate spike-in normalization that confirms our theory and provides accurate occupancy levels over an MNase digestion time course. As with human cells, we observe no overall differences in nucleosome occupancies between Drosophila euchromatin and heterochromatin, which implies that heterochromatic compaction does not reduce MNase accessibility of linker DNA.
Collapse
Affiliation(s)
- Răzvan V Chereji
- Division of Developmental Biology, Eunice Kennedy Shriver National Institute for Child Health and Human Development, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Terri D Bryson
- Howard Hughes Medical Institute and Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA, 98109, USA
| | - Steven Henikoff
- Howard Hughes Medical Institute and Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA, 98109, USA.
| |
Collapse
|
25
|
Bhardwaj V, Semplicio G, Erdogdu NU, Manke T, Akhtar A. MAPCap allows high-resolution detection and differential expression analysis of transcription start sites. Nat Commun 2019; 10:3219. [PMID: 31363093 PMCID: PMC6667505 DOI: 10.1038/s41467-019-11115-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2019] [Accepted: 06/20/2019] [Indexed: 01/06/2023] Open
Abstract
The position, shape and number of transcription start sites (TSS) are critical determinants of gene regulation. Most methods developed to detect TSSs and study promoter usage are, however, of limited use in studies that demand quantification of expression changes between two or more groups. In this study, we combine high-resolution detection of transcription start sites and differential expression analysis using a simplified TSS quantification protocol, MAPCap (Multiplexed Affinity Purification of Capped RNA) along with the software icetea. Applying MAPCap on developing Drosophila melanogaster embryos and larvae, we detected stage and sex-specific promoter and enhancer activity and quantify the effect of mutants of maleless (MLE) helicase at X-chromosomal promoters. We observe that MLE mutation leads to a median 1.9 fold drop in expression of X-chromosome promoters and affects the expression of several TSSs with a sexually dimorphic expression on autosomes. Our results provide quantitative insights into promoter activity during dosage compensation. The position, shape and number of transcription start sites (TSS) regulate gene expression. Here authors present MAPCap, a method for high-resolution detection and differential expression analysis of TSS, and apply MAPCap to early fly development, detecting stage and sex-specific promoter and enhancer activity.
Collapse
Affiliation(s)
- Vivek Bhardwaj
- Max Planck Institute for Immunobiology and Epigenetics, 79108, Freiburg, Germany.,Faculty of Biology, University of Freiburg, 79104, Freiburg, Germany
| | - Giuseppe Semplicio
- Max Planck Institute for Immunobiology and Epigenetics, 79108, Freiburg, Germany
| | - Niyazi Umut Erdogdu
- Max Planck Institute for Immunobiology and Epigenetics, 79108, Freiburg, Germany.,Faculty of Biology, University of Freiburg, 79104, Freiburg, Germany
| | - Thomas Manke
- Max Planck Institute for Immunobiology and Epigenetics, 79108, Freiburg, Germany
| | - Asifa Akhtar
- Max Planck Institute for Immunobiology and Epigenetics, 79108, Freiburg, Germany.
| |
Collapse
|
26
|
Wu AC, Van Werven FJ. Transcribe this way: Rap1 confers promoter directionality by repressing divergent transcription. Transcription 2019; 10:164-170. [PMID: 31057041 PMCID: PMC6602560 DOI: 10.1080/21541264.2019.1608716] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Revised: 04/11/2019] [Accepted: 04/13/2019] [Indexed: 02/07/2023] Open
Abstract
In eukaryotes, divergent transcription is a major source of noncoding RNAs. Recent studies have uncovered that in yeast, the transcription factor Rap1 restricts transcription in the divergent direction and thereby controls promoter directionality. Here, we summarize these findings, propose regulatory principles, and discuss the implications for eukaryotic gene regulation.
Collapse
Affiliation(s)
- Andrew C.K. Wu
- Cell Fate and Gene Regulation Laboratory, The Francis Crick Institute, London, UK
| | | |
Collapse
|
27
|
Chu T, Wang Z, Chou SP, Danko CG. Discovering Transcriptional Regulatory Elements From Run-On and Sequencing Data Using the Web-Based dREG Gateway. CURRENT PROTOCOLS IN BIOINFORMATICS 2019; 66:e70. [PMID: 30589513 PMCID: PMC6584046 DOI: 10.1002/cpbi.70] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Transcription is a chromatin mark that can be used effectively to identify the location of active enhancers and promoters, collectively known as transcriptional regulatory elements (TREs). We recently introduced dREG, a tool for the identification of TREs using run-on and sequencing (RO-seq) assays, including global run-on and sequencing (GRO-seq), precision run-on and sequencing (PRO-seq), and chromatin run-on and sequencing (ChRO-seq). In this protocol, we present step-by-step instructions for running dREG on an arbitrary run-on and sequencing dataset. Users provide dREG with bigWig files (in which each read is represented by a single base) representing the location of RNA polymerase in a cell or tissue sample of interest, and dREG returns a list of genomic regions that are predicted to be active TREs. Finally, we demonstrate the use of dREG regions in discovering transcription factors controlling response to a stimulus and predicting their target genes. Together, this protocol provides detailed instructions for running dREG on arbitrary run-on and sequencing data. © 2018 by John Wiley & Sons, Inc.
Collapse
Affiliation(s)
- Tinyi Chu
- Baker Institute for Animal Health,Cornell University. Hungerford Hill Rd,Ithaca, NY 14853
| | - Zhong Wang
- Baker Institute for Animal Health,Cornell University. Hungerford Hill Rd,Ithaca, NY 14853
| | - Shao-Pei Chou
- Baker Institute for Animal Health,Cornell University. Hungerford Hill Rd,Ithaca, NY 14853
| | - Charles G. Danko
- Baker Institute for Animal Health,Cornell University. Hungerford Hill Rd,Ithaca, NY 14853
| |
Collapse
|
28
|
Abstract
In this review, Core et al. discuss the recent advances in our understanding of the early steps in Pol II transcription, highlighting the events and factors involved in the establishment and release of paused Pol II. They also discuss a number of unanswered questions about the regulation and function of Pol II pausing. Precise spatio–temporal control of gene activity is essential for organismal development, growth, and survival in a changing environment. Decisive steps in gene regulation involve the pausing of RNA polymerase II (Pol II) in early elongation, and the controlled release of paused polymerase into productive RNA synthesis. Here we describe the factors that enable pausing and the events that trigger Pol II release into the gene. We also discuss open questions in the field concerning the stability of paused Pol II, nucleosomes as obstacles to elongation, and potential roles of pausing in defining the precision and dynamics of gene expression.
Collapse
Affiliation(s)
- Leighton Core
- Department of Molecular and Cell Biology, Institute of Systems Genomics, University of Connecticut, Storrs, Connecticut 06269, USA
| | - Karen Adelman
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, Massachusetts 02115, USA
| |
Collapse
|
29
|
Jänes J, Dong Y, Schoof M, Serizay J, Appert A, Cerrato C, Woodbury C, Chen R, Gemma C, Huang N, Kissiov D, Stempor P, Steward A, Zeiser E, Sauer S, Ahringer J. Chromatin accessibility dynamics across C. elegans development and ageing. eLife 2018; 7:37344. [PMID: 30362940 PMCID: PMC6231769 DOI: 10.7554/elife.37344] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2018] [Accepted: 10/25/2018] [Indexed: 12/21/2022] Open
Abstract
An essential step for understanding the transcriptional circuits that control development and physiology is the global identification and characterization of regulatory elements. Here, we present the first map of regulatory elements across the development and ageing of an animal, identifying 42,245 elements accessible in at least one Caenorhabditis elegans stage. Based on nuclear transcription profiles, we define 15,714 protein-coding promoters and 19,231 putative enhancers, and find that both types of element can drive orientation-independent transcription. Additionally, more than 1000 promoters produce transcripts antisense to protein coding genes, suggesting involvement in a widespread regulatory mechanism. We find that the accessibility of most elements changes during development and/or ageing and that patterns of accessibility change are linked to specific developmental or physiological processes. The map and characterization of regulatory elements across C. elegans life provides a platform for understanding how transcription controls development and ageing.
Collapse
Affiliation(s)
- Jürgen Jänes
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Yan Dong
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Michael Schoof
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Jacques Serizay
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Alex Appert
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Chiara Cerrato
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Carson Woodbury
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Ron Chen
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Carolina Gemma
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Ni Huang
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Djem Kissiov
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Przemyslaw Stempor
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Annette Steward
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Eva Zeiser
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| | - Sascha Sauer
- Max Delbrück Center for Molecular Medicine, Berlin, Germany.,Max Planck Institute for Molecular Genetics, Otto-Warburg Laboratories, Berlin, Germany
| | - Julie Ahringer
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom.,The Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
30
|
Ibrahim MM, Karabacak A, Glahs A, Kolundzic E, Hirsekorn A, Carda A, Tursun B, Zinzen RP, Lacadie SA, Ohler U. Determinants of promoter and enhancer transcription directionality in metazoans. Nat Commun 2018; 9:4472. [PMID: 30367057 PMCID: PMC6203779 DOI: 10.1038/s41467-018-06962-z] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2017] [Accepted: 09/11/2018] [Indexed: 12/16/2022] Open
Abstract
Divergent transcription from promoters and enhancers is pervasive in many species, but it remains unclear if it is a general feature of all eukaryotic cis regulatory elements. To address this, here we define cis regulatory elements in C. elegans, D. melanogaster and H. sapiens and investigate the determinants of their transcription directionality. In all three species, we find that divergent transcription is initiated from two separate core promoter sequences and promoter regions display competition between histone modifications on the + 1 and −1 nucleosomes. In contrast, promoter directionality, sequence composition surrounding promoters, and positional enrichment of chromatin states, are different across species. Integrative models of H3K4me3 levels and core promoter sequence are highly predictive of promoter and enhancer directionality and support two directional classes, skewed and balanced. The relative importance of features to these models are clearly distinct for promoters and enhancers. Differences in regulatory architecture within and between metazoans are therefore abundant, arguing against a unified eukaryotic model. Divergent transcription from promoters and enhancers occurs in many species, but it is unclear if it is a general feature of all eukaryotic cis regulatory elements. Here the authors define cis regulatory elements in worms, flies, and human; and identify several differences in regulatory architecture among metazoans.
Collapse
Affiliation(s)
- Mahmoud M Ibrahim
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany.,Department of Biology, Humboldt Universitaet zu Berlin, 10115, Berlin, Germany.,Department of Nephrology and Immunology, Faculty of Medicine, RWTH Aachen University, Pauwelstraat 30, 52074, Aachen, Germany
| | - Aslihan Karabacak
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany.,Department of Biology, Humboldt Universitaet zu Berlin, 10115, Berlin, Germany
| | - Alexander Glahs
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany.,Department of Biology, Humboldt Universitaet zu Berlin, 10115, Berlin, Germany
| | - Ena Kolundzic
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany.,Department of Biology, Humboldt Universitaet zu Berlin, 10115, Berlin, Germany
| | - Antje Hirsekorn
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany
| | - Alexa Carda
- Department of Biostatistics & Bioinformatics, Duke University Medical Center, Durham, 27710, NC, USA
| | - Baris Tursun
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany
| | - Robert P Zinzen
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany
| | - Scott A Lacadie
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany. .,Berlin Institute of Health (BIH), Berlin, 10178, Germany.
| | - Uwe Ohler
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, 13125, Berlin, Germany. .,Department of Biology, Humboldt Universitaet zu Berlin, 10115, Berlin, Germany. .,Department of Biostatistics & Bioinformatics, Duke University Medical Center, Durham, 27710, NC, USA. .,Berlin Institute of Health (BIH), Berlin, 10178, Germany.
| |
Collapse
|
31
|
Khan A, Mathelier A, Zhang X. Super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers. Epigenetics 2018; 13:910-922. [PMID: 30169995 DOI: 10.1080/15592294.2018.1514231] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Super-enhancers and stretch enhancers represent classes of transcriptional enhancers that have been shown to control the expression of cell identity genes and carry disease- and trait-associated variants. Specifically, super-enhancers are clusters of enhancers defined based on the binding occupancy of master transcription factors, chromatin regulators, or chromatin marks, while stretch enhancers are large chromatin-defined regulatory regions of at least 3,000 base pairs. Several studies have characterized these regulatory regions in numerous cell types and tissues to decipher their functional importance. However, the differences and similarities between these regulatory regions have not been fully assessed. We integrated genomic, epigenomic, and transcriptomic data from ten human cell types to perform a comparative analysis of super and stretch enhancers with respect to their chromatin profiles, cell type-specificity, and ability to control gene expression. We found that stretch enhancers are more abundant, more distal to transcription start sites, cover twice as much the genome, and are significantly less conserved than super-enhancers. In contrast, super-enhancers are significantly more enriched for active chromatin marks and cohesin complex, and more transcriptionally active than stretch enhancers. Importantly, a vast majority of super-enhancers (85%) overlap with only a small subset of stretch enhancers (13%), which are enriched for cell type-specific biological functions, and control cell identity genes. These results suggest that super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers, and importantly, most of the stretch enhancers that are distinct from super-enhancers do not show an association with cell identity genes, are less active, and more likely to be poised enhancers.
Collapse
Affiliation(s)
- Aziz Khan
- a Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership , University of Oslo , Oslo , Norway.,b Key Lab of Bioinformatics/Bioinformatics Division, BNRIST (Beijing National Research Center for Information Science and Technology), Department of Automation , Tsinghua University , Beijing , China
| | - Anthony Mathelier
- a Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership , University of Oslo , Oslo , Norway.,c Department of Cancer Genetics, Institute for Cancer Research , Oslo University Hospital Radiumhospitalet , Oslo , Norway
| | - Xuegong Zhang
- b Key Lab of Bioinformatics/Bioinformatics Division, BNRIST (Beijing National Research Center for Information Science and Technology), Department of Automation , Tsinghua University , Beijing , China.,d School of Life Sciences , Tsinghua University , Beijing , China
| |
Collapse
|
32
|
Dao LTM, Spicuglia S. Transcriptional regulation by promoters with enhancer function. Transcription 2018; 9:307-314. [PMID: 29889606 PMCID: PMC6150634 DOI: 10.1080/21541264.2018.1486150] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2018] [Accepted: 06/01/2018] [Indexed: 12/31/2022] Open
Abstract
Promoters with enhancer activity have been described recently. In this point of view, we will discuss current findings highlighting the commonality of this type of regulatory elements, their genetic and epigenetic characteristics, their potential biological roles in the regulation of gene expression and the underlining molecular mechanisms. ABBREVIATIONS TSS: transcription start site; IFN: interferon; STARR-seq: Self-Transcribing Active Regulatory Region sequencing; MPRA: Massively Parallel Reporter Assay; ChIP: chromatin immunoprecipitation; CRISPR: Clustered Regularly Interspaced Short Palindromic Repeats; lncRNA: long non-coding RNA.
Collapse
Affiliation(s)
- Lan T. M. Dao
- Vinmec Research Institute of Stem cell and Gene technology (VRISG), Hanoi, Vietnam
| | - Salvatore Spicuglia
- Aix-Marseille University, INSERM, TAGC, UMR, Marseille, France
- Equipe Labéllisée Ligue Contre le Cancer
| |
Collapse
|