Lysogeny destabilizes computationally simulated microbiomes

Microbiomes are ecosystems, and their stability can impact the health of their hosts. Theory predicts that predators influence ecosystem stability. Phages are key predators of bacteria in microbiomes, but phages are unusual predators because many have lysogenic life cycles. It has been hypothesized that lysogeny can destabilize microbiomes, but lysogeny has no direct analog in classical ecological theory, and no formal theory exists. We studied the stability of computationally simulated microbiomes with different numbers of temperate (lysogenic) and virulent (obligate lytic) phage species. Bacterial populations were more likely to fluctuate over time when there were more temperate phages species. After disturbances, bacterial populations returned to their pre-disturbance densities more slowly when there were more temperate phage species, but cycles engendered by disturbances dampened more slowly when there were more virulent phage species. Our work offers the first formal theory linking lysogeny to microbiome stability.

Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Population genetics of Artemia urmiana species complex (Crustacea, Anostraca): A group with asymmetrical dispersal and gene flow mediated by migratory waterfowl

Bird-mediated dispersal of resting eggs is the main mechanism for Artemia dispersal among catchments. The bisexual populations of Artemia urmiana species complex, which is here considered to be a collection of Artemia genetically close to the so-called "Western Asian Lineage", are mostly distributed in central and western Asia (i.e., in regions falling into the Central Asian Flyway of migratory birds) and live in diversified habitats. Little is known about the genetic relationships among these populations. Aiming to understand the population genetic characteristics and the roles of migratory birds on the dispersal and gene flow of this Artemia group, we evaluated the genetic diversity, genetic differentiation, and gene flow among 14 populations, with their altitudes ranging from 540 to 4870 m above sea level, using 13 microsatellite markers. Almost all populations exhibited high genetic diversity and heterozygote excess, which may be a consequence of combined effects of dispersal and hybridization. The global genetic differentiation (FST) value was 0.092, the pairwise FST values were 0.003-0.246. Discriminant analysis of principal components identified three genetic clusters, consisting of Urmia Lake (Iran), Zhundong (Xinjiang, China), and 12 Qinghai-Tibet Plateau populations, respectively. The among-population genetic differentiation seems to be a consequence of isolation by distance and adaptation to diversified habitats induced by altitudinal gradient. Historical gene flows are asymmetrical, and show an evolutionary source-sink dynamics, with Jingyu Lake (Xinjiang, China) population being the major source. These results support our hypothesis that in Qinghai-Tibet Plateau and surrounding areas the bird-mediated dispersal of Artemia may be biased towards from north to south and/or from higher altitude to lower altitude.

Understanding the joint evolution of dispersal and host specialisation using phytophagous arthropods as a model group

Theory generally predicts that host specialisation and dispersal should evolve jointly. Indeed, many models predict that specialists should be poor dispersers to avoid landing on unsuitable hosts while generalists will have high dispersal abilities. Phytophagous arthropods are an excellent group to test this prediction, given extensive variation in their host range and dispersal abilities. Here, we explore the degree to which the empirical literature on this group is in accordance with theoretical predictions. We first briefly outline the theoretical reasons to expect such a correlation. We then report empirical studies that measured both dispersal and the degree of specialisation in phytophagous arthropods. We find a correlation between dispersal and levels of specialisation in some studies, but with wide variation in this result. We then review theoretical attributes of species and environment that may blur this correlation, namely environmental grain, temporal heterogeneity, habitat selection, genetic architecture, and coevolution between plants and herbivores. We argue that theoretical models fail to account for important aspects, such as phenotypic plasticity and the impact of selective forces stemming from other biotic interactions, on both dispersal and specialisation. Next, we review empirical caveats in the study of this interplay. We find that studies use different measures of both dispersal and specialisation, hampering comparisons. Moreover, several studies do not provide independent measures of these two traits. Finally, variation in these traits may occur at scales that are not being considered. We conclude that this correlation is likely not to be expected from large-scale comparative analyses as it is highly context dependent and should not be considered in isolation from the factors that modulate it, such as environmental scale and heterogeneity, intrinsic traits or biotic interactions. A stronger crosstalk between theoretical and empirical studies is needed to understand better the prevalence and basis of the correlation between dispersal and specialisation.

Host-Parasite Coevolution in Continuous Space Leads to Variation in Local Adaptation across Spatial Scales

AbstractPrevious host-parasite coevolutionary theory has focused on understanding the determinants of local adaptation using spatially discrete models. However, these studies fall short of describing patterns of host-parasite local adaptation across spatial scales. In contrast, empirical work demonstrates that patterns of adaptation depend on the scale at which they are measured. Here, we propose a mathematical model of host-parasite coevolution in continuous space that naturally leads to a scale-dependent definition of local adaptation. In agreement with empirical research, we find that patterns of adaptation vary across spatial scales. In some cases, not only the magnitude of local adaptation but also the identity of the locally adapted species will depend on the spatial scale at which measurements are taken. Building on our results, we suggest a way to consistently measure parasite local adaptation when continuous space is the driver of cross-scale variation. We also describe a way to test whether continuous space is driving cross-scale variation. Taken together, our results provide a new perspective that can be used to understand empirical observations previously unexplained by theoretical expectations and deepens our understanding of the mechanics of host-parasite local adaptation.

The maintenance of genetic diversity under host-parasite coevolution in finite, structured populations

As a corollary to the Red Queen hypothesis, host-parasite coevolution has been hypothesized to maintain genetic variation in both species. Recent theoretical work, however, suggests that reciprocal natural selection alone is insufficient to maintain variation at individual loci. As highlighted by our brief review of the theoretical literature, models of host-parasite coevolution often vary along multiple axes (e.g. inclusion of ecological feedbacks or abiotic selection mosaics), complicating a comprehensive understanding of the effects of interacting evolutionary processes on diversity. Here we develop a series of comparable models to explore the effect of interactions between spatial structures and antagonistic coevolution on genetic diversity. Using a matching alleles model in finite populations connected by migration, we find that, in contrast to panmictic populations, coevolution in a spatially structured environment can maintain genetic variation relative to neutral expectations with migration alone. These results demonstrate that geographic structure is essential for understanding the effect of coevolution on biological diversity.

How do fluctuating ecological dynamics impact the evolution of hosts and parasites?

Theoretical models of the evolution of parasites and their hosts have shaped our understanding of infectious disease dynamics for over 40 years. Many theoretical models assume that the underlying ecological dynamics are at equilibrium or constant, yet we know that in a great many systems there are fluctuations in the ecological dynamics owing to a variety of intrinsic or extrinsic factors. Here, we discuss the challenges presented when modelling evolution in systems with fluctuating ecological dynamics and summarize the main approaches that have been developed to study host-parasite evolution in such systems. We provide an in-depth guide to one of the methods by applying it to two worked examples of host evolution that have not previously been studied in the literature: when cycles occur owing to seasonal forcing in competition, and when the presence of a free-living parasite causes cycles, with accompanying interactive Python code provided. We review the findings of studies that have explored host-parasite evolution when ecological dynamics fluctuate, and point to areas of future research. Throughout we stress the importance of feedbacks between the ecological and evolutionary dynamics in driving the outcomes of infectious disease systems. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

The impact of habitat loss on molecular signatures of coevolution between an iconic butterfly (Alcon blue) and its host plant (Marsh gentian)

Habitat loss is threatening natural communities worldwide. Small and isolated populations suffer from inbreeding and genetic drift, which jeopardize their long-term survival and adaptive capacities. However, the consequences of habitat loss for reciprocal coevolutionary interactions remain poorly studied. In this study, we investigated the effects of decreasing habitat patch size and connectivity associated with habitat loss on molecular signatures of coevolution in the Alcon blue butterfly (Phengaris alcon) and its most limited host, the marsh gentian (Gentiana pneumonanthe). Because reciprocal coevolution is characterized by negative frequency-dependent selection as a particular type of balancing selection, we investigated how signatures of balancing selection vary along a gradient of patch size and connectivity, using single nucleotide polymorphisms (SNPs). We found that signatures of coevolution were unaffected by patch characteristics in the host plants. On the other hand, more pronounced signatures of coevolution were observed in both spatially isolated and in large Alcon populations, together with pronounced spatial variation in SNPs that are putatively involved in coevolution. These findings suggest that habitat loss can facilitate coevolution in large butterfly populations through limiting swamping of locally beneficial alleles by maladaptive ones. We also found that allelic richness (Ar) of the coevolutionary SNPs is decoupled from neutral Ar in the butterfly, indicating that habitat loss has different effects on coevolutionary as compared with neutral processes. We conclude that this specialized coevolutionary system requires particular conservation interventions aiming at generating a spatial mosaic of both connected and of isolated habitat to maintain coevolutionary dynamics.

Effects of habitat destruction on coevolving metacommunities

Habitat destruction is a growing threat to biodiversity and ecosystem services. The ecological consequences of habitat loss and fragmentation involve reductions in species abundance and even the extinction of species and their interactions. However, we do not yet understand how habitat loss alters the coevolutionary trajectories of the remaining species or how coevolution, in turn, affects their response to habitat loss. To investigate this, we develop a spatially explicit model which couples metacommunity and coevolutionary dynamics. We show that, by changing the size, composition and structure of local networks, habitat destruction increases the diversity of coevolutionary trajectories of mutualists across the landscape. Conversely, in antagonistic communities, some species increase while others reduce their spatial trait heterogeneity. Furthermore, we show that while coevolution dampens the negative effects of habitat destruction in mutualistic networks, its effects on the persistence of antagonistic communities tend to be smaller and less predictable.

Large-scale and small-scale population genetic structure of the medically important gastropod species Bulinus truncatus (Gastropoda, Heterobranchia)

Background

Gastropod snails remain strongly understudied, despite their important role in transmitting parasitic diseases. Knowledge of their distribution and population dynamics increases our understanding of the processes driving disease transmission. We report the first study to use high-throughput sequencing (HTS) to elucidate the population genetic structure of the hermaphroditic snail Bulinus truncatus (Gastropoda, Heterobranchia) on a regional (17–150 km) and inter-regional (1000–5400 km) scale. This snail species acts as an intermediate host of Schistosoma haematobium and Schistosoma bovis, which cause human and animal schistosomiasis respectively.

Methods

Bulinus truncatus snails were collected in Senegal, Cameroon, Egypt and France and identified through DNA barcoding. A single-end genotyping-by-sequencing (GBS) library, comprising 87 snail specimens from the respective countries, was built and sequenced on an Illumina HiSeq 2000 platform. Reads were mapped against S. bovis and S. haematobium reference genomes to identify schistosome infections, and single nucleotide polymorphisms (SNPs) were scored using the Stacks pipeline. These SNPs were used to estimate genetic diversity, assess population structure and construct phylogenetic trees of B. truncatus.

Results

A total of 10,750 SNPs were scored and used in downstream analyses. The phylogenetic analysis identified five clades, each consisting of snails from a single country but with two distinct clades within Senegal. Genetic diversity was low in all populations, reflecting high selfing rates, but varied between locations due to habitat variability. Significant genetic differentiation and isolation by distance patterns were observed at both spatial scales, indicating that gene flow is not strong enough to counteract the effects of population bottlenecks, high selfing rates and genetic drift. Remarkably, the population genetic differentiation on a regional scale (i.e. within Senegal) was as large as that between populations on an inter-regional scale. The blind GBS technique was able to pick up parasite DNA in snail tissue, demonstrating the potential of HTS techniques to further elucidate the role of snail species in parasite transmission.

Conclusions

HTS techniques offer a valuable toolbox to further investigate the population genetic patterns of intermediate schistosome host snails and the role of snail species in parasite transmission.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05445-x.

Phylogeography and Cryptic Species Structure of a Locally Adapted Parasite in New Zealand

Phylogeographic patterns of many taxa on New Zealand's South Island are characterised by disjunct distributions that have been attributed to Pleistocene climatic cycles and the formation of the Southern Alps. Pleistocene glaciation has also been implicated in shaping the contemporary genetic differentiation between populations of the aquatic snail Potamopyrgus antipodarum. We investigated whether similar phylogeographic patterns exist for the snail's locally adapted trematode parasite, Atriophallophorus winterbourni. We found evidence for a barrier to gene-flow in sympatry between cryptic, but ecologically divergent species. When focusing on the most common of these species, disjunct geographic distributions are found for mitochondrial lineages that diverged during the Pleistocene. The boundary between these distributions is found in the central part of the South Island and reinforced by low cross-alpine migration. Further support for a vicariant origin of the phylogeographic pattern was found when assessing nuclear multilocus SNP data. Nuclear and mitochondrial population differentiation was concordant in pattern, except for populations in a potential secondary contact zone. Additionally, we found larger than expected differentiation between nuclear- and mitochondrial-based empirical Bayes F_ST estimates (global F_ST : 0.02 versus 0.39 for nuclear and mitochondrial data, respectively). Population subdivision is theoretically expected to be stronger for mitochondrial genomes due to a smaller effective population size, but the strong difference here, together with mito-nuclear discordance in a putative contact zone, is potentially indicative of divergent gene flow of nuclear and mitochondrial genomes.

Antagonistic Coevolution Limits the Range of Host Defense in C. elegans Populations

Host populations often evolve defenses against parasites due to the significant fitness costs imposed by infection. However, adaptation to a specific parasite may alter the effectiveness of the host’s defenses in general. Consequently, the specificity of host defense may be influenced by a host population’s evolutionary history with parasites. Further, the degree of reciprocal change within an interaction may profoundly alter the range of host defense, given that antagonistic coevolutionary interactions are predicted to favor defense against specific parasite genotypes. Here, we examined the effect of host evolutionary history on host defense range by assessing the mortality rates of Caenorhabditis elegans host populations exposed to an array of Serratia marcescens bacterial parasite strains. Importantly, each of the host populations were derived from the same genetic background but have different experimental evolution histories with parasites. Each of these histories (exposure to either heat-killed, fixed genotype, or coevolving parasites) carries a different level of evolutionary reciprocity. Overall, we observed an effect of host evolutionary history in that previously coevolved host populations were generally the most susceptible to novel parasite strains. This data demonstrates that host evolutionary history can have a significant impact on host defense, and that host-parasite coevolution can increase host susceptibility to novel parasites.

Coevolutionary theory of hosts and parasites

Host and parasite evolution are closely intertwined, with selection for adaptations and counter‐adaptations forming a coevolutionary feedback loop. Coevolutionary dynamics are often difficult to intuit due to these feedbacks and are hard to demonstrate empirically in most systems. Theoretical models have therefore played a crucial role in shaping our understanding of host–parasite coevolution. Theoretical models vary widely in their assumptions, approaches and aims, and such variety makes it difficult, especially for non‐theoreticians and those new to the field, to: (1) understand how model approaches relate to one another; (2) identify key modelling assumptions; (3) determine how model assumptions relate to biological systems; and (4) reconcile the results of different models with contrasting assumptions. In this review, we identify important model features, highlight key results and predictions and describe how these pertain to model assumptions. We carry out a literature survey of theoretical studies published since the 1950s (n = 219 papers) to support our analysis. We identify two particularly important features of models that tend to have a significant qualitative impact on the outcome of host–parasite coevolution: population dynamics and the genetic basis of infection. We also highlight the importance of other modelling features, such as stochasticity and whether time proceeds continuously or in discrete steps, that have received less attention but can drastically alter coevolutionary dynamics. We finish by summarizing recent developments in the field, specifically the trend towards greater model complexity, and discuss likely future directions for research.

Co-structure analysis and genetic associations reveal insights into pinworms (Trypanoxyuris) and primates (Alouatta palliata) microevolutionary dynamics

BACKGROUND

In parasitism arm race processes and red queen dynamics between host and parasites reciprocally mold many aspects of their genetics and evolution. We performed a parallel assessment of population genetics and demography of two species of pinworms with different degrees of host specificity (Trypanoxyuris multilabiatus, species-specific; and T. minutus, genus-specific) and their host, the mantled howler monkey (Alouatta palliata), based on mitochondrial DNA sequences and microsatellite loci (these only for the host). Given that pinworms and primates have a close co-evolutionary history, covariation in several genetic aspects of their populations is expected.

RESULTS

Mitochondrial DNA revealed two genetic clusters (West and East) in both pinworm species and howler monkeys, although population structure and genetic differentiation were stronger in the host, while genetic diversity was higher in pinworms than howler populations. Co-divergence tests showed no congruence between host and parasite phylogenies; nonetheless, a significant correlation was found between both pinworms and A. palliata genetic pairwise distances suggesting that the parasites' gene flow is mediated by the host dispersal. Moreover, the parasite most infective and the host most susceptible haplotypes were also the most frequent, whereas the less divergent haplotypes tended to be either more infective (for pinworms) or more susceptible (for howlers). Finally, a positive correlation was found between pairwise p-distance of host haplotypes and that of their associated pinworm haplotypes.

CONCLUSION

The genetic configuration of pinworm populations appears to be molded by their own demography and life history traits in conjunction with the biology and evolutionary history of their hosts, including host genetic variation, social interactions, dispersal and biogeography. Similarity in patterns of genetic structure, differentiation and diversity is higher between howler monkeys and T. multilabiatus in comparison with T. minutus, highlighting the role of host-specificity in coevolving processes. Trypanoxyuris minutus exhibits genetic specificity towards the most frequent host haplotype as well as geographic specificity. Results suggest signals of potential local adaptation in pinworms and further support the notion of correlated evolution between pinworms and their primate hosts.

Spatial structure and dispersal dynamics in a house sparrow metapopulation

The effects of spatial structure on metapopulation dynamics depend upon the interaction between local population dynamics and dispersal, and how this relationship is affected by the geographical isolation and spatial heterogeneity in habitat characteristics. Our aim is to examine how emigration and immigration of house sparrows Passer domesticus in a Norwegian archipelagic metapopulation are affected by key factors predicted by classic metapopulation models to affect dispersal-spatial and temporal variation in population size, inter-island distance, local demography and habitat characteristics. This metapopulation can be divided into two major habitat types: (a) islands closer to the mainland where sparrows breed in colonies on farms, and (b) islands without farms, situated farther away from the mainland where sparrows are exposed to harsher environmental conditions. Dispersal was spatially structured within the metapopulation; there was proportionally and numerically less emigration and immigration involving farm islands, as compared to non-farm islands. Furthermore, emigration and immigration occurred mostly between nearby islands. Moreover, emigration in response to spatial differences in mean population size differed between the habitat types, but populations with large mean received more immigrants in both habitat types. The number of emigrants and immigrants was negatively related to long-term recruit production, which was not the case in non-farm islands. The proportion and number of emigrants was positively related to temporal increases in recruit production on farm islands, however not on non-farm islands. Our results demonstrate that spatial heterogeneity in environmental conditions influences how spatial variation in long-term mean population size, and temporal and spatial variation in recruit production, affects dispersal dynamics. The spatial structure of this metapopulation is therefore best described by a spatially explicit model in which the exchange of individuals within each habitat type is strongly affected by the degree of geographical isolation, population size and recruit production. However, these relationships differed between the two habitat types; non-farm islands showing similarities to a mainland-island model type of structure, whereas farm islands showed features more associated with source-sink or balanced dispersal models. Such differential dispersal dynamics between habitat types are expected to have important consequences for the ecological and evolutionary dynamics within this metapopulation.

Metapopulation Structure Predicts Population Dynamics in the Cakile maritima-Alternaria brassicicola Host-Pathogen Interaction

AbstractIn symbiotic interactions, spatiotemporal variation in the distribution or population dynamics of one species represents spatial and temporal heterogeneity of the landscape for the other. Such interdependent demographic dynamics result in situations where the relative importance of biotic and abiotic factors in determining ecological processes is complicated to decipher. Using a detailed survey of three metapopulations of the succulent plant Cakile maritima and the necrotrophic fungus Alternaria brassicicola located along the southeastern Australian coast, we developed a series of statistical analyses-namely, synchrony analysis, patch occupancy dynamics, and a spatially explicit metapopulation model-to understand how habitat quality, weather conditions, dispersal, and spatial structure determine metapopulation dynamics. Climatic conditions are important drivers, likely explaining the high synchrony among populations. Host availability, landscape features facilitating dispersal, and habitat conditions also impact the occurrence and spread of disease. Overall, we show that the collection of extensive data on host and pathogen population dynamics, in combination with spatially explicit epidemiological modeling, makes it possible to accurately predict disease dynamics-even when there is extreme variability in host population dynamics. Finally, we discuss the importance of genetic information for predicting demographic dynamics in this pathosystem.

Feedback between coevolution and epidemiology can help or hinder the maintenance of genetic variation in host-parasite models

Antagonistic coevolution has long been suggested to help maintain host genetic variation. Although ecological and epidemiological feedbacks are known to have important consequences on coevolutionary allele-frequency dynamics, their effects on the maintenance of genetic variation remains poorly understood. Here, we extend previous work on the maintenance of genetic variation in a classic matching alleles coevolutionary model by exploring the effects of ecological and epidemiological feedbacks, where both allele frequencies and population sizes are allowed to vary over time. We find that coevolution rarely maintains more host genetic variation than expected under neutral genetic drift alone. When and if coevolution maintains or depletes genetic variation relative to neutral drift is determined, predominantly, by two factors: the deterministic stability of the Red Queen allele-frequency cycles and the chance of allele fixation in the pathogen, as this results in directional selection and depletion of genetic variation in the host. Compared to purely coevolutionary models with constant host and pathogen population sizes, ecological and epidemiological feedbacks stabilize Red Queen cycles deterministically, but population fluctuations in the pathogen increase the rate of allele fixation in the pathogen, especially in epidemiological models. Our results illustrate the importance of considering the ecological and epidemiological context in which coevolution occurs when examining the impact of Red Queen cycles on genetic variation.

The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey

Coevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.

Coevolution fails to maintain genetic variation in a host-parasite model with constant finite population size

Coevolutionary negative frequency-dependent selection has been hypothesized to maintain genetic variation in host and parasites. Despite the extensive literature pertaining to host-parasite coevolution, the temporal dynamics of genetic variation have not been examined in a matching-alleles model (MAM) with a finite population size relative to the expectation under neutral genetic drift alone. The dynamics of the MA coevolution in an infinite population, in fact, suggests that genetic variation in these coevolving populations behaves neutrally. By comparing host heterozygosity to the expectation in a single-species model of neutral genetic drift we find that while this is also largely true in finite populations two additional phenomena arise. First, reciprocal natural selection acting on stochastic perturbations in host and pathogen allele frequencies results in a slight increase or decrease in genetic variation depending on the parameter conditions. Second, following the fixation of an allele in the parasite, selection in the MAM becomes directional, which then rapidly erodes genetic variation in the host. Hence, rather than maintain it, we find that, on average, matching-alleles coevolution depletes genetic variation.

Strong Gene Flow Undermines Local Adaptations in a Host Parasite System

Simple Summary

The co-evolution of hosts and parasites depends on their ability to adapt to each other’s defense and counter-defense mechanisms. The strength of selection on those mechanisms may vary among populations, resulting in a geographical mosaic of co-evolution. The boreo-montane paper wasp Polistes biglumis and its parasite Polistes atrimandibularis exemplify this type of co-evolutionary system. Here, we used genetic markers to examine the genetic population structures of these wasps in the western Alps. We found that both host and parasite populations displayed similar levels of genetic variation. In the host species, populations located near to each other were genetically similar; in both the host and the parasite species populations farther apart were significantly different. Thus, apparent dispersal barriers (i.e., high mountains) did not seem to restrict gene flow across populations as expected. Furthermore, there were no major differences in gene flow between the two species, perhaps because P. atrimandibularis parasitizes both alpine and lowland host species and annually migrates between alpine and lowland populations. The presence of strong gene flow in a system where local populations experience variable levels of selection pressure challenges the classical hypothesis that restricted gene flow is required for local adaptations to evolve.

Abstract

The co-evolutionary pathways followed by hosts and parasites strongly depend on the adaptive potential of antagonists and its underlying genetic architecture. Geographically structured populations of interacting species often experience local differences in the strength of reciprocal selection pressures, which can result in a geographic mosaic of co-evolution. One example of such a system is the boreo-montane social wasp Polistes biglumis and its social parasite Polistes atrimandibularis, which have evolved local defense and counter-defense mechanisms to match their antagonist. In this work, we study spatial genetic structure of P. biglumis and P. atrimandibularis populations at local and regional scales in the Alps, by using nuclear markers (DNA microsatellites, AFLP) and mitochondrial sequences. Both the host and the parasite populations harbored similar amounts of genetic variation. Host populations were not genetically structured at the local scale, but geographic regions were significantly differentiated from each other in both the host and the parasite in all markers. The net dispersal inferred from genetic differentiation was similar in the host and the parasite, which may be due to the annual migration pattern of the parasites between alpine and lowland populations. Thus, the apparent dispersal barriers (i.e., high mountains) do not restrict gene flow as expected and there are no important gene flow differences between the species, which contradict the hypothesis that restricted gene flow is required for local adaptations to evolve.

Quantification of Outcrossing Events in Haploid Fungi Using Microsatellite Markers

Species in genera of the fungal family Ceratocystidaceae are known to have different mating strategies, including heterothallism and homothallism. Of these, species of Ceratocystis, typified by the pathogen Ceratocystis fimbriata all undergo unidirectional mating-type switching. This implies that the pathogens possess the ability to self, but also to undergo sexual outcrossing between isolates of different mating types. In this study, we extended the recently developed microsatellite-based technique to determine the extent to which outcrossing occurs in ascospore masses of haploid fungi to two field collections of Ceratocystis albifundus. In this way, the role of reproductive strategies in shaping population structure and diversity could be better understood. Results showed that a high frequency of outcrossing occurs in isolates of the pathogen from both non-native and native areas. This explains the high level of genetic diversity previously observed in this population despite the fact that this pathogen has the ability to self.

Parasite-host specificity: A cross-infection study of the parasite Ophryocystis elektroscirrha

Many parasites are constrained to only one or a few hosts, showing host specificity. It remains unclear why some parasites are specialists and other parasites are generalists. The parasite Ophryocystis elektroscirrha (OE) is a neogregarine protozoan thought to be restricted to monarch butterflies, Danaus plexippus (Nymphaliae) and D. gilippus. Recently, we found OE-like spores in other Lepidoptera, specifically in three noctuid moths: Helicoverpa armigera, H. assulta and H. punctigera, as well as another nymphalid, Parthenos sylvia. To our knowledge, this is the first report of OE-like parasite infections in species other than the genus Danaus. In sequencing 558 bp of 18S rRNA, we found the genetic similarity between OE from D. plexippus and OE-like parasite from the moths H. armigera and H. punctigera to be 95.2%. When we conducted cross-species infection experiments, we could not infect the moths with OE from D. plexippus, but OE-like parasite from H. armigera did infect D. plexippus and a closely related moth species Heliothis virescens. Interestingly, we did not find the OE-like parasite in the H. armigera population from Spain. Inter-population infection experiments with H. armigera demonstrated a higher sensitivity to OE-like infection in the population from Spain compared to the populations from Australia and China. These results suggest geographic variation in OE-like susceptibility and coevolution between parasite and host. Our findings give important new insights into the prevalence and host specificity of OE and OE-like parasites, and provide opportunities to study parasite transmission over spatial and temporal scales.

Parasite-mediated selection in a natural metapopulation of Daphnia magna

Parasite-mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease-related traits. However, nonadaptive processes like migration and extinction-(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life-history and disease-related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite-mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life-history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST -F_ST -like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction-(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.

A look into the genetic diversity of Lecanosticta acicola in northern Europe

For northern Europe Lecanosticta acicola is an emerging pine needle pathogen. This study gives a first look into the population genetics of the pathogen in Estonia, the first population documented in that region. The main aim of this study was to investigate the genetic diversity and population structure of the pathogen in this new region for the fungus. For this purpose, 104 isolates from 2010 to 2017 were analysed with 11 microsatellite and mating type markers. The stand where the pathogen's jump from an exotic host to the native Scots pine was recorded was also involved in this analysis. The analysis revealed low genetic diversity and a high number of clones that indicated L. acicola is an invasive species in northern Europe. Results suggest that several separate introductions have taken place and anthropogenic activity has apparently affected the spread of the pathogen. Clonal reproduction is dominating and although sexual reproduction is possible, it probably takes place infrequently.

Sex differences in dispersal syndrome are modulated by environment and evolution

Dispersal syndromes (i.e. suites of phenotypic correlates of dispersal) are potentially important determinants of local adaptation in populations. Species that exhibit sexual dimorphism in their life history or behaviour may exhibit sex-specific differences in their dispersal syndromes. Unfortunately, there is little empirical evidence of sex differences in dispersal syndromes and how they respond to environmental change or dispersal evolution. We investigated these issues using two same-generation studies and a long-term (greater than 70 generations) selection experiment on laboratory populations of Drosophila melanogaster There was a marked difference between the dispersal syndromes of males and females, the extent of which was modulated by nutrition availability. Moreover, dispersal evolution via spatial sorting reversed the direction of dispersal×sex interaction in one trait (desiccation resistance), while eliminating the sex difference in another trait (body size). Thus, we show that sex differences obtained through same-generation trait-associations ('ecological dispersal syndromes') are probably environment-dependent. Moreover, even under constant environments, they are not good predictors of the sex differences in 'evolutionary dispersal syndrome' (i.e. trait-associations shaped during dispersal evolution). Our findings have implications for local adaptation in the context of sex-biased dispersal and habitat-matching, as well as for the use of dispersal syndromes as a proxy of dispersal.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Role of host genetic diversity for susceptibility-to-infection in the evolution of virulence of a plant virus

Predicting viral emergence is difficult due to the stochastic nature of the underlying processes and the many factors that govern pathogen evolution. Environmental factors affecting the host, the pathogen and the interaction between both are key in emergence. In particular, infectious disease dynamics are affected by spatiotemporal heterogeneity in their environments. A broad knowledge of these factors will allow better estimating where and when viral emergence is more likely to occur. Here, we investigate how the population structure for susceptibility-to-infection genes of the plant Arabidopsis thaliana shapes the evolution of Turnip mosaic virus (TuMV). For doing so we have evolved TuMV lineages in two radically different host population structures: (1) a metapopulation subdivided into six demes (subpopulations); each one being composed of individuals from only one of six possible A. thaliana ecotypes and (2) a well-mixed population constituted by equal number of plants from the same six A. thaliana ecotypes. These two populations were evolved for twelve serial passages. At the end of the experimental evolution, we found faster adaptation of TuMV to each ecotype in the metapopulation than in the well-mixed heterogeneous host populations. However, viruses evolved in well-mixed populations were more pathogenic and infectious than viruses evolved in the metapopulation. Furthermore, the viruses evolved in the demes showed stronger signatures of local specialization than viruses evolved in the well-mixed populations. These results illustrate how the genetic diversity of hosts in an experimental ecosystem favors the evolution of virulence of a pathogen.

The return of the frogs: The importance of habitat refugia in maintaining diversity during a disease outbreak

Recent decades have seen the emergence and spread of numerous infectious diseases, often with severe negative consequences for wildlife populations. Nevertheless, many populations survive the initial outbreaks, and even undergo recoveries. Unfortunately, the long-term effects of these outbreaks on host population genetics are poorly understood; to increase this understanding, we examined the population genetics of two species of rainforest frogs (Litoria nannotis and Litoria serrata) that have largely recovered from a chytridiomycosis outbreak at two national parks in the Wet Tropics of northern Australia. At the wetter, northern park there was little evidence of decreased genetic diversity in either species, and all of the sampled sites had high minor allele frequencies (mean MAF = 0.230-0.235), high heterozygosity (0.318-0.325), and few monomorphic markers (1.4%-4.0%); however, some recovered L. nannotis populations had low N_e values (59.3-683.8) compared to populations that did not decline during the outbreak (1,537.4-1,756.5). At the drier, southern park, both species exhibited lower diversity (mean MAF = 0.084-0.180; heterozygosity = 0.126-0.257; monomorphic markers = 3.7%-43.5%; N_e  = 18.4-676.1). The diversity patterns in this park matched habitat patterns, with both species having higher diversity levels and fewer closely related individuals at sites with higher quality habitat. These patterns were more pronounced for L. nannotis, which has lower dispersal rates than L. serrata. These results suggest that refugia with high quality habitat are important for retaining genetic diversity during disease outbreaks, and that gene flow following disease outbreaks is important for re-establishing diversity in populations where it was reduced.

Can a Terrestrial Ectoparasite Disperse with Its Marine Host?

One of the most extreme examples of parasite adaptation comes from terrestrial ectoparasites exploiting marine hosts. Despite the ubiquity of such ectoparasitism and its ecological and evolutionary importance, investigations of the responses of ectoparasites to conditions encountered on their hosts are rare. In the case of penguins and their ticks, current understanding suggests that ticks freely parasitize their hosts on land but are incapable of surviving extended oceanic journeys. We examined this conjecture by assessing the physiological capacity of little penguin ticks to endure at-sea foraging and dispersal events of their hosts. Survival in penguins ticks was not significantly compromised by exposure to depths commonly associated with host dives (40 and 60 m), repeated seawater exposure relevant to the most common (30 s) and longest (120 s) recorded host dives, or extended (48 h) exposure to seawater. Mean (±SD) closed-phase durations in adult and nymphal ticks exhibiting discontinuous gas exchange ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>339</mml:mn><mml:mo>±</mml:mo><mml:mn>237</mml:mn></mml:mrow></mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>240</mml:mn><mml:mo>±</mml:mo><mml:mn>295</mml:mn></mml:mrow></mml:math> s, respectively) exceeded that of the maximum recorded host dive duration (120 s). Normoxic-anoxic-normoxic respirometry also confirmed spiracle closure. Mean metabolic rates ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>0.354</mml:mn><mml:mo>±</mml:mo><mml:mn>0.220</mml:mn></mml:mrow></mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>4.853</mml:mn><mml:mo>±</mml:mo><mml:mn>4.930</mml:mn></mml:mrow></mml:math> μL/h at 25°C for unfed and fed adult females, respectively) were significantly influenced by temperature; optimal and LT₅₀ temperatures for adult ticks and fed nymphal ticks were typically higher than swimming penguin body temperatures. These findings suggest that marine host dispersal is unlikely to present an insurmountable barrier to long-distance tick dispersal. Such dispersal has important implications for evolutionary theory, conservation, and epidemiology.

Foliar fungal endophyte communities are structured by environment but not host ecotype in Panicum virgatum (switchgrass)

Experimental tests of community assembly mechanisms for host-associated microbiomes in nature are lacking. Asymptomatic foliar fungal endophytes are a major component of the plant microbiome and are increasingly recognized for their impacts on plant performance, including pathogen defense, hormonal manipulation, and drought tolerance. However, it remains unclear whether fungal endophytes preferentially colonize certain host ecotypes or genotypes, reflecting some degree of biotic adaptation in the symbioses, or whether colonization is simply a function of spore type and abundance within the local environment. Whether host ecotype, local environment, or some combination of both controls the pattern of microbiome formation across hosts represents a new dimension to the age-old debate of nature versus nurture. Here, we used a reciprocal transplant design to explore the extent of host specificity and biotic adaptation in the plant microbiome, as evidenced by differential colonization of host genetic types by endophytes. Specifically, replicate plants from three locally-adapted ecotypes of the native grass Panicum virgatum (switchgrass) were transplanted at three geographically distinct field sites (one home and two away) in the Midwestern US. At the end of the growing season, plant leaves were harvested and the fungal microbiome characterized using culture-dependent sequencing techniques. Our results demonstrated that fungal endophyte community structure was determined by local environment (i.e., site), but not by host ecotype. Fungal richness and diversity also strongly differed by site, with lower fungal diversity at a riparian field site, whereas host ecotype had no effect. By contrast, there were significant differences in plant phenotypes across all ecotypes and sites, indicating ecotypic differentiation of host phenotype. Overall, our results indicate that environmental factors are the primary drivers of community structure in the switchgrass fungal microbiome.

Adaptation of a Fungal Pathogen to Host Quantitative Resistance

Impact of host quantitative resistance on pathogen evolution is still poorly documented. In our study, we characterized the adaptation of the pathogenic fungus Colletotrichum gloeosporioides, to the quantitative resistance of its host, the water yam (Dioscorea alata). Genetic and pathogenic diversities of C. gloeosporioides populations were specified at the field scale. We used nuclear markers to describe fungal population structuring within and between six fields of three cultivars differently susceptible to the fungus. Strain aggressiveness was then quantified in the laboratory through cross-inoculation tests. The high level of genetic diversity and significant linkage disequilibrium revealed a significant influence of clonal reproduction in the C. gloeosporioides evolution. The recorded fungal migration between fields was weak (evidence for a dispersion mode via tubers rather than splashing dispersal), which provides the first molecular evidence for limited C. gloeosporioides migration via yam tuber exchanges. C. gloeosporioides's populations are adapted to their host resistance. The aggressiveness of the fungal clones seems to have evolved toward an accumulation of components specific to each host cultivar. Despite the remaining marks of adaptation to the former widely cultivated host, adaptation to current cultivars was clearly depicted.

Immune-Related Transcriptional Responses to Parasitic Infection in a Naturally Inbred Fish: Roles of Genotype and Individual Variation

Parasites are strong drivers of evolutionary change and the genetic variation of both host and parasite populations can co-evolve as a function of parasite virulence and host resistance. The role of transcriptome variation in specific interactions between host and parasite genotypes has been less studied and can be confounded by differences in genetic variation. We employed two naturally inbred lines of a self-fertilizing fish to estimate the role of host genotype in the transcriptome response to parasite infection using RNA-seq. In addition, we targeted several differentially expressed immune-related genes to further investigate the relative role of individual variation in the immune response using RT-qPCR, taking advantage of the genomic uniformity of the self-fertilizing lines. We found significant differences in gene expression between lines in response to infection both in the transcriptome and in individual gene RT-qPCR analyses. Individual RT-qPCR analyses of gene expression identified significant variance differences between lines for six genes but only for three genes between infected and control fish. Our results indicate that although the genetic background plays an important role in the transcriptome response to parasites, it cannot fully explain individual differences within genetically homogeneous lines, which can be important for determining the response to parasites.

Comparison of population genetic structures of the plant Silene stellata and its obligate pollinating seed predator moth Hadena ectypa

BACKGROUND AND AIMS

Population genetic structures and patterns of gene flow of interacting species provide important insights into the spatial scale of their interactions and the potential for local co-adaptation. We analysed the genetic structures of the plant Silene stellata and the nocturnal moth Hadena ectypa. Hadena ectypa acts as one of the important pollinators of S. stellata as well as being an obligate seed parasite on the plant. Although H. ectypa provides a substantial pollination service to S. stellata, this system is largely considered parasitic due to the severe seed predation by the Hadena larvae. Previous research on this system has found variable interaction outcomes across space, indicating the potential for a geographical selection mosaic.

METHODS

Using 11 microsatellite markers for S. stellata and nine markers for H. ectypa, we analysed the population genetic structure and the patterns and intensity of gene flow within and among three local populations in the Appalachians.

KEY RESULTS

We found no spatial genetic structure in the moth populations, while significant differentiation was detected among the local plant populations. Additionally, we observed that gene flow rates among H. ectypa populations were more uniform and that the mean gene flow rate in H. ectypa was twice as large as that in S. stellata.

CONCLUSIONS

Our results suggest that although the moths move frequently among populations, long-distance pollen carryover only happens occasionally. The difference in gene flow rates between S. stellata and H. ectypa could prevent strict local co-adaptation. Furthermore, higher gene flow rates in H. ectypa could also increase resistance of the local S. stellata populations to the parasitic effect of H. ectypa and therefore help to stabilize the Silene-Hadena interaction dynamics.

Evolution of dispersal syndrome and its corresponding metabolomic changes

Dispersal is one of the strategies for organisms to deal with climate change and habitat degradation. Therefore, investigating the effects of dispersal evolution on natural populations is of considerable interest to ecologists and conservation biologists. Although it is known that dispersal itself can evolve due to selection, the behavioral, life-history and metabolic consequences of dispersal evolution are not well understood. Here, we explore these issues by subjecting four outbred laboratory populations of Drosophila melanogaster to selection for increased dispersal. The dispersal-selected populations had similar values of body size, fecundity, and longevity as the nonselected lines (controls), but evolved significantly greater locomotor activity, exploratory tendency, and aggression. Untargeted metabolomic fingerprinting through NMR spectroscopy suggested that the selected flies evolved elevated cellular respiration characterized by greater amounts of glucose, AMP, and NAD. Concurrent evolution of higher level of Octopamine and other neurotransmitters indicate a possible mechanism for the behavioral changes in the selected lines. We discuss the generalizability of our findings in the context of observations from natural populations. To the best of our knowledge, this is the first report of the evolution of metabolome due to selection for dispersal and its connection to dispersal syndrome evolution.

Positive selection and climatic effects on MHC class II gene diversity in hares (Lepus capensis) from a steep ecological gradient

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

Gene flow and adaptive potential in a generalist ectoparasite

Background

In host-parasite systems, relative dispersal rates condition genetic novelty within populations and thus their adaptive potential. Knowledge of host and parasite dispersal rates can therefore help us to understand current interaction patterns in wild populations and why these patterns shift over time and space. For generalist parasites however, estimates of dispersal rates depend on both host range and the considered spatial scale. Here, we assess the relative contribution of these factors by studying the population genetic structure of a common avian ectoparasite, the hen flea Ceratophyllus gallinae, exploiting two hosts that are sympatric in our study population, the great tit Parus major and the collared flycatcher Ficedula albicollis. Previous experimental studies have indicated that the hen flea is both locally maladapted to great tit populations and composed of subpopulations specialized on the two host species, suggesting limited parasite dispersal in space and among hosts, and a potential interaction between these two structuring factors.

Results

C. gallinae fleas were sampled from old nests of the two passerine species in three replicate wood patches and were genotyped at microsatellite markers to assess population genetic structure at different scales (among individuals within a nest, among nests and between host species within a patch and among patches). As expected, significant structure was found at all spatial scales and between host species, supporting the hypothesis of limited dispersal in this parasite. Clustering analyses and estimates of relatedness further suggested that inbreeding regularly occurs within nests. Patterns of isolation by distance within wood patches indicated that flea dispersal likely occurs in a stepwise manner among neighboring nests. From these data, we estimated that gene flow in the hen flea is approximately half that previously described for its great tit hosts.

Conclusion

Our results fall in line with predictions based on observed patterns of adaptation in this host-parasite system, suggesting that parasite dispersal is limited and impacts its adaptive potential with respect to its hosts. More generally, this study sheds light on the complex interaction between parasite gene flow, local adaptation and host specialization within a single host-parasite system.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1205-2) contains supplementary material, which is available to authorized users.

Host-parasite coevolution can promote the evolution of seed banking as a bet-hedging strategy

Seed (egg) banking is a common bet-hedging strategy maximizing the fitness of organisms facing environmental unpredictability by the delayed emergence of offspring. Yet, this condition often requires fast and drastic stochastic shifts between good and bad years. We hypothesize that the host seed banking strategy can evolve in response to coevolution with parasites because the coevolutionary cycles promote a gradually changing environment over longer times than seed persistence. We study the evolution of host germination fraction as a quantitative trait using both pairwise competition and multiple mutant competition methods, while the germination locus can be genetically linked or unlinked with the host locus under coevolution. In a gene-for-gene model of coevolution, hosts evolve a seed bank strategy under unstable coevolutionary cycles promoted by moderate to high costs of resistance or strong disease severity. Moreover, when assuming genetic linkage between coevolving and germination loci, the resistant genotype always evolves seed banking in contrast to susceptible hosts. Under a matching-allele interaction, both hosts' genotypes exhibit the same seed banking strategy irrespective of the genetic linkage between loci. We suggest host-parasite coevolution as an additional hypothesis for the evolution of seed banking as a temporal bet-hedging strategy.

Effect of spatial connectivity on host resistance in a highly fragmented natural pathosystem

Both theory and experimental evolution studies predict migration to influence the outcome of antagonistic coevolution between hosts and their parasites, with higher migration rates leading to increased diversity and evolutionary potential. Migration rates are expected to vary in spatially structured natural pathosystems, yet how spatial structure generates variation in coevolutionary trajectories across populations occupying the same landscape has not been tested. Here, we studied the effect of spatial connectivity on host evolutionary potential in a natural pathosystem characterized by a stable Plantago lanceolata host network and a highly dynamic Podosphaera plantaginis parasite metapopulation. We designed a large inoculation experiment to test resistance of five isolated and five well‐connected host populations against sympatric and allopatric pathogen strains, over 4 years. Contrary to our expectations, we did not find consistently higher resistance against sympatric pathogen strains in the well‐connected populations. Instead, host local adaptation varied considerably among populations and through time with greater fluctuations observed in the well‐connected populations. Jointly, our results suggest that in populations where pathogens have successfully established, they have the upper hand in the coevolutionary arms race, but hosts may be better able to respond to pathogen‐imposed selection in the well‐connected than in the isolated populations. Hence, the ongoing and extensive fragmentation of natural habitats may increase vulnerability to diseases.

Genome-Wide Detection of Genes Under Positive Selection in Worldwide Populations of the Barley Scald Pathogen

Coevolution between hosts and pathogens generates strong selection pressures to maintain resistance and infectivity, respectively. Genomes of plant pathogens often encode major effect loci for the ability to successfully infect specific host genotypes. Hence, spatial heterogeneity in host genotypes coupled with abiotic factors could lead to locally adapted pathogen populations. However, the genetic basis of local adaptation is poorly understood. Rhynchosporium commune, the pathogen causing barley scald disease, interacts at least partially in a gene-for-gene manner with its host. We analyzed global field populations of 125 R. commune isolates to identify candidate genes for local adaptation. Whole genome sequencing data showed that the pathogen is subdivided into three genetic clusters associated with distinct geographic and climatic regions. Using haplotype-based selection scans applied independently to each genetic cluster, we found strong evidence for selective sweeps throughout the genome. Comparisons of loci under selection among clusters revealed little overlap, suggesting that ecological differences associated with each cluster led to variable selection regimes. The strongest signals of selection were found predominantly in the two clusters composed of isolates from Central Europe and Ethiopia. The strongest selective sweep regions encoded protein functions related to biotic and abiotic stress responses. Selective sweep regions were enriched in genes encoding functions in cellular localization, protein transport activity, and DNA damage responses. In contrast to the prevailing view that a small number of gene-for-gene interactions govern plant pathogen evolution, our analyses suggest that the evolutionary trajectory is largely determined by spatially heterogeneous biotic and abiotic selection pressures.

Non-Mendelian segregation influences the infection biology and genetic structure of the African tree pathogen Ceratocystis albifundus

The African fungal tree pathogen, Ceratocystis albifundus, undergoes uni-directional mating type switching, giving rise to either self-fertile or self-sterile progeny. Self-sterile isolates lack the MAT1-2-1 gene and have reduced fitness such as slower growth and reduced pathogenicity, relative to self-fertile isolates. While it has been hypothesized that there is a 1:1 ratio of self-fertile to self-sterile ascospore progeny in relatives of C. albifundus, some studies have reported a significant bias in this ratio. This could be due to the fact that either fewer self-sterile ascospores are produced or that self-sterile ascospores have low viability. We quantified the percentage of self-sterile and self-fertile ascospores from ascospore masses in C. albifundus using real-time PCR. Primers were designed to distinguish between spores that contained the MAT1-2-1 gene and those where this gene had been deleted. A significant bias towards the self-fertile mating type was observed in all single ascospore masses taken from sexual structures produced in haploid-selfed cultures. The same result was observed from a disease outbreak situation in an intensively managed field of cultivated native trees, and this was coupled with very low population diversity in the pathogen. This was in contrast to the results obtained from ascospore masses taken from the crosses performed under laboratory conditions or ascomata on native trees in a non-disease situation, where either self-fertile or self-sterile ascospores were dominant. The results suggest that reproductive strategies play a significant role in the infection biology and genetic structure of C. albifundus populations.

Evidence for parasite-mediated selection during short-lasting toxic algal blooms

Parasites play a role in the control of transient algal blooms, but it is not known whether parasite-mediated selection results in coevolution of the host and the parasites over this short time span. We investigated the presence of coevolution between the toxic dinoflagellate Alexandrium minutum and two naturally occurring endoparasites during blooms lasting a month in two river estuaries, using cross-inoculation experiments across time and space. Higher parasite abundance was associated with a large daily reduction in relative A. minutum abundances, demonstrating strong parasite-mediated selection. There was genetic variability in infectivity in both parasite species, and in resistance in the host. We found no evidence for coevolution in one estuary; however, in the other estuary, we found high genetic diversity in the two parasite species, fluctuations in infectivity and suggestion that the two parasites are well adapted to their host, as in ‘Red Queen’ dynamics. Thus, coevolution is possible over the short time span of a bloom, but geographically variable, and may feedback on community dynamics.

Dynamic, spatial models of parasite transmission in wildlife: Their structure, applications and remaining challenges

Individual differences in contact rate can arise from host, group and landscape heterogeneity and can result in different patterns of spatial spread for diseases in wildlife populations with concomitant implications for disease control in wildlife of conservation concern, livestock and humans. While dynamic disease models can provide a better understanding of the drivers of spatial spread, the effects of landscape heterogeneity have only been modelled in a few well-studied wildlife systems such as rabies and bovine tuberculosis. Such spatial models tend to be either purely theoretical with intrinsic limiting assumptions or individual-based models that are often highly species- and system-specific, limiting the breadth of their utility. Our goal was to review studies that have utilized dynamic, spatial models to answer questions about pathogen transmission in wildlife and identify key gaps in the literature. We begin by providing an overview of the main types of dynamic, spatial models (e.g., metapopulation, network, lattice, cellular automata, individual-based and continuous-space) and their relation to each other. We investigate different types of ecological questions that these models have been used to explore: pathogen invasion dynamics and range expansion, spatial heterogeneity and pathogen persistence, the implications of management and intervention strategies and the role of evolution in host-pathogen dynamics. We reviewed 168 studies that consider pathogen transmission in free-ranging wildlife and classify them by the model type employed, the focal host-pathogen system, and their overall research themes and motivation. We observed a significant focus on mammalian hosts, a few well-studied or purely theoretical pathogen systems, and a lack of studies occurring at the wildlife-public health or wildlife-livestock interfaces. Finally, we discuss challenges and future directions in the context of unprecedented human-mediated environmental change. Spatial models may provide new insights into understanding, for example, how global warming and habitat disturbance contribute to disease maintenance and emergence. Moving forward, better integration of dynamic, spatial disease models with approaches from movement ecology, landscape genetics/genomics and ecoimmunology may provide new avenues for investigation and aid in the control of zoonotic and emerging infectious diseases.

Galba truncatula and Fasciola hepatica: Genetic costructures and interactions with intermediate host dispersal

Antagonistic interactions between hosts and parasites are key structuring forces in natural populations. Demographic factors like extinction, migration and the effective population size shape host-parasite metapopulational dynamics. Therefore, to understand the evolution of host-parasite systems it is necessary to study the distribution of the genetic variation of both entities simultaneously. In this paper, we investigate the population genetics co-structure of parasites and hosts within a metapopulation of the liver fluke, Fasciola hepatica, and two of its intermediate hosts, the main intermediate host in Europe, Galba truncatula, and a new intermediate host, Omphiscola glabra, in Central France. Our results reveal an absence of specificity of flukes as regard to the two alternative hosts though O. glabra shows higher prevalence of F. hepatica. Host and parasites displayed contrasting population genetics structure with very small, highly inbred (selfing) and strongly isolated G. truncatula populations and much bigger, panmictic and more dispersive F. hepatica. This could indicate a local adaptation of the parasite and a local maladaptation of the host. We also unveil a parasite-mediated biased population genetics structure suggesting that infected G. truncatula disperse more; have higher dispersal survival than uninfected snails or, more likely, that immigrant snails are infected more often than local snails (local parasites are less adapted to local hosts). Finally, an absence, or at least an ambiguous signature of isolation by distance was observed in both host and parasite population. A very weak migration rate for G. truncatula provides a reasonable explanation for this ambiguous result. Alternatively, smaller sample sizes combined with modest migration rates might explain the difficulties to unveil the signal in F. hepatica.

Rethinking Conventional Wisdom: Are Locally Adapted Parasites Ahead in the Coevolutionary Race?

The metaphors of the Red Queen and the arms race have inspired a large amount of research aimed at understanding the process of antagonistic coevolution between hosts and parasites. One approach that has been heavily used is to estimate the strength of parasite local adaptation using a reciprocal cross infection or transplant study. These studies frequently conclude that the locally adapted species is ahead in the coevolutionary race. Here, I use mathematical models to decompose parasite infectivity into components attributable to local versus global adaptation and to develop a robust index of which species is ahead in the coevolutionary race, which I term coevolutionary advantage. Computer simulations of coevolving host-parasite interactions demonstrate that because the magnitudes of local and global adaptation are largely independent, the link between the sign of local adaptation and coevolutionary advantage is tenuous. A consequence of the weak coupling between local adaptation and coevolutionary advantage is that the bulk of existing empirical studies do not sample enough populations for any reliable conclusions to be drawn. Together, these results suggest that the long-standing conventional wisdom holding that locally adapted parasites are ahead in the coevolutionary race should be reconsidered.

Impact of life stage-dependent dispersal on the colonization dynamics of host patches by ticks and tick-borne infectious agents

BACKGROUND

When colonization and gene flow depend on host-mediated dispersal, a key factor affecting vector dispersal potential is the time spent on the host for the blood meal, a characteristic that can vary strongly among life history stages. Using a 2-patch vector-pathogen population model and seabird ticks as biological examples, we explore how vector colonization rates and the spread of infectious agents may be shaped by life stage-dependent dispersal. We contrast hard (Ixodidae) and soft (Argasidae) tick systems, which differ strongly in blood- feeding traits.

RESULTS

We find that vector life history characteristics (i.e. length of blood meal) and demographic constraints (Allee effects) condition the colonization potential of ticks; hard ticks, which take a single, long blood meal per life stage, should have much higher colonization rates than soft ticks, which take repeated short meals. Moreover, this dispersal potential has direct consequences for the spread of vector-borne infectious agents, in particular when transmission is transovarial.

CONCLUSIONS

These results have clear implications for predicting the dynamics of vector and disease spread in the context of large-scale environmental change. The findings highlight the need to include life-stage dispersal in models that aim to predict species and disease distributions, and provide testable predictions related to the population genetic structure of vectors and pathogens along expansion fronts.