Spatially structured eco-evolutionary dynamics in a host-pathogen interaction render isolated populations vulnerable to disease –

Spatially structured eco-evolutionary dynamics in a host-pathogen interaction render isolated populations vulnerable to disease –

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Nature Communications volume 13, Article number: 6018 (2022)
14 Altmetric
Metrics details
While the negative effects that pathogens have on their hosts are well-documented in humans and agricultural systems, direct evidence of pathogen-driven impacts in wild host populations is scarce and mixed. Here, to determine how the strength of pathogen-imposed selection depends on spatial structure, we analyze growth rates across approximately 4000 host populations of a perennial plant through time coupled with data on pathogen presence-absence. We find that infection decreases growth more in the isolated than well-connected host populations. Our inoculation study reveals isolated populations to be highly susceptible to disease while connected host populations support the highest levels of resistance diversity, regardless of their disease history. A spatial eco-evolutionary model predicts that non-linearity in the costs to resistance may be critical in determining this pattern. Overall, evolutionary feedbacks define the ecological impacts of disease in spatially structured systems with host gene flow being more important than disease history in determining the outcome.
According to coevolutionary theory, hosts may evolve resistance under pathogen-imposed negative frequency-dependent selection (NFDS), whereby rare host genotypes have an advantage over the common ones1,2. The underlying assumptions of coevolutionary theory are the strong negative fitness effect of infection, with disease-free individuals outperforming infected ones3, and costs of resistance that are central to maintenance of polymorphism within populations4. While consistent negative effects of pathogens on their host populations are well documented in humans and agricultural systems5,6, direct evidence of pathogen-driven ecological and evolutionary change in the wild is scarce and mixed3,7,8,9,10,11. The theoretical expectation is that the selective importance of disease is directly correlated with the frequency and severity of epidemics12. However, our ability to quantify the strength of pathogen-imposed selection in natural populations is limited by few available systematic spatio-temporal data on pathogen occurrence across a sufficient number of host populations.
Spatial structure and heterogeneity supported by natural host populations is in stark contrast to human-managed systems that are typically highly conductive to disease transmission due to large population sizes, high densities and low genetic variability13. Not surprisingly, studies focusing on wild pathosystems have revealed highly variable disease prevalence levels. Moreover, local pathogen populations are typically ephemeral, persisting regionally as metapopulations through extinction and colonization events of local host populations14,15,16,17. Even when infection takes place, the fitness consequences—and the coevolutionary outcomes18—may vary depending on the genetic composition of the host and pathogen populations and their environment, either directly or via GenotypeHOST × GenotypePATHOGEN × Environment—interactions19,20. Moreover, hosts in wild populations may suffer increased mortality or reduced reproduction irrespective of their infection status due to other factors such as extreme weather21. Hence, remarkably little is understood of how pathogens impact the fitness of their host populations in the wild.
There is increasing evidence that host-pathogen dynamics, both epidemiological and evolutionary, may be shaped by the spatial structure of the interaction13,22,23,24. Encounter rates between hosts and their pathogens are expected to be heavily influenced by connectivity to other populations, and the key metapopulation processes—gene flow, extinction, and colonization dynamics—are expected to contribute to the genetic structure of both the colonization dynamics, and the arrival of novel genetic variation into local populations13. As long as rates of migration are low enough to not homogenize local populations, increasing immigration is expected to increase the diversity and evolutionary potential of both host and pathogen populations25. While measuring migration rates in natural populations is difficult26, population connectivity, measured as the Euclidian distances separating populations and calibrated by the species dispersal capacity, provides a powerful proxy for migration rates27. Consequently, spatially structured eco-evolutionary feedback dynamics may emerge, with diversity accumulating in the well-connected populations. In line with this, there is evidence of spatial structure strongly influencing how resistance is distributed, with higher resistance observed in host populations that experience higher rates of gene flow16,28,29. To date, it has not been established what the relative roles of gene flow vs. pathogen-imposed selection are—and how they may vary in space—in generating spatially variable patterns of resistance that have been empirically observed16,28,29.
Here, we combine a spatial analysis of wild host-pathogen populations with an inoculation experiment, and a simulation model to understand how the ecological and evolutionary impacts of disease on host resistance vary in spatially structured populations. Specifically, we ask: (1) Is there evidence of pathogen-imposed selection on its host populations across a large, naturally fragmented host-pathogen metapopulation; (2) Does host population resistance structure, measured through an inoculation assay, reflect variable selection pressure indicated by the spatial analysis; and (3) Using a coevolutionary metapopulation model we explore how gene flow, selection, and costs of resistance contribute to the spatial structure of resistance detected with our empirical approach.
Our analysis is focused on annually recorded population size data (measured visually as coverage; m2) from some ~4000 locations of host plant Plantago lanceolata, and the presence-absence dynamics of its obligate fungal pathogen, Podosphaera plantaginis, in this host population network in the Åland islands, South-Western Finland. Plantago lanceolata is a perennial that produces wind-dispersed pollen, while seeds typically drop close to the mother plant. During the epidemic season, P. plantaginis disperses via clonally produced conidial spores that typically land within close proximity of the infected source plant16. The visually conspicuous symptoms caused by P. plantaginis enable accurate tracking of infection in the wild. Long-term epidemiological data have demonstrated this pathogen to occur as a highly dynamic metapopulation with frequent extinctions and (re)colonizations of local populations, typically persisting in any given host population only for a few years16. The host population spatial structure is a critical determinant of pathogen extinction-colonization dynamics: large host populations are more likely to become colonized and to sustain inifection16. In contrast to predictions of the metapopulation theory27, host population connectivity has a negative impact on pathogen colonization and persistence, suggesting these populations to vary in their suitability for the pathogen16. The host population network does not occur as a metapopulation30, but is characterized by strong fluctuations in population size31,32. These data allow us to study whether the extent of pathogen-imposed selection depends on host population connectivity (SH) and hence, evolutionary potential governed by gene flow, and whether resistance level and diversity vary among host populations depending on their degree of connectivity and disease history. Previous metapopulation models33,34 have demonstrated the existence of overall higher resistance in well-connected populations. To better understand the mechanisms that lead to the significant interaction between population connectivity, infection history and resistance in our inoculation study, we built a host-pathogen coevolutionary metapopulation model, where we examine how different trade-off relationships impact the outcome.
Jointly our results show that the strength of pathogen-imposed selection depends on host population spatial structure. Low disease resistance in the isolated populations renders them vulnerable to pathogen attack. In the well-connected host populations high rates of gene flow associate with high resistance diversity irrespective of population disase history.
We used Spatial Bayesian modelling (Integrated Nested Laplace Approximation; INLA35) to assess how changes in host population size are influenced by the pathogen. We analyzed the relative change in host population size (m2) (defined as population size (t) − population size (t−1))/population size (t−1)) between consecutive years utilizing data from 2001 to 2008, i.e., eight transitions in host population size in response to pathogen presence–absence status at t−1. To assess whether this depends on host population connectivity, we estimated the separate effects of pathogen presence/absence in the previous year for connectivity categories—high-, low, and intermediate—that were based on the 0.2 and 0.8 quantiles of the host-connectivity values (Supplementary Fig. 1). Earlier studies have demonstrated P. lanceolata populations in Åland to be sensitive to drought31,32 and hence, to reliably estimate the effect of the pathogen on host population growth rates, we included data on precipitation and field-estimated drought symptoms in our model. The model also controls for spatio-temporal autocorrelation characteristic of spatial ecological data, that may be due to unmeasured variables (e.g., habitat quality, prevailing wind-direction or other unmeasured biotic of abiotic variation), thereby providing a conservative estimate of the model parameters (Supplementary Table 1)35.
Infection by P. plantaginis had a negative effect on the growth of its host populations. Across all connectivity categories, the estimated mean effects of pathogen presence on host population growth were smaller than the effects with pathogen absence, suggesting an overall negative effect of the pathogen on host-population change (Fig. 1A and Supplementary Table 1). Furthermore, the estimated mean effects of the pathogen within the connectivity categories supports the interpretation that the relative effect of the pathogen on population growth is most negative in the isolated host populations (Fig. 1A and Supplementary Table 1). The posterior uncertainty in the effects of pathogen on the population growth (indicated by the confidence intervals in Fig. 1A) are due to the nature of observational data: pathogen infections were rare at the metapopulation level in studied years, thus there is considerably more pathogen absence observations in these data (See Supplementary Table 2). The temporal autocorrelation in growth in P. lanceolata populations between consecutive years was estimated to be negative (Supplementary Table 1), indicating that local populations exhibit oscillatory dynamics, such that growth in one year is typically followed by a decline in the next year and vice versa. As many of the populations are well-established, these fluctuations could result from populations oscillating around their carrying capacities, dictated by the space and resources available for their growth. The estimated median effects for rainfall in July and August suggest that host population changes are not strongly driven by these effects, although the August rainfall had a slight positive effect on population growth (posterior mean effect 0.03, confidence interval −0.06, 0.12, Fig. 1B, Supplementary Table 1). The proportion of plants expressing drought symptoms in the previous year was significantly associated with a decline in host population size (posterior mean effect −0.38, confidence interval −0.43, −0.33, Fig. 1C, Supplementary Table 1).
The estimated median effects for host population growth with 95% credibility intervals, shown with lines, of the fixed effects of the Bayesian INLA model: A The effect of pathogen presence and absence in the host populations in the three connectivity categories, B the effect of rainfall in July and August; and C the effect of detected drought symptoms in the host populations in the previous and current year. In the model the relative change in host population size (m2) is defined as population size (t) − population size (t1))/population size (t1)) between consecutive years utilizing data from 2001 to 2008 in response to pathogen presence-absence status at t1. Source data are provided as a Source Data file.
To examine whether the diversity and level of resistance vary among host populations depending on their degree of connectivity (SH) and disease history (measured as infection status in years 2001–2014), we performed an inoculation assay to characterize resistance phenotypes in plants sampled from 19 natural P. lanceolata populations. These populations occur in different locations of the host network, and were selected to represent both isolated and well-connected populations. Each plant was inoculated with four strains of P. plantaginis yielding resistance phenotypes values ranging between 0000 and 1111, with one depicting a resistant response and zero a susceptible response (the 16 possible resistance phenotype profiles are shown on the x-axis in Fig. 2A). Our inoculation study confirmed that host plants varied in their resistance against the tested powdery mildew strains (Table 1 and Fig. 2A). We were able to identify all 16 possible resistance phenotypes in the sample of 190 plants (Fig. 2A). In the connected populations, we found a greater diversity of different phenotypes, while isolated populations hosted fewer resistance phenotypes (Fig. 2A). Both the Shannon diversity index (Table 1 and Fig. 2B), and the average level of resistance (Table 1 and Fig. 2C), were higher in the well-connected than in the isolated host populations (Table 1 and Fig. 2B, C).
The phenotype composition of the19 study populations was defined by individual plant (n = 10 per population) responses to the four pathogen strains, resulting in a total of N = 190 observations visualized in panel (A). The matrix of detected resistance phenotypes in the inoculation study shows clustering of similar phenotypic profiles detected in populations in each of the four connectivity (SH)–infection history categories. The columns of the matrix correspond to resistance phenotypes, where the i’th element of the vector is 1, if resistance to pathogen strain I was detected, and zero otherwise. The rows of the matrix encode the observed frequencies of resistance phenotypes within the studied populations. The dendrogram visualizes the similarity structure between the populations, distance along the tree encoding for the degree of similarity between the populations. It is based on a hierarchical clustering (implemented with complete linkage method, aiming to find similar clusters), applied to Euclidean distances between the phenotype profiles within the populations. In panel B the diversity of detected 16 resistance phenotypes for the 19 study populations was characterized with a Shannon diversity index of host populations, shown separately for each connectivity (SH)-disease history category, and C the average resistance (%) of the same populations in each category. The centre lines of the boxplots B, C show the medians, box limits show the 25 and 75% quantiles, and the whiskers span to the data extremes. Purple colours depict isolated populations, and green colours well-connected populations. Source data are provided as a Source Data file.
However, while disease history had no direct effect on phenotypic diversity nor the level of resistance, we found a significant interaction between population connectivity and infection history for both Shannon’s diversity index and level of resistance (Table 1 and Fig. 2B, C). The highest diversity of phenotypes and highest resistance was measured in well-connected populations without any history of disease. In contrast, in isolated populations, we found greater diversity of resistance phenotypes and higher resistance in populations with a history of infection (Fig. 2B, C).
We modelled both the ecological and coevolutionary dynamics of host and pathogen metapopulations by constructing the network in two stages to account for relatively well and poorly connected demes (see “Methods”). We modelled the genetics of the system using a multilocus gene-for-gene framework36 with haploid host and pathogen genotypes characterized by (L) biallelic loci, where 0 and 1 represent the presence and absence, respectively, of resistance and infectivity alleles. Hosts and pathogen with more resistance or infectivity alleles are assumed to pay higher fitness costs, as defined in the methods. We ran 200 stochastic simulations using the tau-leap method37 for each of the parameter sets described in Supplementary Table 3 (example simulation dynamics are shown in Fig. 3D–F). On average, disease prevalence (D), resistance (R), and infectivity (I) were always higher in well-connected than in poorly connected populations regardless of metapopulation structure, transmissibility of the pathogen, or the nature of the trade-offs (Supplementary Table 3). However, the difference between well and poorly connected populations was generally greater when: (1) the metapopulation structure was assortative (i.e., well connected populations are more likely to be connected to other well connected populations than by chance) than random; (2) the pathogen was more transmissible; or (3) host resistance was associated with fitness costs that diminish as resistance increases (i.e., costs of resistance decelerate, ({c}_{H}^{2} ; < ; 0)) (Supplementary Table 4). Overall, we found that the pattern of the empirical results shown in Fig. 2C was most likely to occur when host resistance is associated with diminishing fitness costs and is more likely for transient (Fig. 3B) than long-term dynamics (Fig. 3C).
A Example snapshot of the simulation dynamics at t = 10,000 across a metapopulation with assortative connectivity, highlighting well (green) and poorly (purple) connected populations (unshaded populations are neither well nor poorly connected) that are currently infected (squares) and uninfected (circles). The size of each node corresponds to the mean resistance of the local population. Proportion of simulations which qualitatively match the empirical results as the shape of the host and pathogen cost functions are varied for transient (B) and long-term (C) dynamics: (strong decel. (decelerating): ({c}_{H}^{2},; {c}_{P}^{2}=-10); weak decel.: ({c}_{H}^{2},; {c}_{P}^{2}=-3); weak accel. (accelerating): ({c}_{H}^{2},; {c}_{P}^{2}=3); strong accel.: ({c}_{H}^{2},; {c}_{P}^{2}=10)). DF Example simulation results, showing mean (bold line) and standard deviations (shading) for disease prevalence, i.e., proportion of infected hosts (D), resistance (E), and infectivity, i.e., the average proportion of loci with infectivity alleles in the parasite population (F) in well (green) and poorly (purple) connected populations (({c}_{H}^{2}=-3), ({c}_{P}^{2}=10), (beta=0.01), with assortative network structure). Fixed parameters as defined in Supplementary Table 3.
Here, we show that the negative effect of pathogens on their wild host populations depends on spatial structure. This finding suggests that the strength of pathogen-imposed selection may vary across space in a predictable manner. Overall, finding a consistent negative effect of infection on host population growth is noteworthy given the myriad ecological factors that may hamper our ability to quantify costs of infection in wild populations38. The effect of infection on host population growth was the least negative in well-connected host populations, while isolated host populations were most vulnerable to infection, suggesting that they lack resistance diversity to effectively counter pathogen attack. Indeed, results of the inoculation study confirmed that both the diversity and the average level of resistance were higher in the well-connected than in the isolated host populations. When the interaction is characterized by strain-specific resistance such as in the interaction between P. lanceolata and P. plantaginis, resistance diversity will reduce the probability of establishment by an immigrant pathogen strain, and slow down the spread of established strains due to a mismatch between the specific avirulence alleles of pathogen and resistance alleles of host39. In agriculture, even slight additions of diversity to monocultures have been shown to reduce disease levels significantly40,41.
Theory predicts that pathogens maintain resistance polymorphism in their host populations42,43,44. As described above, our spatial statistical population model demonstrated that the isolated populations went through the strongest reductions in size—most likely through increased mortality of infected individuals31—which could lead to selection increasing in the frequency of resistant phenotypes locally32. Accordingly, in the isolated populations we measured higher resistance diversity in host populations with a history of infection than in host populations that had not been infected in the past. The effect of infection on host population growth rates in the well-connected populations was much weaker, and hence, may explain why we did not detect signs of past selection in these populations. Moreover, high rates of gene flow into the well-connected populations may swamp signatures of pathogen-imposed selection. The resulting differences in resistance among host populations is in line with previous studies that have measured higher resistance levels in well-connected host populations16,28,29. Jointly our results reveal that this pattern is generated by eco-evolutionary feedback resulting from spatial differences in how gene flow vs. selection drive host-pathogen dynamics in the wild. In well-connected populations, gene flow appears more important than pathogen-imposed selection in maintaining resistance diversity. In support of gene flow varying according to population connectivity, we found that population growth rates were lower in the intermediate- and low-connectivity host populations than in the well-connected host populations also in the absence of the pathogen. This suggests that increasing population isolation may also have other genetic consequences, such as lower genetic diversity and higher inbreeding depression, both of which may impact population growth rates45,46.
In theory, polymorphism in resistance within populations is maintained by costs of resistance in the absence of the pathogen, whereas under pathogen attack, the resistant hosts outperform the susceptible ones4. Hence, finding high levels of resistance diversity where pathogen impact has recently been negligible may appear contrary to expectations, and suggests dispersal to be critical for maintaining variation within host populations. Our metapopulation model explored scenarios under which spatial structure, disease dynamics, and life-history trade-offs could yield similar outcomes. We find that the shape of the host trade-off was the critical predictor of whether the simulations would qualitatively match the empirical results. Our results suggest that the costs of resistance are most likely to diminish as resistance increases. Diminishing costs mean that there is an initial large cost associated with resistance and therefore it is less beneficial when disease is rare. While fitness costs associated with resistance have been widely observed, determining the shape of trade-offs from empirical data is challenging, especially when trade-offs are close to linear or vary with environment, and it is impossible to determine trade-off shapes when only two host phenotypes are compared (as is often the case). However, experimental evolution of bacteria and phages has demonstrated that decelerating costs of resistance are possible4. In addition, our simulations suggest that the pattern detected in the empirical results is most likely to occur prior to the system reaching equilibrium and when metapopulation connectivity is assortative. The fact that the transient simulations dynamics tend to provide a better qualitative match to the empirical results does not imply that the resistance patterns detected in the archipelago will necessarily fade in the long-term (many simulations were qualitative matches at equilibrium), although our model indicates that this is a possibility. We think that it is interesting to note that the patterns we see are found for a wider range of parameter values under transient dynamics, but we get the same inference of the key characteristics that lead to the patterns we see. Whether or not the patterns are only transient is an empirical question.
Together, our results show how spatial fragmentation leading to the isolation of host populations drives the loss of diversity and increases host vulnerability to infectious diseases. By combining field population data with a controlled inoculation assay and a simulation model we demonstrate how spatial structure generates variation in the strength of pathogen-imposed selection, and thus provides a compelling example of how landscape fragmentation drives epidemiological and coevolutionary processes in nature.
Plantago lanceolata L. is a perennial monoecious ribwort plantain that reproduces both clonally via the production side rosettes, and sexually via wind pollination. Seeds drop close to the mother plant and usually form a long-term seed bank47. Podospharea plantaginis (Castagne; U. Braun and S. Takamatsu) (Erysiphales, Ascomycota) is an obligate biotrophic powdery mildew that infects only P. lanceolata and requires living host tissue through its life cycle48. It completes its life cycle as localized lesions on host leaves, only the haustorial feeding roots penetrating the leaf tissue to feed nutrients from its host. Infection causes significant stress for host plant and may increase the host mortality31. The interaction between P. lanceolata and P. plantaginis is strain-specific, whereby the same host genotype may be susceptible to some pathogen genotypes while being resistant to others49. The putative resistance mechanism includes two steps. First, resistance occurs when the host plant first recognizes the attacking pathogen and blocks its growth. When the first step fails and infection takes place, the host may mitigate infection development. Both resistance traits vary among host genotypes49.
Approximately 4000 P. lanceolata populations form a network covering an area of 50 × 70 km in the Åland Islands, SW of Finland. Disease incidence (0/1) in these populations has been recorded systematically every year in early September since 2001 by approximately 40 field assistants, who record the occurrence of the fungus P. plantaginis in the local P. lanceolata populations30. At this time, disease symptoms are conspicuous as infected plants are covered by white mycelia and conidia. The coverage (m2) of P. lanceolata in the meadows was recorded between 2001 and 2008 and is used as an estimate of host population size. In the field survey two technicians estimate Plantago population size by visually estimating how much ground/other vegetation P. lanceolata foliage covers (m2) in each meadow. The proportion of P. lanceolata plants in each population suffering from drought is also estimated annually in the survey. Data on average rainfall (mm) in July and August was estimated separately for each population using detailed radar-measured rainfall (obtained by Finnish Meteorological Institute) and it was available for years 2001–2008.
Host population connectivity (SH)27 for each local population i was computed with the formula that takes into account the area of host coverage (m²) of all host populations surveyed, denoted with (Aj), and their spatial location compared to other host populations. We assume that the distribution of dispersal distances from a location are described by negative exponential distribution. Under this assumption, the following formula (1) quantifies for a focal population i, the effect of all other host populations, taking into account their population sizes and how strongly they are connected through immigration to it:
here, dij is the Euclidian distance between populations i and j and 1/α equals the mean dispersal distance, which was set to be two kilometres based on results from a previous study16.
The annual survey data has demonstrated that P. plantaginis infects annually 2–16% of all host populations and persists as a highly dynamic metapopulation through extinctions and re-colonizations of local populations16. The number of host populations has remained relatively stable over the study period49. The first visible symptoms of P. plantaginis infection appear in late June as white-greyish lesions consisting of mycelium supporting the dispersal spores (conidia) that are carried by wind to the same or new host plants. Six to eight clonally produced generations follow one another in rapid succession, often leading to local epidemic with substantial proportion of the infected hosts by late summer within the host local population. Podosphaera plantaginis produces resting structures, chasmothecia, that appear towards the end of growing season in August–September31. Between 20% and 90% of the local pathogen populations go extinct during the winter, and thus the recolonization events play an important role in the persistence of the pathogen regionally16.
To examine whether the diversity and level of resistance vary among host populations depending on their degree of connectivity (SH) and disease history, we selected 20 P. lanceolata populations for an inoculation assay. These populations occur in different locations in the host network, and were selected based on their connectivity values (S H of selected populations was 37–110 in isolated and 237–336 in highly connected category, Fig. 1). We did not include host populations in the intermediate connectivity category that was used in the population dynamic analyses in the inoculation assay due to logistic constraints. Podosphaera plantaginis is an obligate biotrophic pathogen that requires living host tissue throughout its life cycle, and obtaining sufficient inoculum for experiments is extremely time and space consuming. In both isolated and highly connected categories, half of the populations (IDs 193, 260, 311, 313, 337, 507, 1821, 1999, 2818 and 5206) were healthy during the study years 2001–2014, while half of the populations (IDs 271, 294, 309, 321, 490, 609, 1553, 1556, 1676 and 1847) were infected by P. plantaginis for several years during the same period. We collected P. lanceolata seeds from randomly selected ten individual plants around the patch area from each host population in August 2014.
To acquire inoculum for the assay, we collected the pathogen strains as infected leaves, one leaf from ten plant individuals from four additional host populations (IDs 3301, 4684, 1784, and 3108) in August 2014. None of the pathogen populations were same as the sampled host populations and hence, the strains used in the assay all represent allopatric combinations. Both host and pathogen populations selected for the study were separated by at least two kilometres. The collected leaves supporting infection were placed in Petri dishes on moist filter paper and stored at room temperature until later use.
Seeds from ten mother plants from each population were sown in 2:1 mixture of potting soil and sand, and grown in greenhouse conditions at 20 ± 2 °C (day) and 16 ± 2 °C (night) with 16:8 L:D photoperiod. Due to the low germination rate of collected seeds, population 260 (isolated and healthy population) was excluded from the study. Seedlings of ten different mother plants were randomly selected among the germinated plants for each population (n = 190), and grown in individual pots until the plants were eight weeks old.
The pathogen strains were purified through three cycles of single colony inoculations and maintained on live, susceptible leaves on Petri dishes in a growth chamber 20 ± 2 °C with 16:8 L:D photoperiod. Every two weeks, the strains were transferred to fresh P. lanceolata leaves. Purified powdery mildew strains (M1–M4), one representing each allopatric population (3301, 4684, 1784 and 3108), were used for the inoculation assay. To produce enough sporulating fungal material, repeated cycles of inoculations were performed before the assay.
In order to study how the phenotypic resistance of hosts varies depending on population connectivity and infection history, we scored the resistance of 190 host genotypes, ten individuals from each study populations (n = 19), in an inoculation assay. Here, one detached leaf from each plant was exposed to a single pathogen strain (M1–M4) by brushing spores gently with a fine paintbrush onto the leaf. Leaves were placed on moist filter paper in Petri dishes and kept in a growth chamber at 20 ± 2 with a 16/8D photoperiod. All the inoculations were repeated on two individual Petri plates, leading to 760 host genotype—pathogen genotype combinations and a total of 1520 inoculations (19 populations * 10 plant genotypes * 4 pathogen strains * 2 replicates). We then observed and scored the pathogen infection on day 12 post inoculation, under dissecting microscope. The resulting plant phenotypic response was scored as 0 = susceptible (infection) when mycelium and conidia were observed on the leaf surface, and as 1 = resistance (no infection), when no developing lesions could be detected under a dissecting microscope. A genotype was defined resistant only if both inoculated replicates showed similar response (1), and susceptible if one or both replicates became infected (0).
To study how the pathogen infection influences on host population growth, we analyzed the relative change in host population size (m2) (defined as population size (t) − population size (t−1))/population size (t−1)) between consecutive years utilizing data from 2001 to 2008 in response to pathogen presence-absence status at t−1 (Supplementary Table 2). To assess whether this depends on host population connectivity, we estimated the separate effects of pathogen presence/absence in the previous year for connectivity categories—high-, low, and intermediate—that were based on the 0.2 and 0.8 quantiles of the host-connectivity values (Fig. 1A and Supplementary Figs. 1, 2). This allowed us to directly assess and compare the effect of the pathogen on host population growth in the extreme categories between isolated and highly connected host populations which were represented in the sampling for the inoculation study (Fig. 2).
As covariates, we included the proportion (0–100%) of dry host plants measured each year within each local population as well as data on the amount of rainfall at the summer months (June, July, and August) obtained from the satellite images, as these were suggested be relevant for this pathosystem in an earlier analysis16. Observations where the change in host population size, or the host population coverage had absolute values larger than their 0.99 quantiles in the whole data, were regarded as outliers and omitted from the analysis. Before the analyses, all the continuous covariates were scaled and centred, and the categorical variables were transformed into binary variables.
The relative changes in local host population size between consecutive years was analyzed by a Bayesian spatio-temporal statistical model that simultaneously considers the effects of a set of biologically meaningful predictors. The linear predictor thus consists of two parts (2,3):
where (beta) represents the correlation coefficients corresponding to the effects of environmental covariates, ({z}_{t}) corresponds to the spatiotemporal random effect, and ({X}_{t}) and ({A}_{t}) project these to the observation locations. For ({z}_{t}) we assume that the observations from a location in consecutive time points (t−1) and t are described by 1st order autoregressive process:
where ({w}_{t}) corresponds to spatially structured zero-mean random noise, for which a Matern covariance function is assumed. Statistical inference then targets jointly the covariate effects (beta), the temporal autocorrelation (varphi), and the hyperparameters describing the spatial autocorrelation in wt. From these the overall variance, as well as spatial range—a distance after which spatial autocorrelation ceases to be significant—can be inferred (Supplementary Fig. 3). For more detailed description of the structure of the statistical model and how to do efficient inference with it using R-INLA, we refer to refs. 16,50.
The phenotype composition of each study population was defined by individual plant responses to the four pathogen strains, where each response could be “susceptible = 0” or “resistant = 1”. For example, a phenotype “1111” refers to a plant resistant to all four pathogen strains. The diversity of distinct resistance phenotypes within populations was estimated using the Shannon diversity index as implemented in the vegan software package51. The Shannon diversity index for all four study groups was then analyzed using a linear model with class predictors population type (well-connected or isolated), infection history (healthy or infected), and their interaction.
To test whether host population resistance varied depending on connectivity (SH) and infection history, we analyzed the inoculation responses (0 = susceptible, 1=resistant) of each host-pathogen combination by using a logit mixed-effect model in the lme4 package52. The model included the binomial dependent variable (resistance-susceptible; 1/0), and class predictors population type (well-connected or isolated), infection history (healthy or infected), mildew strain (M1, M2, M3, and M4) and their interactions. Plant individual and population were defined as random effects, with plant genotype (sample) hierarchically nested under population. Model fit was assessed using chi-square tests on the log-likelihood values to compare different models and significant interactions, and the best model was selected based on AIC-values. P-values for regression coefficients were obtained by using the car package53. We ran all the analyses in R software54.
We model the ecological and co-evolutionary dynamics of host and pathogen metapopulations to understand key features of the experimental system that impact on the qualitative patterns observed. The structure and parameters in our model are therefore not estimated using experimental data, but rather are chosen to cover a range of possibilities (e.g., low vs high transmission rates, variation in trade-off shapes for fitness costs). We construct the metapopulations in two stages to account for relatively well and poorly connected demes. All demes are identical in quality (i.e., no differences in intrinsic birth or death rates between demes) and only differ in their connectivity. Our metapopulation consists of an outer network of 20 demes, equally spaced around the unit square (0.2 units apart), and a 7×7 inner lattice of demes at a minimum distance of 0.2 units from the outer network (Fig. 3A), giving a total of 69 demes. Demes that are separated by a Euclidean distance of at most 0.2 are then connected to each other. This means that populations near the centre of the metapopulation are highly connected, while those on the boundary of the metapopulation are poorly connected. This also has the effect of making connections between well and poorly connected demes assortative (i.e., well/poorly connected demes tend to be connected to well/poorly connected demes). We relax the assumption of assortativity in a second type of network by randomly reassigning connections between demes, while maintaining the same degree distribution. (i.e., the probability of two demes being connected is proportionate to their degree). While well connected demes still have more connections to other well connected demes than to poorly connected demes, they are not more likely to be connected to a well connected deme than by chance based on the degree distribution. In both types of network structure, we classify a deme as well-connected if it is in the top 20% of the degree distribution and poorly connected if it is in the bottom 20%.
We model the genetics using a multilocus gene-for-gene framework with haploid host and pathogen genotypes characterized by (L) biallelic loci, where 0 and 1 represent the presence and absence, respectively, of resistance and infectivity alleles. Host genotype (i) and pathogen genotype (j) are represented by binary strings: ({x}_{i}^{1}{x}_{i}^{2}ldots {x}_{i}^{L}) and ({y}_{j}^{1}{y}_{j}^{2}ldots {y}_{j}^{L}). Resistance acts multiplicatively such that the probability of host (i) being infected when challenged by pathogen (j) is ({Q}_{{ij}}={sigma }^{{d}_{{ij}}}), where (sigma) is the reduction in infectivity per effective resistance allele and ({d}_{{ij}}={sum }_{k=1}^{L}{x}_{i}^{k}big(1-{y}_{j}^{k}big)) is the number of effective resistance alleles (i.e., the number of loci where hosts have a resistance allele but pathogens do not have a corresponding infectivity allele). Hosts and pathogens with more resistance or infectivity alleles are assumed to pay higher fitness costs, ({c}_{H}left(iright)) eq. (4) and ({c}_{P}left(jright)) eq. (5) with:
where (0 , < , {c}_{H}^{1},; {c}_{P}^{1},le, 1) control the overall strength of the costs (i.e., the maximum proportional reduction in reproduction (hosts) or transmission rate (pathogens)) and ({c}_{H}^{2},; {c}_{P}^{2}in {{mathbb{R}}}_{ne 0}) control the shape of the trade-off. When ({c}_{H}^{2},; {c}_{P}^{2}, < , 0) the costs decelerate (increasing returns) and when ({c}_{H}^{2},; {c}_{P}^{2}, > , 0) the costs accelerate the costs accelerate (decreasing returns) (Supplementary Fig. 4). This formulation, therefore, allows for a wide-range of trade-off shapes that may occur in nature.
The dynamics of the (finite) host and pathogen populations are modelled stochastically using the tau-leap method with a fixed step size of (tau=1). For population (p), the mean host birth rate at time (t) for host (i) (6) is
where (a) is the maximum per-capita birth rate, (q) is the strength of density-dependent competition on births, ({N}_{p}left(tright)={S}_{i}^{p}left(tright)+{I}_{icirc }^{p}left(tright)) is the local host population size, ({S}_{i}^{p}left(tright)) and ({I}_{icirc }^{p}left(tright)={sum }_{j=1}^{n}{I}_{{ij}}^{p}left(tright)) are the local sizes of susceptible and infected individuals of genotype (i), and ({I}_{{ij}}^{p}left(tright)) is the local size of hosts of genotype (i) infected by pathogen (j). Host mutations occur at an average rate of ({mu }_{H}) per loci (limited to at most one mutation per time step), so that the mean number of mutations from host type (i) to ({i}^{{prime} }) is ({mu }_{H}{m}_{i{i}^{{prime} }}{B}_{i}^{p}left(tright)), where ({m}_{i{i}^{{prime} }}=1) if genotypes (i) and ({i}^{{prime} }) differ at exactly one locus, and is 0 otherwise.
The mean local mortalities for susceptible and infected individuals are (b{S}_{i}^{p}left(tright)) and (left(b+alpha right){I}_{{ij}}^{p}left(tright)), respectively, where (b) is the natural mortality rate and (alpha) is the disease-associated mortality rate. The average number of infected hosts that recover is (gamma {I}_{{ij}}^{p}left(tright)), where (gamma) is the recovery rate.
The mean number of new local infections of susceptible host type (i) by pathogen (j) eq. (7) is:
where (beta) is the baseline transmission rate and ({Y}_{j}^{p}left(tright)) is the local number of pathogen propagules following mutation and dispersal. Pathogen mutations occur in a similar manner to host mutations, with mutations from type (j) to ({j}^{{prime} }) occurring at rate ({mu }_{P}{m}_{j{j}^{{prime} }}{I}_{circ j}^{p}left(tright)) where ({mu }_{P}) is the mutation rate per loci (limited to at most one mutation per timestep) and ({I}_{circ j}^{p}left(tright)={sum }_{i=1}^{n}{I}_{{ij}}^{p}left(tright)) is the local number of pathogen (j.) Following mutation, the local number of pathogens of type (j) eq. (8) is:
Pathogen dispersal occurs following mutation at a rate of (rho) between connected demes, given by the adjacency matrix ({G}_{{pr}}), with ({G}_{varSigma p}) the total number of connections for deme (p). The mean local number of pathogen propagules following mutation and dispersal eq. (9) is therefore:
We focus our parameter sweep on: (i) the structure of the network (assortative or random connections); (ii) the strength (left({c}_{H}^{1},; {c}_{P}^{1}right)) and shape (left({c}_{H}^{2},; {c}_{P}^{2}right)) of the trade-offs; (iii) the transmission rate (left(beta right)); and (iv) the dispersal rate (left(rho right)), fixing the remaining parameters as described in Supplementary Table 1 (preliminary investigations suggested they had less of an impact on the qualitative outcome) and conducting 100 simulations per parameter set. For each simulation we initially seed all populations with the most susceptible host type and place the least infective pathogen type in one of the well-connected populations to minimize the risk of early extinction. We then solve the dynamics for 10,000 time steps (preliminary investigations indicated this was a sufficient period for the metapopulations to reach a quasi-equilibrium in terms of overall resistance). We calculate the average level of resistance (proportion of loci with a resistance allele) between time steps 4001 and 5000 (transient dynamics) and over the final 1000 time steps (long-term dynamics) for well and poorly connected demes, categorized according to whether the disease is present in (infected) or absent from (uninfected) the local population at a given time point and discarding simulations where the pathogen is driven globally extinct.
We compare the mean level of resistance in infected/uninfected poorly/well-connected populations across all simulations to the empirical results. We say that a simulation is a qualitative ‘match’ for the empirical findings if: (i) in poorly connected demes, the infected populations are on average at least 5% more resistant than uninfected populations; and (ii) in well-connected demes, the uninfected populations are on average at least 5% more resistant than infected populations. In other words, if ({R}_{{CS}}) is the mean resistance for a population with connectivity (C) ((C=W) and (C=P) for well and poorly connected demes, respectively) and infection status (S) ((S=U) and (S=I) for uninfected and infected populations, respectively), then a parameter set is a qualitative ‘match’ for the empirical findings if ({R}_{{WU}} > 1.05{R}_{{WI}}) and (1.05{R}_{{PI}}, > , 1.05{R}_{{PU}}). If these criteria are not met, then the parameter set is a qualitative ‘mismatch’ for the empirical findings. The model is not intended to be a replica of an empirical metapopulation, but rather is used to reveal the key factors which lead to qualitatively similar distributions of resistance and disease incidences observed in the study of the Åland islands. Hence, the purpose of the model is to determine which biological factors are likely to be crucial to the patterns observed herein.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
All data underlying this study are available at data are provided with this paper.
Code for population growth model is available at and code for the simulation model is available at
Haldane, J. Disease and evolution. Ricerca Sci. Suppl. A 19, 68–75 (1949).
Google Scholar 
Hamilton, W. D. Sex versus Non-Sex versus Parasite. Oikos 35, 282–290 (1980).
MATH  Article  Google Scholar 
Horns, F. & Hood, M. E. The evolution of disease resistance and tolerance in spatially structured populations. Ecol. Evol. 2, 1705–1711 (2012).
PubMed  PubMed Central  Article  Google Scholar 
Antonovics, J. & Thrall, P. H. The cost of resistance and the maintenance of genetic polymorphism in host-pathogen systems. Proc. R. Soc. B 257, 105–110 (1994).
ADS  Article  Google Scholar 
Anderson, R. M. & May, R. M. Infectious Diseases of Humans: Dynamics and Control (Oxford University Press, 1992).
Savary, S. et al. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 3, 430–439 (2019).
PubMed  Article  Google Scholar 
Fowler, N. L. & Clay, K. Environmental heterogeneity, fungal parasitism and the demography of the grass Stipa leucotricha. Oecologia 103, 55–62 (1995).
ADS  PubMed  Article  Google Scholar 
Gulland, F. M. The role of nematode parasites in Soay sheep (Ovis aries L.) mortality during a population crash. Parasitology 105, 493–503 (1992).
PubMed  Article  Google Scholar 
Hudson, P. J., Dobson, A. P. & Newborn, D. Prevention of population cycles by parasite removal. Science 282, 2256–2258 (1998).
ADS  CAS  PubMed  Article  Google Scholar 
Penczykowski, R. M., Walker, E., Soubeyrand, S. & Laine, A.-L. Linking winter conditions to regional disease dynamics in a wild plant-pathogen metapopulation. N. Phytol. 205, 1142–1152 (2015).
Article  Google Scholar 
Prendeville, H. R., Tenhumberg, B. & Pilson, D. Effects of virus on plant fecundity and population dynamics. N. Phytol. 202, 1346–1356 (2014).
Article  Google Scholar 
Gillespie, J. H. Natural selection for resistance to epidemics. Ecology 56, 493–495 (1975).
Article  Google Scholar 
Parratt, S. R., Numminen, E. & Laine, A.-L. Infectious disease dynamics in heterogeneous landscapes. Annu. Rev. Ecol. Evol. Syst. 47, 283–306 (2016).
Article  Google Scholar 
Ericson, L., Müller, W. J. & Burdon, J. J. 28-year temporal sequence of epidemic dynamics in a natural rust–host plant metapopulation. J. Ecol. 105, 701–713 (2017).
Article  Google Scholar 
Grenfell, B. & Harwood, J. (Meta)population dynamics of infectious diseases. Trends Ecol. Evol. 12, 395–399 (1997).
CAS  PubMed  Article  Google Scholar 
Jousimo, J. et al. Ecological and evolutionary effects of fragmentation on infectious disease dynamics. Science 344, 1289–1293 (2014).
ADS  CAS  PubMed  Article  Google Scholar 
Keeling, M. J. & Gilligan, C. A. Metapopulation dynamics of bubonic plague. Nature 407, 903–906 (2000).
ADS  CAS  PubMed  Article  Google Scholar 
MacPherson, A., Keeling, M. J. & Otto, S. P. Feedback between coevolution and epidemiology can help or hinder the maintenance of genetic variation in host-parasite models. Evolution 75, 582–599 (2021).
PubMed  Article  Google Scholar 
Mitchell, S. E., Rogers, E. S., Little, T. J. & Read, A. F. Host-parasite and genotype-by-environment interactions: temperature modifies potential for selection by a sterilizing pathogen. Evolution 59, 70–80 (2005).
PubMed  Article  Google Scholar 
Wolinska, J. & King, K. C. Environment can alter selection in host-parasite interactions. Trends Parasitol. 25, 236–244 (2009).
PubMed  Article  Google Scholar 
van Bergen, E. et al. The effect of summer drought on the predictability of local extinctions in a butterfly metapopulation. Conserv. Biol. 34, 1503–1511 (2020).
PubMed  Article  Google Scholar 
Boots, M. & Sasaki, A. ‘Small worlds’ and the evolution of virulence: Infection occurs locally and at a distance. Proc. Biol. Sci. 266, 1933–1938 (1999).
CAS  PubMed  PubMed Central  Article  Google Scholar 
Gandon, S. & Michalakis, Y. Local adaptation, evolutionary potential and host-parasite coevolution: Interactions between migration, mutation, population size and generation time. J. Evol. Biol. 15, 451–462 (2002).
Article  Google Scholar 
Real, L. A. & Biek, R. Spatial dynamics and genetics of infectious diseases on heterogeneous landscapes. J. R. Soc. Interface 4, 935–948 (2007).
PubMed  PubMed Central  Article  Google Scholar 
Lion, S. & Gandon, S. Evolution of spatially structured host-parasite interactions. J. Evol. Biol. 28, 10–28 (2015).
CAS  PubMed  Article  Google Scholar 
Bonte, D. et al. Costs of dispersal. Biol. Rev. Camb. Philos. Soc. 87, 290–312 (2012).
PubMed  Article  Google Scholar 
Hanski, I. Metapopulation Ecology (Oxford University Press, 1999).
Carlsson-Granér, U. & Thrall, P. H. Host resistance and pathogen infectivity in host populations with varying connectivity. Evolution 69, 926–938 (2015).
PubMed  Article  Google Scholar 
Höckerstedt, L. M., Siren, J. P. & Laine, A. Effect of spatial connectivity on host resistance in a highly fragmented natural pathosystem. J. Evol. Biol. 31, 844–852 (2018).
PubMed  PubMed Central  Article  Google Scholar 
Ojanen, S. P., Nieminen, M., Meyke, E., Pöyry, J. & Hanski, I. Long-term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia): Survey methods, data management, and long-term population trends. Ecol. Evol. 3, 3713–3737 (2013).
PubMed  PubMed Central  Article  Google Scholar 
Laine, A.-L. Resistance variation within and among host populations in a plant–pathogen metapopulation: Implications for regional pathogen dynamics. J. Ecol. 92, 990–1000 (2004).
Article  Google Scholar 
Laine, A.-L. Evolution of host resistance: Looking for coevolutionary hotspots at small spatial scales. Proc. R. Soc. B 273, 267–273 (2006).
PubMed  Article  Google Scholar 
Damgaard, C. Coevolution of a plant host-pathogen gene-for-gene system in a metapopulation model without cost of resistance or cost of virulence. J. Theor. Biol. 201, 1–12 (1999).
Thrall, P. H. & Burdon, J. J. Evolution of gene-for-gene systems in metapopulations: The effect of spatial scale of host and pathogen dispersal. Plant Pathol. 51, 169–184 (2002).
Article  Google Scholar 
Lindgren, F. & Rue, H. Bayesian spatial modelling with R-INLA. J. Stat. Softw. 63, 1–25 (2015).
Article  Google Scholar 
Ashby, B. & Boots, M. Multi-mode fluctuating selection in host-parasite coevolution. Ecol. Lett. 20, 357–365 (2017).
PubMed  Article  Google Scholar 
Gillespie, D. T. Approximate accelerated stochastic simulation of chemically reacting systems. J. Chem. Phys. 115, 1716–1733 (2001).
ADS  CAS  Article  Google Scholar 
Hudson, P., Rizzoli, A., Grenfell, B., Heesterbeek, H. & Dobson, A. The Ecology of Wildlife Diseases (Oxford University Press, 2002).
Dodds, P. & Thrall, P. Recognition events and host-pathogen co-evolution in gene-for-gene resistance to flax rust. Funct. Plant Biol. 36, 395–408 (2009).
CAS  PubMed  PubMed Central  Article  Google Scholar 
Mundt, C. C. Use of multiline cultivars and cultivar mixtures for disease management. Annu. Rev. Phytopathol. 40, 381–410 (2002).
CAS  PubMed  Article  Google Scholar 
Zhu, Y. et al. Genetic diversity and disease control in rice. Nature 406, 718–722 (2000).
ADS  CAS  PubMed  Article  Google Scholar 
Boots, M. & Haraguchi, Y. The evolution of costly resistance in host-parasite systems. Am. Nat. 153, 359–370 (1999).
PubMed  Article  Google Scholar 
Boots, M., White, A., Best, A. & Bowers, R. How specificity and epidemiology drive the coevolution of static trait diversity in hosts and parasites. Evolution 68, 1594–1606 (2014).
PubMed  PubMed Central  Article  Google Scholar 
Taylor, R. L. The genus Lithophragma (Saxifragaceae). Univ. Calif. Publ. Bot. 37, 1–89 (1965).
Google Scholar 
Hughes, A. R., Inouye, B. D., Johnson, M. T. J., Underwood, N. & Vellend, M. Ecological consequences of genetic diversity. Ecol. Lett. 11, 609–623 (2008).
PubMed  Article  Google Scholar 
Nonaka, E. et al. Scaling up the effects of inbreeding depression from individuals to metapopulations. J. Anim. Ecol. 88, 1202–1214 (2019).
PubMed  Article  Google Scholar 
Bos, M. Plantago: A Multidisciplinary Study Vol. 89 Ecological Studies 222–231 (Springer-Verlag, 1992).
Bushnell, W. R. The Powdery Mildews – A Comprehensive Treatise 1–12 (The American Phytological Society, 2002).
Laine, A.-L. Pathogen fitness components and genotypes differ in their sensitivity to nutrient and temperature variation in a wild plant-pathogen association. J. Evol. Biol. 20, 2371–2378 (2007).
PubMed  Article  Google Scholar 
Cameletti, M., Lindgren, F., Simpson, D. & Rue, H. Spatio-temporal modeling of particulate matter concentration through the SPDE approach. Adv. Stat. Anal. 97, 109–131 (2012).
MathSciNet  MATH  Article  Google Scholar 
Oksanen, J., Blanchet, G. & Kindt, R. Vegan: Community 711 Ecology Package. R package version 2.4-3. (2017).
Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Article  Google Scholar 
Fox, J. & Weisberg, S. An R Companion to Applied Regression, Second Edition (SAGE Publications, 2011).
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. (2016).
Höckersted, L. et al. Spatially structured eco-evolutionary dynamics in a host-pathogen interaction render isolated populations vulnerable to disease., (2022).
Höckersted, L. et al. ecoevotheory/Hockerstedt_et_al_2022: Final published version. (2022).
Download references
We would like to acknowledge Krista Raveala and Niko Vilenius for their assistance during the experimental work and all students who participated in annual metapopulation surveys. This work was funded by grants from the Academy of Finland (334276), and the European Research Council (Consolidator Grant RESISTANCE 724508) and SNF (310030_192770/1) to A.-L.L., and LUOVA Doctoral Programme funding to L.H. M.B. acknowledges the Natural Environment Research Council (NE/J009784/1), NIH/R01-GM122061-03 and NSF-DEB- 2011109 for support. B.A. is supported by the Natural Environment Research Council (grant no. NE/N014979/1).
Layla Höckerstedt
Present address: Finnish Meteorological Institute, FI-00101, Helsinki, Finland
These authors contributed equally: Layla Höckerstedt, Elina Numminen, Ben Ashby.
Organismal and Evolutionary Biology Research Program, Faculty of Biological and Environmental Sciences, 00014 University of Helsinki, Helsinki, Finland
Layla Höckerstedt, Elina Numminen & Anna-Liisa Laine
Department of Integrative Biology, University of California, Berkeley, CA, 94720, USA
Ben Ashby & Mike Boots
Department of Mathematics, Simon Fraser University, Burnaby, BC, V3H 5J5, Canada
Ben Ashby
Biosciences, University of Exeter, Penryn, TR10 9EZ, UK
Mike Boots
Department of Evolutionary Biology and Environmental Studies, University of Zürich, CH-8057, Zurich, Switzerland
Anna Norberg & Anna-Liisa Laine
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
A.-L.L., E.N., M.B., and L.H. conceived the ideas and designed the assay; L.H. conducted the experimental work and E.N., L.H., and A.N. analyzed the data. M.B. and B.A. developed and analyzed the simulation model. A.-L.L. and M.B. wrote the first draft of the manuscript. All the authors contributed to the writing of the manuscript and approved the final draft.
Correspondence to Anna-Liisa Laine.
The authors declare no competing interests.
Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit
Reprints and Permissions
Höckerstedt, L., Numminen, E., Ashby, B. et al. Spatially structured eco-evolutionary dynamics in a host-pathogen interaction render isolated populations vulnerable to disease. Nat Commun 13, 6018 (2022).
Download citation
Received: 04 February 2022
Accepted: 27 September 2022
Published: 13 October 2022
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Nature Communications (Nat Commun) ISSN 2041-1723 (online)
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.


About the Author

Leave a Reply

error: Konten dilindungi !!
adana eskort - eskişehir eskort - eskort mersin