Host Genotype and Coinfection Modify the Relationship of within and between Host Transmission
Variation in individual-level disease transmission is well documented, but the underlying causes of this variation are challenging to disentangle in natural epidemics. In general, within-host replication is critical in determining the extent to which infected hosts shed transmission propagules, but which factors cause variation in this relationship are poorly understood. Here, using a plant host, Plantago lanceolata, and the powdery mildew fungus Podosphaera plantaginis, we quantify how the distinct stages of within-host spread (autoinfection), spore release, and successful transmission to new hosts (alloinfection) are influenced by host genotype, pathogen genotype, and the coinfection status of the host. We find that within-host spread alone fails to predict transmission rates, as this relationship is modified by genetic variation in hosts and pathogens. Their contributions change throughout the course of the epidemic. Host genotype and coinfection had particularly pronounced effects on the dynamics of spore release from infected hosts. Confidently predicting disease spread from local levels of individual transmission, therefore, requires a more nuanced understanding of genotype-specific infection outcomes. This knowledge is key to better understanding the drivers of epidemiological dynamics and the resulting evolutionary trajectories of infectious disease.
Understanding the determinants of disease emergence and spread is one of the major challenges in disease biology (May et al. 2001; Yates et al. 2006; Gandon et al. 2013; Lively et al. 2014). Like most biological invasions, pathogen spread is characterized by the movement of small numbers of individuals through a spatially heterogeneous environment (Schreiber and Lloyd-Smith 2009; Hatcher et al. 2012). Small numbers of invading individuals mean that the initial phases of pathogen epidemics are highly stochastic, making the likelihood of successful transmission greatly influenced by the host a pathogen happens to infect (Schreiber and Lloyd-Smith 2009; Gandon et al. 2013; Hartfield and Alizon 2013). From the pathogen’s perspective, not all hosts are equal, and some hosts contribute more than others to the development of epidemics (Lloyd-Smith et al. 2005; Paull et al. 2012). Individual host heterogeneity in disease transmission is well recognized and has been addressed both theoretically (Matthews and Woolhouse 2005; Lloyd-Smith et al. 2006; Bifolchi et al. 2013) and empirically (Woolhouse et al. 1997; Kilpatrick et al. 2006; Cronin et al. 2010; Courtenay et al. 2014; VanderWaal et al. 2014). However, while it is clear that host heterogeneity in transmission is common, we currently know little about the genetic and environmental conditions that may lead some hosts to contribute disproportionately more to transmission (Ferrari et al. 2004; Fellous et al. 2012; Paull et al. 2012; Lass et al. 2013; Vale et al. 2013).
One facet of pathogen transmission that is currently poorly studied is the relationship between within-host colonization and between-host transmission, particularly how this relationship varies among individual hosts (Mideo et al. 2008). Once a pathogen has successfully infected a host, the process of within-host replication is expected to be critical for how much that infected host sheds transmission propagules that potentially lead to new infections. This link is a general feature of all host-parasite systems (Anderson and May 1982), regardless of the idiosyncrasies affecting exactly how within-host infection progresses and how transmission is achieved (Ebert and Weisser 1997; Quinn et al. 2000; Bowen and Walker 2005; Hughes et al. 2011). In plant-pathogen epidemics, infection dynamics can be divided into two distinct phases. During autoinfection, the pathogen spreads from the focal infection to the surrounding leaves within the same host plant. Variation in autoinfection leads to plants with varying numbers of infected leaves, and therefore autoinfection may also be seen as a measure of infection severity. During and following local spread within a host plant, alloinfection occurs, whereby the pathogen is transmitted to other hosts (Robinson 1976; Mundt 2009). A potentially important distinction between the autoinfection and alloinfection processes we describe here is that autoinfection involves the same pathogen transmission stage as alloinfection, while in most animal infections, they result from two different infection stages (e.g., in malaria, merozoites spread infection within hosts, and gametocytes transmit infection to other hosts; Bannister and Mitchell 2003; Schmidt-Hempel 2011).
Given that the window of time for a pathogen epidemic to develop is usually limited to certain environmental conditions (temperature, humidity) and host availability (Garrett et al. 2009), the optimal timing of auto- and alloinfection is crucial to the fate of the pathogen. Most epidemiological studies follow the rate of alloinfections, usually referred to as the transmission rate. Rates of autoinfection have typically received less attention, but some work has shown that host autoinfection can also vary (Lannou et al. 2008), with consequences for pathogen life-history and virulence evolution (van den Berg et al. 2013). Measuring autoinfection is generally easier than quantifying the extent of alloinfection, and traditionally, within-host pathogen replication is assumed to equate with transmission potential (Anderson and May 1982). Given that in most infections the symptoms arise as a consequence of the production of transmission propagules during autoinfection, we might expect hosts with severe disease symptoms to produce higher levels of transmission propagules (autoinfection) and ultimately cause more disease transmission (during alloinfection; Beldomenico and Begon 2010). However, there are many reasons why this relationship may not always hold, and understanding them would help pinpoint which diseased individuals should be targeted in order to maximize the efficacy of disease control.
To date, the effect of host and pathogen genetic diversity on the spread of infectious disease has received little experimental exploration and is currently one of the most pressing questions in disease ecology and evolution (Lively et al. 2014). Given that genetic variation (both of hosts and the pathogens that infect them) is well known to affect all aspects of infection outcomes (Bergelson et al. 2001; Wolinska and King 2009; Tack et al. 2012), the relationship between autoinfection and alloinfection is also likely to vary depending on the genetic composition of both hosts and pathogens. First, pathogen populations are highly variable in their ability to infect and cause disease (Tack et al. 2012), and trade-offs between the different life-history stages of infection may also vary with genotype and genotype-by-genotype interactions (Laine and Barrès 2013). Hence, there may be considerable variation among strains in their ability to grow within an infected host and transmit to new ones. Second, variation among hosts in how they resist becoming infected and mitigate infection development—as well as possible trade-offs between these traits—may further generate variation in autoinfection-alloinfection dynamics (Bergelson and Purrington 1996; Susi and Laine 2015). Finally, hosts are often found to harbor mixed infections, where the infection outcome cannot be predicted from the outcomes of single infections (Petney and Andrews 1998; Pedersen and Fenton 2007). Recognizing the variable nature of pathogen infections is important because mixed-genotype infections often result in increased parasite competition and virulence (Choisy and de Roode 2010; Alizon et al. 2013). Variation in the level of coinfection among individual plants is therefore likely to affect the extent of pathogen transmission propagule production during autoinfection and influence between-host transmission. Therefore, in addition to commonly studied genetic variation in host defenses (Simms and Triplett 1994; Roy and Kirchner 2000; Read et al. 2008; Råberg et al. 2009; Laine et al. 2011; Ayres and Schneider 2012), the level of coinfection may also generate individual host heterogeneity in infectiousness (Lass et al. 2013), potentially modifying the dynamics of infectious disease (Streicker et al. 2013; Susi et al. 2015a).
A clear understanding of what generates variation in the relationship between auto- and alloinfection calls for controlled experiments with two major features: (1) the level of autoinfection can be monitored in a nondestructive manner throughout the course of an infection across a large numbers of individual hosts, and (2) the resulting pathogen propagule shedding and spreading can be precisely measured for each individual host. In addition, a realistic understanding of auto- and alloinfection requires studying both single and mixed infections, as these are the conditions hosts will commonly experience in the wild. Fulfilling all these requirements is challenging, and so experiments such as the one we describe are understandably rare.
Here, we take a common garden approach using a plant-pathogen system to test how host genotype, pathogen genotype, coinfection, and time affect (1) the level of autoinfection within the host and (2) the relationship between autoinfection and alloinfection. The experimental work was carried out with powdery mildew, Podosphaera plantaginis, infecting the host plant, Plantago lanceolata. Our experiment used multiple host genotypes that were cloned into replicates and inoculated with two pathogen strains, either singly or as a coinfection. This host-pathogen interaction is highly amenable to ecological studies, as infection is visually conspicuous on the host surface and the disease cycle lacks extended latency periods. Hence, autoinfection can be visually directly quantified. With a combination of two spore trapping methods—spore traps and live susceptible leaves—we were able to disentangle the relationship between the shedding of transmission propagules and the actual spread of the infection. We monitored both autoinfection and alloinfection dynamics over the course of an epidemic to assess how this relationship changes over time. By genotyping the resulting infections, we were also able to identify the transmitted pathogen genotypes from coinfected host plants.
Material and Methods
Podosphaera plantaginis is a specialist powdery mildew naturally infecting Plantago lanceolata in the Åland archipelago, southwestern Finland. Plantago lanceolata is an obligate outcrossing perennial herb that reproduces both sexually and clonally via side rosettes. The epidemiological dynamics of Podosphaera in its large host population network have been studied since 2001 in the Åland Islands, southwest of Finland. These studies have demonstrated that this pathogen persists as a highly dynamic metapopulation (Jousimo et al. 2014). Visible signs of infection appear in late June in those host populations in which the pathogen has successfully overwintered as resting spores. The epidemic begins from these initial disease foci as the pathogen is transmitted by wind both within and among hosts via clonally produced dispersal spores, or conidia (Laine and Hanski 2006; Ovaskainen and Laine 2006). Some 6–8 clonally produced generations follow one another in quick succession, and as a consequence, infection spreads within (Ovaskainen and Laine 2006) and between host populations (Jousimo et al. 2014). By September, weather conditions turn unfavorable to disease transmission, and the epidemic spread ceases.
In the interaction between Plantago and Podosphaera, disease resistance is strain specific, with the same host genotype blocking infection by some strains of the pathogen while being susceptible to others. There is considerable variation among host individuals and populations in their degree of resistance. The obligate pathogen can establish only on susceptible hosts; hence, variation in resistance plays a fundamental role in determining disease dynamics (Jousimo et al. 2014). Local pathogen populations also support considerable genetic (Tollenaere et al. 2012) and phenotypic (Susi and Laine 2013) diversity. Coinfections, whereby two or more strains of Podosphaera simultaneously infect the same host, are common in the Åland metapopulation (Tollenaere et al. 2012).
Spore Trapping Experiment
To quantify the number of fungal spores released from each host plant and the number of successfully established new infections (a measure of transmission), we carried out a spore trapping experiment. We used eight generally susceptible (i.e., a wide range of Podosphaera strains, including the strains used in this study, are able to infect and sporulate in them) Plantago genotypes originating from five populations in Åland (IDs 4 [three plants], 511 [two plants], 1413 [one plant], 2220 [one plant], and 9031 [1 plant]) as focal plants. Host plants were collected as seed in August 2010 and stored in paper envelopes at room temperature. Seeds were germinated by placing them in 0.8-L pots in 50∶50 sand–potting soil mixture in greenhouse conditions with a 16L∶8D photoperiod at +22°C. Plants were cloned in the greenhouse according to the protocol described in Laine (2004), producing up to 24 ramets. Eight-week-old ramets were placed outside for 2 weeks of acclimation until the experiment was set up. We used two Podosphaera strains originating from Åland (strain 10 from population 2821 and strain 3 from population 877) that were infective on all host genotypes used in this experiment. The strains were purified and maintained on fresh susceptible Plantago leaves on Petri dishes in a growth chamber with a 16L∶8D photoperiod at .
To determine how host genetic background and pathogen treatment (single or coinfection) affect the relationship between autoinfection and alloinfection, we performed a spore trapping experiment under seminatural conditions. The experiment was set up in 2013 at the Lammi Biological Station (61°05′28″N, 25°03′90″E), where neither Plantago nor Podosphaera occur naturally; hence, environmental contamination by fungal spores was highly unlikely. In mid-July, the experimental plants were potted in -cm pots placed at a 1-meter radius from each other and inoculated with strain 3, strain 10, or a coinoculation of strains 3 and 10. The amount of inoculum (all spores brushed off from a 1-cm2 10-day-old sporulating lesion onto one leaf of the receiving plant) was the same for all plants, with the coinfected plants receiving half of the dose of the single genotype inoculum (i.e., all spores brushed from a 0.5-cm2 lesion of strain 3 and a 0.5-cm2 lesion of strain 10 onto one leaf of the receiving plant). We also included a control treatment with no pathogen spores to ensure there was no contamination between the plants. Four replicates of each plant treatment were used, and two replicates of each plant genotype and control treatment were used, resulting in 112 plants in total. The autoinfection rate of the infected plants was measured as the number of leaves on a plant that were infected with powdery mildew.
To quantify spore release and alloinfection, we conducted five spore trapping sessions at 20, 30, 40, 50, and 60 days postinfection (DPI). Our study was focused on quantifying short-distance transmission at a 5-cm distance from the infected host, as this is a relevant distance in the high-density populations of Plantago in Åland (Laine 2004), with most spores landing very close to the infection source (<10 cm; Tack et al. 2014). We used two types of traps: petroleum jelly–coated microscope slides (to quantify the number of spores released) and live detached Plantago leaves (to measure the number of spores that landed on and later germinated on a host, a measure of transmission potential). Four petroleum jelly–coated microscope slides were attached to wooden sticks at a 5-cm distance from the ground and placed between the infected leaves in a radial design to quantify spore release. Sixteen detached live leaves of known susceptible genotypes were attached to moist floral foam at a 5-cm distance from the focal plant to quantify alloinfection. These trap leaves were from five generally susceptible genotypes used in strain maintenance in the laboratory. The trapping period lasted 24 h, after which the traps were removed. The petroleum jelly traps were then kept in 5°C and subsequently examined under a microscope using four transect lines to count the released spores. The live leaves were placed on moist filter paper in a Petri dish and kept in a growth chamber. After 14 days, their infection status was monitored and the infected leaves were collected for subsequent genotyping. The level of autoinfection was assessed during each trapping session by counting the number of infected leaves. All control plants remained uninfected throughout the experiment. The infection status of the plants was monitored at 20 DPI on the leaf that had received the inoculation treatment. In total, 75 plants (79.2%) become infected. Plants that did not show visible signs of infection at 20 DPI (20.8%) were excluded from the analyses. To avoid unnecessary handling of the plants during the experiment, the total number of leaves in each plant was counted at the end of the experiment.
We genotyped the trap leaves infected at time points 40–60 DPI to ensure that there was no cross-contamination between plants and to determine which of the pathogen strains—or whether both—had successfully infected the trap leaves in the coinfection treatments. From each infected leaf, we cut the lesions, consisting of both host tissue and fungal material, into a 1.5-mL tube that was kept at −20°C until DNA was extracted using an E.Z.N.A. Plant Mini Kit (Omega Bio-Tek, Norcross, GA). Samples were genotyped using the 27 single nucleotide polymorphism panel, as in Tollenaere et al. (2012). The pathogen lines used in the experiment can be distinguished from each other as they differ at eight loci used in the genotyping panel. We classified leaves as coinfected if they showed the presence of two alleles in polymorphic loci (Tollenaere et al. 2012).
To test (1) how host genotype, pathogen genotype, coinfection, and time affect autoinfection, (2) how autoinfection correlates with spore release, (3) how autoinfection correlates with infection establishment, and (4) how spore release correlates with infection establishment, we performed four different analyses using generalized linear mixed models in SAS, version 9.2 (SAS Institute 2011), using the GLIMMIX procedure (Littell et al. 2006). Data used in the analyses are deposited in the Dryad Digital Repository:
To understand the relationship between autoinfection and alloinfection, we analyzed the proportion of infected trap leaves (as numerator) to all trap leaves (as denominator) around the focal plant as a response variable with a binomial distribution of errors. There was overdispersion in the data, and we chose a complementary log-log link function to fit the model. The level of autoinfection was included as a covariate; otherwise, the model variables are the same as described above. Finally, to understand the relationship between transmission potential and infection establishment, we analyzed how spore release affects the establishment of infection by having the proportion of infected trap leaves (as numerator) to all trap leaves (denominator) as a response variable, spore release (number of spores trapped on the four microscope slides) as a covariate, plant genotype and pathogen treatment as explanatory categorical variables, and time as DPI as a covariate. We used a binomial distribution of errors and a complementary log-log link function. To identify differences within significant main effects, we used post hoc comparisons by computing least squares means of the main effects in SAS, version 9.2, Proc Glimmix (Littell et al. 2006). In all analyses, nonsignificant interactions were excluded from the final models.
Factors Determining Autoinfection
The fraction of infected leaves increased with time, and there was weak negative correlation between the number of leaves in a plant and the fraction of diseased leaves (fig. A1; table A1). Plant genotype and pathogen treatment had a significant effect on the fraction of the diseased leaves (fig. A1; table A1). By the end of the experiment, none of the plants were saturated with infection: 40% or fewer of the leaves were infected (fig. A1).
The Dynamics of Pathogen Spore Release
Spore release (estimated as the number of spores landing on the petroleum jelly traps) changed throughout the experiment and across all treatments, peaking at 50 DPI (fig. 1; table 1). The number of spores released depended on whether the hosts were infected singly or coinfected, with coinfected hosts shedding more pathogen spores (fig. 1A; table 1). Post hoc comparisons revealed that while the overall number of spores released in the single infections did not differ, there were significant differences between the strains at 40 and 50 DPI (40 DPI, ; 50 DPI, ). This is especially clear at 50 DPI, where strain 3 released more spores than strain 10. The number of spores released from coinfected plants was highest at 50 DPI (fig. 1A). At 60 DPI, the overall number of spores released decreased and the differences between treatments diminished (fig. 1A).
The amount of spores released also depended on the host genotype. Averaged across all single infections and coinfections, some host genotypes (G1, G3, and G6) produced more spores than others throughout the epidemic (fig. 1B; table 1). The effect of host genetic variation on infectiousness is especially pronounced during the peak spore release at 50 DPI (fig. 1B), when the difference in spore release was the largest. Overall, the effect of pathogen treatment on spore release changed through time (significant interaction between time and treatment in table 1; fig. 1A). The amount of spores released from coinfected plants was highest at 30 and 50 DPI, while at 40 DPI, strain 10 had higher spore release than strain 3 and the coinfection treatment (fig. 1A). At the peak of epidemics, strain 3 released more spores than strain 10 in single infections (fig. 1A). At 60 DPI, the differences between treatments diminished (fig. 1A).
The Relationship between Autoinfection and Spore Release
One aim of our experiment was to explore how host genotype and coinfection status contributed to variation in the relationship between autoinfection and spore release. As expected, high autoinfection generally led to high levels of spore release (significant autoinfection term in table 1). However, while autoinfection and spore release are clearly correlated, the strength of this relationship changed over the course of the epidemic and was affected by both host and pathogen genotypes (fig. 2; table 1). Specifically, in hosts infected singly with strain 10, we observed a positive correlation between the rates of autoinfection and spore release, while in plants infected with strain 3, this correlation was still positive but weaker (fig. 2A–2C; table 1). In host plants coinfected with both strains, the relationship between autoinfection and spore release was highly variable (fig. 2A–2C; table 1). Across all single infection and coinfection treatments, there is also variation in the relationship that arises from different host genetic backgrounds. In some host genotypes (e.g., G1 and G2), high autoinfection yielded high spore release, whereas in other genotypes (e.g., G4, G5, and G6), the relationship was weaker (fig. 2D–2F; table 1).
The Establishment of New Infections
Consistent with what we observed for spore release, the establishment of infection by spores landing on trap leaves also peaked at 50 DPI (fig. 3A; table 2). We found that high autoinfection rates led to high rates of new infections becoming established, but we found no significant effects of host genotype and pathogen treatment (table 2). The relationship between the number of spores that were released from the infected plant and the establishment of new infections varied according to which pathogen strain was involved (; fig. 3B; table 3). Generally, when hosts were infected with strain 10, shedding high numbers of spores led to the establishment of many new infections, while this relationship was noticeably weaker when hosts were coinfected with strains 3 and 10 (fig. 3B; table 3). Genotyping of the infected trap leaves at 40–60 DPI from the coinfection treatment revealed that 40% were infected by strain 3, 40% were infected by strain 10, and 20% of new infections consisted of both strains.
Taken together, our results highlight the host’s heterogeneous contribution to the temporal dynamics of epidemics. We found that the contribution of a host individual to epidemics varies over time and correlates with its level of autoinfection as expected but, importantly, that it also depends on the host genetic background and whether it is coinfected or singly infected. We also tested a common assumption that rapid colonization of the host by the pathogen (autoinfection) is associated with the ability to spread infection between hosts (e.g., Robinson 1976). Our results indicate that host and pathogen genetic background can mediate the temporal dynamics of within- and between-host pathogen spread.
The Effect of Pathogen Strain on Transmission Dynamics
In this study, we found that after controlling for the level of autoinfection in our models, the overall spore release and infection establishment rates between the singly infecting strains did not differ. However, we observed significant differences between the strains at individual time points during the epidemics (fig. 1A) and, importantly, differences in how many spores were released at different levels of autoinfection (fig. 2). In plants infected with pathogen strain 10, there was a stronger correlation than in the case of plants infected with strain 3. This may indicate a difference in the latent period of the two strains: the weaker correlation for strain 3 spores may result from a longer period between autoinfection and spore shedding. This suggests that different pathogen strains may have different optimization strategies between within-host and between-host levels. The timing of transmission is essential in disease dynamics, as the time window for infection spread is usually limited by the environment, host life span, or host behavior (Hartfield and Alizon 2013). The differences in transmission at the early stages of the epidemics may have profound consequences for disease epidemiology and pathogen evolution, as rapid transmission is expected to lead to a greater share of the prevailing pathogen population (Day 2003) and therefore higher pathogen fitness (Elena 2001). We also found that in hosts infected with strain 10, the amount of spores released accurately predicted infection establishment, while in strain 3, this trend was not as clear, suggesting that there are differences between the strains in their spore quality (fig. 3B).
Host-Mediated Variation in Pathogen Dynamics within and between Hosts
Autoinfection and alloinfection dynamics have been approached theoretically in epidemiological modeling (Mundt and Leonard 1986; Willocquet and Savary 2004; Mideo et al. 2008), but empirical studies remain rare, and the host’s role in mediating the dynamics has been largely unexplored. We found heterogeneity in the contributions of different host genotypes to pathogen replication (fig. 1B), consistent with the idea that host genotype is critical for infection development (Salvaudon et al. 2008). Not only did host genotypes differ in the number of spores released, the relationship between spore release and infection establishment was also significantly determined by the host genetic background. One reason the relationship between autoinfection and spore release may vary is because hosts may differ in the types of defenses employed during infection and in the level of nutrients available for the pathogen (Laine 2007). Hosts have evolved various strategies to reduce the harm caused by the pathogen (Simms and Triplett 1994; Read et al. 2008; Råberg et al. 2009; Dodds and Rathjen 2010; Ayres and Schneider 2012). One strategy is to resist the pathogen by mitigating its growth, leading to clearance and recovery (Alexander 1992; Schmidt-Hempel 2011). Another strategy is to minimize the negative fitness effects of infection through damage limitation mechanisms that allow hosts to tolerate the presence of pathogens (Roy and Kirchner 2000; Ayres and Schneider 2012; Medzhitov et al. 2012). It remains unclear how the interaction between different host defense strategies may shape the relationship between autoinfection and alloinfection spread. For example, hosts investing heavily in costly resistance mechanisms could clear pathogens, leading to reduced transmission, but could also experience severe symptoms due to immunopathology (Graham et al. 2005; Jones and Dangl 2006). Conversely, hosts investing heavily in damage limitation mechanisms may show only mild disease symptoms (Poland et al. 2009) but would still tolerate the production of transmission propagules during autoinfection and contribute considerably to transmission during alloinfection (Vale et al. 2014). We used host genotypes where mildew was able to establish and grow, but it is possible that host genotypes differed in their ability to tolerate infection. Infection tolerance has not been systematically explored in this system but is an important line of future research. In general, an open question for future studies is how the relationship between autoinfection and spore release may vary between host genotypes with differential investment in resistance or tolerance strategies.
The Effect of Coinfection on Transmission Dynamics
We found that transmission dynamics were altered under coinfection; the overall highest spore shedding was observed in coinfected plants (fig. 1A). This result is in line with a study on mice where coinfection leads to higher helminth oviposition (Lass et al. 2013). Interestingly, in the early phase of the epidemic at 40 DPI, autoinfection rate was not a strong predictor of spore release in the coinfected plants—high spore release was observed at low autoinfection levels, while host plants with high levels of autoinfection released a low number of spores (fig. 2). As the epidemic progressed, a positive correlation between autoinfection and spore release emerged. It has also been suggested previously that the relationship between host disease severity and pathogen transmission may be altered in coinfection (Bremermann and Pickering 1983). Here, we found that the correlation between spore release and alloinfection was weaker under coinfection than in the singly infecting strains. This finding suggests that the accelerated spore production we observe under coinfection comes at a cost of spore quality. These effects of coinfection are likely to be relevant more broadly, as hosts are frequently infected by more than a single pathogen strain (Lopez-Villavicencio et al. 2007; Telfer et al. 2010; Balmer and Tanner 2011), and in some cases, it has been shown that coinfection may lead to increased reproduction of the pathogen (Lass et al. 2013).
Hosts are bottlenecks of pathogen genetic diversity and can therefore shape pathogen epidemics (Fellous et al. 2012; Tack et al. 2014) and evolution (Cisarovsky and Schmid-Hempel 2014). High levels of mixed infections can therefore maintain pathogen genetic diversity at the level of individual hosts, but it is poorly understood whether this diversity spreads to other hosts (Lively et al. 2014). Our genotyping revealed that 20% of the infections resulting from coinfected plants were also coinfected, suggesting the possibility of cotransmission, while the single infections established at equal probability from the coinfected source (both strains found singly in 40% of the infected leaf traps). Whether pathogens spread as single strains or are cotransmitted may have consequences for epidemiology, as infection success has been reported to increase in cotransmission (Karvonen et al. 2012) and is assumed to have further consequences for the evolution of virulence (Alizon 2013). In the case of vector-transmitted pathogens, simultaneous transmission of different strains and species is well documented (Pirone and Blanc 1996), but our study is among the few to investigate the frequency of cotransmission in airborne pathogens. Applied to natural populations, where multihost-multipathogen interaction networks are the norm (Pedersen and Fenton 2007), our results lend insight into a currently outstanding question in disease ecology and evolution (Lively et al. 2014): how pathogen interactions within hosts may regulate pathogen diversity at the population level.
Limitations of the Current Study
In this study, we found that the relationship between spore release and infection establishment was mediated only by pathogen genetic background, but we detected no effect of the focal host genotype (fig. 3). However, it is important to emphasize that the overall proportion of infected trap leaves was low, 0.014. While this is comparable to the estimated success of single bouts of transmission (Ovaskainen and Laine 2006), it may have reduced our power to detect the effects of these potential sources of variation. The size of the current experimental design was determined by the need to study spore release, autoinfection, and alloinfection. In order to increase statistical power to fully disentangle the sources of variation on successful disease spread during alloinfection, future studies should focus on the realized transmission on large numbers of host genotypes representing different resistance backgrounds and large numbers of pathogen stains with different infection profiles. Moreover, as previous studies have found evidence of pathogen local adaptation to host resistance (Laine 2005, 2008), it would be highly relevant to determine how coevolutionary dynamics may alter the relationship between auto- and alloinfection.
We detected important effects of host genotype and coinfection status on the temporal dynamics of spore release, which reinforces the need to identify coinfection in wild hosts as a potential risk factor of disease spread (Mideo et al. 2008; Susi et al. 2015a). Our results also suggest that the degree to which a host plant sheds fungal spores onto its own leaves is not necessarily a good predictor of spread to neighboring host plants: we found the correlation between autoinfection and alloinfection to change over the course of our experimental epidemic and to vary particularly in strength according to the host coinfection status. This result has important implications for the control of disease. In trying to delay the spread of infection, it would seem logical to remove or treat hosts showing the clearest signs of disease. While this may be true generally, our results suggest that timing is crucial, as a strong positive relationship between autoinfection and spore release became apparent only as the epidemic progressed. When designing disease management efforts, it is therefore important to understand the relationship between autoinfection, spore shedding, and alloinfection and how the relationship between them may vary according to the genetic and environmental context experienced by hosts.
The Lammi Biological Station is acknowledged for hosting the experimental study, and S. Hartikainen and P. Hyttinen are acknowledged for their help in the field. L. Ruokolainen is acknowledged for statistical advice. N. Mideo and R. Penczykowski are acknowledged for helpful comments on an early version of the manuscript. B. Barrès helped analyze the genetic samples, and the Institute of Biotechnology and the Institute for Molecular Medicine at the University of Helsinki are acknowledged for carrying out the DNA extractions and for genotyping samples, respectively. H.S. is funded by the LUOVA Doctoral School in Wildlife Biology Research, and P.F.V. is funded by a strategic award from the Wellcome Trust for the Centre for Immunity, Infection and Evolution (
Appendix. Variation in Autoinfection of Podosphaera plantaginis
Alexander, H. M. 1992. Evolution of disease resistance in natural plant populations. Pages 326–344 in R. S. Fritz and E. L. Simms, eds. Plant resistance to herbivores and pathogens ecology, evolution, and genetics. University of Chicago Press, Chicago. Alizon, S. 2013. Parasite co-transmission and the evolutionary epidemiology of virulence. Evolution 67:921–933. Alizon, S., J. C. de Roode, and Y. Michalakis. 2013. Multiple infections and the evolution of virulence. Ecology Letters 16:556–567. Anderson, R. M., and R. M. May. 1982. Coevolution of hosts and parasites. Parasitology 85:411–426. Ayres, J. S., and D. S. Schneider. 2012. Tolerance of infections. Annual Review of Immunology 30:271–294. Balmer, O., and M. Tanner. 2011. Prevalence and implications of multiple-strain infections. Lancet Infectious Diseases 11:868–887. Bannister, L., and G. Mitchell. 2003. The ins, outs and roundabouts of malaria. Trends in Parasitology 19:209–213. Beldomenico, P. M., and M. Begon. 2010. Disease spread, susceptibility and infection intensity: vicious circles? Trends in Ecology and Evolution 25:21–27. Bergelson, J., M. Kreitman, E. A. Stahl, and D. C. Tian. 2001. Evolutionary dynamics of plant R-genes. Science 292:2281–2285. Bergelson, J., and C. B. Purrington. 1996. Surveying patterns in the cost of resistance in plants. American Naturalist 148:536–558. Bifolchi, N., R. Deardon, and Z. Feng. 2013. Spatial approximations of network-based individual level infectious disease models. Spatial and Spatio-temporal Epidemiology 6:59–70. Bowen, D. G., and C. M. Walker. 2005. Adaptive immune responses in acute and chronic hepatitis C virus infection. Nature 436:946–952. Bremermann, H. J., and J. Pickering. 1983. A game-theoretical model of parasite virulence. Journal of Theoretical Biology 100:411–426. Choisy, M., and J. C. de Roode. 2010. Mixed infections and the evolution of virulence: effects of resource competition, parasite plasticity, and impaired host immunity. American Naturalist 175:E105–E118. Cisarovsky, G., and P. Schmid-Hempel. 2014. Few colonies of the host Bombus terrestris disproportionately affect the genetic diversity of its parasite, Crithidia bombi. Infection, Genetics and Evolution 21:192–197. Courtenay, O., C. Carson, L. Calvo-Bado, L. M. Garcez, and R. J. Quinnell. 2014. Heterogeneities in Leishmania infantum infection: using skin parasite burdens to identify highly infectious dogs. PLoS Neglected Tropical Diseases 8:26. Cronin, J. P., M. E. Welsh, M. G. Dekkers, S. T. Abercrombie, and C. E. Mitchell. 2010. Host physiological phenotype explains pathogen reservoir potential. Ecology Letters 13:1221–1232. Day, T. 2003. Virulence evolution and the timing of disease life-history events. Trends in Ecology and Evolution 18:113–118. Dodds, P. N., and J. P. Rathjen. 2010. Plant immunity: towards an integrated view of plant-pathogen interactions. Nature Reviews Genetics 11:539–548. Ebert, D., and W. W. Weisser. 1997. Optimal killing for obligate killers: the evolution of life histories and virulence of semelparous parasites. Proceedings of the Royal Society B: Biological Sciences 264:985–991. Elena, S. F. 2001. Evolutionary history conditions the timing of transmission in vesicular stomatitis virus. Infection, Genetics and Evolution 1:151–159. Fellous, S., A. B. Duncan, E. Quillery, P. F. Vale, and O. Kaltz. 2012. Genetic influence on disease spread following arrival of infected carriers. Ecology Letters 15:186–192. Ferrari, N., I. M. Cattadori, J. Nespereira, A. Rizzoli, and P. J. Hudson. 2004. The role of host sex in parasite dynamics: field experiments on the yellow-necked mouse Apodemus flavicollis. Ecology Letters 7:88–94. Gandon, S., M. E. Hochberg, R. D. Holt, and T. Day. 2013. What limits the evolutionary emergence of pathogens? Philosophical Transactions of the Royal Society B: Biological Sciences 368:20120086. Garrett, K. A., L. N. Zuniga, E. Roncal, G. A. Forbes, C. C. Mundt, Z. Su, and R. J. Nelson. 2009. Intraspecific functional diversity in hosts and its effect on disease risk across a climatic gradient. Ecological Applications 19:1868–1883. Graham, A. L., J. E. Allen, and A. F. Read. 2005. Evolutionary causes and consequences of immunopathology. Annual Review of Ecology, Evolution, and Systematics 36:373–397. Hartfield, M., and S. Alizon. 2013. Introducing the outbreak threshold in epidemiology. PLoS Pathogens 9:e1003277. Hatcher, M. J., J. T. A. Dick, and A. M. Dunn. 2012. Disease emergence and invasions. Functional Ecology 26:1275–1287. Hughes, D. P., S. B. Andersen, N. L. Hywel-Jones, W. Himaman, J. Billen, and J. J. Boomsma. 2011. Behavioral mechanisms and morphological symptoms of zombie ants dying from fungal infection. BMC Ecology 11:13. Jones, J. D. G., and J. L. Dangl. 2006. The plant immune system. Nature 444:323–329. Jousimo, J., A. J. M. Tack, O. Ovaskainen, T. Mononen, H. Susi, C. Tollenaere, and A.-L. Laine. 2014. Ecological and evolutionary effects of fragmentation on infectious disease dynamics. Science 344:1289–1293. Karvonen, A., C. Rellstab, K.-R. Louhi, and J. Jokela. 2012. Synchronous attack is advantageous: mixed genotype infections lead to higher infection success in trematode parasites. Proceedings of the Royal Society B: Biological Sciences 279:171–176. Kilpatrick, A. M., P. Daszak, M. J. Jones, P. P. Marra, and L. D. Kramer. 2006. Host heterogeneity dominates West Nile virus transmission. Proceedings of the Royal Society B: Biological Sciences 273:2327–2333. Laine, A.-L. 2004. Resistance variation within and among host populations in a plant-pathogen metapopulation: implications for regional pathogen dynamics. Journal of Ecology 92:990–1000. ———. 2005. Spatial scale of local adaptation in a plant-pathogen metapopulation. Journal of Evolutionary Biology 18:930–938. ———. 2007. Pathogen fitness components and genotypes differ in their sensitivity to nutrient and temperature variation in a wild plant-pathogen association. Journal of Evolutionary Biology 20:2371–2378. ———. 2008. Temperature-mediated patterns of local adaptation in a natural plant-pathogen metapopulation. Ecology Letters 11:327–337. Laine, A.-L., and B. Barrès. 2013. Epidemiological and evolutionary consequences of life-history trade-offs in pathogens. Plant Pathology 62:96–105. Laine, A.-L., J. J. Burdon, P. N. Dodds, and P. H. Thrall. 2011. Spatial variation in disease resistance: from molecules to metapopulations. Journal of Ecology 99:96–112. Laine, A.-L., and I. Hanski. 2006. Large-scale spatial dynamics of a specialist plant pathogen in a fragmented landscape. Journal of Ecology 94:217–226. Lannou, C., S. Soubeyrand, L. Frezal, and J. Chadœuf. 2008. Autoinfection in wheat leaf rust epidemics. New Phytologist 177:1001–1011. Lass, S., P. J. Hudson, J. Thakar, J. Saric, E. Harvill, R. Albert, and S. E. Perkins. 2013. Generating super-shedders: co-infection increases bacterial load and egg production of a gastrointestinal helminth. Journal of the Royal Society Interface 10:20120588. Littell, R. C., G. A. Milliken, W. W. Stroup, R. D. Wolfinger, and O. Schabenberger. 2006. SAS for mixed models. SAS Institute, Cary, NC. Lively, C. M., J. C. De Roode, M. A. Duffy, A. L. Graham, and B. Koskella. 2014. Interesting open questions in disease ecology and evolution. American Naturalist 184(suppl.):S1–S8. Lloyd-Smith, J. O., S. J. Schreiber, and W. M. Getz. 2006. Moving beyond averages: individual-level variation in disease transmission. Contemporary Mathematics 410:235–258. Lloyd-Smith, J. O., S. J. Schreiber, P. E. Kopp, and W. M. Getz. 2005. Superspreading and the effect of individual variation on disease emergence. Nature 438:355–359. Lopez-Villavicencio, M., O. Jonot, A. Coantic, M. E. Hood, J. Enjalbert, and T. Giraud. 2007. Multiple infections by the anther smut pathogen are frequent and involve related strains. PLoS Pathogens 3:e176. Matthews, L., and M. Woolhouse. 2005. New approaches to quantifying the spread of infection. Nature Reviews Microbiology 3:529–536. May, R. M., S. Gupta, and A. R. McLean. 2001. Infectious disease dynamics: what characterizes a successful invader? Philosophical Transactions of the Royal Society B: Biological Sciences 356:901–910. Medzhitov, R., D. S. Schneider, and M. P. Soares. 2012. Disease tolerance as a defense strategy. Science 335:936–941. Mideo, N., S. Alizon, and T. Day. 2008. Linking within- and between-host dynamics in the evolutionary epidemiology of infectious diseases. Trends in Ecology and Evolution 23:511–517. Mundt, C. C. 2009. Importance of autoinfection to the epidemiology of polycyclic foliar disease. Phytopathology 99:1116–1120. Mundt, C. C., and K. J. Leonard. 1986. Analysis of factors affecting disease increase and spread in mixtures of immune and susceptible plants in computer-simulated epidemics. Phytopathology 76:832–840. Ovaskainen, O., and A.-L. Laine. 2006. Inferring evolutionary signals from ecological data in a plant-pathogen metapopulation. Ecology 87:880–891. Paull, S. H., S. Song, K. M. McClure, L. C. Sackett, A. M. Kilpatrick, and P. T. J. Johnson. 2012. From superspreaders to disease hotspots: linking transmission across hosts and space. Frontiers in Ecology and the Environment 10:75–82. Pedersen, A. B., and A. Fenton. 2007. Emphasizing the ecology in parasite community ecology. Trends in Ecology and Evolution 22:133–139. Petney, T. N., and R. H. Andrews. 1998. Multiparasite communities in animals and humans: frequency, structure and pathogenic significance. International Journal for Parasitology 28:377–393. Pirone, T. P., and S. Blanc. 1996. Helper-dependent vector transmission of plant viruses. Annual Review of Phytopathology 34:227–247. Poland, J. A., P. J. Balint-Kurti, R. J. Wisser, R. C. Pratt, and R. J. Nelson. 2009. Shades of gray: the world of quantitative disease resistance. Trends in Plant Science 14:21–29. Quinn, T. C., M. J. Wawer, N. Sewankambo, D. Serwadda, C. J. Li, F. Wabwire-Mangen, M. O. Meehan, T. Lutalo, and R. H. Gray. 2000. Viral load and heterosexual transmission of human immunodeficiency virus type 1. New England Journal of Medicine 342:921–929. Råberg, L., A. L. Graham, and A. F. Read. 2009. Decomposing health: tolerance and resistance to parasites in animals. Philosophical Transactions of the Royal Society B: Biological Sciences 364:37–49. Read, A. F., A. L. Graham, and L. Råberg. 2008. Animal defenses against infectious agents: is damage control more important than pathogen control? PLoS Biology 6:2638–2641. Robinson, R. R. 1976. Plant pathosystems. Advanced series in agricultural sciences. Vol. 3. Springer, Berlin. Roy, B. A., and J. W. Kirchner. 2000. Evolutionary dynamics of pathogen resistance and tolerance. Evolution 54:51–63. Salvaudon, L., T. Giraud, and J. A. Shykoff. 2008. Genetic diversity in natural populations: a fundamental component of plant-microbe interactions. Current Opinion in Plant Biology 11:135–143. SAS Institute. 2011. The SAS system for Windows. Release 9.2. SAS Institute, Cary, NC. Schmidt-Hempel, P. 2011. Evolutionary parasitology: the integrated study of infections, immunology, ecology, and genetics. Oxford University Press, Oxford. Schreiber, S. J., and J. O. Lloyd-Smith. 2009. Invasion dynamics in spatially heterogeneous environments. American Naturalist 174:490–505. Simms, E. L., and J. Triplett. 1994. Costs and benefits of plant-responses to disease—resistance and tolerance. Evolution 48:1973–1985. Streicker, D. G., A. Fenton, and A. B. Pedersen. 2013. Differential sources of host species heterogeneity influence the transmission and control of multihost parasites. Ecology Letters 16:975–984. Susi, H., B. Barres, P. F. Vale, and A.-L. Laine. 2015a. Coinfection alters population dynamics of infectious disease. Nature Communications 6:5975. Susi, H., and A.-L. Laine. 2013. Pathogen life-history trade-offs revealed in allopatry. Evolution 67:3362–3370. ———. 2015. The effectiveness and costs of pathogen resistance strategies in a perennial plant. Journal of Ecology 103:303–315. Susi, H., P. F. Vale, and A.-L. Laine. 2015b. Data from: Host genotype and coinfection modify the relationship of within and between host transmission. Dryad Digital Repository, http://dx.doi.org/10.5061/dryad.v150f. Tack, A. J. M., J. Hakala, T. Petäjä, M. Kulmala, and A.-L. Laine. 2014. Genotype and spatial structure shape pathogen dispersal and disease dynamics at small spatial scales. Ecology 95:703–714. Tack, A. J. M., P. H. Thrall, L. G. Barrett, J. J. Burdon, and A.-L. Laine. 2012. Variation in infectivity and aggressiveness in space and time in wild host-pathogen systems: causes and consequences. Journal of Evolutionary Biology 25:1918–1936. Telfer, S., X. Lambin, R. Birtles, P. Beldomenico, S. Burthe, S. Paterson, and M. Begon. 2010. Species interactions in a parasite community drive infection risk in a wildlife population. Science 330:243–246. Tollenaere, C., H. Susi, J. Nokso-Koivisto, P. Koskinen, A. Tack, P. Auvinen, L. Paulin, M. J. Frilander, R. Lehtonen, and A.-L. Laine. 2012. SNP design from 454 sequencing of Podosphaera plantaginis transcriptome reveals a genetically diverse pathogen metapopulation with high levels of mixed-genotype infection. PLoS ONE 7:e52492. Vale, P. F., M. Choisy, and T. J. Little. 2013. Host nutrition alters the variance in parasite transmission potential. Biology Letters 9:20121145. Vale, P. F., A. Fenton, and S. Brown. 2014. Limiting damage during infection: lessons from infection tolerance for novel therapeutics. PLoS Biology 12:e1001769. van den Berg, F., S. Gaucel, C. Lannou, C. A. Gilligan, and F. van den Bosch. 2013. High levels of auto-infection in plant pathogens favour short latent periods: a theoretical approach. Evolutionary Ecology 27:409–428. VanderWaal, K. L., E. R. Atwill, L. A. Isbell, and B. McCowan. 2014. Linking social and pathogen transmission networks using microbial genetics in giraffe (Giraffa camelopardalis). Journal of Animal Ecology 83:406–414. Willocquet, L., and S. Savary. 2004. Modelling plant disease epidemics in a spatial hierarchy. Phytopathology 94:S110. Wolinska, J., and K. C. King. 2009. Environment can alter selection in host-parasite interactions. Trends in Parasitology 25:236–244. Woolhouse, M. E. J., C. Dye, J. F. Etard, T. Smith, J. D. Charlwood, G. P. Garnett, P. Hagan, et al. 1997. Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proceedings of the National Academy of Sciences of the USA 94:338–342. Yates, A., R. Antia, and R. R. Regoes. 2006. How do pathogen evolution and host heterogeneity interact in disease emergence? Proceedings of the Royal Society B: Biological Sciences 273:3075–3083.
Associate Editor: Scott L. Nuismer
Editor: Judith L. Bronstein