# How did hosts domestication modify life history traits of its pathogens?

Marie de Gracia, Mathilde Cascales, Pascale Expert, Marie-Noëlle Bellanger,

Bruno Le Cam, Christophe Lemaire

► To cite this version:

Marie de Gracia, Mathilde Cascales, Pascale Expert, Marie-Noëlle Bellanger, Bruno Le Cam, et al. How did hosts domestication modify life history traits of its pathogens?. PLoS ONE, 2015, 10 (6), pp.1-12. 10.1371/journal.pone.0122909. [hal-01392699](https://hal.science/hal-01392699)

HAL Id: hal-01392699

https://hal.science/hal-01392699

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# How Did Host Domestication Modify Life History Traits of Its Pathogens?

Marie De Gracia¹,²,³, Mathilde Cascales¹,²,³, Pascale Expert¹,²,³, Marie-Noelle Bellanger¹,²,³, Bruno Le Cam¹,²,³‡*, Christophe Lemaire¹,²,³‡*

¹ INRA, IRHS, SFR QUASAV, Beaucouzé, 49071, France, ² Université d'Angers, IRHS, PRES UNAM, SFR QUASAV, Angers, 49045, France, ³ Agrocampus Ouest, IRHS, SFR QUASAV, Angers, 49045, France

‡ These authors are joint senior authors on this work.

* christophe.lemaire@univ-angers.fr (CL); bruno.lecam@angers.inra.fr (BLC)

## Abstract

Understanding evolutionary dynamics of pathogens during domestication of their hosts and rise of agro-ecosystems is essential for durable disease management. Here, we investigated changes in life-history traits of the fungal pathogen *Venturia inaequalis* during domestication of the apple. Life traits linked to fungal dispersal were compared between 60 strains that were sampled in domestic and wild habitats in Kazakhstan, the center of origin of both host and pathogen. Our two main findings are that transition from wild to agro-ecosystems was associated with an increase of both spore size and sporulation capacity; and that distribution of quantitative traits of the domestic population mostly overlapped with those of the wild population. Our results suggest that apple domestication had a considerable impact on fungal characters linked to its dispersal through selection from standing phenotypic diversity. We showed that pestification of *V. inaequalis* in orchards led to an enhanced allocation in colonization ability from standing variation in the wild area. This study emphasizes the potential threat that pathogenic fungal populations living in wild environments represent for durability of resistance in agro-ecosystems.

## Introduction

While the impact of domestication on many organisms has been well described [1-4], little attention has been paid to the evolutionary dynamics of pathogens during domestication of their hosts and rise of agro-ecosystems [5, 6]. Understanding how life history traits of pathogens might have been affected through the crucial process of domestication of their host is particularly relevant. Domestication of plants and animals is one of the most important advances in agriculture in the past 13,000 years. It can be defined as an evolutionary process triggered by the anthropogenic selection of individuals from wild populations based on agronomic criteria [7]. Continuous selection of wild species to new anthropic environments like agro-ecosystems has governed our current diet, the rise of civilizations and the world population demography. Anthropic selection based on some phenotypic traits resulted in a strong accumulation of fixation of alleles, some of them being adaptive, in the genomes of domesticated plants [1, 4, 8],

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** De Gracia M, Cascales M, Expert P, Bellanger M-N, Le Cam B, Lemaire C (2015) How Did Host Domestication Modify Life History Traits of Its Pathogens? PLoS ONE 10(6): e0122909.
doi:10.1371/journal.pone.0122909

**Academic Editor:** Richard A Wilson, University of Nebraska-Lincoln, UNITED STATES

**Received:** October 27, 2014

**Accepted:** February 16, 2015

**Published:** June 19, 2015

**Copyright:** © 2015 De Gracia et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** The PhD thesis of Marie De Gracia was funded by both the Département Santé des Plantes et Environnement de l'Institut National de Recherche Agronomique (INRA) (http://www.spe.inra.fr/) and the Région Pays de Loire (http://www.paysdelaloire.fr/). The project was granted by Agence Nationale Pour la Recherche (ANR): Projet Bioadapt ANR Gandalf ANR-12-ADAP-0009 (http://www.agence-nationale-recherche.fr/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
**Competing Interests:** The authors have declared that no competing interests exist.

animals [9] but also microorganisms [10]. The set of morphological and/or genetic differences between domesticated species and their wild progenitors defines the domestication syndrome [11, 12]. Among pathogens attacking our basic food resources, fungi cause more than a third of emerging diseases [13, 14]. The question here is how much does host domestication change the adaptive landscape for pathogens? Understanding how pathogens living on wild hosts adapt to agricultural constraints would allow us to manage emerging diseases and improve the sustainability of host genetic resistances.

Beside the proper phenotypic evolution of the plant during domestication, the influence of transition from natural to domestic environments on host-pathogen co-evolution is not negligible. Phenotypic values of plants living in agro-ecosystems are expected to be quite different from those measured in the wild. Indeed, natural environments are characterized by a random distribution of host plants, low plant density, high genetic diversity and uneven distribution of nutrients [6]. In contrast, the development of agriculture favored the ability of plants to grow in new habitats with higher densities, but less genetic diversity, increased and uniform nutrient availability and less interspecific competition. The ability of a pathogen to adapt to this new environment therefore depends on its life history and potential for evolution of traits [14, 15]. Recently, the term pestification was introduced for the first time by Saleh et al. [16] to describe the influence of rice domestication on its associated fungal pathogen *Magnaporthe oryzae*. In this study, the authors revealed modifications in the genetic structure of the pathogen's population during its evolutionary history without any investigation on life history traits. Here, we define the term pestification as the process by which a pathogen originally living in the wild become adapted to biotic and abiotic crop conditions. Among traits affected by pestification, two categories are particularly relevant in epidemiology: growth and dispersal [17]. Intrinsic hyphal growth is a classical measurement of fitness [17, 18]. Fungal pathogens are often characterized by a dispersive phase mediated by asexual (mitotic) spores also called conidia. Cateria like the number of conidia produced, their size, and their germination rate are good predictors of fitness [17] and are relevant to determine the ability for dispersal. Determining how growth and dispersal of fungal pathogens have evolved during domestication of their host would help us better understanding their adaptive dynamics.

The pathosystem *Malus spp.-Venturia inaequalis* is particularly relevant to investigate changes in life history traits of pathogens that are adapted to domesticated hosts and agro-ecosystems. Indeed, the life history of *V. inaequalis* confers to this fungus a significant potential for evolution [15]. *V. inaequalis* produces both ascospores, spores issued from sexual reproduction, and conidia. Intrinsic dispersal ability of both types of spores is low, the most favorable way for the pathogen to achieve long distance dispersal being linked to human activities, such as transport of infected fruits and plants [19, 20]. It is assumed that the population structure of a pathogen is linked to the history of the spread of its host [21, 22]. Cultivated apple was domesticated in Central Asia from *Malus sieversii* [23–25]. This species makes large forests from the Tien Shan mountains (border of Kazakhstan, Kirghizstan and Xinjiang Province in China) to the Caspian Sea [26, 27]. Migration westwards have involved contact and hybridization with close wild relatives growing along the Silk Road, including *M. sylvestris* in Europe and *M. orientalis* in Caucasus [28]. These wild apple species have thus also contributed to the genetic mixing of domesticated apple. Previous genetic studies on *V. inaequalis* have provided many clues about its origin. A global analysis of diversity showed that *V. inaequalis* shares a common origin with its host in Central Asia [22]. Impact of apple domestication on genetic structure of the pathogen has then been investigated [29]. Focusing their study in Central Asia, the authors identified a population sampled on *M. sieversii* in natural forests of Kazakhstan, which was considered to be a relic of the ancestral population from which have been derived populations currently present in agro-ecosystems worldwide. More recently, a study focusing on
pathogenicity traits of this ancestral population showed that apple domestication had markedly modified pathogenicity of *V. inaequalis* [30]. Indeed, the domestication process was associated with acquisition of virulence and increase of aggressiveness of *V. inaequalis*. The availability of valuable genetic resources and knowledge on this pathosystem provide an excellent framework to investigate the impact of host domestication on evolution of pathogen life history traits.

Here, we compared several life history traits of *V. inaequalis* isolates that were sampled in the ancestral and domestic populations in Kazakhstan. We chose traits that are known to be involved in dispersal and colonization success of fungal pathogens: mycelium growth, asexual spore germination rate, number of asexual spores and size of asexual spores. We addressed the following questions: (i) Are strains sampled in the domestic habitat phenotypically different from those isolated in the wild habitat? If so, which traits are involved in these differences? (ii) Did domestication of apple trees drive evolution of *V. inaequalis* to a new adaptive peak rather than to a standing phenotypic combination selected for life into agroecosystems? (iii) What are the implications of an eventual phenotypic change on the life history of *V. inaequalis* in agroecosystems?

## Material and Methods

We state that no permit was necessary to sample apple tree leaves containing the fungus *Venturia inaequalis* in both wild and agricultural locations in Kazakhstan. The field study did not involve endangered nor protected species.

### Fungal isolates

This study was based on a total of 60 isolates of *Venturia inaequalis* sampled on the wild apple *M. sieversii* in Kazakhstan. Thirty were collected in the wild forest area (43.23°N, 77.28°E) and 30 in the suburbs or agricultural area around Almaty (43.25°N, 76.91°E), the former capital of Kazakhstan. All these isolates were collected in 2006. Throughout the study, the terms wild population (CAM for Central Asian Mountains) and domestic population (CAP, Central Asian Plains) are used to make reference to a collection of individuals sampled in non-anthropic area and in anthropic area, respectively, as defined in [29].

### Phenotypic traits measurements

**Hyphal growth.** Estimation of the mycelium growth capacity of isolates was performed on malt agar medium. For each isolate, eight plugs of mycelium from a two-week-old culture were grown on a malt agar medium. Cultures were kept at 18°C. Considering that *V. inaequalis* growth is circular, we measured colony diameters using a ruler 14 days after mycelium plugs deposit. Two measurements were performed for each plug, which represents 16 measurements per isolate.

**Spore traits measurements.** Spore traits analysis included measuring spore germination rate, sporulation capacity and spore size. Initial inocula were obtained by growing single spores of each isolate on cellophane sheets deposited onto malt agar medium.

Germination rates were assessed for each isolate on malt agar plates covered with a cellophane sheet. Plates were kept in darkness at 18°C to allow conidia germination. After 24 hours, germinated and non-germinated spores were counted using a microscope. The *in vitro* germination rate was then defined as the percentage of germinated spores after 24 hours. For each isolate, a single germination rate was calculated.

Concerning sporulation capacity and spore size measurements, 10 µL of the initial inoculum was calibrated to 105 conidia.mL⁻¹ and deposited on four malt agar Petri dishes covered by cellophane sheets of 3 cm². Petri dishes were then incubated for seven days at 18°C under the
light. After this period, cellophane sheets were water rinsed for measurements of spore size and spore concentration using a spore counter (Counter Multisizer, Beckman). Sixteen measurements were made for each isolate. The counter was calibrated to count only particles with a diameter in the range from 6.5 to 12.5 µm (range size of *Venturia inaequalis* conidia).

## Statistical analyses of the phenotypic data

All statistical analyses were performed using R Version 3.0.0 [31]. In order to test for independence between life history traits and for any trade-off change within each population, we performed pairwise Pearson's correlation tests. Pearson's correlation coefficients were estimated using the "corrgram" package [32].

For each life history trait, distributions were checked for normality using the Kolmogorov-Smirnov test. To test for differences between populations in each phenotypic variable, we used different procedures according to the availability of replicates. For differences between populations in sporulation capacity and spore size, an analysis of variance (ANOVA) was performed with a linear mixed-effect model (LME) implemented in the "nlme" package [33]. The "individual" factor was treated as random effect and nested in "population". Germination rate for each isolate were analyzed using Student's tests. Percentages were transformed by a logit transformation. Population effect on variation in hyphal growth was analyzed using an ANOVA. In all ANOVAs, normality of residuals was checked.

## Results

### Quantitative distribution of life history trait variables

The distribution of each trait measurement is presented in Fig 1. Analyses of hyphal growth were based on diameter measurements 24 days after mycelium plugs were deposited. The overall hyphal growth was between 1.29 and 3.53 cm in diameter. Data ranged from 1.28 to 2.94 cm (mean ± SD: 2.22 ± 0.41) in the wild population and from 1.39 to 3.53 cm (mean ± SD: 2.14± 0.50) in the domestic population. Spore analyses compared spore germination rates, sporulation capacity and spore size.

For the overall sample, germination rate varied between 45 and 96%: this rate ranged from 54 to 96% (mean ± SD: 81 ± 10.82) in the wild population, and from 45 to 96% (mean ± SD: 79.6 ± 14.08) in the domestic one. Sporulation capacity fluctuated much more than the other characters. The lowest sporulation capacity (3,250 spores mL⁻¹) was observed in a strain from the wild population and the highest in a strain originating from the domestic population (534,400 spores mL⁻¹).

Distribution for each population ranged from 3,200 to 298,300 spores mL⁻¹ (mean ± SD: 122,400 ± 80,000) and from 39,600 to 534,400 spores mL⁻¹ (mean ± SD: 200,000 ± 108,100) for wild and domestic populations, respectively. Concerning the spore size, isolates from the wild population ranged from 8.12 to 10.15 µm (mean ± SD: 9.26 ± 0.44), while those from domestic isolates ranged from 8.45 to 10.45 µm (mean ± SD: 9.55 ± 0.36).

### Correlation between life history traits

Pearson's correlations were tested between life history traits and pairwise correlation analyses were performed. Spore size was significantly correlated with sporulation capacity in each population (Fig 2). No other correlation with hyphal growth or germination rate was detected. The germination rate of conidia was not correlated with any other trait, meaning that the ability to germinate does not determine the size of the spores, nor the sporulation capacity.
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Quantitative distribution of life history trait variables. Distribution of quantitative variables: hyphal growth 24 days after deposits on agar (cm), germination rate (in logit units), sporulation capacity (in Log₁₀ units), and spore size (µm). Blue histograms correspond to the wild population and green ones correspond to the domestic population.</figcaption></figure>

doi:10.1371/journal.pone.0122909.g001

Spore size was significantly correlated with sporulation capacity, with Pearson's correlation coefficients of *R* = 0.64 in the wild population (*P*<0.001) and *R* = 0.48 in the domestic pone (*P*<0.01) (Fig 2).

## Life traits differences between wild and domestic populations

The germination rate of the wild population was not significantly different from the rate of the domestic one (81% and 79.6%, respectively) (Student test; *t* = 0.29; *P* = 0.77) (Fig 3).
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

**Fig 2. Pairwise correlations between phenotypic traits.** Scatter plots indicating the degree of correlation between traits (below diagonal) and pairwise Pearson's correlation R values (above diagonal). Pearson's correlation coefficients were calculated from averaged data from each strain for sporulation capacity, spore size and hyphal growth. Statistical significant correlation coefficients are indicated as following: ** P-value < 0.01 and *** P-value < 0.001. Calculations were performed for domestic (A panel) and wild (B panel) populations, respectively.

doi:10.1371/journal.pone.0122909.g002
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Quantitative distribution of life history trait variables between wild (blue) and domestic population (green). The horizontal line in boxplots represents the median. A significant difference between populations is observed for the sporulation capacity and spore size. Sporulation average in wild population: 122,400 spores mL⁻¹; and in domestic population: 199,955 spores mL⁻¹. Spore size average is 9.2 µm in the wild population and is 9.5 µm in the domestic population.</figcaption></figure>

doi:10.1371/journal.pone.0122909.g003

Hyphal growth of the wild population was not significantly different from that one of the domestic population. Twenty-four days after deposition, average hyphal growth was 2.21 and 2.13 for the wild and domestic populations, respectively (ANOVA, *F* = 0.40; *P* = 0.52) (Fig 3).

A significant difference in sporulation capacity between the two populations was detected (ANOVA; *F* = 188; *P*<0.0001). The domestic population exhibited a sporulation rate higher than that of the wild population. The average rates in the domestic and wild populations were 200,000 spores mL⁻¹ and 122,400 spores mL⁻¹, respectively (Fig 3). The spore size of isolates from the wild population was significantly smaller than the spore size of isolates from the domestic population (ANOVA; *F* = 19085; *P*<0.0001). Average size of the domestic population was 9.5 µm and that one of the wild population was 9.2 µm (Fig 3). For these two phenotypic analyses, we noted that distributions of both populations largely overlapped, with nearly all values of the wild isolates encompassed by the distribution of values of the domestic isolates. This result suggests that isolates sampled in the anthropic area were actually a sampling from the largest and most sporulating isolates of the wild population (Fig 3).

Our results demonstrate that apple domestication promoted modifications in traits related to dispersal of the apple scab fungus, *Venturia inaequalis*. Among the four life history traits that were analyzed, both spore production and size showed significant difference between the two tested sampling sites. Our results show that diversity in spore size and production is significantly higher in wild habitats than in agricultural ones.
## Discussion

Impact of plant domestication on pathogenic organisms has rarely been investigated despite its major importance in understanding consequences of modern human activities on disease emergence. Previous identification of fungal genetic resources that did not experience agriculture among all other human activities allowed us addressing this issue of great interest in evolutionary ecology. We used several phenotypic comparisons to investigate changes of different life history traits involved in dispersal of the apple scab fungus, associated with the domestication of its host. We chose strains that had been sampled on *M. sieversii*, the main wild Asian ancestor of the domestic apple tree, located in a wild forest and in the suburbs of Almaty. As previously mentioned, strains from the latter location, although sampled on *M. sieversii* trees, have been shown to belong to a population that is also present on domestic apple trees [29, 30]. Those strains were probably imported with domestic apple in Kazakhstan in the late XIXth century. Indeed, because Kazakhstan was primarily inhabited by nomads, agriculture was not developed in this area before the annexation by the Russian Empire and colonization by Russian settlers. Selecting strains living on the same host species prevents from confounding effects due to adaptation to the host genotype. Our two main findings are that (1) transition from wild to agro-ecosystem was associated with an increase of the average of both spore size and sporulation capacity; and (2) that distribution of each quantitative trait of the domestic population mostly overlapped with those of the wild population. Our results suggest that transition from wild to orchard had a considerable impact on fungal characters linked to its potential for dispersal by selection from standing phenotypic diversity. In the present study, we showed that pestification of *V. inaequalis* in orchards led to an enhanced allocation in colonization ability for which potential capacity did already exist in the wild area.

## Correlation between fitness traits

The analysis of fitness-related traits is complex because fitness is a synthetic function of both survival and reproductive success. Disentangling between these two fitness components implies the measurement of several traits and the analysis of their pairwise correlations.

Our results clearly showed that spore production and spore size on the one hand were not correlated with mycelial growth on the other hand (Fig 2). A likely explanation would be that growth of *Venturia inaequalis* is independent of its ability to produce asexual spores. In other words, the quantity of mycelium produced per time unit does not impact the net production of conidia. This result is at odds with observations made on *Aspergillus niger* for which the mycelium area was significantly correlated with spore number [34]. Mycelium of *Venturia inaequalis* grows under the cuticle of the leaf and is therefore not visible. Aggressiveness is usually determined by the size of lesions harboring conidiospores on the leaf surface. Our results suggest that this way of measuring aggressiveness might not reflect the growth of *Venturia inaequalis* within the leaf tissues. Nevertheless, our results were obtained *in vitro* and need to be validated *in planta*, for instance using thermography or fluorescence-based methods [35]. In addition in our study, hyphal growth was analyzed from measurements of the diameter of the mycelium surface of each isolate. Further analyses would be useful because a mycelium is a network formed by a succession of hyphae with similar branching patterns at different scales that fit fractal geometry [36, 37]. A description and modeling of the fractal geometry of mycelium may be more appropriate to evaluate mycelium growth. Spore size is known to be a trait related to dispersal abilities of fungi. Indeed, smaller fungal spores are currently considered as fitter than large spores because they are more effectively dispersed [17]. Modeling the dynamics of airborne propagules clearly showed that spore size is negatively correlated with dispersal ability [38, 39]. However, in their recent theoretical study that looked at particle size ranging from 1
to 10 µm (the actual size of conidia in our study), Norros et al. [40] claimed that "*the dispersiveness of especially the smallest spores is so high that landing may in fact be a bigger challenge for them than flying*". Larger spores are more likely to be deposited on hosts [39–41]. Thus, even if smaller spores disperse farther than larger spores, the efficiency of dispersal in terms of probability of landing on a host would favor larger spores.

Interestingly, both spore production and size are strongly positively correlated in wild and agricultural habitats. Unlike plants and animals, the more spores produced, the bigger their size. This strong correlation probably reflects a linkage, even pleiotropy, between genes controlling spore production and spore size. Taking into account that larger spores maximize their ability to be deposited on canopy [40, 42], our results suggest that strains that produce more asexual spores have a better efficiency of dispersal. Why did selection not favor strains with higher spore production and larger spores? The variability for both spore size and production in our dataset probably account for some sort of trade-off as previously mentioned by Meerts [41] on Basidomyota. Optimal spore size across species may vary according to life history. Thus, as *Venturia inaequalis* disperses asexual spores on leaves and fruits, mostly by splashing events over short distances [43–45], one might expect that larger spores should be favored in this species. We will discuss this point further when comparing strains from the wild and agricultural habitats.

Surprisingly, the germination rate of conidia did not correlate with any other trait. Larger spores were expected to have greater germination rates [17], but to our knowledge there is no clear demonstration of this correlation in fungi. Our results showed that spore size does not determine the ability to germinate. However, there is no data on survival rate of germinated spores. Indeed, observations showed that mortality could occur from several hours to days after germination. Such a correlation study between spore size and survival should deserve attention in the future.

## Traits related to the efficiency of pathogen's dispersal were impacted during wild-to-crop transition

Our results demonstrated that apple domestication promoted modifications in traits related to dispersal of the apple scab fungus *Venturia inaequalis*. Among the four life history traits that were analyzed, both spore production and spore size showed significant differences between the two tested sampling sites. This means that agro-ecosystems might have selected a high production of larger spores for a more efficient short-distance dispersal. Indeed, wild apple trees, *M. sieversii*, in the Tien-Shan mountains of Kazakhstan make dense and large forests extending along valleys but also form patches in meadows at the peripheral of forests. Wild apple trees exhibit high phenotypic diversity and thus represent highly heterogeneous habitats for *Venturia inaequalis*. In contrast, *M. sieversii* trees around Almaty grow in a landscape dominated by agriculture, especially apple orchards that are characterized by a high homogeneity of genotypes and density. Higher spore size probably increases the probability to be deposited on neighbor apple trees in orchards. If a strain produces spores that are already adapted to the genotype of its host, likely this strain is also adapted to the hosts all around. Thus, selecting for high sporulation rate would also select for large spores, which potentially gives an advantage in competition in a homogeneous habitat. Conversely in wild heterogeneous habitats, a strain adapted to its host might not be adapted to the neighboring hosts. The potential for dispersal of relatively smaller spores might ensure reaching another canopy or distant compatible habitats [40, 41].

## Pestification by selection from standing phenotypic variation in the wild

The domestication process of plants and animals has often been associated to a strong disruptive selection favoring certain traits in each environment. In plants, transition between wild
and cultivated forms was generally associated with different selection-targeted traits related to harvesting conditions, seed production or seedling competition [46].

We showed that phenotypic traits in agro-ecosystems were selected from standing variation in the wild. Indeed, rather than a strong disruptive selection, we observed that the range of each measured trait in the anthropized population mostly overlapped with the range of those traits in the wild population. We showed that phenotypic diversity within the domestic population preexisted in the wild population at least for the traits we have analyzed (Fig 3). The most trivial reason is that humans have not selected the pathogen for any agricultural or biotechnical trait of interest compared to plants or even domesticated fungi such as *Penicillium* [47] or *Saccharomyces* species [48]. As host tracking is a continuous process, we might suppose that plant modifications were progressive and therefore did not exert a strong selection on its pathogen. In the same line, it is noteworthy that *Malus x domestica* remains genetically very close to its wild relative species *M. sieversii*, *M. sylvestris* and *M. orientalis* [24]. We previously showed that moving this fungus outside its native range had not induced shifts in its reproductive mode, which is temperature dependent [22]. We also showed that limited gene flows occur between domestic and wild populations [29]. At this time, we cannot predict whether such low extent of genetic exchanges will maintain the observed phenotypic differences.

We presented here one of the first study of evolution of fitness-related traits in a pathogenic fungus triggered by the domestication of its host. As previously introduced by Saleh *et al.* [16], we called this process pestification. We found that this process led to an increase of the ability for efficient colonization of a homogeneous environment. In *Venturia inaequalis*, pestification mostly proceeded by selection from standing phenotypic variation in the wild. However, we do not know to which extent our conclusions can be transferred to other pathogens because evolution of fitness-related traits is highly dependent of the different life histories like predominance of sexual over asexual reproduction [17]. Our study also focused on conidiospores that ensure clonal dispersal. We have no information about phenotypic variation on ascospores. Dispersal of ascospores by air and rain splashes is crucial for the life cycle of *Venturia inaequalis* because it permits the very first colonization of buds and young leaves [43, 49, 50]. Nevertheless we found a high phenotypic diversity in the wild strains for the four traits we have analyzed.

We cannot exclude such a diversity for other traits like virulence or aggressiveness as suggested by Lê Van *et al.* [30]. This study emphasizes the potential threat for the durability of crop resistances that is represented by populations of pathogens living in wild environments [30, 51, 52].

## Acknowledgments

We are grateful to all the people who helped us collecting strains: AD Dzhangaliev and collaborators, Catherine Peix and François Laurens. We would like to thank V. Caffier for all advices concerning experimental procedures, and Jérôme Collemare for all criticisms. We also thank the French group "REID Champignons" for helpful comments.

## Author Contributions

Conceived and designed the experiments: MDG CL BLC PE MNB. Performed the experiments: MC MDG PE MNB. Analyzed the data: MC MDG CL. Contributed reagents/materials/analysis tools: BLC CL PE MNB. Wrote the paper: MDG BLC CL.

## References

1. Tenaillon MI, U'Ren J, Tenaillon O, Gaut BS. Selection versus demography: A multilocus investigation of the domestication process in maize. Mol Biol and Evol. 2004; 21: 1214–1225. PMID: 15014173
2. Liu A, Burke JM. Patterns of nucleotide diversity in wild and cultivated sunflower. Genetics. 2006; 173: 321-330. PMID: 16322511

3. Haudry A, Cenci A, Ravel C, Bataillon T, Brunel D, Poncet C, et al. Grinding up wheat: A massive loss of nucleotide diversity since domestication. Mol Biol Evol. 2007; 24: 1506-1517. PMID: 17443011

4. Glémin S, Bataillon T. A comparative view of the evolution of grasses under domestication: Tansley review. New Phytol. 2009; 183: 273-290. doi: 10.1111/j.1469-8137.2009.02884.x PMID: 19515223

5. Buckler ES, Thornsberry JM, Kresovich S. Molecular diversity, structure and domestication of grasses. Gen Res. 2001; 77: 213-218.

6. Stukenbrock EH, McDonald BA. The origins of plant pathogens in agroecosystems. Ann Rev Phytopathol. 2008; 46: 75-100. doi: 10.1146/annurev.phyto.010708.154114 PMID: 18680424

7. Clement CR. 1492 and the loss of Amazonian crop genetic resources. I. The relation between domestication and human population decline. Econ Bot; 1999; 53: 188-202.

8. Hyten DL, Song Q, Zhu Y, Choi IY, Nelson RL, Costa JM, et al. Impacts of genetic bottlenecks on soybean genome diversity. P Natl Acad Sci USA. 2006; 103: 16666-16671. PMID: 17068128

9. Rubin CJ, Zody MC, Eriksson J, Meadows JRS, Sherwood E, Webster MT, et al. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010; 464: 587-591. doi: 10.1038/nature08832 PMID: 20220755

10. Liti G, Carter DM, Moses AM, Warringer J, Parts L, James SA, et al. Population genomics of domestic and wild yeasts. Nature. 2009; 458: 337-341. doi: 10.1038/nature07743 PMID: 19212322

11. Hammer K. Das Domestikations syndrom. Kulturpflanze. 1984; 32: 11-34.

12. Li CB, Zhou AL, Sang T. Genetic analysis of rice domestication syndrome with the wild annual species, *Oryza nivara*. New Phytol. 2006; 170: 185-193. PMID: 16539615

13. Anderson PK, Cunningham AA, Patel NG, Morales FJ, Epstein PR, Daszak P. Emerging infectious diseases of plants: pathogen pollution, climate change and agrotechnology drivers. Trends Ecol Evol. 2004; 19: 535-544. PMID: 16701319

14. Giraud T, Gladieux P, Gavrilets S. Linking the emergence of fungal plant diseases with ecological speciation. Trends Ecol Evol. 2010; 25: 387-395. doi: 10.1016/j.tree.2010.03.006 PMID: 20434790

15. McDonald BA, Linde C. Pathogen population genetics, evolutionary potential, and durable resistance. Ann Rev Phytopathol. 2002; 40: 349-379. PMID: 12147764

16. Saleh D, Milazzo J, Adreit H, Fournier E, Tharreau D. South-East Asia is the center of origin, diversity and dispersion of the rice blast fungus, *Magnaporthe oryzae*. New Phytol. 2014; 201: 1440-1456. doi: 10.1111/nph.12627 PMID: 24320224

17. Pringle A, Taylor JW. The fitness of filamentous fungi. Trends Microbiol. 2002; 10: 474-481. PMID: 12377558

18. Brasier CM. Fitness, continuous variation and selection in fungal populations: an ecological perspective. Gen Res. 1999; 77: 213-218.

19. McHardy WE. Apple Scab: Biology, Epidemiology and Management. The American Phytopathological Society. 1996

20. Holb I, Heijine B, Withagen J, Jeger M. Dispersal of *Venturia inaequalis* ascospores and disease gradients from a defined Inoculum source. J Phytopathol. 2004; 152: 639-646.

21. Wirth T, Meyer A, Achtman M. Deciphering host migrations and origins by means of their microbes. Mol Ecol. 2005; 14: 3289-3306. PMID: 16156803

22. Gladieux P, Zhang XG, Afoufa-Bastien D, Sanhueza RMV, Sbaghi M, Le Cam B. On the origin and spread of the scab disease of apple: Out of Central Asia. PloS One. 2008; 1: e1455. doi: 10.1371/journal.pone.0001455 PMID: 18197265

23. Velasco R, Zharkikh A, Affourtit J, Dingra A, Cestaro A, Kalyanaraman A, et al. The genome of the domesticated apple (Malus x domestica Borkh.). Nature Genet. 2010; 42: 833-939. doi: 10.1038/ng.654 PMID: 20802477

24. Cornille A, Gladieux P, Smulders MJ, Roldan-Ruiz I, Laurens F, Le Cam B, et al. New insight into the history of domesticated apple: secondary contribution of the European wild apple to the genome of cultivated varieties. PLoS Genet. 2012; 8: e1002703. doi: 10.1371/journal.pgen.1002703 PMID: 22589740

25. Cornille A, Giraud T, Smulders MJ, Roldan-Ruiz I, Gladieux P. The domestication and evolutionary ecology of apples. Trends Genet. 2014; 30: 57-65. doi: 10.1016/j.tig.2013.10.002 PMID: 24290193

26. Morgan J, Richards A. The Book of Apples. The Ebury Press. 1993

27. Harris SA, Robinson JP, Juniper BE. Genetic clues to the origin of the apple. Trends Genet. 2002; 18: 426-430. PMID: 12142012
28. Cornille A, Giraud T, Bellard C, Tellier A, Le Cam B, Smulders MJM, et al. Postglacial recolonization history of the European crabapple (*Malus sylvestris* Mill.), a wild contributor to the domesticated apple. Mol Ecol. 2013; 22: 2249–2263. doi: 10.1111/mec.12231 PMID: 23402276

29. Gladieux P, Zhang XG, Roldan-Ruiz I, Caffier V, Leroy T, Devaux M, et al. Evolution of the population structure of *Venturia inaequalis*, the apple scab fungus, associated with the domestication of its host. Mol Ecol. 2010; 19: 658–674. doi: 10.1111/j.1365-294X.2009.04498.x PMID: 20088887

30. Lê Van A, Gladieux P, Lemaire C, Cornille A, Giraud T, Durel CE, et al. Evolution of pathogenicity traits in the apple scab fungal pathogen in response to the domestication of its host. Evol Appl. 2012; 5: 694–704. doi: 10.1111/j.1752-4571.2012.00246.x PMID: 23144656

31. R-Development-Core-Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing [http://www.Rproject.org]. 2013 doi: 10.3758/s13428-013-0330-5 PMID: 23519455

32. Wright K. Corrgram: Calculates correlation of variables and displays the results graphically. R package version 1.5. 2013.

33. Pinheiro J, Bates D, Debroy S, Sarkar D. Nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1–117. 2014.

34. de Visser JAGM, Hoekstra RF, Van Den Ende V. Test of interaction between genetic markers that affect fitness in *Aspergillus niger*. Evolution. 1997; 51: 1499–1501

35. Belin E, Rousseau D, Boureau T, Caffier V. Thermography versus chlorophyll fluorescence imaging for detection and quantification of apple scab. Comput Electron Agr. 2013; 90: 159–163.

36. Papagianni M. Quantification of the fractal nature of mycelial aggregation in *Aspergillus niger* sub-merged cultures. Microb Cell Fact. 2006; 5: doi: 10.1186/1475-2859-5-5

37. Obert M, Pfeifer P, Sernetz M. Microbial growth patterns described by fractal geometry. J Bacteriol. 1990; 172: 1180–1185. PMID: 2106504

38. Whitehead R. Wind pollination in the Angiosperms: evolutionary and environmental considerations. Evolution. 1969; 23: 28–35.

39. Dix NJ, Webster J. Fungal Ecology. London, UK: Chapman et Hall. 1995.

40. Norros V, Rannik U, Hussein T, Petäjä T, Vesala T, Ovaskainen O. Do small spores disperse further than large spores? Ecology. 2014; 95: 1612–1621. PMID: 25039225

41. Meerts P. The evolution of spores size in Agarics: do big mushrooms have big spores? J Evol Biol. 1999; 12: 161–165.

42. Ingold CT. Fungal spores: Their liberation and dispersal. Clarendon press. 1971.

43. Aylor DE. The aerobiology of apple scab. Plant Dis. 1998; 82: 838–848.

44. McHardy WE, Gadoury DM, Gessler C. Parasitic and biological fitness of *Venturia inaequalis*: Relationship to disease management strategies. Plant Dis. 2001; 85: 1036–1051.

45. Leroy T. Lemaire C, Dunemann F, Le Cam B. The genetic structure of a Venturia inaequalis population in a heterogeneous host population composed of different Malus species. BMC Evol Biol. 2013; 13: doi: 10.1186/1471-2148-13-64

46. Harlan J, de Wet J, Price E. Comparative evolution of cereals. Evolution. 1973; 27: 311–325.

47. Ropars J, López-Villavicencio M, Dupont J, Snirc A, Gillot G, Coton M, et al. Induction of sexual reproduction and genetic diversity in the cheese fungus *Penicillium roqueforti*. Evol Appl. 2014; 7: 433–441. doi: 10.1111/eva.12140 PMID: 24822078

48. Hittinger CT. Saccharomyces diversity and evolution: a budding model genus. Trends Genet. 2013; 29: 309–317. doi: 10.1016/j.tig.2013.01.002 PMID: 23395329

49. Aylor DE, Sutton TB. Release of Venturia inaequalis ascospores during unsteady rain: relationship to spore transport and deposition. Phytopathology. 1992; 82:532–540.

50. Aylor DE, Ducharme KM. Wind fluctuations near the ground during rain. Agr Forest Meteorol. 1995; 76:59–73.

51. Monteil CL, Cai R, Liu H, Mechan Llontop ME, Leman S, Studholme DJ, et al. Nonagricultural reservoirs contribute to emergence and evolution of *Pseudomonas syringae* crop pathogens. New Phytol. 2013; 199: 800–811. doi: 10.1111/nph.12316 PMID: 23692644

52. Leroy T, Le Cam B, Lemaire C. 2014. When virulence originates from non-agricultural hosts: New insights into plant breeding. Infect Genet Evol. 2014; doi: 10.1016/j.meegid.2013.12.022