key: cord-0992653-5cno180i
authors: Susi, Hanna; Sallinen, Suvi; Laine, Anna‐Liisa
title: Coinfection with a virus constrains within‐host infection load but increases transmission potential of a highly virulent fungal plant pathogen
date: 2022-03-08
journal: Ecol Evol
DOI: 10.1002/ece3.8673
sha: 0af743a0ac4d19141dac48202bfa2c54b6db43f4
doc_id: 992653
cord_uid: 5cno180i

The trade‐off between within‐host infection rate and transmission to new hosts is predicted to constrain pathogen evolution, and to maintain polymorphism in pathogen populations. Pathogen life‐history stages and their correlations that underpin infection development may change under coinfection with other parasites as they compete for the same limited host resources. Cross‐kingdom interactions are common among pathogens in both natural and cultivated systems, yet their impacts on disease ecology and evolution are rarely studied. The host plant Plantago lanceolata is naturally infected by both Phomopsis subordinaria, a seed killing fungus, as well as Plantago lanceolata latent virus (PlLV) in the Åland Islands, SW Finland. We performed an inoculation assay to test whether coinfection with PlLV affects performance of two P. subordinaria strains, and the correlation between within‐host infection rate and transmission potential. The strains differed in the measured life‐history traits and their correlations. Moreover, we found that under virus coinfection, within‐host infection rate of P. subordinaria was smaller but transmission potential was higher compared to strains under single infection. The negative correlation between within‐host infection rate and transmission potential detected under single infection became positive under coinfection with PlLV. To understand whether within‐host and between‐host dynamics are correlated in wild populations, we surveyed 260 natural populations of P. lanceolata for P. subordinaria infection occurrence. When infections were found, we estimated between‐hosts dynamics by determining pathogen population size as the proportion of infected individuals, and within‐host dynamics by counting the proportion of infected flower stalks in 10 infected plants. In wild populations, the proportion of infected flower stalks was positively associated with pathogen population size. Jointly, our results suggest that the trade‐off between within‐host infection load and transmission may be strain specific, and that the pathogen life‐history that underpin epidemics may change depending on the diversity of infection, generating variation in disease dynamics.

To grow, multiply, and transmit, pathogens obtain resources from their host, and theoretically the resulting within-host infection rate is expected to be proportional to the gained resources. However, too high within-host infection rate may come at the cost of increased virulence, that is measurable damage caused to the host (Surico, 2013) , thereby incurring a cost for the pathogen as decreased transmission to other hosts (Blanquart et al., 2016) . Thus, evolutionary theory and experimental studies (de Roode et al., 2008) have established that selection should favor intermediate levels of withinhost infection rates. This trade-off between within-host infection rate and transmission has been proposed to maintain polymorphism in pathogen populations, and to prevent the rise of highly virulent pathogens (Frank, 1992) . Trade-offs have been sought as an evolutionary solution to limit disease epidemics and the emergence of pathogen strains with extremely high within-host infection rate (Zhan et al., 2015) . However, insight on how pathogen within-host infection rate links to transmission during epidemics where pathogens may encounter variation in both biotic and abiotic conditions (Blanquart et al., 2016; Dutta et al., 2021; Susi & Laine, 2013) has remained limited (Acevedo et al., 2019) . In the wild, the limited evidence for trade-offs may be explained by spatial (Osnas et al., 2015) and host-mediated processes (Kubinak et al., 2012) . The trade-offs restraining within-host infection rate may also occur between other traits, that is, adaptation to abiotic conditions (Mboup et al., 2012) or be context-dependent and become evident in stressful environments (Susi & Laine, 2013) .

The drivers of disease evolution and epidemics are rarely limited to the interplay of one host and one pathogen, as in the wild most infections occur as coinfections whereby multiple pathogens are simultaneously infecting the same host (Telfer et al., 2010; Tollenaere et al., 2016) . Coinfection may fundamentally change pathogen's host exploitation strategy in order to outcompete other pathogens sharing the same limited resource (Alizon & van Baalen, 2008; Alizon et al., 2009; de Roode et al., 2005) . Thus, it has been suggested that coinfection is an important driver of disease evolution (Alizon & van Baalen, 2008; Alizon et al., 2009) . Experimental approaches have measured increased within-host infection rate (Bell et al., 2006) and transmission under coinfection but there are also exceptions to this trend (Orton & Brown, 2016) .

Overall, it is well established that the pathogen within-host infection rate may change under coinfection, but studies explicitly testing trade-offs between within-host infection rate and transmission under coinfection are rare, and evidence remains mixed (Sacristan & Garcia-Arenal, 2008; Suffert et al., 2016) . Furthermore, coinfection experiments have often been conducted using strains of the same pathogen species, although interspecific interactions among pathogen species are likely to play an important role, as individual hosts often support diverse pathogen assemblages (Dallas et al., 2019; Susi et al., 2019; Telfer et al., 2010) .

While intraspecific coinfection is a pre-requisite for outcrossing for many pathogens (Suffert et al., 2016) , theory predicts the intensity of competition to increase as relatedness decreases (Alizon et al., 2013) . Across plants (Tollenaere et al., 2016 (Tollenaere et al., , 2017 , animals (Telfer et al., 2010) and humans (Chen et al., 2020; Lawn et al., 2006) interkingdom coinfections are common, and they are often suggested to have serious consequences for disease epidemics and severity. In particular, it is becoming increasingly clear that viruses are ubiquitous in nature Munson-McGee et al., 2018) , although their true diversity and prevalence in natural populations has been under-estimated for a long time (Roossinck et al., 2015; Wren et al., 2006) . The ecological roles of viruses are still poorly understood (Alexander et al., 2014; Roossinck, 2010) , but they have the potential to interact with other pathogen species via competition for shared host resources, and via shared effects on host immunity (Huang et al., 2019; Uehling et al., 2017) . Thus, it is vital to test how coinfection with pathogens from distant taxa may influence withinhost infection rate and transmission, and their potential trade-offs.

Phomopsis subordinaria is a castrating pathogen that infects its hosts through seed stalks. Potential trade-offs between within-host infection rate and transmission would translate in field epidemics as reduced transmission in the populations with high within-host infection rates. Here, we investigate the association between within-host dynamics and between-host dynamics of P. subordinaria by surveying 260 host plant (Plantago lanceolata) populations in the Åland Islands, south-west Finland. We combined data from the field measured disease with laboratory trials. In the laboratory, we challenged P. subordinaria strains with Plantago lanceolata latent virus (PlLV) to understand how cross-kingdom interactions affect within-host infection rate and transmission potential as well as their potential trade-off. Specifically, we ask (1) How common is P. subordinaria in the Åland Islands, and is there natural variation in the within-host infection dynamics in natural P. lanceolata populations, and (2) Is there an association between within-host and between-host dynamics in P. subordinaria? In a laboratory experiment, we tested: (3) Is there a trade-off between within-host infection rate and transmission potential in P. subordinaria? We hypothesize that high within-host K E Y W O R D S cross-kingdom interactions, epidemiology, evolution, host-pathogen interactions, life-history traits, parasites

Disease ecology; Landscape ecology; Life history ecology; Microbial ecology; Population ecology; Spatial ecology infection rate is costly in transmission potential. (4) Does coinfection with PlLV alter P. subordinaria within-host infection rate and transmission potential? We hypothesize that coinfection increases within host infection rate and transmission potential.

Plantago lanceolata is a perennial wind pollinated rosette-forming herb (Sagar & Harper, 1964) with a worldwide distribution. In the Åland Islands, it occurs as a network consisting of ca. 4000 meadows in a highly fragmented landscape (Ojanen et al., 2013) . Annually in September, the 4000 P. lanceolata populations' size, area, and location are monitored and the presence of Podosphaera plantaginis fungus and Melitaea cinxia butterfly are surveyed. (Jousimo et al., 2014; Ojanen et al., 2013) . Due to the fragmented landscape in the Åland Islands, P. lanceolata populations are discrete patches surrounded by unsuitable matrix consisting of rocky outcrops, dense forests, waterbodies, or agricultural fields where the plant cannot grow (Ojanen et al., 2013) . The patches are separated by minimum distance of 20 meters of non-suitable habitat or 50 meters of suitable area without P. lanceolata (Ojanen et al., 2013) . Each patch is GPS delineated and their average size is 0.5 hectares (Ojanen et al., 2013) . In the Åland Islands P. lanceolata has been discovered to host two fungal pathogens, P. plantaginis and Phomopsis subordinaria (Jousimo et al., 2014; Laine, 2003) and five viruses (Susi et al., 2019) thus far. Phomopsis subordinaria (telemorph Diaporthe adunca (Rob.) Niessl.) is a specialist fungal pathogen of P. lanceolata (Laine, 2003) transmitted by weevil Trichosirocalus troglodytes (Nieminen & Vikberg, 2015; de Nooji & van der Aa, 1987) . The pathogen infects its host plant through a wound under the inflorescence of the plant causing the developing seeds to dry out (de Nooji & van der Aa, 1987) . The pathogen kills the plant cells and feeds on dead tissue, subsequently causing death of the whole plant. It produces pycnidia fruiting bodies in which its spores are formed. In our laboratory trial, we use pycnidia density as a measure of transmission potential of the pathogen, more abundant pycnidia formation can be expected to also increase spore load carried by the weevil vector. Plantago lanceolata latent virus (PlLV) is a DNA virus belonging to Capulaviruses in Geminiviridae (Susi et al., 2017) . The virus has been recently characterized (Susi et al., 2017 (Susi et al., , 2019 and it is relatively common in populations of P. lanceolata in the Åland Islands (33% of populations infected) (Susi & Laine, 2021) . Mode of transmission and potential host range of PlLV are currently unknown.

To characterize the distribution and drivers of P. subordinaria, and to measure the relationship between within-host and betweenhost dynamics, we surveyed 260 P. lanceolata populations in early September 2018 for infection by P. subordinaria. These populations were selected to represent different areas in the Åland Islands. The field identification of P. subordinaria was confirmed by microscopy of field collected samples consisting of an infected flower stalk. To understand possible correlations between within-host and betweenhost dynamics, we examined two measures in the field populations.

First, we used proportion of flower stalks infected within a plant as a measure of within-host dynamics because it directly quantifies the harm caused to a susceptible host as the death of the seeds.

Phomopsis subordinaria infection spreads from an infected flower stalk through rosette and eventually kills its host plant within 7-10 weeks (de Nooij & Damme, 1988a) . Highly virulent infection will thus prevent the host from producing any seeds whereas lower virulence of the pathogen would allow the plant to produce greater proportion of healthy stalks and seed. Secondly, we meas- (Chiang & Bock, 2021) . Host population size is expected to increase infection risk of populations (Parratt et al., 2016) and thus, we estimated the host population size as coverage of the host plant in square meters. Host population connectivity was measured in order to understand how distances between host populations may impact pathogen prevalence (Jousimo et al., 2014) . Host population connectivity was calculated as where d ij is the Euclidian distance between patches j and i, and α is the parameter of the negative exponential dispersal kernel, which was set to 1 km −1 (see Jousimo et al. (2014) for more details). A j, is the square root transformation of area (m 2 ) of habitat patch j.

The Åland archipelago is highly fragmented and consists of the main island and smaller islands that jointly form 16 regional districts of similar size. To capture possible spatial variation among the regions within the islands, we used regional district as a spatial unit.

In our sampling area, there were nine regional districts (Eckerö (46), with the number of visited populations within each regional district given in parenthesis).

To investigate potential trade-offs between within-host infection rate and transmission potential, and whether coinfection with To measure transmission potential, we observed pycnidia density.

Pycnidia are the fruiting bodies of the fungus that contain the conidial spores that spread the fungus vectored by the weevil T. troglodytes (de Nooji & van der Aa, 1987) . Pycnidia density measures the capacity of the fungus to produce transmission propagules to be transmitted from plant to plant by the weevil vector, and hence, we consider it as a reasonable proxy for transmission potential. At the end of the experiment, 4 weeks past P. subordinaria inoculation,

we counted the pycnidia density within a square centimeter on the flower stalk that produced the first pycnidia from each plant using a microscope.

To analyze drivers of P. subordinaria epidemics in P. lanceolata populations in the Åland Islands, we ran Generalized linear models in SAS Proc Glimmix software (SAS Institute Inc.). We first fit a model with infection occurrence in the host population (0 = no infection, 1 = infection) as a binary response variable, and host population size and connectivity as covariates, and regional district as categorical explanatory variable. Including host population connectivity allows controlling for spatial variation in putative gene flow among populations (Hanski, 1999) , and including regional district allows controlling for possible regional variation in these data. To under- included the within population prevalence as average proportion of infected stalks in a population as explanatory variable, and, host population connectivity, host population size, and regional district as covariates. A Gamma distribution of errors was assumed. We did not include interactions between covariates in the models as they were not biologically meaningful.

We then analyzed the results from the laboratory trial measuring the performance of two P. subordinaria strains alone and in coinfec- 

We found that P. subordinaria was widely spread in P. lanceolata populations in the Åland Islands as in nearly half of the (124) of the sampled 260 populations; one or more infected plants were found ( Figure 1a ). Our first model testing infection occurrence across the populations showed that large populations were more commonly infected than small, and host population connectivity was positively associated with P. subordinaria infection (Table 1 , Figure 1b ,c). We also found regional spatial variation in infection prevalence that was the highest in Hammarland, where 25 of the 37 surveyed host populations were infected, and lowest in Saltvik where six of the 19 surveyed host populations were infected (Table 1, Figure 1a ).

Within infected host populations, P. subordinaria population size was typically low, but the proportion of infected flower stalks was high (Figure 1d ). In nearly half of the populations, we found one to 10 infected plants (Figure 1e ). In one-third of the infected populations, 11-100 infected plants were found, and in only small fraction of populations we found more than 1000 infected individuals ( Figure 1e ). In most infected host populations less than 10% of all plants were infected, and in only small fraction of the populations we observed more than 25% of infected plants (Figure 1f) . At the plant level, we found that P. subordinaria infections were highly virulent. Our second model testing which factors affect P. subordinaria population size, regional districts differed in the pathogen population size, and the proportion of infected stalks had small but significant positive effect on pathogen population size (Table 1; Figure 1g ).

Neither host population connectivity nor size influenced P. subordinaria population size (Table 1 ). None of the tested variables-regional district, host population connectivity and size, P. subordinaria population size, or within-population prevalence-had a significant effect on proportion of infected flower stalks (within-host infection rate)

tested in our third model (Table 1) .

We investigated the effect of PlLV coinfection on P. subordinaria within-host infection rate (measured as lesion growth), and transmission potential (measured as pycnidia density) in an inoculation experiment with two P. subordinaria strains (P29 and P43). The two strains differed significantly in their lesion growth with strain P43 outperforming strain P29 ( Figure 2c ).

Finally, we tested whether there are trade-offs between lesion growth and pycnidia density, and whether P. subordinaria strain or PlLV coinfection influence trade-offs. In this model, pathogen strain had significant impact on pycnidia density, but PlLV coinfection did not (Table1).

We found that while pycnidia density was not directly correlated with lesion growth, the relationship between pycnidia density and lesion growth was mediated by strain identity (significant interaction lesion Figure 3a ) and by coinfection with PlLV (significant interaction lesion growth ×PlLV coinfection; Table 1 ; Figure 3b ).

There was a negative correlation between lesion growth and pycnidia density in strain P29 suggesting a trade-off, whereas in strain P43

there was no evidence of a trade-off (Figure 3a) . Similarly, in mock inoculated plants trade-offs between pycnidia density and lesion growth was observed, but not in PlLV infected plants (Figure 3b ).

Although the trade-off between within-host infection rate and transmission is a central tenet of pathogen evolution (Alizon et al., 2009 Phomopsis subordinaria was detected in nearly half of the surveyed natural host populations. The other fungal pathogen studied in this same host population network-the powdery mildew fungus P. plantaginis-infects annually 2-20% of P. lanceolata populations (Jousimo et al., 2014) , and hence, by comparison P. subordinaria is relatively common. In few other wild plant pathosystems, similar disease incidence rates have been observed (e.g., Triphragmium ulmariae rust infected 29-69% Filipendula ulmaria host populations (Zhan et al., 2018) and Uromyces valerianae rust infected 43-73%

of Valeriana salina populations (Ericson et al., 1999) ). In our study, spatial structure was the main determinant for pathogen incidence as both host population connectivity and regional district explained variation in P. subordinaria distribution. Positive correlation between host connectivity and infection occurrence suggests that host population connectivity increases between population transmission, as predicted by metapopulation theory (Hanski, 1999) . Large populations were more likely to be infected than small populations, which is also in line with theoretical predictions (Hanski, 1999) . Unlike in P. plantaginis (Jousimo et al., 2014) , host population size did not have an effect on pathogen population size. While airborne pathogens such as P. plantaginis are expected to be sensitive to host population size through density-dependent transmission, P. subordinaria is vector-transmitted and hence, it is expected to be less responsive to variation in host population size (Thrall et al., 1995) . This result highlights the importance of the mode of transmission for epidemics.

Our field survey revealed that in infected host plants the majority inflorescences were infected (on average 85% of inflorescences), while in most infected host populations less than 10% of all plants were infected. We did not find evidence of a trade-off between within-host infection rate and transmission limiting the spread of P.

subordinaria within its host populations. Instead, we found a positive correlation between within-host rate and between-host dynamics.

The observed low within-host populations prevalence of P. subordinaria could also result from a range of other factors unrelated to lifehistory correlations, including variation in host plant susceptibility or pathogen infectivity (de Nooij & Damme, 1988a) , constrained vector transmission (de Nooij, 1988; Pleydell et al., 2018) , and interactions with other pathogens than PlLV (Susi, Barres, et al., 2015; .

Cross-kingdom coinfections are common (Chen et al., 2020; Lawn et al., 2006) (Telfer et al., 2010; Tollenaere et al., 2016 Tollenaere et al., , 2017 , and they are often suggested to have serious consequences in disease epidemics, disease severity, and pathogen evolution.

Cross-kingdom coinfections may drive virulence evolution through virulence transmission trade-offs (Alizon et al., 2009 correlations to be depending on the pathogen genotype (Bruns et al., 2014; Clement et al., 2012) . The negative relationship between lesion growth and pycnidia density became positive under coinfection. This is in line with an earlier study testing coinfection with two strains of powdery mildew fungus P. plantaginis where the strains had higher performance and positive life-history trait correlations when challenged with a competing strain (Laine & Mäkinen, 2018) .

Contrary to other studies on coinfection where response to coinfection has been found pathogen strain specific, here the response of the P. subordinaria strains was similar. In bacteria-fluke coinfections on salmons (Louhi et al., 2015) , genotype specific responses on coinfection were observed. Similarly, virus strain combinations resulted in different within-host growth rates of bacterial pathogen on rice (Tollenaere et al., 2017) .

This study increases our understanding on the factors generating diversity in epidemics in natural populations. Furthermore, this is one of the very first reports addressing the knowledge gap on how pathogen life-history traits correlate in realized epidemics. By

showing that pathogen coinfection and strain identity may alter life history correlations this study contributes to better understanding of disease evolution and epidemiology.

We thank Elise Schuurman for her help in performing the experiment and field survey. We are thankful for Alma Oksanen, Heini Karvinen 

The authors have no conflicts of interest to declare. Funding acquisition (lead); Supervision (supporting).

The data underlying the results presented in the study are available from Dryad (https://doi.org/10.5061/dryad.4mw6m 90c9).

https://orcid.org/0000-0001-6814-9205

Suvi Sallinen https://orcid.org/0000-0001-7125-7186

Anna-Liisa Laine https://orcid.org/0000-0002-0703-5850

Virulence-driven trade-offs in disease transmission: A meta-analysis

Plant-virus interactions and the agroecological interface

Multiple infections and the evolution of virulence

Virulence evolution and the trade-off hypothesis: History, current state of affairs and the future

Multiple infections, immune dynamics, and the evolution of virulence

Within-host competition in genetically diverse malaria infections: Parasite virulence and competitive success

Geometagenomics illuminates the impact of agriculture on the distribution and prevalence of plant viruses at the ecosystem scale

A transmission-virulence evolutionary trade-off explains attenuation of HIV-1 in Uganda. eLife, 5, e20492

THE jack of all trades is master of none: A pathogen's ability to infect a greater number of host genotypes comes at a cost of delayed reproduction

The microbial coinfection in COVID-19

Understanding the ramifications of quantitative ordinal scales on accuracy of estimates of disease severity and data analysis in plant pathology

To be or not to be solitary: Phytophthora infestans' dilemma for optimizing its reproductive fitness in multiple infections

Detecting parasite associations within multi-species host and parasite communities

The role of weevils in the infection process of the fungus Phomopsis subordinaria in Plantago lanceolata

Variation in host susceptibility among and within populations of Plantago lanceolata L. infected by the fungus Phomopsis subordinaria

Variation in pathogenicity among and within populations of the fungus Phomopsis subordinaria infecting Plantago lanceolata

Phomopsis subordinaria and related stalk disease in natural populations of Plantago lanceolata L

Virulence and competitive ability in genetically diverse malaria infections

Virulence-transmission trade-offs and population divergence in virulence in a naturally occuring butterfly parasite

Maintenance of variation in virulence and reproduction in populations of an agricultural plant pathogen

Spatial and temporal dynamics of epidemics of the rust fungus Uromyces valerianae on populations of its host Valeriana salina

Models of plant pathogen coevolution

Metapopulation ecology

Small RNAs -Big players in plant-microbe interactions

Ecological and evolutionary effects of fragmentation on infectious disease dynamics

Experimental viral evolution to specific host MHC genotypes reveals fitness and virulence trade-offs in alternative MHC types

First report of Phomopsis subordinaria in a natural population of Plantago lanceolata in south-west Finland

Life-history correlations change under coinfection leading to higher pathogen load

Impact of HIV infection on the epidemiology of tuberculosis in a peri-urban community in South Africa: The need for age-specific interventions

Interactions among bacterial strains and fluke genotypes shape virulence of co-infection

Genetic structure and local adaptation of European wheat yellow rust populations: The role of temperaturespecific adaptation

A virus or more in (nearly) every cell: Ubiquitous networks of virus-host interactions in extreme environments

The insect community of Plantago lanceolata spikes in the Åland Islands

Long-term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia): Survey methods, data management, and long-term population trends

Reduction of growth and reproduction of the biotrophic fungus Blumeria graminis in the presence of a necrotrophic pathogen

Evolution of pathogen virulence across space during an epidemic

Infectious disease dynamics in heterogeneous landscapes

Estimation of the dispersal distances of an aphid-borne virus in a patchy landscape

Lifestyles of plant viruses

Plant virus metagenomics: Advances in virus discovery

The evolution of virulence and pathogenicity in plant pathogen populations

Plantago major L., P. media L. and P. lanceolata L

The area under the disease progress stairs: Calculation, advantage, and application. Analytical and Theoretical Plant Pathology

Fashionably late partners have more fruitful encounters: Impact of the timing of co-infection and pathogenicity on sexual reproduction in Zymoseptoria tritici

The concepts of plant pathogenicity, virulence/avirulence and effector proteins by a teacher of plant pathology

Co-infection alters population dynamics of infectious disease

Diverse and variable virus communities in wild plant populations revealed by metagenomic tools

Pathogen life-history trade-offs revealed in allopatry

Agricultural land use disrupts biodiversity mediation of disease

Genome sequences of a capulavirus infecting Plantago lanceolata in the land archipelago of Finland

Host genotype and coinfection modify the relationship of within and between host transmission

Variation in infectivity and aggressiveness in space and time in wild host-pathogen systems: Causes and consequences

Species interactions in a parasite community drive infection risk in a wildlife population

Frequency-dependent disease transmission and the dynamics of the silene-ustilago hostpathogen system

Virus-bacteria rice co-infection in Africa: Field estimation, reciprocal effects, molecular mechanisms, and evolutionary implications

Evolutionary and epidemiological implications of multiple infections in plants

Do fungi have an innate immune response? An NLR-based comparison to plant and animal immune systems

Plant virus biodiversity and ecology

Climate change accelerates local disease extinction rates in a long-term wild host-pathogen association

Playing on a pathogen's weakness: Using evolution to guide sustainable plant disease control strategies

Coinfection with a virus constrains within-host infection load but increases transmission potential of a highly virulent fungal plant pathogen