1
|
Vellnow N, Gossmann TI, Waxman D. The pseudoentropy of allele frequency trajectories, the persistence of variation, and the effective population size. Biosystems 2024; 238:105176. [PMID: 38479654 DOI: 10.1016/j.biosystems.2024.105176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Revised: 03/01/2024] [Accepted: 03/01/2024] [Indexed: 03/24/2024]
Abstract
To concisely describe how genetic variation, at individual loci or across whole genomes, changes over time, and to follow transitory allelic changes, we introduce a quantity related to entropy, that we term pseudoentropy. This quantity emerges in a diffusion analysis of the mean time a mutation segregates in a population. For a neutral locus with an arbitrary number of alleles, the mean time of segregation is generally proportional to the pseudoentropy of initial allele frequencies. After the initial time point, pseudoentropy generally decreases, but other behaviours are possible, depending on the genetic diversity and selective forces present. For a biallelic locus, pseudoentropy and entropy coincide, but they are distinct quantities with more than two alleles. Thus for populations with multiple biallelic loci, the language of entropy suffices. Then entropy, combined across loci, serves as a concise description of genetic variation. We used individual based simulations to explore how this entropy behaves under different evolutionary scenarios. In agreement with predictions, the entropy associated with unlinked neutral loci decreases over time. However, deviations from free recombination and neutrality have clear and informative effects on the entropy's behaviour over time. Analysis of publicly available data of a natural D. melanogaster population, that had been sampled over seven years, using a sliding-window approach, yielded considerable variation in entropy trajectories of different genomic regions. These mostly follow a pattern that suggests a substantial effective population size and a limited effect of positive selection on genome-wide diversity over short time scales.
Collapse
Affiliation(s)
- Nikolas Vellnow
- TU Dortmund University, Computational Systems Biology, Faculty of Biochemical and Chemical Engineering, Emil-Figge-Str. 66, 44227 Dortmund, Germany.
| | - Toni I Gossmann
- TU Dortmund University, Computational Systems Biology, Faculty of Biochemical and Chemical Engineering, Emil-Figge-Str. 66, 44227 Dortmund, Germany.
| | - David Waxman
- Fudan University, Centre for Computational Systems Biology, ISTBI, 220 Handan Road, Shanghai 200433, People's Republic of China.
| |
Collapse
|
2
|
Ord J, Gossmann TI, Adrian-Kalchhauser I. High nucleotide diversity accompanies differential DNA methylation in naturally diverging populations. Mol Biol Evol 2023; 40:7083720. [PMID: 36947101 PMCID: PMC10139703 DOI: 10.1093/molbev/msad068] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 02/14/2023] [Accepted: 03/09/2023] [Indexed: 03/23/2023] Open
Abstract
Epigenetic mechanisms such as DNA methylation (DNAme) are thought to comprise an invaluable adaptive toolkit in the early stages of local adaptation, especially when genetic diversity is constrained. However, the link between genetic diversity and DNAme has been scarcely examined in natural populations, despite its potential to shed light on the evolutionary forces acting on methylation state. Here, we analysed reduced-representation bisulfite sequencing and whole genome pool-seq data from marine and freshwater stickleback populations to examine the relationship between DNAme variation (between- and within-population), and nucleotide diversity in the context of freshwater adaptation. We find that sites that are differentially methylated between populations have higher underlying standing genetic variation, with diversity higher among sites that gained methylation in freshwater than those that lost it. Strikingly, while nucleotide diversity is generally lower in the freshwater population as expected from a population bottleneck, this is not the case for sites which lost methylation which instead have elevated nucleotide diversity in freshwater compared to marine. Subsequently, we show that nucleotide diversity is higher among sites with ancestrally variable methylation and also positively correlates with the sensitivity to environmentally induced methylation change. The results suggest that as selection on the control of methylation state becomes relaxed, so too does selection against mutations at the sites themselves. Increased epigenetic variance in a population is therefore likely to precede genetic diversification.
Collapse
Affiliation(s)
- James Ord
- Institute for Fish and Wildlife Health, University of Bern, Länggassstrasse 122, 3012 Bern, Switzerland
- Faculty of Biological and Environmental Sciences, University of Helsinki, Viikinkaari 1, 00790 Helsinki, Finland
| | - Toni I Gossmann
- Computational Systems Biology, Faculty of Biochemical and Chemical Engineering, TU Dortmund University, Emil-Figge-Str. 66, 44227 Dortmund, Germany
| | - Irene Adrian-Kalchhauser
- Institute for Fish and Wildlife Health, University of Bern, Länggassstrasse 122, 3012 Bern, Switzerland
| |
Collapse
|
3
|
Mavreas K, Gossmann TI, Waxman D. Loss and fixation of strongly favoured new variants: Understanding and extending Haldane's result via the Wright-Fisher model. Biosystems 2022; 221:104759. [PMID: 35998748 DOI: 10.1016/j.biosystems.2022.104759] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Revised: 08/04/2022] [Accepted: 08/05/2022] [Indexed: 11/02/2022]
Abstract
The implications of Haldane's analysis for the fixation of beneficial alleles lies at the heart of much of 'population genetic thinking' and underlies many approaches that have been tailored to the detection of positive selection. Within the framework of a branching process, Haldane gave an approximation for the probability that fixation ultimately occurs when the selective advantage of a beneficial allele is small (≪1). Here, we make no use of branching processes. Rather, we work solely within a finite-population Wright-Fisher framework. We use this framework to analyse where Haldane's result applies, and extend Haldane's analysis. In particular, we present results for: (i) the domain of applicability of Haldane's analysis; (ii) the probability that loss occurs up to a given time; (iii) the probabilities that loss and fixation ultimately occur; (iv) an analytic approximation associated with the probability of loss and fixation ultimately occurring; (v) quantification of the crossover from weak to strong selection; (vi) determination of the number of invasive alleles that have a significant probability (>0.95) of invading a novel population. We note that the results obtained for (ii), (iii) and (iv) hold for an arbitrary initial number of mutations, and for selection that can be arbitrarily strong. Our results have fundamental implications for population and conservation genetics, and open up new avenues to identify traces of historically beneficial alleles through comparative genomics.
Collapse
Affiliation(s)
- K Mavreas
- Centre for Computational Systems Biology, ISTBI, Fudan University, 220 Handan Road, Shanghai 200433, China
| | - T I Gossmann
- Department of Evolutionary, Genetics (Animal Behaviour), Bielefeld University, Morgenbreede 45, D-33615 Bielefeld, Germany
| | - D Waxman
- Centre for Computational Systems Biology, ISTBI, Fudan University, 220 Handan Road, Shanghai 200433, China.
| |
Collapse
|
4
|
Vendrami DLJ, Gossmann TI, Chakarov N, Paijmans AJ, Eyre-Walker A, Forcada J, Hoffman JI. Signatures of selection on mitonuclear integrated genes uncover hidden mitogenomic variation in fur seals. Genome Biol Evol 2022; 14:6637498. [PMID: 35809042 PMCID: PMC9338431 DOI: 10.1093/gbe/evac104] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/02/2022] [Indexed: 11/13/2022] Open
Abstract
Nuclear copies of mitochondrial genes (numts) are commonplace in vertebrate genomes and have been characterized in many species. However, relatively little attention has been paid to understanding their evolutionary origins and to disentangling alternative sources of insertions. Numts containing genes with intact mitochondrial reading frames represent good candidates for this purpose. The sequences of the genes they contain can be compared to their mitochondrial homologs to characterize synonymous to non-synonymous substitution rates, which can shed light on the selection pressures these genes have been subjected to. Here, we characterise 25 numts in the Antarctic fur seal (Arctocephalus gazella) genome. Among those containing genes with intact mitochondrial reading frames, three carry multiple substitutions in comparison to their mitochondrial homologs. Our analyses reveal that one represents a historic insertion subjected to strong purifying selection since it colonized the Otarioidea in a genomic region enriched in retrotransposons. By contrast, the other two numts appear to be more recent and their large number of substitutions can be attributed to non-canonical insertions, either the integration of heteroplasmic mtDNA or hybridization. Our study sheds new light on the evolutionary history of pinniped numts and uncovers the presence of hidden sources of mitonuclear variation.
Collapse
Affiliation(s)
- David L J Vendrami
- Department of Animal Behaviour, Bielefeld University, 33501 Bielefeld, Germany
| | - Toni I Gossmann
- Department of Animal Behaviour, Bielefeld University, 33501 Bielefeld, Germany
| | - Nayden Chakarov
- Department of Animal Behaviour, Bielefeld University, 33501 Bielefeld, Germany
| | - Anneke J Paijmans
- Department of Animal Behaviour, Bielefeld University, 33501 Bielefeld, Germany
| | - Adam Eyre-Walker
- School of Life Science, University of Sussex, Brighton, BN1 9QG, UK
| | - Jaume Forcada
- British Antarctic Survey, High Cross, Madingley Road, Cambridge CB3 OET, UK
| | - Joseph I Hoffman
- Department of Animal Behaviour, Bielefeld University, 33501 Bielefeld, Germany.,British Antarctic Survey, High Cross, Madingley Road, Cambridge CB3 OET, UK
| |
Collapse
|
5
|
Gossmann TI, Waxman D. OUP accepted manuscript. Genome Biol Evol 2022; 14:6555663. [PMID: 35349695 PMCID: PMC9016752 DOI: 10.1093/gbe/evac047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/23/2022] [Indexed: 11/30/2022] Open
Abstract
There are many problems in biology and related disciplines involving stochasticity, where a signal can only be detected when it lies above a threshold level, while signals lying below threshold are simply not detected. A consequence is that the detected signal is conditioned to lie above threshold, and is not representative of the actual signal. In this work, we present some general results for the conditioning that occurs due to the existence of such an observational threshold. We show that this conditioning is relevant, for example, to gene-frequency trajectories, where many loci in the genome are simultaneously measured in a given generation. Such a threshold can lead to severe biases of allele frequency estimates under purifying selection. In the analysis presented, within the context of Markov chains such as the Wright–Fisher model, we address two key questions: (1) “What is a natural measure of the strength of the conditioning associated with an observation threshold?” (2) “What is a principled way to correct for the effects of the conditioning?”. We answer the first question in terms of a proportion. Starting with a large number of trajectories, the relevant quantity is the proportion of these trajectories that are above threshold at a later time and hence are detected. The smaller the value of this proportion, the stronger the effects of conditioning. We provide an approximate analytical answer to the second question, that corrects the bias produced by an observation threshold, and performs to reasonable accuracy in the Wright–Fisher model for biologically plausible parameter values.
Collapse
Affiliation(s)
- Toni I. Gossmann
- Department of Evolutionary Genetics, Bielefeld University, Konsequenz 45, 33501 Bielefeld, Germany
- Berlin Institute for Advanced Study, Wallotstrasse 19, 14193 Berlin, Germany
- Corresponding authors: E-mails: ;
| | - David Waxman
- Centre for Computational Systems Biology, ISTBI, Fudan University, 220 Handan Road, Shanghai 20433, People’s Republic of China
- Corresponding authors: E-mails: ;
| |
Collapse
|
6
|
Waterman JM, Gossmann TI, Brandler O, Koprowski JL. Editorial: Ecological, Behavioral and Genomic Consequences in the Rodent Family Sciuridae: Why Are Squirrels So Diverse? Front Ecol Evol 2021. [DOI: 10.3389/fevo.2021.765558] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
|
7
|
Junker N, Gossmann TI. Adaptation-Driven Evolution of Sirtuin 1 (SIRT1), a Key Regulator of Metabolism and Aging, in Marmot Species. Front Ecol Evol 2021. [DOI: 10.3389/fevo.2021.666564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The sirtuin protein family plays a role in the lifespan of various species and is involved in numerous key metabolic processes. To understand the evolutionary role of sirtuins in marmots, a long-living rodent species group with remarkable metabolic shutdown during hibernation, we conducted a phylogeny-based substitution rate analysis of coding genes based on genetic information of seven marmot species. We show that sirtuin 1 (SIRT1) has evolved under positive selection in the marmot lineage. We pinpoint three amino acid changes in four different marmot species that underlie the signal of positive selection and that may favor increased longevity in marmots. Based on a computational structural analysis we can show that all three substitutions affect the secondary structure of the same region in human SIRT1. We propose that the identified region is close to the catalytic domain and that the potential structural changes may impact the catalytic activity of the enzyme and therefore might be playing a functional role in marmot's extended lifespan and metabolic shutdown.
Collapse
|
8
|
Hildebrand F, Gossmann TI, Frioux C, Özkurt E, Myers PN, Ferretti P, Kuhn M, Bahram M, Nielsen HB, Bork P. Dispersal strategies shape persistence and evolution of human gut bacteria. Cell Host Microbe 2021; 29:1167-1176.e9. [PMID: 34111423 PMCID: PMC8288446 DOI: 10.1016/j.chom.2021.05.008] [Citation(s) in RCA: 53] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Revised: 03/19/2021] [Accepted: 05/14/2021] [Indexed: 02/07/2023]
Abstract
Human gut bacterial strains can co-exist with their hosts for decades, but little is known about how these microbes persist and disperse, and evolve thereby. Here, we examined these processes in 5,278 adult and infant fecal metagenomes, longitudinally sampled in individuals and families. Our analyses revealed that a subset of gut species is extremely persistent in individuals, families, and geographic regions, represented often by locally successful strains of the phylum Bacteroidota. These “tenacious” bacteria show high levels of genetic adaptation to the human host but a high probability of loss upon antibiotic interventions. By contrast, heredipersistent bacteria, notably Firmicutes, often rely on dispersal strategies with weak phylogeographic patterns but strong family transmissions, likely related to sporulation. These analyses describe how different dispersal strategies can lead to the long-term persistence of human gut microbes with implications for gut flora modulations. Bacterial strains may persist within family members through transfer Bacteria adapt dispersal strategies: heredipersistent, spatiopersistent, and tenacious Dispersal strategies correlate with genetic bottlenecks and effective population size
Collapse
Affiliation(s)
- Falk Hildebrand
- Gut Microbes and Health, Quadram Institute Bioscience, NR4 7UQ Norwich, UK; Digital Biology, Earlham Institute, NR4 7UZ Norwich, UK; European Molecular Biology Laboratory, Structural and Computational Biology Unit, 69117 Heidelberg, Germany.
| | - Toni I Gossmann
- Department of Animal Behaviour, Bielefeld University, Bielefeld DE-33501, Germany
| | - Clémence Frioux
- Gut Microbes and Health, Quadram Institute Bioscience, NR4 7UQ Norwich, UK; Inria, INRAE, CNRS, Univ. Bordeaux, 33405 Talence, France
| | - Ezgi Özkurt
- Gut Microbes and Health, Quadram Institute Bioscience, NR4 7UQ Norwich, UK; Digital Biology, Earlham Institute, NR4 7UZ Norwich, UK
| | - Pernille Neve Myers
- Clinical Microbiomics A/S, Copenhagen, Denmark; Department of Biotechnology and Biomedicine, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Pamela Ferretti
- European Molecular Biology Laboratory, Structural and Computational Biology Unit, 69117 Heidelberg, Germany
| | - Michael Kuhn
- European Molecular Biology Laboratory, Structural and Computational Biology Unit, 69117 Heidelberg, Germany
| | - Mohammad Bahram
- Department of Ecology, Swedish University of Agricultural Sciences, Ulls väg 16, 750 07 Uppsala, Sweden; Institute of Ecology and Earth Sciences, University of Tartu, Vanemuise 46, 51014 Tartu, Estonia
| | | | - Peer Bork
- European Molecular Biology Laboratory, Structural and Computational Biology Unit, 69117 Heidelberg, Germany; Max Delbrück Center for Molecular Medicine, Berlin, Germany; Yonsei Frontier Lab (YFL), Yonsei University, Seoul 03722, South Korea; Department of Bioinformatics, Biocenter, University of Würzburg, Würzburg, Germany.
| |
Collapse
|
9
|
Yusuf L, Heatley MC, Palmer JPG, Barton HJ, Cooney CR, Gossmann TI. Noncoding regions underpin avian bill shape diversification at macroevolutionary scales. Genome Res 2020; 30:553-565. [PMID: 32269134 PMCID: PMC7197477 DOI: 10.1101/gr.255752.119] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2019] [Accepted: 03/17/2020] [Indexed: 12/18/2022]
Abstract
Recent progress has been made in identifying genomic regions implicated in trait evolution on a microevolutionary scale in many species, but whether these are relevant over macroevolutionary time remains unclear. Here, we directly address this fundamental question using bird beak shape, a key evolutionary innovation linked to patterns of resource use, divergence, and speciation, as a model trait. We integrate class-wide geometric-morphometric analyses with evolutionary sequence analyses of 10,322 protein-coding genes as well as 229,001 genomic regions spanning 72 species. We identify 1434 protein-coding genes and 39,806 noncoding regions for which molecular rates were significantly related to rates of bill shape evolution. We show that homologs of the identified protein-coding genes as well as genes in close proximity to the identified noncoding regions are involved in craniofacial embryo development in mammals. They are associated with embryonic stem cell pathways, including BMP and Wnt signaling, both of which have repeatedly been implicated in the morphological development of avian beaks. This suggests that identifying genotype-phenotype association on a genome-wide scale over macroevolutionary time is feasible. Although the coding and noncoding gene sets are associated with similar pathways, the actual genes are highly distinct, with significantly reduced overlap between them and bill-related phenotype associations specific to noncoding loci. Evidence for signatures of recent diversifying selection on our identified noncoding loci in Darwin finch populations further suggests that regulatory rather than coding changes are major drivers of morphological diversification over macroevolutionary times.
Collapse
Affiliation(s)
- Leeban Yusuf
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, United Kingdom.,Centre for Biological Diversity, School of Biology, University of St. Andrews, Fife, KY16 9TF, United Kingdom
| | - Matthew C Heatley
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, United Kingdom.,Division of Plant and Crop Sciences, School of Biosciences, University of Nottingham, Sutton Bonington LE12 5RD, United Kingdom
| | - Joseph P G Palmer
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, United Kingdom.,School of Biological Sciences, Royal Holloway University of London, Egham, Surrey, TW20 0EX, United Kingdom
| | - Henry J Barton
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, United Kingdom.,Organismal and Evolutionary Biology Research Programme, Viikinkaari 9 (PL 56), University of Helsinki, Helsinki, FI-00014, Finland
| | - Christopher R Cooney
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, United Kingdom
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, United Kingdom.,Department of Animal Behaviour, Bielefeld University, Bielefeld, DE-33501, Germany
| |
Collapse
|
10
|
Affiliation(s)
- Toni I Gossmann
- Department of Animal Behaviour, Bielefeld University, 33501 Bielefeld, Germany
| | - Markus Ralser
- Molecular Biology of Metabolism Laboratory, The Francis Crick Institute, London, NW1 1AT, UK; Department of Biochemistry, Charité, 10117, Berlin, Germany.
| |
Collapse
|
11
|
Kinsella CM, Ruiz-Ruano FJ, Dion-Côté AM, Charles AJ, Gossmann TI, Cabrero J, Kappei D, Hemmings N, Simons MJP, Camacho JPM, Forstmeier W, Suh A. Programmed DNA elimination of germline development genes in songbirds. Nat Commun 2019; 10:5468. [PMID: 31784533 PMCID: PMC6884545 DOI: 10.1038/s41467-019-13427-4] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2019] [Accepted: 11/08/2019] [Indexed: 02/08/2023] Open
Abstract
In some eukaryotes, germline and somatic genomes differ dramatically in their composition. Here we characterise a major germline–soma dissimilarity caused by a germline-restricted chromosome (GRC) in songbirds. We show that the zebra finch GRC contains >115 genes paralogous to single-copy genes on 18 autosomes and the Z chromosome, and is enriched in genes involved in female gonad development. Many genes are likely functional, evidenced by expression in testes and ovaries at the RNA and protein level. Using comparative genomics, we show that genes have been added to the GRC over millions of years of evolution, with embryonic development genes bicc1 and trim71 dating to the ancestor of songbirds and dozens of other genes added very recently. The somatic elimination of this evolutionarily dynamic chromosome in songbirds implies a unique mechanism to minimise genetic conflict between germline and soma, relevant to antagonistic pleiotropy, an evolutionary process underlying ageing and sexual traits. Songbirds have extensive germline–soma genome differences due to developmental elimination of a germline-specific chromosome (GRC). Here, the authors show that the GRC contains dozens of expressed developmental genes, some of which have been on the GRC since the ancestor of all songbirds.
Collapse
Affiliation(s)
- Cormac M Kinsella
- Department of Ecology and Genetics - Evolutionary Biology, Evolutionary Biology Centre (EBC), Science for Life Laboratory, Uppsala University, SE-752 36, Uppsala, Sweden.,Laboratory of Experimental Virology, Department of Medical Microbiology, Amsterdam UMC, University of Amsterdam, 1105 AZ, Amsterdam, The Netherlands
| | - Francisco J Ruiz-Ruano
- Department of Ecology and Genetics - Evolutionary Biology, Evolutionary Biology Centre (EBC), Science for Life Laboratory, Uppsala University, SE-752 36, Uppsala, Sweden. .,Department of Genetics, University of Granada, E-18071, Granada, Spain. .,Department of Organismal Biology - Systematic Biology, Evolutionary Biology Centre (EBC), Science for Life Laboratory, Uppsala University, SE-752 36, Uppsala, Sweden.
| | - Anne-Marie Dion-Côté
- Department of Ecology and Genetics - Evolutionary Biology, Evolutionary Biology Centre (EBC), Science for Life Laboratory, Uppsala University, SE-752 36, Uppsala, Sweden.,Department of Molecular Biology & Genetics, Cornell University, Ithaca, NY, 14853, USA.,Département de Biologie, Université de Moncton, Moncton, NB, E1A 3E9, Canada
| | - Alexander J Charles
- Department of Animal and Plant Sciences, University of Sheffield, S10 2TN, Sheffield, UK
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, S10 2TN, Sheffield, UK.,Department of Animal Behaviour, Bielefeld University, D-33501, Bielefeld, Germany
| | - Josefa Cabrero
- Department of Genetics, University of Granada, E-18071, Granada, Spain
| | - Dennis Kappei
- Cancer Science Institute of Singapore, National University of Singapore, 117599, Singapore, Singapore.,Department of Biochemistry, Yong Loo Lin School of Medicine, National University of Singapore, 117596, Singapore, Singapore
| | - Nicola Hemmings
- Department of Animal and Plant Sciences, University of Sheffield, S10 2TN, Sheffield, UK
| | - Mirre J P Simons
- Department of Animal and Plant Sciences, University of Sheffield, S10 2TN, Sheffield, UK
| | | | | | - Alexander Suh
- Department of Ecology and Genetics - Evolutionary Biology, Evolutionary Biology Centre (EBC), Science for Life Laboratory, Uppsala University, SE-752 36, Uppsala, Sweden. .,Department of Organismal Biology - Systematic Biology, Evolutionary Biology Centre (EBC), Science for Life Laboratory, Uppsala University, SE-752 36, Uppsala, Sweden.
| |
Collapse
|
12
|
Gossmann TI, Shanmugasundram A, Börno S, Duvaux L, Lemaire C, Kuhl H, Klages S, Roberts LD, Schade S, Gostner JM, Hildebrand F, Vowinckel J, Bichet C, Mülleder M, Calvani E, Zelezniak A, Griffin JL, Bork P, Allaine D, Cohas A, Welch JJ, Timmermann B, Ralser M. Ice-Age Climate Adaptations Trap the Alpine Marmot in a State of Low Genetic Diversity. Curr Biol 2019; 29:1712-1720.e7. [PMID: 31080084 PMCID: PMC6538971 DOI: 10.1016/j.cub.2019.04.020] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Revised: 02/16/2019] [Accepted: 04/09/2019] [Indexed: 12/30/2022]
Abstract
Some species responded successfully to prehistoric changes in climate [1, 2], while others failed to adapt and became extinct [3]. The factors that determine successful climate adaptation remain poorly understood. We constructed a reference genome and studied physiological adaptations in the Alpine marmot (Marmota marmota), a large ground-dwelling squirrel exquisitely adapted to the "ice-age" climate of the Pleistocene steppe [4, 5]. Since the disappearance of this habitat, the rodent persists in large numbers in the high-altitude Alpine meadow [6, 7]. Genome and metabolome showed evidence of adaptation consistent with cold climate, affecting white adipose tissue. Conversely, however, we found that the Alpine marmot has levels of genetic variation that are among the lowest for mammals, such that deleterious mutations are less effectively purged. Our data rule out typical explanations for low diversity, such as high levels of consanguineous mating, or a very recent bottleneck. Instead, ancient demographic reconstruction revealed that genetic diversity was lost during the climate shifts of the Pleistocene and has not recovered, despite the current high population size. We attribute this slow recovery to the marmot's adaptive life history. The case of the Alpine marmot reveals a complicated relationship between climatic changes, genetic diversity, and conservation status. It shows that species of extremely low genetic diversity can be very successful and persist over thousands of years, but also that climate-adapted life history can trap a species in a persistent state of low genetic diversity.
Collapse
Affiliation(s)
- Toni I Gossmann
- University of Sheffield, Department of Animal and Plant Sciences, Sheffield S10 2TN, UK; Bielefeld University, Department of Animal Behaviour, 33501 Bielefeld, Germany
| | - Achchuthan Shanmugasundram
- Molecular Biology of Metabolism Laboratory, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK; Centre for Genomic Research, Institute of Integrative Biology, University of Liverpool, Biosciences Building, Crown Street, Liverpool L69 7ZB, UK
| | - Stefan Börno
- Max Planck Institute for Molecular Genetics, Sequencing Core Facility, Ihnestrasse 73, 14195 Berlin, Germany
| | - Ludovic Duvaux
- IRHS, Université d'Angers, INRA, Agrocampus-Ouest, SFR 4207 QuaSaV, 49071 Beaucouzé, France; BIOGECO, INRA, Université de Bordeaux, 69 Route d'Arcachon, 33612 Cestas, France
| | - Christophe Lemaire
- IRHS, Université d'Angers, INRA, Agrocampus-Ouest, SFR 4207 QuaSaV, 49071 Beaucouzé, France
| | - Heiner Kuhl
- Max Planck Institute for Molecular Genetics, Sequencing Core Facility, Ihnestrasse 73, 14195 Berlin, Germany; Department of Ecophysiology and Aquaculture, Leibniz-Institute of Freshwater Ecology and Inland Fisheries, 12587 Berlin, Germany
| | - Sven Klages
- Max Planck Institute for Molecular Genetics, Sequencing Core Facility, Ihnestrasse 73, 14195 Berlin, Germany
| | - Lee D Roberts
- Department of Biochemistry and Cambridge Systems Biology Centre, University of Cambridge, 80 Tennis Court Road, Cambridge CB2 1GA, UK; Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds LS2 9JT, UK
| | - Sophia Schade
- Max Planck Institute for Molecular Genetics, Sequencing Core Facility, Ihnestrasse 73, 14195 Berlin, Germany
| | - Johanna M Gostner
- Division of Medical Biochemistry, Medical University of Innsbruck, 6020 Innsbruck, Austria
| | - Falk Hildebrand
- European Molecular Biology Laboratory (EMBL), 69117 Heidelberg, Germany; Earlham Institute, Norwich Research Park, Norwich NR4 7UZ, UK; Gut Health and Microbes Programme, Quadram Institute, Norwich Research Park, Norwich NR4 7UQ, UK
| | - Jakob Vowinckel
- Department of Biochemistry and Cambridge Systems Biology Centre, University of Cambridge, 80 Tennis Court Road, Cambridge CB2 1GA, UK
| | | | - Michael Mülleder
- Department of Biochemistry and Cambridge Systems Biology Centre, University of Cambridge, 80 Tennis Court Road, Cambridge CB2 1GA, UK; Department of Biochemistry, Charitè, Am Chariteplatz 1, 10117 Berlin, Germany
| | - Enrica Calvani
- Molecular Biology of Metabolism Laboratory, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK; Department of Biochemistry and Cambridge Systems Biology Centre, University of Cambridge, 80 Tennis Court Road, Cambridge CB2 1GA, UK
| | - Aleksej Zelezniak
- Molecular Biology of Metabolism Laboratory, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK; Department of Biology and Biological Engineering, Chalmers University of Technology, 412 96 Göteborg, Sweden; Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm 171 65, Sweden
| | - Julian L Griffin
- Department of Biochemistry and Cambridge Systems Biology Centre, University of Cambridge, 80 Tennis Court Road, Cambridge CB2 1GA, UK
| | - Peer Bork
- European Molecular Biology Laboratory (EMBL), 69117 Heidelberg, Germany; Max-Delbrück-Centre for Molecular Medicine, 13092 Berlin, Germany; Molecular Medicine Partnership Unit, 69120 Heidelberg, Germany
| | - Dominique Allaine
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR 5558, Laboratoire de Biométrie et Biologie Evolutive, 69622 Villeurbanne, France
| | - Aurélie Cohas
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR 5558, Laboratoire de Biométrie et Biologie Evolutive, 69622 Villeurbanne, France
| | - John J Welch
- Department of Genetics, University of Cambridge, Cambridge CB2 3EH, UK
| | - Bernd Timmermann
- Max Planck Institute for Molecular Genetics, Sequencing Core Facility, Ihnestrasse 73, 14195 Berlin, Germany
| | - Markus Ralser
- Molecular Biology of Metabolism Laboratory, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK; Department of Biochemistry and Cambridge Systems Biology Centre, University of Cambridge, 80 Tennis Court Road, Cambridge CB2 1GA, UK; Department of Biochemistry, Charitè, Am Chariteplatz 1, 10117 Berlin, Germany.
| |
Collapse
|
13
|
Karr TL, Southern H, Rosenow MA, Gossmann TI, Snook RR. The Old and the New: Discovery Proteomics Identifies Putative Novel Seminal Fluid Proteins in Drosophila. Mol Cell Proteomics 2019; 18:S23-S33. [PMID: 30760537 PMCID: PMC6427231 DOI: 10.1074/mcp.ra118.001098] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2018] [Revised: 02/11/2019] [Indexed: 12/11/2022] Open
Abstract
Seminal fluid proteins (SFPs), the nonsperm component of male ejaculates produced by male accessory glands, are viewed as central mediators of reproductive fitness. SFPs effect both male and female post-mating functions and show molecular signatures of rapid adaptive evolution. Although Drosophila melanogaster, is the dominant insect model for understanding SFP evolution, understanding of SFP evolutionary causes and consequences require additional comparative analyses of close and distantly related taxa. Although SFP identification was historically challenging, advances in label-free quantitative proteomics expands the scope of studying other systems to further advance the field. Focused studies of SFPs has so far overlooked the proteomes of male reproductive glands and their inherent complex protein networks for which there is little information on the overall signals of molecular evolution. Here we applied label-free quantitative proteomics to identify the accessory gland proteome and secretome in Drosophila pseudoobscura,, a close relative of D. melanogaster,, and use the dataset to identify both known and putative novel SFPs. Using this approach, we identified 163 putative SFPs, 32% of which overlapped with previously identified D. melanogaster, SFPs and show that SFPs with known extracellular annotation evolve more rapidly than other proteins produced by or contained within the accessory gland. Our results will further the understanding of the evolution of SFPs and the underlying male accessory gland proteins that mediate reproductive fitness of the sexes.
Collapse
Affiliation(s)
- Timothy L Karr
- From the ‡Center for Mechanisms of Evolution, The Biodesign Institute, Arizona State University, Tempe, Arizona;.
| | - Helen Southern
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | | | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | - Rhonda R Snook
- Department of Zoology, Stockholm University, Stockholm, Sweden.
| |
Collapse
|
14
|
Laine VN, Gossmann TI, van Oers K, Visser ME, Groenen MAM. Exploring the unmapped DNA and RNA reads in a songbird genome. BMC Genomics 2019; 20:19. [PMID: 30621573 PMCID: PMC6323668 DOI: 10.1186/s12864-018-5378-2] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2018] [Accepted: 12/16/2018] [Indexed: 01/19/2023] Open
Abstract
BACKGROUND A widely used approach in next-generation sequencing projects is the alignment of reads to a reference genome. Despite methodological and hardware improvements which have enhanced the efficiency and accuracy of alignments, a significant percentage of reads frequently remain unmapped. Usually, unmapped reads are discarded from the analysis process, but significant biological information and insights can be uncovered from these data. We explored the unmapped DNA (normal and bisulfite treated) and RNA sequence reads of the great tit (Parus major) reference genome individual. From the unmapped reads we generated de novo assemblies, after which the generated sequence contigs were aligned to the NCBI non-redundant nucleotide database using BLAST, identifying the closest known matching sequence. RESULTS Many of the aligned contigs showed sequence similarity to different bird species and genes that were absent in the great tit reference assembly. Furthermore, there were also contigs that represented known P. major pathogenic species. Most interesting were several species of blood parasites such as Plasmodium and Trypanosoma. CONCLUSIONS Our analyses revealed that meaningful biological information can be found when further exploring unmapped reads. For instance, it is possible to discover sequences that are either absent or misassembled in the reference genome, and sequences that indicate infection or sample contamination. In this study we also propose strategies to aid the capture and interpretation of this information from unmapped reads.
Collapse
Affiliation(s)
- Veronika N Laine
- Department of Animal Ecology, NIOO-KNAW, Wageningen, The Netherlands.
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, The University of Sheffield, Sheffield, UK
| | - Kees van Oers
- Department of Animal Ecology, NIOO-KNAW, Wageningen, The Netherlands
| | - Marcel E Visser
- Department of Animal Ecology, NIOO-KNAW, Wageningen, The Netherlands.,Department of Animal Sciences, Wageningen University, Wageningen, The Netherlands
| | - Martien A M Groenen
- Department of Animal Sciences, Wageningen University, Wageningen, The Netherlands
| |
Collapse
|
15
|
Gossmann TI, Bockwoldt M, Diringer L, Schwarz F, Schumann VF. Evidence for Strong Fixation Bias at 4-fold Degenerate Sites Across Genes in the Great Tit Genome. Front Ecol Evol 2018. [DOI: 10.3389/fevo.2018.00203] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
16
|
Corcoran P, Gossmann TI, Barton HJ, Slate J, Zeng K. Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species. Genome Biol Evol 2018; 9:2987-3007. [PMID: 29045655 PMCID: PMC5714183 DOI: 10.1093/gbe/evx213] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/16/2017] [Indexed: 02/06/2023] Open
Abstract
Population genetic theory predicts that selection should be more effective when the effective population size (Ne) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger Ne. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.
Collapse
Affiliation(s)
- Pádraic Corcoran
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | - Henry J Barton
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | | | - Jon Slate
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | - Kai Zeng
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| |
Collapse
|
17
|
Afanasyeva A, Bockwoldt M, Cooney CR, Heiland I, Gossmann TI. Human long intrinsically disordered protein regions are frequent targets of positive selection. Genome Res 2018; 28:975-982. [PMID: 29858274 PMCID: PMC6028134 DOI: 10.1101/gr.232645.117] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Accepted: 06/01/2018] [Indexed: 12/20/2022]
Abstract
Intrinsically disordered regions occur frequently in proteins and are characterized by a lack of a well-defined three-dimensional structure. Although these regions do not show a higher order of structural organization, they are known to be functionally important. Disordered regions are rapidly evolving, largely attributed to relaxed purifying selection and an increased role of genetic drift. It has also been suggested that positive selection might contribute to their rapid diversification. However, for our own species, it is currently unknown whether positive selection has played a role during the evolution of these protein regions. Here, we address this question by investigating the evolutionary pattern of more than 6600 human proteins with intrinsically disordered regions and their ordered counterparts. Our comparative approach with data from more than 90 mammalian genomes uses a priori knowledge of disordered protein regions, and we show that this increases the power to detect positive selection by an order of magnitude. We can confirm that human intrinsically disordered regions evolve more rapidly, not only within humans but also across the entire mammalian phylogeny. They have, however, experienced substantial evolutionary constraint, hinting at their fundamental functional importance. We find compelling evidence that disordered protein regions are frequent targets of positive selection and estimate that the relative rate of adaptive substitutions differs fourfold between disordered and ordered protein regions in humans. Our results suggest that disordered protein regions are important targets of genetic innovation and that the contribution of positive selection in these regions is more pronounced than in other protein parts.
Collapse
Affiliation(s)
- Arina Afanasyeva
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S102TN, United Kingdom.,Institute of Nanobiotechnologies, Peter the Great St. Petersburg Polytechnic University, Saint-Petersburg 195251, Russia.,Petersburg Nuclear Physics Institute, B.P. Konstantinov NRC Kurchatov Institute, Gatchina, Leningrad District 188300, Russia.,National Institutes of Biomedical Innovation, Health and Nutrition, Ibaraki City, Osaka 567-0085, Japan
| | - Mathias Bockwoldt
- Department of Arctic and Marine Biology, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Christopher R Cooney
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S102TN, United Kingdom
| | - Ines Heiland
- Department of Arctic and Marine Biology, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S102TN, United Kingdom
| |
Collapse
|
18
|
Shafiey H, Gossmann TI, Waxman D. Evolutionary control: Targeted change of allele frequencies in natural populations using externally directed evolution. J Theor Biol 2017; 419:362-374. [PMID: 28130097 DOI: 10.1016/j.jtbi.2017.01.023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2016] [Revised: 11/13/2016] [Accepted: 01/16/2017] [Indexed: 11/30/2022]
Abstract
Random processes in biology, in particular random genetic drift, often make it difficult to predict the fate of a particular mutation in a population. Using principles of theoretical population genetics, we present a form of biological control that ensures a focal allele's frequency, at a given locus, achieves a prescribed probability distribution at a given time. This control is in the form of an additional evolutionary force that acts on a population. We provide the mathematical framework that determines the additional force. Our analysis indicates that generally the additional force depends on the frequency of the focal allele, and it may also depend on the time. We argue that translating this additional force into an externally controlled process, which has the possibility of being implemented in a number of different ways corresponding to selection, migration, mutation, or a combination of these, may provide a flexible instrument for targeted change of traits of interest in natural populations. This framework may be applied, or used as an informed form of guidance, in a variety of different biological scenarios including: yield and pesticide optimisation in crop production, biofermentation, the local regulation of human-associated natural populations, such as parasitic animals, or bacterial communities in hospitals.
Collapse
Affiliation(s)
- Hassan Shafiey
- Centre for Computational Systems Biology, Fudan University, 220 Handan Road, Shanghai 20433, People's Republic of China
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S102TN, UK
| | - David Waxman
- Centre for Computational Systems Biology, Fudan University, 220 Handan Road, Shanghai 20433, People's Republic of China.
| |
Collapse
|
19
|
Gossmann TI, Saleh D, Schmid MW, Spence MA, Schmid KJ. Transcriptomes of Plant Gametophytes Have a Higher Proportion of Rapidly Evolving and Young Genes than Sporophytes. Mol Biol Evol 2016; 33:1669-78. [PMID: 26956888 PMCID: PMC4915351 DOI: 10.1093/molbev/msw044] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Reproductive traits in plants tend to evolve rapidly due to various causes that include plant-pollinator coevolution and pollen competition, but the genomic basis of reproductive trait evolution is still largely unknown. To characterize evolutionary patterns of genome wide gene expression in reproductive tissues in the gametophyte and to compare them to developmental stages of the sporophyte, we analyzed evolutionary conservation and genetic diversity of protein-coding genes using microarray-based transcriptome data from three plant species, Arabidopsis thaliana, rice (Oryza sativa), and soybean (Glycine max). In all three species a significant shift in gene expression occurs during gametogenesis in which genes of younger evolutionary age and higher genetic diversity contribute significantly more to the transcriptome than in other stages. We refer to this phenomenon as "evolutionary bulge" during plant reproductive development because it differentiates the gametophyte from the sporophyte. We show that multiple, not mutually exclusive, causes may explain the bulge pattern, most prominently reduced tissue complexity of the gametophyte, a varying extent of selection on reproductive traits during gametogenesis as well as differences between male and female tissues. This highlights the importance of plant reproduction for understanding evolutionary forces determining the relationship of genomic and phenotypic variation in plants.
Collapse
Affiliation(s)
- Toni I Gossmann
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, Stuttgart, Germany Department of Animal and Plant Sciences, University of Sheffield, Sheffield, United Kingdom
| | - Dounia Saleh
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, Stuttgart, Germany
| | - Marc W Schmid
- Institute for Plant Biology and Zurich-Basel Plant Science Center, University of Zurich, Zurich, Switzerland
| | - Michael A Spence
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, United Kingdom
| | - Karl J Schmid
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, Stuttgart, Germany
| |
Collapse
|
20
|
Laine VN, Gossmann TI, Schachtschneider KM, Garroway CJ, Madsen O, Verhoeven KJF, de Jager V, Megens HJ, Warren WC, Minx P, Crooijmans RPMA, Corcoran P, Sheldon BC, Slate J, Zeng K, van Oers K, Visser ME, Groenen MAM. Evolutionary signals of selection on cognition from the great tit genome and methylome. Nat Commun 2016; 7:10474. [PMID: 26805030 PMCID: PMC4737754 DOI: 10.1038/ncomms10474] [Citation(s) in RCA: 128] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2015] [Accepted: 12/12/2015] [Indexed: 12/30/2022] Open
Abstract
For over 50 years, the great tit (Parus major) has been a model species for research in evolutionary, ecological and behavioural research; in particular, learning and cognition have been intensively studied. Here, to provide further insight into the molecular mechanisms behind these important traits, we de novo assemble a great tit reference genome and whole-genome re-sequence another 29 individuals from across Europe. We show an overrepresentation of genes related to neuronal functions, learning and cognition in regions under positive selection, as well as increased CpG methylation in these regions. In addition, great tit neuronal non-CpG methylation patterns are very similar to those observed in mammals, suggesting a universal role in neuronal epigenetic regulation which can affect learning-, memory- and experience-induced plasticity. The high-quality great tit genome assembly will play an instrumental role in furthering the integration of ecological, evolutionary, behavioural and genomic approaches in this model species. The great tit (Parus major) is known for its complex social-cognitive behaviour. Here, the authors sequence genomes of the great tit and show genes related to learning and cognition in regions under positive selection, as well as neuronal non-CpG methylation patterns similar to those observed in mammals.
Collapse
Affiliation(s)
- Veronika N Laine
- Department of Animal Ecology, Netherlands Institute of Ecology (NIOO-KNAW), P.O. Box 50, 6700AB Wageningen, The Netherlands
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, UK
| | - Kyle M Schachtschneider
- Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700AH Wageningen, The Netherlands.,Department of Animal Sciences, University of Illinois, Urbana, Illinois 61801, USA
| | - Colin J Garroway
- Edward Grey Institute, Department of Zoology, University of Oxford, Oxford OX1 3PS, UK
| | - Ole Madsen
- Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700AH Wageningen, The Netherlands
| | - Koen J F Verhoeven
- Department of Terrestrial Ecology, Netherlands Institute of Ecology (NIOO-KNAW), P.O. Box 50, 6700AB Wageningen, The Netherlands
| | - Victor de Jager
- Department of Animal Ecology, Netherlands Institute of Ecology (NIOO-KNAW), P.O. Box 50, 6700AB Wageningen, The Netherlands
| | - Hendrik-Jan Megens
- Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700AH Wageningen, The Netherlands
| | - Wesley C Warren
- The Genome Institute, Washington University School of Medicine, St Louis, Missouri 63108, USA
| | - Patrick Minx
- The Genome Institute, Washington University School of Medicine, St Louis, Missouri 63108, USA
| | - Richard P M A Crooijmans
- Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700AH Wageningen, The Netherlands
| | - Pádraic Corcoran
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, UK
| | | | - Ben C Sheldon
- Edward Grey Institute, Department of Zoology, University of Oxford, Oxford OX1 3PS, UK
| | - Jon Slate
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, UK
| | - Kai Zeng
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield S10 2TN, UK
| | - Kees van Oers
- Department of Animal Ecology, Netherlands Institute of Ecology (NIOO-KNAW), P.O. Box 50, 6700AB Wageningen, The Netherlands
| | - Marcel E Visser
- Department of Animal Ecology, Netherlands Institute of Ecology (NIOO-KNAW), P.O. Box 50, 6700AB Wageningen, The Netherlands.,Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700AH Wageningen, The Netherlands
| | - Martien A M Groenen
- Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700AH Wageningen, The Netherlands
| |
Collapse
|
21
|
Gossmann TI, Santure AW, Sheldon BC, Slate J, Zeng K. Highly variable recombinational landscape modulates efficacy of natural selection in birds. Genome Biol Evol 2015; 6:2061-75. [PMID: 25062920 PMCID: PMC4231635 DOI: 10.1093/gbe/evu157] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Determining the rate of protein evolution and identifying the causes of its variation across the genome are powerful ways to understand forces that are important for genome evolution. By using a multitissue transcriptome data set from great tit (Parus major), we analyzed patterns of molecular evolution between two passerine birds, great tit and zebra finch (Taeniopygia guttata), using the chicken genome (Gallus gallus) as an outgroup. We investigated whether a special feature of avian genomes, the highly variable recombinational landscape, modulates the efficacy of natural selection through the effects of Hill-Robertson interference, which predicts that selection should be more effective in removing deleterious mutations and incorporating beneficial mutations in high-recombination regions than in low-recombination regions. In agreement with these predictions, genes located in low-recombination regions tend to have a high proportion of neutrally evolving sites and relaxed selective constraint on sites subject to purifying selection, whereas genes that show strong support for past episodes of positive selection appear disproportionally in high-recombination regions. There is also evidence that genes located in high-recombination regions tend to have higher gene expression specificity than those located in low-recombination regions. Furthermore, more compact genes (i.e., those with fewer/shorter introns or shorter proteins) evolve faster than less compact ones. In sum, our results demonstrate that transcriptome sequencing is a powerful method to answer fundamental questions about genome evolution in nonmodel organisms.
Collapse
Affiliation(s)
- Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, United Kingdom
| | - Anna W Santure
- Department of Animal and Plant Sciences, University of Sheffield, United KingdomSchool of Biological Sciences, University of Auckland, New Zealand
| | - Ben C Sheldon
- Edward Grey Institute, Department of Zoology, University of Oxford, United Kingdom
| | - Jon Slate
- Department of Animal and Plant Sciences, University of Sheffield, United Kingdom
| | - Kai Zeng
- Department of Animal and Plant Sciences, University of Sheffield, United Kingdom
| |
Collapse
|
22
|
Gossmann TI, Ziegler M. Sequence divergence and diversity suggests ongoing functional diversification of vertebrate NAD metabolism. DNA Repair (Amst) 2014; 23:39-48. [PMID: 25084685 PMCID: PMC4248024 DOI: 10.1016/j.dnarep.2014.07.005] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2014] [Revised: 05/22/2014] [Accepted: 07/09/2014] [Indexed: 12/04/2022]
Abstract
NAD is not only an important cofactor in redox reactions but has also received attention in recent years because of its physiological importance in metabolic regulation, DNA repair and signaling. In contrast to the redox reactions, these regulatory processes involve degradation of NAD and therefore necessitate a constant replenishment of its cellular pool. NAD biosynthetic enzymes are common to almost all species in all clades, but the number of NAD degrading enzymes varies substantially across taxa. In particular, vertebrates, including humans, have a manifold of NAD degrading enzymes which require a high turnover of NAD. As there is currently a lack of a systematic study of how natural selection has shaped enzymes involved in NAD metabolism we conducted a comprehensive evolutionary analysis based on intraspecific variation and interspecific divergence. We compare NAD biosynthetic and degrading enzymes in four eukaryotic model species and subsequently focus on human NAD metabolic enzymes and their orthologs in other vertebrates. We find that the majority of enzymes involved in NAD metabolism are subject to varying levels of purifying selection. While NAD biosynthetic enzymes appear to experience a rather high level of evolutionary constraint, there is evidence for positive selection among enzymes mediating NAD-dependent signaling. This is particularly evident for members of the PARP family, a diverse protein family involved in DNA damage repair and programmed cell death. Based on haplotype information and substitution rate analysis we pinpoint sites that are potential targets of positive selection. We also link our findings to a three dimensional structure, which suggests that positive selection occurs in domains responsible for DNA binding and polymerization rather than the NAD catalytic domain. Taken together, our results indicate that vertebrate NAD metabolism is still undergoing functional diversification.
Collapse
Affiliation(s)
- Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, Alfred Denny Building, S10 2TN Sheffield, United Kingdom.
| | - Mathias Ziegler
- Department of Molecular Biology, University of Bergen, Postbox 7803, 5020 Bergen, Norway
| |
Collapse
|
23
|
Abstract
It is likely that the strength of selection acting upon a mutation varies through time due to changes in the environment. However, most population genetic theory assumes that the strength of selection remains constant. Here we investigate the consequences of fluctuating selection pressures on the quantification of adaptive evolution using McDonald-Kreitman (MK) style approaches. In agreement with previous work, we show that fluctuating selection can generate evidence of adaptive evolution even when the expected strength of selection on a mutation is zero. However, we also find that the mutations, which contribute to both polymorphism and divergence tend, on average, to be positively selected during their lifetime, under fluctuating selection models. This is because mutations that fluctuate, by chance, to positive selected values, tend to reach higher frequencies in the population than those that fluctuate towards negative values. Hence the evidence of positive adaptive evolution detected under a fluctuating selection model by MK type approaches is genuine since fixed mutations tend to be advantageous on average during their lifetime. Never-the-less we show that methods tend to underestimate the rate of adaptive evolution when selection fluctuates.
Collapse
Affiliation(s)
- Toni I. Gossmann
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, United Kingdom
| | - David Waxman
- Centre for Computational Systems Biology, Fudan University, Shanghai, China
| | - Adam Eyre-Walker
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
| |
Collapse
|
24
|
Gossmann TI, Schmid MW, Grossniklaus U, Schmid KJ. Selection-driven evolution of sex-biased genes is consistent with sexual selection in Arabidopsis thaliana. Mol Biol Evol 2013; 31:574-83. [PMID: 24273323 DOI: 10.1093/molbev/mst226] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Sex-biased genes are genes with a preferential or specific expression in one sex and tend to show an accelerated rate of evolution in animals. Various hypotheses--which are not mutually exclusive--have been put forth to explain observed patterns of rapid evolution. One possible explanation is positive selection, but this has been shown only in few animal species and mostly for male-specific genes. Here, we present a large-scale study that investigates evolutionary patterns of sex-biased genes in the predominantly self-fertilizing plant Arabidopsis thaliana. Unlike most animal species, A. thaliana does not possess sex chromosomes, its flowers develop both male and female sexual organs, and it is characterized by low outcrossing rates. Using cell-specific gene expression data, we identified genes whose expression is enriched in comparison with all other tissues in the male and female gametes (sperm, egg, and central cell), as well as in synergids, pollen, and pollen tubes, which also play an important role in reproduction. Genes specifically expressed in gametes and synergids show higher rates of protein evolution compared with the genome-wide average and no evidence for positive selection. In contrast, pollen- and pollen tube-specific genes not only have lower rates of protein evolution but also exhibit a higher proportion of adaptive amino acid substitutions. We show that this is the result of increased levels of purifying and positive selection among genes with pollen- and pollen tube-specific expression. The increased proportion of adaptive substitutions cannot be explained by the fact that pollen- and pollen tube-expressed genes are enriched in segmental duplications, are on average older, or have a larger effective population size. Our observations are consistent with prezygotic sexual selection as a result of interactions during pollination and pollen tube growth such as pollen tube competition.
Collapse
Affiliation(s)
- Toni I Gossmann
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, Stuttgart, Germany
| | | | | | | |
Collapse
|
25
|
Abstract
NAD is best known as an electron carrier and a cosubstrate of various redox reactions. However, over the past 20 years, NAD(+) has been shown to be a key signaling molecule that mediates post-translational protein modifications and serves as precursor of ADP-ribose-containing messenger molecules, which are involved in calcium mobilization. In contrast to its role as a redox carrier, NAD(+)-dependent signaling processes involve the release of nicotinamide (Nam) and require constant replenishment of cellular NAD(+) pools. So far, very little is known about the evolution of NAD(P) synthesis in eukaryotes. In the present study, genes involved in NAD(P) metabolism in 45 species were identified and analyzed with regard to similarities and differences in NAD(P) synthesis. The results show that the Preiss-Handler pathway and NAD(+) kinase are present in all organisms investigated, and thus seem to be ancestral routes. Additionally, two pathways exist that convert Nam to NAD(+); we identified several species that have apparently functional copies of both biosynthetic routes, which have been thought to be mutually exclusive. Furthermore, our findings suggest the parallel phylogenetic appearance of Nam N-methyltransferase, Nam phosphoribosyl transferase, and poly-ADP-ribosyltransferases.
Collapse
Affiliation(s)
- Toni I Gossmann
- Department of Molecular Biology, University of Bergen, Norway
| | | | | | | | | | | |
Collapse
|
26
|
Gossmann TI, Keightley PD, Eyre-Walker A. The effect of variation in the effective population size on the rate of adaptive molecular evolution in eukaryotes. Genome Biol Evol 2012; 4:658-67. [PMID: 22436998 PMCID: PMC3381672 DOI: 10.1093/gbe/evs027] [Citation(s) in RCA: 107] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The role of adaptation is a fundamental question in molecular evolution. Theory predicts that species with large effective population sizes should undergo a higher rate of adaptive evolution than species with low effective population sizes if adaptation is limited by the supply of mutations. Previous analyses have appeared to support this conjecture because estimates of the proportion of nonsynonymous substitutions fixed by adaptive evolution, α, tend to be higher in species with large Ne. However, α is a function of both the number of advantageous and effectively neutral substitutions, either of which might depend on Ne. Here, we investigate the relationship between Ne and ωa, the rate of adaptive evolution relative to the rate of neutral evolution, using nucleotide polymorphism and divergence data from 13 independent pairs of eukaryotic species. We find a highly significant positive correlation between ωa and Ne. We also find some evidence that the rate of adaptive evolution varies between groups of organisms for a given Ne. The correlation between ωa and Ne does not appear to be an artifact of demographic change or selection on synonymous codon use. Our results suggest that adaptation is to some extent limited by the supply of mutations and that at least some adaptation depends on newly occurring mutations rather than on standing genetic variation. Finally, we show that the proportion of nearly neutral nonadaptive substitutions declines with increasing Ne. The low rate of adaptive evolution and the high proportion of effectively neutral substitution in species with small Ne are expected to combine to make it difficult to detect adaptive molecular evolution in species with small Ne.
Collapse
Affiliation(s)
- Toni I Gossmann
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
| | | | | |
Collapse
|
27
|
Abstract
Gene duplications are one of the most important mechanisms for the origin of evolutionary novelties. Even though various models of the fate of duplicated genes have been established, current knowledge about the role of divergent selection after gene duplication is rather limited. In this study, we analyzed sequence divergence in response to neo- and subfunctionalization of segmentally duplicated genes in the genome of Arabidopsis thaliana. We compared the genomes of A. thaliana and the poplar Populus trichocarpa to identify orthologous pairs of genes and their corresponding inparalogs. Maximum-likelihood analyses of the nonsynonymous and synonymous substitution rate ratio [Formula: see text] of pairs of A. thaliana inparalogs were used to detect differences in the evolutionary rates of protein coding sequences. We analyzed 1,924 A. thaliana paralogous pairs and our results indicate that around 6.9% show divergent ω values between the lineages for a fraction of sites. We observe an enrichment of regulatory sequences, a reduced level of co-expression and an increased number of substitutions that can be attributed to positive selection based on an McDonald-Kreitman type of analysis. Taken together, these results show that selection after duplication contributes substantially to gene novelties and hence functional divergence in plants.
Collapse
Affiliation(s)
- Toni I Gossmann
- Leibniz-Institute of Plant Genetics and Crop Plant Research, Gatersleben, Germany.
| | | |
Collapse
|
28
|
Lau C, Dölle C, Gossmann TI, Agledal L, Niere M, Ziegler M. Isoform-specific targeting and interaction domains in human nicotinamide mononucleotide adenylyltransferases. J Biol Chem 2010; 285:18868-76. [PMID: 20388704 DOI: 10.1074/jbc.m110.107631] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Several important signaling pathways require NAD as substrate, thereby leading to significant consumption of the molecule. Because NAD is also an essential redox carrier, its continuous resynthesis is vital. In higher eukaryotes, maintenance of compartmentalized NAD pools is critical, but so far rather little is known about the regulation and subcellular distribution of NAD biosynthetic enzymes. The key step in NAD biosynthesis is the formation of the dinucleotide by nicotinamide/nicotinic acid mononucleotide adenylyltransferases (NMNATs). The three human isoforms were localized to the nucleus, the Golgi complex, and mitochondria. Here, we show that their genes contain unique exons that encode isoform-specific domains to mediate subcellular targeting and post-translational modifications. These domains are dispensable for catalytic activity, consistent with their absence from NMNATs of lower organisms. We further demonstrate that the Golgi-associated NMNAT is palmitoylated at two adjacent cysteine residues of its isoform-specific domain and thereby anchored at the cytoplasmic surface, a potential mechanism to regulate the cytosolic NAD pool. Insertion of unique domains thus provides a yet unrecognized enzyme targeting mode, which has also been adapted to modulate subcellular NAD supply.
Collapse
Affiliation(s)
- Corinna Lau
- Department of Molecular Biology, University of Bergen, N-5008 Bergen, Norway
| | | | | | | | | | | |
Collapse
|
29
|
Gossmann TI, Song BH, Windsor AJ, Mitchell-Olds T, Dixon CJ, Kapralov MV, Filatov DA, Eyre-Walker A. Genome wide analyses reveal little evidence for adaptive evolution in many plant species. Mol Biol Evol 2010; 27:1822-32. [PMID: 20299543 DOI: 10.1093/molbev/msq079] [Citation(s) in RCA: 203] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The relative contribution of advantageous and neutral mutations to the evolutionary process is a central problem in evolutionary biology. Current estimates suggest that whereas Drosophila, mice, and bacteria have undergone extensive adaptive evolution, hominids show little or no evidence of adaptive evolution in protein-coding sequences. This may be a consequence of differences in effective population size. To study the matter further, we have investigated whether plants show evidence of adaptive evolution using an extension of the McDonald-Kreitman test that explicitly models slightly deleterious mutations by estimating the distribution of fitness effects of new mutations. We apply this method to data from nine pairs of species. Altogether more than 2,400 loci with an average length of approximately 280 nucleotides were analyzed. We observe very similar results in all species; we find little evidence of adaptive amino acid substitution in any comparison except sunflowers. This may be because many plant species have modest effective population sizes.
Collapse
Affiliation(s)
- Toni I Gossmann
- Centre for the Study of Evolution, School of Life Sciences, University of Sussex, Brighton, United Kingdom
| | | | | | | | | | | | | | | |
Collapse
|