1
|
Lievens EJP, Kühn S, Horas EL, Le Pennec G, Peter S, Petrosky AD, Künzel S, Feulner PGD, Becks L. High parasite diversity maintained after an alga-virus coevolutionary arms race. J Evol Biol 2024; 37:795-806. [PMID: 38699979 DOI: 10.1093/jeb/voae053] [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: 10/01/2023] [Revised: 02/27/2024] [Accepted: 05/02/2024] [Indexed: 05/05/2024]
Abstract
Arms race dynamics are a common outcome of host-parasite coevolution. While they can theoretically be maintained indefinitely, realistic arms races are expected to be finite. Once an arms race has ended, for example due to the evolution of a generalist-resistant host, the system may transition into coevolutionary dynamics that favour long-term diversity. In microbial experiments, host-parasite arms races often transition into a stable coexistence of generalist-resistant hosts, (semi-)susceptible hosts, and parasites. While long-term host diversity is implicit in these cases, parasite diversity is usually overlooked. In this study, we examined parasite diversity after the end of an experimental arms race between a unicellular alga (Chlorella variabilis) and its lytic virus (PBCV-1). First, we isolated virus genotypes from multiple time points from two replicate microcosms. A time-shift experiment confirmed that the virus isolates had escalating host ranges, i.e., that arms races had occurred. We then examined the phenotypic and genetic diversity of virus isolates from the post-arms race phase. Post-arms race virus isolates had diverse host ranges, survival probabilities, and growth rates; they also clustered into distinct genetic groups. Importantly, host range diversity was maintained throughout the post-arms race phase, and the frequency of host range phenotypes fluctuated over time. We hypothesize that this dynamic polymorphism was maintained by a combination of fluctuating selection and demographic stochasticity. Together with previous work in prokaryotic systems, our results link experimental observations of arms races to natural observations of long-term host and parasite diversity.
Collapse
Affiliation(s)
- Eva J P Lievens
- Aquatic Ecology and Evolution Group, Department of Biology, University of Konstanz, Konstanz, Germany
| | - Samuel Kühn
- Aquatic Ecology and Evolution Group, Department of Biology, University of Konstanz, Konstanz, Germany
| | - Elena L Horas
- Aquatic Ecology and Evolution Group, Department of Biology, University of Konstanz, Konstanz, Germany
| | - Guénolé Le Pennec
- Department of Fish Ecology and Evolution, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
| | - Sarah Peter
- Aquatic Ecology and Evolution Group, Department of Biology, University of Konstanz, Konstanz, Germany
| | - Azade D Petrosky
- Aquatic Ecology and Evolution Group, Department of Biology, University of Konstanz, Konstanz, Germany
| | - Sven Künzel
- Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Philine G D Feulner
- Department of Fish Ecology and Evolution, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
- Department of Biology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Lutz Becks
- Aquatic Ecology and Evolution Group, Department of Biology, University of Konstanz, Konstanz, Germany
| |
Collapse
|
2
|
Soni V, Jensen JD. Temporal challenges in detecting balancing selection from population genomic data. G3 (BETHESDA, MD.) 2024; 14:jkae069. [PMID: 38551137 DOI: 10.1093/g3journal/jkae069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 12/21/2023] [Accepted: 03/19/2024] [Indexed: 04/28/2024]
Abstract
The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).
Collapse
Affiliation(s)
- Vivak Soni
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, AZ 85281, USA
| | - Jeffrey D Jensen
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, AZ 85281, USA
| |
Collapse
|
3
|
Korfmann K, Temple-Boyer M, Sellinger T, Tellier A. Determinants of rapid adaptation in species with large variance in offspring production. Mol Ecol 2024; 33:e16982. [PMID: 37199145 DOI: 10.1111/mec.16982] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Revised: 04/26/2023] [Accepted: 05/02/2023] [Indexed: 05/19/2023]
Abstract
The speed of population adaptation to changing biotic and abiotic environments is determined by the interaction between genetic drift, positive selection and linkage effects. Many marine species (fish, crustaceans), invertebrates and pathogens of humans and crops, exhibit sweepstakes reproduction characterized by the production of a very large amount of offspring (fecundity phase) from which only a small fraction may survive to the next generation (viability phase). Using stochastic simulations, we investigate whether the occurrence of sweepstakes reproduction affects the efficiency of a positively selected unlinked locus, and thus, the speed of adaptation since fecundity and/or viability have distinguishable consequences on mutation rate, probability and fixation time of advantageous alleles. We observe that the mean number of mutations at the next generation is always the function of the population size, but the variance increases with stronger sweepstakes reproduction when mutations occur in the parents. On the one hand, stronger sweepstakes reproduction magnifies the effect of genetic drift thus increasing the probability of fixation of neutral allele and decreasing that of selected alleles. On the other hand, the time to fixation of advantageous (as well as neutral) alleles is shortened by stronger sweepstakes reproduction. Importantly, fecundity and viability selection exhibit different probabilities and times to fixation of advantageous alleles under intermediate and weak sweepstakes reproduction. Finally, alleles under both strong fecundity and viability selection display a synergistic efficiency of selection. We conclude that measuring and modelling accurately fecundity and/or viability selection are crucial to predict the adaptive potential of species with sweepstakes reproduction.
Collapse
Affiliation(s)
- Kevin Korfmann
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
| | - Marie Temple-Boyer
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
| | - Thibaut Sellinger
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
- Department of Environment and Biodiversity, Paris Lodron University of Salzburg, Salzburg, Austria
| | - Aurélien Tellier
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, Freising, Germany
| |
Collapse
|
4
|
Eldon B, Stephan W. Sweepstakes reproduction facilitates rapid adaptation in highly fecund populations. Mol Ecol 2024; 33:e16903. [PMID: 36896794 DOI: 10.1111/mec.16903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Revised: 02/21/2023] [Accepted: 02/23/2023] [Indexed: 03/11/2023]
Abstract
Adaptation enables natural populations to survive in a changing environment. Understanding the mechanics of adaptation is therefore crucial for learning about the evolution and ecology of natural populations. We focus on the impact of random sweepstakes on selection in highly fecund haploid and diploid populations partitioned into two genetic types, with one type conferring selective advantage. For the diploid populations, we incorporate various dominance mechanisms. We assume that the populations may experience recurrent bottlenecks. In random sweepstakes, the distribution of individual recruitment success is highly skewed, resulting in a huge variance in the number of offspring contributed by the individuals present in any given generation. Using computer simulations, we investigate the joint effects of random sweepstakes, recurrent bottlenecks and dominance mechanisms on selection. In our framework, bottlenecks allow random sweepstakes to have an effect on the time to fixation, and in diploid populations, the effect of random sweepstakes depends on the dominance mechanism. We describe selective sweepstakes that are approximated by recurrent sweeps of strongly beneficial allelic types arising by mutation. We demonstrate that both types of sweepstakes reproduction may facilitate rapid adaptation (as defined based on the average time to fixation of a type conferring selective advantage conditioned on fixation of the type). However, whether random sweepstakes cause rapid adaptation depends also on their interactions with bottlenecks and dominance mechanisms. Finally, we review a case study in which a model of recurrent sweeps is shown to essentially explain population genomic data from Atlantic cod.
Collapse
Affiliation(s)
- Bjarki Eldon
- Institute of Evolution and Biodiversity Science, Natural History Museum Berlin, Berlin, Germany
| | | |
Collapse
|
5
|
Soni V, Terbot JW, Jensen JD. Population genetic considerations regarding the interpretation of within-patient SARS-CoV-2 polymorphism data. Nat Commun 2024; 15:3240. [PMID: 38627371 PMCID: PMC11021480 DOI: 10.1038/s41467-024-46261-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Accepted: 01/29/2024] [Indexed: 04/19/2024] Open
Affiliation(s)
- Vivak Soni
- Center for Evolution & Medicine, Arizona State University, School of Life Sciences, Tempe, AZ, USA
| | - John W Terbot
- Center for Evolution & Medicine, Arizona State University, School of Life Sciences, Tempe, AZ, USA
- Division of Biological Sciences, University of Montana, Missoula, MT, USA
| | - Jeffrey D Jensen
- Center for Evolution & Medicine, Arizona State University, School of Life Sciences, Tempe, AZ, USA.
| |
Collapse
|
6
|
Soni V, Pfeifer SP, Jensen JD. The Effects of Mutation and Recombination Rate Heterogeneity on the Inference of Demography and the Distribution of Fitness Effects. Genome Biol Evol 2024; 16:evae004. [PMID: 38207127 PMCID: PMC10834165 DOI: 10.1093/gbe/evae004] [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: 08/31/2023] [Revised: 12/12/2023] [Accepted: 01/07/2024] [Indexed: 01/13/2024] Open
Abstract
Disentangling the effects of demography and selection has remained a focal point of population genetic analysis. Knowledge about mutation and recombination is essential in this endeavor; however, despite clear evidence that both mutation and recombination rates vary across genomes, it is common practice to model both rates as fixed. In this study, we quantify how this unaccounted for rate heterogeneity may impact inference using common approaches for inferring selection (DFE-alpha, Grapes, and polyDFE) and/or demography (fastsimcoal2 and δaδi). We demonstrate that, if not properly modeled, this heterogeneity can increase uncertainty in the estimation of demographic and selective parameters and in some scenarios may result in mis-leading inference. These results highlight the importance of quantifying the fundamental evolutionary parameters of mutation and recombination before utilizing population genomic data to quantify the effects of genetic drift (i.e. as modulated by demographic history) and selection; or, at the least, that the effects of uncertainty in these parameters can and should be directly modeled in downstream inference.
Collapse
Affiliation(s)
- Vivak Soni
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, AZ, USA
| | - Susanne P Pfeifer
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, AZ, USA
| | - Jeffrey D Jensen
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, AZ, USA
| |
Collapse
|
7
|
Terbot JW, Cooper BS, Good JM, Jensen JD. A Simulation Framework for Modeling the Within-Patient Evolutionary Dynamics of SARS-CoV-2. Genome Biol Evol 2023; 15:evad204. [PMID: 37950882 PMCID: PMC10664409 DOI: 10.1093/gbe/evad204] [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/2023] [Revised: 10/31/2023] [Accepted: 11/07/2023] [Indexed: 11/13/2023] Open
Abstract
The global impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for rarely acting positive selection are best performed via comparison of empirical data with simulated data wherein commonly acting evolutionary factors, including mutation and recombination, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. Although there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intrahost evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them with existing empirical data. Of these, 592 models (∼5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intrahost SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed toward strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.
Collapse
Affiliation(s)
- John W Terbot
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, Arizona, USA
- Division of Biological Sciences, University of Montana, Missoula, Montana, USA
| | - Brandon S Cooper
- Division of Biological Sciences, University of Montana, Missoula, Montana, USA
| | - Jeffrey M Good
- Division of Biological Sciences, University of Montana, Missoula, Montana, USA
| | - Jeffrey D Jensen
- School of Life Sciences, Center for Evolution & Medicine, Arizona State University, Tempe, Arizona, USA
| |
Collapse
|
8
|
Sexton NR, Cline PJ, Gallichotte EN, Fitzmeyer E, Young MC, Janich AJ, Pabilonia KL, Ehrhart N, Ebel GD. SARS-CoV-2 entry into and evolution within a skilled nursing facility. Sci Rep 2023; 13:11657. [PMID: 37468595 DOI: 10.1038/s41598-023-38544-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 07/10/2023] [Indexed: 07/21/2023] Open
Abstract
SARS-CoV-2 belongs to the family Coronaviridae which includes multiple human pathogens that have an outsized impact on aging populations. As a novel human pathogen, SARS-CoV-2 is undergoing continuous adaptation to this new host species and there is evidence of this throughout the scientific and public literature. However, most investigations of SARS-CoV-2 evolution have focused on large-scale collections of data across diverse populations and/or living environments. Here we investigate SARS-CoV-2 evolution in epidemiologically linked individuals within a single outbreak at a skilled nursing facility beginning with initial introduction of the pathogen. The data demonstrate that SARS-CoV-2 was introduced to the facility multiple times without establishing an interfacility transmission chain, followed by a single introduction that infected many individuals within a week. This large-scale introduction by a single genotype then persisted in the facility. SARS-CoV-2 sequences were investigated at both the consensus and intra-host variation levels. Understanding the variability in SARS-CoV-2 during transmission chains will assist in understanding the spread of this disease and can ultimately inform best practices for mitigation strategies.
Collapse
Affiliation(s)
- Nicole R Sexton
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
- Nebraska Center for Virology, School of Biological Sciences, University of Nebraska-Lincoln, Lincoln, NE, 68504, USA
| | - Parker J Cline
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
| | - Emily N Gallichotte
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
| | - Emily Fitzmeyer
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
| | - Michael C Young
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
| | - Ashley J Janich
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
| | - Kristy L Pabilonia
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
| | - Nicole Ehrhart
- Columbine Health Systems Center for Healthy Aging and Department of Clinical Sciences, Colorado State University, Fort Collins, CO, 80523, USA
| | - Gregory D Ebel
- Department of Microbiology, Immunology and Pathology, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, 80523, USA.
| |
Collapse
|
9
|
Terbot JW, Cooper BS, Good JM, Jensen JD. A simulation framework for modeling the within-patient evolutionary dynamics of SARS-CoV-2. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.13.548462. [PMID: 37503016 PMCID: PMC10370031 DOI: 10.1101/2023.07.13.548462] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
The global impact of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for positive selection are best performed via comparison of empirical data to simulated data wherein evolutionary factors, including mutation and recombination rates, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. While there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intra-host evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them to existing empirical data. Of these, 592 models (~5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intra-host SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed towards strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.
Collapse
Affiliation(s)
- John W Terbot
- Arizona State University, School of Life Sciences, Center for Evolution & Medicine, Tempe, Arizona, United States of America
- University of Montana, Division of Biological Sciences, Missoula, Montana, United States of America
| | - Brandon S. Cooper
- University of Montana, Division of Biological Sciences, Missoula, Montana, United States of America
| | - Jeffrey M. Good
- University of Montana, Division of Biological Sciences, Missoula, Montana, United States of America
| | - Jeffrey D. Jensen
- Arizona State University, School of Life Sciences, Center for Evolution & Medicine, Tempe, Arizona, United States of America
| |
Collapse
|
10
|
Terbot JW, Johri P, Liphardt SW, Soni V, Pfeifer SP, Cooper BS, Good JM, Jensen JD. Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples. PLoS Pathog 2023; 19:e1011265. [PMID: 37018331 PMCID: PMC10075409 DOI: 10.1371/journal.ppat.1011265] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/06/2023] Open
Abstract
Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.
Collapse
Affiliation(s)
- John W Terbot
- University of Montana, Division of Biological Sciences, Missoula, Montana, United States of America
- Arizona State University, School of Life Sciences, Center for Evolution & Medicine, Tempe, Arizona, United States of America
| | - Parul Johri
- Arizona State University, School of Life Sciences, Center for Evolution & Medicine, Tempe, Arizona, United States of America
| | - Schuyler W Liphardt
- University of Montana, Division of Biological Sciences, Missoula, Montana, United States of America
| | - Vivak Soni
- Arizona State University, School of Life Sciences, Center for Evolution & Medicine, Tempe, Arizona, United States of America
| | - Susanne P Pfeifer
- Arizona State University, School of Life Sciences, Center for Evolution & Medicine, Tempe, Arizona, United States of America
| | - Brandon S Cooper
- University of Montana, Division of Biological Sciences, Missoula, Montana, United States of America
| | - Jeffrey M Good
- University of Montana, Division of Biological Sciences, Missoula, Montana, United States of America
| | - Jeffrey D Jensen
- Arizona State University, School of Life Sciences, Center for Evolution & Medicine, Tempe, Arizona, United States of America
| |
Collapse
|
11
|
Árnason E, Koskela J, Halldórsdóttir K, Eldon B. Sweepstakes reproductive success via pervasive and recurrent selective sweeps. eLife 2023; 12:80781. [PMID: 36806325 PMCID: PMC9940914 DOI: 10.7554/elife.80781] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 12/28/2022] [Indexed: 02/22/2023] Open
Abstract
Highly fecund natural populations characterized by high early mortality abound, yet our knowledge about their recruitment dynamics is somewhat rudimentary. This knowledge gap has implications for our understanding of genetic variation, population connectivity, local adaptation, and the resilience of highly fecund populations. The concept of sweepstakes reproductive success, which posits a considerable variance and skew in individual reproductive output, is key to understanding the distribution of individual reproductive success. However, it still needs to be determined whether highly fecund organisms reproduce through sweepstakes and, if they do, the relative roles of neutral and selective sweepstakes. Here, we use coalescent-based statistical analysis of population genomic data to show that selective sweepstakes likely explain recruitment dynamics in the highly fecund Atlantic cod. We show that the Kingman coalescent (modelling no sweepstakes) and the Xi-Beta coalescent (modelling random sweepstakes), including complex demography and background selection, do not provide an adequate fit for the data. The Durrett-Schweinsberg coalescent, in which selective sweepstakes result from recurrent and pervasive selective sweeps of new mutations, offers greater explanatory power. Our results show that models of sweepstakes reproduction and multiple-merger coalescents are relevant and necessary for understanding genetic diversity in highly fecund natural populations. These findings have fundamental implications for understanding the recruitment variation of fish stocks and general evolutionary genomics of high-fecundity organisms.
Collapse
Affiliation(s)
- Einar Árnason
- Institute of Life- and environmental Sciences, University of IcelandReykjavikIceland,Department of Organismal and Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Jere Koskela
- Department of Statistics, University of WarwickCoventryUnited Kingdom
| | - Katrín Halldórsdóttir
- Institute of Life- and environmental Sciences, University of IcelandReykjavikIceland
| | - Bjarki Eldon
- Leibniz Institute for Evolution and Biodiversity Science, Museum für NaturkundeBerlinGermany
| |
Collapse
|
12
|
Abstract
We discuss the genetic, demographic, and selective forces that are likely to be at play in restricting observed levels of DNA sequence variation in natural populations to a much smaller range of values than would be expected from the distribution of census population sizes alone-Lewontin's Paradox. While several processes that have previously been strongly emphasized must be involved, including the effects of direct selection and genetic hitchhiking, it seems unlikely that they are sufficient to explain this observation without contributions from other factors. We highlight a potentially important role for the less-appreciated contribution of population size change; specifically, the likelihood that many species and populations may be quite far from reaching the relatively high equilibrium diversity values that would be expected given their current census sizes.
Collapse
Affiliation(s)
- Brian Charlesworth
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Jeffrey D Jensen
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| |
Collapse
|
13
|
Sabin S, Morales-Arce AY, Pfeifer SP, Jensen JD. The impact of frequently neglected model violations on bacterial recombination rate estimation: a case study in Mycobacterium canettii and Mycobacterium tuberculosis. G3 (BETHESDA, MD.) 2022; 12:jkac055. [PMID: 35253851 PMCID: PMC9073693 DOI: 10.1093/g3journal/jkac055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Accepted: 02/28/2022] [Indexed: 12/04/2022]
Abstract
Mycobacterium canettii is a causative agent of tuberculosis in humans, along with the members of the Mycobacterium tuberculosis complex. Frequently used as an outgroup to the M. tuberculosis complex in phylogenetic analyses, M. canettii is thought to offer the best proxy for the progenitor species that gave rise to the complex. Here, we leverage whole-genome sequencing data and biologically relevant population genomic models to compare the evolutionary dynamics driving variation in the recombining M. canettii with that in the nonrecombining M. tuberculosis complex, and discuss differences in observed genomic diversity in the light of expected levels of Hill-Robertson interference. In doing so, we highlight the methodological challenges of estimating recombination rates through traditional population genetic approaches using sequences called from populations of microorganisms and evaluate the likely mis-inference that arises owing to a neglect of common model violations including purifying selection, background selection, progeny skew, and population size change. In addition, we compare performance when full within-host polymorphism data are utilized, versus the more common approach of basing analyses on within-host consensus sequences.
Collapse
Affiliation(s)
- Susanna Sabin
- Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
| | - Ana Y Morales-Arce
- Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
| | - Susanne P Pfeifer
- Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
| | - Jeffrey D Jensen
- Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
| |
Collapse
|
14
|
Morales-Arce AY, Johri P, Jensen JD. Inferring the distribution of fitness effects in patient-sampled and experimental virus populations: two case studies. Heredity (Edinb) 2022; 128:79-87. [PMID: 34987185 PMCID: PMC8728706 DOI: 10.1038/s41437-021-00493-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2021] [Revised: 12/12/2021] [Accepted: 12/13/2021] [Indexed: 11/19/2022] Open
Abstract
We here propose an analysis pipeline for inferring the distribution of fitness effects (DFE) from either patient-sampled or experimentally-evolved viral populations, that explicitly accounts for non-Wright-Fisher and non-equilibrium population dynamics inherent to pathogens. We examine the performance of this approach via extensive power and performance analyses, and highlight two illustrative applications - one from an experimentally-passaged RNA virus, and the other from a clinically-sampled DNA virus. Finally, we discuss how such DFE inference may shed light on major research questions in virus evolution, ranging from a quantification of the population genetic processes governing genome size, to the role of Hill-Robertson interference in dictating adaptive outcomes, to the potential design of novel therapeutic approaches to eradicate within-patient viral populations via induced mutational meltdown.
Collapse
Affiliation(s)
- Ana Y. Morales-Arce
- grid.215654.10000 0001 2151 2636Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ USA
| | - Parul Johri
- grid.215654.10000 0001 2151 2636Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ USA
| | - Jeffrey D. Jensen
- grid.215654.10000 0001 2151 2636Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ USA
| |
Collapse
|
15
|
Retel C, Kowallik V, Becks L, Feulner PGD. Strong Selection and High Mutation Supply Characterize Experimental Chlorovirus Evolution. Virus Evol 2022; 8:veac003. [PMID: 35169490 PMCID: PMC8838748 DOI: 10.1093/ve/veac003] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Revised: 01/14/2022] [Accepted: 01/25/2022] [Indexed: 11/24/2022] Open
Abstract
Characterizing how viruses evolve expands our understanding of the underlying fundamental processes, such as mutation, selection and drift. One group of viruses whose evolution has not yet been extensively studied is the Phycodnaviridae, a globally abundant family of aquatic large double-stranded (ds) DNA viruses. Here we studied the evolutionary change of Paramecium bursaria chlorella virus 1 during experimental coevolution with its algal host. We used pooled genome sequencing of six independently evolved populations to characterize genomic change over five time points. Across six experimental replicates involving either strong or weak demographic fluctuations, we found single nucleotide polymorphisms (SNPs) at sixty-seven sites. The occurrence of genetic variants was highly repeatable, with just two of the SNPs found in only a single experimental replicate. Three genes A122/123R, A140/145R and A540L showed an excess of variable sites, providing new information about potential targets of selection during Chlorella–Chlorovirus coevolution. Our data indicated that the studied populations were not mutation-limited and experienced strong positive selection. Our investigation highlighted relevant processes governing the evolution of aquatic large dsDNA viruses, which ultimately contributes to a better understanding of the functioning of natural aquatic ecosystems.
Collapse
Affiliation(s)
- Cas Retel
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Bio-geochemistry, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Seestrasse 79, Kastanienbaum 6047, Switzerland
- Division of Aquatic Ecology, Institute of Ecology and Evolution, University of Bern, Baltzerstrasse 6, Bern 3012, Switzerland
| | | | | | | |
Collapse
|
16
|
Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context. EPMA J 2022; 13:149-175. [PMID: 35013687 PMCID: PMC8732965 DOI: 10.1007/s13167-021-00267-w] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Accepted: 12/04/2021] [Indexed: 12/12/2022]
Abstract
Aims The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM. Methods In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein–protein interface. Results Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike “boat-shaped” receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the “stern” or “bow” regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the “hull” region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities. Conclusion In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease. Supplementary Information The online version contains supplementary material available at 10.1007/s13167-021-00267-w.
Collapse
|
17
|
Baumdicker F, Bisschop G, Goldstein D, Gower G, Ragsdale AP, Tsambos G, Zhu S, Eldon B, Ellerman EC, Galloway JG, Gladstein AL, Gorjanc G, Guo B, Jeffery B, Kretzschmar WW, Lohse K, Matschiner M, Nelson D, Pope NS, Quinto-Cortés CD, Rodrigues MF, Saunack K, Sellinger T, Thornton K, van Kemenade H, Wohns AW, Wong Y, Gravel S, Kern AD, Koskela J, Ralph PL, Kelleher J. Efficient ancestry and mutation simulation with msprime 1.0. Genetics 2021; 220:6460344. [PMID: 34897427 PMCID: PMC9176297 DOI: 10.1093/genetics/iyab229] [Citation(s) in RCA: 91] [Impact Index Per Article: 30.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Accepted: 12/03/2021] [Indexed: 11/13/2022] Open
Abstract
Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this, a large number of specialized simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and the tskit library. We summarize msprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialized alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.
Collapse
Affiliation(s)
- Franz Baumdicker
- Cluster of Excellence "Controlling Microbes to Fight Infections", Mathematical and Computational Population Genetics, University of Tübingen, 72076 Tübingen, Germany
| | - Gertjan Bisschop
- Institute of Evolutionary Biology,The University of Edinburgh, EH9 3FL, UK
| | - Daniel Goldstein
- Khoury College of Computer Sciences, Northeastern University, MA 02115, USA.,No affiliation
| | - Graham Gower
- Lundbeck GeoGenetics Centre, Globe Institute, University of Copenhagen, 1350 Copenhagen K, Denmark
| | - Aaron P Ragsdale
- Department of Integrative Biology, University of Wisconsin-Madison, WI 53706, USA
| | - Georgia Tsambos
- Melbourne Integrative Genomics, School of Mathematics and Statistics, University of Melbourne, Victoria, 3010, Australia
| | - Sha Zhu
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, OX3 7LF, UK
| | - Bjarki Eldon
- Leibniz Institute for Evolution and Biodiversity Science,Museum für Naturkunde Berlin, 10115, Germany
| | | | - Jared G Galloway
- Institute of Ecology and Evolution, Department of Biology, University of Oregon, OR 97403-5289, USA.,Computational Biology Program, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
| | - Ariella L Gladstein
- Department of Genetics, University of North Carolina at Chapel Hill, NC 27599-7264, USA.,Embark Veterinary, Inc., Boston, MA 02111, USA
| | - Gregor Gorjanc
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, EH25 9RG, UK
| | - Bing Guo
- Institute for Genome Sciences,University of Maryland School of Medicine, Baltimore, MD, 21201, USA
| | - Ben Jeffery
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, OX3 7LF, UK
| | - Warren W Kretzschmar
- Center for Hematology and Regenerative Medicine, Karolinska Institute, 141 83 Huddinge, Sweden
| | - Konrad Lohse
- Institute of Evolutionary Biology,The University of Edinburgh, EH9 3FL, UK
| | | | - Dominic Nelson
- Department of Human Genetics, McGill University, Montréal, QC H3A 0C7, Canada
| | - Nathaniel S Pope
- Department of Entomology, Pennsylvania State University, PA 16802, USA
| | - Consuelo D Quinto-Cortés
- National Laboratory of Genomics for Biodiversity (LANGEBIO), Unit of Advanced Genomics, CINVESTAV, Irapuato, Mexico
| | - Murillo F Rodrigues
- Institute of Ecology and Evolution, Department of Biology, University of Oregon, OR 97403-5289, USA
| | - Kumar Saunack
- IIT Bombay, Powai, Mumbai 400 076, Maharashtra, India
| | - Thibaut Sellinger
- Professorship for Population Genetics, Department of Life Science Systems, Technical University of Munich, 85354 Freising, Germany
| | - Kevin Thornton
- Ecology and Evolutionary Biology, University of California, Irvine, CA 92697, USA
| | | | - Anthony W Wohns
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, OX3 7LF, UK.,Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Yan Wong
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, OX3 7LF, UK
| | - Simon Gravel
- Department of Human Genetics, McGill University, Montréal, QC H3A 0C7, Canada
| | - Andrew D Kern
- Institute of Ecology and Evolution, Department of Biology, University of Oregon, OR 97403-5289, USA
| | - Jere Koskela
- Department of Statistics, University of Warwick, CV4 7AL, UK
| | - Peter L Ralph
- Institute of Ecology and Evolution, Department of Biology, University of Oregon, OR 97403-5289, USA.,Department of Mathematics, University of Oregon, OR 97403-5289 USA
| | - Jerome Kelleher
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, OX3 7LF, UK
| |
Collapse
|
18
|
Charlesworth B, Jensen JD. Effects of Selection at Linked Sites on Patterns of Genetic Variability. ANNUAL REVIEW OF ECOLOGY, EVOLUTION, AND SYSTEMATICS 2021; 52:177-197. [PMID: 37089401 PMCID: PMC10120885 DOI: 10.1146/annurev-ecolsys-010621-044528] [Citation(s) in RCA: 46] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Patterns of variation and evolution at a given site in a genome can be strongly influenced by the effects of selection at genetically linked sites. In particular, the recombination rates of genomic regions correlate with their amount of within-population genetic variability, the degree to which the frequency distributions of DNA sequence variants differ from their neutral expectations, and the levels of adaptation of their functional components. We review the major population genetic processes that are thought to lead to these patterns, focusing on their effects on patterns of variability: selective sweeps, background selection, associative overdominance, and Hill–Robertson interference among deleterious mutations. We emphasize the difficulties in distinguishing among the footprints of these processes and disentangling them from the effects of purely demographic factors such as population size changes. We also discuss how interactions between selective and demographic processes can significantly affect patterns of variability within genomes.
Collapse
Affiliation(s)
- Brian Charlesworth
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3FL, United Kingdom
| | - Jeffrey D. Jensen
- School of Life Sciences, Arizona State University, Tempe, Arizona 85281, USA
| |
Collapse
|
19
|
Menardo F, Gagneux S, Freund F. Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis. Mol Biol Evol 2021; 38:290-306. [PMID: 32667991 PMCID: PMC8480183 DOI: 10.1093/molbev/msaa179] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders). Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.
Collapse
Affiliation(s)
- Fabrizio Menardo
- Department of Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Basel, Switzerland.,University of Basel, Basel, Switzerland
| | - Sébastien Gagneux
- Department of Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Basel, Switzerland.,University of Basel, Basel, Switzerland
| | - Fabian Freund
- Department of Plant Biodiversity and Breeding Informatics, Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, Stuttgart, Germany
| |
Collapse
|
20
|
Freund F, Siri-Jégousse A. The impact of genetic diversity statistics on model selection between coalescents. Comput Stat Data Anal 2021. [DOI: 10.1016/j.csda.2020.107055] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
21
|
Abstract
AbstractWe construct a multitype constant-size population model allowing for general selective interactions as well as extreme reproductive events. Our multidimensional model aims for the generality of adaptive dynamics and the tractability of population genetics. It generalises the idea of Krone and Neuhauser [39] and González Casanova and Spanò [29], who represented the selection by allowing individuals to sample several potential parents in the previous generation before choosing the ‘strongest’ one, by allowing individuals to use any rule to choose their parent. The type of the newborn can even not be one of the types of the potential parents, which allows modelling mutations. Via a large population limit, we obtain a generalisation of
$\Lambda$
-Fleming–Viot processes, with a diffusion term and a general frequency-dependent selection, which allows for non-transitive interactions between the different types present in the population. We provide some properties of these processes related to extinction and fixation events, and give conditions for them to be realised as unique strong solutions of multidimensional stochastic differential equations with jumps. Finally, we illustrate the generality of our model with applications to some classical biological interactions. This framework provides a natural bridge between two of the most prominent modelling frameworks of biological evolution: population genetics and eco-evolutionary models.
Collapse
|
22
|
Jensen JD, Stikeleather RA, Kowalik TF, Lynch M. Imposed mutational meltdown as an antiviral strategy. Evolution 2020; 74:2549-2559. [PMID: 33047822 PMCID: PMC7993354 DOI: 10.1111/evo.14107] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2020] [Revised: 09/30/2020] [Accepted: 10/10/2020] [Indexed: 12/25/2022]
Abstract
Following widespread infections of the most recent coronavirus known to infect humans, SARS‐CoV‐2, attention has turned to potential therapeutic options. With no drug or vaccine yet approved, one focal point of research is to evaluate the potential value of repurposing existing antiviral treatments, with the logical strategy being to identify at least a short‐term intervention to prevent within‐patient progression, while long‐term vaccine strategies unfold. Here, we offer an evolutionary/population‐genetic perspective on one approach that may overwhelm the capacity for pathogen defense (i.e., adaptation) – induced mutational meltdown – providing an overview of key concepts, review of previous theoretical and experimental work of relevance, and guidance for future research. Applied with appropriate care, including target specificity, induced mutational meltdown may provide a general, rapidly implemented approach for the within‐patient eradication of a wide range of pathogens or other undesirable microorganisms.
Collapse
Affiliation(s)
- Jeffrey D Jensen
- School of Life Sciences, Arizona State University, Tempe, Arizona, 85281.,Center for Evolution & Medicine, Arizona State University, Tempe, Arizona, 85281
| | - Ryan A Stikeleather
- Biodesign Center for Mechanisms of Evolution, Arizona State University, Tempe, Arizona, 85281
| | - Timothy F Kowalik
- Department of Microbiology and Physiological Systems, University of Massachusetts Medical School, Worcester, Massachusetts, 01655
| | - Michael Lynch
- School of Life Sciences, Arizona State University, Tempe, Arizona, 85281.,Biodesign Center for Mechanisms of Evolution, Arizona State University, Tempe, Arizona, 85281
| |
Collapse
|
23
|
Morales-Arce AY, Sabin SJ, Stone AC, Jensen JD. The population genomics of within-host Mycobacterium tuberculosis. Heredity (Edinb) 2020; 126:1-9. [PMID: 33060846 DOI: 10.1038/s41437-020-00377-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Revised: 10/02/2020] [Accepted: 10/03/2020] [Indexed: 11/09/2022] Open
Abstract
Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.
Collapse
Affiliation(s)
- Ana Y Morales-Arce
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA.
| | - Susanna J Sabin
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA
| | - Anne C Stone
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA.,School of Human Evolution and Social Change, Arizona State University, Tempe, AZ, USA
| | - Jeffrey D Jensen
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA. .,School of Life Sciences, Arizona State University, Tempe, AZ, USA.
| |
Collapse
|
24
|
Abstract
Natural highly fecund populations abound. These range from viruses to gadids. Many highly fecund populations are economically important. Highly fecund populations provide an important contrast to the low-fecundity organisms that have traditionally been applied in evolutionary studies. A key question regarding high fecundity is whether large numbers of offspring are produced on a regular basis, by few individuals each time, in a sweepstakes mode of reproduction. Such reproduction characteristics are not incorporated into the classical Wright-Fisher model, the standard reference model of population genetics, or similar types of models, in which each individual can produce only small numbers of offspring relative to the population size. The expected genomic footprints of population genetic models of sweepstakes reproduction are very different from those of the Wright-Fisher model. A key, immediate issue involves identifying the footprints of sweepstakes reproduction in genomic data. Whole-genome sequencing data can be used to distinguish the patterns made by sweepstakes reproduction from the patterns made by population growth in a population evolving according to the Wright-Fisher model (or similar models). If the hypothesis of sweepstakes reproduction cannot be rejected, then models of sweepstakes reproduction and associated multiple-merger coalescents will become at least as relevant as the Wright-Fisher model (or similar models) and the Kingman coalescent, the cornerstones of mathematical population genetics, in further discussions of evolutionary genomics of highly fecund populations.
Collapse
Affiliation(s)
- Bjarki Eldon
- Leibniz Institute for Evolution and Biodiversity Science, Museum für Naturkunde, D-10115 Berlin, Germany;
| |
Collapse
|
25
|
Harris RB, Jensen JD. Considering Genomic Scans for Selection as Coalescent Model Choice. Genome Biol Evol 2020; 12:871-877. [PMID: 32396636 PMCID: PMC7313662 DOI: 10.1093/gbe/evaa093] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/06/2020] [Indexed: 12/17/2022] Open
Abstract
First inspired by the seminal work of Lewontin and Krakauer (1973. Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics 74(1):175-195.) and Maynard Smith and Haigh (1974. The hitch-hiking effect of a favourable gene. Genet Res. 23(1):23-35.), genomic scans for positive selection remain a widely utilized tool in modern population genomic analysis. Yet, the relative frequency and genomic impact of selective sweeps have remained a contentious point in the field for decades, largely owing to an inability to accurately identify their presence and quantify their effects-with current methodologies generally being characterized by low true-positive rates and/or high false-positive rates under many realistic demographic models. Most of these approaches are based on Wright-Fisher assumptions and the Kingman coalescent and generally rely on detecting outlier regions which do not conform to these neutral expectations. However, previous theoretical results have demonstrated that selective sweeps are well characterized by an alternative class of model known as the multiple-merger coalescent. Taken together, this suggests the possibility of not simply identifying regions which reject the Kingman, but rather explicitly testing the relative fit of a genomic window to the multiple-merger coalescent. We describe the advantages of such an approach, which owe to the branching structure differentiating selective and neutral models, and demonstrate improved power under certain demographic scenarios relative to a commonly used approach. However, regions of the demographic parameter space continue to exist in which neither this approach nor existing methodologies have sufficient power to detect selective sweeps.
Collapse
|
26
|
Morales-Arce AY, Harris RB, Stone AC, Jensen JD. Evaluating the contributions of purifying selection and progeny-skew in dictating within-host Mycobacterium tuberculosis evolution. Evolution 2020; 74:992-1001. [PMID: 32233086 DOI: 10.1111/evo.13954] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2019] [Accepted: 03/08/2020] [Indexed: 12/28/2022]
Abstract
The within-host evolutionary dynamics of tuberculosis (TB) remain unclear, and underlying biological characteristics render standard population genetic approaches based upon the Wright-Fisher model largely inappropriate. In addition, the compact genome combined with an absence of recombination is expected to result in strong purifying selection effects. Thus, it is imperative to establish a biologically relevant evolutionary framework incorporating these factors in order to enable an accurate study of this important human pathogen. Further, such a model is critical for inferring fundamental evolutionary parameters related to patient treatment, including mutation rates and the severity of infection bottlenecks. We here implement such a model and infer the underlying evolutionary parameters governing within-patient evolutionary dynamics. Results demonstrate that the progeny skew associated with the clonal nature of TB severely reduces genetic diversity and that the neglect of this parameter in previous studies has led to significant mis-inference of mutation rates. As such, our results suggest an underlying de novo mutation rate that is considerably faster than previously inferred, and a progeny distribution differing significantly from Wright-Fisher assumptions. This inference represents a more appropriate evolutionary null model, against which the periodic effects of positive selection, associated with drug-resistance for example, may be better assessed.
Collapse
Affiliation(s)
- Ana Y Morales-Arce
- Center for Evolution and Medicine, Arizona State University, Tempe, Arizona, USA
| | - Rebecca B Harris
- Center for Evolution and Medicine, Arizona State University, Tempe, Arizona, USA
| | - Anne C Stone
- Center for Evolution and Medicine, Arizona State University, Tempe, Arizona, USA.,School of Human Evolution and Social Change, Arizona State University, Tempe, Arizona, USA
| | - Jeffrey D Jensen
- Center for Evolution and Medicine, Arizona State University, Tempe, Arizona, USA.,School of Life Sciences, Arizona State University, Tempe, Arizona, USA
| |
Collapse
|
27
|
Cannings models, population size changes and multiple-merger coalescents. J Math Biol 2020; 80:1497-1521. [PMID: 32008102 PMCID: PMC7052052 DOI: 10.1007/s00285-020-01470-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2019] [Revised: 01/10/2020] [Indexed: 11/17/2022]
Abstract
Multiple-merger coalescents, e.g. \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\varLambda $$\end{document}Λ-n-coalescents, have been proposed as models of the genealogy of n sampled individuals for a range of populations whose genealogical structures are not captured well by Kingman’s n-coalescent. \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\varLambda $$\end{document}Λ-n-coalescents can be seen as the limit process of the discrete genealogies of Cannings models with fixed population size, when time is rescaled and population size \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$N\rightarrow \infty $$\end{document}N→∞. As established for Kingman’s n-coalescent, moderate population size fluctuations in the discrete population model should be reflected by a time-change of the limit coalescent. For \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\varLambda $$\end{document}Λ-n-coalescents, this has been explicitly shown for only a limited subclass of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\varLambda $$\end{document}Λ-n-coalescents and exponentially growing populations. This article gives a more general construction of time-changed \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\varLambda $$\end{document}Λ-n-coalescents as limits of specific Cannings models with rather arbitrary time changes.
Collapse
|
28
|
Guermouche H, Burrel S, Mercier-Darty M, Kofman T, Rogier O, Pawlotsky JM, Boutolleau D, Rodriguez C. Characterization of the dynamics of human cytomegalovirus resistance to antiviral drugs by ultra-deep sequencing. Antiviral Res 2019; 173:104647. [PMID: 31706899 DOI: 10.1016/j.antiviral.2019.104647] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2019] [Revised: 09/30/2019] [Accepted: 11/04/2019] [Indexed: 12/12/2022]
Abstract
Prophylactic or preemptive treatment strategies are required to prevent human cytomegalovirus (CMV) infections in transplant recipients. However, treatment failure occurs when CMV resistant-associated variants (RAVs) are selected. Although the diversity of CMV is lower than that of RNA viruses, CMV appears to show some genetic instability, with possible minor emerging resistance that may be undetectable by Sanger sequencing. We aimed to examine CMV-resistance mutations over time by ultra-deep sequencing (UDS) and Sanger sequencing in a kidney transplant recipient experiencing CMV infection. This patient showed a transient response to three different antiviral drugs (valganciclovir, foscarnet, and maribavir) and four episodes of CMV resistance over two years. The full-length UL97 (2.3kpb) and partial UL54 (2.4kpb) CMV genes were studied by UDS and Sanger sequencing and linkage mutations calculated to determine RAVs. We detected four major and five minor resistance mutations. Minor resistant variants (2-20%) were detected by UDS, whereas major resistance substitutions (>20%) were identified by both UDS and Sanger method. We detected cross-resistance to three drugs, despite high CMV loads, suggesting that the fitness of the viral mutants was not impaired. In conclusion, CMV showed complex dynamic of resistance under antiviral drug pressure, as described for highly variable viruses. The emergence of successive RAVs constitutes a clinically challenging complication and contributes to the difficulty of therapeutic management of patients.
Collapse
Affiliation(s)
- Hélène Guermouche
- Laboratoire de Virologie, CHU Henri Mondor (AP-HP), INSERM U955 Eq18, Plateforme « Génomiques », IMRB, UPEC, Créteil, France
| | - Sonia Burrel
- Centre National de Référence Herpèsvirus (laboratoire associé), Laboratoire de Virologie, Hôpital Universitaire La Pitié-Salpêtrière, GHU AP-PH. Sorbonne Université (AP-HP), INSERM U1136, iPLESP, Sorbonne Université, Paris, France
| | - Mélanie Mercier-Darty
- Laboratoire de Virologie, CHU Henri Mondor (AP-HP), INSERM U955 Eq18, Plateforme « Génomiques », IMRB, UPEC, Créteil, France
| | - Thomas Kofman
- Service de Néphrologie, Hôpital Universitaire Henri Mondor (AP-HP), Créteil, France
| | - Olivier Rogier
- Laboratoire de Virologie, CHU Henri Mondor (AP-HP), INSERM U955 Eq18, Plateforme « Génomiques », IMRB, UPEC, Créteil, France
| | - Jean-Michel Pawlotsky
- Laboratoire de Virologie, CHU Henri Mondor (AP-HP), INSERM U955 Eq18, Plateforme « Génomiques », IMRB, UPEC, Créteil, France
| | - David Boutolleau
- Centre National de Référence Herpèsvirus (laboratoire associé), Laboratoire de Virologie, Hôpital Universitaire La Pitié-Salpêtrière, GHU AP-PH. Sorbonne Université (AP-HP), INSERM U1136, iPLESP, Sorbonne Université, Paris, France
| | - Christophe Rodriguez
- Laboratoire de Virologie, CHU Henri Mondor (AP-HP), INSERM U955 Eq18, Plateforme « Génomiques », IMRB, UPEC, Créteil, France.
| |
Collapse
|
29
|
Retel C, Kowallik V, Huang W, Werner B, Künzel S, Becks L, Feulner PGD. The feedback between selection and demography shapes genomic diversity during coevolution. SCIENCE ADVANCES 2019; 5:eaax0530. [PMID: 31616788 PMCID: PMC6774728 DOI: 10.1126/sciadv.aax0530] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/18/2019] [Accepted: 09/06/2019] [Indexed: 06/10/2023]
Abstract
Species interactions and coevolution are integral to ecological communities, but we lack empirical information on when and how these interactions generate and purge genetic diversity. Using genomic time series data from host-virus experiments, we found that coevolution occurs through consecutive selective sweeps in both species, with temporal consistency across replicates. Sweeps were accompanied by phenotypic change (resistance or infectivity increases) and expansions in population size. In the host, population expansion enabled rapid generation of genetic diversity in accordance with neutral processes. Viral molecular evolution was, in contrast, confined to few genes, all putative targets of selection. This study demonstrates that molecular evolution during species interactions is shaped by both eco-evolutionary feedback dynamics and interspecific differences in how genetic diversity is generated and maintained.
Collapse
Affiliation(s)
- Cas Retel
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
- Division of Aquatic Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Vienna Kowallik
- Community Dynamics Group, Department of Evolutionary Ecology, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Weini Huang
- Group of Theoretical Biology, The State Key Laboratory of Biocontrol, School of Life Science, Sun Yat-sen University, Guangzhou, China
- Complex Systems and Networks Research Group, School of Mathematical Sciences, Queen Mary University of London, London, UK
| | - Benjamin Werner
- Evolutionary Genomics and Modelling Lab, Centre for Evolution and Cancer, The Institute of Cancer Research, London, UK
| | - Sven Künzel
- Department Evolutionary Genetics, Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Lutz Becks
- Community Dynamics Group, Department of Evolutionary Ecology, Max Planck Institute for Evolutionary Biology, Plön, Germany
- Kiel Evolution Center, Biologiezentrum, Kiel, Germany
| | - Philine G. D. Feulner
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
- Division of Aquatic Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| |
Collapse
|
30
|
Filatov DA. Extreme Lewontin's Paradox in Ubiquitous Marine Phytoplankton Species. Mol Biol Evol 2019; 36:4-14. [PMID: 30351418 DOI: 10.1093/molbev/msy195] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Larger populations are expected to have larger genetic diversity. However, as pointed out by Lewontin in 1974, the range of population sizes exceeds the range of genetic diversity by many orders of magnitude (a.k.a. "Lewontin's paradox," LP). The reasons for LP remain obscure. Here, This paper reports an extreme case of LP in astronomically large populations of the ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta)-the species that accounts for 10-20% of primary productivity in the oceans and its blooms are so extensive that they are visible from space. This study demonstrates that despite the wide distribution and enormous population size, the world-wide sample of E. huxleyi strains with sequenced genomes represents a single cohesive species and contains surprisingly limited genetic diversity (π ∼ 0.006 per silent site). The patterns of polymorphism reveal even larger populations in the past, and frequent recombination (ρ ∼ 0.006) throughout the genome, ruling out demographic history and asexual reproduction as possible causes of low polymorphism in E. huxleyi. Natural selection wiping out genetic diversity at linked sites (a.k.a. "genetic draft") must be strong and frequent to account for low polymorphism in E. huxleyi. This study sheds the first light on poorly understood evolutionary genetic processes in astronomically large populations of marine microplankton.
Collapse
Affiliation(s)
- Dmitry A Filatov
- Department of Plant Sciences, University of Oxford, South Parks Road, Oxford OX1 3RB, United Kingdom
| |
Collapse
|
31
|
Ecological and Evolutionary Processes Shaping Viral Genetic Diversity. Viruses 2019; 11:v11030220. [PMID: 30841497 PMCID: PMC6466605 DOI: 10.3390/v11030220] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Revised: 02/22/2019] [Accepted: 02/27/2019] [Indexed: 02/07/2023] Open
Abstract
The contemporary genomic diversity of viruses is a result of the continuous and dynamic interaction of past ecological and evolutionary processes. Thus, genome sequences of viruses can be a valuable source of information about these processes. In this review, we first describe the relevant processes shaping viral genomic variation, with a focus on the role of host–virus coevolution and its potential to give rise to eco-evolutionary feedback loops. We further give a brief overview of available methodology designed to extract information about these processes from genomic data. Short generation times and small genomes make viruses ideal model systems to study the joint effect of complex coevolutionary and eco-evolutionary interactions on genetic evolution. This complexity, together with the diverse array of lifetime and reproductive strategies in viruses ask for extensions of existing inference methods, for example by integrating multiple information sources. Such integration can broaden the applicability of genetic inference methods and thus further improve our understanding of the role viruses play in biological communities.
Collapse
|
32
|
Inferring Demography and Selection in Organisms Characterized by Skewed Offspring Distributions. Genetics 2019; 211:1019-1028. [PMID: 30651284 DOI: 10.1534/genetics.118.301684] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2018] [Accepted: 01/15/2019] [Indexed: 01/01/2023] Open
Abstract
The recent increase in time-series population genomic data from experimental, natural, and ancient populations has been accompanied by a promising growth in methodologies for inferring demographic and selective parameters from such data. However, these methods have largely presumed that the populations of interest are well-described by the Kingman coalescent. In reality, many groups of organisms, including viruses, marine organisms, and some plants, protists, and fungi, typified by high variance in progeny number, may be best characterized by multiple-merger coalescent models. Estimation of population genetic parameters under Wright-Fisher assumptions for these organisms may thus be prone to serious mis-inference. We propose a novel method for the joint inference of demography and selection under the Ψ-coalescent model, termed Multiple-Merger Coalescent Approximate Bayesian Computation, or MMC-ABC. We first demonstrate mis-inference under the Kingman, and then exhibit the superior performance of MMC-ABC under conditions of skewed offspring distributions. In order to highlight the utility of this approach, we reanalyzed previously published drug-selection lines of influenza A virus. We jointly inferred the extent of progeny-skew inherent to viral replication and identified putative drug-resistance mutations.
Collapse
|
33
|
Ormond L, Liu P, Matuszewski S, Renzette N, Bank C, Zeldovich K, Bolon DN, Kowalik TF, Finberg RW, Jensen JD, Wang JP. The Combined Effect of Oseltamivir and Favipiravir on Influenza A Virus Evolution. Genome Biol Evol 2017; 9:1913-1924. [PMID: 28854600 PMCID: PMC5570085 DOI: 10.1093/gbe/evx138] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/18/2017] [Indexed: 01/14/2023] Open
Abstract
Influenza virus inflicts a heavy death toll annually and resistance to existing antiviral drugs has generated interest in the development of agents with novel mechanisms of action. Favipiravir is an antiviral drug that acts by increasing the genome-wide mutation rate of influenza A virus (IAV). Potential synergistic benefits of combining oseltamivir and favipiravir have been demonstrated in animal models of influenza, but the population-level effects of combining the drugs are unknown. In order to elucidate the underlying evolutionary processes at play, we performed genome-wide sequencing of IAV experimental populations subjected to serial passaging in vitro under a combined protocol of oseltamivir and favipiravir. We describe the interplay between mutation, selection, and genetic drift that ultimately culminates in population extinction. In particular, selective sweeps around oseltamivir resistance mutations reduce genome-wide variation while deleterious mutations hitchhike to fixation given the increased mutational load generated by favipiravir. This latter effect reduces viral fitness and accelerates extinction compared with IAV populations treated with favipiravir alone, but risks spreading both established and newly emerging mutations, including possible drug resistance mutations, if transmission occurs before the viral populations are eradicated.
Collapse
Affiliation(s)
- Louise Ormond
- École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.,Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland
| | - Ping Liu
- Department of Medicine, University of Massachusetts Medical School
| | - Sebastian Matuszewski
- École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.,Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland
| | - Nicholas Renzette
- Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland.,Department of Microbiology and Physiological Systems, University of Massachusetts Medical School
| | - Claudia Bank
- École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.,Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland.,Instituto Gulbenkian de Ciência, Oeiras, Portugal
| | - Konstantin Zeldovich
- Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School
| | - Daniel N Bolon
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School
| | - Timothy F Kowalik
- Department of Microbiology and Physiological Systems, University of Massachusetts Medical School
| | - Robert W Finberg
- Department of Medicine, University of Massachusetts Medical School
| | - Jeffrey D Jensen
- École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.,Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland.,School of Life Sciences, Center for Evolution & Medicine, Arizona State University
| | - Jennifer P Wang
- Department of Medicine, University of Massachusetts Medical School
| |
Collapse
|
34
|
Coalescent Processes with Skewed Offspring Distributions and Nonequilibrium Demography. Genetics 2017; 208:323-338. [PMID: 29127263 DOI: 10.1534/genetics.117.300499] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2017] [Accepted: 10/30/2017] [Indexed: 11/18/2022] Open
Abstract
Nonequilibrium demography impacts coalescent genealogies leaving detectable, well-studied signatures of variation. However, similar genomic footprints are also expected under models of large reproductive skew, posing a serious problem when trying to make inference. Furthermore, current approaches consider only one of the two processes at a time, neglecting any genomic signal that could arise from their simultaneous effects, preventing the possibility of jointly inferring parameters relating to both offspring distribution and population history. Here, we develop an extended Moran model with exponential population growth, and demonstrate that the underlying ancestral process converges to a time-inhomogeneous psi-coalescent. However, by applying a nonlinear change of time scale-analogous to the Kingman coalescent-we find that the ancestral process can be rescaled to its time-homogeneous analog, allowing the process to be simulated quickly and efficiently. Furthermore, we derive analytical expressions for the expected site-frequency spectrum under the time-inhomogeneous psi-coalescent, and develop an approximate-likelihood framework for the joint estimation of the coalescent and growth parameters. By means of extensive simulation, we demonstrate that both can be estimated accurately from whole-genome data. In addition, not accounting for demography can lead to serious biases in the inferred coalescent model, with broad implications for genomic studies ranging from ecology to conservation biology. Finally, we use our method to analyze sequence data from Japanese sardine populations, and find evidence of high variation in individual reproductive success, but few signs of a recent demographic expansion.
Collapse
|
35
|
Evolution of highly fecund haploid populations. Theor Popul Biol 2017; 119:48-56. [PMID: 29111301 DOI: 10.1016/j.tpb.2017.10.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2017] [Revised: 10/14/2017] [Accepted: 10/17/2017] [Indexed: 11/24/2022]
Abstract
We consider a model of viability selection in a highly fecund haploid population with sweepstakes reproduction. We use simulations to estimate the time until the allelic type with highest fitness has reached high frequency in a finite population. We compare the time between two reproduction modes of high and low fecundity. We also consider the probability that the allelic type with highest fitness is lost from the population before reaching high frequency. Our simulation results indicate that highly fecund populations can evolve faster (in some cases much faster) than populations of low fecundity. However, high fecundity and sweepstakes reproduction also confer much higher risk of losing the allelic type with highest fitness from the population by chance. The impact of selection on driving alleles to high frequency varies depending on the trait value conferring highest fitness; in some cases the effect of selection can hardly be detected.
Collapse
|
36
|
Matuszewski S, Ormond L, Bank C, Jensen JD. Two sides of the same coin: A population genetics perspective on lethal mutagenesis and mutational meltdown. Virus Evol 2017; 3:vex004. [PMID: 29977604 PMCID: PMC6007402 DOI: 10.1093/ve/vex004] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
The extinction of RNA virus populations upon application of a mutagenic drug is frequently referred to as evidence for the existence of an error threshold, above which the population cannot sustain the mutational load. To explain the extinction process after reaching this threshold, models of lethal mutagenesis have been proposed, in which extinction is described as a deterministic (and thus population size-independent) process. As a separate body of literature, the population genetics community has developed models of mutational meltdown, which focus on the stochastic (and thus population-size dependent) processes governing extinction. However, recent extensions of both models have blurred these boundaries. Here, we first clarify definitions in terms of assumptions, expectations, and relevant parameter spaces, and then assess similarities and differences. As concepts from both fields converge, we argue for a unified theoretical framework that is focused on the evolutionary processes at play, rather than dispute over terminology.
Collapse
Affiliation(s)
- Sebastian Matuszewski
- School of Life Sciences, École Polytechnique Fédérale de Lausanne, Lausanne
1015, Switzerland
| | - Louise Ormond
- School of Life Sciences, École Polytechnique Fédérale de Lausanne, Lausanne
1015, Switzerland
| | - Claudia Bank
- Instituto Gulbenkian de Ciência, Oeiras 2780-156, Portugal
| | - Jeffrey D. Jensen
- School of Life Sciences, École Polytechnique Fédérale de Lausanne, Lausanne
1015, Switzerland
- Center for Evolution and Medicine, School of Life Sciences, Arizona State
University, Tempe, AZ 85287, USA
| |
Collapse
|
37
|
Bierne N, Bonhomme F, Arnaud-Haond S. Dedicated population genomics for the silent world: the specific questions of marine population genetics. Curr Zool 2016; 62:545-550. [PMID: 29491944 PMCID: PMC5804263 DOI: 10.1093/cz/zow107] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Nicolas Bierne
- Université de Montpellier, Montpellier, France
- CNRS – Institut des Sciences de l’Evolution Montpellier, UMR 5554 UM – CNRS
– IRD – EPHE, Station Marine OREME, Sète, France
| | - François Bonhomme
- Université de Montpellier, Montpellier, France
- CNRS – Institut des Sciences de l’Evolution Montpellier, UMR 5554 UM – CNRS
– IRD – EPHE, Station Marine OREME, Sète, France
| | - Sophie Arnaud-Haond
- Université de Montpellier, Montpellier, France
- Ifremer – MARine Biodiversity, Exploitation and Conservation, UMR 9190 IRD –
IFREMER – UM – CNRS, Sète, France
| |
Collapse
|