Population size changes and selection drive patterns of parallel evolution in a host-virus system

High parasite diversity maintained after an alga-virus coevolutionary arms race

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.

Incomplete immunity in a natural animal-microbiota interaction selects for higher pathogen virulence

Incomplete immunity in recovered hosts is predicted to favor more virulent pathogens upon re-infection in the population.1 The microbiota colonizing animals can generate a similarly long-lasting, partial immune response, allowing for infection but dampened disease severity.2 We tracked the evolutionary trajectories of a widespread pathogen (Pseudomonas aeruginosa), experimentally passaged through populations of nematodes immune-primed by a natural microbiota member (P. berkeleyensis). This bacterium can induce genes regulated by a mitogen-activated protein kinase (MAPK) signaling pathway effective at conferring protection against pathogen-induced death despite infection.3 Across host populations, this incomplete immunity selected for pathogens more than twice as likely to kill as those evolved in non-primed (i.e., naive) or immune-compromised (mutants with a knockout of the MAPK ortholog) control populations. Despite the higher virulence, pathogen molecular evolution in immune-primed hosts was slow and constrained. In comparison, evolving pathogens in immune-compromised hosts were characterized by substantial genomic differentiation and attenuated virulence. These findings directly attribute the incomplete host immunity induced from microbiota as a significant force shaping the virulence and evolutionary dynamics of novel infectious diseases.

Species and condition shape the mutational spectrum in experimentally evolved biofilms

IMPORTANCE

Biofilm formation is a vital factor for the survival and adaptation of bacteria in diverse environmental niches. Experimental evolution combined with the advancement of whole-population genome sequencing provides us a powerful tool to understand the genomic dynamic of evolutionary adaptation to different environments, such as during biofilm development. Previous studies described the genetic and phenotypic changes of selected clones from experimentally evolved Bacillus thuringiensis and Bacillus subtilis that were adapted under abiotic and biotic biofilm conditions. However, the full understanding of the dynamic evolutionary landscapes was lacking. Furthermore, the differences and similarities of adaptive mechanisms in B. thuringiensis and B. subtilis were not identified. To overcome these limitations, we performed longitudinal whole-population genome sequencing to study the underlying genetic dynamics at high resolution. Our study provides the first comprehensive mutational landscape of two bacterial species' biofilms that is adapted to an abiotic and biotic surface.

Evolutionary repeatability of emergent properties of ecological communities

Most species belong to ecological communities where their interactions give rise to emergent community-level properties, such as diversity and productivity. Understanding and predicting how these properties change over time has been a major goal in ecology, with important practical implications for sustainability and human health. Less attention has been paid to the fact that community-level properties can also change because member species evolve. Yet, our ability to predict long-term eco-evolutionary dynamics hinges on how repeatably community-level properties change as a result of species evolution. Here, we review studies of evolution of both natural and experimental communities and make the case that community-level properties at least sometimes evolve repeatably. We discuss challenges faced in investigations of evolutionary repeatability. In particular, only a handful of studies enable us to quantify repeatability. We argue that quantifying repeatability at the community level is critical for approaching what we see as three major open questions in the field: (i) Is the observed degree of repeatability surprising? (ii) How is evolutionary repeatability at the community level related to repeatability at the level of traits of member species? (iii) What factors affect repeatability? We outline some theoretical and empirical approaches to addressing these questions. Advances in these directions will not only enrich our basic understanding of evolution and ecology but will also help us predict eco-evolutionary dynamics. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Characterizing viral within-host diversity in fast and non-equilibrium demo-genetic dynamics

High-throughput sequencing has opened the route for a deep assessment of within-host genetic diversity that can be used, e.g., to characterize microbial communities and to infer transmission links in infectious disease outbreaks. The performance of such characterizations and inferences cannot be analytically assessed in general and are often grounded on computer-intensive evaluations. Then, being able to simulate within-host genetic diversity across time under various demo-genetic assumptions is paramount to assess the performance of the approaches of interest. In this context, we built an original model that can be simulated to investigate the temporal evolution of genotypes and their frequencies under various demo-genetic assumptions. The model describes the growth and the mutation of genotypes at the nucleotide resolution conditional on an overall within-host viral kinetics, and can be tuned to generate fast non-equilibrium demo-genetic dynamics. We ran simulations of this model and computed classic diversity indices to characterize the temporal variation of within-host genetic diversity (from high-throughput amplicon sequences) of virus populations under three demographic kinetic models of viral infection. Our results highlight how demographic (viral load) and genetic (mutation, selection, or drift) factors drive variations in within-host diversity during the course of an infection. In particular, we observed a non-monotonic relationship between pathogen population size and genetic diversity, and a reduction of the impact of mutation on diversity when a non-specific host immune response is activated. The large variation in the diversity patterns generated in our simulations suggests that the underlying model provides a flexible basis to produce very diverse demo-genetic scenarios and test, for instance, methods for the inference of transmission links during outbreaks.

Predicting cancer evolution for patient benefit: Renal cell carcinoma paradigm

Evolutionary features of cancer have important clinical implications, but their evaluation in the clinic is currently limited. Here, we review current approaches to reconstruct tumour subclonal structure and discuss tumour sampling method and experimental design influence. We describe clear-cell renal cell carcinoma (ccRCC) as an exemplar for understanding and predicting cancer evolutionary dynamics. Finally, we discuss how understanding cancer evolution can benefit patients.

Neural networks enable efficient and accurate simulation-based inference of evolutionary parameters from adaptation dynamics

The rate of adaptive evolution depends on the rate at which beneficial mutations are introduced into a population and the fitness effects of those mutations. The rate of beneficial mutations and their expected fitness effects is often difficult to empirically quantify. As these 2 parameters determine the pace of evolutionary change in a population, the dynamics of adaptive evolution may enable inference of their values. Copy number variants (CNVs) are a pervasive source of heritable variation that can facilitate rapid adaptive evolution. Previously, we developed a locus-specific fluorescent CNV reporter to quantify CNV dynamics in evolving populations maintained in nutrient-limiting conditions using chemostats. Here, we use CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based likelihood-free inference approaches. We tested the suitability of 2 evolutionary models: a standard Wright-Fisher model and a chemostat model. We evaluated 2 likelihood-free inference algorithms: the well-established Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC) algorithm, and the recently developed Neural Posterior Estimation (NPE) algorithm, which applies an artificial neural network to directly estimate the posterior distribution. By systematically evaluating the suitability of different inference methods and models, we show that NPE has several advantages over ABC-SMC and that a Wright-Fisher evolutionary model suffices in most cases. Using our validated inference framework, we estimate the CNV formation rate at the GAP1 locus in the yeast Saccharomyces cerevisiae to be 10-4.7 to 10-4 CNVs per cell division and a fitness coefficient of 0.04 to 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our inference-based estimates using 2 distinct experimental methods-barcode lineage tracking and pairwise fitness assays-which provide independent confirmation of the accuracy of our approach. Our results are consistent with a beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining the outsized importance of CNVs in rapid adaptive evolution. More generally, our study demonstrates the utility of novel neural network-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data with possible applications ranging from tumor to viral evolution.

Towards an integrative view of virus phenotypes

Understanding how phenotypes emerge from genotypes is a foundational goal in biology. As challenging as this task is when considering cellular life, it is further complicated in the case of viruses. During replication, a virus as a discrete entity (the virion) disappears and manifests itself as a metabolic amalgam between the virus and the host (the virocell). Identifying traits that unambiguously constitute a virus's phenotype is straightforward for the virion, less so for the virocell. Here, we present a framework for categorizing virus phenotypes that encompasses both virion and virocell stages and considers functional and performance traits of viruses in the context of fitness. Such an integrated view of virus phenotype is necessary for comprehensive interpretation of viral genome sequences and will advance our understanding of viral evolution and ecology.

Functional comparison of MERS-coronavirus lineages reveals increased replicative fitness of the recombinant lineage 5

Middle East respiratory syndrome coronavirus (MERS-CoV) is enzootic in dromedary camels across the Middle East and Africa. Virus-induced pneumonia in humans results from animal contact, with a potential for limited onward transmission. Phenotypic changes have been suspected after a novel recombinant clade (lineage 5) caused large nosocomial outbreaks in Saudi Arabia and South Korea in 2016. However, there has been no functional assessment. Here we perform a comprehensive in vitro and ex vivo comparison of viruses from parental and recombinant virus lineages (lineage 3, n = 7; lineage 4, n = 8; lineage 5, n = 9 viruses) from Saudi Arabia, isolated immediately before and after the shift toward lineage 5. Replication of lineage 5 viruses is significantly increased. Transcriptional profiling finds reduced induction of immune genes IFNB1, CCL5, and IFNL1 in lung cells infected with lineage 5 strains. Phenotypic differences may be determined by IFN antagonism based on experiments using IFN receptor knock out and signaling inhibition. Additionally, lineage 5 is more resilient against IFN pre-treatment of Calu-3 cells (ca. 10-fold difference in replication). This phenotypic change associated with lineage 5 has remained undiscovered by viral sequence surveillance, but may be a relevant indicator of pandemic potential.

MERS-CoV is enzootic in dromedary camels, can spread to humans but undergoes limited onward transmission. Here, Schroeder et al. compare clinical isolates of MERS-CoV in vitro and show that the predominantly circulating recombinant lineage 5 possess a fitness advantage over parental lineage 3 and 4 due to reduced activation of innate immune signaling.

Next-generation cophylogeny: unravelling eco-evolutionary processes

A fundamental question in evolutionary biology is how microevolutionary processes translate into species diversification. Cophylogeny provides an appropriate framework to address this for symbiotic associations, but historically has been primarily limited to unveiling patterns. We argue that it is essential to integrate advances from ecology and evolutionary biology into cophylogeny, to gain greater mechanistic insights and transform cophylogeny into a platform to advance understanding of interspecific interactions and diversification more widely. We discuss key directions, such as incorporating trait reconstruction and considering multiple scales of network organization, and highlight recent developments for implementation. A new quantitative framework is proposed to allow integration of relevant information, such as quantitative traits and assessment of the contribution of individual mechanisms to cophylogenetic patterns.

Novel genomic approaches to study antagonistic coevolution between hosts and parasites

Host-parasite coevolution is ubiquitous, shaping genetic and phenotypic diversity and the evolutionary trajectory of interacting species. With the advances of high throughput sequencing technologies applicable to model and non-model organisms alike, it is now feasible to study in greater detail (a) the genetic underpinnings of coevolution, (b) the speed and type of dynamics at coevolving loci, and (c) the genomic consequences of coevolution. This review focuses on three recently developed approaches that leverage information from host and parasite full genome data simultaneously to pinpoint coevolving loci and draw inference on the coevolutionary history. First, co-genome-wide association study (co-GWAS) methods allow pinpointing the loci underlying host-parasite interactions. These methods focus on detecting associations between genetic variants and the outcome of experimental infection tests or on correlations between genomes of naturally infected hosts and their infecting parasites. Second, extensions to population genomics methods can detect genes under coevolution and infer the coevolutionary history, such as fitness costs. Third, correlations between host and parasite population size in time are indicative of coevolution, and polymorphism levels across independent spatially distributed populations of hosts and parasites can reveal coevolutionary loci and infer coevolutionary history. We describe the principles of these three approaches and discuss their advantages and limitations based on coevolutionary theory. We present recommendations for their application to various host (prokaryotes, fungi, plants, and animals) and parasite (viruses, bacteria, fungi, and macroparasites) species. We conclude by pointing out methodological and theoretical gaps to be filled to extract maximum information from full genome data and thereby to shed light on the molecular underpinnings of coevolution.

Neutral genomic signatures of host-parasite coevolution

Background

Coevolution is a selective process of reciprocal adaptation in hosts and parasites or in mutualistic symbionts. Classic population genetics theory predicts the signatures of selection at the interacting loci of both species, but not the neutral genome-wide polymorphism patterns. To bridge this gap, we build an eco-evolutionary model, where neutral genomic changes over time are driven by a single selected locus in hosts and parasites via a simple biallelic gene-for-gene or matching-allele interaction. This coevolutionary process may lead to cyclic changes in the sizes of the interacting populations.

Results

We investigate if and when these changes can be observed in the site frequency spectrum of neutral polymorphisms from host and parasite full genome data. We show that changes of the host population size are too smooth to be observable in its polymorphism pattern over the course of time. Conversely, the parasite population may undergo a series of strong bottlenecks occurring on a slower relative time scale, which may lead to observable changes in a time series sample. We also extend our results to cases with 1) several parasites per host accelerating relative time, and 2) multiple parasite generations per host generation slowing down rescaled time.

Conclusions

Our results show that time series sampling of host and parasite populations with full genome data are crucial to understand if and how coevolution occurs. This model provides therefore a framework to interpret and draw inference from genome-wide polymorphism data of interacting species.

The feedback between selection and demography shapes genomic diversity during coevolution

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.

The phylogenomics of evolving virus virulence

How virulence evolves after a virus jumps to a new host species is central to disease emergence. Our current understanding of virulence evolution is based on insights drawn from two perspectives that have developed largely independently: long-standing evolutionary theory based on limited real data examples that often lack a genomic basis, and experimental studies of virulence-determining mutations using cell culture or animal models. A more comprehensive understanding of virulence mutations and their evolution can be achieved by bridging the gap between these two research pathways through the phylogenomic analysis of virus genome sequence data as a guide to experimental study.

Population size changes and extinction risk of populations driven by mutant interactors

Spontaneous random mutations are an important source of variation in populations. Many evolutionary models consider mutants with a fixed fitness, chosen from a fitness distribution without considering microscopic interactions among the residents and mutants. Here, we go beyond this and consider "mutant interactors," which lead to new interactions between the residents and invading mutants that can affect the population size and the extinction risk of populations. We model microscopic interactions between individuals by using a dynamic interaction matrix, the dimension of which increases with the emergence of a new mutant and decreases with extinction. The new interaction parameters of the mutant follow a probability distribution around the payoff entries of its ancestor. These new interactions can drive the population away from the previous equilibrium and lead to changes in the population size. Thus, the population size is an evolving property rather than an externally controlled variable. We calculate the average population size of our stochastic system over time and quantify the extinction risk of the population by the mean time to extinction.

Ecological and Evolutionary Processes Shaping Viral Genetic Diversity

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.

Construction and Characterization of Synthetic Bacterial Community for Experimental Ecology and Evolution

Experimental microbial ecology and evolution have yielded foundational insights into ecological and evolutionary processes using simple microcosm setups and phenotypic assays with one- or two-species model systems. The fields are now increasingly incorporating more complex systems and exploration of the molecular basis of observations. For this purpose, simplified, manageable and well-defined multispecies model systems are required that can be easily investigated using culturing and high-throughput sequencing approaches, bridging the gap between simpler and more complex synthetic or natural systems. Here we address this need by constructing a completely synthetic 33 bacterial strain community that can be cultured in simple laboratory conditions. We provide whole-genome data for all the strains as well as metadata about genomic features and phenotypic traits that allow resolving individual strains by amplicon sequencing and facilitate a variety of envisioned mechanistic studies. We further show that a large proportion of the strains exhibit coexistence in co-culture over serial transfer for 48 days in the absence of any experimental manipulation to maintain diversity. The constructed bacterial community can be a valuable resource in future experimental work.

Population size changes and selection drive patterns of parallel evolution in a host-virus system

Predicting the repeatability of evolution remains elusive. Theory and empirical studies suggest that strong selection and large population sizes increase the probability for parallel evolution at the phenotypic and genotypic levels. However, selection and population sizes are not constant, but rather change continuously and directly affect each other even on short time scales. Here, we examine the degree of parallel evolution shaped through eco-evolutionary dynamics in an algal host population coevolving with a virus. We find high degrees of parallelism at the level of population size changes (ecology) and at the phenotypic level between replicated populations. At the genomic level, we find evidence for parallelism, as the same large genomic region was duplicated in all replicated populations, but also substantial novel sequence divergence between replicates. These patterns of genome evolution can be explained by considering population size changes as an important driver of rapid evolution.