# Genetic and Pathogenicity Diversity of Aphanomyces euteiches Populations From Pea-Growing Regions in France

Anne Quillévéré-Hamard, Gwenola Le Roy, Anne Moussart, Alain Baranger, Didier Andrivon, Marie-Laure Pilet-Nayel, Christophe Le May

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

## ► To cite this version:

Anne Quillévéré-Hamard, Gwenola Le Roy, Anne Moussart, Alain Baranger, Didier Andrivon, et al. Genetic and Pathogenicity Diversity of Aphanomyces euteiches Populations From Pea-Growing Regions in France. Frontiers in Plant Science, 2018, 9, pp.1673. 10.3389/fpls.2018.01673. hal-01943017

HAL Id: hal-01943017

https://institut-agro-rennes-angers.hal.science/hal-01943017

Submitted on 3 Dec 2018

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_2.png" /></figure>

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

Check for
updates

# Genetic and Pathogenicity Diversity of *Aphanomyces euteiches* Populations From Pea-Growing Regions in France

Anne Quillévéré-Hamard<fnref n="1" />,², Gwenola Le Roy<fnref n="1" />,², Anne Moussart<fnref n="1" />,²,³, Alain Baranger<fnref n="1" />,², Didier Andrivon<fnref n="1" />, Marie-Laure Pilet-Nayel<fnref n="1" />,²†† and Christophe Le May²,⁴*†

OPEN ACCESS

**Edited by:**
Silvia Matesanz,
Universidad Rey Juan Carlos, Spain

**Reviewed by:**
Alfredo García-Fernández,
Universidad Rey Juan Carlos, Spain
Yessica Rico,
Instituto de Ecología (INECOL),
Mexico

*Correspondence:
Christophe Le May
lemay@agrocampus-ouest.fr

†These authors have contributed
equally to this work

‡Marie-Laure Pilet-Nayel
orcid.org/0000-0002-4505-1907

**Specialty section:**
This article was submitted to
Evolutionary and Population Genetics,
a section of the journal
Frontiers in Plant Science

**Received:** 05 April 2018
**Accepted:** 26 October 2018
**Published:** 19 November 2018

**Citation:**
Quillévéré-Hamard A, Le Roy G, Moussart A, Baranger A, Andrivon D, Pilet-Nayel M-L and Le May C (2018) Genetic and Pathogenicity Diversity of Aphanomyces euteiches Populations From Pea-Growing Regions in France. *Front. Plant Sci.* 9:1673.
doi: 10.3389/fpls.2018.01673

*Aphanomyces euteiches* is an oomycete pathogen with a broad host-range on legumes that causes devastating root rot disease in many pea-growing countries and especially in France. Genetic resistance is a promising way to manage the disease since consistent QTL controlling partial resistance have been identified in near isogenic lines of pea. However, there are still no resistant pea varieties cultivated in France. This study aimed to evaluate the phenotypic and genetic diversity of *A. euteiches* populations from the major pea-growing regions in France. A collection of 205 isolates, from soil samples collected in infested pea fields located in five French regions, was established and genotyped using 20 SSR markers. Thirteen multilocus genotypes were found among the 205 isolates which displayed a low genotypic richness (ranged from 0 to 0.333). Two main clusters of isolates were identified using PCoA and STRUCTURE, including a predominant group comprising 88% of isolates and another group representing 12% of isolates mainly from the Bourgogne region. A subset of 34 isolates, representative of the fields sampled, was phenotyped for aggressiveness on a set of resistant and susceptible varieties of four legume hosts (pea, faba bean, vetch, alfalfa). Significant differences in disease severity were found among isolates and three groups of aggressiveness comprising 16, 17, and 2 isolates, respectively, were identified using HCA analysis. A higher diversity in pathogen aggressiveness was observed among isolates from Bourgogne, which included different legumes in its crop history. Little relationship was observed between genetic clusters and pathogenicity in the subset of 34 isolates, as expected using neutral markers. This study provides useful knowledge on the current state of low to moderate diversity among *A. euteiches* populations before resistant pea varieties are grown in France. New insights and hypotheses about the major factors shaping the diversity and evolution of *A. euteiches* are also discussed.

**Keywords:** Pisum sativum, legume crops, Aphanomyces root rot, genetic structure, host specificity

<footnote n="1">¹IGEPP, INRA, Agrocampus Ouest, Université de Rennes 1, Le Rheu, France, ²UMT PISOM INRA/Terres Inovia, Le Rheu, France, ³Terres Inovia, Thiverval Grignon, France, ⁴IGEPP, Agrocampus Ouest, INRA, Université de Rennes 1, Université Bretagne-Loire, Rennes, France</footnote>
# INTRODUCTION

*Aphanomyces euteiches* Drechsler is an oomycete pathogen of legumes, which causes the devastating root rot disease of pea (Pisum sativum L.) worldwide (Kraft and Pfleger, 2001). In Europe, *A. euteiches* was first observed in Norway in 1925 (Sundheim, 1972), and was reported a few years later in France (Labrousse, 1933), where it has been considered to be highly damaging in infested pea areas since 1993 (Didelot and Chaillet, 1995). The increased frequency of pea crops in French crop rotations since 1978 is believed to have favored the development of the disease in all pea-producing areas (Wicker and Rouxel, 2001). To date, the disease has been recorded in pea-growing regions such as the Bassin-Parisien, Bretagne, Rhônes-Alpes, Pyrénées, and Charente-Maritime (Wicker, 2001).

*A. euteiches* is a diploid, homothallic pathogen, i.e., having the capacity to complete its sexual cycle in the absence of any other individual (Lin and Heitman, 2007), that produces both oospores (sexual reproduction) and zoospores (asexual reproduction). Oospores can survive in the soil for more than 10 years (Papavizas and Ayers, 1974), and can resist unfavorable conditions (e.g., desiccation, freezing). *A. euteiches* was also reported to infect other legume species including alfalfa (Medicago sativa) (Delwiche et al., 1987), green bean (Phaseolus vulgaris) (Pfender and Hagedorn, 1982), faba bean (Vicia faba) (Lamari and Bernier, 1985) and common vetch (Vicia sativa) (Tsvetkova and Kotova, 1985). Currently, the main methods for managing the disease in France include avoidance of highly infested fields, which are diagnosed using an inoculum potential test, measuring disease severity on susceptible plants grown in field soil samples within a bioassay (Moussart et al., 2009), and crop rotations with non-host or resistant legume crops (Moussart et al., 2013). There are still no cultivated resistant varieties in France but consistent Quantitative Trait Loci (QTL) controlling partial resistance have been identified, validated in pea (Lavaud et al., 2015; Desgroux et al., 2016) and are being used in French pea breeding programs for the development of future resistant varieties. To support breeding for efficient resistance against *A. euteiches* populations present in French pea-growing regions more knowledge is required about the genetic diversity and the adaptive capacities of *A. euteiches* populations. A description of the state of *A. euteiches* population diversity before resistant varieties are grown is necessary to provide reference information for future comparisons with pathogenicity diversity after the development of resistant varieties.

Previous studies have described the diversity of *A. euteiches* populations worldwide. In the USA, results, obtained with dominant markers, showed high genetic diversity at the field and/or regional scales (Malvick et al., 1998; Grünwald and Hoheisel, 2006). Different genetic subpopulations were identified, which were distinguished based on the host of origin or host preference. However, all pea-infecting populations of *A. euteiches* showed significant linkage disequilibrium between markers, which suggested that selfing had played an important role in shaping the genetic structure of these populations (Grünwald and Hoheisel, 2006). In France, a preliminary study of the genetic diversity using AFLP markers found unstructured

populations among 56 pea *A. euteiches* isolates collected in infested fields (unpublished data). Recently, Mieuzet et al. (2016) and Le May et al. (2018) used SSR and Sequence Related Amplified Polymorphisms (SRAP) markers, respectively, to show the absence of genetic structure among French pea isolates sampled from research disease nurseries. However, in the same study, Le May also found significant genetic structure among American isolates sampled from four American nurseries.

Prior studies also analyzed phenotypic diversity among pea-infecting *A. euteiches* isolates based on their pathogenicity on legume species and pea lines. In our study, the term “pathogenicity” included two principal characteristics generally evaluated in pathogens, including virulence and aggressiveness. Virulence refers to the ability of the pathogen to cause a susceptible response on a host plant (Parlevliet, 2002). Aggressiveness refers to the quantitative variation measured on susceptible hosts. In practice, aggressiveness can be measured through a variety of quantitative traits, including infection efficiency, latent period, spore production rate, infectious period and lesion size, expressed during the host-pathogen interaction (Sackett and Mundt, 2005). In the USA, Malvick and Percich (1998) observed field-dependent variations in pathogenicity on pea, alfalfa, and bean hosts among 114 isolates from pea fields in Minnesota, Wisconsin and Oregon. In France, Wicker et al. (2001) identified four types of pathogenicity among 91 French isolates, according to their aggressiveness on five legume species. Most of the isolates were pathogenic on at least three of the five species. Wicker and Rouxel (2001) analyzed the pathogenic variability of 88 French pea-infecting isolates compared to 21 foreign isolates on a differential set of six pea genotypes. All the French isolates belonged to one major pathotype, i.e., virulence group (pathotype I) and a wide range of aggressiveness was observed between these isolates. Another pathotype (pathotype III), which showed lower aggressiveness toward the pea genotype MN313, was identified among American isolates.

Thus, low genetic and pathogenicity diversity among pea-infecting *A. euteiches* isolates was previously reported in France. However, these findings were obtained from collections of small to medium size isolates (<100 individuals) partially isolated from nurseries, using mostly dominant markers that could not reveal heterozygosity on the *A. euteiches* genome, and testing the strain aggressiveness on a limited set of genotypes or legume species. However, we can hypothesize that the diversity of *A. euteiches* populations in France is wider than previously described, because pea is not grown at the same frequency and density in the different French regions and other host legume crops can be grown as intercrops. In addition, *A. euteiches* isolates analyzed in previous studies do not represent the diversity of agronomic and pedo-climatic situations existing in the different French regions grown in peas. Therefore, this study addressed the following question: what is the genetic and pathogenicity diversity of *A. euteiches* populations in pea-growing regions of France? To answer this question, a large novel collection of pea-infecting isolates from major French pea-growing regions was established, as well as a set of differential genotypes of four legume species used in crop rotations in France. The collection was characterized using co-dominant and neutral SSR markers
and genetic structure of the French pathogen populations was analyzed. Pathogenicity diversity was evaluated on a subset of isolates representing the genetic diversity within the populations.

# MATERIALS AND METHODS

## Creation of the Collection of *A. euteiches* Isolates

A total of 205 *A. euteiches* isolates were obtained from soils sampled between 2011 and 2013 in 17 infested fields, distributed within eight departments corresponding to the main French pea-growing regions (Table 1). These fields belong to a network of plots for epidemic monitoring in the main production areas of pea. The choice of fields was based on the positive results of the inoculum potential tests (Moussart et al., 2009). Fifteen to twenty soil sub-samples were collected in each field at a 20–25 cm depth across the field, according to a spatial distribution in “W” (Campbell and Madden, 1990). Soils were sampled at regular intervals along the spatial pattern in “W.” The soil sub-samples from each field were mixed and 3 liters of soil per field were stored at 5°C. The baiting method used to sample isolates from each 3 liter-soil sample was adapted from Moussart et al. (2001), using the soil indexing method developed by Sherwood and Hagedorn (1958). Soils of each field were placed in pots (4 pots/field) and five plants of the pea cv. Lumina were sown in each pot (20 plants/field). A total of 68 pots was used to harvest the 205 *A.euteiches* isolates from the 17 infested fields. Four to twenty-two isolates were isolated from necrotic roots observed on 20 plants per soil, grown for 14 days, in a climatic chamber under favorable conditions to the disease (photoperiod 8–16 h, temperature 23–25°C, saturating moisture). The heterogeneous number of isolates per soil sample was due to variation in the number of symptomatic plants and level of plant infection between fields. No isolates were isolated from single-zoospore progeny since low level of polymorphism was previously reported between single-zoospore progeny and parental strains (Grünwald and Hoheisel, 2006).

## DNA Extraction and SSR Genotyping

DNA was extracted from mycelium samples of the 205 *A. euteiches* isolates. Mycelial explants were grown for 6 days at 25°C on Corn Meal Agar medium (CMA). Seven to ten agar discs (3 mm diameter) per culture were then transferred to peptone-glucose broth and grown for 6 days at 25°C. Mycelial mats were vacuum-filtered on Whatman paper, rinsed three times with sterile water, and transferred to Eppendorf tubes. The harvested mycelia were lyophilized and stored at –20°C. Twenty mg of mycelium were ground with a FastPrep® grinder and extracted with the Macherey-Nagel NucleoSpin® Plant II kit (Macherey-Nagel GmbH & Co. KG, France) according to the manufacturer's recommendations. The quality and quantity of DNA were evaluated using spectrophotometry and adjusted (Nanodrop ND-100; Nanodrop Technologies) to 10 ng/µL. The extracted DNAs were then stored at –20°C.

Twenty SSR primers, previously developed to amplify *A. euteiches* DNA (Mieuzet et al., 2016), were used in this study. The PCR reaction in simplex was carried out in 10 µL

containing 1X Green GoTaq® Flexi Buffer, 2 mM MgCl2, 0.2 mM each dNTP (Promega, France), 1 µM each forward and reverse primer and 0.25 µM fluorescent-labeled M13 primer (FAM, Applied Biosystem), 1 U of GoTaq Flexi DNA polymerase (final concentration) (Promega, France) and 20 ng template DNA. Volumes were adjusted to 10 µL with sterile distilled water. Amplifications were conducted on a S1000 Thermal Cycler (Bio-Rad). The cycling conditions for PCR amplification were as described in Mieuzet et al. (2016). The amplified products were separated on 1.5% agarose (Lonza) gels, stained in 1X Tris-Borate EDTA (TBE) buffer containing 1X Sybr® Safe (Invitrogen) and visualized under UV light. They were diluted 1:75 then 1:25 in sterile distilled water to decrease the signal intensity and to perform a clear analysis of peaks on GeneMapper 3.7 (Applied Biosystems) Two microliters of this dilution was then added to 10 µL HiDi formamide (Applied Biosystems) and 0.1 µl GeneScan 500 LIZ Size Standard (Applied Biosystems) then run on an ABI 3130xl Genetic Analyzer (Applied Biosystems). The presence of alleles and their sizes were assigned using GeneMapper 3.7 (Applied Biosystems). To check for the absence of genotyping error, the reference isolate RB84 and water were used as positive and negative control, respectively, in three replicates. Samples displaying non-obvious peaks or specific molecular patterns were duplicated.

## Genetic Diversity Analysis

SSR data were used to define MultiLocus Genotypes (MLGs) and checked for repeated MLGs. The number of repeated MLGs was identified using GENCLONE 2.0 (Arnaud-Haond and Belkhir, 2007). Genotypic evenness was evaluated using the index $R = (G-1)/(N-1)$, with G the number of distinct multilocus genotypes and N the number of isolates (Grünwald et al., 2003). The allelic richness (Ar) was estimated using the rarefaction method implemented in POPULATIONS 1.2.32 software (Langella, 1999), which estimated the number of alleles per locus for a reduced sample size. This analysis was carried out based on the smallest population size ($n = 4$) in a first estimation. In a second analysis small-size populations ($n < 10$) were rejected and an analysis was carried out only with the smallest remaining population size ($n = 12$). Observed (Ho) and expected (He) heterozygosity were computed using GenAlEx 6.5.2 (Peakall and Smouse, 2012). MICRO-CHECKER version 2.2.3 (Van Oosterhout et al., 2004) was used to estimate the presence of null alleles at the SSR markers. Fis (Fis=1-Ho/He) (Weir and Cockerham, 1984) was calculated using GENEPOP (Raymond and Rousset, 1995) with and without markers that displayed null alleles for each population. Clonality was assessed with the index of association ($I_A$) and $\bar{r}_d$ statistic, a measure of the multilocus linkage disequilibrium, calculated using Multilocus software version 3.1b (Agapow and Burt, 2001). The index of association between the scored alleles was estimated by comparing the variance of the genetic distances among alleles in the current dataset to the mean variance of 1,000 artificial re-sampled datasets. The $\bar{r}_d$ statistic is much less dependent on the number of loci than the index of association (Agapow and Burt, 2001). The $I_A$ and $\bar{r}_d$ are expected to be zero if populations are freely recombining and greater than zero if
<table><caption>TABLE 1 | Origin of the 205 *Aphanomyces euteiches* isolates of the collection, collected in different French fields between 2011 and 2013.</caption><thead><tr><th>Region</th><th>French department (area code)</th><th>Field (population code)</th><th>N</th><th>Latitude</th><th>Longitude</th></tr></thead><tbody><tr><td rowspan="2">Bretagne</td><td>Finistère (29)</td><td>Riec/Belon (RB)</td><td>6</td><td>47°52′40.69″N</td><td>3°42′28.57″O</td></tr><tr><td>Morbihan (56)</td><td>Bignan (BI)</td><td>6</td><td>47°52′51.30″N</td><td>2°44′21.70″O</td></tr><tr><td rowspan="4">Bourgogne</td><td rowspan="4">Côte d'Or (21)</td><td>Bretenière SO (BRO)</td><td>8</td><td>–</td><td>–</td></tr><tr><td>Bretenière S1 (BR1)</td><td>16</td><td>47°14′17.68″N</td><td>5°05′39.37″E</td></tr><tr><td>Bretenière S2 (BR2)</td><td>4</td><td>47°14′18.96″N</td><td>5°05′42.15″E</td></tr><tr><td>Bretenière S3 (BR3)</td><td>8</td><td>47°14′34.44″N</td><td>5°05′44.90″E</td></tr><tr><td rowspan="5">Center</td><td rowspan="5">Eure-et-Loir (28)</td><td>Houville La Branche SO (HBO)</td><td>12</td><td>48°27′03.55″N</td><td>1°38′25.93″E</td></tr><tr><td>Houville La Branche S1 (HB1)</td><td>22</td><td>48°27′23.30″N</td><td>1°38′20.98″E</td></tr><tr><td>Pierres (P)</td><td>20</td><td>48°34′21.54″N</td><td>1°32′07.79″E</td></tr><tr><td>Fresnay l'Evêque (F)</td><td>12</td><td>48°14′45.00″N</td><td>1°48′32.00″E</td></tr><tr><td>Bleville (BL)</td><td>13</td><td>48°19′47.80″N</td><td>2°22′96.20″E</td></tr><tr><td rowspan="2">Nord-Pas-de-Calais</td><td>Loiret (45)</td><td>Ennevelin SO (ENO)</td><td>21</td><td>50°31′50.12″N</td><td>3°07′45.10″E</td></tr><tr><td>Nord (59)</td><td>Ennevelin S1 (EN1)</td><td>8</td><td>50°31′50.71″N</td><td>3°08′09.90″E</td></tr><tr><td rowspan="3">Ile de France</td><td rowspan="3">Essonne (91)</td><td>Nainvilles Les Roches (NR)</td><td>16</td><td>48°51′96.00″N</td><td>2°49′00.00″E</td></tr><tr><td>Mondeville (MO)</td><td>9</td><td>48°28′57.00″N</td><td>2°25′15.00″E</td></tr><tr><td>Boigneville (BO)</td><td>11</td><td>48°32′94.90″N</td><td>2°38′72.20″E</td></tr><tr><td>Seine et Marne (77)</td><td></td><td>Crisenoy (CR)</td><td>13</td><td>48°35′43.13″N</td><td>2°45′10.58″E</td></tr></tbody></table>

N, Number of isolates sampled in each field.

linkage disequilibrium between alleles is maintained through selfing (clonality).

## Population Structure Analysis

Partition of molecular diversity among and within regions, as well as among and within populations, was studied using analysis of molecular variance (AMOVA) (Excoffier et al., 1992; Lynch and Milligan, 1994). AMOVA was performed using Arlequin 3.1 (Excoffier et al., 2005).

A principal coordinate analysis (PCoA) of the mean pairwise population genetic distance matrix was performed using the standardized genetic distance in GenALEx 6.5 (Peakall and Smouse, 2012). *A. euteiches* isolates were clustered with and without Bourgogne population using the Bayesian Clustering approach implemented in STRUCTURE 2.3.3 (Pritchard et al., 2000). The analysis was performed using 5 × 10⁵ burn-in replicates and a run length of 2 × 10⁵ Markov chain Monte Carlo (MCMC) replicates, adopting the admixed model and the correlated allele frequencies option. The number of genetic groups (K-value) was estimated using the model developed by Evanno et al. (2005), which provides an estimate of the posterior probability of the data for a given K, Pr (X/K). We used the height of the modal value of the distribution as an indicator of the strength of the signal detected by STRUCTURE software. Five independent runs were performed for K values between one and 10 in order to verify the convergence of parameter estimates. A Discriminant Analysis of Principal Components (DAPC) was also performed using the “adegenet” package of R software (Jombart et al., 2018), to confirm and describe the genetic clusters identified. Finally, a Minimum spanning network (MSN) using Nei’s distance for the different populations of *A. euteiches* was calculated with the “poppr” package (Kamvar et al., 2018) from R software, version 3.2.2 (R Core Team, 2014). The interactive

tool “imsn()” was used to create minimum spanning networks. Multilocus genotypes (MLG) were collapsed to multilocus haplotypes based on the minimum genetic distance at which two individuals would be considered from different clonal lineages. Haplotypes were represented by circles containing the number of isolates, and sized in proportion to haplotype frequency.

## Pathogenicity Testing

### Pathogen and Plant Material

The strain aggressiveness level was evaluated for 34 of the 205 *A. euteiches* isolates, including two isolates randomly sampled in each of the 17 fields used for the collection. The RB84 French reference strain (pathotype I, Moussart et al., 2007) was also used as a control.

Pathogenicity tests were conducted on eight genotypes from four leguminous species, which previously showed various levels of resistance vs. susceptibility to the RB84 strain of *A. euteiches* (Moussart et al., 2008). The plant genotypes included: pea (*Pisum sativum*) cv. Lumina (susceptible), gm. MN313 (susceptible to pathotype I isolates, but partially resistant to pathotype III isolates, Wicker and Rouxel, 2001) and gm. PI180693 (partially resistant); alfalfa (*Medicago sativa*) cv. Zenith (susceptible); vetch (*Vicia sativa*) cv. Amethyste (susceptible) and cv. Topaze (resistant); faba bean (*Vicia faba*) cv. Baraca (moderately susceptible) and cv. Melodie (resistant).

### Pathogenicity Tests and Disease Severity Assessment

The 34 isolates were assayed in an experiment comprising eight pathogenicity tests conducted in controlled conditions (Additional File 1). Each test included four to five isolates and the RB84 strain, individually inoculated on the eight genotypes. Each test comprised four replicates of four to five plants per
genotype and isolate tested. The experiment was performed twice, in 2016 and 2017, respectively, in the same climatic chamber, in order to confirm the results. A modified version of the standardized test developed for evaluating pea resistance to *A. euteiches* was used (Moussart et al., 2001). Seeds were sown in plastic pots (9*9*9.5 cm) containing unsterilized vermiculite (VERMEX, M). Faba bean seeds were soaked in water for 2 h before sowing. In each pot, five seeds of one pea, vetch or alfalfa genotype, and four seeds of one faba bean genotype were sown. Each pot constituted a replicate. Trays containing the four replicates (pots) were placed in a randomized design in a growth chamber (thermo period: 25/23°C and 16 h photoperiod). Seven days after sowing, each plant was inoculated by applying 5 ml of a zoospore suspension adjusted to 5.10³ spores/ml, as previously described by Moussart et al. (2001). After inoculation, the vermiculite was saturated with water to provide favorable conditions for infection. After 10 days, the plants were carefully removed from the vermiculite, the roots were washed in tap water and disease severity (DS) was scored on each plant using a 0-5 scoring scale (14): 0 = no symptoms; 1 = traces of discoloration on the roots (<25%); 2 = discoloration of 25 to 50% of the roots; 3 = discoloration of 50 to 75% of the roots; 4 = discoloration of >75% of the roots; 5 = plant dead.

## Statistical Analysis of Phenotypic Data

Statistical analysis was performed using R software, version 3.2.2 (R Core Team, 2014).

To check that a given isolate shows consistent results across experiments, correlations were estimated between DS score means obtained for (i) the different isolates in the two experiments and (ii) the RB84 strain in the sixteen pathogenicity tests over the two experiments, using Pearson coefficients (p-value = 0.05).

DS scores were analyzed as ordinal qualitative data using a cumulative link mixed model (CLMM; “clmm” function, “ordinal” package) (Christensen, 2015). In a first analysis, the DS score was considered as the dependant variable, the isolate, legume species, or genotypes as fixed factors and the replicate and experiment as random factors. ANOVA tests using “car” and “RVAideMemoire” packages (Hervé M., 2015; Fox and Weisberg, 2016) were performed to evaluate the legume species and genotype effects (p-value = 0.05) in the CLMM model. In a second analysis conducted for each genotype, the DS score was considered as the dependant variable, the isolate as a fixed factor and the replicate and experiment as random factors. ANOVA tests using “car” and “RVAideMemoire” packages (Hervé M., 2015; Fox and Weisberg, 2016) were carried out to evaluate the isolate effects (p-value = 0.05) in the CLMM model. Least square means (LSMeans) were calculated for each isolate and genotype using the “lsmeans” function of the “lsmeans” package (Lenth and Herve, 2015). Then, for each genotype, LSMeans were compared to the value obtained for the RB84 reference isolate in all tests with a Tukey test (p-value = 0.001), using the “cld” function of the “MultCompView” package (Graves et al., 2015). The probability to obtain each DS scores (0 to 5) for each isolate per genotype was calculated with the “rating lsmeans” function. The probability of DS scores was comprised between 0 and 1.

Virulence is usually associated with the pathogen's capacity to infect a specific host (Van der Plank, 1968) and to multiply in the host (Casadevall and Pirofski, 2001). Virulence phenotypes were defined based on the DS score probability obtained from the CLMM for the four susceptible legume genotypes. An isolate was declared virulent on a legume species (P: pea; V: vetch; F: faba bean; A: alfalfa) when the probability of symptom appearance on the susceptible legume genotype (DS score ≥1, presence of oospores in roots; Wicker et al., 2001) was greater than 0.5.

Isolates were assigned to previously described pea pathotypes (Wicker and Rouxel, 2001), based on their level of aggressiveness on the three pea genotypes used in this study. An isolate was assigned to pathotype I when LSMean scores obtained from CLMM on the three pea genotypes were significantly different (p-value < 0.05), with Lumina, the most susceptible genotype, MN313, the genotype with intermediate behavior and PI180693, the most resistant genotype. An isolate was assigned to pathotype III when the LSMean score obtained on MN313 was both significantly lower than on Lumina and equal to or lower than on PI180693 (Onfroy et al., personal communication).

A principal component analysis (PCA; “pca” function, “FactoMineR” package) (Husson et al., 2007) was performed using DS LSmeans scores for each isolate and genotype, in order to analyze the structure of *A. euteiches* isolates according to their phenotypic variability on the eight genotypes. A hierarchical clustering analysis (HCA) was then carried out to define groups of isolates using the “hclust” function of the “fastcluster” package (Müllner, 2015).

Statistical analysis was carried out in order to evaluate the differences between the DS scores of the different genetic clusters from the Bayesian analysis (STRUCTURE) on each plant genotype. DS scores were analyzed using a cumulative link mixed model (CLMM; “clm” function, “ordinal” package), in which the DS score was considered as the dependant variable, the genetic cluster and the genotypes as fixed factors and the replicate and experiment as random factor. ANOVA tests using “car” and “RVAideMemoire” packages were performed to evaluate the cluster and genotype effects (p-value = 0.05) in the CLMM model. In a second analysis conducted for each genotype, the DS score was considered as the dependent variable, the cluster as a fixed factor and the replicate and experiment as random factors. ANOVA tests using “car” and “RVAideMemoire” packages (Hervé M., 2015; Fox and Weisberg, 2016) were carried out to evaluate the cluster effects (p-value = 0.05) in the CLMM model. Least square means (LSMeans) were calculated for each cluster using the “lsmeans” function of the “lsmeans” package (Lenth and Herve, 2015). Then, for each genotype, LSMeans were compared with a Tukey test (p-value = 0.05), using the “cld” function of the “MultCompView” package (Graves et al., 2015).

# RESULTS

## Genetic Diversity and Structure of
*A. euteiches* Populations

Genotyping of the 205 *A. euteiches* isolates with the 20 SSRs revealed a total of 42 alleles, with one to three alleles
per SSR marker within the collection. Most of the isolates were homozygous at SSR loci. However, some isolates from the Bourgogne region displayed distinct patterns, with a lot of heterozygous loci (Additional File 2). A low level of SSR genotypic diversity was observed among isolates within fields, according to the number of distinct MLGs observed per field ($G \le 3$, except at BR1 and ENO; Table 2). Indeed, among the 205 isolates, GENCLONE analysis detected 19 different MLGs, including six MLGs differing from each other by one allele. Due to scoring errors, distinct MLGs belonging to the same clone can be found (Meirmans and Van Tienderen, 2004). As this was observed in this microsatellite data sets, each MLG pair presenting extremely low distance (i.e., 1–2 allelic differences) was grouped into the same multilocus lineages (Arnaud-Haond et al., 2007) resulting in a total of 13 MLG identified among the isolates in the collection (Additional File 3). MLGs 1 and 7 were the most represented, comprising 54 and 16% of the isolates studied, respectively. Genotypic richness (R) ranged from 0 (BRO, BR2, HBO, F, BO, and CR) to 0.333 (BR1). Only four fields (ENI, BR1, RB, and B1) showed an R index equal or greater than 0.20 (Table 2).

Genetic diversity indexes also showed low variability among isolates within fields, with a low level of heterozygosity and high proportion of allele fixation. Ar values did not exceed 1.84 alleles per population and were similar between populations except for BR1 and BR3 (Table 2). The rate of heterozygosity observed in each population was low, except for BR1 and BR3 (Ho = 0.269 and Ho = 0.788, respectively). Null alleles were present in the BR1, HB1, BL, P, EN0, EN1, NR, MO, and CR populations.

The excess of heterozygosity in the two fields BR1 and BR3 was confirmed by a negative *Fis* value, whether corrected or not with null alleles, whereas complete allele fixation was observed for the other fields except for the HB1 population which showed a quasi-complete allele fixation (Table 2). A majority of the populations had significant $I_A$ and $\bar{r}_d$ (*p*-value < 0.01), suggesting high multi-locus linkage disequilibrium and clonal populations (Table 2).

AMOVA confirmed the lack of genetic subdivision between regions and revealed that 53% (*p*-value < 0.001) and 45% (*p*-value < 0.001) of the total genetic variance was partitioned among populations and within populations, respectively, whereas 2% (*p*-value = 0.025) was attributable to differences among regions.

STRUCTURE analysis showed a maximum log likelihood score at $K=2$ (delta $K=100$) (Additional File 4A), suggesting no differentiation between populations except for most populations from Bourgogne (Figure 1A, Additional File 3). Results showed that 87.8 and 12.2% of the isolates belonged to cluster 1 and cluster 2, respectively, with moderate to high inference probabilities of isolates to cluster 1 (*p* > 0.99) and cluster 2 (*p* > 0.68). Cluster 1 included isolates with homozygous genotypic patterns from the five different French regions sampled. Cluster 2 included 25 isolates, among which 14 isolates from Bourgogne (BR1, BR3) showed heterozygous genotypic profiles. PCoA analysis confirmed distinct genetic structure of isolates from Bourgogne (BR1, BR3) compared to isolates from the other fields, with the first principal axis contributing to 52.11% of the variation (Figure 2A). The distinct genetic

<table><caption>TABLE 2 | SSR genotypic and genetic diversity parameters for each *Aphanomyces euteiches* population.</caption><thead><tr><th rowspan="2">Field<br>(population)<br>code</th><th colspan="3">SSR genotypic diversity</th><th colspan="6">Genetic diversity</th><th colspan="2">Gametic desequilibrium</th></tr><tr><th>N</th><th>G</th><th>R</th><th>Ar (n = 4)</th><th>Ar (n = 12)</th><th>He</th><th>Ho</th><th>Na</th><th>Fis non-<br>corrected</th><th>Fis<br>corrected</th><th>I<sub>A</sub></th><th>r&#773;<sub>d</sub></th></tr></thead><tbody><tr><td>RB</td><td>6</td><td>2</td><td>0.2</td><td>1.15</td><td>-</td><td>0.067</td><td>0</td><td>0</td><td>1</td><td>1</td><td>2.000***</td><td>1.000***</td></tr><tr><td>BI</td><td>6</td><td>2</td><td>0.2</td><td>1.05</td><td>-</td><td>0.022</td><td>0</td><td>0</td><td>1</td><td>1</td><td>-</td><td>-</td></tr><tr><td>BRO</td><td>8</td><td>1</td><td>0</td><td>1</td><td>-</td><td>0</td><td>0</td><td>0</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BR1</td><td>16</td><td>6</td><td>0.333</td><td>1.75</td><td>1.4</td><td>0.267</td><td>0.269</td><td>2</td><td>-0.03</td><td>-0.146</td><td>11.372***</td><td>0.769***</td></tr><tr><td>BR2</td><td>4</td><td>1</td><td>0</td><td>1</td><td>-</td><td>0</td><td>0</td><td>0</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BR3</td><td>8</td><td>2</td><td>0.14</td><td>1.84</td><td>-</td><td>0.403</td><td>0.788</td><td>0</td><td>-0.95</td><td>-0.95</td><td>1.000**</td><td>1.000***</td></tr><tr><td>HBO</td><td>12</td><td>1</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>HB1</td><td>22</td><td>2</td><td>0.048</td><td>1.14</td><td>1</td><td>0.053</td><td>0.005</td><td>3</td><td>0.78</td><td>-0.024</td><td>1.809***</td><td>0.677***</td></tr><tr><td>BL</td><td>13</td><td>2</td><td>0.083</td><td>1.05</td><td>1.05</td><td>0.021</td><td>0</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-</td></tr><tr><td>P</td><td>20</td><td>3</td><td>0.105</td><td>1.11</td><td>1.05</td><td>0.048</td><td>0.003</td><td>2</td><td>0.96</td><td>1</td><td>0.214</td><td>0.114</td></tr><tr><td>F</td><td>12</td><td>1</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>ENO</td><td>21</td><td>4</td><td>0.15</td><td>1.12</td><td>1.1</td><td>0.054</td><td>0.005</td><td>2</td><td>0.93</td><td>1</td><td>0.656***</td><td>0.349***</td></tr><tr><td>EN1</td><td>8</td><td>3</td><td>0.286</td><td>1.09</td><td>-</td><td>0.034</td><td>0</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-</td></tr><tr><td>NR</td><td>16</td><td>2</td><td>0.066</td><td>1.05</td><td>1.05</td><td>0.025</td><td>0</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-</td></tr><tr><td>MO</td><td>9</td><td>2</td><td>0.125</td><td>1.19</td><td>-</td><td>0.069</td><td>0</td><td>4</td><td>1</td><td>1</td><td>3.000***</td><td>1.000***</td></tr><tr><td>BO</td><td>12</td><td>1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>CR</td><td>13</td><td>1</td><td>0</td><td>1.05</td><td>1.05</td><td>0.018</td><td>0</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-</td></tr></tbody></table>

N, number of isolates; G, number of distinct multilocus genotype; R, genotypic richness; Ar, allelic richness corrected for sample size (n); He: expected heterozygosity without biais (Nei, 1978); Ho, observed heterozygosity; Na, number of SSR markers with null alleles; Fis, average deviation from Hardy-Weinberg proportions with (non-corrected) or without (corrected) null alleles; I<sub>A</sub> index of association and r&#773;<sub>d</sub> the standardized index of association (***p-value < 0.001,**p-value < 0.01, *p-value < 0.05, .p-value < 0.1); - missing data.
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | (A) Geographical position of fields distributed in five regions and Bayesian analysis (STRUCTURE) of the corresponding isolates, with Bourgogne (BR) data. Assignment probabilities of *A. euteiches* isolates are presented for *k* = 2 and a delta *k* = 100. Each vertical dotted line separates isolates from each of the 17 sites. (B) Bayesian analysis of the corresponding isolates, without Bourgogne (BR) data. Assignment probabilities of *A. euteiches* isolates are presented for *k* = 3 and a delta *k* = 1,600. Each vertical dotted line separates isolates from each of the 13 sites.</figcaption></figure>

structure of BR1 and BR3 populations was also confirmed by DAPC analysis (Additional File 5A).

Structure analysis excluding isolates from Bourgogne revealed a highest level of hierarchical sub-structure with *K* = 3 (delta *K* = 1600) (Figure 1B, Additional Files 3, 4B). However, this result is not fully consistent with the results obtained with PcoA and DAPC, which more highlighted a structure of *A. euteiches* populations according to their region of origin (Figure 2B, Additional File 5B).

The MSN analysis confirmed results obtained with STRUCTURE and PCoA on the entire collection, since the nineteen identified haplotypes were classified into two main groups separated by more than 0.25 Nei's genetic distance. One genetic group included isolates from Bourgogne in seven poorly represented haplotypes and the other genetic group included most of the isolates in 12 haplotypes. For the latter group there was no relationship between haplotype and geographical origin (Figure 3).

## Pathogenicity Diversity of *A. euteiches* Isolates

DS scores on the eight legume genotypes tested for all 34 isolates were significantly correlated between the two experiments (*r* > 0.84, *P* < 0.05). DS scores for the RB84 strain were significantly

correlated between the sixteen tests in the two experiments (*r* > 0.93, *P* < 0.05).

*A. euteiches* isolates caused significantly different levels of disease severity among the tested legume species (*p*-value = 6.34·10⁻¹⁵), genotypes (*p*-value = 2.07·10⁻⁸) and isolates (*p*-value < 2.10⁻¹⁶). Highest levels of disease were observed on pea, followed by alfalfa, vetch and faba bean. In each legume species, our results confirmed the expected response of resistance or susceptibility of the genotypes tested. In pea, Lumina, MN313 and PI180693 presented high, intermediate and low average DS values, respectively, in response to each of the isolates tested (Table 3).

Significant DS differences were observed between several isolates of the collection and RB84, on the pea genotype PI180693, the vetch genotype Améthyste and the faba bean genotype Baraca (*p*-value < 0.0049) (Table 3). In addition, some isolates had the same effects on different genotypes (low variance) and other had different effects (high variance) depending on the French regions and the plant genotype. Isolates from Bourgogne showed the highest DS variance on almost all the genotypes tested. The vetch genotype Topaze, as well as the faba bean genotype Baraca and mostly Melodie, had the most highly diverse DS scores in response to the isolates sampled from the different regions (Additional File 6). Twenty eight isolates were virulent on the four susceptible legume genotypes (PVFA
<figure><img src="image_5.png" /></figure>

**FIGURE 2** | Principal Coordinates Analysis (PCoA) conducted from SSR genotypic data on the *A.euteiches* isolates of the collection. (A) Analysis including all the fields. (B) Analysis performed without Bourgogne fields. Each color corresponds to regions: yellow Ile de France, green: Center, red: Bretagne, black: Nord Pas-de-Calais, blue: Bourgogne.

virulence phenotype profile). Four isolates (BR2-4, EN1-3, BO5, and BO13) were virulent on all the legumes tested but vetch (PFA profile) and two isolates from Bourgogne (BRO-2 and BRO-6) were virulent only on pea (P profile) (Table 3, Additional File 7). All 34 isolates were assigned to pea pathotype I (Table 3, Additional File 8).

Hierarchical clustering and PCA identified three clusters of pathogenicity (cluster A: 16 isolates, cluster B: 17 isolates, and cluster C: 2 isolates) among the 34 isolates (Figure 4A). Cluster C included two isolates from Bourgogne (BRO-2, BRO-6), isolated from a different field than those identified in genetic group 2 (BR1, BR3). The first PCA axis ($R^2$ = 58.04%) separated cluster C from clusters A and B (Figure 4A). Hierarchical clustering on the PCA analysis performed without isolates from Bourgogne separated isolates belonging to clusters A and B regardless of their region of origin, with a first PCA axis accounting for

39.26% of the total phenotypic variability (Figure 4B). Isolates belonging to cluster A included the reference isolate RB84 and showed a PVFA virulence phenotype profile. Isolates belonging to cluster B displayed PVFA and PFA profiles and the two isolates belonging to cluster C displayed a P virulence phenotype profile (Table 3).

A. *euteiches* isolates caused significantly different levels of disease severity among the genetic cluster and genotypes (*p*-value < 2.10⁻¹⁶). On each plant genotype, significant differences (*p*-value < 0.05) were observed between LSMeans of DS scores of two out of the three genetic clusters obtained from the Bayesian analysis (STRUCTURE), except on the vetch genotypes Amethyste and Topaze. On Amethyste, DS Score LSMeans were significantly different between the three genetic clusters, whereas on Topaze, these were not different between clusters (Additional File 9).
<figure><img src="image_6.png" /></figure>

**FIGURE 3** | Minimum spanning network (MSN) of 19 haplotypes detected in the *A. euteiches* collection of 205 isolates. Each circle represents a unique haplotype and the colors represent the sampling fields. The circle size represents the haplotype frequency and the number of isolates was indicated in circles. Line widths and the shading represent relatedness of the haplotypes based on Nei's genetic distance.

DISCUSSION

A. *euteiches* is a major devastating disease of spring pea in France. Intensification of spring crops and short rotations resulted in the significant development of the disease. This study is the first report analyzing both the genetic and pathogenicity diversity of *A. euteiches* populations from French pea-growing regions using co-dominant markers. This study gives insights into (i) the genetic diversity of French *A. euteiches* populations at the regional scale, (ii) the phenotypic diversity of isolates for pathogenicity on different legume hosts, and (iii) the relationship between genotypic and phenotypic group.

Most French *A. euteiches* Populations
Have Low Genetic Structure and Diversity

Based on the genetic polymorphism of SSR markers, this
study shows that the genetic structure of French *A. euteiches*
populations is low and not related to geographical origin. A
low level of genetic diversity in *A. euteiches* populations was
observed between the different locations in France, as shown
by the moderate number of multilocus genotypes identified (13

different MLGs within 205 isolates) including a predominant
MLG (MLG1) shared by 54% of the isolates studied. This
result is consistent with previous studies conducted on French
and American *A. euteiches* populations, using codominant or
dominant markers (Malvick and Percich, 1998; Wicker, 2001;
Grünwald and Hoheisel, 2006; Mieuzet et al., 2016; Le May et al.,
2018). In this study, the similar level of differentiation observed
within and among populations also confirmed the absence of
a clear genetic structure, as previously reported (Malvick and
Percich, 1998; Wicker et al., 2001; Grünwald and Hoheisel,
2006). The absence of resistance level toward the disease and the
deployment of the same pea susceptible cultivars across several
geographic area in France limit selection by host. This could
explain the low level of neutral genetic differentiation effects.
Moreover, founder effect which can lead to reduce the genetic
variation within the *A. euteiches* populations may not occur due
to the production of a large number of oospores at each cropping
season and the ability of oospores to be conserved in the soil for
many years. The high level of linkage disequilibrium observed
suggests clonal reproduction in the populations studied, in line
with the low level of genotypic diversity observed. These results
<table><caption>**TABLE 3** | Disease severity (DS) LSMean values obtained from CLMM analysis of DS ratings on roots of eight genotypes of four different legume species, in response to inoculation with 34 *A. euteiches* isolates from the collection (two isolates per field).</caption><thead><tr><th rowspan="2">Species<br>Isolates</th><th colspan="3">Pea</th><th colspan="2">Vetch</th><th colspan="2">Faba bean</th><th>Alfafa</th><th rowspan="2">Virulence<br>phenotype</th><th rowspan="2">Pathotype</th><th rowspan="2">MLG</th></tr><tr><th>LSMeans<br>Lumina</th><th>LSMeans<br>MN313</th><th>LSMeans<br>PI180693</th><th>LSMeans<br>Amethyste</th><th>LSMeans<br>Topaze</th><th>LSMeans<br>Baraca</th><th>LSMeans<br>Melodie</th><th>LSMeans<br>Zenith</th></tr></thead><tbody><tr><td>RB84</td><td>12.328</td><td>6.813</td><td>2.236</td><td>4.119</td><td>-3.785</td><td>1.111</td><td>-2.143</td><td>1.274</td><td>PVFA</td><td>I</td><td>3</td></tr><tr><td>RB3</td><td>5.892</td><td>4.830</td><td>-0.537***</td><td>-0.391***</td><td>-1.729</td><td>-1.861***</td><td>-5.595</td><td>-1.150***</td><td>PVFA</td><td>I</td><td>3</td></tr><tr><td>RB5</td><td>8.316</td><td>5.616</td><td>1.089</td><td>0.542***</td><td>-3.856</td><td>-1.899***</td><td>-5.615</td><td>-0.071</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>BI1</td><td>7.285</td><td>6.764</td><td>0.614</td><td>1.612</td><td>-4.255</td><td>-2.261***</td><td>-3.853</td><td>0.067</td><td>PVFA</td><td>I</td><td>4</td></tr><tr><td>BI6</td><td>8.316</td><td>7.845</td><td>1.292</td><td>1.565</td><td>-4.652</td><td>-1.655</td><td>-3.294</td><td>0.168</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>BRO-2</td><td>-1.013</td><td>-3.535***</td><td>-4.583***</td><td>-9.432***</td><td>-2.320</td><td>-5.620***</td><td>-5.623</td><td>-5.468***</td><td>P</td><td>I</td><td>1</td></tr><tr><td>BRO-6</td><td>0.444</td><td>-1.891***</td><td>-2.775***</td><td>-7.566***</td><td>-1.547***</td><td>-5.535***</td><td>-5.603</td><td>-4.034***</td><td>P</td><td>I</td><td>1</td></tr><tr><td>BR1-2</td><td>2.452</td><td>3.888</td><td>-0.753***</td><td>1.121</td><td>-3.426</td><td>0.494</td><td>-2.671</td><td>2.258</td><td>PVFA</td><td>I</td><td>12</td></tr><tr><td>BR1-3</td><td>5.396</td><td>5.975</td><td>-0.762***</td><td>-0.654***</td><td>-2.977</td><td>0.649</td><td>-1.644</td><td>1.650</td><td>PVFA</td><td>I</td><td>11</td></tr><tr><td>BR2-1</td><td>7.105</td><td>5.752</td><td>0.994</td><td>-3.603***</td><td>-5.411</td><td>-1.829***</td><td>-3.510</td><td>-0.027</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>BR2-4</td><td>9.035</td><td>6.675</td><td>0.923</td><td>-6.639***</td><td>-4.555</td><td>-2.508***</td><td>-5.656</td><td>-2.495***</td><td>PFA</td><td>I</td><td>1</td></tr><tr><td>BR3-1</td><td>5.547</td><td>5.901</td><td>-1.939***</td><td>3.124</td><td>-4.271</td><td>1.175</td><td>-1.665</td><td>1.804</td><td>PVFA</td><td>I</td><td>11</td></tr><tr><td>BR3-5</td><td>5.749</td><td>5.978</td><td>-1.003***</td><td>-0.949***</td><td>-3.497</td><td>0.707</td><td>-0.765</td><td>0.009</td><td>PVFA</td><td>I</td><td>10</td></tr><tr><td>HBO-1</td><td>8.175</td><td>4.128</td><td>-0.130</td><td>-1.640***</td><td>-4.317</td><td>-1.401</td><td>-4.028</td><td>1.939</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>HBO-9</td><td>9.012</td><td>4.981</td><td>1.829</td><td>1.991</td><td>-3.003</td><td>-0.018</td><td>-2.898</td><td>-0.617</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>HB1-3</td><td>7.536</td><td>6.546</td><td>-0.954***</td><td>-1.058***</td><td>-4.687</td><td>-1.396</td><td>-3.355</td><td>-1.135</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>HB1-14</td><td>8.982</td><td>5.966</td><td>4.200</td><td>0.919</td><td>-5.411</td><td>-0.305</td><td>-1.313</td><td>0.560</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>BL5</td><td>8.941</td><td>5.617</td><td>2.057</td><td>3.666</td><td>-3.294</td><td>1.002</td><td>-3.260</td><td>0.485</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>BL8</td><td>9.012</td><td>4.371</td><td>1.713</td><td>1.414</td><td>-3.679</td><td>-0.575</td><td>-3.826</td><td>0.969</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>P8</td><td>8.292</td><td>5.737</td><td>0.242</td><td>0.969</td><td>-4.688</td><td>-0.978</td><td>-2.772</td><td>0.181</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>P14</td><td>8.966</td><td>4.687</td><td>0.610</td><td>0.919</td><td>-4.688</td><td>-0.365</td><td>-3.547</td><td>-0.357</td><td>PVFA</td><td>I</td><td>7</td></tr><tr><td>F3</td><td>9.035</td><td>6.081</td><td>3.466</td><td>2.439</td><td>-5.447</td><td>0.377</td><td>-1.589</td><td>-0.264</td><td>PVFA</td><td>I</td><td>8</td></tr><tr><td>F10</td><td>8.959</td><td>7.125</td><td>2.641</td><td>2.133</td><td>-4.634</td><td>0.339</td><td>-1.074</td><td>-0.123</td><td>PVFA</td><td>I</td><td>8</td></tr><tr><td>ENO-12</td><td>8.975</td><td>7.479</td><td>1.317</td><td>2.298</td><td>-3.954</td><td>0.634</td><td>-0.850</td><td>1.663</td><td>PVFA</td><td>I</td><td>4</td></tr><tr><td>ENO-16</td><td>8.064</td><td>5.592</td><td>1.121</td><td>2.466</td><td>-3.570</td><td>0.211</td><td>-2.630</td><td>2.494</td><td>PVFA</td><td>I</td><td>5</td></tr><tr><td>EN1-3</td><td>8.175</td><td>7.312</td><td>2.270</td><td>-8.745***</td><td>-3.109</td><td>-0.861</td><td>-3.242</td><td>1.825</td><td>PFA</td><td>I</td><td>7</td></tr><tr><td>EN1-5</td><td>9.035</td><td>6.196</td><td>0.451</td><td>1.633</td><td>-3.131</td><td>0.495</td><td>-2.399</td><td>0.582</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>NR8</td><td>9.012</td><td>4.660</td><td>1.087</td><td>1.979</td><td>-4.707</td><td>-1.776***</td><td>-2.198</td><td>-0.240</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>NR14</td><td>8.982</td><td>6.262</td><td>2.903</td><td>2.411</td><td>-3.952</td><td>0.605</td><td>-2.097</td><td>0.095</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>MO1</td><td>7.853</td><td>6.783</td><td>1.920</td><td>0.843</td><td>-4.254</td><td>-1.298</td><td>-4.846</td><td>0.510</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>MO5</td><td>9.006</td><td>5.584</td><td>2.487</td><td>0.558***</td><td>-5.395</td><td>-1.065</td><td>-4.914</td><td>0.003</td><td>PVFA</td><td>I</td><td>6</td></tr><tr><td>BO5</td><td>7.291</td><td>5.211</td><td>1.982</td><td>-11.334***</td><td>-4.271</td><td>-2.422***</td><td>-2.769</td><td>-0.673</td><td>PFA</td><td>I</td><td>7</td></tr><tr><td>BO13</td><td>7.105</td><td>5.932</td><td>2.507</td><td>-10.555***</td><td>-3.698</td><td>-1.906***</td><td>-3.890</td><td>-1.883***</td><td>PFA</td><td>I</td><td>7</td></tr><tr><td>CR5</td><td>9.035</td><td>4.925</td><td>0.390</td><td>-1.097***</td><td>-5.410</td><td>-0.168</td><td>-2.871</td><td>0.770</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>CR9</td><td>8.989</td><td>6.161</td><td>1.636</td><td>1.427</td><td>-4.688</td><td>0.252</td><td>-2.435</td><td>0.515</td><td>PVFA</td><td>I</td><td>1</td></tr><tr><td>Mean</td><td>7.580</td><td>5.370</td><td>0.873</td><td>-0.672</td><td>-4.008</td><td>-0.847</td><td>-3.155</td><td>0.037</td><td></td><td></td><td></td></tr><tr><td>Standard Deviation</td><td>2.556</td><td>2.234</td><td>1.780</td><td>4.211</td><td>0.986</td><td>1.622</td><td>1.432</td><td>1.640</td><td></td><td></td><td></td></tr></tbody></table>

LSMean scores highlighted in gray with the isolates tested are significantly different from the LSMean score obtained with the RB84 reference strain for the same genotype (***p-value < 0.001). A Virulence phenotype is defined on pea (P), vetch (V), faba bean (F) and alfalfa (A) when the probability of DS score ≥ 1 on the four susceptible legumes (Pea: Lumina; Vetch: Amethyste; Faba bean: Baraca; Alfalfa: Zenith) is equal or higher than 0.5, according to the CLMM analysis. Pathotype I was attributed for all the isolates since LSMean scores between the three pea genotypes were significantly different (p-value < 0.05), with Lumina, MN313 and PI180693 being the susceptible, intermediate and partially resistant genotypes, respectively (Additional File 8). MLGs (Multilocus groups) are indicated according to the genetic diversity analysis conducted in this study.

are consistent with the homothallic mode of reproduction of the pathogen and its lower dissemination in the soil than that of airborne pathogens in the air. These *A. euteiches* life traits contribute to limit genetic mixing and maintain highly inbred pathogen populations.

However, *A. euteiches* populations from Bourgogne appeared to have unique characteristics compared to the other locations sampled in this study. Indeed, some *A. euteiches* isolates from Bourgogne showed a specific molecular pattern, with a high level of heterozygosity, which contributed to the high level of
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | PCA of LSMean disease severity data obtained on eight legume genotypes in response to inoculation with 34 isolates of *A.euteiches*. (A) Analysis included all the field sites. (B) Analysis was performed without Bourgogne field sites. Each color correspond to group from the hierarchical clustering ascendant analysis: black, field sites from group A, red: field sites from group B and green, field sites from group C.</figcaption></figure>

genetic diversity observed in this region. SSR markers used in our study allowed heterozygosity to be revealed in populations from Bourgogne, whereas no previous study using dominant markers could be able to reveal such genetic profiles. An excess of heterozygous genotypes was observed for BR1 and BR3 populations, which clustered separately from the other populations in the analysis and DAPC analysis. Our results are consistent with those of Grünwald and Hoheisel (2006), which suggested that population structure of the homothallic *A. euteiches* pathogen is mostly determined by regular selfing, but also occasional recombination, indicating a mixed mode of reproduction in *A. euteiches* populations. This complex pattern of sexuality was identified in other oomycetes (Francis and Stclair, 1993; Whisson et al., 1994), suggesting that outcrossing could occur even in pathogen populations with a high level of selfing. Indeed, Grünwald and Hoheisel (2006) suggested that outcrossing and migration occurred, albeit rarely, and contributed to the genetic diversity and differentiation observed in *A. euteiches* populations sampled from two fields in Oregon and Washington in the USA.

In our study, heterozygous isolates were obtained from fields in the Bourgogne region, which included legumes (faba bean, vetch, and alfalfa) other than pea in their cropping history, in contrast to all the other fields sampled, which only had a history of pea production. We could thus hypothesize that the occurrence of these isolates may result from outcrossing between genetically distinct isolates adapted to different legumes. *A. euteiches* is reported to attack other legume species including common bean, broad bean, faba bean, clover, and alfalfa (Pfender and Hagedorn, 1982; Lamari and Bernier, 1985; Tivoli et al., 2006; Moussart et al., 2008). In some regions of the United States, where pea and alfalfa crops were frequently included in cropping systems, populations of pea-infecting *A. euteiches* showed differentiation into sub-populations with differences in genotypes and virulence toward pea and alfalfa

(Holub et al., 1991; Malvick et al., 1998, 2009; Malvick and Grau, 2001). To investigate the presence of such distinct isolate genotypes, influenced by their host of origin in France, it would thus be necessary to increase the number of isolates sampled from the Bourgogne fields studied and from other fields with various legumes in their crop histories.

## Pathogenicity Diversity of French

### *A. euteiches* Populations Depends on Plant Host and Genotypes

Low to moderate pathogenicity diversity and structure was observed among the 34 isolates sampled in the collection from the different locations studied, depending on isolates and plant hosts and genotypes. Most isolates showed high aggressiveness on pea and virulence on vetch, faba bean and alfalfa (PVFA) with variable aggressiveness. Moreover, the level of aggressiveness of the isolates recorded on pea differential set used in this study indicated that all 34 isolates belong to pathotype I. Because of the heaviness of pathological tests and low genetic structure of populations studied, sub-samples of two isolates per field population were chosen to have a good geographic representation of the aggressiveness of French populations. The pathogenicity diversity observed within the French *A. euteiches* populations studied is consistent with previous studies. Using 91 pea-infecting isolates of *A. euteiches*, Wicker et al. (2001) reported that French isolates were able to infect a wide range of legume species with a low diversity of aggressiveness on pea genotypes and a moderate diversity of aggressiveness on vetch and alfalfa. The prevalence of pea crops in France since the 1980s and the susceptibility to *A. euteiches* of the pea cultivars used by the growers may explain the high aggressiveness and small variation of pathogenicity on pea of most of the French isolates tested. Moussart et al. (2008) reported a lower variability of disease response to *A. euteiches* infection in pea than in vetch and faba bean, with most varieties showing high level of susceptibility.
In this study, a poor relationship between genetic clusters and pathogenicity was observed in the subset of 34 isolates analyzed. Despite the low isolate-sample size used in this study to compare genotypic and pathogenicity diversity, the type of markers used could explain the lack of relationship observed between genetic and phenotypic structure of *A. euteiches* populations. Indeed, the neutral nature of SSR markers would probably not allow the *A. euteiches* isolates to be genotyped in genomic regions under selection. Similar results were obtained previously using RAPD or AFLP markers. Malvick and Percich (1998) did not find any correlation between pathogenic and genetic diversity in four American populations of *A. euteiches* from Minnesota, Wisconsin and Oregon, despite pathogenic variation being observed among the populations on a differential pea set. Grünwald and Hoheisel (2006) reported discrepancies between the levels of genetic and pathogenic diversity observed in two *A. euteiches* populations baited on pea from two fields located in Athena, Oregon and Mount Vernon, Washington. The use of Next Generation Sequencing methodologies to sequence the genome of *A. euteiches* isolates (Madoui et al., 2007; Gaulin et al., 2008) and further knowledge of pathogen effectors (Ramirez-Garces et al., 2016) would help to decipher the relationship between genetic and pathogenic diversity in *A. euteiches* populations.

## CONCLUSION

A low level of genetic and phenotypic diversity among French *A. euteiches* populations was observed in this study, based on the analysis of a collection of 205 isolates sampled from the major pea-growing regions in France. However, significant differences in aggressiveness were observed between several isolates on some genotypes. In addition, some isolates originating from fields where other legumes besides pea have been grown, also showed a distinct genetic structure with heterozygous genotypic patterns. This could suggest that crosses may have occurred between isolates which may have evolved on different hosts. Up till now, French *A. euteiches* populations have been subjected to limited selective pressure from hosts, since the diversification of legume species in rotation was restricted and no pea resistant varieties have yet been cultivated in France. Breeding for resistance is making progress, especially due to cumulating stable resistance QTL which were identified and recently confirmed in germplasm (Desgroux et al., 2016; Lavaud et al., 2016). Future partially resistant pea varieties and the increasing diversification

of cropping systems with host legumes will likely modify selection pressures on the pathogen populations and thus their genetic and pathogenic structures, as suggested in this study from some isolates. A better understanding and management of the evolutionary forces affecting *A. euteiches* populations will thus be required, in order to develop integrated control strategies for the durable management of Aphanomyces root rot disease.

## AUTHOR CONTRIBUTIONS

AQ-H and GL generated the phenotypic and genotypic data. AQ-H and CL performed the statistical and genetic analyses and drafted the manuscript. AM contributed to the establishment of the collection. AM, AB, and DA provided scientific expertise on the conception of the study. CL and M-LP-N coordinated and supervised the manuscript writing and the experiments. All the authors approved the final draft of the manuscript.

## FUNDING

This study was supported by INRA, AAP 2010 Gestion Durable des Résistances, the PeaMUST project, which received funding from the French Government managed by the Research National Agency (ANR) under the Investments for the Future call 2011 (ANR-11-BTBR-0002), and the PROGRAILIVE project (grant: RBRE160116CR0530019) funded by the Bretagne region, France, and European FEADER grants.

## ACKNOWLEDGMENTS

The authors greatly thank Dr. C. Lavaud (INRA IGEPP, France) for his helpful advice on statistical analysis, as well as G. Boutet and A. Lesné (INRA IGEPP, France) for their help using GeneMapper software. We also thank M-N. Even and C. Onfroy (Terres Inovia, UMT PISOM) for providing soil and several isolates. We acknowledge the greenhouse and experimental device plat-form of IGEPP for providing and managing equipment for the experiments.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01673/full#supplementary-material

## REFERENCES

* Agapow, P. M., and Burt, A. (2001). Indices of multilocus linkage disequilibrium. *Mol. Ecol.* Notes 1, 101–102. doi: 10.1046/j.1471-8278.2000.00014.x

* Arnaud-Haond, S., and Belkhir, K. (2007). GENCLONE: a computer program to analyse genotypic data, test for clonality and describe spatial clonal organization. *Mol. Ecol.* Notes 7, 15–17. doi: 10.1111/j.1471-8286.2006.01522.x

* Arnaud-Haond, S., Migliaccio, M., Diaz-Almela, E., Teixeira, S., Vand de Vliet, M. J., Alberto, F., et al. (2007). Vicariance patterns in the Mediterranean Sea: east-west cleavage and low dispersal in the endemic seagrass *Posidonia oceanica*. *J. Biogeogr.* 34, 963–976. doi: 10.1111/j.1365-2699.2006.01671.x

* Campbell, C. L., and Madden, L. V. (1990). *Introduction to Plant Disease Epidemiology*. New York, NY: J. Wiley and Sons.

* Casadevall, A., and Pirofski, L. (2001). Host-pathogen interactions: the attributes of virulence. *J. Infect. Dis.* 184, 337–344. doi: 10.1086/322044

* Christensen, R. (2015). *Regression Models For Ordinal Data*. Available Online at: https://cran.r-project.org/package=ordinal

* Delwiche, P. A., Grau, C. R., Holub, E. B., and Perry, J. B. (1987). Characterization of *Aphanomyces euteiches* isolates recovered from alfalfa in Wisconsin. *Plant Dis.* 71, 155–161. doi: 10.1094/PD-71-0155

* Desgroux, A., L'Anthoëne, V., Roux-Duparque, M., Rivière, J. P., Aubert, G., Tayeh, N., et al. (2016). Genome-wide association mapping of partial
resistance to *Aphanomyces euteiches* in pea. *BMC Genomics* 20:124. doi: 10.1186/s12864-016-2429-4

Didelot, D., and Chaillet, I. (1995). Relevance and interest of root disease prediction tests for pea crop in France in AEP (ed) 2nd European Conference of Grain Legumes. July 9-13th, Copenhagen.

Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. *Mol. Ecol.* 14, 2611-2620. doi: 10.1111/j.1365-294X.2005.02553.x

Excoffier, L., Laval, G., and Schneider, S. (2005). Arlequin (version 3.0): an integrated software package for population genetics data analysis. *Evol. Bioinform.* 1, 47-50. doi: 10.1177/117693430500100003

Excoffier, L., Smouse, P. E., and Quattro, J. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. *Genetics* 131, 479-491

Fox, J., and Weisberg, S. (2016). *Companion to Applied Regression*. Available online at: https://cran.r-project.org/package=car

Francis, D. M., and Stclair, D. A. (1993). Outcrossing in the homothallic oomycete, *Phythium ultimum*, detected with molecular markers. *Curr. Genet.* 24, 100-106. doi: 10.1007/BF00324672

Gaulin, E., Madoui, M. A., Bottin, A., Jacquet, C., Mathe, C., Couloux, A., et al. (2008). Transcriptome of *Aphanomyces euteiches*: new oomycete putative pathogenicity factors and metabolic pathways. *PLoS ONE* 3:e1723. doi: 10.1371/journal.pone.0001723

Graves, S., Piepho, H., and Selzer, L. (2015). *Visualizations of Paired Comparisons*. Available online at: https://cran.r-project.org/web/packages/multcompView/index.html

Grünwald, N. J., Goodwin, S. B., Milgroom, M. G., and Fry, W. E. (2003). Analysis of genotypic diversity data for populations of microorganisms. *Phytopathology* 93, 738-746. doi: 10.1094/PHYTO.2003.93.6.738

Grünwald, N. J., and Hoheisel, G. A. (2006). Hierarchical analysis of diversity, selfing, and genetic differentiation in populations of the oomycete *Aphanomyces euteiches*. *Phytopathology* 96,1134-1141. doi: 10.1094/PHYTO-96-1134

Hervé M. (2015). *RVAideMemoire: Diverse Basic Statistical and Graphical Functions*. Available online at: https://cran.r-project.org/package=RVAideMemoire

Holub, E. B., Grau, C. R., and Parke, J. L. (1991). Evaluation of the forma-specialis concept in *Aphanomyces euteiches*. *Mycol. Res.* 95, 147-157. doi: 10.1016/S0953-7562(91)80104-6

Husson, F., Josse, J., Lê S and Mazet J (2007). *Factor Analysis and Data Mining with R*. Available online at: https://cran.r-project.org/package=FactoMineR

Jombart, T., Kamvar, Z. N., Collins, C., Lustrik, R., Beugin, M. P., Knaus, B. J., et al. (2018). *Adegenet: Exploratory Analysis of Genetic And Genomic Data*. Available online at: https://cran.r-project.org/web/packages/adegenet/index.html

Kamvar, Z. N., Tabina, J. F., Everhart, S. E., Brooks, J. C., Krueger-Hadfield, S. A., Sotka, E., et al. (2018). *Genetic Analysis of Populations With Mixed Reproduction*. Available online at: https://cran.r-project.org/web/packages/poppr/index.html

Kraft, J. M., and Pfleger, F. L. (2001). *Compendium of Pea Diseases and Pests*. 2edn. St.Paul, MN: The American Phytopathological Society Press.

Labrousse, F. (1933). Notes de pathologie végétale. *Revue de pathologie végétale et d'Entomologie agricole* 19, 71-84.

Lamari, L., and Bernier, C. C. (1985). Etiology of seedling blight and root rot of Faba bean (Vicia faba) in Manitoba. *Can. J. Plant Pathol.* 7,139-145.

Langella, O. (1999). *Population Genetic Software*. Available online at: http://bioinformatics.org/~tryphon/populations

Lavaud, C., Bavière, M., Roy, G. L., Hervé M., Moussart, A., Delourme, R., et al. (2016). Single and multiple resistance QTL reduce *Aphanomyces euteiches* infection and colonization in roots of pea Near-Isogenic-Lines. *BMC Plant Biol.* 16:166. doi: 10.1186/s12870-016-0822-4

Lavaud, C., Lesne, A., Piriou, C., Le Roy, G., Boutet, G., Moussart, A., et al. (2015). Validation of QTL for resistance to *Aphanomyces euteiches* in different pea genetic backgrounds using near-isogenic lines. *Theor. Appl. Genet.* 128, 2273-2288. doi: 10.1007/s00122-015-2583-0

Le May, C., Onfroy, C., Moussart, A., Andrivon, D., Baranger, A., Pilet-Nayel, M. L., et al. (2018). Genetic structure of *Aphanomyces euteiches* populations sampled from United States and France pea nurseries. *Eur. J. Plant Pathol.* 150, 275-286. doi: 10.1007/s10658-017-1274-x

Lenth, R., and Herve, M. (2015). *Lsmeans: Least-Squares Means*. Available online at: https://cran.r-project.org/package=LSmeans

Lin, X., and Heitman, J. (2007). “Mechanisms of homothallism in fungi-transitions between heterothallism and homothallism,” in *Sex in Fungi: molecular determination and evolutionary implications*, eds J. Heitman, J. W. Kronstad, J. W. Taylor, L. A. Casselton (Washington DC: ASM Press), 35-58.

Lynch, M., and Milligan, B. G. (1994). Analysis of population genetic structure with RAPD markers. *Mol. Ecol.* 3, 91-99. doi: 10.1111/j.1365-294X.1994.tb00109.x

Madoui, M. A., Gaulin, E., Mathé C., San Clemente, H., Couloux, A., Wincker, P., et al. (2007). AphanoDB: a genomic resource for Aphanomyces pathogens. *BMC Genomics* 8:471. doi:10.1186/1471-2164-8-471

Malvick, D., Grunwald, N., and Dyer, A. (2009). Population structure, races, and host range of *Aphanomyces euteiches* from alfalfa production fields in the central USA. *Eur. J. Plant Pathol.* 123, 171-182. doi: 10.1007/s10658-008-9354-6

Malvick, D. K., and Grau, C. R. (2001). Characteristics and frequency of *Aphanomyces euteiches* races 1 and 2 associated with alfalfa in the midwestern United States. *Plant Dis.* 85, 740-744. doi: 10.1094/PDIS.2001.85.7.740

Malvick, D. K., Grau, C. R., and Percich, J. A. (1998). Characterization of *Aphanomyces euteiches* strains based on pathogenicity tests and random amplified polymorphic DNA analyses. *Mycol. Res.* 102,465-475. doi: 10.1017/S0953756297005029

Malvick, D. K., and Percich, J. A. (1998). Genotypic and pathogenic diversity among pea-infecting strains of *Aphanomyces euteiches* from the central and western United States. *Phytopathology* 88, 915-921. doi: 10.1094/PHYTO.1998.88.9.915

Meirmans, P. G., and Van Tienderen, P. H. (2004). GENOTYPE and GENODIVE: two programs for the analysis of genetic diversity of asexual organisms. *Mol. Ecol. Notes* 4, 792-794. doi: 10.1111/j.1471-8286.2004.00770.x

Mieuzet, L., Quillévéré A., Pilet, M. L., and Le May, C. (2016). Development and characterization of microsatellite markers for the oomycete *Aphanomyces euteiches*. *Fungal Genet. Biol.* 91, 1-5. doi: 10.1016/j.fgb.2016.03.001

Moussart, A., Even, M. N., Lesné A., and Tivoli, B. (2013). Successive legumes tested in a greenhouse crop rotation experiment modify the inoculum potential of soils naturally infested by *Aphanomyces euteiches*. *Plant Pathol.* 62, 545-551 doi: 10.1111/j.1365-3059.2012.02679.x

Moussart, A., Even, M. N., and Tivoli, B. (2008). Reaction of genotypes from several species of grain and forage legumes to infection with a French pea isolate of the oomycete *Aphanomyces euteiches*. *Eur. J. Plant Pathol.* 122, 321-333. doi: 10.1007/s10658-008-9297-y

Moussart, A., Onfroy, C., Lesné A., Esquibet, M., Grenier, E., and Tivoli, B. (2007). Host status and reaction of *Medicago truncatula* accessions to infection by three major pathogens of pea (*Pisum sativum*) and alfalfa (*Medicago sativa*). *Eur. J. Plant Pathol.* 117, 57-69. doi: 10.1007/s10658-006-9071-y

Moussart, A., Wicker, E., Duparque, M., and Rouxfel, F. (2001). “Development of an efficient screening test for pea resistance to *Aphanomyces euteiches*,” in AEP, ed 4th European Conference on Grain Legumes. July 8-12th) Cracow, Poland

Moussart, A., Wicker, E., Le Delliou, B., Abelard, J. M., Esnault, R., Lemarchand, E., et al. (2009). Spatial distribution of *Aphanomyces euteiches* inoculums in a naturally infested pea fields. *Eur. J. Plant Pathol.* 123, 153-158 doi: 10.1007/s10658-008-9350-x

Müllner, D. (2015). *Fast Hierarchical Clustering Routines for R and Python*. Available online at: https://cran.r-project.org/package=fastcluster

Nei, M. (1978). Estimation of average heterozygosity and genetic distance from a small number of individuals. *Genetics*. 89, 583-590.

Papavizas, G., and Ayers, W. (1974). *Aphanomyces species and their root diseases on pea and sugarbeet*. US Department of Agricultural Research Technical Bulletin 1484, 1-157

Parlevliet, J. E. (2002). Durability of resistance against fungal, bacterial and viral pathogens: present situation. *Euphytica* 124, 147-156. doi: 10.1023/A:1015601731446

Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5. *Bioinformatics* 28, 2537-2539. doi: 10.1093/bioinformatics/bts460

PFender, W. F., and Hagedorn, D. J. (1982). *Aphanomyces euteiches* f.sp. phaseoli, a causal agent of bean root and hypocotyl rot. *Phytopathology* 72, 306-310. doi: 10.1094/Phyto-72-306

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. *Genetics* 155, 945-959
R Core Team (2014). *R: A Language and Environment for Statistical Computing*. Vienna: R foundation for Statistical Computing

Ramirez-Garces, D., Camborde, L., Pel, M. J. C., Jauneau, A., Martinez, Y., Néant, I., et al. (2016). CRN13 candidate effectors from plant and animal eukaryotic pathogens are DNA-binding proteins which trigger host DNA damage response. *New Phytologist*. 210, 602–617. doi: 10.1111/nph.13774

Raymond, M., and Rousset, F. (1995). GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism. *J. Hered*. 86, 248–249. doi: 10.1093/oxfordjournals.jhered.a11573

Sackett, K. E., and Mundt, C. C. (2005). Primary disease gradients of wheat stripe rust in large field plots. *Phytopathology* 95, 983–991. doi: 10.1094/PHYTO-95-0983

Sherwood, R. T., and Hagedorn, D. J. (1958). Determining common root rot potential of pea fields. *Wisconsin Agricultural Experimental Station Bulletin* 531: 11p

Sundheim, L. (1972). Physiologic specialization in *Aphanomyces euteiches*. *Physiol. Plant Pathol.* 2, 301–306. doi: 10.1016/0048-4059(72)90013-6

Tivoli, B., Baranger, A., Sivasithamparam, K., and Barbetti, M. J. (2006). Annual *Medicago*: from a model crop challenged by a spectrum of necrotrophic pathogens to a model plant to explore the nature of disease resistance. *Ann. Bot.* 98, 1117–1128. doi: 10.1093/aob/mcl132

Tsvetkova, N., and Kotova, V. (1985). Root rots of peas and vetch. *Zashchita Rastenij* 8, 50–51.

Van der Plank, J. E. (1968). *Disease Resistance in Plants*. New York, NY: Academic Press.

Van Oosterhout, C., Hutchinson, W. F., Wills, D. P., and Shipley, P. (2004). MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. *Mol. Ecol. Notes* 4, 535–538. doi: 10.1111/j.1471-8286.2004.00684.x

Weir, B. S., and Cockerham, C. C. (1984). Estimating F-statistics for the analysis of population structure. *Evolution* 38, 1358–1370. doi: 10.2307/2408641

Whisson, S. C., Drenth, A., Maclean, D. J., and Irwin, J. A. (1994). Evidence for outcrossing in *Phytophthora sojae* and likage of a DNA marker to 2 avirulence genes. *Curr. Genet.* 27, 77–82. doi: 10.1007/BF00326582

Wicker, E. (2001). *Diversité des populations françaises d'Aphanomyces euteiches Drechs., agent de la pourriture racinaire du pois: variabilité pathogène et moléculaire*. Thesis. Ecole Nationale Superieure Agronomique de Rennes.

Wicker, E., Hulle, M., and Rouxel, F. (2001). Pathogenic characteristics of isolates of *Aphanomyces euteiches* from pea in France. *Plant Pathol.* 50, 433–442. doi: 10.1046/j.1365-3059.2001.00590.x

Wicker, E., and Rouxel, F. (2001). Specific behaviour of French Aphanomyces euteiches Drechs. Populations for virulence and aggressiveness on pea, related to isolates from Europe, America and New Zealand. *Eur. J. Plant Pathol.* 107, 919–929. doi: 10.1023/A:1013171217610

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer AG-F and handling Editor declared their shared affiliation.

Copyright © 2018 Quillévéré-Hamard, Le Roy, Moussart, Baranger, Andrivon, Pilet-Nayel and Le May. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.