The effects of migration and drift on local adaptation to a heterogeneous environment

Local adaptation experiments are widely used to quantify the levels of adaptation within a heterogeneous environment. However, theoretical studies generally focus on the probability of fixation of alleles or the mean fitness of populations, rather than local adaptation as it is commonly measured experimentally or in field studies. Here, we develop mathematical models and use them to generate analytical predictions for the level of local adaptation as a function of selection, migration and genetic drift. First, we contrast mean fitness and local adaptation measures and show that the latter can be expressed in a simple and general way as a function of the spatial covariance between population mean phenotype and local environmental conditions. Second, we develop several approximations of a population genetics model to show that the system exhibits different behaviours depending on the rate of migration. The main insights are the following: with intermediate migration, both genetic drift and migration decrease local adaptation; with low migration, drift decreases local adaptation but migration speeds up adaptation; with high migration, genetic drift has no effect on local adaptation. Third, we extend this analysis to cases where the trait under selection is continuous using classical quantitative genetics theory. Finally, we discuss these results in the light of recent experimental work on local adaptation.     

Tissue electroporation: quantification and analysis of heterogeneous transport in multicellular environments

Although electroporation is gaining increased attention as a technology to enhance clinical chemotherapy and gene therapy of tissues, direct measurements of electroporation-mediated transport in multicellular environments are lacking. In this study, we used multicellular tumor spheroids of DU145 prostate cancer cells as a model tissue to measure the levels and distribution of molecular uptake in a multicellular environment as a function of electrical and other parameters. These measurements, and subsequent analysis, were used to test the hypothesis that cells in a multicellular environment respond to electroporation in a heterogeneous manner that differs from isolated cells in suspension due to differences in cell state, local solute concentration, and local electric field. In support of the hypothesis, molecular uptake was consistently lower for cells within spheroids than cells in dilute suspension and was spatially heterogeneous, with progressively less uptake observed for cells located deeper within spheroid interiors. Reduced uptake and heterogeneity can be explained quantitatively by accounting for the effects of cell size on transmembrane voltage and cell volume, limited extracellular solute reservoir, heterogeneous field strength due to influence of neighboring cells, and diffusional lag times.     

Parasites promote host gene flow in a metapopulation

Local adaptation is a powerful mechanism to maintain genetic diversity in subdivided populations. It counteracts the homogenizing effect of gene flow because immigrants have an inferior fitness in the new habitat. This picture may be reversed in host populations where parasites influence the success of immigrating hosts. Here we report two experiments testing whether parasite abundance and genetic background influences the success of host migration among pools in a Daphnia magna metapopulation. In 22 natural populations of D. magna, immigrant hosts were found to be on average more successful when the resident populations experienced high prevalences of a local microsporidian parasite. We then determined whether this success is due to parasitism per se, or the genetic background of the parasites. In a common garden competition experiment, we found that parasites reduced the fitness of their local hosts relatively more than the fitness of allopatric host genotypes. Our experiments are consistent with theoretical predictions based on coevolutionary host-parasite models in metapopulations. A direct consequence of the observed mechanism is an elevated effective migration rate for the host in the metapopulation.     

Compartmentalization of surface envelope glycoprotein of human immunodeficiency virus type 1 during acute and chronic infection

Human immunodeficiency virus type 1 is characterized by extensive genetic heterogeneity. Having previously demonstrated that, in the peripheral blood, the initial viral population is more homogeneous than at subsequent stages of infection, we have extended our studies to tissue samples, allowing comparisons between viral populations in peripheral blood and tissues during both the acute and chronic stages of infection. We found that homogeneity in gp120 sequences during the acute infection phase is not just restricted to the peripheral blood but also extends to other tissue compartments. However, in chronically infected individuals, heterogeneous and distinct viral populations were found in different compartments. We therefore conclude that the dominant and homogeneous viral population observed during the acute infection phase is likely to infiltrate lymphoid tissues and form the genetic bases for subsequent diversification. It is therefore likely that the compartmentalization of viral sequences observed in chronically infected patients reflects a gradual diversification of a common dominant viral variant rather than the preferential migration of distinct viral populations to different tissue compartments at the beginning of infection.     

Local adaptation, resistance, and virulence in a hemiparasitic plant-host plant interaction

Abstract.— Coevolution may lead to local adaptation of parasites to their sympatric hosts. Locally adapted parasites are, on average, more infectious to sympatric hosts than to allopatric hosts of the same species or their fitness on the sympatric hosts is superior to that on allopatric hosts. We tested local adaptation of a hemiparasitic plant, Rhinanthus serotinus (Scrophulariaceae), to its host plant, the grass Agrostis capillaris . Using a reciprocal cross-infection experiment, we exposed host plants from four sites to hemiparasites originating from the same four sites in a common environment. The parasites were equally able to establish haustorial connections to sympatric and allopatric hosts, and their performance was similar on both host types. Therefore, these results do not indicate local adaptation of the parasites to their sympatric hosts. However, the parasite populations differed in average biomass and number of flowers per plant and in their effect on host biomass. These results indicate that the virulence of the parasite varied among populations, suggesting genetic variation. Theoretical models suggest that local adaptation is likely to be detected if the host and the parasite have different evolutionary potentials, different migration rates, and the parasite is highly virulent. In the interaction between R. serotinus and A. capillaris all the theoretical prerequisites for local adaptation may not be fulfilled.     

Evolutionary feedback mediated through population density, illustrated with viruses in chemostats

A cornerstone of evolutionary ecology is that population density affects adaptation: r and K selection is the obvious example. The reverse is also appreciated: adaptation impacts population density. Yet, empirically demonstrating a direct connection between population density and adaptation is challenging. Here, we address both evolution and ecology of population density in models of viral (bacteriophage) chemostats. Chemostats supply nutrients for host cell growth, and the hosts are prey for viral reproduction. Two different chemostat designs have profoundly different consequences for viral evolution. If host and virus are confined to the same chamber, as in a predator-prey system, viral regulation of hosts feeds back to maintain low viral density (measured as infections per cell). Viral adaptation impacts host density but has a small effect on equilibrium viral density. More interesting are chemostats that supply the viral population with hosts from a virus-free refuge. Here, a type of evolutionary succession operates: adaptation at low viral density leads to higher density, but high density then favors competitive ability. Experiments support these models with both phenotypic and molecular data. Parallels to these designs exist in many natural systems, so these experimental systems may yield insights to the evolution and regulation of natural populations.     

A synthesis of experimental work on parasite local adaptation

The study of parasite local adaptation, whereby parasites perform better on sympatric hosts than on allopatric hosts and/or better on their own host population than do other parasites, is of great importance to both basic and applied biology. Theoretical examination of host-parasite coevolution predicts that parasite migration rate, generation time and virulence all contribute to the pattern of parasite local adaptation, such that parasites with greater dispersal ability, more frequent reproduction and/or high virulence ought to exhibit increased infectivity on local hosts. Here, we present a meta-analysis of experimental work from 57 host-parasite systems across 54 local adaptation studies to directly test theoretical predictions concerning the effect of each attribute on parasite adaptation. As expected, we find that studies of parasites with higher migration rates than their hosts report local adaptation, as measured by infection success, significantly more often than studies of parasites with relatively low migration rates. Furthermore, this synthesis serves to identify biases in the current body of work and highlight areas with the greatest need for further study. We emphasize the importance of unifying the field with regard to experimental methods, local adaptation definitions and reported statistics for cross-infection studies.     

Viral infection in internally structured hosts. I. Conditions for persistent infection

For a virus population within its host, two important levels of structure can be considered: multiple cell types which can be infected, and tissue types or body compartments which may be coupled via movement. We develop a model with both types of structure. Migration between compartments can create "sources" and "sinks" within the virus population, where realized viral growth rate and abundance is lowered in some compartments compared to what would be observed in isolation. Using both analytical and numerical methods, we investigate how this within-host spatial structure affects the conditions for persistent viral infection. We find that migration between compartments makes the establishment of infection more difficult than it would be in the absence of migration, implying that within-host spatial structure combined with viral movement decreases the likelihood of viral establishment. If migration is symmetrical and compartments are heterogeneous, an increase in migration rates between compartments generally makes establishment less likely. This may help to explain the tissue specificity observed for many viruses. There are, however, important exceptions to this result. These include circumstances where the virus initially invades a compartment that is unfavorable to population growth and migration is necessary to infect other parts of the host body. Stochastic aspects of viral establishment may also favor increased migration as it tends to dampen the amplitude of fluctuations in population size during the initial transient phase of establishment.     

Vesicular stomatitis virus evolution during alternation between persistent infection in insect cells and acute infection in mammalian cells is dominated by the persistence phase

Vesicular stomatitis virus has the potential for very rapid evolution in the laboratory, but like many other arboviruses, it evolves at a relatively slow rate in the natural environment. Previous work showed that alternating replication in different cell types does not promote stasis. In order to determine whether other factors promote stasis, we compared the fitness trajectories of populations evolving during acute infections in mammalian cells, populations evolving during persistent infections in insect cells, and populations evolving during alternating acute and persistent infection cycles. Populations evolving under constant conditions increased in fitness in the environment in which they replicated. An asymmetric trade-off was observed such that acute infection had no cost for persistence but persistent replication had a dramatic cost for acute infection in mammalian cells. After an initial period of increase, fitness remained approximately constant in all the populations that included persistent replication, but fitness continuously increased in populations evolving during acute infections. Determination of the consensus sequence of the genes encoding the N, P, M, and G proteins showed that the pattern of mutation accumulation was coherent with fitness changes during persistence so that once fitness reached a maximum, the rate of mutation accumulation dropped. Persistent replication dominated both the genetic and the phenotypic evolution of the populations that alternated between acute infection of mammalian cells and persistence in insect cells, and fitness loss was observed in the mammalian environment despite periodic replication in mammalian cells. These results show that stasis can be achieved without good levels of adaptation to both the mammalian and the insect environments.     

Do mites evolving in alternating host plants adapt to host switch?

A fluctuating environment may be perceived as a composition of different environments, or as an environment per se, in which it is the fluctuation itself that poses a selection pressure. If so, then organisms may adapt to this alternation. We tested this using experimental populations of spider mites that have been evolving for 45 generations in a homogeneous environment (pepper or tomato plants), or in a heterogeneous environment composed of an alternation of these two plants approximately at each generation. The performance (daily oviposition rate and juvenile survival) of individuals from these populations was tested in each of the homogeneous environments, and in two alternating environments, one every 3 days and the other between generations. To discriminate between potential genetic interactions between alleles conferring adaptation to each host plant and environmental effects of evolving in a fluctuating environment, we compared the performance of all lines with that of a cross between tomato and pepper lines. As a control, two lines within each selection regime were also crossed. We found that crosses between alternating lines and between pepper and tomato lines performed worse than crosses between lines evolving in homogeneous environments when tested in that environment. In contrast, alternating lines performed either better or similarly to lines evolving in homogeneous environments when tested in a fluctuating environment. Our results suggest that fluctuating environments are more than the juxtaposition of two environments. Hence, tests for adaptation of organisms evolving in such environments should be carried out in fluctuating conditions.     

BIO-LGCA: A cellular automaton modelling class for analysing collective cell migration

Collective dynamics in multicellular systems such as biological organs and tissues plays a key role in biological development, regeneration, and pathological conditions. Collective tissue dynamics—understood as population behaviour arising from the interplay of the constituting discrete cells—can be studied with on- and off-lattice agent-based models. However, classical on-lattice agent-based models, also known as cellular automata, fail to replicate key aspects of collective migration, which is a central instance of collective behaviour in multicellular systems. To overcome drawbacks of classical on-lattice models, we introduce an on-lattice, agent-based modelling class for collective cell migration, which we call biological lattice-gas cellular automaton (BIO-LGCA). The BIO-LGCA is characterised by synchronous time updates, and the explicit consideration of individual cell velocities. While rules in classical cellular automata are typically chosen ad hoc, rules for cell-cell and cell-environment interactions in the BIO-LGCA can also be derived from experimental cell migration data or biophysical laws for individual cell migration. We introduce elementary BIO-LGCA models of fundamental cell interactions, which may be combined in a modular fashion to model complex multicellular phenomena. We exemplify the mathematical mean-field analysis of specific BIO-LGCA models, which allows to explain collective behaviour. The first example predicts the formation of clusters in adhesively interacting cells. The second example is based on a novel BIO-LGCA combining adhesive interactions and alignment. For this model, our analysis clarifies the nature of the recently discovered invasion plasticity of breast cancer cells in heterogeneous environments.     

MALADAPTATION AS A SOURCE OF SENESCENCE IN HABITATS VARIABLE IN SPACE AND TIME

In this study, we use a quantitative genetics model of structured populations to investigate the evolution of senescence in a variable environment. Adaptation to local environments depends on phenotypic traits whose optimal values vary with age and with environmental conditions. We study different scenarios of environmental heterogeneity, where the environment changes abruptly, gradually, or cyclically with time and where the environment is heterogeneous in space with different populations connected by migration. The strength of selection decreases with age, which predicts slower adaptation of traits expressed late in the life cycle, potentially generating stronger senescence in habitats where selection changes in space or in time. This prediction is however complicated by the fact that the genetic variance also increases with age. Using numerical calculations, we found that the rate of senescence is generally increased when the environment varies. In particular, migration between different habitats is a source of senescence in heterogeneous landscapes. We also show that the rate of senescence can vary transiently when the population is not at equilibrium, with possible implications for experimental evolution and the study of invasive species. Our results highlight the need to study age‐specific adaptation, as a changing environment can have a different impact on different age classes.     

Eph/ephrin signaling in cancer

The Eph receptor tyrosine kinases family and their membrane bound ligands, the ephrins, represents a complex signaling network of cell communication for cell sorting during tissue patterning in development and in the normal physiology and homeostasis of adult tissues. This molecular family has adapted to evolving tissue complexity in multicellular organisms through the emergence of more members and complex mechanisms of expression and signaling that result in the fine-tuning of cell positioning. Since their initial identification from an erythropoietin producing hepatocellular (Eph) carcinoma cell line in 1987, Eph/ephrin signaling has been a matter of intensive investigation for their plausible role in cancer. Similarly to their context dependent modus operandi in normal tissues, Eph/ephrin signaling in cancer is an intricate and puzzling network of events that tumors “manage” to their benefit in multiple aspects like cell adhesion to substrate, migration, invasion or growth.     

ROLE OF SEX AND MIGRATION IN ADAPTATION TO SINK ENVIRONMENTS

Understanding the effects of sex and migration on adaptation to novel environments remains a key problem in evolutionary biology. Using a single‐cell alga Chlamydomonas reinhardtii, we investigated how sex and migration affected rates of evolutionary rescue in a sink environment, and subsequent changes in fitness following evolutionary rescue. We show that sex and migration affect both the rate of evolutionary rescue and subsequent adaptation. However, their combined effects change as the populations adapt to a sink habitat. Both sex and migration independently increased rates of evolutionary rescue, but the effect of sex on subsequent fitness improvements, following initial rescue, changed with migration, as sex was beneficial in the absence of migration but constraining adaptation when combined with migration. These results suggest that sex and migration are beneficial during the initial stages of adaptation, but can become detrimental as the population adapts to its environment.     

The impact of migration from parasite-free patches on antagonistic host-parasite coevolution

Natural populations of hosts and parasites are often subdivided and patchily distributed such that some regions of a host species' range will be free from a given parasite. Host migration from parasite-free to parasite-containing patches is expected to alter coevolutionary dynamics by changing the evolutionary potential of antagonists. Specifically, host immigration can favor parasites by increasing transmission opportunities, or hosts by introducing genetic variation. We tested these predictions in coevolving populations of Pseudomonas fluorescens and phage Phi2 that received immigrants from phage-free populations. We observed a negative quadratic relationship between sympatric resistance to phage and host immigration rate (highest at intermediate immigration) but a positive quadratic relationship between coevolution rate and host immigration rate (lowest at intermediate immigration). These results indicate that for a wide range of rates, host immigration from parasite-free patches can increase the evolutionary potential of parasites, and increase the coevolutionary rate if parasite adaptation is limiting in the absence of immigration.     

Models of RNA virus evolution and their roles in vaccine design

Viruses are fast evolving pathogens that continuously adapt to the highly variable environments they live and reproduce in. Strategies devoted to inhibit virus replication and to control their spread among hosts need to cope with these extremely heterogeneous populations and with their potential to avoid medical interventions. Computational techniques such as phylogenetic methods have broadened our picture of viral evolution both in time and space, and mathematical modeling has contributed substantially to our progress in unraveling the dynamics of virus replication, fitness, and virulence. Integration of multiple computational and mathematical approaches with experimental data can help to predict the behavior of viral pathogens and to anticipate their escape dynamics. This piece of information plays a critical role in some aspects of vaccine development, such as viral strain selection for vaccinations or rational attenuation of viruses. Here we review several aspects of viral evolution that can be addressed quantitatively, and we discuss computational methods that have the potential to improve vaccine design.     

Viral biogeography revealed by signatures in Sulfolobus islandicus genomes

Viruses are a driving force of microbial evolution. Despite their importance, the evolutionary dynamics that shape diversity in viral populations are not well understood. One of the primary factors that define viral population structure is coevolution with microbial hosts. Experimental models predict that the trajectory of coevolution will be determined by the relative migration rates of viruses and their hosts; however, there are no natural microbial systems in which both have been examined. The biogeographic distribution of viruses that infect Sulfolobus islandicus is investigated using genome comparisons among four newly identified, integrated, Sulfolobus spindle-shaped viruses and previously sequenced viral strains. Core gene sequences show a biogeographic distribution where viral genomes are specifically associated with each local population. In addition, signatures of host–virus interactions recorded in the sequence-specific CRISPR (clustered regularly interspaced short palindromic repeats) system show that hosts have interacted with viral communities that are more closely related to local viral strains than to foreign ones. Together, both proviral and CRISPR sequences show a clear biogeographic structure for Sulfolobus viral populations. Our findings demonstrate that virus–microbe coevolution must be examined in a spatially explicit framework. The combination of host and virus biogeography suggests a model for viral diversification driven by host immunity and local adaptation.     

Estimation of the rate and effect of new beneficial mutations in asexual populations

The rate and effect of available beneficial mutations are key parameters in determining how a population adapts to a new environment. However, these parameters are poorly known, in large part because of the difficulty of designing and interpreting experiments to examine the rare and intrinsically stochastic process of mutation occurrence. We present a new approach to estimate the rate and selective advantage of beneficial mutations that underlie the adaptation of asexual populations. We base our approach on the analysis of experiments that track the effect of newly arising beneficial mutations on the dynamics of a neutral marker in evolving bacterial populations and develop efficient estimators of mutation rate and selective advantage. Using extensive simulations, we evaluate the accuracy of our estimators and conclude that they are quite robust to the use of relatively low experimental replication. To validate the predictions of our model, we compare theoretical and experimentally determined estimates of the selective advantage of the first beneficial mutation to fix in a series of ten replicate populations. We find that our theoretical predictions are not significantly different from experimentally determined selection coefficients. Application of our method to suitably designed experiments will allow estimation of how population evolvability depends on demographic and initial fitness parameters.     

Lack of evolutionary stasis during alternating replication of an arbovirus in insect and mammalian cells

The evolution of vesicular stomatitis virus (VSV) in a constant environment, consisting of either mammalian or insect cells, has been compared to the evolution of the same viral population in changing environments consisting in alternating passages in mammalian and insect cells. Fitness increases were observed in all cases. An initial fitness loss of VSV passaged in insect cells was noted when fitness was measured in BHK-21 cells, but this effect could be attributed to a difference of temperature during VSV replication at 37 degrees C in BHK-21 cells. Sequencing of nucleotides 1-4717 at the 3' end of the VSV genome (N, P, M and G genes) showed that at passage 80 the number of mutations accumulated during alternated passages (seven mutations) is similar or larger than that observed in populations evolving in a constant environment (two to four mutations). Our results indicate that insect and mammalian cells can constitute similar environments for viral replication. Thus, the slow rates of evolution observed in natural populations of arboviruses are not necessarily due to the need for the virus to compromise between adaptation to both arthropod and vertebrate cell types.     

Evolution of polygenic traits under global vs local adaptation

Observations about the number, frequency, effect size, and genomic distribution of alleles associated with complex traits must be interpreted in light of evolutionary process. These characteristics, which constitute a trait’s genetic architecture, can dramatically affect evolutionary outcomes in applications from agriculture to medicine, and can provide a window into how evolution works. Here, I review theoretical predictions about the evolution of genetic architecture under spatially homogeneous, global adaptation as compared with spatially heterogeneous, local adaptation. Due to the tension between divergent selection and migration, local adaptation can favor “concentrated” genetic architectures that are enriched for alleles of larger effect, clustered in a smaller number of genomic regions, relative to expectations under global adaptation. However, the evolution of such architectures may be limited by many factors, including the genotypic redundancy of the trait, mutation rate, and temporal variability of environment. I review the circumstances in which predictions differ for global vs local adaptation and discuss where progress can be made in testing hypotheses using data from natural populations and lab experiments. As the field of comparative population genomics expands in scope, differences in architecture among traits and species will provide insights into how evolution works, and such differences must be interpreted in light of which kind of selection has been operating.