1
|
McQuarrie DWJ, Alizada A, Nicholson BC, Soller M. Rapid evolution of promoters from germline-specifically expressed genes including transposon silencing factors. BMC Genomics 2024; 25:678. [PMID: 38977960 PMCID: PMC11229233 DOI: 10.1186/s12864-024-10584-9] [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: 04/10/2024] [Accepted: 07/01/2024] [Indexed: 07/10/2024] Open
Abstract
BACKGROUND The piRNA pathway in animal gonads functions as an 'RNA-based immune system', serving to silence transposable elements and prevent inheritance of novel invaders. In Drosophila, this pathway relies on three gonad-specific Argonaute proteins (Argonaute-3, Aubergine and Piwi) that associate with 23-28 nucleotide piRNAs, directing the silencing of transposon-derived transcripts. Transposons constitute a primary driver of genome evolution, yet the evolution of piRNA pathway factors has not received in-depth exploration. Specifically, channel nuclear pore proteins, which impact piRNA processing, exhibit regions of rapid evolution in their promoters. Consequently, the question arises whether such a mode of evolution is a general feature of transposon silencing pathways. RESULTS By employing genomic analysis of coding and promoter regions within genes that function in transposon silencing in Drosophila, we demonstrate that the promoters of germ cell-specific piRNA factors are undergoing rapid evolution. Our findings indicate that rapid promoter evolution is a common trait among piRNA factors engaged in germline silencing across insect species, potentially contributing to gene expression divergence in closely related taxa. Furthermore, we observe that the promoters of genes exclusively expressed in germ cells generally exhibit rapid evolution, with some divergence in gene expression. CONCLUSION Our results suggest that increased germline promoter evolution, in partnership with other factors, could contribute to transposon silencing and evolution of species through differential expression of genes driven by invading transposons.
Collapse
Affiliation(s)
- David W J McQuarrie
- School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
- Birmingham Centre for Genome Biology, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
| | - Azad Alizada
- Cancer Research UK Cambridge Institute, University of Cambridge, Li Ka Shing Centre, Cambridge, CB2 0RE, UK
| | - Benjamin Czech Nicholson
- Cancer Research UK Cambridge Institute, University of Cambridge, Li Ka Shing Centre, Cambridge, CB2 0RE, UK
| | - Matthias Soller
- School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK.
- Birmingham Centre for Genome Biology, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK.
| |
Collapse
|
2
|
Farmiloe G, van Bree EJ, Robben SF, Janssen LJM, Mol L, Jacobs FMJ. Structural Evolution of Gene Promoters Driven by Primate-Specific KRAB Zinc Finger Proteins. Genome Biol Evol 2023; 15:evad184. [PMID: 37847041 PMCID: PMC10653712 DOI: 10.1093/gbe/evad184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 09/26/2023] [Accepted: 10/09/2023] [Indexed: 10/18/2023] Open
Abstract
Krüppel-associated box (KRAB) zinc finger proteins (KZNFs) recognize and repress transposable elements (TEs); TEs are DNA elements that are capable of replicating themselves throughout our genomes with potentially harmful consequences. However, genes from this family of transcription factors have a much wider potential for genomic regulation. KZNFs have become integrated into gene-regulatory networks through the control of TEs that function as enhancers and gene promoters; some KZNFs also bind directly to gene promoters, suggesting an additional, more direct layer of KZNF co-option into gene-regulatory networks. Binding site analysis of ZNF519, ZNF441, and ZNF468 suggests the structural evolution of KZNFs to recognize TEs can result in coincidental binding to gene promoters independent of TE sequences. We show a higher rate of sequence turnover in gene promoter KZNF binding sites than neighboring regions, implying a selective pressure is being applied by the binding of a KZNF. Through CRISPR/Cas9 mediated genetic deletion of ZNF519, ZNF441, and ZNF468, we provide further evidence for genome-wide co-option of the KZNF-mediated gene-regulatory functions; KZNF knockout leads to changes in expression of KZNF-bound genes in neuronal lineages. Finally, we show that the opposite can be established upon KZNF overexpression, further strengthening the support for the role of KZNFs as bona-fide gene regulators. With no eminent role for ZNF519 in controlling its TE target, our study may provide a snapshot into the early stages of the completed co-option of a KZNF, showing the lasting, multilayered impact that retrovirus invasions and host response mechanisms can have upon the evolution of our genomes.
Collapse
Affiliation(s)
- Grace Farmiloe
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
- Complex Trait Genetics, Amsterdam Neuroscience, Amsterdam, The Netherlands
| | - Elisabeth J van Bree
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
- Complex Trait Genetics, Amsterdam Neuroscience, Amsterdam, The Netherlands
| | - Stijn F Robben
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
| | - Lara J M Janssen
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
| | - Lisa Mol
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
| | - Frank M J Jacobs
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
- Complex Trait Genetics, Amsterdam Neuroscience, Amsterdam, The Netherlands
| |
Collapse
|
3
|
McQuarrie DWJ, Read AM, Stephens FHS, Civetta A, Soller M. Indel driven rapid evolution of core nuclear pore protein gene promoters. Sci Rep 2023; 13:8035. [PMID: 37198214 DOI: 10.1038/s41598-023-34985-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 05/11/2023] [Indexed: 05/19/2023] Open
Abstract
Nuclear pore proteins (Nups) prominently are among the few genes linked to speciation from hybrid incompatibility in Drosophila. These studies have focused on coding sequence evolution of Nup96 and Nup160 and shown evidence of positive selection driving nucleoporin evolution. Intriguingly, channel Nup54 functionality is required for neuronal wiring underlying the female post-mating response induced by male-derived sex-peptide. A region of rapid evolution in the core promoter of Nup54 suggests a critical role for general transcriptional regulatory elements at the onset of speciation, but whether this is a general feature of Nup genes has not been determined. Consistent with findings for Nup54, additional channel Nup58 and Nup62 promoters also rapidly accumulate insertions/deletions (indels). Comprehensive examination of Nup upstream regions reveals that core Nup complex gene promoters accumulate indels rapidly. Since changes in promoters can drive changes in expression, these results indicate an evolutionary mechanism driven by indel accumulation in core Nup promoters. Compensation of such gene expression changes could lead to altered neuronal wiring, rapid fixation of traits caused by promoter changes and subsequently the rise of new species. Hence, the nuclear pore complex may act as a nexus for species-specific changes via nucleo-cytoplasmic transport regulated gene expression.
Collapse
Affiliation(s)
- David W J McQuarrie
- School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
- Birmingham Centre for Genome Biology, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
| | - Adam M Read
- School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
| | - Frannie H S Stephens
- School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
| | - Alberto Civetta
- Department of Biology, University of Winnipeg, Winnipeg, MB, R3B 2E9, Canada.
| | - Matthias Soller
- School of Biosciences, College of Life and Environmental Sciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK.
- Birmingham Centre for Genome Biology, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK.
| |
Collapse
|
4
|
Ji X, Cai J, Liang L, Shi T, Liu J. Gene expression variability across cells and species shapes the relationship between renal resident macrophages and infiltrated macrophages. BMC Bioinformatics 2023; 24:72. [PMID: 36858955 PMCID: PMC9976410 DOI: 10.1186/s12859-023-05198-z] [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: 08/13/2022] [Accepted: 02/21/2023] [Indexed: 03/03/2023] Open
Abstract
BACKGROUND Two main subclasses of macrophages are found in almost all solid tissues: embryo-derived resident tissue macrophages and bone marrow-derived infiltrated macrophages. These macrophage subtypes show transcriptional and functional divergence, and the programs that have shaped the evolution of renal macrophages and related signaling pathways remain poorly understood. To clarify these processes, we performed data analysis based on single-cell transcriptional profiling of renal tissue-resident and infiltrated macrophages in human, mouse and rat. RESULTS In this study, we (i) characterized the transcriptional divergence among species and (ii) illustrated variability in expression among cells of each subtype and (iii) compared the gene regulation network and (iv) ligand-receptor pairs in human and mouse. Using single-cell transcriptomics, we mapped the promoter architecture during homeostasis. CONCLUSIONS Transcriptionally divergent genes, such as the differentially TF-encoding genes expressed in resident and infiltrated macrophages across the three species, vary among cells and include distinct promoter structures. The gene regulatory network in infiltrated macrophages shows comparatively better species-wide consistency than resident macrophages. The conserved transcriptional gene regulatory network in infiltrated macrophages among species is uniquely enriched in pathways related to kinases, and TFs associated with largely conserved regulons among species are uniquely enriched in kinase-related pathways.
Collapse
Affiliation(s)
- Xiangjun Ji
- grid.284723.80000 0000 8877 7471Guangdong Provincial Key Laboratory of Proteomics, Department of Pathophysiology, School of Basic Medical Sciences, Southern Medical University, Guangzhou, 510515 China
| | - Junwei Cai
- grid.284723.80000 0000 8877 7471Guangdong Provincial Key Laboratory of Proteomics, Department of Pathophysiology, School of Basic Medical Sciences, Southern Medical University, Guangzhou, 510515 China
| | - Lixin Liang
- grid.284723.80000 0000 8877 7471Guangdong Provincial Key Laboratory of Proteomics, Department of Pathophysiology, School of Basic Medical Sciences, Southern Medical University, Guangzhou, 510515 China
| | - Tieliu Shi
- Center for Bioinformatics and Computational Biology, The Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China. .,Beijing Advanced Innovation Center, for Big Data-Based Precision Medicine, Beihang University and Capital Medical University, Beijing, 100083, China.
| | - Jinghua Liu
- Guangdong Provincial Key Laboratory of Proteomics, Department of Pathophysiology, School of Basic Medical Sciences, Southern Medical University, Guangzhou, 510515, China.
| |
Collapse
|
5
|
Fang Y, Deng S, Li C. A generalizable deep learning framework for inferring fine-scale germline mutation rate maps. NAT MACH INTELL 2022. [DOI: 10.1038/s42256-022-00574-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
6
|
Young RS, Talmane L, Marion de Procé S, Taylor MS. The contribution of evolutionarily volatile promoters to molecular phenotypes and human trait variation. Genome Biol 2022; 23:89. [PMID: 35379293 PMCID: PMC8978360 DOI: 10.1186/s13059-022-02634-w] [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] [Received: 05/07/2021] [Accepted: 02/16/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Promoters are sites of transcription initiation that harbour a high concentration of phenotype-associated genetic variation. The evolutionary gain and loss of promoters between species (collectively, termed turnover) is pervasive across mammalian genomes and may play a prominent role in driving human phenotypic diversity. RESULTS We classified human promoters by their evolutionary history during the divergence of mouse and human lineages from a common ancestor. This defined conserved, human-inserted and mouse-deleted promoters, and a class of functional-turnover promoters that align between species but are only active in humans. We show that promoters of all evolutionary categories are hotspots for substitution and often, insertion mutations. Loci with a history of insertion and deletion continue that mode of evolution within contemporary humans. The presence of an evolutionary volatile promoter within a gene is associated with increased expression variance between individuals, but only in the case of human-inserted and mouse-deleted promoters does that correspond to an enrichment of promoter-proximal genetic effects. Despite the enrichment of these molecular quantitative trait loci (QTL) at evolutionarily volatile promoters, this does not translate into a corresponding enrichment of phenotypic traits mapping to these loci. CONCLUSIONS Promoter turnover is pervasive in the human genome, and these promoters are rich in molecularly quantifiable but phenotypically inconsequential variation in gene expression. However, since evolutionarily volatile promoters show evidence of selection, coupled with high mutation rates and enrichment of QTLs, this implicates them as a source of evolutionary innovation and phenotypic variation, albeit with a high background of selectively neutral expression variation.
Collapse
Affiliation(s)
- Robert S Young
- Usher Institute, University of Edinburgh, Teviot Place, Edinburgh, EH8 9AG, UK. .,Zhejiang University - University of Edinburgh Institute, Zhejiang University, 718 East Haizhou Road, 314400, Haining, China. .,MRC Human Genetics Unit, Institute for Genetics and Cancer, University of Edinburgh, Crewe Road, Edinburgh, EH4 2XU, UK.
| | - Lana Talmane
- MRC Human Genetics Unit, Institute for Genetics and Cancer, University of Edinburgh, Crewe Road, Edinburgh, EH4 2XU, UK
| | - Sophie Marion de Procé
- Usher Institute, University of Edinburgh, Teviot Place, Edinburgh, EH8 9AG, UK.,MRC Human Genetics Unit, Institute for Genetics and Cancer, University of Edinburgh, Crewe Road, Edinburgh, EH4 2XU, UK
| | - Martin S Taylor
- MRC Human Genetics Unit, Institute for Genetics and Cancer, University of Edinburgh, Crewe Road, Edinburgh, EH4 2XU, UK.
| |
Collapse
|
7
|
Huang D, Wang X, Liu Y, Huang Z, Hu X, Hu W, Li Q, Chan H, Zou Y, Ho IHT, Wang Y, Cheng ASL, Kang W, To KF, Wang MHT, Wong SH, Yu J, Gin T, Zhang Q, Li Z, Shen J, Zhang L, Chan MTV, Liu X, Wu WKK. Multi-omic analysis suggests tumor suppressor genes evolved specific promoter features to optimize cancer resistance. Brief Bioinform 2021; 22:6200210. [PMID: 33783485 DOI: 10.1093/bib/bbab040] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Revised: 12/23/2020] [Accepted: 01/28/2021] [Indexed: 12/31/2022] Open
Abstract
Tumor suppressor genes (TSGs) exhibit distinct evolutionary features. We speculated that TSG promoters could have evolved specific features that facilitate their tumor-suppressing functions. We found that the promoter CpG dinucleotide frequencies of TSGs are significantly higher than that of non-cancer genes across vertebrate genomes, and positively correlated with gene expression across tissue types. The promoter CpG dinucleotide frequencies of all genes gradually increase with gene age, for which young TSGs have been subject to a stronger evolutionary pressure. Transcription-related features, namely chromatin accessibility, methylation and ZNF263-, SP1-, E2F4- and SP2-binding elements, are associated with gene expression. Moreover, higher promoter CpG dinucleotide frequencies and chromatin accessibility are positively associated with the ability of TSGs to resist downregulation during tumorigenesis. These results were successfully validated with independent datasets. In conclusion, TSGs evolved specific promoter features that optimized cancer resistance through achieving high expression in normal tissues and resistance to downregulation during tumorigenesis.
Collapse
Affiliation(s)
- Dan Huang
- Chinese University of Hong Kong and the CUHK-Shenzhen Research Institute, China
| | - Xiansong Wang
- Chinese University of Hong Kong and the CUHK-Shenzhen Research Institute, China
| | - Yingzhi Liu
- Chinese University of Hong Kong and the CUHK-Shenzhen Research Institute, China
| | - Ziheng Huang
- Chinese University of Hong Kong and the CUHK-Shenzhen Research Institute, China
| | - Xiaoxu Hu
- Chinese University of Hong Kong and the CUHK-Shenzhen Research Institute, China
| | - Wei Hu
- Chinese University of Hong Kong, China
| | - Qing Li
- Chinese University of Hong Kong, China
| | - Hung Chan
- Chinese University of Hong Kong, China
| | - Yidan Zou
- Chinese University of Hong Kong, China
| | | | - Yan Wang
- Chinese University of Hong Kong, China
| | | | - Wei Kang
- Chinese University of Hong Kong, China
| | - Ka F To
- Chinese University of Hong Kong, China
| | - Maggie H T Wang
- Chinese University of Hong Kong and the CUHK-Shenzhen Research Institute, China
| | | | - Jun Yu
- Chinese University of Hong Kong, China
| | - Tony Gin
- Chinese University of Hong Kong, China
| | | | - Zheng Li
- Peking Union Medical College Hospital, China
| | | | - Lin Zhang
- Chinese University of Hong Kong, China
| | | | | | - William K K Wu
- Chinese University of Hong Kong and a researcher at the CUHK-Shenzhen Research Institute, China
| |
Collapse
|
8
|
Halachev M, Meynert A, Taylor MS, Vitart V, Kerr SM, Klaric L, Aitman TJ, Haley CS, Prendergast JG, Pugh C, Hume DA, Harris SE, Liewald DC, Deary IJ, Semple CA, Wilson JF. Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions. PLoS Genet 2019; 15:e1008480. [PMID: 31765389 PMCID: PMC6901239 DOI: 10.1371/journal.pgen.1008480] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2019] [Revised: 12/09/2019] [Accepted: 10/15/2019] [Indexed: 01/03/2023] Open
Abstract
Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.
Collapse
Affiliation(s)
- Mihail Halachev
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | - Alison Meynert
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | - Martin S. Taylor
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | - Veronique Vitart
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | - Shona M. Kerr
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | - Lucija Klaric
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | | | - Timothy J. Aitman
- Centre for Genomic and Experimental Medicine, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | - Chris S. Haley
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, United Kingdom
| | - James G. Prendergast
- The Roslin Institute, University of Edinburgh, Easter Bush, Midlothian, United Kingdom
| | - Carys Pugh
- Centre for Clinical Brain Sciences, Division of Psychiatry, University of Edinburgh, Royal Edinburgh Hospital, Edinburgh, United Kingdom
| | - David A. Hume
- Mater Research Institute, University of Queensland, Woolloongabba, Australia
| | - Sarah E. Harris
- Centre for Cognitive Ageing and Cognitive Epidemiology, Department of Psychology, School of Philosophy, Psychology and Language Sciences, University of Edinburgh, George Square, Edinburgh, United Kingdom
| | - David C. Liewald
- Centre for Cognitive Ageing and Cognitive Epidemiology, Department of Psychology, School of Philosophy, Psychology and Language Sciences, University of Edinburgh, George Square, Edinburgh, United Kingdom
| | - Ian J. Deary
- Centre for Cognitive Ageing and Cognitive Epidemiology, Department of Psychology, School of Philosophy, Psychology and Language Sciences, University of Edinburgh, George Square, Edinburgh, United Kingdom
| | - Colin A. Semple
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
| | - James F. Wilson
- MRC Human Genetics Unit, MRC IGMM, University of Edinburgh, Crewe Road, Edinburgh, United Kingdom
- Centre for Global Health Research, Usher Institute of Population Health Sciences and Informatics, University of Edinburgh, Teviot Place, Edinburgh, United Kingdom
| |
Collapse
|
9
|
Nute M, Saleh E, Warnow T. Evaluating Statistical Multiple Sequence Alignment in Comparison to Other Alignment Methods on Protein Data Sets. Syst Biol 2019; 68:396-411. [PMID: 30329135 PMCID: PMC6472439 DOI: 10.1093/sysbio/syy068] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 09/27/2018] [Accepted: 10/11/2018] [Indexed: 01/15/2023] Open
Abstract
The estimation of multiple sequence alignments of protein sequences is a basic step in many bioinformatics pipelines, including protein structure prediction, protein family identification, and phylogeny estimation. Statistical coestimation of alignments and trees under stochastic models of sequence evolution has long been considered the most rigorous technique for estimating alignments and trees, but little is known about the accuracy of such methods on biological benchmarks. We report the results of an extensive study evaluating the most popular protein alignment methods as well as the statistical coestimation method BAli-Phy on 1192 protein data sets from established benchmarks as well as on 120 simulated data sets. Our study (which used more than 230 CPU years for the BAli-Phy analyses alone) shows that BAli-Phy has better precision and recall (with respect to the true alignments) than the other alignment methods on the simulated data sets but has consistently lower recall on the biological benchmarks (with respect to the reference alignments) than many of the other methods. In other words, we find that BAli-Phy systematically underaligns when operating on biological sequence data but shows no sign of this on simulated data. There are several potential causes for this change in performance, including model misspecification, errors in the reference alignments, and conflicts between structural alignment and evolutionary alignments, and future research is needed to determine the most likely explanation. We conclude with a discussion of the potential ramifications for each of these possibilities. [BAli-Phy; homology; multiple sequence alignment; protein sequences; structural alignment.]
Collapse
Affiliation(s)
- Michael Nute
- Department of Statistics, University of Illinois at Urbana-Champaign, 725 S Wright St #101, Champaign, IL 61820, USA
| | - Ehsan Saleh
- Department of Computer Science, University of Illinois at Urbana-Champaign, 201 N. Goodwin Ave, Urbana, IL 61801, USA
| | - Tandy Warnow
- Department of Computer Science, University of Illinois at Urbana-Champaign, 201 N. Goodwin Ave, Urbana, IL 61801, USA.,Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, 1205 W. Clark St., Urbana, IL 61801, USA.,National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
10
|
Neininger K, Marschall T, Helms V. SNP and indel frequencies at transcription start sites and at canonical and alternative translation initiation sites in the human genome. PLoS One 2019; 14:e0214816. [PMID: 30978217 PMCID: PMC6461226 DOI: 10.1371/journal.pone.0214816] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2018] [Accepted: 03/20/2019] [Indexed: 11/30/2022] Open
Abstract
Single-nucleotide polymorphisms (SNPs) are the most common form of genetic variation in humans and drive phenotypic variation. Due to evolutionary conservation, SNPs and indels (insertion and deletions) are depleted in functionally important sequence elements. Recently, population-scale sequencing efforts such as the 1000 Genomes Project and the Genome of the Netherlands Project have catalogued large numbers of sequence variants. Here, we present a systematic analysis of the polymorphisms reported by these two projects in different coding and non-coding genomic elements of the human genome (intergenic regions, CpG islands, promoters, 5’ UTRs, coding exons, 3’ UTRs, introns, and intragenic regions). Furthermore, we were especially interested in the distribution of SNPs and indels in direct vicinity to the transcription start site (TSS) and translation start site (CSS). Thereby, we discovered an enrichment of dinucleotides CpG and CpA and an accumulation of SNPs at base position −1 relative to the TSS that involved primarily CpG and CpA dinucleotides. Genes having a CpG dinucleotide at TSS position -1 were enriched in the functional GO terms “Phosphoprotein”, “Alternative splicing”, and “Protein binding”. Focusing on the CSS, we compared SNP patterns in the flanking regions of canonical and alternative AUG and near-cognate start sites where we considered alternative starts previously identified by experimental ribosome profiling. We observed similar conservation patterns of canonical and alternative translation start sites, which underlines the importance of alternative translation mechanisms for cellular function.
Collapse
Affiliation(s)
- Kerstin Neininger
- Center for Bioinformatics, Saarland University, 66123 Saarbrücken, Germany
- Graduate School of Computer Science, Saarland University, 66123 Saarbrücken, Germany
| | - Tobias Marschall
- Center for Bioinformatics, Saarland University, 66123 Saarbrücken, Germany
- Max Planck Institute for Informatics, 66123 Saarbrücken, Germany
| | - Volkhard Helms
- Center for Bioinformatics, Saarland University, 66123 Saarbrücken, Germany
- * E-mail:
| |
Collapse
|
11
|
Promoter analysis and prediction in the human genome using sequence-based deep learning models. Bioinformatics 2019; 35:2730-2737. [DOI: 10.1093/bioinformatics/bty1068] [Citation(s) in RCA: 60] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2018] [Revised: 12/03/2018] [Accepted: 12/27/2018] [Indexed: 12/14/2022] Open
Abstract
Abstract
Motivation
Computational identification of promoters is notoriously difficult as human genes often have unique promoter sequences that provide regulation of transcription and interaction with transcription initiation complex. While there are many attempts to develop computational promoter identification methods, we have no reliable tool to analyze long genomic sequences.
Results
In this work, we further develop our deep learning approach that was relatively successful to discriminate short promoter and non-promoter sequences. Instead of focusing on the classification accuracy, in this work we predict the exact positions of the transcription start site inside the genomic sequences testing every possible location. We studied human promoters to find effective regions for discrimination and built corresponding deep learning models. These models use adaptively constructed negative set, which iteratively improves the model’s discriminative ability. Our method significantly outperforms the previously developed promoter prediction programs by considerably reducing the number of false-positive predictions. We have achieved error-per-1000-bp rate of 0.02 and have 0.31 errors per correct prediction, which is significantly better than the results of other human promoter predictors.
Availability and implementation
The developed method is available as a web server at http://www.cbrc.kaust.edu.sa/PromID/.
Collapse
|
12
|
Li C, Lenhard B, Luscombe NM. Integrated analysis sheds light on evolutionary trajectories of young transcription start sites in the human genome. Genome Res 2018; 28:676-688. [PMID: 29618487 PMCID: PMC5932608 DOI: 10.1101/gr.231449.117] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2017] [Accepted: 03/21/2018] [Indexed: 01/06/2023]
Abstract
Understanding the molecular mechanisms and evolution of the gene regulatory system remains a major challenge in biology. Transcription start sites (TSSs) are especially interesting because they are central to initiating gene expression. Previous studies revealed widespread transcription initiation and fast turnover of TSSs in mammalian genomes. Yet, how new TSSs originate and how they evolve over time remain poorly understood. To address these questions, we analyzed ∼200,000 human TSSs by integrating evolutionary (inter- and intra-species) and functional genomic data, particularly focusing on evolutionarily young TSSs that emerged in the primate lineage. TSSs were grouped according to their evolutionary age using sequence alignment information as a proxy. Comparisons of young and old TSSs revealed that (1) new TSSs emerge through a combination of intrinsic factors, like the sequence properties of transposable elements and tandem repeats, and extrinsic factors such as their proximity to existing regulatory modules; (2) new TSSs undergo rapid evolution that reduces the inherent instability of repeat sequences associated with a high propensity of TSS emergence; and (3) once established, the transcriptional competence of surviving TSSs is gradually enhanced, with evolutionary changes subject to temporal (fewer regulatory changes in younger TSSs) and spatial constraints (fewer regulatory changes in more isolated TSSs). These findings advance our understanding of how regulatory innovations arise in the genome throughout evolution and highlight the genomic robustness and evolvability in these processes.
Collapse
Affiliation(s)
- Cai Li
- The Francis Crick Institute, London NW1 1AT, United Kingdom
| | - Boris Lenhard
- Computational Regulatory Genomics, MRC London Institute of Medical Sciences, London W12 0NN, United Kingdom.,Institute of Clinical Sciences, Faculty of Medicine, Imperial College London, London W12 0NN, United Kingdom.,Sars International Centre for Marine Molecular Biology, University of Bergen, N-5008 Bergen, Norway
| | - Nicholas M Luscombe
- The Francis Crick Institute, London NW1 1AT, United Kingdom.,UCL Genetics Institute, University College London, London WC1E 6BT, United Kingdom.,Okinawa Institute of Science & Technology Graduate University, Okinawa, 904-0495, Japan
| |
Collapse
|
13
|
Sequence and functional characterization of MIRNA164 promoters from Brassica shows copy number dependent regulatory diversification among homeologs. Funct Integr Genomics 2018. [DOI: 10.1007/s10142-018-0598-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
14
|
Zhai Y, Alexandre BC. A Poissonian Model of Indel Rate Variation for Phylogenetic Tree Inference. Syst Biol 2018; 66:698-714. [PMID: 28204784 DOI: 10.1093/sysbio/syx033] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2015] [Accepted: 01/27/2017] [Indexed: 01/22/2023] Open
Abstract
While indel rate variation has been observed and analyzed in detail, it is not taken into account by current indel-aware phylogenetic reconstruction methods. In this work, we introduce a continuous time stochastic process, the geometric Poisson indel process, that generalizes the Poisson indel process by allowing insertion and deletion rates to vary across sites. We design an efficient algorithm for computing the probability of a given multiple sequence alignment based on our new indel model. We describe a method to construct phylogeny estimates from a fixed alignment using neighbor joining. Using simulation studies, we show that ignoring indel rate variation may have a detrimental effect on the accuracy of the inferred phylogenies, and that our proposed method can sidestep this issue by inferring latent indel rate categories. We also show that our phylogenetic inference method may be more stable to taxa subsampling than methods that either ignore indels or indel rate variation. [evolutionary stochastic process; indel rate variation; Poisson indel process; TKF91.].
Collapse
Affiliation(s)
- Yongliang Zhai
- Department of Statistics, University of British Columbia, Vancouver, British Columbia, V6T 1Z4, Canada
| | - Bouchard-Côté Alexandre
- Department of Statistics, University of British Columbia, Vancouver, British Columbia, V6T 1Z4, Canada
| |
Collapse
|
15
|
Hardingham GE, Pruunsild P, Greenberg ME, Bading H. Lineage divergence of activity-driven transcription and evolution of cognitive ability. Nat Rev Neurosci 2017; 19:9-15. [PMID: 29167525 DOI: 10.1038/nrn.2017.138] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Excitation-transcription coupling shapes network formation during brain development and controls neuronal survival, synaptic function and cognitive skills in the adult. New studies have uncovered differences in the transcriptional responses to synaptic activity between humans and mice. These differences are caused both by the emergence of lineage-specific activity-regulated genes and by the acquisition of signal-responsive DNA elements in gene regulatory regions that determine whether a gene can be transcriptionally induced by synaptic activity or alter the extent of its inducibility. Such evolutionary divergence may have contributed to lineage-related advancements in cognitive abilities.
Collapse
Affiliation(s)
- Giles E Hardingham
- UK Dementia Research Institute at The University of Edinburgh, Edinburgh Medical School, 47 Little France Crescent, Edinburgh EH16 4TJ, UK
| | - Priit Pruunsild
- Department of Neurobiology, Interdisciplinary Center for Neurosciences (IZN), Heidelberg University, 69120 Heidelberg, Germany
| | - Michael E Greenberg
- Department of Neurobiology, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Hilmar Bading
- Department of Neurobiology, Interdisciplinary Center for Neurosciences (IZN), Heidelberg University, 69120 Heidelberg, Germany
| |
Collapse
|
16
|
Tokizawa M, Kusunoki K, Koyama H, Kurotani A, Sakurai T, Suzuki Y, Sakamoto T, Kurata T, Yamamoto YY. Identification of Arabidopsis genic and non-genic promoters by paired-end sequencing of TSS tags. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2017; 90:587-605. [PMID: 28214361 DOI: 10.1111/tpj.13511] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2016] [Revised: 02/02/2017] [Accepted: 02/06/2017] [Indexed: 06/06/2023]
Abstract
Information about transcription start sites (TSSs) provides baseline data for the analysis of promoter architecture. In this paper we used paired- and single-end deep sequencing to analyze Arabidopsis TSS tags from several libraries prepared from roots, shoots, flowers and etiolated seedlings. The clustering of approximately 33 million mapped TSS tags led to the identification of 324 461 promoters that covered 79.7% (21 672/27 206) of protein-coding genes in the Arabidopsis genome. In addition we identified intragenic, antisense and orphan promoters that were not associated with any gene models. Of these, intragenic promoters exhibited unique characteristics regarding dinucleotide sequences at TSSs and core promoter element composition, suggesting that these promoters use different mechanisms of transcriptional initiation. An analysis of base composition with regard to promoter position revealed a low GC content throughout the promoter region and several local strand biases that were evident for TATA-type promoters, but not for Coreless-type promoters. Most observed strand biases coincided with strand biases of single nucleotide polymorphism rate. Our analysis also revealed that transcription of a gene is supported by an average of 2.7 genic promoters, among which one specific promoter, designated as a top promoter, substantially determines the expression level of the gene.
Collapse
Affiliation(s)
- Mutsutomo Tokizawa
- United Graduate School of Agriculture, Gifu University, Yanagido 1-1, Gifu City, Gifu, 501-1193, Japan
| | - Kazutaka Kusunoki
- United Graduate School of Agriculture, Gifu University, Yanagido 1-1, Gifu City, Gifu, 501-1193, Japan
| | - Hiroyuki Koyama
- United Graduate School of Agriculture, Gifu University, Yanagido 1-1, Gifu City, Gifu, 501-1193, Japan
- Faculty of Applied Biological Sciences, Gifu University, Yanagido 1-1, Gifu City, Gifu, 501-1193, Japan
| | - Atsushi Kurotani
- RIKEN Center for Sustainable Resource Science, Suehiro-cho 1-7-22, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan
| | - Tetsuya Sakurai
- RIKEN Center for Sustainable Resource Science, Suehiro-cho 1-7-22, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan
| | - Yutaka Suzuki
- Institute of Medical Science, University of Tokyo, Shiroganedai 4-6-1, Minato-ku, Tokyo, 108-8639, Japan
| | - Tomoaki Sakamoto
- Plant Global Education Project, Graduate School of Biological Sciences, Nara Institute of Science and Technology, Takayam-cho 8916-5, Ikoma, Nara, 630-0192, Japan
| | - Tetsuya Kurata
- Plant Global Education Project, Graduate School of Biological Sciences, Nara Institute of Science and Technology, Takayam-cho 8916-5, Ikoma, Nara, 630-0192, Japan
| | - Yoshiharu Y Yamamoto
- United Graduate School of Agriculture, Gifu University, Yanagido 1-1, Gifu City, Gifu, 501-1193, Japan
- Faculty of Applied Biological Sciences, Gifu University, Yanagido 1-1, Gifu City, Gifu, 501-1193, Japan
- RIKEN Center for Sustainable Resource Science, Suehiro-cho 1-7-22, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan
- JST ALCA, Tokyo, Japan
| |
Collapse
|
17
|
Yella VR, Bansal M. DNA structural features of eukaryotic TATA-containing and TATA-less promoters. FEBS Open Bio 2017; 7:324-334. [PMID: 28286728 PMCID: PMC5337902 DOI: 10.1002/2211-5463.12166] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Accepted: 11/16/2016] [Indexed: 01/26/2023] Open
Abstract
Eukaryotic genes can be broadly classified as TATA‐containing and TATA‐less based on the presence of TATA box in their promoters. Experiments on both classes of genes have revealed a disparity in the regulation of gene expression and cellular functions between the two classes. In this study, we report characteristic differences in promoter sequences and associated structural properties of the two categories of genes in six different eukaryotes. We have analyzed three structural features, DNA duplex stability, bendability, and curvature along with the distribution of A‐tracts, G‐quadruplex motifs, and CpG islands. The structural feature analyses reveal that while the two classes of gene promoters are distinctly different from each other, the properties are also distinguishable across the six organisms.
Collapse
Affiliation(s)
- Venkata Rajesh Yella
- Molecular Biophysics Unit Indian Institute of Science Bangalore India; Present address: Department of Biotechnology K L University, Vaddeswaram Guntur 522502 India
| | - Manju Bansal
- Molecular Biophysics Unit Indian Institute of Science Bangalore India
| |
Collapse
|
18
|
Dong X, Wang X, Zhang F, Tian W. Genome-Wide Identification of Regulatory Sequences Undergoing Accelerated Evolution in the Human Genome. Mol Biol Evol 2016; 33:2565-75. [PMID: 27401230 PMCID: PMC5026254 DOI: 10.1093/molbev/msw128] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Accelerated evolution of regulatory sequence can alter the expression pattern of target genes, and cause phenotypic changes. In this study, we used DNase I hypersensitive sites (DHSs) to annotate putative regulatory sequences in the human genome, and conducted a genome-wide analysis of the effects of accelerated evolution on regulatory sequences. Working under the assumption that local ancient repeat elements of DHSs are under neutral evolution, we discovered that ∼0.44% of DHSs are under accelerated evolution (ace-DHSs). We found that ace-DHSs tend to be more active than background DHSs, and are strongly associated with epigenetic marks of active transcription. The target genes of ace-DHSs are significantly enriched in neuron-related functions, and their expression levels are positively selected in the human brain. Thus, these lines of evidences strongly suggest that accelerated evolution on regulatory sequences plays important role in the evolution of human-specific phenotypes.
Collapse
Affiliation(s)
- Xinran Dong
- State Key Laboratory of Genetic Engineering, Collaborative Innovation Center of Genetics and Development, Department of Biostatistics and Computational Biology, School of Life Sciences, Fudan University, Shanghai, P.R. China
| | - Xiao Wang
- State Key Laboratory of Genetic Engineering, Collaborative Innovation Center of Genetics and Development, Department of Biostatistics and Computational Biology, School of Life Sciences, Fudan University, Shanghai, P.R. China
| | - Feng Zhang
- State Key Laboratory of Genetic Engineering, Collaborative Innovation Center of Genetics and Development, Department of Biostatistics and Computational Biology, School of Life Sciences, Fudan University, Shanghai, P.R. China
| | - Weidong Tian
- State Key Laboratory of Genetic Engineering, Collaborative Innovation Center of Genetics and Development, Department of Biostatistics and Computational Biology, School of Life Sciences, Fudan University, Shanghai, P.R. China Children's Hospital of Fudan University, Shanghai, P.R. China
| |
Collapse
|
19
|
Yang H, He J, Wei W, Chu W, Yu S, Tian Y, Peng F, Liu H, Zhang Z, Chen J. The c.-360 T>C mutation affects PGAM2 transcription activity and is linked with the water holding capacity of the longissimus lumborum muscle in pigs. Meat Sci 2016; 122:139-144. [PMID: 27538264 DOI: 10.1016/j.meatsci.2016.07.023] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2016] [Revised: 07/26/2016] [Accepted: 07/28/2016] [Indexed: 11/18/2022]
Abstract
The phosphoglycerate mutase 2 (PGAM2) gene encodes a key enzyme in the glycolytic process. This study examined a functional mutation in the PGAM2 gene and evaluated its relationship with water holding capacity (WHC). RT-qPCR analysis showed the PGAM2 mRNA level was significantly higher in the low-WHC group than in the high-WHC group (P<0.05). The c.-360 T>C mutation was identified through sequencing and found to have opposite allele distributions in the two groups. The allele was further genotyped in 170 Duroc×Large White×Yorkshire crossbred pigs using allele-specific PCR. The CC genotype was associated with lower WHC and higher PGAM2 mRNA levels, whereas the TT genotype corresponded to a higher WHC and lower PGAM2 mRNA levels (P<0.05). A luciferase activity assay also showed that the CC-genotype promoter had higher activity than the TT-genotype promoter (P<0.05). In conclusion, we discovered the c.-360 T>C mutation in the PGAM2 gene, which is a promising marker for improving pork WHC.
Collapse
Affiliation(s)
- Haoxin Yang
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Jiawen He
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Wei Wei
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Weiwei Chu
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Shigang Yu
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Ye Tian
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Fengyi Peng
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Hongcheng Liu
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Zengkai Zhang
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China
| | - Jie Chen
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095, PR China.
| |
Collapse
|
20
|
Zhang Y, Sun T, Liu S, Dong L, Liu C, Song W, Liu J, Gai S. MYC cis-Elements in PsMPT Promoter Is Involved in Chilling Response of Paeonia suffruticosa. PLoS One 2016; 11:e0155780. [PMID: 27228117 PMCID: PMC4882030 DOI: 10.1371/journal.pone.0155780] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2016] [Accepted: 05/04/2016] [Indexed: 11/25/2022] Open
Abstract
The MPT transports Pi to synthesize ATP. PsMPT, a chilling-induced gene, was previously reported to promote energy metabolism during bud dormancy release in tree peony. In this study, the regulatory elements of PsMPT promoter involved in chilling response were further analyzed. The PsMPT transcript was detected in different tree peony tissues and was highly expressed in the flower organs, including petal, stigma and stamen. An 1174 bp of the PsMPT promoter was isolated by TAIL-PCR, and the PsMPT promoter::GUS transgenic Arabidopsis was generated and analyzed. GUS staining and qPCR showed that the promoter was active in mainly the flower stigma and stamen. Moreover, it was found that the promoter activity was enhanced by chilling, NaCl, GA, ACC and NAA, but inhibited by ABA, mannitol and PEG. In transgenic plants harboring 421 bp of the PsMPT promoter, the GUS gene expression and the activity were significantly increased by chilling treatment. When the fragment from -421 to -408 containing a MYC cis-element was deleted, the chilling response could not be observed. Further mutation analysis confirmed that the MYC element was one of the key motifs responding to chilling in the PsMPT promoter. The present study provides useful information for further investigation of the regulatory mechanism of PsMPT during the endo-dormancy release.
Collapse
Affiliation(s)
- Yuxi Zhang
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
| | - Tingzhao Sun
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
| | - Shaoqing Liu
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
| | - Lei Dong
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
| | - Chunying Liu
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
| | - Wenwen Song
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
| | - Jingjing Liu
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
| | - Shupeng Gai
- College of Life Sciences, Qingdao Agricultural University, Key Lab of Plant Biotechnology in Universities of Shandong Province, Changcheng Road 700, Qingdao, China
- * E-mail:
| |
Collapse
|
21
|
Abstract
As a species, we possess unique biological features that distinguish us from other primates. Here, we review recent efforts to identify changes in gene regulation that drove the evolution of novel human phenotypes. We discuss genotype-directed comparisons of human and nonhuman primate genomes to identify human-specific genetic changes that may encode new regulatory functions. We also review phenotype-directed approaches, which use comparisons of gene expression or regulatory function in homologous human and nonhuman primate cells and tissues to identify changes in expression levels or regulatory activity that may be due to genetic changes in humans. Together, these studies are beginning to reveal the landscape of regulatory innovation in human evolution and point to specific regulatory changes for further study. Finally, we highlight two novel strategies to model human-specific regulatory functions in vivo: primate induced pluripotent stem cells and the generation of humanized mice by genome editing.
Collapse
Affiliation(s)
- Steven K Reilly
- Department of Genetics, Yale School of Medicine, New Haven, Connecticut 06510;
| | - James P Noonan
- Department of Genetics, Yale School of Medicine, New Haven, Connecticut 06510; .,Department of Ecology and Evolutionary Biology, Yale University, New Haven, Connecticut 06511.,Kavli Institute for Neuroscience, Yale School of Medicine, New Haven, Connecticut 06510
| |
Collapse
|
22
|
Young RS, Hayashizaki Y, Andersson R, Sandelin A, Kawaji H, Itoh M, Lassmann T, Carninci P, Bickmore WA, Forrest AR, Taylor MS. The frequent evolutionary birth and death of functional promoters in mouse and human. Genome Res 2015; 25:1546-57. [PMID: 26228054 PMCID: PMC4579340 DOI: 10.1101/gr.190546.115] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2015] [Accepted: 07/28/2015] [Indexed: 12/04/2022]
Abstract
Promoters are central to the regulation of gene expression. Changes in gene regulation are thought to underlie much of the adaptive diversification between species and phenotypic variation within populations. In contrast to earlier work emphasizing the importance of enhancer evolution and subtle sequence changes at promoters, we show that dramatic changes such as the complete gain and loss (collectively, turnover) of functional promoters are common. Using quantitative measures of transcription initiation in both humans and mice across 52 matched tissues, we discriminate promoter sequence gains from losses and resolve the lineage of changes. We also identify expression divergence and functional turnover between orthologous promoters, finding only the latter is associated with local sequence changes. Promoter turnover has occurred at the majority (>56%) of protein-coding genes since humans and mice diverged. Tissue-restricted promoters are the most evolutionarily volatile where retrotransposition is an important, but not the sole, source of innovation. There is considerable heterogeneity of turnover rates between promoters in different tissues, but the consistency of these in both lineages suggests that the same biological systems are similarly inclined to transcriptional rewiring. The genes affected by promoter turnover show evidence of adaptive evolution. In mice, promoters are primarily lost through deletion of the promoter containing sequence, whereas in humans, many promoters appear to be gradually decaying with weak transcriptional output and relaxed selective constraint. Our results suggest that promoter gain and loss is an important process in the evolutionary rewiring of gene regulation and may be a significant source of phenotypic diversification.
Collapse
Affiliation(s)
- Robert S Young
- MRC Human Genetics Unit, MRC Institute for Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| | - Yoshihide Hayashizaki
- RIKEN Preventive Medicine and Diagnosis Innovation Program, Wako, Saitama, 351-0198, Japan
| | - Robin Andersson
- Department of Biology and Biotech Research and Innovation Centre, Copenhagen University, 2200 Copenhagen N, Denmark
| | - Albin Sandelin
- Department of Biology and Biotech Research and Innovation Centre, Copenhagen University, 2200 Copenhagen N, Denmark
| | - Hideya Kawaji
- RIKEN Preventive Medicine and Diagnosis Innovation Program, Wako, Saitama, 351-0198, Japan; RIKEN Center for Life Science Technologies, Division of Genomic Technologies, Tsurumi-ku, Yokohama, 230-0045, Japan
| | - Masayoshi Itoh
- RIKEN Preventive Medicine and Diagnosis Innovation Program, Wako, Saitama, 351-0198, Japan; RIKEN Center for Life Science Technologies, Division of Genomic Technologies, Tsurumi-ku, Yokohama, 230-0045, Japan
| | - Timo Lassmann
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, Tsurumi-ku, Yokohama, 230-0045, Japan
| | - Piero Carninci
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, Tsurumi-ku, Yokohama, 230-0045, Japan
| | | | - Wendy A Bickmore
- MRC Human Genetics Unit, MRC Institute for Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| | - Alistair R Forrest
- RIKEN Center for Life Science Technologies, Division of Genomic Technologies, Tsurumi-ku, Yokohama, 230-0045, Japan; Systems Biology and Genomics, Harry Perkins Institute of Medical Research, QEII Medical Centre, Nedlands, Western Australia 6009, Australia
| | - Martin S Taylor
- MRC Human Genetics Unit, MRC Institute for Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| |
Collapse
|
23
|
Kashkin K, Chernov I, Stukacheva E, Monastyrskaya G, Uspenskaya N, Kopantzev E, Sverdlov E. Cancer specificity of promoters of the genes controlling cell proliferation. J Cell Biochem 2015; 116:299-309. [PMID: 25187488 DOI: 10.1002/jcb.24968] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2014] [Accepted: 08/29/2014] [Indexed: 12/20/2022]
Abstract
Violation of proliferation control is a common feature of cancer cells. We put forward the hypothesis that promoters of genes involved in the control of cell proliferation should possess intrinsic cancer specific activity. We cloned promoter regions of CDC6, POLD1, CKS1B, MCM2, and PLK1 genes into pGL3 reporter vector and studied their ability to drive heterologous gene expression in transfected cancer cells of different origin and in normal human fibroblasts. Each promoter was cloned in short (335-800 bp) and long (up to 2.3 kb) variants to cover probable location of core and whole promoter regulatory elements. Cloned promoters were significantly more active in cancer cells than in normal fibroblasts that may indicate their cancer specificity. Both versions of CDC6 promoters were shown to be most active while the activities of others were close to that of BIRC5 gene (survivin) gene promoter. Long and short variants of each cloned promoter demonstrated very similar cancer specificity with the exception of PLK1-long promoter that was substantially more specific than its short variant and other promoters under study. The data indicate that most of the important cis-regulatory transcription elements responsible for intrinsic cancer specificity are located in short variants of the promoters under study. CDC6 short promoter may serve as a promising candidate for transcription targeted cancer gene therapy.
Collapse
Affiliation(s)
- Kirill Kashkin
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, ul. Miklukho-Maklaya 16/10, Moscow, 117997, Russia
| | | | | | | | | | | | | |
Collapse
|
24
|
Scala G, Affinito O, Miele G, Monticelli A, Cocozza S. Evidence for evolutionary and nonevolutionary forces shaping the distribution of human genetic variants near transcription start sites. PLoS One 2014; 9:e114432. [PMID: 25474578 PMCID: PMC4256220 DOI: 10.1371/journal.pone.0114432] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2014] [Accepted: 11/09/2014] [Indexed: 11/19/2022] Open
Abstract
The regions surrounding transcription start sites (TSSs) of genes play a critical role in the regulation of gene expression. At the same time, current evidence indicates that these regions are particularly stressed by transcription-related mutagenic phenomena. In this work we performed a genome-wide analysis of the distribution of single nucleotide polymorphisms (SNPs) inside the 10 kb region flanking human TSSs by dividing SNPs into four classes according to their frequency (rare, two intermediate classes, and common). We found that, in this 10 kb region, the distribution of variants depends on their frequency and on their localization relative to the TSS. We found that the distribution of variants is generally different for TSSs located inside or outside of CpG islands. We found a significant relationship between the distribution of rare variants and nucleosome occupancy scores. Furthermore, our analysis suggests that evolutionary (purifying selection) and nonevolutionary (biased gene conversion) forces both play a role in determining the relative SNP frequency around TSSs. Finally, we analyzed the potential pathogenicity of each class of variant using the Combined Annotation Dependent Depletion score. In conclusion, this study provides a novel and detailed view of the distribution of genomic variants around TSSs, providing insight into the forces that instigate and maintain variability in such critical regions.
Collapse
Affiliation(s)
- Giovanni Scala
- Gruppo Interdipartimentale di Bioinformatica e Biologia Computazionale, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Dipartimento di Fisica, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Naples, Italy
- * E-mail:
| | - Ornella Affinito
- Gruppo Interdipartimentale di Bioinformatica e Biologia Computazionale, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Dipartimento di Medicina Molecolare e Biotecnologie Mediche, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Istituto di Endocrinologia ed Oncologia Sperimentale (IEOS), CNR, Naples, Italy
| | - Gennaro Miele
- Gruppo Interdipartimentale di Bioinformatica e Biologia Computazionale, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Dipartimento di Fisica, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Naples, Italy
| | - Antonella Monticelli
- Gruppo Interdipartimentale di Bioinformatica e Biologia Computazionale, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Istituto di Endocrinologia ed Oncologia Sperimentale (IEOS), CNR, Naples, Italy
| | - Sergio Cocozza
- Gruppo Interdipartimentale di Bioinformatica e Biologia Computazionale, Università degli Studi di Napoli “Federico II”, Naples, Italy
- Dipartimento di Medicina Molecolare e Biotecnologie Mediche, Università degli Studi di Napoli “Federico II”, Naples, Italy
| |
Collapse
|
25
|
Clark BS, Blackshaw S. Long non-coding RNA-dependent transcriptional regulation in neuronal development and disease. Front Genet 2014; 5:164. [PMID: 24936207 PMCID: PMC4047558 DOI: 10.3389/fgene.2014.00164] [Citation(s) in RCA: 118] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2014] [Accepted: 05/18/2014] [Indexed: 01/17/2023] Open
Abstract
Comprehensive analysis of the mammalian transcriptome has revealed that long non-coding RNAs (lncRNAs) may make up a large fraction of cellular transcripts. Recent years have seen a surge of studies aimed at functionally characterizing the role of lncRNAs in development and disease. In this review, we discuss new findings implicating lncRNAs in controlling development of the central nervous system (CNS). The evolution of the higher vertebrate brain has been accompanied by an increase in the levels and complexities of lncRNAs expressed within the developing nervous system. Although a limited number of CNS-expressed lncRNAs are now known to modulate the activity of proteins important for neuronal differentiation, the function of the vast majority of neuronal-expressed lncRNAs is still unknown. Topics of intense current interest include the mechanism by which CNS-expressed lncRNAs might function in epigenetic and transcriptional regulation during neuronal development, and how gain and loss of function of individual lncRNAs contribute to neurological diseases.
Collapse
Affiliation(s)
- Brian S Clark
- Solomon Snyder Department of Neuroscience, Johns Hopkins University School of Medicine Baltimore, MD, USA
| | - Seth Blackshaw
- Solomon Snyder Department of Neuroscience, Johns Hopkins University School of Medicine Baltimore, MD, USA ; Department of Ophthalmology, Johns Hopkins University School of Medicine Baltimore, MD, USA ; Department of Neurology, Johns Hopkins University School of Medicine Baltimore, MD, USA ; Center for High-Throughput Biology, Johns Hopkins University School of Medicine Baltimore, MD, USA ; Institute for Cell Engineering, Johns Hopkins University School of Medicine Baltimore, MD, USA
| |
Collapse
|
26
|
Talhaoui I, Couve S, Gros L, Ishchenko AA, Matkarimov B, Saparbaev MK. Aberrant repair initiated by mismatch-specific thymine-DNA glycosylases provides a mechanism for the mutational bias observed in CpG islands. Nucleic Acids Res 2014; 42:6300-13. [PMID: 24692658 PMCID: PMC4041421 DOI: 10.1093/nar/gku246] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2013] [Revised: 03/09/2014] [Accepted: 03/13/2014] [Indexed: 12/13/2022] Open
Abstract
The human thymine-DNA glycosylase (TDG) initiates the base excision repair (BER) pathway to remove spontaneous and induced DNA base damage. It was first biochemically characterized for its ability to remove T mispaired with G in CpG context. TDG is involved in the epigenetic regulation of gene expressions by protecting CpG-rich promoters from de novo DNA methylation. Here we demonstrate that TDG initiates aberrant repair by excising T when it is paired with a damaged adenine residue in DNA duplex. TDG targets the non-damaged DNA strand and efficiently excises T opposite of hypoxanthine (Hx), 1,N(6)-ethenoadenine, 7,8-dihydro-8-oxoadenine and abasic site in TpG/CpX context, where X is a modified residue. In vitro reconstitution of BER with duplex DNA containing Hx•T pair and TDG results in incorporation of cytosine across Hx. Furthermore, analysis of the mutation spectra inferred from single nucleotide polymorphisms in human population revealed a highly biased mutation pattern within CpG islands (CGIs), with enhanced mutation rate at CpA and TpG sites. These findings demonstrate that under experimental conditions used TDG catalyzes sequence context-dependent aberrant removal of thymine, which results in TpG, CpA→CpG mutations, thus providing a plausible mechanism for the putative evolutionary origin of the CGIs in mammalian genomes.
Collapse
Affiliation(s)
- Ibtissam Talhaoui
- Groupe Réparation de l’ADN, Université Paris Sud, Laboratoire Stabilité Génétique et Oncogenèse CNRS, UMR 8200, Gustave Roussy, F-94805 Villejuif Cedex, France
| | - Sophie Couve
- Laboratoire de Génétique Oncologique EPHE, INSERM U753, Gustave Roussy, F-94805 Villejuif, France
| | - Laurent Gros
- Groupe Réparation de l’ADN, Université Paris Sud, Laboratoire Stabilité Génétique et Oncogenèse CNRS, UMR 8200, Gustave Roussy, F-94805 Villejuif Cedex, France
- AB Science SA, 75008 Paris, France
| | - Alexander A. Ishchenko
- Groupe Réparation de l’ADN, Université Paris Sud, Laboratoire Stabilité Génétique et Oncogenèse CNRS, UMR 8200, Gustave Roussy, F-94805 Villejuif Cedex, France
| | - Bakhyt Matkarimov
- Nazarbayev University Research and Innovation System, Astana 010000, Kazakhstan
| | - Murat K. Saparbaev
- Groupe Réparation de l’ADN, Université Paris Sud, Laboratoire Stabilité Génétique et Oncogenèse CNRS, UMR 8200, Gustave Roussy, F-94805 Villejuif Cedex, France
| |
Collapse
|
27
|
Abstract
Evolutionary conservation has been an accurate predictor of functional elements across the first decade of metazoan genomics. More recently, there has been a move to define functional elements instead from biochemical annotations. Evolutionary methods are, however, more comprehensive than biochemical approaches can be and can assess quantitatively, especially for subtle effects, how biologically important--how injurious after mutation--different types of elements are. Evolutionary methods are thus critical for understanding the large fraction (up to 10%) of the human genome that does not encode proteins and yet might convey function. These methods can also capture the ephemeral nature of much noncoding functional sequence, with large numbers of functional elements having been gained and lost rapidly along each mammalian lineage. Here, we review how different strengths of purifying selection have impacted on protein-coding and non-protein-coding loci and on transcription factor binding sites in mammalian and fruit fly genomes.
Collapse
Affiliation(s)
- Wilfried Haerty
- MRC Functional Genomics Unit, Department of Physiology, Anatomy, and Genetics, University of Oxford, Oxford OX1 3PT, United Kingdom; ,
| | | |
Collapse
|
28
|
Huang S, Li J, Xu A, Huang G, You L. Small Insertions Are More Deleterious than Small Deletions in Human Genomes. Hum Mutat 2013; 34:1642-9. [DOI: 10.1002/humu.22435] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2013] [Accepted: 08/22/2013] [Indexed: 11/09/2022]
Affiliation(s)
- Shengfeng Huang
- State Key Laboratory of Biocontrol; Guangdong Key Laboratory of Pharmaceutical Functional Genes; College of Life Sciences, Sun Yat-sen University; Guangzhou 510275 People's Republic of China
| | - Jie Li
- State Key Laboratory of Biocontrol; Guangdong Key Laboratory of Pharmaceutical Functional Genes; College of Life Sciences, Sun Yat-sen University; Guangzhou 510275 People's Republic of China
| | - Anlong Xu
- State Key Laboratory of Biocontrol; Guangdong Key Laboratory of Pharmaceutical Functional Genes; College of Life Sciences, Sun Yat-sen University; Guangzhou 510275 People's Republic of China
- Beijing University of Chinese Medicine, Chao-yang District; Beijing 100029 People's Republic of China
| | - Guangrui Huang
- State Key Laboratory of Biocontrol; Guangdong Key Laboratory of Pharmaceutical Functional Genes; College of Life Sciences, Sun Yat-sen University; Guangzhou 510275 People's Republic of China
| | - Leiming You
- State Key Laboratory of Biocontrol; Guangdong Key Laboratory of Pharmaceutical Functional Genes; College of Life Sciences, Sun Yat-sen University; Guangzhou 510275 People's Republic of China
| |
Collapse
|
29
|
Abstract
Recently articles have been published disputing the main finding of the ENCODE project that the majority of the human genome exhibits biochemical indices of function, based primarily on low sequence conservation and the existence of larger genomes in some ostensibly simpler organisms (the C-value enigma), indicating the likely presence of significant amounts of junk. Here we challenge these arguments, showing that conservation is a relative measure based on circular assumptions of the non-functionality of transposon-derived sequences and uncertain comparison sets, and that regulatory sequence evolution is subject to different and much more plastic structure-function constraints than protein-coding sequences, as well as positive selection for adaptive radiation. We also show that polyploidy accounts for the higher than expected genome sizes in some eukaryotes, compounded by variable levels of repetitive sequences of unknown significance. We argue that the extent of precise dynamic and differential cell- and tissue-specific transcription and splicing observed from the majority of the human genome is a more reliable indicator of genetic function than conservation, although the unexpectedly large amount of regulatory RNA presents a conceptual challenge to the traditional protein-centric view of human genetic programming. Finally, we suggest that resistance to these findings is further motivated in some quarters by the use of the dubious concept of junk DNA as evidence against intelligent design.
Collapse
|
30
|
Qu W, Hashimoto SI, Shimada A, Nakatani Y, Ichikawa K, Saito TL, Ogoshi K, Matsushima K, Suzuki Y, Sugano S, Takeda H, Morishita S. Genome-wide genetic variations are highly correlated with proximal DNA methylation patterns. Genome Res 2012; 22:1419-25. [PMID: 22689467 PMCID: PMC3409255 DOI: 10.1101/gr.140236.112] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
5-methyl-cytosines at CpG sites frequently mutate into thymines, accounting for a large proportion of spontaneous point mutations. The repair system would leave substantial numbers of errors in neighboring regions if the synthesis of erased gaps around deaminated 5-methyl-cytosines is error-prone. Indeed, we identified an unexpected genome-wide role of the CpG methylation state as a major determinant of proximal natural genetic variation. Specifically, 507 Mbp (∼18%) of the human genome was within 10 bp of a CpG site; in these regions, the single nucleotide polymorphism (SNP) rate significantly increased by ∼50% (P < 10−566 by a two-proportion z-test) if the neighboring CpG sites are methylated. To reconfirm this finding in another vertebrate, we compared six single-base resolution methylomes in two inbred medaka (Oryzias latipes) strains with sufficient genetic divergence (3.4%). We found that the SNP rate also increased by ∼50% (P < 10−2170), and the substitution rates in all dinucleotides increased simultaneously (P < 10−441) around methylated CpG sites. In the hypomethylated regions, the “CGCG” motif was significantly enriched (P < 10−680) and evolutionarily conserved (P = ∼ 0.203%), and slow CpG deamination rather than fast CpG gain was seen, indicating a possible role of CGCG as a candidate cis-element for the hypomethylation state. In regions that were hypermethylated in germline-like tissues but were hypomethylated in somatic liver cells, the SNP rate was significantly smaller than that in hypomethylated regions in both tissues, suggesting a positive selective pressure during DNA methylation reprogramming. This is the first report of findings showing that the CpG methylation state is significantly correlated with the characteristics of evolutionary change in neighboring DNA.
Collapse
Affiliation(s)
- Wei Qu
- Department of Computational Biology, Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa 277-0882, Japan.
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
31
|
Conservation and divergence in Toll-like receptor 4-regulated gene expression in primary human versus mouse macrophages. Proc Natl Acad Sci U S A 2012; 109:E944-53. [PMID: 22451944 DOI: 10.1073/pnas.1110156109] [Citation(s) in RCA: 235] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Evolutionary change in gene expression is generally considered to be a major driver of phenotypic differences between species. We investigated innate immune diversification by analyzing interspecies differences in the transcriptional responses of primary human and mouse macrophages to the Toll-like receptor (TLR)-4 agonist lipopolysaccharide (LPS). By using a custom platform permitting cross-species interrogation coupled with deep sequencing of mRNA 5' ends, we identified extensive divergence in LPS-regulated orthologous gene expression between humans and mice (24% of orthologues were identified as "divergently regulated"). We further demonstrate concordant regulation of human-specific LPS target genes in primary pig macrophages. Divergently regulated orthologues were enriched for genes encoding cellular "inputs" such as cell surface receptors (e.g., TLR6, IL-7Rα) and functional "outputs" such as inflammatory cytokines/chemokines (e.g., CCL20, CXCL13). Conversely, intracellular signaling components linking inputs to outputs were typically concordantly regulated. Functional consequences of divergent gene regulation were confirmed by showing LPS pretreatment boosts subsequent TLR6 responses in mouse but not human macrophages, in keeping with mouse-specific TLR6 induction. Divergently regulated genes were associated with a large dynamic range of gene expression, and specific promoter architectural features (TATA box enrichment, CpG island depletion). Surprisingly, regulatory divergence was also associated with enhanced interspecies promoter conservation. Thus, the genes controlled by complex, highly conserved promoters that facilitate dynamic regulation are also the most susceptible to evolutionary change.
Collapse
|
32
|
Abstract
Copy number variation has recently received considerable attention, and copy number variants (CNVs) have been shown to be both common in mammalian genomes and important for understanding genetic and phenotypic variation. As empirical knowledge and detection methods are quickly advancing, evolutionary theories about CNVs are rapidly updated and often revised. Here, we review recent progress on understanding CNVs, and we discuss some key issues for future research. In essence, we discuss four major forces in population genetics, recombination, mutation, selection, and demography, in relation to CNVs.
Collapse
|
33
|
Zhao X, Valen E, Parker BJ, Sandelin A. Systematic clustering of transcription start site landscapes. PLoS One 2011; 6:e23409. [PMID: 21887249 PMCID: PMC3160847 DOI: 10.1371/journal.pone.0023409] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2011] [Accepted: 07/15/2011] [Indexed: 12/27/2022] Open
Abstract
Genome-wide, high-throughput methods for transcription start site (TSS) detection have shown that most promoters have an array of neighboring TSSs where some are used more than others, forming a distribution of initiation propensities. TSS distributions (TSSDs) vary widely between promoters and earlier studies have shown that the TSSDs have biological implications in both regulation and function. However, no systematic study has been made to explore how many types of TSSDs and by extension core promoters exist and to understand which biological features distinguish them. In this study, we developed a new non-parametric dissimilarity measure and clustering approach to explore the similarities and stabilities of clusters of TSSDs. Previous studies have used arbitrary thresholds to arrive at two general classes: broad and sharp. We demonstrated that in addition to the previous broad/sharp dichotomy an additional category of promoters exists. Unlike typical TATA-driven sharp TSSDs where the TSS position can vary a few nucleotides, in this category virtually all TSSs originate from the same genomic position. These promoters lack epigenetic signatures of typical mRNA promoters and a substantial subset of them are mapping upstream of ribosomal protein pseudogenes. We present evidence that these are likely mapping errors, which have confounded earlier analyses, due to the high similarity of ribosomal gene promoters in combination with known G addition bias in the CAGE libraries. Thus, previous two-class separations of promoter based on TSS distributions are motivated, but the ultra-sharp TSS distributions will confound downstream analyses if not removed.
Collapse
Affiliation(s)
- Xiaobei Zhao
- Department of Biology and Biotech Research and Innovation Centre, The Bioinformatics Centre, Copenhagen University, Copenhagen, Denmark
| | - Eivind Valen
- Department of Biology and Biotech Research and Innovation Centre, The Bioinformatics Centre, Copenhagen University, Copenhagen, Denmark
| | - Brian J. Parker
- Department of Biology and Biotech Research and Innovation Centre, The Bioinformatics Centre, Copenhagen University, Copenhagen, Denmark
| | - Albin Sandelin
- Department of Biology and Biotech Research and Innovation Centre, The Bioinformatics Centre, Copenhagen University, Copenhagen, Denmark
| |
Collapse
|
34
|
Variation in the IGF2 gene promoter region is associated with intramuscular fat content in porcine skeletal muscle. Mol Biol Rep 2011; 39:4101-10. [PMID: 21779802 DOI: 10.1007/s11033-011-1192-5] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2011] [Accepted: 07/11/2011] [Indexed: 12/16/2022]
Abstract
Intramuscular fat (IMF) and subcutaneous fat (back fat-BF) are two of the major fat depots in livestock. A QTN located in the insulin-like growth factor 2 gene (IGF2) has been associated with a desirable reduction in BF depth in pigs. Given that the lipid metabolism of intramuscular adipocytes differs from that of subcutaneous fat adipocytes, this study aimed to search for genetic variation in the IGF2 gene that may be associated with IMF, as well as BF, in diverse pig breeds. Four proximal promoter regions of the IGF2 gene were characterised and the association of IGF2 genetic variation with IMF and BF was assessed. Six promoter SNPs were identified in four promoter regions (P1-P4; sequence coverage 945, 866, 784 and 864 bp, respectively) in phenotypically diverse F1 cross populations. Three promoter SNPs were subsequently genotyped in three pure breeds (Pietrain = 98, Duroc = 99 and Large White = 98). All three SNPs were >95% monomorphic in the Pietrain and Duroc breeds but minor alleles were at moderate frequencies in the Large White breed. These SNPs were linked and one was located in a putative transcription factor binding site. Five haplotypes were inferred and three combined diplotypes tested for association with IMF and BF in the Large White. As expected haplotype 1 (likely in LD with the beneficial QTN allele) was superior for BF level. In contrast, the heterozygote diplotype of the most common haplotypes (1 and 2) was associated with higher IMF and marbling scores compared to either homozygote. Gene expression analysis of divergent animals showed that IGF2 was 1.89 fold up-regulated in muscle with higher compared to lower IMF content. These findings suggest that genetic variation in the promoter region of the IGF2 gene is associated with IMF content in porcine skeletal muscle and that greater expression of the IGF2 gene is associated with higher IMF content.
Collapse
|
35
|
Association between promoter polymorphisms in a key cytoskeletal gene (Ankyrin 1) and intramuscular fat and water-holding capacity in porcine muscle. Mol Biol Rep 2011; 39:3903-14. [DOI: 10.1007/s11033-011-1169-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2010] [Accepted: 06/30/2011] [Indexed: 10/18/2022]
|
36
|
Yamamoto YY, Yoshioka Y, Hyakumachi M, Obokata J. Characteristics of core promoter types with respect to gene structure and expression in Arabidopsis thaliana. DNA Res 2011; 18:333-42. [PMID: 21745829 PMCID: PMC3190954 DOI: 10.1093/dnares/dsr020] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It is now well known that vertebrates use multiple types of core promoter to accomplish differentiated tasks in Pol II-dependent transcription. Several transcriptional characteristics are known to be associated with core types, including distribution patterns of transcription start sites (TSSs) and selection between tissue-specific and constitutive expression profiles. However, their relationship to gene structure is poorly understood. In this report, we carried a comparative analysis of three Arabidopsis core types, TATA, GA, and Coreless, with regard to gene structure. Our genome-wide investigation was based on the peak TSS positions in promoters that had been identified in a large-scale experimental analysis. This analysis revealed that the types of core promoter are related with the room for promoters that is measured as the distance from the TSS to the end of the upstream gene, the distance from the TSS to the start position of the coding sequence (CDS), and the number and species of the cis-regulatory elements. Of these, it was found that the distance from the TSS to the CDS has a tight, inverse correlation to the expression level, and thus the observed relationship to the core type appears to be indirect. However, promoter length and preference of cis-elements are thought to be a direct reflection of core type-specific transcriptional initiation mechanisms.
Collapse
Affiliation(s)
- Yoshiharu Y Yamamoto
- Faculty of Applied Biological Sciences, Gifu University, Yanagido 1-1, Gifu City, Gifu 501-1193, Japan.
| | | | | | | |
Collapse
|
37
|
Behrens S, Vingron M. Studying the evolution of promoter sequences: a waiting time problem. J Comput Biol 2011; 17:1591-606. [PMID: 21128851 DOI: 10.1089/cmb.2010.0084] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
To gain a better understanding of the evolutionary dynamics of regulatory DNA sequences, we address the following questions: (1) How long does it take until a given transcription factor (TF) binding site emerges at random in a promoter sequence? and (2) How does the composition of a TF binding site affect this waiting time? Using two different probabilistic models (an i.i.d. model and a neighbor dependent model), we can compute the expected waiting time for every k-mer, k ranging from 5 to 10, until it appears in a promoter of a species. Our findings indicate that new TF binding sites can be created on a short evolutionary time scale, i.e. in a time span below the speciation time of human and chimp. Furthermore, one can conclude that the composition of a TF binding site plays a crucial role concerning the waiting time until it appears and that the CpG methylation-deamination substitution process probably accelerates the creation of new TF binding sites. A screening of existing TF binding sites moreover reveals that k-mers predicted to have short waiting times occur more frequently than others. Supplementary Material is available at www.libertonline.com/cmb .
Collapse
Affiliation(s)
- Sarah Behrens
- Computational Molecular Biology, Max Planck Institute for Molecular Genetics, Berlin, Germany.
| | | |
Collapse
|
38
|
Bullaughey K. Changes in selective effects over time facilitate turnover of enhancer sequences. Genetics 2011; 187:567-82. [PMID: 21098721 PMCID: PMC3030497 DOI: 10.1534/genetics.110.121590] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2010] [Accepted: 11/10/2010] [Indexed: 11/18/2022] Open
Abstract
Correct gene expression is often critical and consequently stabilizing selection on expression is widespread. Yet few genes possess highly conserved regulatory DNA, and for the few enhancers that have been carefully characterized, substantial functional reorganization has often occurred. Given that natural selection removes mutations of even very small deleterious effect, how can transcription factor binding evolve so readily when it underlies a conserved phenotype? As a first step toward addressing this question, I combine a computational model for regulatory function that incorporates many aspects of our present biological knowledge with a model for the fitness effects of misexpression. I then use this model to study the evolution of enhancers. Several robust behaviors emerge: First, the selective effects of mutations at a site change dramatically over time due to substitutions elsewhere in the enhancer, and even the overall degree of constraint across the enhancer can change considerably. Second, many of the substitutions responsible for changes in binding occur at sites where previously the mutation would have been strongly deleterious, suggesting that fluctuations in selective effects at a site are important for functional turnover. Third, most substitutions contributing to the repatterning of binding and constraint are effectively neutral, highlighting the importance of genetic drift-even for enhancers underlying conserved phenotypes. These findings have important implications for phylogenetic inference of function and for interpretations of selection coefficients estimated for regulatory DNA.
Collapse
Affiliation(s)
- Kevin Bullaughey
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois 60637, USA.
| |
Collapse
|
39
|
Harris EE. Nonadaptive processes in primate and human evolution. AMERICAN JOURNAL OF PHYSICAL ANTHROPOLOGY 2011; 143 Suppl 51:13-45. [PMID: 21086525 DOI: 10.1002/ajpa.21439] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Evolutionary biology has tended to focus on adaptive evolution by positive selection as the primum mobile of evolutionary trajectories in species while underestimating the importance of nonadaptive evolutionary processes. In this review, I describe evidence that suggests that primate and human evolution has been strongly influenced by nonadaptive processes, particularly random genetic drift and mutation. This is evidenced by three fundamental effects: a relative relaxation of selective constraints (i.e., purifying selection), a relative increase in the fixation of slightly deleterious mutations, and a general reduction in the efficacy of positive selection. These effects are observed in protein-coding, regulatory regions, and in gene expression data, as well as in an augmentation of fixation of large-scale mutations, including duplicated genes, mobile genetic elements, and nuclear mitochondrial DNA. The evidence suggests a general population-level explanation such as a reduction in effective population size (N(e)). This would have tipped the balance between the evolutionary forces of natural selection and random genetic drift toward genetic drift for variants having small selective effects. After describing these proximate effects, I describe the potential consequences of these effects for primate and human evolution. For example, an increase in the fixation of slightly deleterious mutations could potentially have led to an increase in the fixation rate of compensatory mutations that act to suppress the effects of slightly deleterious substitutions. The potential consequences of compensatory evolution for the evolution of novel gene functions and in potentially confounding the detection of positively selected genes are explored. The consequences of the passive accumulation of large-scale genomic mutations by genetic drift are unclear, though evidence suggests that new gene copies as well as insertions of transposable elements into genes can potentially lead to adaptive phenotypes. Finally, because a decrease in selective constraint at the genetic level is expected to have effects at the morphological level, I review studies that compare rates of morphological change in various mammalian and island populations where N(e) is reduced. Furthermore, I discuss evidence that suggests that craniofacial morphology in the Homo lineage has shifted from an evolutionary rate constrained by purifying selection toward a neutral evolutionary rate.
Collapse
Affiliation(s)
- Eugene E Harris
- Department of Biological Sciences and Geology, Queensborough Community College, City University of New York, Bayside, NY 10364, USA.
| |
Collapse
|
40
|
Abstract
The medaka fish, Oryzias latipes, is an emerging vertebrate model and now has a high quality draft genome and a number of unique mutants. The long history of medaka research in Japan has provided medaka with unique features, which are complementary to other vertebrate models. A large collection of spontaneous mutants collected over a century, the presence of highly polymorphic inbred lines established over decades, and the recently completed genome sequence all give the medaka a big boost. This review focuses on the state of the art in medaka genetics and genomics, such as the first isolation of active transposons in vertebrates, the influence of chromatin structure on sequence variation, fine quantitative trait locus (QTL) analysis, and versatile mutants as human disease models.
Collapse
Affiliation(s)
- Hiroyuki Takeda
- Department of Biological Sciences, Graduate School of Science, University of Tokyo, Tokyo 113-0033, Japan.
| | | |
Collapse
|
41
|
Aslan O, Sweeney T, Mullen AM, Hamill RM. Regulatory polymorphisms in the bovine Ankyrin 1 gene promoter are associated with tenderness and intramuscular fat content. BMC Genet 2010; 11:111. [PMID: 21159195 PMCID: PMC3022666 DOI: 10.1186/1471-2156-11-111] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2010] [Accepted: 12/15/2010] [Indexed: 12/01/2022] Open
Abstract
Background Recent QTL and gene expression studies have highlighted ankyrins as positional and functional candidate genes for meat quality. Our objective was to characterise the promoter region of the bovine ankyrin 1 gene and to test polymorphisms for association with sensory and technological meat quality measures. Results Seven novel promoter SNPs were identified in a 1.11 kb region of the ankyrin 1 promoter in Angus, Charolais and Limousin bulls (n = 15 per breed) as well as 141 crossbred beef animals for which meat quality data was available. Eighteen haplotypes were inferred with significant breed variation in haplotype frequencies. The five most frequent SNPs and the four most frequent haplotypes were subsequently tested for association with sensory and technological measures of meat quality in the crossbred population. SNP1, SNP3 and SNP4 (which were subsequently designated regulatory SNPs) and SNP5 were associated with traits that contribute to sensorial and technological measurements of tenderness and texture; Haplotype 1 and haplotype 4 were oppositely correlated with traits contributing to tenderness (P < 0.05). While no single SNP was associated with intramuscular fat (IMF), a clear association with increased IMF and juiciness was observed for haplotype 2. Conclusion The conclusion from this study is that alleles defining haplotypes 2 and 4 could usefully contribute to marker SNP panels used to select individuals with improved IMF/juiciness or tenderness in a genome-assisted selection framework.
Collapse
Affiliation(s)
- Ozlem Aslan
- Teagasc, Food Research Centre, Ashtown, Dublin, Ireland
| | | | | | | |
Collapse
|
42
|
Sverdlov ED, Vinogradova TV. Core promoters as an example of the effect of whole-genome information on the evolution of views on molecular mechanisms of vital activity. Mol Biol 2010. [DOI: 10.1134/s002689331005002x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
43
|
Charoensawan V, Adryan B, Martin S, Söllner C, Thisse B, Thisse C, Wright GJ, Teichmann SA. The impact of gene expression regulation on evolution of extracellular signaling pathways. Mol Cell Proteomics 2010; 9:2666-77. [PMID: 20935258 PMCID: PMC3101855 DOI: 10.1074/mcp.m110.003020] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Extracellular protein interactions are crucial to the development of multicellular organisms because they initiate signaling pathways and enable cellular recognition cues. Despite their importance, extracellular protein interactions are often under-represented in large scale protein interaction data sets because most high throughput assays are not designed to detect low affinity extracellular interactions. Due to the lack of a comprehensive data set, the evolution of extracellular signaling pathways has remained largely a mystery. We investigated this question using a combined data set of physical pairwise interactions between zebrafish extracellular proteins, mainly from the immunoglobulin superfamily and leucine-rich repeat families, and their spatiotemporal expression profiles. We took advantage of known homology between proteins to estimate the relative rates of changes of four parameters after gene duplication, namely extracellular protein interaction, expression pattern, and the divergence of extracellular and intracellular protein sequences. We showed that change in expression profile is a major contributor to the evolution of signaling pathways followed by divergence in intracellular protein sequence, whereas extracellular sequence and interaction profiles were relatively more conserved. Rapidly evolving expression profiles will eventually drive other parameters to diverge more quickly because differentially expressed proteins get exposed to different environments and potential binding partners. This allows homologous extracellular receptors to attain specialized functions and become specific to tissues and/or developmental stages.
Collapse
Affiliation(s)
- Varodom Charoensawan
- Medical Research Council Laboratory of Molecular Biology, Cambridge CB20QH, United Kingdom.
| | | | | | | | | | | | | | | |
Collapse
|
44
|
Sharif J, Endo TA, Toyoda T, Koseki H. Divergence of CpG island promoters: A consequence or cause of evolution? Dev Growth Differ 2010; 52:545-54. [DOI: 10.1111/j.1440-169x.2010.01193.x] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
45
|
Chodroff RA, Goodstadt L, Sirey TM, Oliver PL, Davies KE, Green ED, Molnár Z, Ponting CP. Long noncoding RNA genes: conservation of sequence and brain expression among diverse amniotes. Genome Biol 2010; 11:R72. [PMID: 20624288 PMCID: PMC2926783 DOI: 10.1186/gb-2010-11-7-r72] [Citation(s) in RCA: 194] [Impact Index Per Article: 13.9] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2010] [Revised: 05/17/2010] [Accepted: 07/12/2010] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Long considered to be the building block of life, it is now apparent that protein is only one of many functional products generated by the eukaryotic genome. Indeed, more of the human genome is transcribed into noncoding sequence than into protein-coding sequence. Nevertheless, whilst we have developed a deep understanding of the relationships between evolutionary constraint and function for protein-coding sequence, little is known about these relationships for non-coding transcribed sequence. This dearth of information is partially attributable to a lack of established non-protein-coding RNA (ncRNA) orthologs among birds and mammals within sequence and expression databases. RESULTS Here, we performed a multi-disciplinary study of four highly conserved and brain-expressed transcripts selected from a list of mouse long intergenic noncoding RNA (lncRNA) loci that generally show pronounced evolutionary constraint within their putative promoter regions and across exon-intron boundaries. We identify some of the first lncRNA orthologs present in birds (chicken), marsupial (opossum), and eutherian mammals (mouse), and investigate whether they exhibit conservation of brain expression. In contrast to conventional protein-coding genes, the sequences, transcriptional start sites, exon structures, and lengths for these non-coding genes are all highly variable. CONCLUSIONS The biological relevance of lncRNAs would be highly questionable if they were limited to closely related phyla. Instead, their preservation across diverse amniotes, their apparent conservation in exon structure, and similarities in their pattern of brain expression during embryonic and early postnatal stages together indicate that these are functional RNA molecules, of which some have roles in vertebrate brain development.
Collapse
Affiliation(s)
- Rebecca A Chodroff
- Department of Physiology, Anatomy, and Genetics, Le Gros Clark Building South Parks Road, University of Oxford, Oxford OX1 3QX, UK
| | | | | | | | | | | | | | | |
Collapse
|
46
|
Jung CH, Makunin IV, Mattick JS. Identification of conserved Drosophila-specific euchromatin-restricted non-coding sequence motifs. Genomics 2010; 96:154-66. [PMID: 20595017 DOI: 10.1016/j.ygeno.2010.05.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2010] [Revised: 05/25/2010] [Accepted: 05/26/2010] [Indexed: 01/19/2023]
Abstract
Non-protein-coding DNA comprises the majority of animal genomes but its functions are largely unknown. We identified over 17,000 different tetranucleotide pairs in the Drosophila melanogaster genome that are over-represented at distances up to 100nt in conserved non-exonic sequences. Those exhibiting the highest information content in surrounding nucleotides were classified into five groups: tRNAs, motifs associated with histone genes, Suppressor-of-Hairy-wing binding sites, and two sets of previously unrecognized motifs (DLM3 and DLM4). There are hundreds to thousands of copies of DLM3 and DLM4, respectively, in the genome, located almost exclusively in non-coding regions. They have similar copy numbers among drosophilids, but are largely absent in other insects. DLM3 is likely a cis-regulatory element, whereas DLM4 sequences are capable of forming a short hairpin structure and are expressed as approximately 80nt RNAs. This work reports the existence of Drosophila genus-specific sequence motifs, and suggests that many more novel functional elements may be discovered in genomes using the general approach outlined herein.
Collapse
Affiliation(s)
- Chol-Hee Jung
- Institute for Molecular Bioscience, The University of Queensland, St Lucia QLD, Australia
| | | | | |
Collapse
|
47
|
Abstract
We have produced an evolutionary model for promoters, analogous to the commonly used synonymous/nonsynonymous mutation models for protein-coding sequences. Although our model, called Sunflower, relies on some simple assumptions, it captures enough of the biology of transcription factor action to show clear correlation with other biological features. Sunflower predicts a binding profile of transcription factors to DNA sequences, in which different factors compete for the same potential binding sites. The parametrized model simultaneously estimates a continuous measurement of binding occupancy across the genomic sequence for each factor. We can then introduce a localized mutation, rerun the binding model, and record the difference in binding profiles. A single mutation can alter interactions both upstream and downstream of its position due to potential overlapping binding sites, and our statistic captures this domino effect. Over evolutionary time, we observe a clear excess of low-scoring mutations fixed in promoters, consistent with most changes being neutral. However, this is not consistent across all promoters, and some promoters show more rapid divergence. This divergence often occurs in the presence of relatively constant protein-coding divergence. Interestingly, different classes of promoters show different sensitivity to mutations, with phosphorylation-related genes having promoters inherently more sensitive to mutations than immune genes. Although there have previously been a number of models attempting to handle transcription factor binding, Sunflower provides a richer biological model, incorporating weak binding sites and the possibility of competition. The results show the first clear correlations between such a model and evolutionary processes.
Collapse
Affiliation(s)
- Michael M Hoffman
- EMBL-European Bioinformatics Institute, Wellcome Trust Genome Campus, Cambridge CB10 1SD, United Kingdom
| | | |
Collapse
|
48
|
Sjödin P, Bataillon T, Schierup MH. Insertion and deletion processes in recent human history. PLoS One 2010; 5:e8650. [PMID: 20098729 PMCID: PMC2808225 DOI: 10.1371/journal.pone.0008650] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2009] [Accepted: 12/14/2009] [Indexed: 11/25/2022] Open
Abstract
Background Although insertions and deletions (indels) account for a sizable portion of genetic changes within and among species, they have received little attention because they are difficult to type, are alignment dependent and their underlying mutational process is poorly understood. A fundamental question in this respect is whether insertions and deletions are governed by similar or different processes and, if so, what these differences are. Methodology/Principal Findings We use published resequencing data from Seattle SNPs and NIEHS human polymorphism databases to construct a genomewide data set of short polymorphic insertions and deletions in the human genome (n = 6228). We contrast these patterns of polymorphism with insertions and deletions fixed in the same regions since the divergence of human and chimpanzee (n = 10546). The macaque genome is used to resolve all indels into insertions and deletions. We find that the ratio of deletions to insertions is greater within humans than between human and chimpanzee. Deletions segregate at lower frequency in humans, providing evidence for deletions being under stronger purifying selection than insertions. The insertion and deletion rates correlate with several genomic features and we find evidence that both insertions and deletions are associated with point mutations. Finally, we find no evidence for a direct effect of the local recombination rate on the insertion and deletion rate. Conclusions/Significance Our data strongly suggest that deletions are more deleterious than insertions but that insertions and deletions are otherwise generally governed by the same genomic factors.
Collapse
Affiliation(s)
- Per Sjödin
- Bioinformatics Research Center, C. F. Møllers Alle, Arhus, Denmark.
| | | | | |
Collapse
|
49
|
Dong X, Akalin A, Sharma Y, Lenhard B. Translog, a web browser for studying the expression divergence of homologous genes. BMC Bioinformatics 2010; 11 Suppl 1:S59. [PMID: 20122234 PMCID: PMC3009532 DOI: 10.1186/1471-2105-11-s1-s59] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
Background Increasing amount of data from comparative genomics, and newly developed technologies producing accurate gene expression data facilitate the study of the expression divergence of homologous genes. Previous studies have individually highlighted factors that contribute to the expression divergence of duplicate genes, e.g. promoter changes, exon structure heterogeneity, asymmetric histone modifications and genomic neighborhood conservation. However, there is a lack of a tool to integrate multiple factors and visualize their variety among homologous genes in a straightforward way. Results We introduce Translog (a web-based tool for Transcriptome comparison of homologous genes) that assists in the comparison of homologous genes by displaying the loci in three different views: promoter view for studying the sharing/turnover of transcription initiations, exon structure for displaying the exon-intron structure changes, and genomic neighborhood to show the macro-synteny conservation in a larger scale. CAGE data for transcription initiation are mapped for each transcript and can be used to study transcription turnover and expression changes. Alignment anchors between homologous loci can be used to define the precise homologous transcripts. We demonstrate how these views can be used to visualize the changes of homologous genes during evolution, particularly after the 2R and 3R whole genome duplication. Conclusion We have developed a web-based tool for assisting in the transcriptome comparison of homologous genes, facilitating the study of expression divergence.
Collapse
|
50
|
Mattick JS. Deconstructing the dogma: a new view of the evolution and genetic programming of complex organisms. Ann N Y Acad Sci 2009; 1178:29-46. [PMID: 19845626 DOI: 10.1111/j.1749-6632.2009.04991.x] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Since the birth of molecular biology it has been generally assumed that most genetic information is transacted by proteins, and that RNA plays an intermediary role. This led to the subsidiary assumption that the vast tracts of noncoding sequences in the genomes of higher organisms are largely nonfunctional, despite the fact that they are transcribed. These assumptions have since become articles of faith, but they are not necessarily correct. I propose an alternative evolutionary history whereby developmental and cognitive complexity has arisen by constructing sophisticated RNA-based regulatory networks that interact with generic effector complexes to control gene expression patterns and the epigenetic trajectories of differentiation and development. Environmental information can also be conveyed into this regulatory system via RNA editing, especially in the brain. Moreover, the observations that RNA-directed epigenetic changes can be inherited raises the intriguing question: has evolution learnt how to learn?
Collapse
Affiliation(s)
- John S Mattick
- Institute for Molecular Bioscience, The University of Queensland, St Lucia QLD 4072, Australia.
| |
Collapse
|