Predicting evolution in experimental range expansions of an aquatic model system.

Predicting range expansion dynamics is an important goal of both fundamental and applied research in conservation and global change biology. However, this is challenging if ecological and evolutionary processes occur on the same time scale. Using the freshwater ciliate Paramecium caudatum, we combined experimental evolution and mathematical modeling to assess the predictability of evolutionary change during range expansions. In the experiment, we followed ecological dynamics and trait evolution in independently replicated microcosm populations in range core and front treatments, where episodes of natural dispersal alternated with periods of population growth. These eco-evolutionary conditions were recreated in a predictive mathematical model, parametrized with dispersal and growth data of the 20 founder strains in the experiment. We found that short-term evolution was driven by selection for increased dispersal in the front treatment and general selection for higher growth rates in all treatments. There was a good quantitative match between predicted and observed trait changes. Phenotypic divergence was further mirrored by genetic divergence between range core and front treatments. In each treatment, we found the repeated fixation of the same cytochrome c oxidase I (COI) marker genotype, carried by strains that also were the most likely winners in our model. Long-term evolution in the experimental range front lines resulted in the emergence of a dispersal syndrome, namely a competition-colonization trade-off. Altogether, both model and experiment highlight the potential importance of dispersal evolution as a driver of range expansions. Thus, evolution at range fronts may follow predictable trajectories, at least for simple scenarios, and predicting these dynamics may be possible from knowledge of few key parameters.

Travelling with a parasite: the evolution of resistance and dispersal syndromes during experimental range expansion.

Rapid evolutionary change during range expansions can lead to diverging range core and front populations, with the emergence of dispersal syndromes (coupled responses in dispersal and life-history traits). Besides intraspecific effects, range expansions may be impacted by interspecific interactions such as parasitism. Yet, despite the potentially large impact of parasites imposing additional selective pressures on the host, their role on range expansions remains largely unexplored. Using microcosm populations of the ciliate Paramecium caudatum and its bacterial parasite Holospora undulata, we studied experimental range expansions under parasite presence or absence. We found that the interaction of range expansion and parasite treatments affected the evolution of host dispersal syndromes. Namely, front populations showed different associations of population growth parameters and swimming behaviours than core populations, indicating divergent evolution. Parasitism reshaped trait associations, with hosts evolved in the presence of the parasite exhibiting overall increased resistance and reduced dispersal. Nonetheless, when comparing infected range core and front populations, we found a positive association, suggesting joint evolution of resistance and dispersal at the front. We conclude that host-parasite interactions during range expansions can change evolutionary trajectories; this in turn may feedback on the ecological dynamics of the range expansion and parasite epidemics.

Spatial autocorrelation of local patch extinctions drives recovery dynamics in metacommunities.

Human activities put ecosystems under increasing pressure, often resulting in local extinctions. However, it is unclear how local extinctions affect regional processes, such as the distribution of diversity in space, especially if extinctions show spatial patterns, such as being clustered. Therefore, it is crucial to investigate extinctions and their consequences in a spatially explicit framework. Using highly controlled microcosm experiments and theoretical models, we ask here how the number and spatial autocorrelation of extinctions interactively affect metacommunity dynamics. We found that local patch extinctions increased local diversity (α-diversity) and inter-patch diversity (ß-diversity) by delaying the exclusion of inferior competitors. Importantly, recolonization dynamics depended more strongly on the spatial distribution than on the number of patch extinctions: clustered local patch extinctions resulted in slower recovery, lower α-diversity and higher ß-diversity. Our results highlight that the spatial distribution of perturbations should be taken into account when studying and managing spatially structured communities.

Niche breadth affects bacterial transcription patterns along a salinity gradient.

Understanding the molecular mechanisms that determine a species' life history is important for predicting their susceptibility to environmental change. While specialist species with a narrow niche breadth (NB) maximize their fitness in their optimum habitat, generalists with broad NB adapt to multiple environments. The main objective of this study was to identify general transcriptional patterns that would distinguish bacterial strains characterized by contrasted NBs along a salinity gradient. More specifically, we hypothesized that genes encoding fitness-related traits, such as biomass production, have a higher degree of transcriptional regulation in specialists than in generalists, because the fitness of specialists is more variable under environmental change. By contrast, we expected that generalists would exhibit enhanced transcriptional regulation of genes encoding traits that protect them against cellular damage. To test these hypotheses, we assessed the transcriptional regulation of fitness-related and adaptation-related genes of 11 bacterial strains in relation to their NB and stress exposure under changing salinity conditions. The results suggested that transcriptional regulation levels of fitness- and adaptation-related genes correlated with the NB and/or the stress exposure of the inspected strains. We further identified a shortlist of candidate stress marker genes that could be used in future studies to monitor the susceptibility of bacterial populations or communities to environmental changes.

Parasitism and host dispersal plasticity in an aquatic model system.

Dispersal is a central determinant of spatial dynamics in communities and ecosystems, and various ecological factors can shape the evolution of constitutive and plastic dispersal behaviours. One important driver of dispersal plasticity is the biotic environment. Parasites, for example, influence the internal condition of infected hosts and define external patch quality. Thus, state-dependent dispersal may be determined by infection status and context-dependent dispersal by the abundance of infected hosts in the population. A prerequisite for such dispersal plasticity to evolve is a genetic basis on which natural selection can act. Using interconnected microcosms, we investigated dispersal in experimental populations of the freshwater protist Paramecium caudatum in response to the bacterial parasite Holospora undulata. For a collection of 20 natural host strains, we found substantial variation in constitutive dispersal and to a lesser degree in dispersal plasticity. First, infection tended to increase or decrease dispersal relative to uninfected controls, depending on strain identity, indicative of state-dependent dispersal plasticity. Infection additionally decreased host swimming speed compared to the uninfected counterparts. Second, for certain strains, there was a weak negative association between dispersal and infection prevalence, such that uninfected hosts dispersed less when infection was more frequent in the population, indicating context-dependent dispersal plasticity. Future experiments may test whether the observed differences in dispersal plasticity are sufficiently strong to be picked up by natural selection. The evolution of dispersal plasticity as a strategy to mitigate parasite effects spatially may have important implications for epidemiological dynamics.

An evolutionary trade-off between parasite virulence and dispersal at experimental invasion fronts.

Exploitative parasites are predicted to evolve in highly connected populations or in expanding epidemics. However, many parasites rely on host dispersal to reach new populations, potentially causing conflict between local transmission and global spread. We performed experimental range expansions in interconnected microcosms of the protozoan Paramecium caudatum, allowing natural dispersal of hosts infected with the bacterial parasite Holospora undulata. Parasites from range front treatments facilitated host dispersal and were less virulent, but also invested less in horizontal transmission than parasites from range cores. These differences were consistent with parameter estimates derived from an epidemiological model fitted on population-level time-series data. Our results illustrate how dispersal selection can have profound consequences for the evolution of parasite life history and virulence. Decrypting the eco-evolutionary processes that shape parasite 'dispersal syndromes' may be important for the management of spreading epidemics in changing environments, biological invasions or in other spatial non-equilibrium settings.

Phage steering of antibiotic-resistance evolution in the bacterial pathogen, Pseudomonas aeruginosa.

BACKGROUND AND OBJECTIVES: Antimicrobial resistance is a growing global concern and has spurred increasing efforts to find alternative therapeutics. Bacteriophage therapy has seen near constant use in Eastern Europe since its discovery over a century ago. One promising approach is to use phages that not only reduce bacterial pathogen loads but also select for phage resistance mechanisms that trade-off with antibiotic resistance-so called 'phage steering'. METHODOLOGY: Recent work has shown that the phage OMKO1 can interact with efflux pumps and in so doing select for both phage resistance and antibiotic sensitivity of the pathogenic bacterium Pseudomonas aeruginosa. We tested the robustness of this approach to three different antibiotics in vitro (tetracycline, erythromycin and ciprofloxacin) and one in vivo (erythromycin). RESULTS: We show that in vitro OMKO1 can reduce antibiotic resistance of P. aeruginosa (Washington PAO1) even in the presence of antibiotics, an effect still detectable after ca.70 bacterial generations in continuous culture with phage. Our in vivo experiment showed that phage both increased the survival times of wax moth larvae (Galleria mellonella) and increased bacterial sensitivity to erythromycin. This increased antibiotic sensitivity occurred both in lines with and without the antibiotic. CONCLUSIONS AND IMPLICATIONS: Our study supports a trade-off between antibiotic resistance and phage sensitivity. This trade-off was maintained over co-evolutionary time scales even under combined phage and antibiotic pressure. Similarly, OMKO1 maintained this trade-off in vivo, again under dual phage/antibiotic pressure. Our findings have implications for the future clinical use of steering in phage therapies. Lay Summary: Given the rise of antibiotic-resistant bacterial infection, new approaches to treatment are urgently needed. Bacteriophages (phages) are bacterial viruses. The use of such viruses to treat infections has been in near-continuous use in several countries since the early 1900s. Recent developments have shown that these viruses are not only effective against routine infections but can also target antibiotic resistant bacteria in a novel, unexpected way. Similar to other lytic phages, these so-called 'steering phages' kill the majority of bacteria directly. However, steering phages also leave behind bacterial variants that resist the phages, but are now sensitive to antibiotics. Treatment combinations of these phages and antibiotics can now be used to greater effect than either one independently. We evaluated the impact of steering using phage OMKO1 and a panel of three antibiotics on Pseudomonas aeruginosa, an important pathogen in hospital settings and in people with cystic fibrosis. Our findings indicate that OMKO1, either alone or in combination with antibiotics, maintains antibiotic sensitivity both in vitro and in vivo, giving hope that phage steering will be an effective treatment option against antibiotic-resistant bacteria.

Transient negative effects of antibiotics on phages do not jeopardise the advantages of combination therapies.

Phages, the viruses of bacteria, have been proposed as antibacterial agents to complement or replace antibiotics due to the growing problem of resistance. In nature and in the clinic, antibiotics are ubiquitous and may affect phages indirectly via impacts on bacterial hosts. Even if the synergistic association of phages and antibiotics has been shown in several studies, the focus is often on bacteria with little known about the impact on phages. Evolutionary studies have demonstrated that time scale is an important factor in understanding the consequences of antimicrobial strategies, but this perspective is generally overlooked in phage-antibiotic combination studies. Here, we explore the effects of antibiotics on phages targeting the opportunistic pathogen Pseudomonas aeruginosa. We go beyond previous studies by testing the interaction between several types of antibiotics and phages, and evaluate the effects on several important phage parameters during 8 days of experimental co-evolution with bacteria. Our study reveals that antibiotics had a negative effect on phage density and efficacy early on, but not in the later stages of the experiment. The results indicate that antibiotics can affect phage adaptation, but that phages can nevertheless contribute to managing antibiotic resistance levels.

Evolutionary rescue and local adaptation under different rates of temperature increase: a combined analysis of changes in phenotype expression and genotype frequency in Paramecium microcosms.

Evolutionary rescue (ER) occurs when populations, which have declined due to rapid environmental change, recover through genetic adaptation. The success of this process and the evolutionary trajectory of the population strongly depend on the rate of environmental change. Here we investigated how different rates of temperature increase (from 23 to 32 °C) affect population persistence and evolutionary change in experimental microcosms of the protozoan Paramecium caudatum. Consistent with theory on ER, we found that those populations experiencing the slowest rate of temperature increase were the least likely to become extinct and tended to be the best adapted to the new temperature environment. All high-temperature populations were more tolerant to severe heat stress (35, 37 °C), indicating a common mechanism of heat protection. High-temperature populations also had superior growth rates at optimum temperatures, leading to the absence of a pattern of local adaptation to control (23 °C) and high-temperature (32 °C) environments. However, high-temperature populations had reduced growth at low temperatures (5-9 °C), causing a shift in the temperature niche. In part, the observed evolutionary change can be explained by selection from standing variation. Using mitochondrial markers, we found complete divergence between control and high-temperature populations in the frequencies of six initial founder genotypes. Our results confirm basic predictions of ER and illustrate how adaptation to an extreme local environment can produce positive as well as negative correlated responses to selection over the entire range of the ecological niche.

Network structure and local adaptation in co-evolving bacteria-phage interactions.

Numerous theoretical and experimental studies have investigated antagonistic co-evolution between parasites and their hosts. Although experimental tests of theory from a range of biological systems are largely concordant regarding the influence of several driving processes, we know little as to how mechanisms acting at the smallest scales (individual molecular and phenotypic changes) may result in the emergence of structures at larger scales, such as co-evolutionary dynamics and local adaptation. We capitalized on methods commonly employed in community ecology to quantify how the structure of community interaction matrices, so-called bipartite networks, reflected observed co-evolutionary dynamics, and how phages from these communities may or may not have adapted locally to their bacterial hosts. We found a consistent nested network structure for two phage types, one previously demonstrated to exhibit arms race co-evolutionary dynamics and the other fluctuating co-evolutionary dynamics. Both phages increased their host ranges through evolutionary time, but we found no evidence for a trade-off with impact on bacteria. Finally, only bacteria from the arms race phage showed local adaptation, and we provide preliminary evidence that these bacteria underwent (sometimes different) molecular changes in the wzy gene associated with the LPS receptor, while bacteria co-evolving with the fluctuating selection phage did not show local adaptation and had partial deletions of the pilF gene associated with type IV pili. We conclude that the structure of phage-bacteria interaction networks is not necessarily specific to co-evolutionary dynamics, and discuss hypotheses for why only one of the two phages was, nevertheless, locally adapted.

Antibiotic stress selects against cooperation in the pathogenic bacterium Pseudomonas aeruginosa.

Cheats are a pervasive threat to public goods production in natural and human communities, as they benefit from the commons without contributing to it. Although ecological antagonisms such as predation, parasitism, competition, and abiotic environmental stress play key roles in shaping population biology, it is unknown how such stresses generally affect the ability of cheats to undermine cooperation. We used theory and experiments to address this question in the pathogenic bacterium, Pseudomonas aeruginosa Although public goods producers were selected against in all populations, our competition experiments showed that antibiotics significantly increased the advantage of nonproducers. Moreover, the dominance of nonproducers in mixed cultures was associated with higher resistance to antibiotics than in either monoculture. Mathematical modeling indicates that accentuated costs to producer phenotypes underlie the observed patterns. Mathematical analysis further shows how these patterns should generalize to other taxa with public goods behaviors. Our findings suggest that explaining the maintenance of cooperative public goods behaviors in certain natural systems will be more challenging than previously thought. Our results also have specific implications for the control of pathogenic bacteria using antibiotics and for understanding natural bacterial ecosystems, where subinhibitory concentrations of antimicrobials frequently occur.

Size evolution in microorganisms masks trade-offs predicted by the growth rate hypothesis.

Adaptation to local resource availability depends on responses in growth rate and nutrient acquisition. The growth rate hypothesis (GRH) suggests that growing fast should impair competitive abilities for phosphorus and nitrogen due to high demand for biosynthesis. However, in microorganisms, size influences both growth and uptake rates, which may mask trade-offs and instead generate a positive relationship between these traits (size hypothesis, SH). Here, we evolved a gradient of maximum growth rate (µmax) from a single bacterium ancestor to test the relationship among µmax, competitive ability for nutrients and cell size, while controlling for evolutionary history. We found a strong positive correlation between µmax and competitive ability for phosphorus, associated with a trade-off between µmax and cell size: strains selected for high µmax were smaller and better competitors for phosphorus. Our results strongly support the SH, while the trade-offs expected under GRH were not apparent. Beyond plasticity, unicellular populations can respond rapidly to selection pressure through joint evolution of their size and maximum growth rate. Our study stresses that physiological links between these traits tightly shape the evolution of competitive strategies.

Predicting biotic interactions and their variability in a changing environment.

Global environmental change is altering the patterns of biodiversity worldwide. Observation and theory suggest that species' distributions and abundances depend on a suite of processes, notably abiotic filtering and biotic interactions, both of which are constrained by species' phylogenetic history. Models predicting species distribution have historically mostly considered abiotic filtering and are only starting to integrate biotic interaction. However, using information on present interactions to forecast the future of biodiversity supposes that biotic interactions will not change when species are confronted with new environments. Using bacterial microcosms, we illustrate how biotic interactions can vary along an environmental gradient and how this variability can depend on the phylogenetic distance between interacting species.

Long-term selection experiment produces breakdown of horizontal transmissibility in parasite with mixed transmission mode.

Evolutionary transitions from parasitism toward beneficial or mutualistic associations may encompass a change from horizontal transmission to (strict) vertical transmission. Parasites with both vertical and horizontal transmission are amendable to study factors driving such transitions. In a long-term experiment, microcosm populations of the protozoan Paramecium caudatum and its bacterial parasite Holospora undulata were exposed to three growth treatments, manipulating vertical transmission opportunities over ca. 800 host generations. In inoculation tests, horizontal transmission propagules produced by parasites from a "high-growth" treatment, with elevated host division rates increasing levels of parasite vertical transmission, showed a near-complete loss of infectivity. A similar reduction was observed for parasites from a treatment alternating between high growth and low growth (i.e., low levels of population turn-over). Parasites from a low-growth treatment had the highest infectivity on all host genotypes tested. Our results complement previous findings of reduced investment in horizontal transmission and increased vertical transmissibility of high-growth parasites. We explain the loss of horizontal transmissibility by epidemiological feedbacks and resistance evolution, reducing the frequency of susceptible hosts in the population and thereby decreasing the selective advantage of horizontal transmission. This illustrates how environmental conditions may push parasites with a mixed transmission mode toward becoming vertically transmitted nonvirulent symbionts.

Competition-colonization dynamics in experimental bacterial metacommunities.

One of the simplest hypotheses used to explain species coexistence is the competition-colonization trade-off, that is, species can stably coexist in a landscape if they show a trade-off between competitive and colonization abilities. Despite extensive theory, the dynamics predicted to result from competition-colonization trade-offs are largely untested. Landscape change, such as habitat destruction, is thought to greatly influence coexistence under competition-colonization dynamics, although there is no formal test of this prediction. Here we present the first illustration of competition-colonization dynamics that fully transposes theory into a controlled experimental metacommunity of two Pseudomonas bacterial strains. The competition-colonization dynamics were achieved by directly manipulating trade-off strength and colonization rates to generate the full range of coexistence conditions and responses to habitat destruction. Our study successfully generates competition-colonization dynamics matching theoretical predictions, and our results further reveal a negative relationship between diversity and productivity when scaling up to entire metacommunities.

Phylogenetic constraints on ecosystem functioning.

There is consensus that biodiversity losses will result in declining ecosystem functioning if species have different functional traits. Phylogenetic diversity has recently been suggested as a predictor of ecosystem functioning because it could approximate the functional complementarity among species. Here we describe an experiment that takes advantage of the rapid evolutionary response of bacteria to disentangle the role of phylogenetic and species diversity. We impose a strong selection regime on marine bacterial lineages and assemble the ancestral and evolved lines in microcosms of varying lineage and phylogenetic diversity. We find that the relationship between phylogenetic diversity and productivity is strong for the ancestral lineages but brakes down for the evolved lineages. Our results not only emphasize the potential of using phylogeny to evaluate ecosystem functioning, but also they warn against using phylogenetics as a proxy for functional diversity without good information on species evolutionary history.

Evolutionary dynamics of separate and combined exposure of Pseudomonas fluorescens SBW25 to antibiotics and bacteriophage.

The use of bacteriophages against pathogenic bacteria in health care and in the food industry is now being advocated as an alternative to the use of antibiotics. But what is the evolutionary response for a bacterial population if both antibiotics and phages are used in combination? We employ an experimental evolution approach to address these questions and exposed Pseudomonas fluorescens SBW25 and a related hypermutator strain (mutS-) to the action of the antibiotic rifampicin and the lytic bacteriophage SBW25ϕ2. We then compared the densities, growth rates, and the mutations at the rpoB locus leading to rifampicin resistance of the evolved bacterial populations. We observed that the evolutionary response of populations under different treatments varied depending on the order in which the antimicrobials were added and whether the bacterium was a hypermutator. We found that wild-type rifampicin-resistant populations involved in biofilm formation often reverted to rifampicin sensitivity when stresses were added sequentially. In contrast, when the mortality agents were added simultaneously, phage populations frequently went extinct and the bacteria evolved antibiotic resistance. However, populations of the hypermutator mutS- converged to a single genotype at the rpoB locus. Future investigation on other bacteria and using different antibiotics and bacteriophage are needed to evaluate the generality of our findings.

Contrasted effects of diversity and immigration on ecological insurance in marine bacterioplankton communities.

The ecological insurance hypothesis predicts a positive effect of species richness on ecosystem functioning in a variable environment. This effect stems from temporal and spatial complementarity among species within metacommunities coupled with optimal levels of dispersal. Despite its importance in the context of global change by human activities, empirical evidence for ecological insurance remains scarce and controversial. Here we use natural aquatic bacterial communities to explore some of the predictions of the spatial and temporal aspects of the ecological insurance hypothesis. Addressing ecological insurance with bacterioplankton is of strong relevance given their central role in fundamental ecosystem processes. Our experimental set up consisted of water and bacterioplankton communities from two contrasting coastal lagoons. In order to mimic environmental fluctuations, the bacterioplankton community from one lagoon was successively transferred between tanks containing water from each of the two lagoons. We manipulated initial bacterial diversity for experimental communities and immigration during the experiment. We found that the abundance and production of bacterioplankton communities was higher and more stable (lower temporal variance) for treatments with high initial bacterial diversity. Immigration was only marginally beneficial to bacterial communities, probably because microbial communities operate at different time scales compared to the frequency of perturbation selected in this study, and of their intrinsic high physiologic plasticity. Such local "physiological insurance" may have a strong significance for the maintenance of bacterial abundance and production in the face of environmental perturbations.

Ecological strategies shape the insurance potential of biodiversity.

Biodiversity is thought to provide insurance for ecosystem functioning under heterogeneous environments; however, such insurance potential is under serious threat following unprecedented loss of biodiversity. One of the key mechanism underlying ecological insurance is that niche differentiation allows asynchronous responses to fluctuating environments, although the role of different ecological strategies (e.g., specialists vs. generalists) has yet to be formally evaluated. We present here a simple experimental study that illustrates how different ecological strategies (i.e., generalists vs. specialists) can shape the biodiversity-insurance relationship. We assembled microcosm of generalists and specialist bacteria over a gradient of salinity and found that, bacterial communities made up of generalists were more productive and more stable over time under environmental fluctuations. We discuss our results in context with simple theoretical predictions and propose future directions for biological insurance theory. We argue that beyond species richness itself, it is essential to incorporate the distribution of ecological strategies across relevant environmental gradients as predictors of the insurance potential of biodiversity in natural ecosystems.

Terminal investment induced by a bacteriophage in a rhizosphere bacterium.

Despite knowledge about microbial responses to abiotic stress, few studies have investigated stress responses to antagonistic species, such as competitors, predators and pathogens. While it is often assumed that interacting populations of bacteria and phage will coevolve resistance and exploitation strategies, an alternative is that individual bacteria tolerate or evade phage predation through inducible responses to phage presence. Using the microbial model Pseudomonas fluorescens SBW25 and its lytic DNA phage SBW25Φ2, we demonstrate the existence of an inducible response in the form of a transient increase in population growth rate, and found that the response was induced by phage binding. This response was accompanied by a decrease in bacterial cell size, which we propose to be an associated cost. We discuss these results in the context of bacterial ecology and phage-bacteria co-evolution.