1
|
Münzbergová Z, Šurinová M, Biscarini F, Níčová E. Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity. J Evol Biol 2024; 37:704-716. [PMID: 38761114 DOI: 10.1093/jeb/voae060] [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/29/2023] [Revised: 04/15/2024] [Accepted: 05/17/2024] [Indexed: 05/20/2024]
Abstract
The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.
Collapse
Affiliation(s)
- Zuzana Münzbergová
- Department of Botany, Faculty of Science, Charles University, Benátská 2, Prague, Czech Republic
- Department of Population Ecology, Institute of Botany, Czech Academy of Sciences, Zámek 1, Průhonice, Czech Republic
| | - Maria Šurinová
- Department of Botany, Faculty of Science, Charles University, Benátská 2, Prague, Czech Republic
- Department of Population Ecology, Institute of Botany, Czech Academy of Sciences, Zámek 1, Průhonice, Czech Republic
| | - Filippo Biscarini
- Institute of Agricultural Biology and Biotechnology, National Research Council (IBBA-CNR), Milan, Italy
| | - Eva Níčová
- Department of Population Ecology, Institute of Botany, Czech Academy of Sciences, Zámek 1, Průhonice, Czech Republic
| |
Collapse
|
2
|
Amin MR, Hasan M, Arnab SP, DeGiorgio M. Tensor Decomposition-based Feature Extraction and Classification to Detect Natural Selection from Genomic Data. Mol Biol Evol 2023; 40:msad216. [PMID: 37772983 PMCID: PMC10581699 DOI: 10.1093/molbev/msad216] [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: 03/02/2023] [Revised: 08/10/2023] [Accepted: 09/14/2023] [Indexed: 09/30/2023] Open
Abstract
Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under nonconvex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data although preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx, which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.
Collapse
Affiliation(s)
- Md Ruhul Amin
- Department of Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, FL 33431, USA
| | - Mahmudul Hasan
- Department of Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, FL 33431, USA
| | - Sandipan Paul Arnab
- Department of Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, FL 33431, USA
| | - Michael DeGiorgio
- Department of Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, FL 33431, USA
| |
Collapse
|
3
|
Tavares H, Readshaw A, Kania U, de Jong M, Pasam RK, McCulloch H, Ward S, Shenhav L, Forsyth E, Leyser O. Artificial selection reveals complex genetic architecture of shoot branching and its response to nitrate supply in Arabidopsis. PLoS Genet 2023; 19:e1010863. [PMID: 37616321 PMCID: PMC10482290 DOI: 10.1371/journal.pgen.1010863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2022] [Revised: 09/06/2023] [Accepted: 07/08/2023] [Indexed: 08/26/2023] Open
Abstract
Quantitative traits may be controlled by many loci, many alleles at each locus, and subject to genotype-by-environment interactions, making them difficult to map. One example of such a complex trait is shoot branching in the model plant Arabidopsis, and its plasticity in response to nitrate. Here, we use artificial selection under contrasting nitrate supplies to dissect the genetic architecture of this complex trait, where loci identified by association mapping failed to explain heritability estimates. We found a consistent response to selection for high branching, with correlated responses in other traits such as plasticity and flowering time. Genome-wide scans for selection and simulations suggest that at least tens of loci control this trait, with a distinct genetic architecture between low and high nitrate treatments. While signals of selection could be detected in the populations selected for high branching on low nitrate, there was very little overlap in the regions selected in three independent populations. Thus the regulatory network controlling shoot branching can be tuned in different ways to give similar phenotypes.
Collapse
Affiliation(s)
- Hugo Tavares
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
| | - Anne Readshaw
- Department of Biology, University of York, York, United Kingdom
| | - Urszula Kania
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
| | - Maaike de Jong
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
| | - Raj K. Pasam
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
| | - Hayley McCulloch
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
| | - Sally Ward
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
- Department of Biology, University of York, York, United Kingdom
| | - Liron Shenhav
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
| | - Elizabeth Forsyth
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
| | - Ottoline Leyser
- Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom
- Department of Biology, University of York, York, United Kingdom
| |
Collapse
|
4
|
Busch JW, Bodbyl‐Roels S, Tusuubira S, Kelly JK. Pollinator loss causes rapid adaptive evolution of selfing and dramatically reduces genome-wide genetic variability. Evolution 2022; 76:2130-2144. [PMID: 35852008 PMCID: PMC9543508 DOI: 10.1111/evo.14572] [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: 11/12/2021] [Revised: 03/23/2022] [Accepted: 04/20/2022] [Indexed: 01/22/2023]
Abstract
Although selfing populations harbor little genetic variation limiting evolutionary potential, the causes are unclear. We experimentally evolved large, replicate populations of Mimulus guttatus for nine generations in greenhouses with or without pollinating bees and studied DNA polymorphism in descendants. Populations without bees adapted to produce more selfed seed yet exhibited striking reductions in DNA polymorphism despite large population sizes. Importantly, the genome-wide pattern of variation cannot be explained by a simple reduction in effective population size, but instead reflects the complicated interaction between selection, linkage, and inbreeding. Simulations demonstrate that the spread of favored alleles at few loci depresses neutral variation genome wide in large populations containing fully selfing lineages. It also generates greater heterogeneity among chromosomes than expected with neutral evolution in small populations. Genome-wide deviations from neutrality were documented in populations with bees, suggesting widespread influences of background selection. After applying outlier tests to detect loci under selection, two genome regions were found in populations with bees, yet no adaptive loci were otherwise mapped. Large amounts of stochastic change in selfing populations compromise evolutionary potential and undermine outlier tests for selection. This occurs because genetic draft in highly selfing populations makes even the largest changes in allele frequency unremarkable.
Collapse
Affiliation(s)
- Jeremiah W. Busch
- School of Biological SciencesWashington State UniversityPullmanWashington99164
| | - Sarah Bodbyl‐Roels
- Trefny Innovative Instruction CenterColorado School of MinesGoldenColorado80401
| | - Sharif Tusuubira
- Department of Ecology and Evolutionary BiologyUniversity of KansasLawrenceKansas66045
| | - John K. Kelly
- Department of Ecology and Evolutionary BiologyUniversity of KansasLawrenceKansas66045
| |
Collapse
|
5
|
Genome-wide signatures of synergistic epistasis during parallel adaptation in a Baltic Sea copepod. Nat Commun 2022; 13:4024. [PMID: 35821220 PMCID: PMC9276764 DOI: 10.1038/s41467-022-31622-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2021] [Accepted: 06/27/2022] [Indexed: 01/01/2023] Open
Abstract
The role of epistasis in driving adaptation has remained an unresolved problem dating back to the Evolutionary Synthesis. In particular, whether epistatic interactions among genes could promote parallel evolution remains unexplored. To address this problem, we employ an Evolve and Resequence (E&R) experiment, using the copepod Eurytemora affinis, to elucidate the evolutionary genomic response to rapid salinity decline. Rapid declines in coastal salinity at high latitudes are a predicted consequence of global climate change. Based on time-resolved pooled whole-genome sequencing, we uncover a remarkably parallel, polygenic response across ten replicate selection lines, with 79.4% of selected alleles shared between lines by the tenth generation of natural selection. Using extensive computer simulations of our experiment conditions, we find that this polygenic parallelism is consistent with positive synergistic epistasis among alleles, far more so than other mechanisms tested. Our study provides experimental and theoretical support for a novel mechanism promoting repeatable polygenic adaptation, a phenomenon that may be common for selection on complex physiological traits. Using time-series whole-genome sequencing data from a laboratory evolution experiment, along with extensive computer simulations, the authors show that synergistic epistasis could drive rapid parallel freshwater adaptation in a saline copepod.
Collapse
|
6
|
Scossa F, Fernie AR. The evolution of metabolism: How to test evolutionary hypotheses at the genomic level. Comput Struct Biotechnol J 2020; 18:482-500. [PMID: 32180906 PMCID: PMC7063335 DOI: 10.1016/j.csbj.2020.02.009] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2019] [Revised: 02/12/2020] [Accepted: 02/13/2020] [Indexed: 01/21/2023] Open
Abstract
The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.
Collapse
Affiliation(s)
- Federico Scossa
- Max-Planck-Institut für Molekulare Pflanzenphysiologie, 14476 Potsdam-Golm, Germany
- Council for Agricultural Research and Economics (CREA), Research Centre for Genomics and Bioinformatics (CREA-GB), Via Ardeatina 546, 00178 Rome, Italy
| | - Alisdair R. Fernie
- Max-Planck-Institut für Molekulare Pflanzenphysiologie, 14476 Potsdam-Golm, Germany
- Center of Plant Systems Biology and Biotechnology (CPSBB), Plovdiv, Bulgaria
| |
Collapse
|
7
|
Fishman L, McIntosh M. Standard Deviations: The Biological Bases of Transmission Ratio Distortion. Annu Rev Genet 2019; 53:347-372. [DOI: 10.1146/annurev-genet-112618-043905] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The rule of Mendelian inheritance is remarkably robust, but deviations from the equal transmission of alternative alleles at a locus [a.k.a. transmission ratio distortion (TRD)] are also commonly observed in genetic mapping populations. Such TRD reveals locus-specific selection acting at some point between the diploid heterozygous parents and progeny genotyping and therefore can provide novel insight into otherwise-hidden genetic and evolutionary processes. Most of the classic selfish genetic elements were discovered through their biasing of transmission, but many unselfish evolutionary and developmental processes can also generate TRD. In this review, we describe methodologies for detecting TRD in mapping populations, detail the arenas and genetic interactions that shape TRD during plant and animal reproduction, and summarize patterns of TRD from across the genetic mapping literature. Finally, we point to new experimental approaches that can accelerate both detection of TRD and characterization of the underlying genetic mechanisms.
Collapse
Affiliation(s)
- Lila Fishman
- Division of Biological Sciences, University of Montana, Missoula, Montana 59812, USA
| | - Mariah McIntosh
- Division of Biological Sciences, University of Montana, Missoula, Montana 59812, USA
| |
Collapse
|
8
|
Brown KE, Kelly JK. Severe inbreeding depression is predicted by the “rare allele load” in
Mimulus guttatus
*. Evolution 2019; 74:587-596. [DOI: 10.1111/evo.13876] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Revised: 09/30/2019] [Accepted: 10/09/2019] [Indexed: 01/04/2023]
Affiliation(s)
- Keely E. Brown
- Ecology and Evolutionary Biology University of Kansas Lawrence Kansas 66045
| | - John K. Kelly
- Ecology and Evolutionary Biology University of Kansas Lawrence Kansas 66045
| |
Collapse
|
9
|
Hoedjes KM, van den Heuvel J, Kapun M, Keller L, Flatt T, Zwaan BJ. Distinct genomic signals of lifespan and life history evolution in response to postponed reproduction and larval diet in Drosophila. Evol Lett 2019; 3:598-609. [PMID: 31867121 PMCID: PMC6906992 DOI: 10.1002/evl3.143] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2019] [Revised: 09/08/2019] [Accepted: 09/11/2019] [Indexed: 12/14/2022] Open
Abstract
Reproduction and diet are two major factors controlling the physiology of aging and life history, but how they interact to affect the evolution of longevity is unknown. Moreover, although studies of large-effect mutants suggest an important role of nutrient sensing pathways in regulating aging, the genetic basis of evolutionary changes in lifespan remains poorly understood. To address these questions, we analyzed the genomes of experimentally evolved Drosophila melanogaster populations subjected to a factorial combination of two selection regimes: reproductive age (early versus postponed), and diet during the larval stage ("low," "control," "high"), resulting in six treatment combinations with four replicate populations each. Selection on reproductive age consistently affected lifespan, with flies from the postponed reproduction regime having evolved a longer lifespan. In contrast, larval diet affected lifespan only in early-reproducing populations: flies adapted to the "low" diet lived longer than those adapted to control diet. Here, we find genomic evidence for strong independent evolutionary responses to either selection regime, as well as loci that diverged in response to both regimes, thus representing genomic interactions between the two. Overall, we find that the genomic basis of longevity is largely independent of dietary adaptation. Differentiated loci were not enriched for "canonical" longevity genes, suggesting that naturally occurring genic targets of selection for longevity differ qualitatively from variants found in mutant screens. Comparing our candidate loci to those from other "evolve and resequence" studies of longevity demonstrated significant overlap among independent experiments. This suggests that the evolution of longevity, despite its presumed complex and polygenic nature, might be to some extent convergent and predictable.
Collapse
Affiliation(s)
- Katja M Hoedjes
- Department of Ecology and Evolution University of Lausanne Lausanne Switzerland
| | - Joost van den Heuvel
- Laboratory of Genetics, Plant Sciences Group Wageningen University Wageningen The Netherlands.,Institute for Cell and Molecular Biosciences Newcastle University Newcastle Upon Tyne United Kingdom
| | - Martin Kapun
- Department of Ecology and Evolution University of Lausanne Lausanne Switzerland.,Department of Biology University of Fribourg Fribourg Switzerland.,Current Address: Department of Evolutionary Biology and Environmental Studies University of Zurich Zurich Switzerland
| | - Laurent Keller
- Department of Ecology and Evolution University of Lausanne Lausanne Switzerland
| | - Thomas Flatt
- Department of Ecology and Evolution University of Lausanne Lausanne Switzerland.,Department of Biology University of Fribourg Fribourg Switzerland
| | - Bas J Zwaan
- Laboratory of Genetics, Plant Sciences Group Wageningen University Wageningen The Netherlands
| |
Collapse
|
10
|
Vlachos C, Burny C, Pelizzola M, Borges R, Futschik A, Kofler R, Schlötterer C. Benchmarking software tools for detecting and quantifying selection in evolve and resequencing studies. Genome Biol 2019; 20:169. [PMID: 31416462 PMCID: PMC6694636 DOI: 10.1186/s13059-019-1770-8] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2019] [Accepted: 07/22/2019] [Indexed: 12/15/2022] Open
Abstract
BACKGROUND The combination of experimental evolution with whole-genome resequencing of pooled individuals, also called evolve and resequence (E&R) is a powerful approach to study the selection processes and to infer the architecture of adaptive variation. Given the large potential of this method, a range of software tools were developed to identify selected SNPs and to measure their selection coefficients. RESULTS In this benchmarking study, we compare 15 test statistics implemented in 10 software tools using three different scenarios. We demonstrate that the power of the methods differs among the scenarios, but some consistently outperform others. LRT-1, CLEAR, and the CMH test perform best despite LRT-1 and the CMH test not requiring time series data. CLEAR provides the most accurate estimates of selection coefficients. CONCLUSION This benchmark study will not only facilitate the analysis of already existing data, but also affect the design of future data collections.
Collapse
Affiliation(s)
- Christos Vlachos
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Wien, 1210, Austria
- Vienna Graduate School of Population Genetics, Vienna, Austria
| | - Claire Burny
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Wien, 1210, Austria
- Vienna Graduate School of Population Genetics, Vienna, Austria
| | - Marta Pelizzola
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Wien, 1210, Austria
- Vienna Graduate School of Population Genetics, Vienna, Austria
| | - Rui Borges
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Wien, 1210, Austria
| | - Andreas Futschik
- Institute of Applied Statistics, Johannes Kepler University, Linz, 4040, Austria
- Plattform Bioinformatik und Biostatistik, Vetmeduni Vienna, Veterinärplatz 1, Wien, 1210, Austria
| | - Robert Kofler
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Wien, 1210, Austria.
| | - Christian Schlötterer
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Wien, 1210, Austria.
| |
Collapse
|
11
|
Flagel LE, Blackman BK, Fishman L, Monnahan PJ, Sweigart A, Kelly JK. GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity. PLoS Comput Biol 2019; 15:e1006949. [PMID: 30986215 PMCID: PMC6483263 DOI: 10.1371/journal.pcbi.1006949] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2018] [Revised: 04/25/2019] [Accepted: 03/15/2019] [Indexed: 11/18/2022] Open
Abstract
Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.
Collapse
Affiliation(s)
- Lex E. Flagel
- Bayer Crop Science, Chesterfield, MO, United States of America
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, MN, United States of America
- * E-mail: (LEF); (JKK)
| | - Benjamin K. Blackman
- Department of Plant and Microbial Biology, University of California—Berkeley, Berkeley, CA, United States of America
| | - Lila Fishman
- Division of Biological Sciences, University of Montana, Missoula, MT, United States of America
| | - Patrick J. Monnahan
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, United States of America
- Department of Ecology, Evolution, and Behavior, University of Minnesota, St. Paul, MN, United States of America
| | - Andrea Sweigart
- Department of Genetics, University of Georgia, Athens, GA, United States of America
| | - John K. Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, United States of America
- * E-mail: (LEF); (JKK)
| |
Collapse
|
12
|
Kelly JK, Hughes KA. Pervasive Linked Selection and Intermediate-Frequency Alleles Are Implicated in an Evolve-and-Resequencing Experiment of Drosophila simulans. Genetics 2019; 211:943-961. [PMID: 30593495 PMCID: PMC6404262 DOI: 10.1534/genetics.118.301824] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2018] [Accepted: 12/15/2018] [Indexed: 11/18/2022] Open
Abstract
We develop analytical and simulation tools for evolve-and-resequencing experiments and apply them to a new study of rapid evolution in Drosophila simulans Likelihood test statistics applied to pooled population sequencing data suggest parallel evolution of 138 SNPs across the genome. This number is reduced by orders of magnitude from previous studies (thousands or tens of thousands), owing to differences in both experimental design and statistical analysis. Whole genome simulations calibrated from Drosophila genetic data sets indicate that major features of the genome-wide response could be explained by as few as 30 loci under strong directional selection with a corresponding hitchhiking effect. Smaller effect loci are likely also responding, but are below the detection limit of the experiment. Finally, SNPs showing strong parallel evolution in the experiment are intermediate in frequency in the natural population (usually 30-70%) indicative of balancing selection in nature. These loci also exhibit elevated differentiation among natural populations of D. simulans, suggesting environmental heterogeneity as a potential balancing mechanism.
Collapse
Affiliation(s)
- John K Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas 66045
| | - Kimberly A Hughes
- Department of Biological Science, Florida State University, Tallahassee, Florida 32306
| |
Collapse
|
13
|
Nelson TC, Monnahan PJ, McIntosh MK, Anderson K, MacArthur-Waltz E, Finseth FR, Kelly JK, Fishman L. Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers. Mol Ecol 2018; 28:1460-1475. [PMID: 30346101 DOI: 10.1111/mec.14904] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2018] [Revised: 10/03/2018] [Accepted: 10/08/2018] [Indexed: 12/15/2022]
Abstract
Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness-relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate-frequency three-copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over-represented in late-flowering and large-flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250-300 copies of RLG1a totalling ~5.7 Mb (20-40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR-based copy counts indicate that high+ is a single tandem array unlinked to the single-copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single-copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress-responsive life-history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.
Collapse
Affiliation(s)
- Thomas C Nelson
- Division of Biological Sciences, University of Montana, Missoula, Montana
| | - Patrick J Monnahan
- Department of Ecology and Evolution, University of Kansas, Lawrence, Kansas
| | - Mariah K McIntosh
- Division of Biological Sciences, University of Montana, Missoula, Montana
| | - Kayli Anderson
- Division of Biological Sciences, University of Montana, Missoula, Montana
| | | | - Findley R Finseth
- Division of Biological Sciences, University of Montana, Missoula, Montana
| | - John K Kelly
- Department of Ecology and Evolution, University of Kansas, Lawrence, Kansas
| | - Lila Fishman
- Division of Biological Sciences, University of Montana, Missoula, Montana
| |
Collapse
|
14
|
Signor SA, New FN, Nuzhdin S. A Large Panel of Drosophila simulans Reveals an Abundance of Common Variants. Genome Biol Evol 2018; 10:189-206. [PMID: 29228179 PMCID: PMC5767965 DOI: 10.1093/gbe/evx262] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/07/2017] [Indexed: 01/03/2023] Open
Abstract
The rapidly expanding availability of large NGS data sets provides an opportunity to investigate population genetics at an unprecedented scale. Drosophila simulans is the sister species of the model organism Drosophila melanogaster, and is often presumed to share similar demographic history. However, previous population genetic and ecological work suggests very different signatures of selection and demography. Here, we sequence a new panel of 170 inbred genotypes of a North American population of D. simulans, a valuable complement to the DGRP and other D. melanogaster panels. We find some unexpected signatures of demography, in the form of excess intermediate frequency polymorphisms. Simulations suggest that this is possibly due to a recent population contraction and selection. We examine the outliers in the D. simulans genome determined by a haplotype test to attempt to parse the contribution of demography and selection to the patterns observed in this population. Untangling the relative contribution of demography and selection to genomic patterns of variation is challenging, however, it is clear that although D. melanogaster was thought to share demographic history with D. simulans different forces are at work in shaping genomic variation in this population of D. simulans.
Collapse
Affiliation(s)
- Sarah A Signor
- Department of Molecular and Computational Biology, University of Southern California
| | - Felicia N New
- Department of Molecular Genetics and Microbiology, University of Florida College of Medicine
| | - Sergey Nuzhdin
- Department of Molecular and Computational Biology, University of Southern California
| |
Collapse
|
15
|
Troth A, Puzey JR, Kim RS, Willis JH, Kelly JK. Selective trade-offs maintain alleles underpinning complex trait variation in plants. SCIENCE (NEW YORK, N.Y.) 2018; 361:475-478. [PMID: 30072534 DOI: 10.1126/science.aat5760] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Received: 03/13/2018] [Accepted: 06/22/2018] [Indexed: 11/02/2022]
Abstract
To understand evolutionary factors that maintain complex trait variation, we sequenced genomes from a single population of the plant Mimulus guttatus, identifying hundreds of nucleotide variants associated with morphological and life history traits. Alleles that delayed flowering also increased size at reproduction, which suggests pervasive antagonistic pleiotropy in this annual plant. The "large and slow" alleles, which were less common in small, rapidly flowering populations, became more abundant in populations with greater plant size. Furthermore, natural selection within the field population favored alternative alleles from year to year. Our results suggest that environmental fluctuations and selective trade-offs maintain polygenic trait variation within populations and also contribute to the geographic divergence in this wildflower species.
Collapse
Affiliation(s)
- Ashley Troth
- Department of Biology, Duke University, Durham, NC 27708, USA
| | - Joshua R Puzey
- Department of Biology, Duke University, Durham, NC 27708, USA.,Department of Biology, College of William and Mary, Williamsburg, VA 23187, USA
| | - Rebecca S Kim
- Department of Biology, Duke University, Durham, NC 27708, USA.,Department of Environmental Medicine, NYU Langone, New York, NY 10016, USA
| | - John H Willis
- Department of Biology, Duke University, Durham, NC 27708, USA
| | - John K Kelly
- Department of Ecology and Evolution, University of Kansas, Lawrence, KS 27708, USA.
| |
Collapse
|
16
|
Zuellig MP, Sweigart AL. Gene duplicates cause hybrid lethality between sympatric species of Mimulus. PLoS Genet 2018; 14:e1007130. [PMID: 29649209 PMCID: PMC5896889 DOI: 10.1371/journal.pgen.1007130] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2017] [Accepted: 11/27/2017] [Indexed: 01/30/2023] Open
Abstract
Hybrid incompatibilities play a critical role in the evolution and maintenance of species. We have discovered a simple genetic incompatibility that causes lethality in hybrids between two closely related species of yellow monkeyflower (Mimulus guttatus and M. nasutus). This hybrid incompatibility, which causes one sixteenth of F2 hybrid seedlings to lack chlorophyll and die shortly after germination, occurs between sympatric populations that are connected by ongoing interspecific gene flow. Using complimentary genetic mapping and gene expression analyses, we show that lethality occurs in hybrids that lack a functional copy of the critical photosynthetic gene pTAC14. In M. guttatus, this gene was duplicated, but the ancestral copy is no longer expressed. In M. nasutus, the duplication is missing altogether. As a result, hybrids die when they are homozygous for the nonfunctional M. guttatus copy and missing the duplicate from M. nasutus, apparently due to misregulated transcription of key photosynthetic genes. Our study indicates that neutral evolutionary processes may play an important role in the evolution of hybrid incompatibilities and opens the door to direct investigations of their contribution to reproductive isolation among naturally hybridizing species. The evolution of hybrid incompatibilities (gene interactions that cause hybrids to be sterile or inviable) is a common outcome of genomic divergence between lineages. However, evaluating the importance of hybrid incompatibilities for speciation requires that we identify the causal genes and evolutionary forces in recently diverged, wild species. We discovered that hybrid seedling lethality between two closely related sister species of yellow monkeyflower is caused by duplicate copies of a gene critical for chloroplast development. Because each lineage carries its one functional gene copy in a distinct genomic location, some hybrids inherit only inactive (or missing) alleles. We infer that hybrid lethality in this young species pair has arisen through divergent resolution of gene duplicates by degenerative mutations and (likely) genetic drift. These findings are an important step toward understanding the evolutionary dynamics of hybrid incompatibility genes in nature, as well as the role of such genes in species divergence.
Collapse
Affiliation(s)
- Matthew P. Zuellig
- Department of Genetics, University of Georgia, Athens, Georgia, United States of America
- * E-mail:
| | - Andrea L. Sweigart
- Department of Genetics, University of Georgia, Athens, Georgia, United States of America
| |
Collapse
|
17
|
Zhang L, Xu P, Cai Y, Ma L, Li S, Li S, Xie W, Song J, Peng L, Yan H, Zou L, Ma Y, Zhang C, Gao Q, Wang J. The draft genome assembly of Rhododendron delavayi Franch. var. delavayi. Gigascience 2017; 6:1-11. [PMID: 29020749 PMCID: PMC5632301 DOI: 10.1093/gigascience/gix076] [Citation(s) in RCA: 49] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Revised: 05/28/2017] [Accepted: 08/04/2017] [Indexed: 01/16/2023] Open
Abstract
Rhododendron delavayi Franch. is globally famous as an ornamental plant. Its distribution in southwest China covers several different habitats and environments. However, not much research had been conducted on Rhododendron spp. at the molecular level, which hinders understanding of its evolution, speciation, and synthesis of secondary metabolites, as well as its wide adaptability to different environments. Here, we report the genome assembly and gene annotation of R. delavayi var. delavayi (the second genome sequenced in the Ericaceae), which will facilitate the study of the family. The genome assembly will have further applications in genome-assisted cultivar breeding. The final size of the assembled R. delavayi var. delavayi genome (695.09 Mb) was close to the 697.94 Mb, estimated by k-mer analysis. A total of 336.83 gigabases (Gb) of raw Illumina HiSeq 2000 reads were generated from 9 libraries (with insert sizes ranging from 170 bp to 40 kb), achieving a raw sequencing depth of ×482.6. After quality filtering, 246.06 Gb of clean reads were obtained, giving ×352.55 coverage depth. Assembly using Platanus gave a total scaffold length of 695.09 Mb, with a contig N50 of 61.8 kb and a scaffold N50 of 637.83 kb. Gene prediction resulted in the annotation of 32 938 protein-coding genes. The genome completeness was evaluated by CEGMA and BUSCO and reached 95.97% and 92.8%, respectively. The gene annotation completeness was also evaluated by CEGMA and BUSCO and reached 97.01% and 87.4%, respectively. Genome annotation revealed that 51.77% of the R. delavayi genome is composed of transposable elements, and 37.48% of long terminal repeat elements (LTRs). The de novo assembled genome of R. delavayi var. delavayi (hereinafter referred to as R. delavayi) is the second genomic resource of the family Ericaceae and will provide a valuable resource for research on future comparative genomic studies in Rhododendron species. The availability of the R. delavayi genome sequence will hopefully provide a tool for scientists to tackle open questions regarding molecular mechanisms underlying environmental interactions in the genus Rhododendron, more accurately understand the evolutionary processes and systematics of the genus, facilitate the identification of genes encoding pharmaceutically important compounds, and accelerate molecular breeding to release elite varieties.
Collapse
Affiliation(s)
- Lu Zhang
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Pengwei Xu
- BGI-Shenzhen, BGI Park, No. 21 Hongan 3rd Street, Yantian District, Shenzhen 518083, China
| | - Yanfei Cai
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Lulin Ma
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Shifeng Li
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Shufa Li
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Weijia Xie
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Jie Song
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Lvchun Peng
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Huijun Yan
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Ling Zou
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| | - Yongpeng Ma
- Kunming Botanical Garden, Kunming Institute of Botany, Chinese Academy of Sciences, No. 132 Lanhei Road, Panlong District, Kunming, Yunnan 650201, China
| | - Chengjun Zhang
- Germplasm Bank of Wild species, Kunming Institute of Botany, Chinese Academy of Sciences, No. 132 Lanhei Road, Panlong District, Kunming, Yunnan 650201, China
| | - Qiang Gao
- BGI-Shenzhen, BGI Park, No. 21 Hongan 3rd Street, Yantian District, Shenzhen 518083, China
| | - Jihua Wang
- Flower Research Institute of Yunnan Academy of Agricultural Sciences, National Engineering Research Center For Ornamental Horticulture, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
- Key Lab of Yunnan Flower Breeding, No. 2238 Beijing Road, Panlong District, Kunming 650205, China
| |
Collapse
|
18
|
Monnahan PJ, Kelly JK. The Genomic Architecture of Flowering Time Varies Across Space and Time in Mimulus guttatus. Genetics 2017; 206:1621-1635. [PMID: 28455350 PMCID: PMC5500155 DOI: 10.1534/genetics.117.201483] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2017] [Accepted: 04/23/2017] [Indexed: 11/18/2022] Open
Abstract
The degree to which genomic architecture varies across space and time is central to the evolution of genomes in response to natural selection. Bulked-segregant mapping combined with pooled sequencing provides an efficient means to estimate the effect of genetic variants on quantitative traits. We develop a novel likelihood framework to identify segregating variation within multiple populations and generations while accommodating estimation error on a sample- and SNP-specific basis. We use this method to map loci for flowering time within natural populations of Mimulus guttatus, collecting the early- and late-flowering plants from each of three neighboring populations and two consecutive generations. Structural variants, such as inversions, and genes from multiple flowering-time pathways exhibit the strongest associations with flowering time. We find appreciable variation in genetic effects on flowering time across both time and space; the greatest differences evident between populations, where numerous factors (environmental variation, genomic background, and private polymorphisms) likely contribute to heterogeneity. However, the changes across years within populations clearly identify genotype-by-environment interactions as an important influence on flowering time variation.
Collapse
Affiliation(s)
- Patrick J Monnahan
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas 66045
| | - John K Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas 66045
| |
Collapse
|
19
|
Hoban S, Kelley JL, Lotterhos KE, Antolin MF, Bradburd G, Lowry DB, Poss ML, Reed LK, Storfer A, Whitlock MC. Finding the Genomic Basis of Local Adaptation: Pitfalls, Practical Solutions, and Future Directions. Am Nat 2016; 188:379-97. [PMID: 27622873 PMCID: PMC5457800 DOI: 10.1086/688018] [Citation(s) in RCA: 431] [Impact Index Per Article: 53.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Uncovering the genetic and evolutionary basis of local adaptation is a major focus of evolutionary biology. The recent development of cost-effective methods for obtaining high-quality genome-scale data makes it possible to identify some of the loci responsible for adaptive differences among populations. Two basic approaches for identifying putatively locally adaptive loci have been developed and are broadly used: one that identifies loci with unusually high genetic differentiation among populations (differentiation outlier methods) and one that searches for correlations between local population allele frequencies and local environments (genetic-environment association methods). Here, we review the promises and challenges of these genome scan methods, including correcting for the confounding influence of a species' demographic history, biases caused by missing aspects of the genome, matching scales of environmental data with population structure, and other statistical considerations. In each case, we make suggestions for best practices for maximizing the accuracy and efficiency of genome scans to detect the underlying genetic basis of local adaptation. With attention to their current limitations, genome scan methods can be an important tool in finding the genetic basis of adaptive evolutionary change.
Collapse
Affiliation(s)
- Sean Hoban
- Morton Arboretum, Lisle, Illinois 60532; and National Institute for Mathematical and Biological Synthesis (NIMBioS), Knoxville, Tennessee 37966
| | - Joanna L. Kelley
- School of Biological Sciences, Washington State University, Pullman, Washington 99164
| | - Katie E. Lotterhos
- Department of Marine and Environmental Sciences, Northeastern University Marine Science Center, Nahant, Massachusetts 01908
| | - Michael F. Antolin
- Department of Biology, Colorado State University, Fort Collins, Colorado 80523
| | - Gideon Bradburd
- Museum of Vertebrate Zoology and Department of Environmental Science, Policy, and Management, University of California, Berkeley, California 94720
| | - David B. Lowry
- Department of Plant Biology, Michigan State University, East Lansing, Michigan 48824
| | - Mary L. Poss
- Department of Biology and Veterinary and Biomedical Sciences, Penn State University, University Park, Pennsylvania 16802
| | - Laura K. Reed
- Department of Biological Sciences, University of Alabama, Tuscaloosa, Alabama 35406
| | - Andrew Storfer
- School of Biological Sciences, Washington State University, Pullman, Washington 99164
| | | |
Collapse
|
20
|
Lee YW, Fishman L, Kelly JK, Willis JH. A Segregating Inversion Generates Fitness Variation in Yellow Monkeyflower (Mimulus guttatus). Genetics 2016; 202:1473-84. [PMID: 26868767 PMCID: PMC4905537 DOI: 10.1534/genetics.115.183566] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2015] [Accepted: 02/01/2016] [Indexed: 11/18/2022] Open
Abstract
Polymorphic chromosomal rearrangements can bind hundreds of genes into single genetic loci with diverse effects. Rearrangements are often associated with local adaptation and speciation and may also be an important component of genetic variation within populations. We genetically and phenotypically characterize a segregating inversion (inv6) in the Iron Mountain (IM) population of Mimulus guttatus (yellow monkeyflower). We initially mapped inv6 as a region of recombination suppression in three F2 populations resulting from crosses among IM plants. In each case, the F1 parent was heterozygous for a derived haplotype, homogenous across markers spanning over 5 Mb of chromsome 6. In the three F2 populations, inv6 reduced male and female fitness components. In addition,i nv6 carriers suffered an ∼30% loss of pollen viability in the field. Despite these costs, inv6 exists at moderate frequency (∼8%) in the natural population, suggesting counterbalancing fitness benefits that maintain the polymorphism. Across 4 years of monitoring in the field, inv6 had an overall significant positive effect on seed production (lifetime female fitness) of carriers. This benefit was particularly strong in harsh years and may be mediated (in part) by strong positive effects on flower production. These data suggest that opposing fitness effects maintain an intermediate frequency, and as a consequence, inv6 generates inbreeding depression and high genetic variance. We discuss these findings in relation to the theory of inbreeding depression and the maintenance of fitness variation.
Collapse
Affiliation(s)
- Young Wha Lee
- Biology Department, Duke University, Durham, North Carolina 27708
| | - Lila Fishman
- Division of Biological Sciences, University of Montana, Missoula, Montana 59812
| | - John K Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas 66045
| | - John H Willis
- Biology Department, Duke University, Durham, North Carolina 27708
| |
Collapse
|
21
|
Mikheyev AS, Tin MMY, Arora J, Seeley TD. Museum samples reveal rapid evolution by wild honey bees exposed to a novel parasite. Nat Commun 2015; 6:7991. [PMID: 26246313 PMCID: PMC4918369 DOI: 10.1038/ncomms8991] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2015] [Accepted: 07/06/2015] [Indexed: 11/26/2022] Open
Abstract
Understanding genetic changes caused by novel pathogens and parasites can reveal mechanisms of adaptation and genetic robustness. Using whole-genome sequencing of museum and modern specimens, we describe the genomic changes in a wild population of honey bees in North America following the introduction of the ectoparasitic mite, Varroa destructor. Even though colony density in the study population is the same today as in the past, a major loss of haplotypic diversity occurred, indicative of a drastic mitochondrial bottleneck, caused by massive colony mortality. In contrast, nuclear genetic diversity did not change, though hundreds of genes show signs of selection. The genetic diversity within each bee colony, particularly as a consequence of polyandry by queens, may enable preservation of genetic diversity even during population bottlenecks. These findings suggest that genetically diverse honey bee populations can recover from introduced diseases by evolving rapid tolerance, while maintaining much of the standing genetic variation.
Collapse
Affiliation(s)
- Alexander S. Mikheyev
- Okinawa Institute of Science and Technology, 1919-1 Tancha, Onna-son, Kunigami-gun 904-0412, Japan
- Research School of Biology, Australian National University, Canberra, Australian Capital Territory 0200, Australia
| | - Mandy M. Y. Tin
- Okinawa Institute of Science and Technology, 1919-1 Tancha, Onna-son, Kunigami-gun 904-0412, Japan
| | - Jatin Arora
- Okinawa Institute of Science and Technology, 1919-1 Tancha, Onna-son, Kunigami-gun 904-0412, Japan
- Department of Paediatrics and Adolescent Medicine, Medical University of Vienna, Lazarettgasse 14, AKH BT 25.3, 1090 Vienna, Austria
| | - Thomas D. Seeley
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
22
|
Colicchio JM, Miura F, Kelly JK, Ito T, Hileman LC. DNA methylation and gene expression in Mimulus guttatus. BMC Genomics 2015; 16:507. [PMID: 26148779 PMCID: PMC4492170 DOI: 10.1186/s12864-015-1668-0] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2015] [Accepted: 05/29/2015] [Indexed: 11/13/2022] Open
Abstract
Background The presence of methyl groups on cytosine nucleotides across an organism’s genome (methylation) is a major regulator of genome stability, crossing over, and gene regulation. The capacity for DNA methylation to be altered by environmental conditions, and potentially passed between generations, makes it a prime candidate for transgenerational epigenetic inheritance. Here we conduct the first analysis of the Mimulus guttatus methylome, with a focus on the relationship between DNA methylation and gene expression. Results We present a whole genome methylome for the inbred line Iron Mountain 62 (IM62). DNA methylation varies across chromosomes, genomic regions, and genes. We develop a model that predicts gene expression based on DNA methylation (R2 = 0.2). Post hoc analysis of this model confirms prior relationships, and identifies novel relationships between methylation and gene expression. Additionally, we find that DNA methylation is significantly depleted near gene transcriptional start sites, which may explain the recently discovered elevated rate of recombination in these same regions. Conclusions The establishment here of a reference methylome will be a useful resource for the continued advancement of M. guttatus as a model system. Using a model-based approach, we demonstrate that methylation patterns are an important predictor of variation in gene expression. This model provides a novel approach for differential methylation analysis that generates distinct and testable hypotheses regarding gene expression. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1668-0) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Jack M Colicchio
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, 66045, USA.
| | - Fumihito Miura
- Department of Medical Biochemistry, Department of Biochemistry, Fukuoka 812-8581, Fukuoka 812-8582, Japan
| | - John K Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, 66045, USA
| | - Takashi Ito
- Department of Medical Biochemistry, Department of Biochemistry, Fukuoka 812-8581, Fukuoka 812-8582, Japan
| | - Lena C Hileman
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, 66045, USA
| |
Collapse
|
23
|
Monnahan PJ, Kelly JK. Epistasis Is a Major Determinant of the Additive Genetic Variance in Mimulus guttatus. PLoS Genet 2015; 11:e1005201. [PMID: 25946702 PMCID: PMC4422649 DOI: 10.1371/journal.pgen.1005201] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2015] [Accepted: 04/08/2015] [Indexed: 11/21/2022] Open
Abstract
The influence of genetic interactions (epistasis) on the genetic variance of quantitative traits is a major unresolved problem relevant to medical, agricultural, and evolutionary genetics. The additive genetic component is typically a high proportion of the total genetic variance in quantitative traits, despite that underlying genes must interact to determine phenotype. This study estimates direct and interaction effects for 11 pairs of Quantitative Trait Loci (QTLs) affecting floral traits within a single population of Mimulus guttatus. With estimates of all 9 genotypes for each QTL pair, we are able to map from QTL effects to variance components as a function of population allele frequencies, and thus predict changes in variance components as allele frequencies change. This mapping requires an analytical framework that properly accounts for bias introduced by estimation errors. We find that even with abundant interactions between QTLs, most of the genetic variance is likely to be additive. However, the strong dependency of allelic average effects on genetic background implies that epistasis is a major determinant of the additive genetic variance, and thus, the population’s ability to respond to selection. Complex traits are influenced not only by the effects of individual genes but also by the myriad ways that these genes interact with one another, commonly referred to as epistasis. Theory suggests that epistasis could have important population-level implications in terms of the genetic variance components that govern evolution in response to natural or artificial selection. Unfortunately, empirical examples extending from observed interactions between genes to genetic variances are scant, particularly for natural populations. Here, we characterize epistasis between naturally segregating polymorphisms in M. guttatus and determine the cumulative effect of epistasis on population genetic variance components. To do this, we first elaborate the necessary statistical theory to accommodate estimation error in genetic effects, as failing to do so will upwardly bias variance predictions. We find that gene interactions have a net positive effect on both the total and additive genetic variance for most traits; however, the contribution of individual loci to the additive variance depends heavily on the genotype frequencies at other loci. Therefore, the effect of epistasis extends beyond the individual’s phenotype to influence how both populations and their component alleles respond to selection.
Collapse
Affiliation(s)
- Patrick J. Monnahan
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Lawrence, Kansas, United States of America
- * E-mail:
| | - John K. Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Lawrence, Kansas, United States of America
| |
Collapse
|
24
|
Schlötterer C, Kofler R, Versace E, Tobler R, Franssen SU. Combining experimental evolution with next-generation sequencing: a powerful tool to study adaptation from standing genetic variation. Heredity (Edinb) 2015; 114:431-40. [PMID: 25269380 PMCID: PMC4815507 DOI: 10.1038/hdy.2014.86] [Citation(s) in RCA: 158] [Impact Index Per Article: 17.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2014] [Revised: 07/01/2014] [Accepted: 07/14/2014] [Indexed: 12/20/2022] Open
Abstract
Evolve and resequence (E&R) is a new approach to investigate the genomic responses to selection during experimental evolution. By using whole genome sequencing of pools of individuals (Pool-Seq), this method can identify selected variants in controlled and replicable experimental settings. Reviewing the current state of the field, we show that E&R can be powerful enough to identify causative genes and possibly even single-nucleotide polymorphisms. We also discuss how the experimental design and the complexity of the trait could result in a large number of false positive candidates. We suggest experimental and analytical strategies to maximize the power of E&R to uncover the genotype-phenotype link and serve as an important research tool for a broad range of evolutionary questions.
Collapse
Affiliation(s)
- C Schlötterer
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria
| | - R Kofler
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria
| | - E Versace
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria
- Center for Mind/Brain Sciences, University of Trento, Rovereto, Italy
| | - R Tobler
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria
- Vienna Graduate School of Population Genetics, Vienna, Austria
| | - S U Franssen
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria
| |
Collapse
|
25
|
Beissinger TM, Rosa GJM, Kaeppler SM, Gianola D, de Leon N. Defining window-boundaries for genomic analyses using smoothing spline techniques. Genet Sel Evol 2015; 47:30. [PMID: 25928167 PMCID: PMC4404117 DOI: 10.1186/s12711-015-0105-9] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2014] [Accepted: 02/04/2015] [Indexed: 01/29/2023] Open
Abstract
Background High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome. Results Simulations applying this method were performed to identify selection signatures from pooled sequencing FST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing FST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach. Conclusions We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on FST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.
Collapse
Affiliation(s)
| | - Guilherme J M Rosa
- Department of Animal Sciences, University of Wisconsin, Madison, 53706, USA. .,Department of Biostatistics and Medical Informatics, University of Wisconsin, Madison, 53792, USA.
| | - Shawn M Kaeppler
- Department of Agronomy, University of Wisconsin, Madison, 53706, USA. .,Department of Energy Great Lakes Bioenergy Research Center, University of Wisconsin, Madison, 53706, USA.
| | - Daniel Gianola
- Department of Animal Sciences, University of Wisconsin, Madison, 53706, USA. .,Department of Biostatistics and Medical Informatics, University of Wisconsin, Madison, 53792, USA. .,Department of Dairy Science, University of Wisconsin, Madison, 53706, USA.
| | - Natalia de Leon
- Department of Agronomy, University of Wisconsin, Madison, 53706, USA. .,Department of Energy Great Lakes Bioenergy Research Center, University of Wisconsin, Madison, 53706, USA.
| |
Collapse
|
26
|
Sequencing pools of individuals — mining genome-wide polymorphism data without big funding. Nat Rev Genet 2014; 15:749-63. [DOI: 10.1038/nrg3803] [Citation(s) in RCA: 512] [Impact Index Per Article: 51.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
27
|
Abstract
Local adaptation and adaptive clines are pervasive in natural plant populations, yet the effects of these types of adaptation on genomic diversity are not well understood. With a data set of 202 accessions of Medicago truncatula genotyped at almost 2 million single nucleotide polymorphisms, we used mixed linear models to identify candidate loci responsible for adaptation to three climatic gradients-annual mean temperature (AMT), precipitation in the wettest month (PWM), and isothermality (ITH)-representing the major axes of climate variation across the species' range. Loci with the strongest association to these climate gradients tagged genome regions with high sequence similarity to genes with functional roles in thermal tolerance, drought tolerance, or resistance to herbivores of pathogens. Genotypes at these candidate loci also predicted the performance of an independent sample of plant accessions grown in climate-controlled conditions. Compared to a genome-wide sample of randomly drawn reference SNPs, candidates for two climate gradients, AMT and PWM, were significantly enriched for genic regions, and genome segments flanking genic AMT and PWM candidates harbored less nucleotide diversity, elevated differentiation between haplotypes carrying alternate alleles, and an overrepresentation of the most common haplotypes. These patterns of diversity are consistent with a history of soft selective sweeps acting on loci underlying adaptation to climate, but not with a history of long-term balancing selection.
Collapse
|
28
|
A genome-wide scan for evidence of selection in a maize population under long-term artificial selection for ear number. Genetics 2013; 196:829-40. [PMID: 24381334 DOI: 10.1534/genetics.113.160655] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
A genome-wide scan to detect evidence of selection was conducted in the Golden Glow maize long-term selection population. The population had been subjected to selection for increased number of ears per plant for 30 generations, with an empirically estimated effective population size ranging from 384 to 667 individuals and an increase of more than threefold in the number of ears per plant. Allele frequencies at >1.2 million single-nucleotide polymorphism loci were estimated from pooled whole-genome resequencing data, and FST values across sliding windows were employed to assess divergence between the population preselection and the population postselection. Twenty-eight highly divergent regions were identified, with half of these regions providing gene-level resolution on potentially selected variants. Approximately 93% of the divergent regions do not demonstrate a significant decrease in heterozygosity, which suggests that they are not approaching fixation. Also, most regions display a pattern consistent with a soft-sweep model as opposed to a hard-sweep model, suggesting that selection mostly operated on standing genetic variation. For at least 25% of the regions, results suggest that selection operated on variants located outside of currently annotated coding regions. These results provide insights into the underlying genetic effects of long-term artificial selection and identification of putative genetic elements underlying number of ears per plant in maize.
Collapse
|