1
|
Gutiérrez-Valencia J, Zervakis PI, Postel Z, Fracassetti M, Losvik A, Mehrabi S, Bunikis I, Soler L, Hughes PW, Désamoré A, Laenen B, Abdelaziz M, Pettersson OV, Arroyo J, Slotte T. Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum. Mol Biol Evol 2024; 41:msae087. [PMID: 38709782 PMCID: PMC11114476 DOI: 10.1093/molbev/msae087] [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: 11/17/2023] [Revised: 03/22/2024] [Accepted: 04/29/2024] [Indexed: 05/08/2024] Open
Abstract
Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.
Collapse
Affiliation(s)
- Juanita Gutiérrez-Valencia
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Panagiotis-Ioannis Zervakis
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Zoé Postel
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Marco Fracassetti
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Aleksandra Losvik
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Sara Mehrabi
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Ignas Bunikis
- Department of Immunology, Genetics and Pathology, Uppsala Genome Center, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Lucile Soler
- Department of Medical Biochemistry and Microbiology, Uppsala University, National Bioinformatics Infrastructure Sweden (NBIS), Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - P William Hughes
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Aurélie Désamoré
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Benjamin Laenen
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | | | - Olga Vinnere Pettersson
- Department of Immunology, Genetics and Pathology, Uppsala Genome Center, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Juan Arroyo
- Department of Plant Biology and Ecology, University of Seville, Seville, Spain
| | - Tanja Slotte
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| |
Collapse
|
2
|
Martínez Villar M, Bechsgaard J, Bilde T, Albo MJ, Tomasco IH. Impact of pre-copulatory sexual cannibalism on genetic diversity and efficacy of selection. Biol Lett 2024; 20:20230505. [PMID: 38746981 DOI: 10.1098/rsbl.2023.0505] [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/30/2023] [Accepted: 03/21/2024] [Indexed: 07/31/2024] Open
Abstract
Factors that increase reproductive variance among individuals act to reduce effective population size (Ne), which accelerates the loss of genetic diversity and decreases the efficacy of purifying selection. These factors include sexual cannibalism, offspring investment and mating system. Pre-copulatory sexual cannibalism, where the female consumes the male prior to mating, exacerbates this effect. We performed comparative transcriptomics in two spider species, the cannibalistic Trechaleoides biocellata and the non-cannibalistic T. keyserlingi, to generate genomic evidence to support these predictions. First, we estimated heterozygosity and found that genetic diversity is relatively lower in the cannibalistic species. Second, we calculated dN/dS ratios as a measure of purifying selection; a higher dN/dS ratio indicated relaxed purifying selection in the cannibalistic species. These results are consistent with the hypothesis that sexual cannibalism impacts operational sex ratio and demographic processes, which interact with evolutionary forces to shape the genetic structure of populations. However, other factors such as the mating system and life-history traits contribute to shaping Ne. Comparative analyses across multiple contrasting species pairs would be required to disentangle these effects. Our study highlights that extreme behaviours such as pre-copulatory cannibalism may have profound eco-evolutionary effects.
Collapse
Affiliation(s)
- Mauro Martínez Villar
- Departamento de Ecología y Evolución, Facultad de Ciencias, Universidad de la República , Montevideo, 11400, Uruguay
| | | | - Trine Bilde
- Departament of Biology, Aarhus University , Aarhus 8000, Denmark
| | - Maria Jose Albo
- Departament of Biology, Aarhus University , Aarhus 8000, Denmark
- Laboratory of Ethology, Ecology and Evolution, Clemente Estable Biological Research Institute , Montevideo, 11600, Uruguay
| | - Ivanna H Tomasco
- Departamento de Ecología y Evolución, Facultad de Ciencias, Universidad de la República , Montevideo, 11400, Uruguay
| |
Collapse
|
3
|
Zhao W, Gao J, Hall D, Andersson BA, Bruxaux J, Tomlinson KW, Drouzas AD, Suyama Y, Wang XR. Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow. THE NEW PHYTOLOGIST 2024. [PMID: 38515228 DOI: 10.1111/nph.19694] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Accepted: 03/04/2024] [Indexed: 03/23/2024]
Abstract
Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers withc . $$ c. $$ 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.
Collapse
Affiliation(s)
- Wei Zhao
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
| | - Jie Gao
- CAS Key Laboratory of Tropical Forest Ecology, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Yunnan, 666303, China
| | - David Hall
- Forestry Research Institute of Sweden (Skogforsk), Sävar, SE-91833, Sweden
| | - Bea Angelica Andersson
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
| | - Jade Bruxaux
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
| | - Kyle W Tomlinson
- Center for Integrative Conservation & Yunnan Key Laboratory for Conservation of Tropical Rainforests and Asian Elephant, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Yunnan, 666303, China
| | - Andreas D Drouzas
- Laboratory of Systematic Botany and Phytogeography, School of Biology, Aristotle University of Thessaloniki, Thessaloniki, 54124, Greece
| | - Yoshihisa Suyama
- Graduate School of Agricultural Science, Tohoku University, Miyagi, 989-6711, Japan
| | - Xiao-Ru Wang
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
- National Engineering Research Center of Tree Breeding and Ecological Restoration, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, 100083, China
| |
Collapse
|
4
|
Russo CAM, Eyre-Walker A, Katz LA, Gaut BS. Forty Years of Inferential Methods in the Journals of the Society for Molecular Biology and Evolution. Mol Biol Evol 2024; 41:msad264. [PMID: 38197288 PMCID: PMC10763999 DOI: 10.1093/molbev/msad264] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Accepted: 11/27/2023] [Indexed: 01/11/2024] Open
Abstract
We are launching a series to celebrate the 40th anniversary of the first issue of Molecular Biology and Evolution. In 2024, we will publish virtual issues containing selected papers published in the Society for Molecular Biology and Evolution journals, Molecular Biology and Evolution and Genome Biology and Evolution. Each virtual issue will be accompanied by a perspective that highlights the historic and contemporary contributions of our journals to a specific topic in molecular evolution. This perspective, the first in the series, presents an account of the broad array of methods that have been published in the Society for Molecular Biology and Evolution journals, including methods to infer phylogenies, to test hypotheses in a phylogenetic framework, and to infer population genetic processes. We also mention many of the software implementations that make methods tractable for empiricists. In short, the Society for Molecular Biology and Evolution community has much to celebrate after four decades of publishing high-quality science including numerous important inferential methods.
Collapse
Affiliation(s)
- Claudia A M Russo
- Departamento de Genética, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | | | - Laura A Katz
- Department of Biological Sciences, Smith College, Northampton, MA, USA
| | - Brandon S Gaut
- School of Biological Sciences, University of California, Irvine, CA, USA
| |
Collapse
|
5
|
Wang X, Ingvarsson PK. Quantifying adaptive evolution and the effects of natural selection across the Norway spruce genome. Mol Ecol 2023; 32:5288-5304. [PMID: 37622583 DOI: 10.1111/mec.17106] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Revised: 08/07/2023] [Accepted: 08/09/2023] [Indexed: 08/26/2023]
Abstract
Detecting natural selection is one of the major goals of evolutionary genomics. Here, we sequenced the whole genome of 25 Picea abies individuals and quantified the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we showed that both negative selection and the rate of positively selected substitutions are very limited in coding regions. We found a positive correlation between the rate of adaptive substitutions and recombination rate and a negative correlation between the rate of adaptive substitutions and gene density, suggesting a widespread influence from Hill-Robertson interference on the efficiency of protein adaptation in P. abies. Finally, the distinct population statistics between genomic regions under either positive or balancing selection with that under neutral regions indicated the impact of natural selection on the genomic architecture of Norway spruce. Further gene ontology enrichment analysis for genes located in regions identified as undergoing either positive or long-term balancing selection also highlighted the specific molecular functions and biological processes that appear to be targets of selection in Norway spruce.
Collapse
Affiliation(s)
- Xi Wang
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden
| | - Pär K Ingvarsson
- Linnean Centre for Plant Biology, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
6
|
Andersson BA, Zhao W, Haller BC, Brännström Å, Wang XR. Inference of the distribution of fitness effects of mutations is affected by single nucleotide polymorphism filtering methods, sample size and population structure. Mol Ecol Resour 2023; 23:1589-1603. [PMID: 37340611 DOI: 10.1111/1755-0998.13825] [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: 06/20/2022] [Revised: 06/02/2023] [Accepted: 06/08/2023] [Indexed: 06/22/2023]
Abstract
The distribution of fitness effects (DFE) of new mutations has been of interest to evolutionary biologists since the concept of mutations arose. Modern population genomic data enable us to quantify the DFE empirically, but few studies have examined how data processing, sample size and cryptic population structure might affect the accuracy of DFE inference. We used simulated and empirical data (from Arabidopsis lyrata) to show the effects of missing data filtering, sample size, number of single nucleotide polymorphisms (SNPs) and population structure on the accuracy and variance of DFE estimates. Our analyses focus on three filtering methods-downsampling, imputation and subsampling-with sample sizes of 4-100 individuals. We show that (1) the choice of missing-data treatment directly affects the estimated DFE, with downsampling performing better than imputation and subsampling; (2) the estimated DFE is less reliable in small samples (<8 individuals), and becomes unpredictable with too few SNPs (<5000, the sum of 0- and 4-fold SNPs); and (3) population structure may skew the inferred DFE towards more strongly deleterious mutations. We suggest that future studies should consider downsampling for small data sets, and use samples larger than 4 (ideally larger than 8) individuals, with more than 5000 SNPs in order to improve the robustness of DFE inference and enable comparative analyses.
Collapse
Affiliation(s)
| | - Wei Zhao
- Department of Ecology and Environmental Sciences, Umeå University, Umeå, Sweden
| | - Benjamin C Haller
- Department of Computational Biology, Cornell University, Ithaca, New York, USA
| | - Åke Brännström
- Department of Mathematics and Mathematical Statistics, Umeå University, Umeå, Sweden
- Advancing Systems Analysis Program, International Institute for Applied Systems Analysis, Laxenburg, Austria
- Complexity Science and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, Kunigami, Japan
| | - Xiao-Ru Wang
- Department of Ecology and Environmental Sciences, Umeå University, Umeå, Sweden
| |
Collapse
|
7
|
Fyon F, Berbel‐Filho WM. Influence of the mutation load on the genomic composition of hybrids between outcrossing and self-fertilizing species. Ecol Evol 2023; 13:e10538. [PMID: 37720059 PMCID: PMC10502466 DOI: 10.1002/ece3.10538] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 07/24/2023] [Accepted: 09/04/2023] [Indexed: 09/19/2023] Open
Abstract
Hybridization is a natural process whereby two diverging evolutionary lineages reproduce and create offspring of mixed ancestry. Differences in mating systems (e.g., self-fertilization and outcrossing) are expected to affect the direction and extent of hybridization and introgression in hybrid zones. Among other factors, selfers and outcrossers are expected to differ in their mutation loads. This has been studied both theoretically and empirically; however, conflicting predictions have been made on the effects mutation loads of parental species with different mating systems can have on the genomic composition of hybrids. Here, we develop a multi-locus, selective model to study how the different mutation load built up in selfers and outcrossers as a result of selective interference and homozygosity impact the long-term genetic composition of hybrid populations. Notably, our results emphasize that genes from the parental population with lesser mutation load get rapidly overrepresented in hybrid genomes, regardless of the hybrids own mating system. When recombination tends to be more important than mutation, outcrossers' genomes tend to be of higher quality and prevail. When recombination rates are low, however, selfers' genomes may reach higher quality than outcrossers' genomes and prevail in the hybrids. Taken together, these results provide concrete insights into one of the multiple factors influencing hybrid genome ancestry and introgression patterns in hybrid zones containing species with different mating systems.
Collapse
Affiliation(s)
- Fréderic Fyon
- Department of BiologyRoyal Holloway University of LondonEghamUK
| | | |
Collapse
|
8
|
Whole Transcriptome Sequencing Unveils the Genomic Determinants of Putative Somaclonal Variation in Mint ( Mentha L.). Int J Mol Sci 2022; 23:ijms23105291. [PMID: 35628103 PMCID: PMC9141282 DOI: 10.3390/ijms23105291] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2022] [Revised: 05/06/2022] [Accepted: 05/07/2022] [Indexed: 12/21/2022] Open
Abstract
Mint (Mentha L., Lamiaceae) is a strongly scented herb of the family Lamiaceae that is grown mostly by clonal propagation, making it a valuable species for the study of somaclonal variation and its phenotypic consequences. The recent introduction of a few species of mint in South America, followed by a presumably rampant propagation, make this region particularly ideal for studying the extent of somaclonal genetic diversity. Hence, the objective of this work was to offer a preliminary characterization of somaclonal genetically coding diversity of the mint in the northern Andes in order to address the question of whether somaclonal variants may have emerged despite relatively recent introductions in a region where mint is not native. A total of 29 clonally propagated specimens, collected in mint export farms in the province of Antioquia, a major region for mint production in the northwest Andes of Colombia, were genotyped using RNA sequencing (RNA-Seq). SNP calling was carried out from the leaves’ transcriptome profiles of each plant by combining the GATK4 and TRINITY protocols, obtaining a total of 2033 loci across 912 transcripts with a minimum read depth of 20X and 4% of missing data. Unsupervised machine learning algorithms considered the K-means, AGNES and UPGMA approaches, all of which suggested three genetic clusters for M. spicata and a unique cluster for M. × piperita. The results indicate that at least two different origins of M. spicata reached the eastern region of the Antioquia province, clonally propagated in the locality ever since for local consumption and export. One of these ancestries had more population structure, possibly due to environmental or anthropological pressures that intervened in the fragmentation of this genetic group or to a higher somaclonal mutation rate. This work offers a first step into the study of the accumulation and transmission of presumably quasi-neutral somatic mutations at coding regions in an herbaceous clonally propagated scented species such as mint, likely favored by an expected population expansion after its Andean introduction. These ad hoc hypotheses warrant further study as part of future research.
Collapse
|
9
|
Liang YY, Chen XY, Zhou BF, Mitchell-Olds T, Wang B. Globally Relaxed Selection and Local Adaptation in Boechera stricta. Genome Biol Evol 2022; 14:evac043. [PMID: 35349686 PMCID: PMC9011030 DOI: 10.1093/gbe/evac043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/23/2022] [Indexed: 11/25/2022] Open
Abstract
The strength of selection varies among populations and across the genome, but the determinants of efficacy of selection remain unclear. In this study, we used whole-genome sequencing data from 467 Boechera stricta accessions to quantify the strength of selection and characterize the pattern of local adaptation. We found low genetic diversity on 0-fold degenerate sites and conserved non-coding sites, indicating functional constraints on these regions. The estimated distribution of fitness effects and the proportion of fixed substitutions suggest relaxed negative and positive selection in B. stricta. Among the four population groups, the NOR and WES groups have smaller effective population size (Ne), higher proportions of effectively neutral sites, and lower rates of adaptive evolution compared with UTA and COL groups, reflecting the effect of Ne on the efficacy of natural selection. We also found weaker selection on GC-biased sites compared with GC-conservative (unbiased) sites, suggested that GC-biased gene conversion has affected the strength of selection in B. stricta. We found mixed evidence for the role of the recombination rate on the efficacy of selection. The positive and negative selection was stronger in high-recombination regions compared with low-recombination regions in COL but not in other groups. By scanning the genome, we found different subsets of selected genes suggesting differential adaptation among B. stricta groups. These results show that differences in effective population size, nucleotide composition, and recombination rate are important determinants of the efficacy of selection. This study enriches our understanding of the roles of natural selection and local adaptation in shaping genomic variation.
Collapse
Affiliation(s)
- Yi-Ye Liang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences,
Guangzhou, China
- University of the Chinese Academy of Sciences, Beijing, China
| | - Xue-Yan Chen
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences,
Guangzhou, China
- University of the Chinese Academy of Sciences, Beijing, China
| | - Biao-Feng Zhou
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences,
Guangzhou, China
- University of the Chinese Academy of Sciences, Beijing, China
| | | | - Baosheng Wang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences,
Guangzhou, China
- Center of Conservation Biology, Core Botanical Gardens, Chinese Academy of Sciences, Guangzhou, China
| |
Collapse
|
10
|
Abstract
It is known that methods to estimate the rate of adaptive evolution, which are based on the McDonald–Kreitman test, can be biased by changes in effective population size. Here, we demonstrate theoretically that changes in population size can also generate an artifactual correlation between the rate of adaptive evolution and any factor that is correlated to the strength of selection acting against deleterious mutations. In this context, we have investigated whether several site-level factors influence the rate of adaptive evolution in the divergence of humans and chimpanzees, two species that have been inferred to have undergone population size contraction since they diverged. We find that the rate of adaptive evolution, relative to the rate of mutation, is higher for more exposed amino acids, lower for amino acid pairs that are more dissimilar in terms of their polarity, volume, and lower for amino acid pairs that are subject to stronger purifying selection, as measured by the ratio of the numbers of nonsynonymous to synonymous polymorphisms (pN/pS). All of these correlations are opposite to the artifactual correlations expected under contracting population size. We therefore conclude that these correlations are genuine.
Collapse
Affiliation(s)
- Vivak Soni
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
| | - Ana Filipa Moutinho
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
- Department for Evolutionary Genetics, Max Planck Institute for Evolutionary Biology, Plon, Germany
| | - Adam Eyre-Walker
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
- Corresponding author: E-mail:
| |
Collapse
|
11
|
Horvath R, Josephs EB, Pesquet E, Stinchcombe JR, Wright SI, Scofield D, Slotte T. Selection on Accessible Chromatin Regions in Capsella grandiflora. Mol Biol Evol 2021; 38:5563-5575. [PMID: 34498072 PMCID: PMC8662636 DOI: 10.1093/molbev/msab270] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.
Collapse
Affiliation(s)
- Robert Horvath
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Emily B Josephs
- Department of Plant Biology, Michigan State University, Lansing, MI, USA
| | - Edouard Pesquet
- Department of Ecology, Environment and Plant Sciences, Stockholm University, Stockholm, Sweden
| | - John R Stinchcombe
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - Stephen I Wright
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - Douglas Scofield
- Department of Ecology and Genetics, Uppsala University, Uppsala, Sweden
| | - Tanja Slotte
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| |
Collapse
|
12
|
Xu K. Mutation accumulation in inbreeding populations under evolution of the selfing rate. J Evol Biol 2021; 35:23-39. [PMID: 34860448 DOI: 10.1111/jeb.13968] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Revised: 11/17/2021] [Accepted: 11/23/2021] [Indexed: 11/28/2022]
Abstract
It is theoretically established that self-fertilization can facilitate mutation accumulation, thus increasing extinction risk. However, in previous studies, selfing rates are often set as fixed parameters, but in natural systems, evolution of selfing rates and deleterious mutations may mutually affect each other. I carried out simulations to investigate the dynamics of selfing rates and mutation accumulation, by allowing deleterious mutations to coevolve with alleles that modify the selfing rate (selfing modifiers). I found that selfing rates will often fluctuate over time, due to successive invasion of alleles that increase selfing and outcrossing. Since mutation fixation is mainly caused by Muller's ratchet, its rate is sensitive to the change of the selfing rate mutations will accumulate in a punctuated pattern. The dynamics are influenced by several factors, such as recombination and the selfing rate effects of selfing modifier loci. Also, such temporal variation produces variation of selfing rates and mutation accumulation rates between multiple conspecific populations, which can increase the average fitness across populations. As factors, such as the genomic mutation rate of deleterious mutations, can simultaneously influence the selfing rate and mutation fixation, effects of these factors on mutation accumulation rates can be complicated and non-monotonic.
Collapse
Affiliation(s)
- Kuangyi Xu
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina, USA
| |
Collapse
|
13
|
Wang XJ, Barrett SCH, Zhong L, Wu ZK, Li DZ, Wang H, Zhou W. The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly. Mol Biol Evol 2021; 38:168-180. [PMID: 32761213 PMCID: PMC7782863 DOI: 10.1093/molbev/msaa199] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (dN/dS) and polymorphism (πN/πS) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.
Collapse
Affiliation(s)
- Xin-Jia Wang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,Plant Germplasm and Genomics Center, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Spencer C H Barrett
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Canada
| | - Li Zhong
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Zhi-Kun Wu
- Department of Pharmacy, Guizhou University of Traditional Chinese Medicine, Guiyang, Guizhou, China
| | - De-Zhu Li
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,Plant Germplasm and Genomics Center, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - Hong Wang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - Wei Zhou
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,Plant Germplasm and Genomics Center, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| |
Collapse
|
14
|
Holding ML, Sovic MG, Colston TJ, Gibbs HL. The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey. Biol J Linn Soc Lond 2020. [DOI: 10.1093/biolinnean/blaa192] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Abstract
Coevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.
Collapse
Affiliation(s)
- Matthew L Holding
- Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
- Department of Biological Science, Florida State University, Tallahassee, FL, USA
| | - Michael G Sovic
- Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
| | - Timothy J Colston
- Department of Biological Science, Florida State University, Tallahassee, FL, USA
| | - H Lisle Gibbs
- Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
| |
Collapse
|
15
|
Kayum MA, Nath UK, Park JI, Hossain MR, Kim HT, Kim HR, Nou IS. Glucosinolate profile and Myrosinase gene expression are modulated upon Plasmodiophora brassicae infection in cabbage. FUNCTIONAL PLANT BIOLOGY : FPB 2020; 48:103-118. [PMID: 32780986 DOI: 10.1071/fp19278] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/03/2019] [Accepted: 07/20/2020] [Indexed: 06/11/2023]
Abstract
Clubroot is a devastating disease of Brassicaceae caused by the biotrophic protist Plasmodiophora brassicae. The progression of clubroot disease is modulated by the glucosinolate (GSL) profile of the host plant. GSL is hydrolysed by the enzyme myrosinase upon cell disruption and gives rise to metabolites like isothiocyanate, nitriles, thiocyanates, epithionitriles and oxazolidines. Some of these metabolites play important roles in the plant's defence mechanism. We identified 13 Myrosinase (Myro) and 28 Myrosinase-Binding Protein-like (MBP) genes from Brassica oleracea L. using a comparative genomics approach and characterised them through in silico analyses. We compared the expression patterns of these genes in a clubroot-susceptible line and a resistant line following inoculation with P. brassicae. Two BolMyro and 12 BolMBP genes were highly expressed in the susceptible line, whereas only one BolMyro and five BolMBP genes were highly expressed in the resistant line. Principal component analysis confirmed that specific GSL profiles and gene expression were modulated due to pathogen infection. Plants with higher levels of neoglucobrassicin, glucobrassicin and methooxyglucobrassicin produced disease symptoms and formed galls, whereas, plants with higher levels of sinigrin, hydroxyglucobrassicin and progoitrin produced less symptoms with almost no galls. Our results provide insights into the roles of Myro and MBP genes in GSL hydrolysis during P. brassicae infection, which will help for developing clubroot resistant cabbage lines.
Collapse
Affiliation(s)
- Md Abdul Kayum
- Department of Horticulture, Sunchon National University, 255 Jungang-ro, Suncheon, Jeonnam 57922, South Korea
| | - Ujjal Kumar Nath
- Department of Horticulture, Sunchon National University, 255 Jungang-ro, Suncheon, Jeonnam 57922, South Korea
| | - Jong-In Park
- Department of Horticulture, Sunchon National University, 255 Jungang-ro, Suncheon, Jeonnam 57922, South Korea; and Corresponding authors. ;
| | - Mohammad Rashed Hossain
- Department of Horticulture, Sunchon National University, 255 Jungang-ro, Suncheon, Jeonnam 57922, South Korea
| | - Hoy-Taek Kim
- Department of Horticulture, Sunchon National University, 255 Jungang-ro, Suncheon, Jeonnam 57922, South Korea
| | - Hye-Ran Kim
- Plant Systems Engineering Research Center, Korea Research Institute of Bioscience and Biotechnology, Daejeon, South Korea
| | - Ill-Sup Nou
- Department of Horticulture, Sunchon National University, 255 Jungang-ro, Suncheon, Jeonnam 57922, South Korea; and Corresponding authors. ;
| |
Collapse
|
16
|
Kutschera VE, Poelstra JW, Botero-Castro F, Dussex N, Gemmell NJ, Hunt GR, Ritchie MG, Rutz C, Wiberg RAW, Wolf JBW. Purifying Selection in Corvids Is Less Efficient on Islands. Mol Biol Evol 2020; 37:469-474. [PMID: 31633794 PMCID: PMC6993847 DOI: 10.1093/molbev/msz233] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size Ne. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species.
Collapse
Affiliation(s)
- Verena E Kutschera
- Department of Evolutionary Biology, Uppsala University, Uppsala, Sweden.,Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden
| | | | - Fidel Botero-Castro
- Division of Evolutionary Biology, Faculty of Biology, LMU Munich, Planegg-Martinsried, Germany
| | - Nicolas Dussex
- Department of Anatomy, University of Otago, Dunedin, New Zealand.,Department of Bioinformatics and Genetics, Swedish Museum of Natural History, Stockholm, Sweden
| | - Neil J Gemmell
- Department of Anatomy, University of Otago, Dunedin, New Zealand
| | | | - Michael G Ritchie
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, United Kingdom
| | - Christian Rutz
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, United Kingdom
| | - R Axel W Wiberg
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, United Kingdom.,Department of Environmental Sciences, Evolutionary Biology, University of Basel, Basel, Switzerland
| | - Jochen B W Wolf
- Department of Evolutionary Biology, Uppsala University, Uppsala, Sweden.,Division of Evolutionary Biology, Faculty of Biology, LMU Munich, Planegg-Martinsried, Germany
| |
Collapse
|
17
|
Hämälä T, Tiffin P. Biased Gene Conversion Constrains Adaptation in Arabidopsis thaliana. Genetics 2020; 215:831-846. [PMID: 32414868 PMCID: PMC7337087 DOI: 10.1534/genetics.120.303335] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2020] [Accepted: 05/14/2020] [Indexed: 02/01/2023] Open
Abstract
Reduction of fitness due to deleterious mutations imposes a limit to adaptive evolution. By characterizing features that influence this genetic load we may better understand constraints on responses to both natural and human-mediated selection. Here, using whole-genome, transcriptome, and methylome data from >600 Arabidopsis thaliana individuals, we set out to identify important features influencing selective constraint. Our analyses reveal that multiple factors underlie the accumulation of maladaptive mutations, including gene expression level, gene network connectivity, and gene-body methylation. We then focus on a feature with major effect, nucleotide composition. The ancestral vs. derived status of segregating alleles suggests that GC-biased gene conversion, a recombination-associated process that increases the frequency of G and C nucleotides regardless of their fitness effects, shapes sequence patterns in A. thaliana Through estimation of mutational effects, we present evidence that biased gene conversion hinders the purging of deleterious mutations and contributes to a genome-wide signal of decreased efficacy of selection. By comparing these results to two outcrossing relatives, Arabidopsis lyrata and Capsella grandiflora, we find that protein evolution in A. thaliana is as strongly affected by biased gene conversion as in the outcrossing species. Last, we perform simulations to show that natural levels of outcrossing in A. thaliana are sufficient to facilitate biased gene conversion despite increased homozygosity due to selfing. Together, our results show that even predominantly selfing taxa are susceptible to biased gene conversion, suggesting that it may constitute an important constraint to adaptation among plant species.
Collapse
Affiliation(s)
- Tuomas Hämälä
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, Minnesota 55108
| | - Peter Tiffin
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, Minnesota 55108
| |
Collapse
|
18
|
Josephs EB, Lee YW, Wood CW, Schoen DJ, Wright SI, Stinchcombe JR. The Evolutionary Forces Shaping Cis- and Trans-Regulation of Gene Expression within a Population of Outcrossing Plants. Mol Biol Evol 2020; 37:2386-2393. [DOI: 10.1093/molbev/msaa102] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Abstract
Understanding the persistence of genetic variation within populations has long been a goal of evolutionary biology. One promising route toward achieving this goal is using population genetic approaches to describe how selection acts on the loci associated with trait variation. Gene expression provides a model trait for addressing the challenge of the maintenance of variation because it can be measured genome-wide without information about how gene expression affects traits. Previous work has shown that loci affecting the expression of nearby genes (local or cis-eQTLs) are under negative selection, but we lack a clear understanding of the selective forces acting on variants that affect the expression of genes in trans. Here, we identify loci that affect gene expression in trans using genomic and transcriptomic data from one population of the obligately outcrossing plant, Capsella grandiflora. The allele frequencies of trans-eQTLs are consistent with stronger negative selection acting on trans-eQTLs than cis-eQTLs, and stronger negative selection acting on trans-eQTLs associated with the expression of multiple genes. However, despite this general pattern, we still observe the presence of a trans-eQTL at intermediate frequency that affects the expression of a large number of genes in the same coexpression module. Overall, our work highlights the different selective pressures shaping variation in cis- and trans-regulation.
Collapse
Affiliation(s)
- Emily B Josephs
- Department of Plant Biology, Michigan State University, East Lansing, MI
| | | | - Corlett W Wood
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, PA
| | - Daniel J Schoen
- Department of Biology, McGill University, Montreal, QC, Canada
| | - Stephen I Wright
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - John R Stinchcombe
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| |
Collapse
|
19
|
Uzunović J, Josephs EB, Stinchcombe JR, Wright SI. Transposable Elements Are Important Contributors to Standing Variation in Gene Expression in Capsella Grandiflora. Mol Biol Evol 2020; 36:1734-1745. [PMID: 31028401 DOI: 10.1093/molbev/msz098] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Transposable elements (TEs) make up a significant portion of eukaryotic genomes and are important drivers of genome evolution. However, the extent to which TEs affect gene expression variation on a genome-wide scale in comparison with other types of variants is still unclear. We characterized TE insertion polymorphisms and their association with gene expression in 124 whole-genome sequences from a single population of Capsella grandiflora, and contrasted this with the effects of single nucleotide polymorphisms (SNPs). Population frequency of insertions was negatively correlated with distance to genes, as well as density of conserved noncoding elements, suggesting that the negative effects of TEs on gene regulation are important in limiting their abundance. Rare TE variants strongly influence gene expression variation, predominantly through downregulation. In contrast, rare SNPs contribute equally to up- and down-regulation, but have a weaker individual effect than TEs. An expression quantitative trait loci (eQTL) analysis shows that a greater proportion of common TEs are eQTLs as opposed to common SNPs, and a third of the genes with TE eQTLs do not have SNP eQTLs. In contrast with rare TE insertions, common insertions are more likely to increase expression, consistent with recent models of cis-regulatory evolution favoring enhancer alleles. Taken together, these results imply that TEs are a significant contributor to gene expression variation and are individually more likely than rare SNPs to cause extreme changes in gene expression.
Collapse
Affiliation(s)
- Jasmina Uzunović
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
| | - Emily B Josephs
- Department of Plant Biology, Michigan State University, East Lansing, MI
| | - John R Stinchcombe
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada.,Koffler Scientific Reserve, University of Toronto, Toronto, Ontario, Canada
| | - Stephen I Wright
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada.,Centre for the Analysis of Genome Evolution and Function, University of Toronto, Toronto, Ontario, Canada
| |
Collapse
|
20
|
Abstract
Adaptive mutations play an important role in molecular evolution. However, the frequency and nature of these mutations at the intramolecular level are poorly understood. To address this, we analyzed the impact of protein architecture on the rate of adaptive substitutions, aiming to understand how protein biophysics influences fitness and adaptation. Using Drosophila melanogaster and Arabidopsis thaliana population genomics data, we fitted models of distribution of fitness effects and estimated the rate of adaptive amino-acid substitutions both at the protein and amino-acid residue level. We performed a comprehensive analysis covering genome, gene, and protein structure, by exploring a multitude of factors with a plausible impact on the rate of adaptive evolution, such as intron number, protein length, secondary structure, relative solvent accessibility, intrinsic protein disorder, chaperone affinity, gene expression, protein function, and protein-protein interactions. We found that the relative solvent accessibility is a major determinant of adaptive evolution, with most adaptive mutations occurring at the surface of proteins. Moreover, we observe that the rate of adaptive substitutions differs between protein functional classes, with genes encoding for protein biosynthesis and degradation signaling exhibiting the fastest rates of protein adaptation. Overall, our results suggest that adaptive evolution in proteins is mainly driven by intermolecular interactions, with host-pathogen coevolution likely playing a major role.
Collapse
Affiliation(s)
- Ana Filipa Moutinho
- Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Fernanda Fontes Trancoso
- Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Julien Yann Dutheil
- Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Biology, Plön, Germany.,Unité Mixte de Recherche 5554 Institut des Sciences de l'Evolution, CNRS, IRD, EPHE, Université de Montpellier, Montpellier, France
| |
Collapse
|
21
|
Mattila TM, Laenen B, Slotte T. Population Genomics of Transitions to Selfing in Brassicaceae Model Systems. Methods Mol Biol 2020; 2090:269-287. [PMID: 31975171 DOI: 10.1007/978-1-0716-0199-0_11] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Many plants harbor complex mechanisms that promote outcrossing and efficient pollen transfer. These include floral adaptations as well as genetic mechanisms, such as molecular self-incompatibility (SI) systems. The maintenance of such systems over long evolutionary timescales suggests that outcrossing is favorable over a broad range of conditions. Conversely, SI has repeatedly been lost, often in association with transitions to self-fertilization (selfing). This transition is favored when the short-term advantages of selfing outweigh the costs, primarily inbreeding depression. The transition to selfing is expected to have major effects on population genetic variation and adaptive potential, as well as on genome evolution. In the Brassicaceae, many studies on the population genetic, gene regulatory, and genomic effects of selfing have centered on the model plant Arabidopsis thaliana and the crucifer genus Capsella. The accumulation of population genomics datasets have allowed detailed investigation of where, when and how the transition to selfing occurred. Future studies will take advantage of the development of population genetics theory on the impact of selfing, especially regarding positive selection. Furthermore, investigation of systems including recent transitions to selfing, mixed mating populations and/or multiple independent replicates of the same transition will facilitate dissecting the effects of mating system variation from processes driven by demography.
Collapse
Affiliation(s)
- Tiina M Mattila
- Department of Ecology and Genetics, University of Oulu, Oulu, Finland
| | - Benjamin Laenen
- Department of Ecology, Environment, and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Tanja Slotte
- Department of Ecology, Environment, and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden.
| |
Collapse
|
22
|
Moutinho AF, Bataillon T, Dutheil JY. Variation of the adaptive substitution rate between species and within genomes. Evol Ecol 2019. [DOI: 10.1007/s10682-019-10026-z] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
AbstractThe importance of adaptive mutations in molecular evolution is extensively debated. Recent developments in population genomics allow inferring rates of adaptive mutations by fitting a distribution of fitness effects to the observed patterns of polymorphism and divergence at sites under selection and sites assumed to evolve neutrally. Here, we summarize the current state-of-the-art of these methods and review the factors that affect the molecular rate of adaptation. Several studies have reported extensive cross-species variation in the proportion of adaptive amino-acid substitutions (α) and predicted that species with larger effective population sizes undergo less genetic drift and higher rates of adaptation. Disentangling the rates of positive and negative selection, however, revealed that mutations with deleterious effects are the main driver of this population size effect and that adaptive substitution rates vary comparatively little across species. Conversely, rates of adaptive substitution have been documented to vary substantially within genomes. On a genome-wide scale, gene density, recombination and mutation rate were observed to play a role in shaping molecular rates of adaptation, as predicted under models of linked selection. At the gene level, it has been reported that the gene functional category and the macromolecular structure substantially impact the rate of adaptive mutations. Here, we deliver a comprehensive review of methods used to infer the molecular adaptive rate, the potential drivers of adaptive evolution and how positive selection shapes molecular evolution within genes, across genes within species and between species.
Collapse
|
23
|
Ho EKH, Bartkowska M, Wright SI, Agrawal AF. Population genomics of the facultatively asexual duckweed Spirodela polyrhiza. THE NEW PHYTOLOGIST 2019; 224:1361-1371. [PMID: 31298732 DOI: 10.1111/nph.16056] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2019] [Accepted: 06/18/2019] [Indexed: 06/10/2023]
Abstract
Clonal propagation allows some plant species to achieve massive population sizes quickly but also reduces the evolutionary independence of different sites in the genome. We examine genome-wide genetic diversity in Spirodela polyrhiza, a duckweed that reproduces primarily asexually. We find that this geographically widespread and numerically abundant species has very low levels of genetic diversity. Diversity at nonsynonymous sites relative to synonymous sites is high, suggesting that purifying selection is weak. A potential explanation for this observation is that a very low frequency of sex renders selection ineffective. However, there is a pronounced decay in linkage disequilibrium over 40 kb, suggesting that though sex may be rare at the individual level it is not too infrequent at the population level. In addition, neutral diversity is affected by the physical proximity of selected sites, which would be unexpected if sex was exceedingly rare at the population level. The amount of genetic mixing as assessed by the decay in linkage disequilibrium is not dissimilar from selfing species such as Arabidopsis thaliana, yet selection appears to be much less effective in duckweed. We discuss alternative explanations for the signature of weak purifying selection.
Collapse
Affiliation(s)
- Eddie K H Ho
- Department of Ecology and Evolutionary Biology, University of Toronto, 25 Willcocks Street, Toronto, ON, M5S 3B2, Canada
| | - Magdalena Bartkowska
- Department of Ecology and Evolutionary Biology, University of Toronto, 25 Willcocks Street, Toronto, ON, M5S 3B2, Canada
| | - Stephen I Wright
- Department of Ecology and Evolutionary Biology, University of Toronto, 25 Willcocks Street, Toronto, ON, M5S 3B2, Canada
- Center for Analysis of Genome Evolution and Function, University of Toronto, 25 Willcocks Street, Toronto, ON, M5S 3B2, Canada
| | - Aneil F Agrawal
- Department of Ecology and Evolutionary Biology, University of Toronto, 25 Willcocks Street, Toronto, ON, M5S 3B2, Canada
- Center for Analysis of Genome Evolution and Function, University of Toronto, 25 Willcocks Street, Toronto, ON, M5S 3B2, Canada
| |
Collapse
|
24
|
Mattila TM, Laenen B, Horvath R, Hämälä T, Savolainen O, Slotte T. Impact of demography on linked selection in two outcrossing Brassicaceae species. Ecol Evol 2019; 9:9532-9545. [PMID: 31534673 PMCID: PMC6745670 DOI: 10.1002/ece3.5463] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 06/28/2019] [Accepted: 07/02/2019] [Indexed: 12/13/2022] Open
Abstract
Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.
Collapse
Affiliation(s)
- Tiina M. Mattila
- Department of Ecology and GeneticsUniversity of OuluOuluFinland
- Present address:
Department of Organismal BiologyUppsala UniversityUppsalaSweden
| | - Benjamin Laenen
- Science for Life Laboratory, Department of Ecology, Environment, and Plant SciencesStockholm UniversityStockholmSweden
| | - Robert Horvath
- Science for Life Laboratory, Department of Ecology, Environment, and Plant SciencesStockholm UniversityStockholmSweden
| | - Tuomas Hämälä
- Department of Ecology and GeneticsUniversity of OuluOuluFinland
- Biocenter OuluUniversity of OuluOuluFinland
- Present address:
Department of Plant and Microbial BiologyUniversity of Minnesota Twin CitiesSt. PaulMNUSA
| | - Outi Savolainen
- Department of Ecology and GeneticsUniversity of OuluOuluFinland
- Biocenter OuluUniversity of OuluOuluFinland
| | - Tanja Slotte
- Science for Life Laboratory, Department of Ecology, Environment, and Plant SciencesStockholm UniversityStockholmSweden
| |
Collapse
|
25
|
Braasch J, Barker BS, Dlugosch KM. Expansion history and environmental suitability shape effective population size in a plant invasion. Mol Ecol 2019; 28:2546-2558. [PMID: 30993767 DOI: 10.1111/mec.15104] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2018] [Revised: 04/02/2019] [Accepted: 04/04/2019] [Indexed: 12/19/2022]
Abstract
The margins of an expanding range are predicted to be challenging environments for adaptation. Marginal populations should often experience low effective population sizes (Ne ) where genetic drift is high due to demographic expansion and/or census population size is low due to unfavourable environmental conditions. Nevertheless, invasive species demonstrate increasing evidence of rapid evolution and potential adaptation to novel environments encountered during colonization, calling into question whether significant reductions in Ne are realized during range expansions in nature. Here we report one of the first empirical tests of the joint effects of expansion dynamics and environment on effective population size variation during invasive range expansion. We estimate contemporary values of Ne using rates of linkage disequilibrium among genome-wide markers within introduced populations of the highly invasive plant Centaurea solstitialis (yellow starthistle) in North America (California, USA), and within native Eurasian populations. As predicted, we find that Ne within the invaded range is positively correlated with both expansion history (time since founding) and habitat quality (abiotic climate). History and climate had independent additive effects with similar effect sizes, indicating an important role for both factors in this invasion. These results support theoretical expectations for the population genetics of range expansion, though whether these processes can ultimately arrest the spread of an invasive species remains an unanswered question.
Collapse
Affiliation(s)
- Joseph Braasch
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona
| | - Brittany S Barker
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona.,Integrated Plant Protection Center and Department of Horticulture, Oregon State University, Corvallis, Oregon
| | - Katrina M Dlugosch
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona
| |
Collapse
|
26
|
Deleterious Mutation Burden and Its Association with Complex Traits in Sorghum ( Sorghum bicolor). Genetics 2019; 211:1075-1087. [PMID: 30622134 PMCID: PMC6404259 DOI: 10.1534/genetics.118.301742] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Accepted: 12/22/2018] [Indexed: 11/18/2022] Open
Abstract
Sorghum (Sorghum bicolor (L.) Moench) is a major staple food cereal for millions of people worldwide. Valluru et al. identify putative deleterious mutations among ∼5.5M segregating variants of 229 diverse sorghum... Sorghum (Sorghum bicolor L.) is a major food cereal for millions of people worldwide. The sorghum genome, like other species, accumulates deleterious mutations, likely impacting its fitness. The lack of recombination, drift, and the coupling with favorable loci impede the removal of deleterious mutations from the genome by selection. To study how deleterious variants impact phenotypes, we identified putative deleterious mutations among ∼5.5 M segregating variants of 229 diverse biomass sorghum lines. We provide the whole-genome estimate of the deleterious burden in sorghum, showing that ∼33% of nonsynonymous substitutions are putatively deleterious. The pattern of mutation burden varies appreciably among racial groups. Across racial groups, the mutation burden correlated negatively with biomass, plant height, specific leaf area (SLA), and tissue starch content (TSC), suggesting that deleterious burden decreases trait fitness. Putatively deleterious variants explain roughly one-half of the genetic variance. However, there is only moderate improvement in total heritable variance explained for biomass (7.6%) and plant height (average of 3.1% across all stages). There is no advantage in total heritable variance for SLA and TSC. The contribution of putatively deleterious variants to phenotypic diversity therefore appears to be dependent on the genetic architecture of traits. Overall, these results suggest that incorporating putatively deleterious variants into genomic models slightly improves prediction accuracy because of extensive linkage. Knowledge of deleterious variants could be leveraged for sorghum breeding through either genome editing and/or conventional breeding that focuses on the selection of progeny with fewer deleterious alleles.
Collapse
|
27
|
Abstract
A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.
Collapse
Affiliation(s)
- Sylvain Glémin
- Institut des Sciences de l'Evolution, UMR5554, Université Montpellier II, Montpellier, France
| | - Clémentine M François
- Institut des Sciences de l'Evolution, UMR5554, Université Montpellier II, Montpellier, France
| | - Nicolas Galtier
- Institut des Sciences de l'Evolution, UMR5554, Université Montpellier II, Montpellier, France.
| |
Collapse
|
28
|
Horvath R, Slotte T. The Role of Small RNA-Based Epigenetic Silencing for Purifying Selection on Transposable Elements in Capsella grandiflora. Genome Biol Evol 2018; 9:2911-2920. [PMID: 29036316 PMCID: PMC5737465 DOI: 10.1093/gbe/evx206] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/27/2017] [Indexed: 12/13/2022] Open
Abstract
To avoid negative effects of transposable element (TE) proliferation, plants epigenetically silence TEs using a number of mechanisms, including RNA-directed DNA methylation. These epigenetic modifications can extend outside the boundaries of TE insertions and lead to silencing of nearby genes, resulting in a trade-off between TE silencing and interference with nearby gene regulation. Therefore, purifying selection is expected to remove silenced TE insertions near genes more efficiently and prevent their accumulation within a population. To explore how effects of TE silencing on gene regulation shapes purifying selection on TEs, we analyzed whole genome sequencing data from 166 individuals of a large population of the outcrossing species Capsella grandiflora. We found that most TEs are rare, and in chromosome arms, silenced TEs are exposed to stronger purifying selection than those that are not silenced by 24-nucleotide small RNAs, especially with increasing proximity to genes. An age-of-allele test of neutrality on a subset of TEs supports our inference of purifying selection on silenced TEs, suggesting that our results are robust to varying transposition rates. Our results provide new insights into the processes affecting the accumulation of TEs in an outcrossing species and support the view that epigenetic silencing of TEs results in a trade-off between preventing TE proliferation and interference with nearby gene regulation. We also suggest that in the centromeric and pericentromeric regions, the negative aspects of epigenetic TE silencing are missing.
Collapse
Affiliation(s)
- Robert Horvath
- Science for Life Laboratory, Department of Ecology, Environment and Plant Sciences, Stockholm University, Sweden
| | - Tanja Slotte
- Science for Life Laboratory, Department of Ecology, Environment and Plant Sciences, Stockholm University, Sweden
| |
Collapse
|
29
|
Paape T, Briskine RV, Halstead-Nussloch G, Lischer HEL, Shimizu-Inatsugi R, Hatakeyama M, Tanaka K, Nishiyama T, Sabirov R, Sese J, Shimizu KK. Patterns of polymorphism and selection in the subgenomes of the allopolyploid Arabidopsis kamchatica. Nat Commun 2018; 9:3909. [PMID: 30254374 PMCID: PMC6156220 DOI: 10.1038/s41467-018-06108-1] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2018] [Accepted: 08/10/2018] [Indexed: 12/30/2022] Open
Abstract
Genome duplication is widespread in wild and crop plants. However, little is known about genome-wide selection in polyploids due to the complexity of duplicated genomes. In polyploids, the patterns of purifying selection and adaptive substitutions may be affected by masking owing to duplicated genes or homeologs as well as effective population size. Here, we resequence 25 accessions of the allotetraploid Arabidopsis kamchatica, which is derived from the diploid species A. halleri and A. lyrata. We observe a reduction in purifying selection compared with the parental species. Interestingly, proportions of adaptive non-synonymous substitutions are significantly positive in contrast to most plant species. A recurrent pattern observed in both frequency and divergence–diversity neutrality tests is that the genome-wide distributions of both subgenomes are similar, but the correlation between homeologous pairs is low. This may increase the opportunity of different evolutionary trajectories such as in the HMA4 gene involved in heavy metal hyperaccumulation. Despite the prevalence of genome duplication in plants, little is known about the evolutionary patterns of entire subgenomes. Here the authors resequence allopolyploid Arabidopsis kamchatica genome to estimate diversity, linkage disequilibrium and strengths of both positive and purifying selection.
Collapse
Affiliation(s)
- Timothy Paape
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland. .,Department of Plant and Microbial Biology, University of Zurich, Zollikerstrasse 107, CH-8008, Zurich, Switzerland.
| | - Roman V Briskine
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland.,Department of Environmental Systems Science, ETH Zurich, CH-8092, Zurich, Switzerland.,Functional Genomics Center Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland
| | - Gwyneth Halstead-Nussloch
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland
| | - Heidi E L Lischer
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland.,Swiss Institute of Bioinformatics (SIB), Lausanne, 1015, Switzerland
| | - Rie Shimizu-Inatsugi
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland.,Department of Plant and Microbial Biology, University of Zurich, Zollikerstrasse 107, CH-8008, Zurich, Switzerland
| | - Masaomi Hatakeyama
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland.,Department of Plant and Microbial Biology, University of Zurich, Zollikerstrasse 107, CH-8008, Zurich, Switzerland.,Swiss Institute of Bioinformatics (SIB), Lausanne, 1015, Switzerland.,Functional Genomics Center Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland
| | - Kenta Tanaka
- Sugadaira Montane Research Center, University of Tsukuba, Nagano, Ueda, 386-2204, Japan
| | - Tomoaki Nishiyama
- Advanced Science Research Center, Kanazawa University, 13-1 Takara-machi, Kanazawa, 920-0934, Japan
| | - Renat Sabirov
- Institute of Marine Geology and Geophysics, Far East Branch, Russian Academy of Sciences, Nauki street, 1-B, Yuzhno-Sakhalinsk, 693022, Russian Federation
| | - Jun Sese
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, 135-0064, Japan.,AIST-Tokyo Tech Real World Big-Data Computation Open Innovation Laboratory, Tokyo, 152-8550, Japan
| | - Kentaro K Shimizu
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland. .,Department of Plant and Microbial Biology, University of Zurich, Zollikerstrasse 107, CH-8008, Zurich, Switzerland. .,Kihara Institute for Biological Research, Yokohama City University, 641-12 Maioka, Yokohama, 244-0813, Japan.
| |
Collapse
|
30
|
Ho EKH, Agrawal AF. Mutation accumulation in selfing populations under fluctuating selection. Evolution 2018; 72:1759-1772. [DOI: 10.1111/evo.13553] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2018] [Revised: 06/26/2018] [Accepted: 07/01/2018] [Indexed: 12/13/2022]
Affiliation(s)
- Eddie K. H. Ho
- Department of Ecology and Evolutionary Biology University of Toronto 25 Willcocks Street Toronto ON M5S 3B2 Canada
| | - Aneil F. Agrawal
- Department of Ecology and Evolutionary Biology University of Toronto 25 Willcocks Street Toronto ON M5S 3B2 Canada
| |
Collapse
|
31
|
Petrone Mendoza S, Lascoux M, Glémin S. Competitive ability of Capsella species with different mating systems and ploidy levels. ANNALS OF BOTANY 2018; 121:1257-1264. [PMID: 29471370 PMCID: PMC5946883 DOI: 10.1093/aob/mcy014] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Accepted: 01/19/2018] [Indexed: 05/31/2023]
Abstract
BACKGROUND AND AIMS Capsella is a model genus for studying the transition from outcrossing to selfing, with or without change in ploidy levels. The genomic consequences and changes in reproductive traits (selfing syndrome) associated with these shifts have been studied in depth. However, potential ecological divergence among species of the genus has not been determined. Among ecological traits, competitive ability could be relevant for selfing evolution, as selfing has been shown to be statistically associated with reduced competitiveness in a recent meta-analysis. METHODS We assessed the effect of competition on three Capsella species differing in their mating system and ploidy level. We used an experimental design where fitness related traits were measured in focal individuals with and without competitors. KEY RESULTS The diploid selfer (C. rubella) was most sensitive to competition, whereas the tetraploid selfer (C. bursa-pastoris) performed the best, with the diploid outcrosser (C. grandiflora) being intermediate. CONCLUSIONS These results add to the detailed characterization of Capsella species and highlight the possible roles of ecological context and ploidy in the evolutionary trajectories of selfing species.
Collapse
Affiliation(s)
- Sandra Petrone Mendoza
- Department of Ecology and Genetics, Evolutionary Biology Center and Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Martin Lascoux
- Department of Ecology and Genetics, Evolutionary Biology Center and Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Sylvain Glémin
- Department of Ecology and Genetics, Evolutionary Biology Center and Science for Life Laboratory, Uppsala University, Uppsala, Sweden
- Institut des Sciences de l’Evolution (ISEM - UMR 5554 Université de Montpellier-CNRS-IRD-EPHE), France
| |
Collapse
|
32
|
Grivet D, Avia K, Vaattovaara A, Eckert AJ, Neale DB, Savolainen O, González-Martínez SC. High rate of adaptive evolution in two widespread European pines. Mol Ecol 2017; 26:6857-6870. [PMID: 29110402 DOI: 10.1111/mec.14402] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2015] [Revised: 09/14/2017] [Accepted: 09/25/2017] [Indexed: 12/18/2022]
Abstract
Comparing related organisms with differing ecological requirements and evolutionary histories can shed light on the mechanisms and drivers underlying genetic adaptation. Here, by examining a common set of hundreds of loci, we compare patterns of nucleotide diversity and molecular adaptation of two European conifers (Scots pine and maritime pine) living in contrasted environments and characterized by distinct population genetic structure (low and clinal in Scots pine, high and ecotypic in maritime pine) and demographic histories. We found higher nucleotide diversity in Scots pine than in maritime pine, whereas rates of new adaptive substitutions (ωa ), as estimated from the distribution of fitness effects, were similar across species and among the highest found in plants. Sample size and population genetic structure did not appear to have resulted in significant bias in estimates of ωa . Moreover, population contraction-expansion dynamics for each species did not affect differentially the rate of adaptive substitution in these two pines. Several methodological and biological factors may underlie the unusually high rate of adaptive evolution of Scots pine and maritime pine. By providing two new case studies with contrasting evolutionary histories, we contribute to disentangling the multiple factors potentially affecting adaptive evolution in natural plant populations.
Collapse
Affiliation(s)
- Delphine Grivet
- Department of Forest Ecology and Genetics, Forest Research Centre, INIA-CIFOR, Madrid, Spain.,Sustainable Forest Management Research Institute, INIA - University of Valladolid, Palencia, Spain
| | - Komlan Avia
- Department of Ecology and Genetics and Biocenter Oulu, University of Oulu, Oulu, Finland.,Algal Genetics Group, UMR 8227, CNRS, Sorbonne Universités, UPMC, Station Biologique Roscoff, Roscoff, France.,UMI 3614 Evolutionary Biology and Ecology of Algae, CNRS, Sorbonne Universités, UPMC, Pontificia Universidad Católica de Chile, Universidad Austral de Chile, Station Biologique Roscoff, Roscoff, France
| | - Aleksia Vaattovaara
- Department of Ecology and Genetics and Biocenter Oulu, University of Oulu, Oulu, Finland.,Division of Plant Biology, Department of Biosciences, Viikki Plant Science Centre (ViPS), University of Helsinki, Helsinki, Finland
| | - Andrew J Eckert
- Department of Biology, Virginia Commonwealth University, Richmond, VA, USA
| | - David B Neale
- Department of Plant Sciences, University of California at Davis, Davis, CA, USA
| | - Outi Savolainen
- Department of Ecology and Genetics and Biocenter Oulu, University of Oulu, Oulu, Finland
| | - Santiago C González-Martínez
- Department of Forest Ecology and Genetics, Forest Research Centre, INIA-CIFOR, Madrid, Spain.,Sustainable Forest Management Research Institute, INIA - University of Valladolid, Palencia, Spain.,BIOGECO, INRA, Univ. Bordeaux, Cestas, France
| |
Collapse
|
33
|
Jiang X, Assis R. Natural Selection Drives Rapid Functional Evolution of Young Drosophila Duplicate Genes. Mol Biol Evol 2017; 34:3089-3098. [PMID: 28961791 PMCID: PMC5850746 DOI: 10.1093/molbev/msx230] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Gene duplication is thought to play a major role in phenotypic evolution. Yet the forces involved in the functional divergence of young duplicate genes remain unclear. Here, we use population-genetic inference to elucidate the role of natural selection in the functional evolution of young duplicate genes in Drosophila melanogaster. We find that negative selection acts on young duplicates with ancestral functions, and positive selection on those with novel functions, suggesting that natural selection may determine whether and how young duplicate genes are retained. Moreover, evidence of natural selection is strongest in protein-coding regions and 3' UTRs of young duplicates, indicating that selection may primarily target encoded proteins and regulatory sequences specific to 3' UTRs. Further analysis reveals that natural selection acts immediately after duplication and weakens over time, possibly explaining the observed bias toward the acquisition of new functions by young, rather than old, duplicate gene copies. Last, we find an enrichment of testis-related functions in young duplicates that underwent recent positive selection, but not in young duplicates that did not undergo recent positive selection, or in old duplicates that either did or did not undergo recent positive selection. Thus, our findings reveal that natural selection is a key player in the functional evolution of young duplicate genes, acts rapidly and in a region-specific manner, and may underlie the origin of novel testis-specific phenotypes in Drosophila.
Collapse
Affiliation(s)
- Xueyuan Jiang
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA
| | - Raquel Assis
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA
- Department of Biology, Pennsylvania State University, University Park, PA
| |
Collapse
|
34
|
Hawkins AK, Garza ER, Dietz VA, Hernandez OJ, Hawkins WD, Burrell AM, Pepper AE. Transcriptome Signatures of Selection, Drift, Introgression, and Gene Duplication in the Evolution of an Extremophile Endemic Plant. Genome Biol Evol 2017; 9:3478-3494. [PMID: 29220486 PMCID: PMC5751042 DOI: 10.1093/gbe/evx259] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/04/2017] [Indexed: 12/26/2022] Open
Abstract
Plants on serpentine soils provide extreme examples of adaptation to environment, and thus offer excellent models for the study of evolution at the molecular and genomic level. Serpentine outcrops are derived from ultramafic rock and have extremely low levels of essential plant nutrients (e.g., N, P, K, and Ca), as well as toxic levels of heavy metals (e.g., Ni, Cr, and Co) and low moisture availability. These outcrops provide habitat to a number of endemic plant species, including the annual mustard Caulanthus amplexicaulis var. barbarae (Cab) (Brassicaceae). Its sister taxon, C. amplexicaulis var. amplexicaulis (Caa), is intolerant to serpentine soils. Here, we assembled and annotated comprehensive reference transcriptomes of both Caa and Cab for use in protein coding sequence comparisons. A set of 29,443 reciprocal best Blast hit (RBH) orthologs between Caa and Cab was compared with identify coding sequence variants, revealing a high genome-wide dN/dS ratio between the two taxa (mean = 0.346). We show that elevated dN/dS likely results from the composite effects of genetic drift, positive selection, and the relaxation of negative selection. Further, analysis of paralogs within each taxon revealed the signature of a period of elevated gene duplication (∼10 Ma) that is shared with other species of the tribe Thelypodieae, and may have played a role in the striking morphological and ecological diversity of this tribe. In addition, distribution of the synonymous substitution rate, dS, is strongly bimodal, indicating a history of reticulate evolution that may have contributed to serpentine adaptation.
Collapse
|
35
|
De La Torre AR, Li Z, Van de Peer Y, Ingvarsson PK. Contrasting Rates of Molecular Evolution and Patterns of Selection among Gymnosperms and Flowering Plants. Mol Biol Evol 2017; 34:1363-1377. [PMID: 28333233 PMCID: PMC5435085 DOI: 10.1093/molbev/msx069] [Citation(s) in RCA: 111] [Impact Index Per Article: 15.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The majority of variation in rates of molecular evolution among seed plants remains both unexplored and unexplained. Although some attention has been given to flowering plants, reports of molecular evolutionary rates for their sister plant clade (gymnosperms) are scarce, and to our knowledge differences in molecular evolution among seed plant clades have never been tested in a phylogenetic framework. Angiosperms and gymnosperms differ in a number of features, of which contrasting reproductive biology, life spans, and population sizes are the most prominent. The highly conserved morphology of gymnosperms evidenced by similarity of extant species to fossil records and the high levels of macrosynteny at the genomic level have led scientists to believe that gymnosperms are slow-evolving plants, although some studies have offered contradictory results. Here, we used 31,968 nucleotide sites obtained from orthologous genes across a wide taxonomic sampling that includes representatives of most conifers, cycads, ginkgo, and many angiosperms with a sequenced genome. Our results suggest that angiosperms and gymnosperms differ considerably in their rates of molecular evolution per unit time, with gymnosperm rates being, on average, seven times lower than angiosperm species. Longer generation times and larger genome sizes are some of the factors explaining the slow rates of molecular evolution found in gymnosperms. In contrast to their slow rates of molecular evolution, gymnosperms possess higher substitution rate ratios than angiosperm taxa. Finally, our study suggests stronger and more efficient purifying and diversifying selection in gymnosperm than in angiosperm species, probably in relation to larger effective population sizes.
Collapse
Affiliation(s)
- Amanda R De La Torre
- Department of Plant Sciences, University of California-Davis, Davis, CA.,Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden
| | - Zhen Li
- Department of Plant Systems Biology, VIB, Ghent, Belgium.,Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
| | - Yves Van de Peer
- Department of Plant Systems Biology, VIB, Ghent, Belgium.,Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.,Genomics Research Institute, University of Pretoria, Hatfield Campus, Pretoria, South Africa
| | - Pär K Ingvarsson
- Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden.,Department of Plant Biology, Uppsala Biocenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
36
|
Inference of Distribution of Fitness Effects and Proportion of Adaptive Substitutions from Polymorphism Data. Genetics 2017; 207:1103-1119. [PMID: 28951530 PMCID: PMC5676230 DOI: 10.1534/genetics.117.300323] [Citation(s) in RCA: 87] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2017] [Accepted: 09/13/2017] [Indexed: 11/18/2022] Open
Abstract
The distribution of fitness effects (DFE) encompasses the fraction of deleterious, neutral, and beneficial mutations. It conditions the evolutionary trajectory of populations, as well as the rate of adaptive molecular evolution (α). Inferring DFE and α from patterns of polymorphism, as given through the site frequency spectrum (SFS) and divergence data, has been a longstanding goal of evolutionary genetics. A widespread assumption shared by previous inference methods is that beneficial mutations only contribute negligibly to the polymorphism data. Hence, a DFE comprising only deleterious mutations tends to be estimated from SFS data, and α is then predicted by contrasting the SFS with divergence data from an outgroup. We develop a hierarchical probabilistic framework that extends previous methods to infer DFE and α from polymorphism data alone. We use extensive simulations to examine the performance of our method. While an outgroup is still needed to obtain an unfolded SFS, we show that both a DFE, comprising both deleterious and beneficial mutations, and α can be inferred without using divergence data. We also show that not accounting for the contribution of beneficial mutations to polymorphism data leads to substantially biased estimates of the DFE and α. We compare our framework with one of the most widely used inference methods available and apply it on a recently published chimpanzee exome data set.
Collapse
|
37
|
Abstract
Mistranslation errors compromise fitness by wasting resources on nonfunctional proteins. In order to reduce the cost of mistranslations, natural selection chooses the most accurately translated codons at sites that are particularly important for protein structure and function. We investigated the determinants underlying selection for translational accuracy in several species of plants belonging to three clades: Brassicaceae, Fabidae, and Poaceae. Although signatures of translational selection were found in genes from a wide range of species, the underlying factors varied in nature and intensity. Indeed, the degree of synonymous codon bias at evolutionarily conserved sites varied among plant clades while remaining uniform within each clade. This is unlikely to solely reflect the diversity of tRNA pools because there is little correlation between synonymous codon bias and tRNA abundance, so other factors must affect codon choice and translational accuracy in plant genes. Accordingly, synonymous codon choice at a given site was affected not only by the selection pressure at that site, but also its participation in protein domains or mRNA secondary structures. Although these effects were detected in all the species we analyzed, their impact on translation accuracy was distinct in evolutionarily distant plant clades. The domain effect was found to enhance translational accuracy in dicot and monocot genes with a high GC content, but to oppose the selection of more accurate codons in monocot genes with a low GC content.
Collapse
|
38
|
Balao F, Trucchi E, Wolfe TM, Hao B, Lorenzo MT, Baar J, Sedman L, Kosiol C, Amman F, Chase MW, Hedrén M, Paun O. Adaptive sequence evolution is driven by biotic stress in a pair of orchid species (Dactylorhiza) with distinct ecological optima. Mol Ecol 2017; 26:3649-3662. [PMID: 28370647 PMCID: PMC5518283 DOI: 10.1111/mec.14123] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2016] [Revised: 02/24/2017] [Accepted: 03/15/2017] [Indexed: 12/23/2022]
Abstract
The orchid family is the largest in the angiosperms, but little is known about the molecular basis of the significant variation they exhibit. We investigate here the transcriptomic divergence between two European terrestrial orchids, Dactylorhiza incarnata and Dactylorhiza fuchsii, and integrate these results in the context of their distinct ecologies that we also document. Clear signals of lineage-specific adaptive evolution of protein-coding sequences are identified, notably targeting elements of biotic defence, including both physical and chemical adaptations in the context of divergent pools of pathogens and herbivores. In turn, a substantial regulatory divergence between the two species appears linked to adaptation/acclimation to abiotic conditions. Several of the pathways affected by differential expression are also targeted by deviating post-transcriptional regulation via sRNAs. Finally, D. incarnata appears to suffer from insufficient sRNA control over the activity of RNA-dependent DNA polymerase, resulting in increased activity of class I transposable elements and, over time, in larger genome size than that of D. fuchsii. The extensive molecular divergence between the two species suggests significant genomic and transcriptomic shock in their hybrids and offers insights into the difficulty of coexistence at the homoploid level. Altogether, biological response to selection, accumulated during the history of these orchids, appears governed by their microenvironmental context, in which biotic and abiotic pressures act synergistically to shape transcriptome structure, expression and regulation.
Collapse
Affiliation(s)
- Francisco Balao
- Department of Botany and Biodiversity ResearchUniversity of ViennaViennaAustria
- Departamento de Biología Vegetal y EcologíaUniversity of SevilleSevillaSpain
| | - Emiliano Trucchi
- Department of Botany and Biodiversity ResearchUniversity of ViennaViennaAustria
- Department of Life Sciences and BiotechnologiesUniversity of FerraraFerraraItaly
| | - Thomas M. Wolfe
- Department of Botany and Biodiversity ResearchUniversity of ViennaViennaAustria
- Vienna Graduate School of Population GeneticsViennaAustria
| | - Bao‐Hai Hao
- Department of Botany and Biodiversity ResearchUniversity of ViennaViennaAustria
| | - Maria Teresa Lorenzo
- Department of Botany and Biodiversity ResearchUniversity of ViennaViennaAustria
- Departamento de Biología Vegetal y EcologíaUniversity of SevilleSevillaSpain
| | - Juliane Baar
- Department of Botany and Biodiversity ResearchUniversity of ViennaViennaAustria
| | - Laura Sedman
- Gregor Mendel Institute for Plant Molecular BiologyViennaAustria
| | - Carolin Kosiol
- Institut für PopulationsgenetikVetmeduni ViennaViennaAustria
- Centre of Biological DiversitySchool of BiologyUniversity of St AndrewsSt AndrewsUK
| | - Fabian Amman
- Department of Chromosome BiologyUniversity of ViennaViennaAustria
| | - Mark W. Chase
- Royal Botanic Gardens KewRichmondUK
- School of Plant BiologyUniversity of Western AustraliaCrawley, PerthWAAustralia
| | | | - Ovidiu Paun
- Department of Botany and Biodiversity ResearchUniversity of ViennaViennaAustria
| |
Collapse
|
39
|
Ho EKH, Agrawal AF. Aging asexual lineages and the evolutionary maintenance of sex. Evolution 2017; 71:1865-1875. [PMID: 28444897 DOI: 10.1111/evo.13260] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2016] [Revised: 04/13/2017] [Accepted: 04/17/2017] [Indexed: 12/24/2022]
Abstract
Finite populations of asexual and highly selfing species suffer from a reduced efficacy of selection. Such populations are thought to decline in fitness over time due to accumulating slightly deleterious mutations or failing to adapt to changing conditions. These within-population processes that lead nonrecombining species to extinction may help maintain sex and outcrossing through species level selection. Although inefficient selection is proposed to elevate extinction rates over time, previous models of species selection for sex assumed constant diversification rates. For sex to persist, classic models require that asexual species diversify at rates lower than sexual species; the validity of this requirement is questionable, both conceptually and empirically. We extend past models by allowing asexual lineages to decline in diversification rates as they age, that is nonrecombining lineages "senesce" in diversification rates. At equilibrium, senescing diversification rates maintain sex even when asexual lineages, at young ages, diversify faster than their sexual progenitors. In such cases, the age distribution of asexual lineages contains a peak at intermediate values rather than showing the exponential decline predicted by the classic model. Coexistence requires only that the average rate of diversification in asexuals be lower than that of sexuals.
Collapse
Affiliation(s)
- Eddie K H Ho
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, M5S 3B2, Canada
| | - Aneil F Agrawal
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, M5S 3B2, Canada
| |
Collapse
|
40
|
Badouin H, Gladieux P, Gouzy J, Siguenza S, Aguileta G, Snirc A, Le Prieur S, Jeziorski C, Branca A, Giraud T. Widespread selective sweeps throughout the genome of model plant pathogenic fungi and identification of effector candidates. Mol Ecol 2017; 26:2041-2062. [DOI: 10.1111/mec.13976] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2016] [Revised: 12/15/2016] [Accepted: 12/19/2016] [Indexed: 12/11/2022]
Affiliation(s)
- H. Badouin
- Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech; Université Paris-Saclay; 91400 Orsay France
| | - P. Gladieux
- Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech; Université Paris-Saclay; 91400 Orsay France
- UMR BGPI; Campus International de Baillarguet; INRA; 34398 Montpellier France
| | - J. Gouzy
- Laboratoire des Interactions Plantes-Microorganismes (LIPM); UMR441; INRA; 31326 Castanet-Tolosan France
- Laboratoire des Interactions Plantes-Microorganismes (LIPM); UMR2594; CNRS; 31326 Castanet-Tolosan France
| | - S. Siguenza
- Laboratoire des Interactions Plantes-Microorganismes (LIPM); UMR441; INRA; 31326 Castanet-Tolosan France
- Laboratoire des Interactions Plantes-Microorganismes (LIPM); UMR2594; CNRS; 31326 Castanet-Tolosan France
| | - G. Aguileta
- Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech; Université Paris-Saclay; 91400 Orsay France
| | - A. Snirc
- Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech; Université Paris-Saclay; 91400 Orsay France
| | - S. Le Prieur
- Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech; Université Paris-Saclay; 91400 Orsay France
| | - C. Jeziorski
- Genotoul; GeT-PlaGe; INRA Auzeville 31326 Castanet-Tolosan France
- UAR1209; INRA Auzeville 31326 Castanet-Tolosan France
| | - A. Branca
- Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech; Université Paris-Saclay; 91400 Orsay France
| | - T. Giraud
- Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech; Université Paris-Saclay; 91400 Orsay France
| |
Collapse
|
41
|
Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora. Proc Natl Acad Sci U S A 2017; 114:1087-1092. [PMID: 28096395 DOI: 10.1073/pnas.1612561114] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.
Collapse
|
42
|
Braverman JM, Hamilton MB, Johnson BA. Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA. J Mol Evol 2016; 83:97-109. [PMID: 27592229 DOI: 10.1007/s00239-016-9752-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2015] [Accepted: 07/14/2016] [Indexed: 01/09/2023]
Abstract
There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima's 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26-32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.
Collapse
Affiliation(s)
- John M Braverman
- Department of Biology, Saint Joseph's University, Philadelphia, PA, USA.
| | | | - Brent A Johnson
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, NY, USA
| |
Collapse
|
43
|
Limits to Adaptation in Partially Selfing Species. Genetics 2016; 203:959-74. [PMID: 27098913 DOI: 10.1534/genetics.116.188821] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2016] [Accepted: 04/15/2016] [Indexed: 12/26/2022] Open
Abstract
In outcrossing populations, "Haldane's sieve" states that recessive beneficial alleles are less likely to fix than dominant ones, because they are less exposed to selection when rare. In contrast, selfing organisms are not subject to Haldane's sieve and are more likely to fix recessive types than outcrossers, as selfing rapidly creates homozygotes, increasing overall selection acting on mutations. However, longer homozygous tracts in selfers also reduce the ability of recombination to create new genotypes. It is unclear how these two effects influence overall adaptation rates in partially selfing organisms. Here, we calculate the fixation probability of beneficial alleles if there is an existing selective sweep in the population. We consider both the potential loss of the second beneficial mutation if it has a weaker advantage than the first one, and the possible replacement of the initial allele if the second mutant is fitter. Overall, loss of weaker adaptive alleles during a first selective sweep has a larger impact on preventing fixation of both mutations in highly selfing organisms. Furthermore, the presence of linked mutations has two opposing effects on Haldane's sieve. First, recessive mutants are disproportionally likely to be lost in outcrossers, so it is likelier that dominant mutations will fix. Second, with elevated rates of adaptive mutation, selective interference annuls the advantage in selfing organisms of not suffering from Haldane's sieve; outcrossing organisms are more able to fix weak beneficial mutations of any dominance value. Overall, weakened recombination effects can greatly limit adaptation in selfing organisms.
Collapse
|
44
|
Steige KA, Slotte T. Genomic legacies of the progenitors and the evolutionary consequences of allopolyploidy. CURRENT OPINION IN PLANT BIOLOGY 2016; 30:88-93. [PMID: 26943938 DOI: 10.1016/j.pbi.2016.02.006] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2015] [Revised: 02/13/2016] [Accepted: 02/15/2016] [Indexed: 06/05/2023]
Abstract
The formation of an allopolyploid species involves the merger of genomes with separate evolutionary histories and thereby different genomic legacies. Contrary to expectations from theory, genes from one are often lost preferentially in allopolyploids - there is biased fractionation. Here, we provide an overview of two ways in which the genomic legacies of the progenitors may impact the fate of duplicated genes in allopolyploids. Specifically, we discuss the role of homeolog expression biases in setting the stage for biased fractionation, and the evidence for transposable element silencing as a possible mechanism for homeolog expression biases. Finally, we highlight how differences between the progenitors with respect to accumulation of deleterious variation may affect trajectories of duplicate gene evolution in allopolyploids.
Collapse
Affiliation(s)
- Kim A Steige
- Department of Ecology and Genetics, Uppsala University, Sweden; Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Sweden
| | - Tanja Slotte
- Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Sweden.
| |
Collapse
|
45
|
Settepani V, Bechsgaard J, Bilde T. Phylogenetic analysis suggests that sociality is associated with reduced effectiveness of selection. Ecol Evol 2016; 6:469-77. [PMID: 26843931 PMCID: PMC4729245 DOI: 10.1002/ece3.1886] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2015] [Accepted: 11/24/2015] [Indexed: 11/07/2022] Open
Abstract
The evolution of sociality in spiders is associated with female bias, reproductive skew and an inbreeding mating system, factors that cause a reduction in effective population size and increase effects of genetic drift. These factors act to decrease the effectiveness of selection, thereby increasing the fixation probability of deleterious mutations. Comparative studies of closely related species with contrasting social traits and mating systems provide the opportunity to test consequences of low effective population size on the effectiveness of selection empirically. We used phylogenetic analyses of three inbred social spider species and seven outcrossing subsocial species of the genus Stegodyphus, and compared dN/dS ratios and codon usage bias between social Inbreeding and subsocial outcrossing mating systems to assess the effectiveness of selection. The overall results do not differ significantly between the social inbreeding and outcrossing species, but suggest a tendency for lower codon usage bias and higher dN/dS ratios in the social inbreeding species compared with their outcrossing congeners. The differences in dN/dS ratio and codon usage bias between social and subsocial species are modest but consistent with theoretical expectations of reduced effectiveness of selection in species with relatively low effective population size. The modest differences are consistent with relatively recent evolution of social mating systems. Additionally, the short terminal branches and lack of speciation of the social lineages, together with low genetic diversity lend support for the transient state of permanent sociality in spiders.
Collapse
Affiliation(s)
- Virginia Settepani
- Department of BioscienceAarhus UniversityNy Munkegade 116, Building 15408000Aarhus CDenmark
| | - Jesper Bechsgaard
- Department of BioscienceAarhus UniversityNy Munkegade 116, Building 15408000Aarhus CDenmark
| | - Trine Bilde
- Department of BioscienceAarhus UniversityNy Munkegade 116, Building 15408000Aarhus CDenmark
| |
Collapse
|
46
|
Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species. Genetics 2015; 202:1185-200. [PMID: 26721855 PMCID: PMC4788117 DOI: 10.1534/genetics.115.183152] [Citation(s) in RCA: 65] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2015] [Accepted: 12/24/2015] [Indexed: 12/30/2022] Open
Abstract
A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.
Collapse
|
47
|
Association mapping reveals the role of purifying selection in the maintenance of genomic variation in gene expression. Proc Natl Acad Sci U S A 2015; 112:15390-5. [PMID: 26604315 DOI: 10.1073/pnas.1503027112] [Citation(s) in RCA: 64] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The evolutionary forces that maintain genetic variation in quantitative traits within populations remain poorly understood. One hypothesis suggests that variation is under purifying selection, resulting in an excess of low-frequency variants and a negative correlation between minor allele frequency and selection coefficients. Here, we test these predictions using the genetic loci associated with total expression variation (eQTLs) and allele-specific expression variation (aseQTLs) mapped within a single population of the plant Capsella grandiflora. In addition to finding eQTLs and aseQTLs for a large fraction of genes, we show that alleles at these loci are rarer than expected and exhibit a negative correlation between phenotypic effect size and frequency. Overall, our results show that the distribution of frequencies and effect sizes of the loci responsible for local expression variation within a single outcrossing population are consistent with the effects of purifying selection.
Collapse
|
48
|
Hartfield M. Evolutionary genetic consequences of facultative sex and outcrossing. J Evol Biol 2015; 29:5-22. [DOI: 10.1111/jeb.12770] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2015] [Revised: 09/24/2015] [Accepted: 09/28/2015] [Indexed: 01/16/2023]
Affiliation(s)
- M. Hartfield
- Department of Ecology and Evolutionary Biology; University of Toronto; Toronto ON Canada
- Bioinformatics Research Centre; University of Aarhus; Aarhus Denmark
| |
Collapse
|
49
|
Castellano D, Coronado-Zamora M, Campos JL, Barbadilla A, Eyre-Walker A. Adaptive Evolution Is Substantially Impeded by Hill-Robertson Interference in Drosophila. Mol Biol Evol 2015; 33:442-55. [PMID: 26494843 PMCID: PMC4794616 DOI: 10.1093/molbev/msv236] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Hill-Robertson interference (HRi) is expected to reduce the efficiency of natural selection when two or more linked selected sites do not segregate freely, but no attempt has been done so far to quantify the overall impact of HRi on the rate of adaptive evolution for any given genome. In this work, we estimate how much HRi impedes the rate of adaptive evolution in the coding genome of Drosophila melanogaster. We compiled a data set of 6,141 autosomal protein-coding genes from Drosophila, from which polymorphism levels in D. melanogaster and divergence out to D. yakuba were estimated. The rate of adaptive evolution was calculated using a derivative of the McDonald-Kreitman test that controls for slightly deleterious mutations. We find that the rate of adaptive amino acid substitution at a given position of the genome is positively correlated to both the rate of recombination and the mutation rate, and negatively correlated to the gene density of the region. These correlations are robust to controlling for each other, for synonymous codon bias and for gene functions related to immune response and testes. We show that HRi diminishes the rate of adaptive evolution by approximately 27%. Interestingly, genes with low mutation rates embedded in gene poor regions lose approximately 17% of their adaptive substitutions whereas genes with high mutation rates embedded in gene rich regions lose approximately 60%. We conclude that HRi hampers the rate of adaptive evolution in Drosophila and that the variation in recombination, mutation, and gene density along the genome affects the HRi effect.
Collapse
Affiliation(s)
- David Castellano
- Genomics, Bioinformatics and Evolution Group, Institut de Biotecnologia i de Biomedicina (IBB) and Department de Genètica i Microbiologia, Campus Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Spain
| | - Marta Coronado-Zamora
- Genomics, Bioinformatics and Evolution Group, Institut de Biotecnologia i de Biomedicina (IBB) and Department de Genètica i Microbiologia, Campus Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Spain
| | - Jose L Campos
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Antonio Barbadilla
- Genomics, Bioinformatics and Evolution Group, Institut de Biotecnologia i de Biomedicina (IBB) and Department de Genètica i Microbiologia, Campus Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Spain
| | - Adam Eyre-Walker
- Centre for the Study of Evolution, School of Life Sciences, University of Sussex, Brighton, United Kingdom
| |
Collapse
|
50
|
Steige KA, Reimegård J, Koenig D, Scofield DG, Slotte T. Cis-Regulatory Changes Associated with a Recent Mating System Shift and Floral Adaptation in Capsella. Mol Biol Evol 2015; 32:2501-14. [PMID: 26318184 PMCID: PMC4576713 DOI: 10.1093/molbev/msv169] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 ka. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE; however, only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.
Collapse
Affiliation(s)
- Kim A Steige
- Department of Ecology and Genetics, Uppsala University, Uppsala, Sweden
| | - Johan Reimegård
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Daniel Koenig
- Department of Molecular Biology, Max Planck Institute for Developmental Biology, Tübingen, Germany
| | - Douglas G Scofield
- Department of Ecology and Genetics, Uppsala University, Uppsala, Sweden Uppsala Multidisciplinary Center for Advanced Computational Science, Uppsala University, Uppsala Sweden
| | - Tanja Slotte
- Department of Ecology and Genetics, Uppsala University, Uppsala, Sweden Department of Ecology, Environment and Plant Sciences, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| |
Collapse
|