1
|
Anderson NW, Kirk L, Schraiber JG, Ragsdale AP. A Path Integral Approach for Allele Frequency Dynamics Under Polygenic Selection. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.14.599114. [PMID: 38915613 PMCID: PMC11195211 DOI: 10.1101/2024.06.14.599114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/26/2024]
Abstract
Many phenotypic traits have a polygenic genetic basis, making it challenging to learn their genetic architectures and predict individual phenotypes. One promising avenue to resolve the genetic basis of complex traits is through evolve-and-resequence experiments, in which laboratory populations are exposed to some selective pressure and trait-contributing loci are identified by extreme frequency changes over the course of the experiment. However, small laboratory populations will experience substantial random genetic drift, and it is difficult to determine whether selection played a roll in a given allele frequency change. Predicting how much allele frequencies change under drift and selection had remained an open problem well into the 21st century, even those contributing to simple, monogenic traits. Recently, there have been efforts to apply the path integral, a method borrowed from physics, to solve this problem. So far, this approach has been limited to genic selection, and is therefore inadequate to capture the complexity of quantitative, highly polygenic traits that are commonly studied. Here we extend one of these path integral methods, the perturbation approximation, to selection scenarios that are of interest to quantitative genetics. In particular, we derive analytic expressions for the transition probability (i.e., the probability that an allele will change in frequency from x , to y in time t ) of an allele contributing to a trait subject to stabilizing selection, as well as that of an allele contributing to a trait rapidly adapting to a new phenotypic optimum. We use these expressions to characterize the use of allele frequency change to test for selection, as well as explore optimal design choices for evolve-and-resequence experiments to uncover the genetic architecture of polygenic traits under selection.
Collapse
Affiliation(s)
- Nathan W. Anderson
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Lloyd Kirk
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Joshua G. Schraiber
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, 90089, USA
| | - Aaron P. Ragsdale
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, WI, 53706, USA
| |
Collapse
|
2
|
Fine AG, Steinrücken M. A novel expectation-maximization approach to infer general diploid selection from time-series genetic data. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.10.593575. [PMID: 38798346 PMCID: PMC11118272 DOI: 10.1101/2024.05.10.593575] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Detecting and quantifying the strength of selection is a main objective in population genetics. Since selection acts over multiple generations, many approaches have been developed to detect and quantify selection using genetic data sampled at multiple points in time. Such time series genetic data is commonly analyzed using Hidden Markov Models, but in most cases, under the assumption of additive selection. However, many examples of genetic variation exhibiting non-additive mechanisms exist, making it critical to develop methods that can characterize selection in more general scenarios. Thus, we extend a previously introduced expectation-maximization algorithm for the inference of additive selection coefficients to the case of general diploid selection, in which heterozygote and homozygote fitnesses are parameterized independently. We furthermore introduce a framework to identify bespoke modes of diploid selection from given data, as well as a procedure for aggregating data across linked loci to increase power and robustness. Using extensive simulation studies, we find that our method accurately and efficiently estimates selection coefficients for different modes of diploid selection across a wide range of scenarios; however, power to classify the mode of selection is low unless selection is very strong. We apply our method to ancient DNA samples from Great Britain in the last 4,450 years, and detect evidence for selection in six genomic regions, including the well-characterized LCT locus. Our work is the first genome-wide scan characterizing signals of general diploid selection.
Collapse
Affiliation(s)
- Adam G Fine
- Department of Ecology and Evolution, University of Chicago
- Graduate Program in Biophysical Sciences, University of Chicago
| | - Matthias Steinrücken
- Department of Ecology and Evolution, University of Chicago
- Department of Human Genetics, University of Chicago
| |
Collapse
|
3
|
Sabolić I, Mira Ó, Brandt DYC, Lisičić D, Stapley J, Novosolov M, Bakarić R, Cizelj I, Glogoški M, Hudina T, Taverne M, Allentoft ME, Nielsen R, Herrel A, Štambuk A. Plastic and genomic change of a newly established lizard population following a founder event. Mol Ecol 2024; 33:e17255. [PMID: 38133599 DOI: 10.1111/mec.17255] [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: 07/14/2022] [Revised: 12/06/2023] [Accepted: 12/13/2023] [Indexed: 12/23/2023]
Abstract
Understanding how phenotypic divergence arises among natural populations remains one of the major goals in evolutionary biology. As part of competitive exclusion experiment conducted in 1971, 10 individuals of Italian wall lizard (Podarcis siculus (Rafinesque-Schmaltz, 1810)) were transplanted from Pod Kopište Island to the nearby island of Pod Mrčaru (Adriatic Sea). Merely 35 years after the introduction, the newly established population on Pod Mrčaru Island had shifted their diet from predominantly insectivorous towards omnivorous and changed significantly in a range of morphological, behavioural, physiological and ecological characteristics. Here, we combine genomic and quantitative genetic approaches to determine the relative roles of genetic adaptation and phenotypic plasticity in driving this rapid phenotypic shift. Our results show genome-wide genetic differentiation between ancestral and transplanted population, with weak genetic erosion on Pod Mrčaru Island. Adaptive processes following the founder event are indicated by highly differentiated genomic loci associating with ecologically relevant phenotypic traits, and/or having a putatively adaptive role across multiple lizard populations. Diverged traits related to head size and shape or bite force showed moderate heritability in a crossing experiment, but between-population differences in these traits did not persist in a common garden environment. Our results confirm the existence of sufficient additive genetic variance for traits to evolve under selection while also demonstrating that phenotypic plasticity and/or genotype by environment interactions are the main drivers of population differentiation at this early evolutionary stage.
Collapse
Affiliation(s)
- Iva Sabolić
- Department of Biology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Óscar Mira
- Department of Biology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Débora Y C Brandt
- Department of Integrative Biology, University of Berkeley, Berkeley, California, USA
| | - Duje Lisičić
- Department of Biology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Jessica Stapley
- Department of Environmental Sciences, ETH Zurich, Zurich, Switzerland
| | - Maria Novosolov
- Lundbeck Foundation GeoGenetics Centre, GLOBE Institute, University of Copenhagen, Copenhagen, Denmark
| | - Robert Bakarić
- Department of Biology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Ivan Cizelj
- Zoological Garden of Zagreb, Zagreb, Croatia
| | - Marko Glogoški
- Department of Biology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | | | - Maxime Taverne
- C.N.R.S/M.N.H.N., Département d'Ecologie et de Gestion de la Biodiversité, Paris, France
| | - Morten E Allentoft
- Lundbeck Foundation GeoGenetics Centre, GLOBE Institute, University of Copenhagen, Copenhagen, Denmark
- Trace and Environmental DNA (TrEnD) Laboratory, School of Molecular and Life Sciences, Curtin University, Perth, Western Australia, Australia
| | - Rasmus Nielsen
- Department of Integrative Biology, University of Berkeley, Berkeley, California, USA
| | - Anthony Herrel
- C.N.R.S/M.N.H.N., Département d'Ecologie et de Gestion de la Biodiversité, Paris, France
- Department of Biology, Evolutionary Morphology of Vertebrates, Ghent University, Ghent, Belgium
- Department of Biology, University of Antwerp, Wilrijk, Belgium
- Naturhistorisches Museum Bern, Bern, Switzerland
| | - Anamaria Štambuk
- Department of Biology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| |
Collapse
|
4
|
Bernatchez L, Ferchaud AL, Berger CS, Venney CJ, Xuereb A. Genomics for monitoring and understanding species responses to global climate change. Nat Rev Genet 2024; 25:165-183. [PMID: 37863940 DOI: 10.1038/s41576-023-00657-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/29/2023] [Indexed: 10/22/2023]
Abstract
All life forms across the globe are experiencing drastic changes in environmental conditions as a result of global climate change. These environmental changes are happening rapidly, incur substantial socioeconomic costs, pose threats to biodiversity and diminish a species' potential to adapt to future environments. Understanding and monitoring how organisms respond to human-driven climate change is therefore a major priority for the conservation of biodiversity in a rapidly changing environment. Recent developments in genomic, transcriptomic and epigenomic technologies are enabling unprecedented insights into the evolutionary processes and molecular bases of adaptation. This Review summarizes methods that apply and integrate omics tools to experimentally investigate, monitor and predict how species and communities in the wild cope with global climate change, which is by genetically adapting to new environmental conditions, through range shifts or through phenotypic plasticity. We identify advantages and limitations of each method and discuss future research avenues that would improve our understanding of species' evolutionary responses to global climate change, highlighting the need for holistic, multi-omics approaches to ecosystem monitoring during global climate change.
Collapse
Affiliation(s)
- Louis Bernatchez
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, Quebec, Canada
| | - Anne-Laure Ferchaud
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, Quebec, Canada.
- Parks Canada, Office of the Chief Ecosystem Scientist, Protected Areas Establishment, Quebec City, Quebec, Canada.
| | - Chloé Suzanne Berger
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, Quebec, Canada
| | - Clare J Venney
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, Quebec, Canada
| | - Amanda Xuereb
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec City, Quebec, Canada
| |
Collapse
|
5
|
Frachon L, Schiestl FP. Rapid genomic evolution in Brassica rapa with bumblebee selection in experimental evolution. BMC Ecol Evol 2024; 24:7. [PMID: 38195402 PMCID: PMC10775529 DOI: 10.1186/s12862-023-02194-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Accepted: 12/20/2023] [Indexed: 01/11/2024] Open
Abstract
BACKGROUND Insect pollinators shape rapid phenotypic evolution of traits related to floral attractiveness and plant reproductive success. However, the underlying genomic changes remain largely unknown despite their importance in predicting adaptive responses to natural or to artificial selection. Based on a nine-generation experimental evolution study with fast cycling Brassica rapa plants adapting to bumblebees, we investigate the genomic evolution associated with the previously observed parallel phenotypic evolution. In this current evolve and resequencing (E&R) study, we conduct a genomic scan of the allele frequency changes along the genome in bumblebee-pollinated and hand-pollinated plants and perform a genomic principal component analysis (PCA). RESULTS We highlight rapid genomic evolution associated with the observed phenotypic evolution mediated by bumblebees. Controlling for genetic drift, we observe significant changes in allelic frequencies at multiple loci. However, this pattern differs according to the replicate of bumblebee-pollinated plants, suggesting putative non-parallel genomic evolution. Finally, our study underlines an increase in genomic variance implying the putative involvement of multiple loci in short-term pollinator adaptation. CONCLUSIONS Overall, our study enhances our understanding of the complex interactions between pollinator and plants, providing a stepping stone towards unravelling the genetic basis of plant genomic adaptation to biotic factors in the environment.
Collapse
Affiliation(s)
- Léa Frachon
- Department of Systematic and Evolutionary Botany, University of Zürich, Zürich, Switzerland.
| | - Florian P Schiestl
- Department of Systematic and Evolutionary Botany, University of Zürich, Zürich, Switzerland
| |
Collapse
|
6
|
Langmüller AM, Nolte V, Dolezal M, Schlötterer C. The genomic distribution of transposable elements is driven by spatially variable purifying selection. Nucleic Acids Res 2023; 51:9203-9213. [PMID: 37560917 PMCID: PMC10516647 DOI: 10.1093/nar/gkad635] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 07/10/2023] [Accepted: 07/18/2023] [Indexed: 08/11/2023] Open
Abstract
It is widely accepted that the genomic distribution of transposable elements (TEs) mainly reflects the outcome of purifying selection and insertion bias (1). Nevertheless, the relative importance of these two evolutionary forces could not be tested thoroughly. Here, we introduce an experimental system, which allows separating purifying selection from TE insertion bias. We used experimental evolution to study the TE insertion patterns in Drosophila simulans founder populations harboring 1040 insertions of an active P-element. After 10 generations at a large population size, we detected strong selection against P-element insertions. The exception were P-element insertions in genomic regions for which a strong insertion bias has been proposed (2-4). Because recurrent P-element insertions cannot explain this pattern, we conclude that purifying selection, with variable strength along the chromosomes, is the major determinant of the genomic distribution of P-elements. Genomic regions with relaxed purifying selection against P-element insertions exhibit normal levels of purifying selection against base substitutions. This suggests that different types of purifying selection operate on base substitutions and P-element insertions. Our results highlight the power of experimental evolution to understand basic evolutionary processes, which are difficult to infer from patterns of natural variation alone.
Collapse
Affiliation(s)
- Anna M Langmüller
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, 1210 Wien, Austria
- Vienna Graduate School of Population Genetics, Vetmeduni Vienna, Veterinärplatz 1, 1210 Vienna, Austria
| | - Viola Nolte
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, 1210 Wien, Austria
| | - Marlies Dolezal
- Plattform Bioinformatik und Biostatistik, Vetmeduni Vienna, Veterinärplatz 1, 1210 Vienna, Austria
| | - Christian Schlötterer
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, 1210 Wien, Austria
| |
Collapse
|
7
|
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
|
8
|
Souto-Maior C, Serrano Negron YL, Harbison ST. Nonlinear expression patterns and multiple shifts in gene network interactions underlie robust phenotypic change in Drosophila melanogaster selected for night sleep duration. PLoS Comput Biol 2023; 19:e1011389. [PMID: 37561813 PMCID: PMC10443883 DOI: 10.1371/journal.pcbi.1011389] [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: 02/15/2023] [Revised: 08/22/2023] [Accepted: 07/25/2023] [Indexed: 08/12/2023] Open
Abstract
All but the simplest phenotypes are believed to result from interactions between two or more genes forming complex networks of gene regulation. Sleep is a complex trait known to depend on the system of feedback loops of the circadian clock, and on many other genes; however, the main components regulating the phenotype and how they interact remain an unsolved puzzle. Genomic and transcriptomic data may well provide part of the answer, but a full account requires a suitable quantitative framework. Here we conducted an artificial selection experiment for sleep duration with RNA-seq data acquired each generation. The phenotypic results are robust across replicates and previous experiments, and the transcription data provides a high-resolution, time-course data set for the evolution of sleep-related gene expression. In addition to a Hierarchical Generalized Linear Model analysis of differential expression that accounts for experimental replicates we develop a flexible Gaussian Process model that estimates interactions between genes. 145 gene pairs are found to have interactions that are different from controls. Our method appears to be not only more specific than standard correlation metrics but also more sensitive, finding correlations not significant by other methods. Statistical predictions were compared to experimental data from public databases on gene interactions. Mutations of candidate genes implicated by our results affected night sleep, and gene expression profiles largely met predicted gene-gene interactions.
Collapse
Affiliation(s)
- Caetano Souto-Maior
- Laboratory of Systems Genetics, Systems Biology Center, National Heart Lung and Blood Institute, Bethesda, Maryland, United States of America
| | - Yazmin L. Serrano Negron
- Laboratory of Systems Genetics, Systems Biology Center, National Heart Lung and Blood Institute, Bethesda, Maryland, United States of America
| | - Susan T. Harbison
- Laboratory of Systems Genetics, Systems Biology Center, National Heart Lung and Blood Institute, Bethesda, Maryland, United States of America
| |
Collapse
|
9
|
Whitehouse LS, Schrider DR. Timesweeper: accurately identifying selective sweeps using population genomic time series. Genetics 2023; 224:iyad084. [PMID: 37157914 PMCID: PMC10324941 DOI: 10.1093/genetics/iyad084] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 07/25/2022] [Accepted: 04/25/2023] [Indexed: 05/10/2023] Open
Abstract
Despite decades of research, identifying selective sweeps, the genomic footprints of positive selection, remains a core problem in population genetics. Of the myriad methods that have been developed to tackle this task, few are designed to leverage the potential of genomic time-series data. This is because in most population genetic studies of natural populations, only a single period of time can be sampled. Recent advancements in sequencing technology, including improvements in extracting and sequencing ancient DNA, have made repeated samplings of a population possible, allowing for more direct analysis of recent evolutionary dynamics. Serial sampling of organisms with shorter generation times has also become more feasible due to improvements in the cost and throughput of sequencing. With these advances in mind, here we present Timesweeper, a fast and accurate convolutional neural network-based tool for identifying selective sweeps in data consisting of multiple genomic samplings of a population over time. Timesweeper analyzes population genomic time-series data by first simulating training data under a demographic model appropriate for the data of interest, training a one-dimensional convolutional neural network on said simulations, and inferring which polymorphisms in this serialized data set were the direct target of a completed or ongoing selective sweep. We show that Timesweeper is accurate under multiple simulated demographic and sampling scenarios, identifies selected variants with high resolution, and estimates selection coefficients more accurately than existing methods. In sum, we show that more accurate inferences about natural selection are possible when genomic time-series data are available; such data will continue to proliferate in coming years due to both the sequencing of ancient samples and repeated samplings of extant populations with faster generation times, as well as experimentally evolved populations where time-series data are often generated. Methodological advances such as Timesweeper thus have the potential to help resolve the controversy over the role of positive selection in the genome. We provide Timesweeper as a Python package for use by the community.
Collapse
Affiliation(s)
- Logan S Whitehouse
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27514, USA
| | - Daniel R Schrider
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27514, USA
| |
Collapse
|
10
|
Burke MK. Embracing Complexity: Yeast Evolution Experiments Featuring Standing Genetic Variation. J Mol Evol 2023; 91:281-292. [PMID: 36752827 PMCID: PMC10276092 DOI: 10.1007/s00239-023-10094-4] [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: 09/26/2022] [Accepted: 01/16/2023] [Indexed: 02/09/2023]
Abstract
The yeast Saccharomyces cerevisiae has a long and esteemed history as a model system for laboratory selection experiments. The majority of yeast evolution experiments begin with an isogenic ancestor, impose selection as cells divide asexually, and track mutations that arise and accumulate over time. Within the last decade, the popularity of S. cerevisiae as a model system for exploring the evolution of standing genetic variation has grown considerably. As a facultatively sexual microbe, it is possible to initiate experiments with populations that harbor diversity and also to maintain that diversity by promoting sexual recombination as the experiment progresses. These experimental choices expand the scope of evolutionary hypotheses that can be tested with yeast. And, in this review, I argue that yeast is one of the best model systems for testing such hypotheses relevant to eukaryotic species. Here, I compile a list of yeast evolution experiments that involve standing genetic variation, initially and/or by implementing protocols that induce sexual recombination in evolving populations. I also provide an overview of experimental methods required to set up such an experiment and discuss the unique challenges that arise in this type of research. Throughout the article, I emphasize the best practices emerging from this small but growing niche of the literature.
Collapse
Affiliation(s)
- Molly K Burke
- Department of Integrative Biology, Oregon State University, Corvallis, OR, 97333, USA.
| |
Collapse
|
11
|
Schlötterer C. How predictable is adaptation from standing genetic variation? Experimental evolution in Drosophila highlights the central role of redundancy and linkage disequilibrium. Philos Trans R Soc Lond B Biol Sci 2023; 378:20220046. [PMID: 37004724 PMCID: PMC10067264 DOI: 10.1098/rstb.2022.0046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2023] Open
Abstract
Experimental evolution is well-suited to test the predictability of evolution without the confounding effects of inaccurate forecasts about future environments. Most of the literature about parallel (and thus predictable) evolution has been carried out in asexual microorganisms, which adapt by de novo mutations. Nevertheless, parallel evolution has also been studied in sexual species at the genomic level. Here, I review the evidence for parallel evolution in Drosophila, the best-studied obligatory outcrossing model for adaptation from standing genetic variation in the laboratory. Similar to asexual microorganisms, evidence for parallel evolution varies between the focal hierarchical levels. Selected phenotypes consistently respond in a very predicable way, but the underlying allele frequency changes are much less predictable. The most important insight is that the predictability of the genomic selection response for polygenic traits depends highly on the founder population and to a much lesser extent on the selection regime. This implies that predicting adaptive genomic response is challenging and requires a good understanding of the adaptive architecture (including linkage disequilibrium) in the ancestral populations. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.
Collapse
Affiliation(s)
- Christian Schlötterer
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, 1210 Wien, Austria
| |
Collapse
|
12
|
Garcia-Costoya G, Williams CE, Faske TM, Moorman JD, Logan ML. Evolutionary constraints mediate extinction risk under climate change. Ecol Lett 2023; 26:529-539. [PMID: 36756845 DOI: 10.1111/ele.14173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Revised: 12/14/2022] [Accepted: 01/13/2023] [Indexed: 02/10/2023]
Abstract
Mounting evidence suggests that rapid evolutionary adaptation may rescue some organisms from the impacts of climate change. However, evolutionary constraints might hinder this process, especially when different aspects of environmental change generate antagonistic selection on genetically correlated traits. Here, we use individual-based simulations to explore how genetic correlations underlying the thermal physiology of ectotherms might influence their responses to the two major components of climate change-increases in mean temperature and thermal variability. We found that genetic correlations can influence population dynamics under climate change, with declines in population size varying three-fold depending on the type of correlation present. Surprisingly, populations whose thermal performance curves were constrained by genetic correlations often declined less rapidly than unconstrained populations. Our results suggest that accurate forecasts of the impact of climate change on ectotherms will require an understanding of the genetic architecture of the traits under selection.
Collapse
Affiliation(s)
| | | | | | - Jacob D Moorman
- University of California, Los Angeles, Los Angeles, California, USA
| | | |
Collapse
|
13
|
Johnson SE, Tittes S, Franks SJ. Rapid, nonparallel genomic evolution of Brassica rapa (field mustard) under experimental drought. J Evol Biol 2023; 36:550-562. [PMID: 36721268 DOI: 10.1111/jeb.14152] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Revised: 09/22/2022] [Accepted: 12/08/2022] [Indexed: 02/02/2023]
Abstract
While we know that climate change can potentially cause rapid phenotypic evolution, our understanding of the genetic basis and degree of genetic parallelism of rapid evolutionary responses to climate change is limited. In this study, we combined the resurrection approach with an evolve-and-resequence design to examine genome-wide evolutionary changes following drought. We exposed genetically similar replicate populations of the annual plant Brassica rapa derived from a field population in southern California to four generations of experimental drought or watered conditions in a greenhouse. Genome-wide sequencing of ancestral and descendant population pools identified hundreds of SNPs that showed evidence of rapidly evolving in response to drought. Several of these were in stress response genes, and two were identified in a prior study of drought response in this species. However, almost all genetic changes were unique among experimental populations, indicating that the evolutionary changes were largely nonparallel, despite the fact that genetically similar replicates of the same founder population had experienced controlled and consistent selection regimes. This nonparallelism of evolution at the genetic level is potentially because of polygenetic adaptation allowing for multiple different genetic routes to similar phenotypic outcomes. Our findings help to elucidate the relationship between rapid phenotypic and genomic evolution and shed light on the degree of parallelism and predictability of genomic evolution to environmental change.
Collapse
Affiliation(s)
- Stephen E Johnson
- Department of Biological Sciences and Louis Calder Center, Fordham University, Bronx, New York, USA
| | - Silas Tittes
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, USA
| | - Steven J Franks
- Department of Biological Sciences and Louis Calder Center, Fordham University, Bronx, New York, USA
| |
Collapse
|
14
|
Resurrection genomics provides molecular and phenotypic evidence of rapid adaptation to salinization in a keystone aquatic species. Proc Natl Acad Sci U S A 2023; 120:e2217276120. [PMID: 36730191 PMCID: PMC9963159 DOI: 10.1073/pnas.2217276120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Ecologists and evolutionary biologists are increasingly cognizant of rapid adaptation in wild populations. Rapid adaptation to anthropogenic environmental change is critical for maintaining biodiversity and ecosystems services into the future. Anthropogenic salinization of freshwater ecosystems is quickly emerging as a primary threat, which is well documented in the northern temperate ecoregion. Specifically, many northern temperate lakes have undergone extensive salinization because of urbanization and the associated increase in impervious surfaces causing runoff, and the extensive use of road deicing salts (e.g., NaCl). It remains unclear whether increasing salinization will lead to extirpation of species from these systems. Using a "resurrection genomics" approach, we investigated whether the keystone aquatic herbivore, Daphnia pulicaria, has evolved increased salinity tolerance in a severely salinized lake located in Minnesota, USA. Whole-genome resequencing of 54 Daphnia clones from the lake and hatched from resting eggs that represent a 25-y temporal contrast demonstrates that many regions of the genome containing genes related to osmoregulation are under selection in the study population. Tolerance assays of clones revealed that the most recent clones are more tolerant to salinity than older clones; this pattern is concomitant with the temporal pattern of stabilizing salinity in this lake. Together, our results demonstrate that keystone species such as Daphnia can rapidly adapt to increasing freshwater salinization. Further, our results indicate that rapid adaptation to salinity may allow lake Daphnia populations to persist in the face of anthropogenic salinization maintaining the food webs and ecosystem services they support despite global environmental change.
Collapse
|
15
|
Stuart KC, Sherwin WB, Edwards RJ, Rollins LA. Evolutionary genomics: Insights from the invasive European starlings. Front Genet 2023; 13:1010456. [PMID: 36685843 PMCID: PMC9845568 DOI: 10.3389/fgene.2022.1010456] [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: 08/03/2022] [Accepted: 11/23/2022] [Indexed: 01/06/2023] Open
Abstract
Two fundamental questions for evolutionary studies are the speed at which evolution occurs, and the way that this evolution may present itself within an organism's genome. Evolutionary studies on invasive populations are poised to tackle some of these pressing questions, including understanding the mechanisms behind rapid adaptation, and how it facilitates population persistence within a novel environment. Investigation of these questions are assisted through recent developments in experimental, sequencing, and analytical protocols; in particular, the growing accessibility of next generation sequencing has enabled a broader range of taxa to be characterised. In this perspective, we discuss recent genetic findings within the invasive European starlings in Australia, and outline some critical next steps within this research system. Further, we use discoveries within this study system to guide discussion of pressing future research directions more generally within the fields of population and evolutionary genetics, including the use of historic specimens, phenotypic data, non-SNP genetic variants (e.g., structural variants), and pan-genomes. In particular, we emphasise the need for exploratory genomics studies across a range of invasive taxa so we can begin understanding broad mechanisms that underpin rapid adaptation in these systems. Understanding how genetic diversity arises and is maintained in a population, and how this contributes to adaptability, requires a deep understanding of how evolution functions at the molecular level, and is of fundamental importance for the future studies and preservation of biodiversity across the globe.
Collapse
Affiliation(s)
- Katarina C. Stuart
- Evolution & Ecology Research Centre, School of Biological, Earth and Environmental Sciences, UNSW Sydney, Sydney, NSW, Australia,*Correspondence: Katarina C. Stuart,
| | - William B. Sherwin
- Evolution & Ecology Research Centre, School of Biological, Earth and Environmental Sciences, UNSW Sydney, Sydney, NSW, Australia
| | - Richard J. Edwards
- Evolution & Ecology Research Centre, School of Biotechnology and Biomolecular Sciences, UNSW Sydney, Sydney, NSW, Australia
| | - Lee A Rollins
- Evolution & Ecology Research Centre, School of Biological, Earth and Environmental Sciences, UNSW Sydney, Sydney, NSW, Australia
| |
Collapse
|
16
|
Aggarwal DD, Mishra P, Singh M. An analysis of direct and indirect effects in Drosophila melanogaster undergoing a few cycles of experimental evolution for stress-related traits. Comp Biochem Physiol B Biochem Mol Biol 2023; 263:110795. [PMID: 35970341 DOI: 10.1016/j.cbpb.2022.110795] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Revised: 08/04/2022] [Accepted: 08/08/2022] [Indexed: 10/15/2022]
Abstract
The physiological mechanisms underpinning adaptations to starvation and cold stresses have been extensively studied in Drosophila, yet the understanding of correlated changes in stress-related and life-history traits, as well as the energetics of stress tolerance, still remains elusive. To answer the questions empirically in this context, we allowed D. melanogaster to evolve for either increased starvation or cold tolerance (24-generations / regime) in an experimental evolution system, and examined whether selection of either trait affects un-selected stress trait, as well as the impacts potential changes in life-history and mating success-related traits. Our results revealed remarkable changes in starvation/cold tolerance (up to 1.5-fold) as a direct effect of selection, while cold tolerance had been dramatically reduced (1.26-fold) in the starvation tolerant (ST) lines compared to control counterparts, although no such changes were evident in cold-tolerant (CT) lines. ST lines exhibited a higher level of body lipids and a reduced level of trehalose content, while CT lines accumulated a greater levels of body lipid and trehalose contents. Noticeably, we found that selection for starvation or cold tolerance positively correlates with larval development time, longevity, and copulation duration, indicating that these traits are among the most common targets of selection trajectories shaping stress tolerance. Altogether, this study highlights the complexity of mechanisms evolved in ST lines that contribute to enhanced starvation tolerance, but also negatively impact cold tolerance. Nevertheless, mechanisms foraging enhanced cold tolerance in CT lines appear not to target starvation tolerance. Moreover, the parallel changes in life history/mating success traits across stress regimes could indicate some generic pathways evolved in stressful environments, targeting life-history and mating success characteristics to optimize fitness.
Collapse
Affiliation(s)
- Dau Dayal Aggarwal
- Department of Biochemistry, University of Delhi South Campus, New Delhi 110021, India.
| | - Prachi Mishra
- Department of Biochemistry, University of Delhi South Campus, New Delhi 110021, India
| | - Manvender Singh
- Department of Biotechnology, University Institute of Technology, Maharshi Dayanand University, Rohtak 124001, India
| |
Collapse
|
17
|
Barata C, Borges R, Kosiol C. Bait-ER: A Bayesian method to detect targets of selection in Evolve-and-Resequence experiments. J Evol Biol 2023; 36:29-44. [PMID: 36544394 PMCID: PMC10108205 DOI: 10.1111/jeb.14134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Revised: 11/09/2022] [Accepted: 11/11/2022] [Indexed: 12/24/2022]
Abstract
For over a decade, experimental evolution has been combined with high-throughput sequencing techniques. In so-called Evolve-and-Resequence (E&R) experiments, populations are kept in the laboratory under controlled experimental conditions where their genomes are sampled and allele frequencies monitored. However, identifying signatures of adaptation in E&R datasets is far from trivial, and it is still necessary to develop more efficient and statistically sound methods for detecting selection in genome-wide data. Here, we present Bait-ER - a fully Bayesian approach based on the Moran model of allele evolution to estimate selection coefficients from E&R experiments. The model has overlapping generations, a feature that describes several experimental designs found in the literature. We tested our method under several different demographic and experimental conditions to assess its accuracy and precision, and it performs well in most scenarios. Nevertheless, some care must be taken when analysing trajectories where drift largely dominates and starting frequencies are low. We compare our method with other available software and report that ours has generally high accuracy even for trajectories whose complexity goes beyond a classical sweep model. Furthermore, our approach avoids the computational burden of simulating an empirical null distribution, outperforming available software in terms of computational time and facilitating its use on genome-wide data. We implemented and released our method in a new open-source software package that can be accessed at https://doi.org/10.5281/zenodo.7351736.
Collapse
Affiliation(s)
- Carolina Barata
- Centre for Biological Diversity, University of St Andrews, St Andrews, UK
| | - Rui Borges
- Institute of Population Genetics, Wien, Austria
| | - Carolin Kosiol
- Centre for Biological Diversity, University of St Andrews, St Andrews, UK.,Institute of Population Genetics, Wien, Austria
| |
Collapse
|
18
|
Cosma BM, Shirali Hossein Zade R, Jordan EN, van Lent P, Peng C, Pillay S, Abeel T. Evaluating long-read de novo assembly tools for eukaryotic genomes: insights and considerations. Gigascience 2022; 12:giad100. [PMID: 38000912 PMCID: PMC10673639 DOI: 10.1093/gigascience/giad100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Revised: 06/18/2023] [Accepted: 10/31/2023] [Indexed: 11/26/2023] Open
Abstract
BACKGROUND Assembly algorithm choice should be a deliberate, well-justified decision when researchers create genome assemblies for eukaryotic organisms from third-generation sequencing technologies. While third-generation sequencing by Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio) has overcome the disadvantages of short read lengths specific to next-generation sequencing (NGS), third-generation sequencers are known to produce more error-prone reads, thereby generating a new set of challenges for assembly algorithms and pipelines. However, the introduction of HiFi reads, which offer substantially reduced error rates, has provided a promising solution for more accurate assembly outcomes. Since the introduction of third-generation sequencing technologies, many tools have been developed that aim to take advantage of the longer reads, and researchers need to choose the correct assembler for their projects. RESULTS We benchmarked state-of-the-art long-read de novo assemblers to help readers make a balanced choice for the assembly of eukaryotes. To this end, we used 12 real and 64 simulated datasets from different eukaryotic genomes, with different read length distributions, imitating PacBio continuous long-read (CLR), PacBio high-fidelity (HiFi), and ONT sequencing to evaluate the assemblers. We include 5 commonly used long-read assemblers in our benchmark: Canu, Flye, Miniasm, Raven, and wtdbg2 for ONT and PacBio CLR reads. For PacBio HiFi reads , we include 5 state-of-the-art HiFi assemblers: HiCanu, Flye, Hifiasm, LJA, and MBG. Evaluation categories address the following metrics: reference-based metrics, assembly statistics, misassembly count, BUSCO completeness, runtime, and RAM usage. Additionally, we investigated the effect of increased read length on the quality of the assemblies and report that read length can, but does not always, positively impact assembly quality. CONCLUSIONS Our benchmark concludes that there is no assembler that performs the best in all the evaluation categories. However, our results show that overall Flye is the best-performing assembler for PacBio CLR and ONT reads, both on real and simulated data. Meanwhile, best-performing PacBio HiFi assemblers are Hifiasm and LJA. Next, the benchmarking using longer reads shows that the increased read length improves assembly quality, but the extent to which that can be achieved depends on the size and complexity of the reference genome.
Collapse
Affiliation(s)
- Bianca-Maria Cosma
- Delft Bioinformatics Lab, Intelligent Systems, Delft University of Technology, 2628 XE, Delft, The Netherlands
| | - Ramin Shirali Hossein Zade
- Delft Bioinformatics Lab, Intelligent Systems, Delft University of Technology, 2628 XE, Delft, The Netherlands
| | - Erin Noel Jordan
- Delft Bioinformatics Lab, Intelligent Systems, Delft University of Technology, 2628 XE, Delft, The Netherlands
- Technical Biochemistry, TU Dortmund University, 44227, Dortmund, Germany
| | - Paul van Lent
- Delft Bioinformatics Lab, Intelligent Systems, Delft University of Technology, 2628 XE, Delft, The Netherlands
| | - Chengyao Peng
- Delft Bioinformatics Lab, Intelligent Systems, Delft University of Technology, 2628 XE, Delft, The Netherlands
| | - Stephanie Pillay
- Delft Bioinformatics Lab, Intelligent Systems, Delft University of Technology, 2628 XE, Delft, The Netherlands
| | - Thomas Abeel
- Delft Bioinformatics Lab, Intelligent Systems, Delft University of Technology, 2628 XE, Delft, The Netherlands
- Infectious Disease and Microbiome Program, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| |
Collapse
|
19
|
Hoedjes KM, Kostic H, Keller L, Flatt T. Natural alleles at the Doa locus underpin evolutionary changes in Drosophila lifespan and fecundity. Proc Biol Sci 2022; 289:20221989. [PMID: 36350205 PMCID: PMC9653240 DOI: 10.1098/rspb.2022.1989] [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] [Indexed: 11/11/2022] Open
Abstract
'Evolve and resequence' (E&R) studies in Drosophila melanogaster have identified many candidate loci underlying the evolution of ageing and life history, but experiments that validate the effects of such candidates remain rare. In a recent E&R study we have identified several alleles of the LAMMER kinase Darkener of apricot (Doa) as candidates for evolutionary changes in lifespan and fecundity. Here, we use two complementary approaches to confirm a functional role of Doa in life-history evolution. First, we used transgenic RNAi to study the effects of Doa at the whole-gene level. Ubiquitous silencing of expression in adult flies reduced both lifespan and fecundity, indicating pleiotropic effects. Second, to characterize segregating variation at Doa, we examined four candidate single nucleotide polymorphisms (SNPs; Doa-1, -2, -3, -4) using a genetic association approach. Three candidate SNPs had effects that were qualitatively consistent with expectations based on our E&R study: Doa-2 pleiotropically affected both lifespan and late-life fecundity; Doa-1 affected lifespan (but not fecundity); and Doa-4 affected late-life fecundity (but not lifespan). Finally, the last candidate allele (Doa-3) also affected lifespan, but in the opposite direction from predicted.
Collapse
Affiliation(s)
- Katja M. Hoedjes
- Department of Ecology and Evolution, University of Lausanne, 1015 Lausanne, Switzerland
| | - Hristina Kostic
- Department of Ecology and Evolution, University of Lausanne, 1015 Lausanne, Switzerland
| | - Laurent Keller
- Department of Ecology and Evolution, University of Lausanne, 1015 Lausanne, Switzerland
| | - Thomas Flatt
- Department of Ecology and Evolution, University of Lausanne, 1015 Lausanne, Switzerland,Department of Biology, University of Fribourg, 1700 Fribourg, Switzerland
| |
Collapse
|
20
|
Burny C, Nolte V, Dolezal M, Schlötterer C. Genome-wide selection signatures reveal widespread synergistic effects of two different stressors in Drosophila melanogaster. Proc Biol Sci 2022; 289:20221857. [PMID: 36259211 PMCID: PMC9579754 DOI: 10.1098/rspb.2022.1857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Experimental evolution combined with whole-genome sequencing (evolve and resequence (E&R)) is a powerful approach to study the adaptive architecture of selected traits. Nevertheless, so far the focus has been on the selective response triggered by a single stressor. Building on the highly parallel selection response of founder populations with reduced variation, we evaluated how the presence of a second stressor affects the genomic selection response. After 20 generations of adaptation to laboratory conditions at either 18°C or 29°C, strong genome-wide selection signatures were observed. Only 38% of the selection signatures can be attributed to laboratory adaptation (no difference between temperature regimes). The remaining selection responses are either caused by temperature-specific effects, or reflect the joint effects of temperature and laboratory adaptation (same direction, but the magnitude differs between temperatures). The allele frequency changes resulting from the combined effects of temperature and laboratory adaptation were more extreme in the hot environment for 83% of the affected genomic regions-indicating widespread synergistic effects of the two stressors. We conclude that E&R with reduced genetic variation is a powerful approach to study genome-wide fitness consequences driven by the combined effects of multiple environmental factors.
Collapse
Affiliation(s)
- Claire Burny
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Vienna 1210, Austria.,Vienna Graduate School of Population Genetics, Vetmeduni Vienna, Vienna 1210, Austria
| | - Viola Nolte
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Vienna 1210, Austria
| | - Marlies Dolezal
- Plattform Bioinformatik und Biostatistik, Vetmeduni Vienna, Vienna 1210, Austria
| | - Christian Schlötterer
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Vienna 1210, Austria
| |
Collapse
|
21
|
Roberts RJV, Pop S, Prieto-Godino LL. Evolution of central neural circuits: state of the art and perspectives. Nat Rev Neurosci 2022; 23:725-743. [DOI: 10.1038/s41583-022-00644-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/22/2022] [Indexed: 11/09/2022]
|
22
|
Visher E, Uricchio L, Bartlett L, DeNamur N, Yarcan A, Alhassani D, Boots M. The evolution of host specialization in an insect pathogen. Evolution 2022; 76:2375-2388. [PMID: 35946063 DOI: 10.1111/evo.14594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Revised: 05/21/2022] [Accepted: 06/06/2022] [Indexed: 01/22/2023]
Abstract
Niche breadth coevolution between biotic partners underpins theories of diversity and co-existence and influences patterns of disease emergence and transmission in host-parasite systems. Despite these broad implications, we still do not fully understand how the breadth of parasites' infectivity evolves, the nature of any associated costs, or the genetic basis of specialization. Here, we serially passage a granulosis virus on multiple inbred populations of its Plodia interpunctella host to explore the dynamics and outcomes of specialization. In particular, we collect time series of phenotypic and genetic data to explore the dynamics of host genotype specialization throughout the course of experimental evolution and examine two fitness components. We find that the Plodia interpunctella granulosis virus consistently evolves and increases in overall specialization, but that our two fitness components evolve independently such that lines can specialize in productivity or infectivity. Furthermore, we find that specialization in our experiment is a highly polygenic trait best explained by a combination of evolutionary mechanisms. These results are important for understanding the evolution of specialization in host-parasite interactions and its broader implications for co-existence, diversification, and infectious disease management.
Collapse
Affiliation(s)
- Elisa Visher
- Department of Integrative Biology, University of California, Berkeley, CA, 94720, USA
| | | | - Lewis Bartlett
- Center for the Ecology of Infectious Diseases, University of Georgia, Athens, GA, 30602, USA
| | | | - Aren Yarcan
- University of California, Berkeley, CA, 94720, USA
| | | | - Mike Boots
- Department of Integrative Biology, University of California, Berkeley, CA, 94720, USA.,Centre for Ecology and Conservation, College of Life and Environmental Sciences, University of Exeter Penryn Campus, Penryn, TR10 9FE, UK
| |
Collapse
|
23
|
Experimental evolution reveals the synergistic genomic mechanisms of adaptation to ocean warming and acidification in a marine copepod. Proc Natl Acad Sci U S A 2022; 119:e2201521119. [PMID: 36095205 PMCID: PMC9499500 DOI: 10.1073/pnas.2201521119] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Metazoan adaptation to global change relies on selection of standing genetic variation. Determining the extent to which this variation exists in natural populations, particularly for responses to simultaneous stressors, is essential to make accurate predictions for persistence in future conditions. Here, we identified the genetic variation enabling the copepod Acartia tonsa to adapt to experimental ocean warming, acidification, and combined ocean warming and acidification (OWA) over 25 generations of continual selection. Replicate populations showed a consistent polygenic response to each condition, targeting an array of adaptive mechanisms including cellular homeostasis, development, and stress response. We used a genome-wide covariance approach to partition the allelic changes into three categories: selection, drift and replicate-specific selection, and laboratory adaptation responses. The majority of allele frequency change in warming (57%) and OWA (63%) was driven by shared selection pressures across replicates, but this effect was weaker under acidification alone (20%). OWA and warming shared 37% of their response to selection but OWA and acidification shared just 1%, indicating that warming is the dominant driver of selection in OWA. Despite the dominance of warming, the interaction with acidification was still critical as the OWA selection response was highly synergistic with 47% of the allelic selection response unique from either individual treatment. These results disentangle how genomic targets of selection differ between single and multiple stressors and demonstrate the complexity that nonadditive multiple stressors will contribute to predictions of adaptation to complex environmental shifts caused by global change.
Collapse
|
24
|
Sherwin WB. Bray-Curtis (AFD) differentiation in molecular ecology: Forecasting, an adjustment ( A A), and comparative performance in selection detection. Ecol Evol 2022; 12:e9176. [PMID: 36110882 PMCID: PMC9465203 DOI: 10.1002/ece3.9176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Revised: 07/04/2022] [Accepted: 07/06/2022] [Indexed: 11/07/2022] Open
Abstract
Geographic genetic differentiation measures are used for purposes such as assessing genetic diversity and connectivity, and searching for signals of selection. Confirmation by unrelated measures can minimize false positives. A popular differentiation measure, Bray‐Curtis, has been used increasingly in molecular ecology, renamed AFD (hereafter called BCAFD). Critically, BCAFD is expected to be partially independent of the commonly used Hill “Q‐profile” measures. BCAFD needs scrutiny for potential biases, by examining limits on its value, and comparing simulations against expectations. BCAFD has two dependencies on within‐population (alpha) variation, undesirable for a between‐population (beta) measure. The first dependency is derived from similarity to GST and FST. The second dependency is that BCAFD cannot be larger than the highest allele proportion in either location (alpha variation), which can be overcome by data‐filtering or by a modified statistic AA or “Adjusted AFD”. The first dependency does not forestall applications such as assessing connectivity or selection, if we know the measure's null behavior under selective neutrality with specified conditions—which is shown in this article for AA, for equilibrium, and nonequilibrium, for the commonly used data type of single‐nucleotide‐polymorphisms (SNPs) in two locations. Thus, AA can be used in tandem with mathematically contrasting differentiation measures, with the aim of reducing false inferences. For detecting adaptive loci, the relative performance of AA and other measures was evaluated, showing that it is best to use two mathematically different measures simultaneously, and that AA is in one of the best such pairwise criteria. For any application, using AA, rather than BCAFD, avoids the counterintuitive limitation by maximum allele proportion within localities.
Collapse
Affiliation(s)
- William B Sherwin
- Evolution and Ecology Research Centre, School of BEES UNSW-Sydney Sydney New South Wales Australia
| |
Collapse
|
25
|
Olivares‐Zambrano D, Daane J, Hyde J, Sandel MW, Aguilar A. Speciation genomics and the role of depth in the divergence of rockfishes (
Sebastes
) revealed through Pool‐seq analysis of enriched sequences. Ecol Evol 2022; 12:e9341. [PMID: 36188524 PMCID: PMC9502067 DOI: 10.1002/ece3.9341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2022] [Revised: 08/22/2022] [Accepted: 08/30/2022] [Indexed: 11/07/2022] Open
Abstract
Speciation in the marine environment is challenged by the wide geographic distribution of many taxa and potential for high rates of gene flow through larval dispersal mechanisms. Depth has recently been proposed as a potential driver of ecological divergence in fishes, and yet it is unclear how adaptation along these gradients' shapes genomic divergence. The genus Sebastes contains numerous species pairs that are depth‐segregated and can provide a better understanding of the mode and tempo of genomic diversification. Here, we present exome data on two species pairs of rockfishes that are depth‐segregated and have different degrees of divergence: S. chlorostictus–S. rosenblatti and S. crocotulus–S. miniatus. We were able to reliably identify “islands of divergence” in the species pair with more recent divergence (S. chlorostictus–S. rosenblatti) and discovered a number of genes associated with neurosensory function, suggesting a role for this pathway in the early speciation process. We also reconstructed demographic histories of divergence and found the best supported model was isolation followed by asymmetric secondary contact for both species pairs. These results suggest past ecological/geographic isolation followed by asymmetric secondary contact of deep to shallow species. Our results provide another example of using rockfish as a model for studying speciation and support the role of depth as an important mechanism for diversification in the marine environment.
Collapse
Affiliation(s)
- Daniel Olivares‐Zambrano
- Department of Biological SciencesCalifornia State University Los AngelesLos AngelesCaliforniaUSA
- Present address:
Department of Marine and Environmental BiologyUniversity of Southern CaliforniaLos AngelesCaliforniaUSA
| | - Jacob Daane
- Department of Biology and BiochemistryUniversity of HoustonHoustonTexasUSA
| | - John Hyde
- National Oceanic and Atmospheric Administration, National Marine Fisheries ServiceNational Marine Fisheries ServiceSouthwest Fisheries Science CenterLa JollaCaliforniaUSA
| | - Michael W. Sandel
- Biological and Environmental SciencesUniversity of West AlabamaLivingstonAlabamaUSA
- Department of WIldlifeFisheries, and Aquaculture, Mississippi State UniversityMississippi StateMississippiUSA
| | - Andres Aguilar
- Department of Biological SciencesCalifornia State University Los AngelesLos AngelesCaliforniaUSA
| |
Collapse
|
26
|
Host-Associated Rhizobial Fitness: Dependence on Nitrogen, Density, Community Complexity, and Legume Genotype. Appl Environ Microbiol 2022; 88:e0052622. [PMID: 35852362 PMCID: PMC9361818 DOI: 10.1128/aem.00526-22] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
The environmental context of the nitrogen-fixing mutualism between leguminous plants and rhizobial bacteria varies over space and time. Variation in resource availability, population density, and composition likely affect the ecology and evolution of rhizobia and their symbiotic interactions with hosts. We examined how host genotype, nitrogen addition, rhizobial density, and community complexity affected selection on 68 rhizobial strains in the Sinorhizobium meliloti–Medicago truncatula mutualism. As expected, host genotype had a substantial effect on the size, number, and strain composition of root nodules (the symbiotic organ). The understudied environmental variable of rhizobial density had a stronger effect on nodule strain frequency than the addition of low nitrogen levels. Higher inoculum density resulted in a nodule community that was less diverse and more beneficial but only in the context of the more selective host genotype. Higher density resulted in more diverse and less beneficial nodule communities with the less selective host. Density effects on strain composition deserve additional scrutiny as they can create feedback between ecological and evolutionary processes. Finally, we found that relative strain rankings were stable across increasing community complexity (2, 3, 8, or 68 strains). This unexpected result suggests that higher-order interactions between strains are rare in the context of nodule formation and development. Our work highlights the importance of examining mechanisms of density-dependent strain fitness and developing theoretical predictions that incorporate density dependence. Furthermore, our results have translational relevance for overcoming establishment barriers in bioinoculants and motivating breeding programs that maintain beneficial plant-microbe interactions across diverse agroecological contexts. IMPORTANCE Legume crops establish beneficial associations with rhizobial bacteria that perform biological nitrogen fixation, providing nitrogen to plants without the economic and greenhouse gas emission costs of chemical nitrogen inputs. Here, we examine the influence of three environmental factors that vary in agricultural fields on strain relative fitness in nodules. In addition to manipulating nitrogen, we also use two biotic variables that have rarely been examined: the rhizobial community's density and complexity. Taken together, our results suggest that (i) breeding legume varieties that select beneficial strains despite environmental variation is possible, (ii) changes in rhizobial population densities that occur routinely in agricultural fields could drive evolutionary changes in rhizobial populations, and (iii) the lack of higher-order interactions between strains will allow the high-throughput assessments of rhizobia winners and losers during plant interactions.
Collapse
|
27
|
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
|
28
|
Khosrovyan A, Doria HB, Kahru A, Pfenninger M. Polyamide microplastic exposure elicits rapid, strong and genome-wide evolutionary response in the freshwater non-biting midge Chironomus riparius. CHEMOSPHERE 2022; 299:134452. [PMID: 35367228 DOI: 10.1016/j.chemosphere.2022.134452] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 03/20/2022] [Accepted: 03/25/2022] [Indexed: 06/14/2023]
Abstract
Susceptibility to hazardous materials and contamination is largely determined by genetic make-up and evolutionary history of affected organisms. Yet evolutionary adaptation and microevolutionary processes triggered by contaminants are rarely considered in ecotoxicology. Using an evolve and resequencing approach, we investigated genome-wide responses of the midge C. riparius exposed to virgin polyamide microplastics (0-180 μm size range, at concentration 1 g kg-1) during seven consecutive generations. The results were integrated to a parallel life-cycle experiment ran under the same exposure conditions. Emergence, life-cycle trait, showed first a substantial reduction in larval survival, followed by a rapid recovery within three generations. On the genomic level, we observed substantial selectively driven allele frequency changes (mean 0.566 ± 0.0879) within seven generations, associated with a mean selection coefficient of 0.322, indicating very strong selection pressure. Putative selection targets were mainly connected to oxidative stress in the microplastics exposed C. riparius population. This is the first multigenerational study on chironomids to provide evidence that upon exposure to polyamide microplastic there are changes on the genomic level, providing basis to rapid adaptation of aquatic organisms to microplastics.
Collapse
Affiliation(s)
- Alla Khosrovyan
- National Institute of Chemical Physics and Biophysics, Laboratory of Environmental Toxicology, 23 Akadeemia Tee, 12618, Tallinn, Estonia.
| | - Halina Binde Doria
- Dept. Molecular Ecology, Senckenberg Biodiversity and Climate Research Centre, Georg-Voigt-Str. 14-16, D-60325, Frankfurt am Main, Germany; LOEWE Centre for Translational Biodiversity Genomics, Senckenberg Biodiversity and Climate Research Centre, Senckenberganlage 25, 60325, Frankfurt am Main, Germany.
| | - Anne Kahru
- National Institute of Chemical Physics and Biophysics, Laboratory of Environmental Toxicology, 23 Akadeemia Tee, 12618, Tallinn, Estonia; Estonian Academy of Sciences, 6 Kohtu, 10130, Tallinn, Estonia
| | - Markus Pfenninger
- Dept. Molecular Ecology, Senckenberg Biodiversity and Climate Research Centre, Georg-Voigt-Str. 14-16, D-60325, Frankfurt am Main, Germany; LOEWE Centre for Translational Biodiversity Genomics, Senckenberg Biodiversity and Climate Research Centre, Senckenberganlage 25, 60325, Frankfurt am Main, Germany; Institute for Molecular and Organismic Evolution, Johannes Gutenberg University, Johann-Joachim-Becher-Weg 7, 55128, Mainz, Germany
| |
Collapse
|
29
|
Lahondère C, Bonizzoni M. Thermal biology of invasive Aedes mosquitoes in the context of climate change. CURRENT OPINION IN INSECT SCIENCE 2022; 51:100920. [PMID: 35421621 DOI: 10.1016/j.cois.2022.100920] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Revised: 02/24/2022] [Accepted: 03/30/2022] [Indexed: 05/06/2023]
Abstract
The increasing incidence of arboviral diseases in tropical endemic areas and their emergence in new temperate countries is one of the most important challenges that Public Health agencies are currently facing. Because mosquitoes are poikilotherms, shifts in temperature influence physiological functions besides egg viability. These traits impact not only vector density, but also their interaction with their hosts and arboviruses. As such the relationship among mosquitoes, arboviral diseases and temperature is complex. Here, we summarize current knowledge on the thermal biology of Aedes invasive mosquitoes, highlighting differences among species. We also emphasize the need to expand knowledge on the variability in thermal sensitivity across populations within a species, especially in light of climate change that encompasses increase not only in mean environmental temperature but also in the frequency of hot and cold snaps. Finally, we suggest a novel experimental approach to investigate the molecular architecture of thermal adaptation in mosquitoes.
Collapse
Affiliation(s)
- Chloé Lahondère
- Department of Biochemistry, USA; The Fralin Life Science Institute, USA; Center of Emerging, Zoonotic and Arthropod-borne Pathogens, USA; The Global Change Center, USA; Department of Entomology at Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA.
| | | |
Collapse
|
30
|
Brennan RS, deMayo JA, Dam HG, Finiguerra MB, Baumann H, Pespeni MH. Loss of transcriptional plasticity but sustained adaptive capacity after adaptation to global change conditions in a marine copepod. Nat Commun 2022; 13:1147. [PMID: 35241657 PMCID: PMC8894427 DOI: 10.1038/s41467-022-28742-6] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Accepted: 02/04/2022] [Indexed: 12/04/2022] Open
Abstract
Adaptive evolution and phenotypic plasticity will fuel resilience in the geologically unprecedented warming and acidification of the earth’s oceans, however, we have much to learn about the interactions and costs of these mechanisms of resilience. Here, using 20 generations of experimental evolution followed by three generations of reciprocal transplants, we investigated the relationship between adaptation and plasticity in the marine copepod, Acartia tonsa, in future global change conditions (high temperature and high CO2). We found parallel adaptation to global change conditions in genes related to stress response, gene expression regulation, actin regulation, developmental processes, and energy production. However, reciprocal transplantation showed that adaptation resulted in a loss of transcriptional plasticity, reduced fecundity, and reduced population growth when global change-adapted animals were returned to ambient conditions or reared in low food conditions. However, after three successive transplant generations, global change-adapted animals were able to match the ambient-adaptive transcriptional profile. Concurrent changes in allele frequencies and erosion of nucleotide diversity suggest that this recovery occurred via adaptation back to ancestral conditions. These results demonstrate that while plasticity facilitated initial survival in global change conditions, it eroded after 20 generations as populations adapted, limiting resilience to new stressors and previously benign environments. Rapid adaptation will facilitate species resilience under global climate change, but its effects on plasticity are less commonly investigated. This study shows that 20 generations of experimental adaptation in a marine copepod drives a rapid loss of plasticity that carries costs and might have impacts on future resilience to environmental change.
Collapse
Affiliation(s)
- Reid S Brennan
- Department of Biology, University of Vermont, Burlington, VT, USA. .,Marine Evolutionary Ecology, GEOMAR Helmholtz Centre for Ocean Research Kiel, Kiel, Germany.
| | - James A deMayo
- Department of Marine Sciences, University of Connecticut, Groton, CT, USA.,Department of Integrative Biology, University of Colorado Denver, Denver, CO, USA
| | - Hans G Dam
- Department of Marine Sciences, University of Connecticut, Groton, CT, USA
| | - Michael B Finiguerra
- Department of Ecology and Evolutionary Biology, University of Connecticut, Groton, CT, USA
| | - Hannes Baumann
- Department of Marine Sciences, University of Connecticut, Groton, CT, USA
| | | |
Collapse
|
31
|
Kasimatis KR, Moerdyk-Schauwecker MJ, Lancaster R, Smith A, Willis JH, Phillips PC. Post-insemination selection dominates pre-insemination selection in driving rapid evolution of male competitive ability. PLoS Genet 2022; 18:e1010063. [PMID: 35157717 PMCID: PMC8880957 DOI: 10.1371/journal.pgen.1010063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Revised: 02/25/2022] [Accepted: 01/28/2022] [Indexed: 11/30/2022] Open
Abstract
Sexual reproduction is a complex process that contributes to differences between the sexes and divergence between species. From a male’s perspective, sexual selection can optimize reproductive success by acting on the variance in mating success (pre-insemination selection) as well as the variance in fertilization success (post-insemination selection). The balance between pre- and post-insemination selection has not yet been investigated using a strong hypothesis-testing framework that directly quantifies the effects of post-insemination selection on the evolution of reproductive success. Here we use experimental evolution of a uniquely engineered genetic system that allows sperm production to be turned off and on in obligate male-female populations of Caenorhabditis elegans. We show that enhanced post-insemination competition increases the efficacy of selection and surpasses pre-insemination sexual selection in driving a polygenic response in male reproductive success. We find that after 10 selective events occurring over 30 generations post-insemination selection increased male reproductive success by an average of 5- to 7-fold. Contrary to expectation, enhanced pre-insemination competition hindered selection and slowed the rate of evolution. Furthermore, we found that post-insemination selection resulted in a strong polygenic response at the whole-genome level. Our results demonstrate that post-insemination sexual selection plays a critical role in the rapid optimization of male reproductive fitness. Therefore, explicit consideration should be given to post-insemination dynamics when considering the population effects of sexual selection. Some of the most dramatic and diverse phenotypes observed in nature––such as head-butting in wild sheep and the elaborate tails of peacocks––are sexually dimorphic. These remarkable phenotypes are a result of sexual selection optimizing reproductive success in females and males independently. For males, total reproductive success is comprised of winning a mating event and then translating that mating event into a fertilization event. Therefore, to understand not only how male reproductive success is comprised, but also how it evolves, we must examine the interaction between pre- and post-insemination sexual selection. We combine environmentally-inducible control of sperm production within a highly reproducible factorial experimental evolution design to directly quantify the contribution of post-insemination selection to male reproductive evolution. We demonstrate that enhanced sperm competition increases the efficacy of selection and enhances the rate of male evolution. Alternatively, we show that enhanced pre-insemination competition slows the evolutionary rate. Using whole-genome approaches, we identify over 60 genes that contribute to male fertilization success. Brought together, our new approaches and results demonstrate that the unseen world of molecular interactions occurring during post-insemination are as fundamentally important as pre-mating factors.
Collapse
Affiliation(s)
- Katja R. Kasimatis
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
- * E-mail: (KRK); (PCP)
| | | | - Ruben Lancaster
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - Alexander Smith
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - John H. Willis
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - Patrick C. Phillips
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
- * E-mail: (KRK); (PCP)
| |
Collapse
|
32
|
Phillips MA, Arnold KR, Vue Z, Beasley HK, Garza-Lopez E, Marshall AG, Morton DJ, McReynolds MR, Barter TT, Hinton A. Combining Metabolomics and Experimental Evolution Reveals Key Mechanisms Underlying Longevity Differences in Laboratory Evolved Drosophila melanogaster Populations. Int J Mol Sci 2022; 23:ijms23031067. [PMID: 35162994 PMCID: PMC8835531 DOI: 10.3390/ijms23031067] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 01/07/2022] [Accepted: 01/11/2022] [Indexed: 12/22/2022] Open
Abstract
Experimental evolution with Drosophila melanogaster has been used extensively for decades to study aging and longevity. In recent years, the addition of DNA and RNA sequencing to this framework has allowed researchers to leverage the statistical power inherent to experimental evolution to study the genetic basis of longevity itself. Here, we incorporated metabolomic data into to this framework to generate even deeper insights into the physiological and genetic mechanisms underlying longevity differences in three groups of experimentally evolved D. melanogaster populations with different aging and longevity patterns. Our metabolomic analysis found that aging alters mitochondrial metabolism through increased consumption of NAD+ and increased usage of the TCA cycle. Combining our genomic and metabolomic data produced a list of biologically relevant candidate genes. Among these candidates, we found significant enrichment for genes and pathways associated with neurological development and function, and carbohydrate metabolism. While we do not explicitly find enrichment for aging canonical genes, neurological dysregulation and carbohydrate metabolism are both known to be associated with accelerated aging and reduced longevity. Taken together, our results provide plausible genetic mechanisms for what might be driving longevity differences in this experimental system. More broadly, our findings demonstrate the value of combining multiple types of omic data with experimental evolution when attempting to dissect mechanisms underlying complex and highly polygenic traits such as aging.
Collapse
Affiliation(s)
- Mark A. Phillips
- Department of Integrative Biology, Oregon State University, Corvallis, OR 97331, USA;
| | - Kenneth R. Arnold
- Department of Ecology and Evolutionary Biology, University of California, Irvine, CA 92697, USA; (K.R.A.); (T.T.B.)
| | - Zer Vue
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA; (Z.V.); (H.K.B.); (A.G.M.)
| | - Heather K. Beasley
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA; (Z.V.); (H.K.B.); (A.G.M.)
- Department of Biochemistry, Cancer Biology, Neuroscience, and Pharmacology, Meharry Medical College, Nashville, TN 37208, USA
| | - Edgar Garza-Lopez
- Hinton and Garza-Lopez Family Consulting Company, Iowa City, IA 52246, USA;
| | - Andrea G. Marshall
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA; (Z.V.); (H.K.B.); (A.G.M.)
| | - Derrick J. Morton
- Department of Biological Sciences, University of Southern California, Los Angeles, CA 90089, USA;
| | - Melanie R. McReynolds
- Department of Biochemistry and Molecular Biology, Huck Institute of the Life Sciences, Pennsylvania State University, University Park, PA 16802, USA;
| | - Thomas T. Barter
- Department of Ecology and Evolutionary Biology, University of California, Irvine, CA 92697, USA; (K.R.A.); (T.T.B.)
| | - Antentor Hinton
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA; (Z.V.); (H.K.B.); (A.G.M.)
- Hinton and Garza-Lopez Family Consulting Company, Iowa City, IA 52246, USA;
- Correspondence:
| |
Collapse
|
33
|
Galewski P, Funk A, McGrath JM. Select and Sequence of a Segregating Sugar Beet Population Provides Genomic Perspective of Host Resistance to Seedling Rhizoctonia solani Infection. FRONTIERS IN PLANT SCIENCE 2022; 12:785267. [PMID: 35095959 PMCID: PMC8793884 DOI: 10.3389/fpls.2021.785267] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Accepted: 11/12/2021] [Indexed: 05/15/2023]
Abstract
Understanding the genetic basis of polygenic traits is a major challenge in agricultural species, especially in non-model systems. Select and sequence (SnS) experiments carried out within existing breeding programs provide a means to simultaneously identify the genomic background of a trait while improving the mean phenotype for a population. Using pooled whole genome sequencing (WGS) of selected and unselected bulks derived from a synthetic outcrossing sugar beet population EL57 (PI 663212), which segregates for seedling rhizoctonia resistance, we identified a putative genomic background involved in conditioning a resistance phenotype. Population genomic parameters were estimated to measure fixation (He), genome divergence (F ST ), and allele frequency changes between bulks (DeltaAF). We report on the genome wide patterns of variation resulting from selection and highlight specific genomic features associated with resistance. Expected heterozygosity (He) showed an increased level of fixation in the resistant bulk, indicating a greater selection pressure was applied. In total, 1,311 biallelic loci were detected as significant FST outliers (p < 0.01) in comparisons between the resistant and susceptible bulks. These loci were detected in 206 regions along the chromosomes and contained 275 genes. We estimated changes in allele frequency between bulks resulting from selection for resistance by leveraging the allele frequencies of an unselected bulk. DeltaAF was a more stringent test of selection and recovered 186 significant loci, representing 32 genes, all of which were also detected using FST. Estimates of population genetic parameters and statistical significance were visualized with respect to the EL10.2 physical map and produced a candidate gene list that was enriched for function in cell wall metabolism and plant disease resistance, including pathogen perception, signal transduction, and pathogen response. Specific variation associated with these genes was also reported and represents genetic markers for validation and prediction of resistance to Rhizoctonia. Select and sequence experiments offer a means to characterize the genetic base of sugar beet, inform selection within breeding programs, and prioritize candidate variation for functional studies.
Collapse
Affiliation(s)
- Paul Galewski
- United States Department of Agriculture – Agricultural Research Service (USDA-ARS) Northwest Irrigation and Soils Research Laboratory, Kimberly, ID, United States
- Department of Plant, Soil, and Microbial Science, Plant Breeding, Genetics, and Biotechnology Program, Michigan State University, East Lansing, MI, United States
| | - Andrew Funk
- Department of Plant, Soil, and Microbial Science, Plant Breeding, Genetics, and Biotechnology Program, Michigan State University, East Lansing, MI, United States
- United States Department of Agriculture – National Institute of Food and Agriculture (USDA-NIFA) Institute of Food Production and Sustainability, Kansas City, MO, United States
| | - J. Mitchell McGrath
- United States Department of Agriculture – Agricultural Research Service (USDA-ARS) Sugar Beet and Bean Research Unit USDA-ARS, East Lansing, MI, United States
| |
Collapse
|
34
|
Saleh D, Chen J, Leplé J, Leroy T, Truffaut L, Dencausse B, Lalanne C, Labadie K, Lesur I, Bert D, Lagane F, Morneau F, Aury J, Plomion C, Lascoux M, Kremer A. Genome-wide evolutionary response of European oaks during the Anthropocene. Evol Lett 2022; 6:4-20. [PMID: 35127134 PMCID: PMC8802238 DOI: 10.1002/evl3.269] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2021] [Revised: 11/26/2021] [Accepted: 12/02/2021] [Indexed: 12/23/2022] Open
Abstract
The pace of tree microevolution during Anthropocene warming is largely unknown. We used a retrospective approach to monitor genomic changes in oak trees since the Little Ice Age (LIA). Allelic frequency changes were assessed from whole-genome pooled sequences for four age-structured cohorts of sessile oak (Quercus petraea) dating back to 1680, in each of three different oak forests in France. The genetic covariances of allelic frequency changes increased between successive time periods, highlighting genome-wide effects of linked selection. We found imprints of parallel linked selection in the three forests during the late LIA, and a shift of selection during more recent time periods of the Anthropocene. The changes in allelic covariances within and between forests mirrored the documented changes in the occurrence of extreme events (droughts and frosts) over the last 300 years. The genomic regions with the highest covariances were enriched in genes involved in plant responses to pathogens and abiotic stresses (temperature and drought). These responses are consistent with the reported sequence of frost (or drought) and disease damage ultimately leading to the oak dieback after extreme events. They provide support for adaptive evolution of long-lived species during recent climatic changes. Although we acknowledge that other sources (e.g., gene flow, generation overlap) may have contributed to temporal covariances of allelic frequency changes, the consistent and correlated response across the three forests lends support to the existence of a systematic driving force such as natural selection.
Collapse
Affiliation(s)
- Dounia Saleh
- UMR BIOGECO, INRAEUniversité de BordeauxCestas33612France
| | - Jun Chen
- College of Life SciencesZhejiang UniversityHangzhou310058China
| | | | - Thibault Leroy
- Department of Botany and Biodiversity ResearchUniversity of ViennaVienna1010Austria
| | - Laura Truffaut
- UMR BIOGECO, INRAEUniversité de BordeauxCestas33612France
| | | | - Céline Lalanne
- UMR BIOGECO, INRAEUniversité de BordeauxCestas33612France
| | - Karine Labadie
- Genoscope, Institut de Biologie François Jacob, Commissariat à l’énergie atomique (CEA)Université de Paris‐SaclayEvry91057France
| | | | - Didier Bert
- UMR BIOGECO, INRAEUniversité de BordeauxCestas33612France
| | | | - François Morneau
- Département Recherche Développement InnovationOffice National des ForêtsBoigny‐Sur‐Bionne45760France,Current Address: Service de l'Information Statistique Forestière et EnvironnementaleInstitut National de l'Information géographique et ForestièreNogent‐sur‐Vernisson45290France
| | - Jean‐Marc Aury
- Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRSUniv Evry, Université Paris‐SaclayEvry91057France
| | | | - Martin Lascoux
- Department of Ecology and Genetics, Evolutionary Biology CentreUppsala UniversityUppsalaSE‐75236Sweden
| | - Antoine Kremer
- UMR BIOGECO, INRAEUniversité de BordeauxCestas33612France
| |
Collapse
|
35
|
Lewis JA, Morran LT. Advantages of laboratory natural selection in the applied sciences. J Evol Biol 2021; 35:5-22. [PMID: 34826161 DOI: 10.1111/jeb.13964] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Revised: 11/22/2021] [Accepted: 11/23/2021] [Indexed: 11/29/2022]
Abstract
In the past three decades, laboratory natural selection has become a widely used technique in biological research. Most studies which have utilized this technique are in the realm of basic science, often testing hypotheses related to mechanisms of evolutionary change or ecological dynamics. While laboratory natural selection is currently utilized heavily in this setting, there is a significant gap with its usage in applied studies, especially when compared to the other selection experiment methodologies like artificial selection and directed evolution. This is despite avenues of research in the applied sciences which seem well suited to laboratory natural selection. In this review, we place laboratory natural selection in context with other selection experiments, identify the characteristics which make it well suited for particular kinds of applied research and briefly cover key examples of the usefulness of selection experiments within applied science. Finally, we identify three promising areas of inquiry for laboratory natural selection in the applied sciences: bioremediation technology, identifying mechanisms of drug resistance and optimizing biofuel production. Although laboratory natural selection is currently less utilized in applied science when compared to basic research, the method has immense promise in the field moving forward.
Collapse
Affiliation(s)
- Jordan A Lewis
- Population Biology, Ecology, and Evolution Graduate Program, Emory University, Atlanta, Georgia, USA
| | - Levi T Morran
- Population Biology, Ecology, and Evolution Graduate Program, Emory University, Atlanta, Georgia, USA.,Department of Biology, Emory University, Atlanta, Georgia, USA
| |
Collapse
|
36
|
Burny C, Nolte V, Dolezal M, Schlötterer C. Highly Parallel Genomic Selection Response in Replicated Drosophila melanogaster Populations with Reduced Genetic Variation. Genome Biol Evol 2021; 13:6409861. [PMID: 34694407 PMCID: PMC8599828 DOI: 10.1093/gbe/evab239] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/21/2021] [Indexed: 12/12/2022] Open
Abstract
Many adaptive traits are polygenic and frequently more loci contributing to the phenotype are segregating than needed to express the phenotypic optimum. Experimental evolution with replicated populations adapting to a new controlled environment provides a powerful approach to study polygenic adaptation. Because genetic redundancy often results in nonparallel selection responses among replicates, we propose a modified evolve and resequence (E&R) design that maximizes the similarity among replicates. Rather than starting from many founders, we only use two inbred Drosophila melanogaster strains and expose them to a very extreme, hot temperature environment (29 °C). After 20 generations, we detect many genomic regions with a strong, highly parallel selection response in 10 evolved replicates. The X chromosome has a more pronounced selection response than the autosomes, which may be attributed to dominance effects. Furthermore, we find that the median selection coefficient for all chromosomes is higher in our two-genotype experiment than in classic E&R studies. Because two random genomes harbor sufficient variation for adaptive responses, we propose that this approach is particularly well-suited for the analysis of polygenic adaptation.
Collapse
Affiliation(s)
- Claire Burny
- Institut für Populationsgenetik, Vetmeduni Vienna, Austria.,Vienna Graduate School of Population Genetics, Vetmeduni Vienna, Austria
| | - Viola Nolte
- Institut für Populationsgenetik, Vetmeduni Vienna, Austria
| | - Marlies Dolezal
- Plattform Bioinformatik und Biostatistik, Vetmeduni Vienna, Wien, Austria
| | | |
Collapse
|
37
|
Acosta K, Appenroth KJ, Borisjuk L, Edelman M, Heinig U, Jansen MAK, Oyama T, Pasaribu B, Schubert I, Sorrels S, Sree KS, Xu S, Michael TP, Lam E. Return of the Lemnaceae: duckweed as a model plant system in the genomics and postgenomics era. THE PLANT CELL 2021; 33:3207-3234. [PMID: 34273173 PMCID: PMC8505876 DOI: 10.1093/plcell/koab189] [Citation(s) in RCA: 76] [Impact Index Per Article: 25.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 06/18/2021] [Indexed: 05/05/2023]
Abstract
The aquatic Lemnaceae family, commonly called duckweed, comprises some of the smallest and fastest growing angiosperms known on Earth. Their tiny size, rapid growth by clonal propagation, and facile uptake of labeled compounds from the media were attractive features that made them a well-known model for plant biology from 1950 to 1990. Interest in duckweed has steadily regained momentum over the past decade, driven in part by the growing need to identify alternative plants from traditional agricultural crops that can help tackle urgent societal challenges, such as climate change and rapid population expansion. Propelled by rapid advances in genomic technologies, recent studies with duckweed again highlight the potential of these small plants to enable discoveries in diverse fields from ecology to chronobiology. Building on established community resources, duckweed is reemerging as a platform to study plant processes at the systems level and to translate knowledge gained for field deployment to address some of society's pressing needs. This review details the anatomy, development, physiology, and molecular characteristics of the Lemnaceae to introduce them to the broader plant research community. We highlight recent research enabled by Lemnaceae to demonstrate how these plants can be used for quantitative studies of complex processes and for revealing potentially novel strategies in plant defense and genome maintenance.
Collapse
Affiliation(s)
- Kenneth Acosta
- Department of Plant Biology, Rutgers the State University of New Jersey, New Brunswick, NJ 08901, USA
| | - Klaus J Appenroth
- Plant Physiology, Matthias Schleiden Institute, University of Jena, Jena 07737, Germany
| | - Ljudmilla Borisjuk
- The Leibniz Institute of Plant Genetics and Crop Plant Research, Gatersleben D-06466, Germany
| | - Marvin Edelman
- Department of Plant and Environmental Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Uwe Heinig
- Department of Plant and Environmental Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Marcel A K Jansen
- School of Biological, Earth and Environmental Sciences, Environmental Research Institute, University College Cork, Cork T23 TK30, Ireland
| | - Tokitaka Oyama
- Department of Botany, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Buntora Pasaribu
- Department of Plant Biology, Rutgers the State University of New Jersey, New Brunswick, NJ 08901, USA
| | - Ingo Schubert
- The Leibniz Institute of Plant Genetics and Crop Plant Research, Gatersleben D-06466, Germany
| | - Shawn Sorrels
- Department of Plant Biology, Rutgers the State University of New Jersey, New Brunswick, NJ 08901, USA
| | - K Sowjanya Sree
- Department of Environmental Science, Central University of Kerala, Periye 671320, India
| | - Shuqing Xu
- Institute for Evolution and Biodiversity, University of Münster, Münster 48149, Germany
| | | | - Eric Lam
- Author for correspondence: (E.L.), (T.P.M.)
| |
Collapse
|
38
|
Abstract
Aging has provided fruitful challenges for evolutionary theory, and evolutionary theory has deepened our understanding of aging. A great deal of genetic and molecular data now exists concerning mortality regulation and there is a growing body of knowledge concerning the life histories of diverse species. Assimilating all relevant data into a framework for the evolution of aging promises to significantly advance the field. We propose extensions of some key concepts to provide greater precision when applying these concepts to age-structured contexts. Secondary or byproduct effects of mutations are proposed as an important factor affecting survival patterns, including effects that may operate in small populations subject to genetic drift, widening the possibilities for mutation accumulation and pleiotropy. Molecular and genetic studies have indicated a diverse array of mechanisms that can modify aging and mortality rates, while transcriptome data indicate a high level of tissue and species specificity for genes affected by aging. The diversity of mechanisms and gene effects that can contribute to the pattern of aging in different organisms may mirror the complex evolutionary processes behind aging.
Collapse
Affiliation(s)
- Stewart Frankel
- Biology Department, University of Hartford, West Hartford, CT, United States
| | - Blanka Rogina
- Genetics and Genome Sciences, Institute for Systems Genomics, School of Medicine, University of Connecticut Health Center, Farmington, CT, United States
| |
Collapse
|
39
|
Shahrestani P, King E, Ramezan R, Phillips M, Riddle M, Thornburg M, Greenspan Z, Estrella Y, Garcia K, Chowdhury P, Malarat G, Zhu M, Rottshaefer SM, Wraight S, Griggs M, Vandenberg J, Long AD, Clark AG, Lazzaro BP. The molecular architecture of Drosophila melanogaster defense against Beauveria bassiana explored through evolve and resequence and quantitative trait locus mapping. G3-GENES GENOMES GENETICS 2021; 11:6371870. [PMID: 34534291 PMCID: PMC8664422 DOI: 10.1093/g3journal/jkab324] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 08/17/2021] [Indexed: 12/02/2022]
Abstract
Little is known about the genetic architecture of antifungal immunity in natural populations. Using two population genetic approaches, quantitative trait locus (QTL) mapping and evolve and resequence (E&R), we explored D. melanogaster immune defense against infection with the fungus Beauveria bassiana. The immune defense was highly variable both in the recombinant inbred lines from the Drosophila Synthetic Population Resource used for our QTL mapping and in the synthetic outbred populations used in our E&R study. Survivorship of infection improved dramatically over just 10 generations in the E&R study, and continued to increase for an additional nine generations, revealing a trade-off with uninfected longevity. Populations selected for increased defense against B. bassiana evolved cross resistance to a second, distinct B. bassiana strain but not to bacterial pathogens. The QTL mapping study revealed that sexual dimorphism in defense depends on host genotype, and the E&R study indicated that sexual dimorphism also depends on the specific pathogen to which the host is exposed. Both the QTL mapping and E&R experiments generated lists of potentially causal candidate genes, although these lists were nonoverlapping.
Collapse
Affiliation(s)
- Parvin Shahrestani
- Department of Biological Science, California State University Fullerton, Fullerton CA, 92831, USA
| | - Elizabeth King
- Division of Biological Sciences, University of Missouri, Columbia MO, 65211, USA
| | - Reza Ramezan
- Department of Statistics and Actuarial Science, University of Waterloo, Waterloo ON, N2L 3G1, Canada
| | - Mark Phillips
- Department of Integrative Biology, Oregon State University, Corvallis OR, 97331, USA
| | - Melissa Riddle
- Department of Biological Science, California State University Fullerton, Fullerton CA, 92831, USA
| | - Marisa Thornburg
- Department of Biological Science, California State University Fullerton, Fullerton CA, 92831, USA
| | - Zachary Greenspan
- Department of Ecology and Evolutionary Biology, University of California Irvine, Irvine CA, 92692, USA
| | | | - Kelly Garcia
- Department of Entomology, Cornell University, Ithaca NY, 14853, USA
| | - Pratik Chowdhury
- Department of Entomology, Cornell University, Ithaca NY, 14853, USA
| | - Glen Malarat
- Department of Entomology, Cornell University, Ithaca NY, 14853, USA
| | - Ming Zhu
- Department of Entomology, Cornell University, Ithaca NY, 14853, USA
| | | | - Stephen Wraight
- USDA ARS Emerging Pets and Pathogens Research Unit, Robert W. Holley Center for Agriculture & Health, Ithaca NY, 14853, USA
| | - Michael Griggs
- USDA ARS Emerging Pets and Pathogens Research Unit, Robert W. Holley Center for Agriculture & Health, Ithaca NY, 14853, USA
| | - John Vandenberg
- USDA ARS Emerging Pets and Pathogens Research Unit, Robert W. Holley Center for Agriculture & Health, Ithaca NY, 14853, USA
| | - Anthony D Long
- Department of Ecology and Evolutionary Biology, University of California Irvine, Irvine CA, 92692, USA
| | - Andrew G Clark
- Department of Molecular Biology and Genetics, Cornell University, Ithaca NY, 14853, USA
| | - Brian P Lazzaro
- Department of Entomology, Cornell University, Ithaca NY, 14853, USA
| |
Collapse
|
40
|
Croze M, Kim Y. Inference of population genetic parameters from an irregular time series of seasonal influenza virus sequences. Genetics 2021; 217:6066165. [PMID: 33724414 DOI: 10.1093/genetics/iyaa039] [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: 11/22/2020] [Accepted: 12/17/2020] [Indexed: 11/12/2022] Open
Abstract
Basic summary statistics that quantify the population genetic structure of influenza virus are important for understanding and inferring the evolutionary and epidemiological processes. However, the sampling dates of global virus sequences in the last several decades are scattered nonuniformly throughout the calendar. Such temporal structure of samples and the small effective size of viral population hampers the use of conventional methods to calculate summary statistics. Here, we define statistics that overcome this problem by correcting for the sampling-time difference in quantifying a pairwise sequence difference. A simple linear regression method jointly estimates the mutation rate and the level of sequence polymorphism, thus providing an estimate of the effective population size. It also leads to the definition of Wright's FST for arbitrary time-series data. Furthermore, as an alternative to Tajima's D statistic or the site-frequency spectrum, a mismatch distribution corrected for sampling-time differences can be obtained and compared between actual and simulated data. Application of these methods to seasonal influenza A/H3N2 viruses sampled between 1980 and 2017 and sequences simulated under the model of recurrent positive selection with metapopulation dynamics allowed us to estimate the synonymous mutation rate and find parameter values for selection and demographic structure that fit the observation. We found that the mutation rates of HA and PB1 segments before 2007 were particularly high and that including recurrent positive selection in our model was essential for the genealogical structure of the HA segment. Methods developed here can be generally applied to population genetic inferences using serially sampled genetic data.
Collapse
Affiliation(s)
- Myriam Croze
- Division of EcoScience, Ewha Womans University, Seoul 03760, Korea
| | - Yuseob Kim
- Division of EcoScience, Ewha Womans University, Seoul 03760, Korea.,Department of Life Science, Ewha Womans University, Seoul 03760, Korea
| |
Collapse
|
41
|
Henry LP, Bruijning M, Forsberg SKG, Ayroles JF. The microbiome extends host evolutionary potential. Nat Commun 2021; 12:5141. [PMID: 34446709 PMCID: PMC8390463 DOI: 10.1038/s41467-021-25315-x] [Citation(s) in RCA: 104] [Impact Index Per Article: 34.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 08/03/2021] [Indexed: 02/07/2023] Open
Abstract
The microbiome shapes many host traits, yet the biology of microbiomes challenges traditional evolutionary models. Here, we illustrate how integrating the microbiome into quantitative genetics can help untangle complexities of host-microbiome evolution. We describe two general ways in which the microbiome may affect host evolutionary potential: by shifting the mean host phenotype and by changing the variance in host phenotype in the population. We synthesize the literature across diverse taxa and discuss how these scenarios could shape the host response to selection. We conclude by outlining key avenues of research to improve our understanding of the complex interplay between hosts and microbiomes.
Collapse
Affiliation(s)
- Lucas P. Henry
- grid.16750.350000 0001 2097 5006Dept. of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ USA ,grid.16750.350000 0001 2097 5006Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ USA
| | - Marjolein Bruijning
- grid.16750.350000 0001 2097 5006Dept. of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ USA
| | - Simon K. G. Forsberg
- grid.16750.350000 0001 2097 5006Dept. of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ USA ,grid.16750.350000 0001 2097 5006Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ USA ,grid.8993.b0000 0004 1936 9457Dept. of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Julien F. Ayroles
- grid.16750.350000 0001 2097 5006Dept. of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ USA ,grid.16750.350000 0001 2097 5006Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ USA
| |
Collapse
|
42
|
Ulrich GF, Zemp N, Vorburger C, Boulain H. Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance. Heredity (Edinb) 2021; 127:219-232. [PMID: 34012059 PMCID: PMC8322320 DOI: 10.1038/s41437-021-00444-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 04/26/2021] [Accepted: 04/27/2021] [Indexed: 02/04/2023] Open
Abstract
Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.
Collapse
Affiliation(s)
- Gabriel F. Ulrich
- grid.418656.80000 0001 1551 0562EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Überlandstrasse 133, 8600 Dübendorf, Switzerland ,grid.5801.c0000 0001 2156 2780Institute of Integrative Biology, ETH Zürich, Universitätsstrasse 16, 8092 Zürich, Switzerland
| | - Niklaus Zemp
- Genetic Diversity Centre, Department of Environmental Systems Sciences, ETH Zürich, 8092 Zürich, Switzerland
| | - Christoph Vorburger
- grid.418656.80000 0001 1551 0562EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Überlandstrasse 133, 8600 Dübendorf, Switzerland ,grid.5801.c0000 0001 2156 2780Institute of Integrative Biology, ETH Zürich, Universitätsstrasse 16, 8092 Zürich, Switzerland
| | - Hélène Boulain
- grid.418656.80000 0001 1551 0562EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Überlandstrasse 133, 8600 Dübendorf, Switzerland ,grid.9851.50000 0001 2165 4204Present Address: Department of Ecology and Evolution, University of Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
43
|
Janzen T, Diaz F. Individual‐based simulations of genome evolution with ancestry: The
GenomeAdmixR
R package. Methods Ecol Evol 2021. [DOI: 10.1111/2041-210x.13612] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Thijs Janzen
- Groningen Institute for Evolutionary Life Sciences University of Groningen Groningen The Netherlands
- Carl von Ossietzky University Oldenburg Germany
| | - Fernando Diaz
- Department of Entomology University of Arizona Tucson AZ USA
| |
Collapse
|
44
|
White NJ, Butlin RK. Multidimensional divergent selection, local adaptation, and speciation. Evolution 2021; 75:2167-2178. [PMID: 34263939 DOI: 10.1111/evo.14312] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Revised: 06/29/2021] [Accepted: 07/05/2021] [Indexed: 12/24/2022]
Abstract
Divergent selection applied to one or more traits drives local adaptation and may lead to ecological speciation. Divergent selection on many traits might be termed "multidimensional" divergent selection. There is a commonly held view that multidimensional divergent selection is likely to promote local adaptation and speciation to a greater extent than unidimensional divergent selection. We disentangle the core concepts underlying dimensionality as a property of the environment, phenotypes, and genome. In particular, we identify a need to separate the overall strength of selection and the number of loci affected from dimensionality per se, and to distinguish divergence dimensionality from dimensionality of stabilizing selection. We then critically scrutinize this commonly held view that multidimensional selection promotes speciation, re-examining the evidence base from theory, experiments, and nature. We conclude that the evidence base is currently weak and generally suffers from confounding of possible causal effects. Finally, we propose several mechanisms by which multidimensional divergent selection and related processes might influence divergence, both as a driver and as a barrier.
Collapse
Affiliation(s)
- Nathan J White
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, S10 2TN, United Kingdom
| | - Roger K Butlin
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, S10 2TN, United Kingdom.,Department of Marine Sciences, University of Gothenburg, Gothenburg, SE-40530, Sweden
| |
Collapse
|
45
|
Bitner K, Rutledge GA, Kezos JN, Mueller LD. The effects of adaptation to urea on feeding rates and growth in Drosophila larvae. Ecol Evol 2021; 11:9516-9529. [PMID: 34306639 PMCID: PMC8293711 DOI: 10.1002/ece3.7770] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2020] [Revised: 04/22/2021] [Accepted: 05/21/2021] [Indexed: 11/18/2022] Open
Abstract
A collection of forty populations were used to study the phenotypic adaptation of Drosophila melanogaster larvae to urea-laced food. A long-term goal of this research is to map genes responsible for these phenotypes. This mapping requires large numbers of populations. Thus, we studied fifteen populations subjected to direct selection for urea tolerance and five controls. In addition, we studied another twenty populations which had not been exposed to urea but were subjected to stress or demographic selection. In this study, we describe the differentiation in these population for six phenotypes: (1) larval feeding rates, (2) larval viability in urea-laced food, (3) larval development time in urea-laced food, (4) adult starvation times, (5) adult desiccation times, and (6) larval growth rates. No significant differences were observed for desiccation resistance. The demographically/stress-selected populations had longer times to starvation than urea-selected populations. The urea-adapted populations showed elevated survival and reduced development time in urea-laced food relative to the control and nonadapted populations. The urea-adapted populations also showed reduced larval feeding rates relative to controls. We show that there is a strong linear relationship between feeding rates and growth rates at the same larval ages feeding rates were measured. This suggests that feeding rates are correlated with food intake and growth. This relationship between larval feeding rates, food consumption, and efficiency has been postulated to involve important trade-offs that govern larval evolution in stressful environments. Our results support the idea that energy allocation is a central organizing theme in adaptive evolution.
Collapse
Affiliation(s)
- Kathreen Bitner
- Department of Ecology and Evolutionary BiologyUniversity of California, IrvineIrvineCAUSA
| | - Grant A. Rutledge
- Department of Ecology and Evolutionary BiologyUniversity of California, IrvineIrvineCAUSA
- USDA HNRCA at Tufts UniversityBostonMAUSA
| | - James N. Kezos
- Department of Ecology and Evolutionary BiologyUniversity of California, IrvineIrvineCAUSA
- Department of Development, Aging, and RegenerationSanford Burnham Prebys Medical Discovery InstituteLa JollaCAUSA
| | - Laurence D. Mueller
- Department of Ecology and Evolutionary BiologyUniversity of California, IrvineIrvineCAUSA
| |
Collapse
|
46
|
Langmüller AM, Dolezal M, Schlötterer C. Fine Mapping without Phenotyping: Identification of Selection Targets in Secondary Evolve and Resequence Experiments. Genome Biol Evol 2021; 13:6311659. [PMID: 34190980 PMCID: PMC8358229 DOI: 10.1093/gbe/evab154] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/24/2021] [Indexed: 12/19/2022] Open
Abstract
Evolve and Resequence (E&R) studies investigate the genomic selection response of populations in an Experimental Evolution setup. Despite the popularity of E&R, empirical studies in sexually reproducing organisms typically suffer from an excess of candidate loci due to linkage disequilibrium, and single gene or SNP resolution is the exception rather than the rule. Recently, so-called "secondary E&R" has been suggested as promising experimental follow-up procedure to confirm putatively selected regions from a primary E&R study. Secondary E&R provides also the opportunity to increase mapping resolution by allowing for additional recombination events, which separate the selection target from neutral hitchhikers. Here, we use computer simulations to assess the effect of different crossing schemes, population size, experimental duration, and number of replicates on the power and resolution of secondary E&R. We find that the crossing scheme and population size are crucial factors determining power and resolution of secondary E&R: A simple crossing scheme with few founder lines consistently outcompetes crossing schemes where evolved populations from a primary E&R experiment are mixed with a complex ancestral founder population. Regardless of the experimental design tested, a population size of at least 4,800 individuals, which is roughly five times larger than population sizes in typical E&R studies, is required to achieve a power of at least 75%. Our study provides an important step toward improved experimental designs aiming to characterize causative SNPs in Experimental Evolution studies.
Collapse
Affiliation(s)
- Anna Maria Langmüller
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria.,Vienna Graduate School of Population Genetics, Vetmeduni Vienna, Vienna, Austria
| | - Marlies Dolezal
- Plattform Bioinformatik und Biostatistik, Vetmeduni Vienna, Vienna, Austria
| | | |
Collapse
|
47
|
Pointer MD, Gage MJG, Spurgin LG. Tribolium beetles as a model system in evolution and ecology. Heredity (Edinb) 2021; 126:869-883. [PMID: 33767370 PMCID: PMC8178323 DOI: 10.1038/s41437-021-00420-1] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Revised: 02/19/2021] [Accepted: 02/19/2021] [Indexed: 02/01/2023] Open
Abstract
Flour beetles of the genus Tribolium have been utilised as informative study systems for over a century and contributed to major advances across many fields. This review serves to highlight the significant historical contribution that Tribolium study systems have made to the fields of ecology and evolution, and to promote their use as contemporary research models. We review the broad range of studies employing Tribolium to make significant advances in ecology and evolution. We show that research using Tribolium beetles has contributed a substantial amount to evolutionary and ecological understanding, especially in the fields of population dynamics, reproduction and sexual selection, population and quantitative genetics, and behaviour, physiology and life history. We propose a number of future research opportunities using Tribolium, with particular focus on how their amenability to forward and reverse genetic manipulation may provide a valuable complement to other insect models.
Collapse
Affiliation(s)
- Michael D Pointer
- School of Biological Sciences, University of East Anglia, Norwich, UK.
| | - Matthew J G Gage
- School of Biological Sciences, University of East Anglia, Norwich, UK
| | - Lewis G Spurgin
- School of Biological Sciences, University of East Anglia, Norwich, UK.
| |
Collapse
|
48
|
Kim NY, Lee HJ, Kim HS, Lee SH, Moon JS, Jeong RD. Identification of Plant Viruses Infecting Pear Using RNA Sequencing. THE PLANT PATHOLOGY JOURNAL 2021; 37:258-267. [PMID: 34111915 PMCID: PMC8200581 DOI: 10.5423/ppj.oa.01.2021.0009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/25/2021] [Revised: 03/30/2021] [Accepted: 04/23/2021] [Indexed: 06/12/2023]
Abstract
Asian pear (Pyrus pyrifolia) is a widely cultivated and commercially important fruit crop, which is occasionally subject to severe economic losses due to latent viral infections. Thus, the aim of the present study was to examine and provide a comprehensive overview of virus populations infecting a major pear cultivar ('Singo') in Korea. From June 2017 to October 2019, leaf samples (n = 110) of pear trees from 35 orchards in five major pear-producing regions were collected and subjected to RNA sequencing. Most virus-associated contigs matched the sequences of known viruses, including apple stem grooving virus (ASGV) and apple stem pitting virus (ASPV). However, some contigs matched the sequences of apple green crinkle-associated virus and cucumber mosaic virus. In addition, three complete or nearly complete genomes were constructed based on transcriptome data and subjected to phylogenetic analyses. Based on the number of virus-associated reads, ASGV and ASPV were identified as the dominant viruses of 'Singo.' The present study describes the virome of a major pear cultivar in Korea, and looks into the diversity of viral communities in this cultivar. This study can provide valuable information on the complexity of genetic variability of viruses infecting pear trees.
Collapse
Affiliation(s)
- Nam-Yeon Kim
- Department of Applied Biology, Institute of Environmentally Friendly Agriculture, Chonnam National University, Gwangju 61185, Korea
| | - Hyo-Jeong Lee
- Department of Applied Biology, Institute of Environmentally Friendly Agriculture, Chonnam National University, Gwangju 61185, Korea
| | - Hong-Sup Kim
- Seed Testing & Research Center, Korea Seed & Variety Service, Gimcheon 39660, Korea
| | - Su-Heon Lee
- School of Applied Biosciences, Kyungpook National University, Daegu 98411, Korea
| | - Jae-Sun Moon
- Plant Genome Research Center, Korea Research Institute of Biosciences & Biotechnology, Daejeon 34141, Korea
| | - Rae-Dong Jeong
- Department of Applied Biology, Institute of Environmentally Friendly Agriculture, Chonnam National University, Gwangju 61185, Korea
| |
Collapse
|
49
|
North HL, McGaughran A, Jiggins CD. Insights into invasive species from whole-genome resequencing. Mol Ecol 2021; 30:6289-6308. [PMID: 34041794 DOI: 10.1111/mec.15999] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Revised: 03/12/2021] [Accepted: 04/30/2021] [Indexed: 12/12/2022]
Abstract
Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.
Collapse
Affiliation(s)
- Henry L North
- Department of Zoology, University of Cambridge, Cambridge, UK
| | - Angela McGaughran
- Te Aka Mātuatua/School of Science, University of Waikato, Hamilton, New Zealand
| | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge, UK
| |
Collapse
|
50
|
Pereira RJ, Lima TG, Pierce-Ward NT, Chao L, Burton RS. Recovery from hybrid breakdown reveals a complex genetic architecture of mitonuclear incompatibilities. Mol Ecol 2021; 30:6403-6416. [PMID: 34003535 DOI: 10.1111/mec.15985] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2020] [Revised: 03/29/2021] [Accepted: 05/11/2021] [Indexed: 01/03/2023]
Abstract
Reproductive isolation is often achieved when genes that are neutral or beneficial in their genomic background become functionally incompatible in a foreign genomic background, causing inviability, sterility or other forms of low fitness in hybrids. Recent studies suggest that mitonuclear interactions are among the initial incompatibilities to evolve at early stages of population divergence across taxa. Yet, the genomic architecture of mitonuclear incompatibilities has rarely been elucidated. We employ an experimental evolution approach starting with low-fitness F2 interpopulation hybrids of the copepod Tigriopus californicus, in which frequencies of compatible and incompatible nuclear alleles change in response to an alternative mitochondrial background. After about nine generations, we observe a generalized increase in population size and in survivorship, suggesting efficiency of selection against maladaptive phenotypes. Whole genome sequencing of evolved populations showed some consistent allele frequency changes across three replicates of each reciprocal cross, but markedly different patterns between mitochondrial backgrounds. In only a few regions (~6.5% of the genome), the same parental allele was overrepresented irrespective of the mitochondrial background. About 33% of the genome showed allele frequency changes consistent with divergent selection, with the location of these genomic regions strongly differing between mitochondrial backgrounds. In 87% and 89% of these genomic regions, the dominant nuclear allele matched the associated mitochondrial background, consistent with mitonuclear co-adaptation. These results suggest that mitonuclear incompatibilities have a complex polygenic architecture that differs between populations, potentially generating genome-wide barriers to gene flow between closely related taxa.
Collapse
Affiliation(s)
- Ricardo J Pereira
- Division of Evolutionary Biology, Faculty of Biology, Ludwig-Maximilians-Universität München, Munich, Germany
| | - Thiago G Lima
- Marine Biology Research Division, Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA, USA
| | - N Tessa Pierce-Ward
- Marine Biology Research Division, Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA, USA
| | - Lin Chao
- Division of Biological Sciences, University of California San Diego, La Jolla, CA, USA
| | - Ronald S Burton
- Marine Biology Research Division, Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA, USA
| |
Collapse
|