# Ten Broad Spectrum Resistances to Downy Mildew Physically Mapped on the Sunflower Genome.

Yann Pecrix, Charlotte Penouilh-Suzette, Stephane Munos, Felicity Vear, Laurence Godiard

► To cite this version:

Yann Pecrix, Charlotte Penouilh-Suzette, Stephane Munos, Felicity Vear, Laurence Godiard. Ten Broad Spectrum Resistances to Downy Mildew Physically Mapped on the Sunflower Genome. Frontiers in Plant Science, 2018, 9 (1780), pp.1-15. 10.3389/fpls.2018.01780. hal-01990160

HAL Id: hal-01990160

https://hal.science/hal-01990160

Submitted on 22 Jan 2019

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>

Check for
updates

# Ten Broad Spectrum Resistances to Downy Mildew Physically Mapped on the Sunflower Genome

Yann Pecrix<fnref n="1" />, Charlotte Penouilh-Suzette<fnref n="1" />, Stéphane Muños<fnref n="1" />, Felicity Vear<fnref n="2" />* and Laurence Godiard<fnref n="1" />*

Resistance to downy mildew (*Plasmopara halstedii*) in sunflower (*Helianthus annuus* L.) is conferred by major resistance genes, denoted *Pl*. Twenty-two *Pl* genes have been identified and genetically mapped so far. However, over the past 50 years, wide-scale presence of only a few of them in sunflower crops led to the appearance of new, more virulent pathotypes (races) so it is important for sunflower varieties to carry as wide a range of resistance genes as possible. We analyzed phenotypically 12 novel resistant sources discovered in breeding pools derived from two wild *Helianthus* species and in eight wild *H. annuus* ecotypes. All were effective against at least 16 downy mildew pathotypes. We mapped their resistance genes on the sunflower reference genome of 3,600 Mb, in intervals that varied from 75 Kb to 32 Mb using an AXIOM® genotyping array of 49,449 SNP. Ten probably new genes were identified according to resistance spectrum, map position, hypersensitive response to the transient expression of a *P. halstedii* RXLR effector, or the ecotype/species from which they originated. The resistance source HAS6 was found to carry the first downy mildew resistance gene mapped on chromosome 11, whereas the other resistances were positioned on chromosomes 1, 2, 4, and 13 carrying already published *Pl* genes that we also mapped physically on the same reference genome. The new genes were designated *Pl*23-*Pl*32 according to the current nomenclature. However, since sunflower downy mildew resistance genes have not yet been sequenced, rules for designation are discussed. This is the first large scale physical mapping of both 10 new and 10 already reported downy mildew resistance genes in sunflower.

**Keywords:** Helianthus annuus, Plasmopara halstedii, resistance gene (R gene), SNP array, oomycete, RXLR effector, hypersensitive response (HR), crop protection

## OPEN ACCESS

**Edited by:**
Jacqueline Batley,
University of Western Australia,
Australia

**Reviewed by:**
Leire Molinero-Ruiz,
Instituto de Agricultura Sostenible
(IAS), Spain
Yang Zhang,
St. Jude Children's Research
Hospital, United States

*Correspondence:
Felicity Vear
felicity.vear@wanadoo.fr
Laurence Godiard
laurence.godiard@inra.fr

**Specialty section:**
This article was submitted to
Plant Breeding,
a section of the journal
Frontiers in Plant Science

Received: 17 July 2018
Accepted: 15 November 2018
Published: 04 December 2018

**Citation:**
Pecrix Y, Penouilh-Suzette C, Muños S, Vear F and Godiard L (2018) Ten Broad Spectrum Resistances to Downy Mildew Physically Mapped on the Sunflower Genome. Front. Plant Sci. 9:1780. doi: 10.3389/fpls.2018.01780

## INTRODUCTION

Downy mildew of sunflower, caused by *Plasmopara halstedii* can occur in most countries where sunflowers (*Helianthus annuus* L.) are grown. It is an obligate parasitic oomycete belonging to the Peronosporales showing physiological races or pathotypes (Gascuel et al., 2015) which interact with sunflower genotypes according to the resistance genes they contain. More than 20 pathotypes have been described worldwide (Viranyi et al., 2015) and 22 major resistance genes (denoted *Pl*) have been identified in sunflower (Qi et al., 2016b).

<footnote n="1">¹ Laboratoire des Interactions Plantes Microorganismes, INRA, CNRS, Université de Toulouse, Castanet-Tolosan, France,</footnote>

<footnote n="2">² INRA, Génétique, Diversité, Ecophysiologie des Céréales, UMR 1095, Clermont-Ferrand, France</footnote>
Sunflower was domesticated by Native Americans from wild *H. annuus* but the first large scale developments of the crop were in Russia between 1900 and 1950. *P. halstedii* also originated in North America (Leppik, 1962) but the first important disease attacks occurred in Russia (Pustovoit, 1966) when sunflower became an important crop. In the period 1960–70, G.V. Pustovoit introduced resistance from *H. tuberosus* into cultivated sunflower and bred open pollinated varieties (OPV) Progress and Novinka which were resistant to downy mildew in Russia. When research on sunflower started in Canada, accidental crosses with wild *H. annuus* in 1962–1963 showed that resistance was available in this weedy species (Putt, 1964). In Western Europe, where the climate is favorable for downy mildew, the genes *Pl*₁ and *Pl*₂, from wild *H. annuus*, were rapidly incorporated into hybrid varieties and it became clear that the pathotype prevalent in Europe (100) was not the same as that in the United States (300) (Zimmer and Kinman, 1972). Between 1978 and 1988, nearly all varieties carried these genes and the appearance, first in the United States and then in Europe, of new pathotypes not controlled by these genes (700, 703, 710) (Gulya et al., 1991) may well have been caused by the selection pressure applied on *P. halstedii* as well as by accidental downy mildew introductions in sunflower seed (Ahmed et al., 2012). Progress and Novinka and also some ecotypes of wild *H. annuus* and other annual *Helianthus* spp. (Miller and Gulya, 1988) provided effective resistance. However, in France, rapid appearance of further new pathotypes (304, 314, 334) provided confirmation that this was favored if all widely grown varieties carried the same resistance genes and that it was important to make a search for additional genes (Tourvieille de Labrouhe, 2003).

Studies have been made of quantitative, perhaps pathotype-
non-specific, resistance (Tourvieille de Labrouhe et al., 2008;
Vear et al., 2008b; Vincourt et al., 2012) but this is not easy
to use in breeding, so control of downy mildew in sunflower
still depends largely on major resistance genes. It is thus most
important to find a wide range of genes giving resistance to
the most prevalent pathotypes in areas where sunflower is an
important crop.

Several countries have developed collections of wild *H. annuus*
and other *Helianthus* species, and research for new downy
mildew resistance genes has been made in crosses of this material
with cultivated sunflower. The genes introduced into cultivated
sunflower up to 2008 were summarized by Vear et al. (2008a).
Since then, *Pl*₁₄ (Bachlava et al., 2011) and *Pl*₁₅ (Bertero de
Romano et al., 2010) were identified in Argentina from local OPV
and *Pl*₁₃ (Mulpuri et al., 2009) and *Pl*₁₆ (Liu et al., 2012) in the
United States from similar origins; *Pl*₁₇ and *Pl*₁₉ originated from
wild *H. annuus* (Qi et al., 2015; Zhang et al., 2017) and *Pl*₁₈ and
*Pl*₂₀ were introduced from *H. argophyllus* (Qi et al., 2016a; Ma
et al., 2017). *Pl*₂₁, identified in many restorer genotypes, probably
originated from the same wild *H. annuus* as *Pl*₂ (Vincourt et al.,
2012). *Pl*₂₂, obtained from the OPV Novinka, should be from
*H. tuberosus* but could be from *H. argophyllus* according to
its mapping position on the reference genome (Pecrix et al.,
2018).

To determine whether different resistant accessions contain
different genes, the classic method is to study phenotypic

reactions of progenies from crosses between resistant lines.
However, additional genotypic analyses and fine genetic mapping
are necessary to obtain more precise conclusions. Most of
the downy mildew resistance genes identified since 1972 have
been localized on molecular maps using either RFLP (Pl₁, Pl₂,
Pl₅, Pl₆/7, Pl₈), SSR (simple sequence repeats) (PlARG, Pl₁₃,
Pl₁₅, Pl₁₈, Pl₁₉, Pl₂₁) and more recently SNP (single nucleotide
polymorphisms) (Pl₁₇, Pl₂₀, Pl₂₂). Pl₂₂ was the first resistance
gene to be physically mapped to a region of 1.7 Mb of
chromosome 13 (Pecrix et al., 2018). Some genes have been
difficult to localize; for example Pl₁₃ (Mulpuri et al., 2009) and
Pl₂₂ could not be localized by Bulk Segregant Analyses with RFLP
and SSR in 2000–2005. In addition, several Pl genes are clustered
in complex genomic regions on chromosomes 1, 2, 4, 8, and
13 carrying many resistance gene candidates (RGCs; Radwan
et al., 2003, 2008) and showing duplications (Bachlava et al.,
2011). It is therefore not easy to localize precisely the different
genes relative to each other, especially if different polymorphic
markers are used to map them. Physical mapping of both
known and novel resistance Pl genes on the reference genome
sequence is therefore a major progress and should be an asset for
breeders.

Oomycete pathogens produce virulence protein factors called effectors, secreted and translocated in plants to facilitate colonization (Fawke et al., 2015). Mestre et al. (2016) showed that *P. halstedii* expressed effectors during sunflower infection and overexpression of some cloned effectors induced hypersensitive reactions (HRs) in resistant sunflower lines and not in near-isogenic susceptible lines (Gascuel et al., 2016). Specific recognition of effectors acting as avirulent factors has been used as a reliable marker of resistance genes in potato (Vleeshouwers et al., 2008) and in sunflower (Gascuel et al., 2016; Pecrix et al., 2018). We used this property to distinguish some closely linked resistances.

This paper reports genetic mapping of the broad spectrum
resistance genes in 12 sunflower lines developed from crosses of
cultivated and wild *Helianthus* spp., and their physical mapping
on the sunflower reference genome using an AXIOM® SNP
genotyping array. Ten putatively new downy mildew resistance
genes were located on the sunflower genome sequence together
with 10 already published *Pl* genes mapped on the same
chromosomes. The question of what constitutes a new resistance
gene is discussed.

MATERIALS AND METHODS

Sunflower Genotypes and Crosses Made

In 2005–2006, downy mildew resistance was identified in
interspecific pools and progenies from crosses of cultivated and
wild *H. annuus* developed using the INRA collection. Inbred
lines were developed from progenies showing resistance to both
pathotypes 710 and 304 and were found to give resistance to all
or most French pathotypes (Vear et al., 2008a). The inbred lines
carrying known downy mildew resistance genes and these new
sources of resistance are listed in **Table 1**. HIS or INTER denotes
origins from interspecific pools, HAS denotes hybridization with
<table><caption>**TABLE 1** | Pedigrees of known and new inbred sunflower lines with downy mildew resistance.</caption><thead><tr><th>Resistant line</th><th>INRA code</th><th>Pedigree</th><th>Downy mildew resistance gene</th><th>Source of resistance</th><th>Breeder</th></tr></thead><tbody><tr><td colspan="6"><strong>New downy mildew resistance sources</strong></td></tr><tr><td>HIS33</td><td>UTQ</td><td>INRA interspecific breeding pool H. resinosus 1451</td><td></td><td>H. resinosus 1451</td><td>INRA</td></tr><tr><td>INTER35</td><td>UUQ</td><td>INRA interspecific breeding pool H. resinosus 1323</td><td></td><td>H. resinosus 1323</td><td>INRA</td></tr><tr><td>HAS62</td><td>HAS62</td><td>(VKQ × (OEG × (H. annuus Utah 775RT1A11)))</td><td></td><td>H. annuus Utah 775</td><td>INRA</td></tr><tr><td>HAS40</td><td>HBQ</td><td>((H. annuus Texas 650 × RT1A11) × OEK)</td><td></td><td>H. annuus Texas 650</td><td>INRA</td></tr><tr><td>HAS103</td><td>HAS103</td><td>(OEG × (H. annuus Kansas 999 × 90R19))</td><td></td><td>H. annuus Kansas 999</td><td>INRA</td></tr><tr><td>IDAHO</td><td>IDAHO</td><td>(HA434 × H. annuus Idaho PI468435)F4</td><td></td><td>H. annuus Idaho</td><td>USDA/INRA</td></tr><tr><td>HIS32</td><td>HIS32</td><td>(VKQ × INRA interspecific breeding pool H. tomentosus 1450)</td><td></td><td>H. tomentosus 1450</td><td>INRA</td></tr><tr><td>HIS36</td><td>HIS36</td><td>INRA interspecific breeding pool H. tomentosus</td><td></td><td>H. tomentosus</td><td>INRA</td></tr><tr><td>HAS85</td><td>HEQ</td><td>(OW × (H. annuus Wyoming 963 × RT1A11))</td><td></td><td>H. annuus Wyoming 963</td><td>INRA</td></tr><tr><td>HAS6</td><td>HAQ</td><td>(OW × (H. annuus Wyoming 358 × RT1A11))</td><td></td><td>H. annuus Wyoming 358</td><td>INRA</td></tr><tr><td>HAS42</td><td>HCQ</td><td>((H. annuus Texas 651RT1A11) × OEU)</td><td></td><td>H. annuus Texas 651</td><td>INRA</td></tr><tr><td>HAS54</td><td>HAS54</td><td>(OZ × (H. annuus Oklahoma 662 × 90R19))</td><td></td><td>H. annuus Oklahoma 662</td><td>INRA</td></tr><tr><td colspan="6"><strong>Lines with known downy mildew resistance genes</strong></td></tr><tr><td>PAA1</td><td>PAA1</td><td>(AR22.A.A.B.D.D (from hybrid with H. argophyllus) × PBP1)</td><td>PI<sub>ARG</sub></td><td>AR22 – H. argophyllus</td><td>INRA</td></tr><tr><td>0QP7</td><td>0QP7</td><td>(PAA1 × RHA345)</td><td>PI<sub>ARG</sub></td><td>PAA1 – H. argophyllus</td><td>INRA</td></tr><tr><td>0QP8</td><td>0QP8</td><td>(PAA1 × RHA345)</td><td>PI<sub>ARG</sub></td><td>PAA1 – H. argophyllus</td><td>INRA</td></tr><tr><td>RHA419</td><td>RHA419</td><td>(RHA373 × H. argophyllus 1575-2)</td><td>PI<sub>ARG</sub></td><td>H. argophyllus 1575</td><td>USDA</td></tr><tr><td>HA458</td><td>COQ</td><td>(HA434 × H. annuus Idaho PI468435)</td><td>PI<sub>17</sub></td><td>H. annuus Idaho</td><td>USDA</td></tr><tr><td>VKQ</td><td>VKQ</td><td>(H52 × YDQ (HA335 H. annuus Texas))</td><td>PI<sub>6</sub></td><td>HA335 – H. annuus</td><td>INRA</td></tr><tr><td>XRQ</td><td>XRQ</td><td>(HA89 × Progress 1816)</td><td>PI<sub>5</sub></td><td>Progress – H. tuberosus</td><td>INRA</td></tr><tr><td>HA460</td><td>U8Q</td><td>(HA434 × RHA340)</td><td>PI<sub>8</sub></td><td>RHA340 – H. argophyllus</td><td>USDA</td></tr></tbody></table>

Inbred lines of cultivated sunflower in the pedigrees of new resistance sources were all bred by INRA (Institut National de la Recherche Agronomique, France). USDA, United States Department of Agriculture.

wild *H. annuus*. The line IDAHO was bred from USDA material from which HA458 was also developed (Hulke et al., 2010). It was crossed with PMI3 in 2006 to study both resistance origins. **Table 2** presents resistances of all the new lines to 15 French pathotypes of *P. halstedii* and a Spanish pathotype 330 not present in France.

To determine the number of genes controlling resistance in these new lines and, at the same time whether some are the same (or at least clustered) or whether they are inherited independently, the new lines were crossed firstly with lines representing the resistance genes clusters on chromosomes 1, 8, and 13 together with COQ (HA458) whose gene had not yet been mapped (now located on chromosome 4 (Qi et al., 2015), and then between each other, for most of the possible combinations (**Supplementary Table 1**). With knowledge of the origin of IDAHO, this line was not used to produce test crosses. The F1 hybrids were crossed with downy mildew susceptible lines to obtain test cross progenies which were tested for their downy mildew resistance to determine segregations indicating allelism or independent inheritance.

In addition, to prepare for molecular mapping, the new lines were crossed with lines susceptible to at least some downy mildew pathotypes. F1 and F2 generations were selfed in order to develop F3 progenies which were phenotyped to determine the genotype of the F2 plants from which they were derived (**Table 3**).

## Downy Mildew Resistance Tests

Classic seed immersion downy mildew tests (Mouzeyar et al., 1993) were used to phenotype both test cross and F3 progenies. These tests were made in confined growth chambers in agreement with European quarantine regulations. Pathotypes (710, 304, 334, 714) were among those known to be present in France, and chosen to permit determination of the presence or absence of the resistance genes studied. As in previous studies (Vear et al., 1997), plants with no sporulation (R) or with sporulation only on cotyledons (RII) were counted as resistant; those with sporulation on cotyledons and leaves, whether profuse (S) or rare (SII) were considered as susceptible.

For each test cross progeny, the number of plants tested varied according to the number of seed available. All the results were summed, such that the total number of plants tested with one pathotype for one progeny varied between 91 and 763, except for one progeny with 51 plants, which showed no segregation for resistance (all numbers of plants tested are given in **Supplementary Table 2**). Chi square tests were made for segregation patterns: either no segregation (same gene or closely linked genes), 3R:1S (two independent genes) or 7R:1S (three independent genes).

For each F3 progeny, 20–25 seedlings were observed and, 72 to 177 F3 progenies for each resistance origin were tested using one pathotype and between 41 and 95 when tests were made with a second pathotype (**Table 3**). Chi-square tests were made for segregation patterns for one gene (1RR : 2SEG : 1SS) or for
<table><caption>**TABLE 2** | Reaction of sunflower resistance sources to 15 *Plasmopara halstedii* pathotypes found in France and *P. halstedii* pathotype 330 from Spain.</caption><thead><tr><th rowspan="2">Sunflower genotype</th><th colspan="16">16 P. halstedii pathotypes</th></tr><tr><th>100</th><th>300</th><th>304-10<sup>a</sup></th><th>304-30<sup>a</sup></th><th>307</th><th>314</th><th>330</th><th>334</th><th>700</th><th>703</th><th>704</th><th>707</th><th>710</th><th>714</th><th>730</th><th>774</th></tr></thead><tbody><tr><td>HIS33</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>INTER35</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HAS62</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HAS40</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HAS103</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>Rll</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>Rll</td></tr><tr><td>HAS85</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HIS32</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HIS36</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>IDAHO</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HAS6</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HAS42</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr><tr><td>HAS54</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td><td>R</td></tr></tbody></table>

Plants with no sporulation (R) or with sporulation only on cotyledons (Rll) were considered as resistant; those with sporulation on cotyledons and leaves were noted as susceptible (S). <sup>a</sup>Samples of 304 collected in 2000 and 2002, named according to the proposition of Tourvieille de Labrouhe et al. (2012) using additional differentials.

<table><caption>**TABLE 3** | Downy mildew reactions for F3 progenies from crosses of new resistance sources with lines susceptible to the pathotypes tested.</caption><thead><tr><th>Cross</th><th>Downy mildew pathotype</th><th>Number of F3 progenies</th><th>Results</th><th>Theoretical segregation</th><th>X<sup>2</sup></th><th>Notes</th></tr></thead><tbody><tr><td>HIS33 × OPC3</td><td>710</td><td>101</td><td>24RR : 53SEG : 24SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>0.207ns</td></tr><tr><td>PAS5 × INTER35</td><td>710</td><td>165</td><td>24RR : 96SEG : 45SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>9.65**</td></tr><tr><td>VKQ × HAS62</td><td>334</td><td>90</td><td>44RR : 37SEG : 9SS</td><td>2 genes</td><td>7RR : 8SEG : 1SS</td><td>2.18ns</td></tr><tr><td>VKQ × HAS62</td><td>714</td><td>43</td><td>12RR : 22SEG : 9SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>0.442ns</td></tr><tr><td>OEH × HAS40</td><td>710</td><td>141</td><td>34RR : 74SEG : 33SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>0.096ns</td></tr><tr><td>VKQ+XRQ × HAS103</td><td>334</td><td>87</td><td>30RR : 39SEG : 18SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>4.24ns</td></tr><tr><td>PMI3 × IDAHO</td><td>710</td><td>177</td><td>51RR : 81SEG : 45SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>1.678ns</td></tr><tr><td>VKQ × HIS32</td><td>334</td><td>117</td><td>50RR : 59SEG : 8SS</td><td>2 genes</td><td>7RR : 8SEG : 1SS</td><td>0.096ns</td></tr><tr><td>VKQ × HAS32</td><td>714</td><td>95</td><td>19RR : 60SEG : 16SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>6.77**<sup>a</sup></td></tr><tr><td>VKQHIS36</td><td>334</td><td>72</td><td>17RR : 33SEG : 22SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>1.194ns</td></tr><tr><td>DI × HAS85</td><td>710</td><td>136</td><td>37RR : 70SEG : 29SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>1.96ns</td></tr><tr><td>HAS6 × OPC2</td><td>334</td><td>130</td><td>38RR : 63SEG : 29SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>0.955ns</td></tr><tr><td>HAS42 × PSM7</td><td>710</td><td>148</td><td>33RR : 82SEG : 33SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>1.73ns</td></tr><tr><td>VKQ × HAS54</td><td>334</td><td>95</td><td>35RR : 53SEG : 7SS</td><td>2 genes</td><td>7RR : 8SEG : 1SS</td><td>1.283ns</td></tr><tr><td>VKQ × HAS54</td><td>714</td><td>41</td><td>12RR : 18SEG : 11SS</td><td>1 gene</td><td>1RR : 2SEG : 1SS</td><td>0.659ns</td></tr></tbody></table>

RR, homozygous resistant; SEG, progeny segregating for resistance; SS, homozygous susceptible. <sup>a</sup>If 4 of the 60 progenies noted SEG were in fact SS, the X² would not be significant. **: significantly different P > 0.01; ns: not significantly different.

two genes (7RR : 8SEG : 1SS). Progenies chosen for molecular analyses were preferentially those homozygous for resistance or susceptibility.

Transient Expression Assays of a Cloned *P. halstedii* RXLR Effector

A construct consisting of the *P. halstedii* PhRXLR-C01 RXLR-type effector fused to YFP (Yellow Fluorescent protein) under the control of the 35S promoter was infiltrated in sunflower leaves of entire plants in *Agrobacterium*-mediated transient expression experiments, as previously described (Gascuel et al., 2016). The construct consisting of YFP under the control of the 35S promoter was used as an infiltration negative control. Two leaves of 12 plants were tested in two independent experiments and scored on a 0 to 5 cell death index scale (Gascuel et al.,

2016). The triggering of an hypersensitive response (scores of the infiltrated leaf areas of 4 to 5) was indicative of the recognition of the effector by the resistance gene carried by the plant (Gascuel et al., 2016).

Genotyping Procedures

From the phenotype segregation patterns, the numbers of F3 progenies chosen for genotyping each resistance source varied from 27 to 72, according to apparent similarity with another line or novelty, and to number of seeds available (Supplementary Table 3). Preference was given to progenies showing homozygous resistance or susceptibility. Leaf disks cut from eight seedlings, raised in a greenhouse, for each progeny, were bulked for DNA extraction. Extraction and hybridization experiments were conducted by the Gentyane platform (Plateforme Gentyane,
UMR INRA/UBP 1095 Génétique Diversité et Ecophysiologie des Céréales, 5 chemin de Beaulieu - Clermont-Ferrand, France) on a GeneTitan® (Affymetrix) according to the manufacturer's instructions. Genotyping was performed using a 50K SNP AXIOM® array (Mangin et al., 2017). Marker quality control and allele calling was carried out using the Axiom Analysis Suite software (v1.1.1.66 Affymetrix<fnref n="1" />), following the Best Practices Genotyping Workflow. Samples with a dish quality control (DQC) value ≥ 0.9 and QC call rate ≥ 0.90 thresholds were considered to have passed the quality control assessment. SNP quality control parameters and thresholds were defined as follows: polyploid species-type; cr-cutoff ≥ 90; fld-cutoff ≥ 3.6; het-so-cutoff ≥ -0.3; het-so-otv-cutoff ≥ -0.3; hom-ro-1-cutoff ≥ 0.6; hom-ro-2-cutoff ≥ 0.3; hom-ro-3-cutoff ≥ -0.9. For mapping, only PolyHighResolution SNPs considered to be reliable and informative were retained. A final SNP filter focused on genotype data of parental lines was applied: for either or both parents, markers that failed SNP calling, or were genotyped as heterozygous, or showed divergences of allele calling between technical replicates were discarded from the analysis.

## Mapping Procedures

The genetic maps were constructed using CarthaGène v1.2.3 software (de Givry et al., 2005). Markers were merged with the *mrkmerges* function and linkage groups were obtained with the *group* 0.3 3 function for all mapping populations except for HAS6xOPC2 (*group* 0.5 3). The genetic distances between markers were calculated using the *lkh* 1 -1 function and the genetic maps were obtained with the *bestprintd* function. Published linked markers of *Pl* genes were positioned on the sequenced sunflower genome (INRA inbred line XRQ; Badouin et al., 2017) using BLAST<fnref n="2" /> and the genetic and physical maps were drawn using Mapchart 2.3 software (Voorrips, 2002).

# RESULTS

## Distinction of Downy Mildew Resistance Genes by Phenotypic Studies

**Supplementary Table 2** presents the results of downy mildew resistance tests on all the test crosses, ordered according to the chromosome on which the downy mildew resistance has been mapped (as described below). The X<fnref n="2" /> were calculated for expected segregations in the case of two or three independent genes (3R:1S or 7R:1S), linkage was not considered. Ninety four of the 122 progeny segregations agreed with those expected (non-significant X<fnref n="2" />). In some cases, SII or RII plants could have caused some misinterpretation. Pathotype 710 was used for the majority of tests and gave the highest percentage of rejection of the null hypothesis (significant X<fnref n="2" />): 33% compared with 27% for pathotype 304 and 12% for pathotype 334. Among the lines used in the crosses, COQ (HA458) and HAS40 gave the largest numbers of significant X<fnref n="2" />, but in the case of HAS40, there were more segregations with unexpectedly low numbers of

susceptible plants than those with excess susceptible plants. All X<fnref n="2" /> agreed with the hypotheses proposed for the line HIS32. From these results, the resistances of HIS33, INTER35 and HAS62 showed linkage to *Pl*ARG on chromosome 1, those of IDAHO, HIS32, HIS36 and HAS85 linkage to *Pl*17 on chromosome 4 and those of HAS42 and HAS54 linkage to *Pl*5 on chromosome 13, whereas resistances of HAS103 and HAS6 appeared independent of all the known *Pl* genes tested. It was only for HAS40 that no conclusion could be drawn; it showed few susceptible plants in crosses with *Pl*ARG but also in crosses involving *Pl*17.

**Table 3** presents a summary of the reactions of F3 progenies. Tests were first made with pathotype 710, or pathotype 334 if the “susceptible” parent was VKQ, carrying *Pl*6 conferring resistance to pathotype 710. Three lines, HAS62, HIS32, and HAS54 were found to show segregations indicating the presence of two genes for resistance to pathotype 334 although *Pl*6 is not effective against this pathotype. The three lines are from quite different origins, but they appear to contain a *Pl*2-type gene in addition to that giving resistance to a large number of pathotypes. When tested with pathotype 714, the progenies showed segregations agreeing with the hypothesis for a single dominant gene. These segregations were used to define the resistance genotype of each F2 plant used for mapping.

Transient expression assays with the PhRXLR-C01 effector of *P. halstedii* showed that only three of the 12 resistant lines responded by an HR: HAS103, HAS6 and HAS85, as well as the positive control line PMI3 (Pecrix et al., 2018) (Figure 1). All the other resistant lines showed no HR to the infiltration of the PhRXLR-C01 effector, suggesting that the resistances they carry do not recognize this specific effector. XRQ showed a discoloration phenotype with both YFP negative control and PhRXLR-C01 constructs, suggesting a stress response to the Agroinfiltration process in this ecotype but no specific recognition of the *P. halstedii* effector.

## Physical Mapping of the Resistances With SNP Markers on XRQ Reference Genome

The mapping interval of a resistance gene is delimited by the two closest markers, one on each side, revealing one recombination event between the resistance gene and the marker. This region may contain several markers co-segregating completely with resistance in the progeny tested. **Supplementary Table 3** presents the numbers of polymorphic markers retained per chromosome and those used to map the resistance genes. For seven susceptible × resistant line progenies, a small number of SNP markers (from 0 to 10 for HAS42) co-segregated with resistance which was located in a physical interval of 75 to 4,222 Kb. Four progenies, on chromosomes 1 and 13 (HAS54), showed between 50 and 100 markers co-segregating with resistance, and larger interval sizes from 8.6 to 32 Mb (Supplementary Tables 3, 4). Mapping was not precise for the gene in the line HAS103 since we found 219 SNP markers of chromosome 2 (on a total of 345) co-segregating with resistance, and a physical interval for the resistance of 114 Mb.

<footnote n="1">¹http://www.affymetrix.com</footnote>

<footnote n="2">²https://www.heliagene.org/</footnote>
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Recognition of *Plasmopara halstedii* PhRXLR-C01 effector in resistant sunflower genotypes by HR-like responses. Phenotypes observed in sunflower leaves of the 12 new resistant genotypes, after *Agrobacterium*-transient expression of the *P. halstedii* effector PhRXLR-C01 fused to YFP, or the YFP control, both driven by the 35S promoter. PMI3 (*Pl*22) and XRQ (*Pl*5) were used as positive and negative control lines to PhRXLR-C01 HR response, respectively. Pictures were taken 11 days post-infiltration (for more details see Gascuel et al., 2016). Sunflower leaf size varied from 5 to 6 cm long and 2 to 3 cm wide, depending on the tested line.</figcaption></figure>

## Four New Resistances Map to Chromosome 1

The resistances in HIS33, INTER35, HAS62, and HAS40 were mapped on chromosome 1 (Figure 2A). The first three co-localize with markers for *Pl*ARG (Wieckhorst et al., 2010) in an area representing 20 to 30 Mb on the reference genome sequence and all overlap (Figure 2B). In contrast, the resistance of HAS40 appears situated in the region carrying *Pl*13, *Pl*14, and *Pl*16 (Liu et al., 2012) but differs from these genes since it is effective against more pathotypes (Figure 2C). Comparison with the *Pl* genes previously mapped was rather difficult as there are differences in marker and gene orders in the genotypes that were used for the crosses. The gene in HAS40 mapped above the SNP marker HT636 and so appears different from *Pl*13 and *Pl*16, which were mapped below this marker (Mulpuri et al., 2009; Liu et al., 2012) (Figure 2C). *Pl*14 also mapped above this marker but the RGCs identified in this region by Bachlava et al. (2011) could not be used in the present study because they mapped by blastn at different locations including on different chromosomes of the reference genome, with high e-values (,<10<sup>-100</sup>). The resistance of HAS40 was quite precisely mapped on the genome sequence in a region of 75 Kb, but surprisingly corresponding to a large genetic region of 16.8 cM. This genomic region does not appear to include any RGC in XRQ, where there may have been a deletion.

## One Resistance Maps to Chromosome 2

The resistance of HAS103, a wild *H. annuus* from Kansas, mapped on chromosome 2 (Figure 3). Although on the genetic map the gene appears to be located in only 3.5 cM, on the physical

map it is somewhere in 114 Mb. This appears to be related to a lack of recombination on a large part of this chromosome. It is therefore not possible to conclude whether this gene colocalizes with markers for *Pl*18 also on chromosome 2 but coming from *H. argophyllus* (Qi et al., 2016a). HAS103 showed HR with PhRXLR-C01 (Figure 1), its resistance would be distinguished from that of *Pl*18, if lines containing this gene did not respond to the effector.

## Four Resistances Map to Chromosome 4

The resistances in IDAHO, HIS32, HIS36, and HAS85 mapped in a small region at one end of chromosome 4 (Figure 4A). While those of HIS32 and IDAHO overlap *Pl*17, those in HIS36 and HAS85 mapped in a distinct region with markers reported to be linked to *Pl*19 (Figure 4B). The resistance of IDAHO, derived from the same *H. annuus* origin as HA458 carrying *Pl*17 and conferring resistance to the same 16 pathotypes, mapped precisely in 1 Mb in part of the area of this gene, so we concluded that it was controlled by the same gene (Qi et al., 2015). In contrast since the resistances of HIS32 and HIS36, both from the *H. tomentosus* breeding pool, mapped to non-overlapping genomic regions of the reference genome, they must, at present, be considered as controlled by different genes. On the physical map they are placed in regions covering *Pl*17 and *Pl*19 respectively but also some distances outside the area of these genes (Figure 4B). The resistance of HAS85 was mapped very close to markers for *Pl*19 (Zhang et al., 2017) and it differs from the three others, IDAHO, HIS32 and HIS36, since it is the
<figure><img src="image_4.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>C</figcaption></figure>

**FIGURE 2** | Four downy mildew resistances map to chromosome 1 of sunflower. (A) Genetic (yellow) and physical (gray) maps of the *PI23_HIS33*, *PI23_INTER35*, *PI24_HAS62*, and *PI25_HAS40* resistance genes on chromosome 1. The cross used to map the resistance is indicated on top of each map, with the resistant parent in bold. Only the non-redundant SNP markers are indicated with their genetic (left map) and corresponding physical positions on the XRG genome (right map).

(B) Enlarged view of the physical map showing the localization of *PI23_HIS33*, *PI23_INTER35*, and *PI24_HAS62* genes. The chromosomal region where *PIARG* was located (Wieckhorst et al., 2010) is indicated in purple on the physical map. The polymorphic SNP markers of the sunflower AXIOM® array denominated AX-xxxxxxx are indicated with their physical positions on the sunflower genomic sequence. (C) Enlarged view of the physical map showing the localization of *PI25_HAS40* gene. *PI14*, *PI13*, and *PI16* resistance genes were localized according to published linked genetic markers (Mulpuri et al., 2009; Bachlava et al., 2011; Liu et al., 2012).
<figure><img src="image_6.png" /></figure>

**FIGURE 3** | One downy mildew resistance maps to chromosome 2 of sunflower. Genetic (yellow) and physical (gray) map of the *Pl*26 HAS103 resistance gene on chromosome 2. The cross used to map the resistance is indicated on top of the map, with the resistant parent in bold. Only the non-redundant SNP markers are indicated with their genetic (left map) and corresponding physical positions on the XRQ genome (right map).

only one of the four to respond by a hypersensitive response to the *P. halstedii* RXLR effector PhRXLR-C01 transiently over-expressed (Figure 1 and Table 4).

## One New Resistance Maps to Chromosome 11

The downy mildew resistance of HAS6 was the first to be mapped on chromosome 11, in a short region spanning 1.3 Mb (Figure 5). This resistance also shows HR with PhRXLR-C01 (Figure 1 and Table 4).

## Two Resistances Map to Chromosome 13

The resistances of HAS42 and HAS54 both mapped in a small area on chromosome 13 where *Pl*5, *Pl*8, *Pl*21, and *Pl*22 had already been mapped, closest to *Pl*8 (Figure 6). However, they differ from all these genes, either by their reaction to certain downy mildew pathotypes (*Pl*5, *Pl*21, and *Pl*22) or by their resistance

phenotypes (*Pl*8 always shows RII type reactions whereas they show R) (Gascuel et al., 2015).

We positioned these *Pl* genes on the physical map of chromosome 13 according to their linked published markers. Radwan et al. (2003), located *Pl*5 2 cM and *Pl*8 4.5 cM from marker NTIR7H3 using two different populations. However, recombination in this area appears to vary between crosses, Qi et al. (2017) reported 5 cM between, ORS316 and *Pl*8 in one cross and 1.2 cM in another. Taking into consideration mapping of *Pl*5 and *Pl*21 below ORS316 but above ORS630 by Vincourt et al. (2012), we propose the map order in **Figure 6B**. *Pl*22 is above ORS316, with *Pl*5 and *Pl*21 between the two SSR markers and then *Pl*8, and the resistance genes in HAS42 and HAS54 in the same zone, covering 8 Mb (**Figure 6B**). Bachlava et al. (2011) reported at least 14 RGC in the region of *Pl*8, so it is perhaps not surprising that many *Pl* genes map in this area.

Using the AXIOM® genotyping array, the physical mapping of 11 of the 12 resistance genes was consistent with the conclusions of the allelism tests performed initially and the only resistance (HAS40) that could not be localized by phenotypic studies was mapped successfully with the array on chromosome 1.

## DISCUSSION

The results presented here describe mapping, on the sunflower reference genome (Badouin et al., 2017), using a high-throughput mapping SNP AXIOM® array, of 12 sources of resistance effective against at least 16 downy mildew pathotypes. The results suggest that these sources may be considered as providing new resistances carried by 10 novel resistance genes.

### Phenotyping

Before the availability of molecular analyses, determination of linkage or independence of downy mildew resistance genes depended on resistance tests on seedlings with adequate pathotypes. The present results show that a large number of crosses and progenies must be tested for analyses of segregations to permit decisions as to whether *Pl* genes are in one or different clusters. Seventy seven percent of tests reported here did not differ from the Mendelian segregations 3R: 1S, 7R:1S or no segregation. Of those with significant X², some were probably due to the wide use of the highly aggressive pathotype 710, which, in breeding programs, has occasionally been found to give symptoms in plants containing resistance genes but with background genotypes that are unfavorable for expression of resistance. Concerning parental lines, COQ (HA458) and HAS40 showed the highest proportion of “abnormal” segregations (6/13 and 8/17 respectively). For COQ, there was an excess of susceptible plants in five tests out of the six, and it may be suggested that this is due to background susceptibility to pathotype 710. In contrast, for HAS40, five of the eight “abnormal” segregations showed an excess of resistant plants. In three cases they were in crosses with lines (PAA1, HIS33, and HAS62) that were confirmed by the SNP genotyping to have resistance genes on the same region of chromosome 1. These segregations thus probably corresponded to loose linkage
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Three downy mildew resistances map to chromosome 4 of sunflower, close to *PI*17. (A) Genetic (yellow) and physical (gray) maps of the *PI*17, *PI*27\_HIS32, *PI*28\_HIS36, and *PI*29\_HAS85, resistance genes on chromosome 4. The cross used to map the resistance is indicated on top of each map, with the resistant parent in bold. Only the non-redundant SNP markers are indicated with their genetic (left map) and corresponding physical positions on the XRO genome (right map).
(B) Enlarged view of the physical map showing the localization of *PI*27\_HIS32, *PI*28\_HIS36, and *PI*29\_HAS85 genes. *PI*17 and *PI*19 resistance genes were localized according to published linked genetic markers (Qi et al., 2015; Zhang et al., 2017). The polymorphic SNP markers of the sunflower AXIOM® array denominated AX-xxxxxxx are indicated with their physical positions on the sunflower XRO genomic sequence.</figcaption></figure>
<table><caption>**TABLE 4** | Hypersensitive reactions (HR) with *P. halstedii* PhRXLR-C01 effector and resistance gene names.</caption><thead><tr><th>Chromosome</th><th>Sunflower genotype</th><th>HR with <i>P. halstedii</i> PhRXLR-C01 effector</th><th>Origin difference from co-localized genes</th><th>PI gene name</th></tr></thead><tbody><tr><td>1</td><td>HIS33</td><td>No HR</td><td>From <i>H. resinosus</i>, PI<sub>ARG</sub> from <i>H. argophyllus</i></td><td>PI<sub>23_HIS33</sub></td></tr><tr><td>1</td><td>INTER35</td><td>No HR</td><td>From <i>H. resinosus</i>, PI<sub>ARG</sub> from <i>H. argophyllus</i></td><td>PI<sub>23_INTER35</sub></td></tr><tr><td>1</td><td>HAS62</td><td>No HR</td><td>From wild <i>H. annuus</i> Utah, PI<sub>ARG</sub> from <i>H. argophyllus</i></td><td>PI<sub>24_HAS62</sub></td></tr><tr><td>1</td><td>HAS40</td><td>No HR</td><td>From wild <i>H. annuus</i> Texas, PI<sub>13</sub>, PI<sub>14</sub>, PI<sub>16</sub> from cultivated sunflower, HAR4, HAR5</td><td>PI<sub>25_HAS40</sub></td></tr><tr><td>2</td><td>HAS103</td><td>HR</td><td>From wild <i>H. annuus</i> Kansas PI<sub>18</sub> from <i>H. argophyllus</i></td><td>PI<sub>26_HAS103</sub></td></tr><tr><td>4</td><td>IDAHO</td><td>No HR</td><td></td><td>PI<sub>17</sub></td></tr><tr><td>4</td><td>HIS32</td><td>No HR</td><td>From <i>H. tomentosus</i>, PI<sub>17</sub> from wild <i>H. annuus</i> Idaho</td><td>PI<sub>27_HIS32</sub></td></tr><tr><td>4</td><td>HIS36</td><td>No HR</td><td>From <i>H. tomentosus</i>, PI<sub>19</sub> from wild <i>H. annuus</i> Texas</td><td>PI<sub>28_HAS36</sub></td></tr><tr><td>4</td><td>HAS85</td><td>HR</td><td>From wild <i>H. annuus</i> Wyoming, PI<sub>19</sub> from wild <i>H. annuus</i> Texas</td><td>PI<sub>29_HAS85</sub></td></tr><tr><td>11</td><td>HAS6</td><td>HR</td><td></td><td>PI<sub>30_HAS6</sub></td></tr><tr><td>13</td><td>HAS42</td><td>No HR</td><td>From wild <i>H. annuus</i> Texas</td><td>PI<sub>31_HAS42</sub></td></tr><tr><td>13</td><td>HAS54</td><td>No HR</td><td>From wild <i>H. annuus</i> Oklahoma</td><td>PI<sub>32_HAS54</sub></td></tr></tbody></table>

between the genes, which could not be estimated precisely with the numbers of plants phenotyped for their downy mildew resistance. It was concluded that phenotyping by resistance tests provided first indications of the genetics of the downy mildew resistances studied but that molecular studies were necessary to confirm distinction and precise location.

## Mapping With SNPs on a Reference Sunflower Genome

The AXIOM® array used for genotyping has theoretically 49,449 SNP. Over the 12 crosses, we found an average of 5953 high-quality polymorphic SNPs, between a minimum of 4138 for HAS40 and a maximum of 8548 for HIS33 (i.e., 17% of the total SNPs present on the array). Although some physical information was lost according to the cross, there were always sufficient markers to map correctly the resistances in regions of 10s of Mb, and sometimes in 10s of kb (HAS40). This was quite surprising since HAS40 had a small number of polymorphic markers on the chromosome concerned (221). However, the F3 population used for mapping was large (70 plants) making it possible to take account of a sufficient number of recombination events. In addition, the segregation ratios observed for resistance were close to theoretical values, suggesting that mapping population size and phenotyping quality are important criteria for precise fine mapping. On the contrary, the cross involving HAS103 had 50% more polymorphic markers on the chromosome concerned, but the mapping population was smaller (45) and segregation for resistance not optimal. Together with a lack of recombination events on chromosome 2 in the HAS103 x XRQ population, the result is poor physical mapping.

**Supplementary Figure 1**, plotting size of intervals in which resistances were mapped against numbers of markers polymorphic between parental lines on the chromosome concerned, suggests that there could be a relation between polymorphism, a measure of distance between lines, and frequency of recombination. In other words, high levels of polymorphism would disfavor recombination events at the DNA level in segregating progenies, rendering fine mapping

more difficult. High levels of polymorphism have already been observed when the parental lines of crosses used to map resistance were of very different origins, with one parent from wild *H. annuus* ecotypes or from different *Helianthus* species. In the first report mapping *PI<sub>ARG</sub>* chromosome 1 appeared much shorter in a cross involving *H. argophyllus* and cultivated sunflower (Dussle et al., 2004) than in crosses between two cultivated sunflower lines (Wieckhorst et al., 2010). Moreover, it has recently been shown in rice that regions with high levels of polymorphism, in some cases associated with clusters of resistance genes, display reduced crossover frequencies (Mieulet et al., 2018); illustrating a negative correlation between genome divergence of the parental lines and recombination rate.

In addition, this work permits, for the first time, comparisons between genetic and physical maps in sunflower, for 5 of the total 17 chromosomes and highlights several interesting features. First, we observed the presence of differences between the two maps, particularly concerning chromosomes 1, 11, and 13. On the physical map of the reference genome, there appear to be some large genetic regions (10s of cM) devoid of SNPs, whereas the SNPs selected for the study provide relatively regular chromosome coverage. Thus, in the populations studied, the recombination frequency may be high in these regions compared with the reference genome. Alternatively, there may have been deletions in the reference genome compared with the genotypes used for mapping. Second, we observed inversions between genetic and physical maps, for example on chromosome 4, again suggesting differences between the parental genotypes used to map the resistances and the reference genome. It would be interesting to have the sequence of more than one genome to allow better comparisons with genotypes of various origins and to understand sunflower genome evolution.

## Physical Mapping of New and Published PI Resistances

Mapping and positioning physically on the reference genome made it possible to compare localizations of new downy mildew resistance genes from different genetic backgrounds and also
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | The downy mildew resistance of HAS6 is the first to be localized on chromosome 11 of sunflower. Genetic (yellow) and physical (gray) map of the *Pl30_HAS6* resistance gene on chromosome 11. The cross used to map the resistance is indicated on top of the map, with the resistant parent in bold. Only the non-redundant SNP markers are indicated with their genetic (left map) and corresponding physical positions on the XRQ genome (right map).</figcaption></figure>

*Pl* genes that were previously genetically mapped in the same regions. The genetic regions defined for formerly mapped *Pl* genes appear to be very large on the physical map, encompassing up to 10 Mb for *PlARG* and up to 20 Mb for *Pl14*, both on chromosome 1. Badouin et al. (2017) estimated that 1 Mb on

the XRQ sequenced genome contains an average of 20 genes. It is therefore highly probable to find closely linked but different resistance genes in such large regions, and one of the reasons why we propose two new genes, close to *PlARG* on chromosome 1. However, it should be kept in mind that gene order in some genotypes could be slightly different from that in the reference genotype, making fine mapping difficult.

## Gene Distinction and Nomenclature

Even with the availability of a sunflower genome sequence, all the questions regarding the genetic and molecular architecture of resistance to downy mildew will not be solved in the near future, so it appears of interest to fix a rule to designate genes in sunflower involved in the interaction with *P. halstedii*.

In the past, if two resistance sources had resistance to the same downy mildew pathotypes and showed no segregation in test crosses, they were considered to carry the same gene. For example two of the pathotype differential lines (Gascuel et al., 2015), PM17 (D5) and 803.1 (D6), showed no segregation when crossed with XRQ, carrying *Pl5*, but they showed slight differences in phenotyping reactions to pathotypes: PM17 is susceptible to a Spanish isolate of pathotype 330 to which XRQ is resistant, so it was considered to have *Pl5* whereas 803.1 is resistant to the US isolate of 330 to which XRQ is susceptible so its gene was defined as *Pl5* (Vear et al., 2008a). If they had been mapped recently, they would probably been given different names. To distinguish resistance sources with closely linked genes by segregation studies would require very large populations and certainty that observations of one or a few susceptible plants were correct. At present, there are many sources that cannot be distinguished by their pathotype reaction, because they are effective against all known pathotypes. Some of them are located in the same genomic regions, but since slightly different markers are polymorphic according to the lines studied, it is difficult to be certain that they are in exactly the same place and to decide whether they are different genes. An example involved in mapping of the resistance of HAS40 concerned *Pl13*, *Pl14* and *Pl16* that gave initially identical reactions to downy mildew pathotypes. These genes were identified in rust resistant pools developed by USDA from Argentinean origins: *Pl13* in HAR5 (Mulpuri et al., 2009), *Pl14* in a line derived from HAR4 (Bachlava et al., 2011) and *Pl16* in HAR4 (Liu et al., 2012). All three genes were mapped to the lower end of chromosome 1, but, surprisingly, and depending on the position of one marker HT636, the two genes from HAR4, *Pl14* and *Pl16* appeared to be at different positions, the latter being mapped in the same region as *Pl13* according to our results. However, recent results confirm that *Pl13* differs from *Pl14* and *Pl16* because pathotype 705, differentiating these resistances, was reported by Garcia-Carneros and Molinero-Ruiz (2017) and all the genes in this cluster may be slightly different, comparable with those on chromosome 13.

Recent studies of presence or absence of host reactions to pathogen effectors provide new evidence that some of the genes situated in the same region are different (Pecrix et al., 2018). In the case of chromosome 1, none of the four sources responded by an HR to the *P. halstedii* effector PhRXLR-C01. However, this test
$Pl_{23}$ and $Pl_{24}$ show that this is not the case. The different clusters could in part be duplications as sunflower is considered to have undergone chromosome doubling (Badouin et al., 2017) but the many genes in one cluster could well be sequences that have been imperfectly duplicated or the result of unequal crossing-overs as has been shown for other major resistance genes (Meziadi et al., 2016).

We performed the colinearity studies of the genetic maps using the genome sequence of the cultivated sunflower line XRQ (Badouin et al., 2017). However, although this resource was a major breakthrough for genetic studies in sunflower such as for the mapping of the resistance genes to downy mildew of our study, this reference sequence does not represent the diversity in cultivated sunflower and even less the diversity in wild *Helianthus*. Plant genomes vary considerably, including within the same species. It is now relatively easy to re-sequence genomes and identify SNP polymorphisms by comparing the re-sequencing data with a reference sequence. However, the structural variations (inversion, translocation, deletion, insertion, or copy number variation) that can range from a few nucleotides to 100s of kb, are more difficult to identify. By aligning genetic maps on the reference genome sequence, we identified such polymorphisms which suggest that some resistances, originating from wild sunflowers, could be due to structural variations. In order to facilitate the sequencing of the downy mildew resistance genes or other loci of interest, it would be very useful to produce new tools or new genomic resources in sunflower. For example, from a pan-genome, it would be interesting to develop a genotyping tool based on SNPs identified on this pan-genome in order to identify large structural variations. At a finer scale, making the optical maps of the genomes from 10s of accessions would be a useful resource too. Finally, on an even finer scale, assembling high quality genome sequences of sunflowers and wild *Helianthus* would considerably speed up research in sunflower, which could become a model plant species in many scientific fields (genetics, evolution, floral development, responses to biotic and abiotic stresses...).

The present results suggest that *Helianthus* species, and more particularly wild *H. annuus* are quite rich in broad spectrum downy mildew resistance genes. Use of combinations of those described here should help to provide some reasonably durable resistance, protecting the sunflower crop against changes in downy mildew pathotypes.

## ETHICS STATEMENT

Sunflower plant material used in this study was grown from seeds made available by INRA as part of INRA sunflower genetic

resources. The lines Idaho, COQ and U8Q were obtained by selfing from USDA sources.

## AUTHOR CONTRIBUTIONS

YP, FV, and LG conceived the project. FV produced the plant material and carried out phenotyping. CP-S and YP produced the plant material for genotyping. SM designed the 50K-SNP AXIOM® array. YP and LG acquired genotyping data. YP mapped the resistance genes. FV, LG, SM, and YP wrote the manuscript. All authors were involved in revising the manuscript critically and approved it.

## FUNDING

This work was supported by the French Laboratory of Excellence project 'TULIP' (Toward a unified theory of biotic interactions: role of environmental perturbations) (ANR-10-LABX-41 and ANR-11-IDEX-0002-02) and by PROMOSOL ("Association pour la promotion de la sélection des plantes oléagineuses"), project "Milvarsunres."

## ACKNOWLEDGMENTS

The authors thank CRB sunflower for providing seeds (LIPM, Castanet-Tolosan, France), Marie-Claude Boniface, Nicolas Blanchet, Nicolas Pouilly, Anaïs Maistre (LIPM, Castanet-Tolosan, France), Sylvie Roche and Frédéric Serre (INRA, UMR1095, Clermont-Ferrand, France) for technical help in plant experiments and Clémence Paris (Syngenta seeds, Toulouse, France) for her help in genotyping analyses. The authors also thank the Gentyane Platform (Clermont-Ferrand, France) for genotyping experiments. They acknowledge Denis Tourvieille de Labrouhe (INRA Clermont-Ferrand, France) for *P. halstedii* pathotypes, Joël Piquemal (Syngenta seeds, Toulouse, France) and Emmanuelle Mestries (Terres Inovia, France) for being members of the "Milvarsunres" project steering committee, and Patrick Vincourt for critical reading of the manuscript.

## SUPPLEMENTARY MATERIAL

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

## REFERENCES

Ahmed, S., de Labrouhe, D. T., and Delmotte, F. (2012). Emerging virulence arising from hybridisation facilitated by multiple introductions of the sunflower downy mildew pathogen *Plasmopara halstedii*. *Fungal Genet. Biol.* 49, 847–855. doi: 10.1016/j.fgb.2012.06.012

Bachlava, E., Radwan, O. E., Abratti, G., Tang, S. X., Gao, W. X., Heesacker, A. F., et al. (2011). Downy mildew (Pl (8) and Pl (14) and rust (R (Adv) ) resistance genes reside in close proximity to tandemly duplicated clusters of non-TIR-like NBS-LRR-encoding genes on sunflower chromosomes 1 and 13. *Theor. Appl. Genet.* 122, 1211–1221. doi: 10.1007/s00122-010-1525-0
Badouin, H., Gouzy, J., Grassa, C. J., Murat, F., Staton, S. E., Cottret, L., et al. (2017). The sunflower genome provides insights into oil metabolism, flowering and asterid evolution. *Nature* 546, 148–152. doi: 10.1038/nature22380

Bertero de Romano, A., Romano, C., Bulos, M., Altieri, E., and Sala, C. (2010). "A new gene for resistance to downy mildew in sunflower," in *Proceedings of the International Symposium on Sunflower Breeding for Resistance to Disease* (Krasnodar: ISA), 141–146.

de Givry, S., Bouchez, M., Chabrier, P., Milan, D., and Schiex, T. (2005). Carthagene: multipopulation integrated genetic and radiation hybrid mapping. *Bioinformatics* 21, 1703–1704. doi: 10.1093/bioinformatics/bti222

Dimitrijevic, A., and Horn, R. (2018). Sunflower hybrid breeding: from markers to genomic selection. *Front. Plant Sci.* 8:20. doi: 10.3389/fpls.2017.02238

Dussle, C. M., Hahn, V., Knapp, S. J., and Bauer, E. (2004). PlArg from *Helianthus argophyllus* is unlinked to other known downy mildew resistance genes in sunflower. *Theor. Appl. Genet.* 109, 1083–1086. doi: 10.1007/s00122-004-1722-9

Fawke, S., Doumane, M., and Schornack, S. (2015). Oomycete interactions with plants: infection strategies and resistance principles. *Microbiol. Mol. Biol. Rev.* 79, 263–280. doi: 10.1128/mmbr.00010-15

Garcia-Carneros, A. B., and Molinero-Ruiz, L. (2017). First report of the highly virulent race 705 of *Plasmopara halstedii* (Downy Mildew of Sunflower) in portugal and in spain. *Plant Dis.* 101, 1555–1555. doi: 10.1094/pdis-02-17-0289-pdn

Gascuel, Q., Buendia, L., Pecrix, Y., Blanchet, N., Muños, S., Vear, F., et al. (2016). RXLR and CRN effectors from the sunflower downy mildew pathogen *Plasmopara halstedii* induce hypersensitive-like responses in resistant sunflower lines. *Front. Plant Sci.* 7:1887. doi: 10.3389/fpls.2016.01887

Gascuel, Q., Martinez, Y., Boniface, M. C., Vear, F., Pichon, M., and Godiard, L. (2015). The sunflower downy mildew pathogen *Plasmopara halstedii*. *Mol. Plant Pathol.* 16, 109–122. doi: 10.1111/mpp.12164

Gulya, T. J., Sackston, W. E., Viranyi, F., Masirevic, S., and Rashid, K. Y. (1991). New races of the sunflower downy mildew pathogen (*Plasmopara halstedii*) in Europe and North and South-America. *J. Phytopathol.* 132, 303–311. doi: 10.1111/j.1439-0434.1991.tb00125.x

Gulya, T. J., Tourvieille de Labrouhe, D., Masirevic, S., Penaud, A., Rashid, K., and Viranyi, F. (1998). "Proposal for standardized nomenclature and identifaction of races of *Plasmopara halstedii* (Sunflower downy mildew)," in *Proceedings of the Symposium on Sunflower Downy Mildew* (Fargo, ND: ISA), 130–136.

Hulke, B. S., Miller, J. F., Gulya, T. J., and Vick, B. A. (2010). Registration of the oilseed sunflower genetic stocks HA458, HA459, and HA460 possessing genes for resistance to Downy Mildew. *J. Plant Regist.* 4, 93–97. doi: 10.3198/jpr2009.08.0426crgs

Leppik, E. E. (1962). Distribution of downy mildew and some other seed-borne pathogens on sunflower. *FAO Plant Prot. Bull.* 10, 126–129.

Liu, Z., Gulya, T. J., Seiler, G. J., Vick, B. A., and Jan, C. C. (2012). Molecular mapping of the Pl(16) downy mildew resistance gene from HA-R4 to facilitate marker-assisted selection in sunflower. *Theor. Appl. Genet.* 125, 121–131. doi: 10.1007/s00122-012-1820-z

Ma, G. J., Markell, S. G., Song, Q. J., and Qi, L. L. (2017). Genotyping-by-sequencing targeting of a novel downy mildew resistance gene Pl (20) from wild *Helianthus argophyllus* for sunflower (*Helianthus annuus* L.). *Theor. Appl. Genet.* 130, 1519–1529. doi: 10.1007/s00122-017-2906-4

Mangin, B., Pouilly, N., Boniface, M. C., Langlade, N. B., Vincourt, P., Vear, F., et al. (2017). Molecular diversity of sunflower populations maintained as genetic resources is affected by multiplication processes and breeding for major traits. *Theor. Appl. Genet.* 130, 1099–1112. doi: 10.1007/s00122-017-2872-x

Mestre, P., Carrere, S., Gouzy, J., Piron, M. C., de Labrouhe, D. T., Vincourt, P., et al. (2016). Comparative analysis of expressed CRN and RXLR effectors from two Plasmopara species causing grapevine and sunflower downy mildew. *Plant Pathol.* 65, 767–781. doi: 10.1111/ppa.12469

Meziadi, C., Richard, M. M. S., Derquennes, A., Thareau, V., Blanchet, S., Gratias, A., et al. (2016). Development of molecular markers linked to disease resistance genes in common bean based on whole genome sequence. *Plant Sci.* 242, 351–357. doi: 10.1016/j.plantsci.2015.09.006

Michelmore, R., and Wong, J. (2008). Classical and molecular genetics of Bremia lactucae, cause of lettuce downy mildew. *Eur. J. Plant Pathol.* 122, 19–30. doi: 10.1007/s10658-008-9305-2

Mieulet, D., Aubert, G., Bres, C., Klein, A., Droc, G., Vieille, E., et al. (2018). Unleashing meiotic crossovers in crops. *bioRxiv* [Preprint]. doi: 10.1101/343509

Miller, J. F., and Gulya, T. J. (1988). Registration of six downy mildew resistant sunflower germplasm lines. *Crop Sci.* 28, 1040–1041. doi: 10.2135/cropsci1988.0011183X002800060073x

Mouzeyar, S., Tourvieille de Labrouhe, D., and Vear, F. (1993). Histopathological studies of resistance of sunflower (*Helianthus annuus* L.) to downy mildew (Plasmopara halstedii). *J. Phytopathol.* 139, 289–297. doi: 10.1111/j.1439-0434.1993.tb01430.x

Mulpuri, S., Liu, Z., Feng, J., Gulya, T. J., and Jan, C. C. (2009). Inheritance and molecular mapping of a downy mildew resistance gene, Pl (13) in cultivated sunflower (*Helianthus annuus* L.). *Theor. Appl. Genet.* 119, 795–803. doi: 10.1007/s00122-009-1089-z

Pecrix, Y., Buendia, L., Penouilh-Suzette, C., Maréchaux, M., Legrand, L., Bouchez, O., et al. (2018). Sunflower resistance to various downy mildew pathotypes revealed by recognition of conserved effectors of the oomycete *Plasmopara halstedii*. *Plant J.* doi: 10.1111/tpj.14157

Pustovoit, G. V. (1966). "Distant (Interspecific) hybridization of sunflowers in the USSR," in *Proceedings of the 2nd International Sunflower Conference* (Morden: ISA).

Putt, E. D. (1964). "Sunflower breeding in Canada," in *Proceedings of the 1st International Sunflower Conference* (College Station, TX: ISA).

Qi, L. L., Foley, M. E., Cai, X. W., and Gulya, T. J. (2016a). Genetics and mapping of a novel downy mildew resistance gene, Pl (18), introgressed from wild *Helianthus argophyllus* into cultivated sunflower (*Helianthus annuus* L.). *Theor. Appl. Genet.* 129, 741–752. doi: 10.1007/s00122-015-2662-2

Qi, L. L., Ma, G., Talukder, Z. I., Seiler, G. J., Hulke, B. S., Jan, C. C., et al. (2016b). "Molecular mapping of the disease resistance gene and its impact on sunflower breeding," in *Proceedings of the 19th Intetnational Sunflower Conference*, Edirne, 20–30.

Qi, L. L., Long, Y. M., Jan, C. C., Ma, G. J., and Gulya, T. J. (2015). Pl(17) is a novel gene independent of known downy mildew resistance genes in the cultivated sunflower (*Helianthus annuus* L.). *Theor. Appl. Genet.* 128, 757–767. doi: 10.1007/s00122-015-2470-8

Qi, L. L., Talukder, Z. I., Hulke, B. S., and Foley, M. E. (2017). Development and dissection of diagnostic SNP markers for the downy mildew resistance genes Pl (Arg) and Pl (8) and marker-assisted gene pyramiding in sunflower (*Helianthus annuus* L.). *Mol. Genet. Genomics* 292, 551–563. doi: 10.1007/s00438-017-1290-8

Radwan, O., Bouzidi, M. F., Nicolas, P., and Mouzeyar, S. (2004). Development of PCR markers for the Pl5/Pl8 locus for resistance to *Plasmopara halstedii* in sunflower, *Helianthus annuus* L. from complete CC-NBS-LRR sequences. *Theor. Appl. Genet.* 109, 176–185. doi: 10.1007/s00122-004-1613-0

Radwan, O., Bouzidi, M. F., Vear, F., Philippon, J., de Labrouhe, D. T., Nicolas, P., et al. (2003). Identification of non-TIR-NBS-LRR markers for the Pl5/Pl8 locus for resistance to downy mildew in sunflower. *Theor. Appl. Genet.* 106, 1438–1446. doi: 10.1007/s00122-003-1196-1

Radwan, O., Gandhi, S., Heesacker, A., Whitaker, B., Taylor, C., Plocik, A., et al. (2008). Genetic diversity and genomic distribution of homologs encoding NBS-LRR disease resistance proteins in sunflower. *Mol. Genet. Genomics* 280, 111–125. doi: 10.1007/s00438-008-0346-1

Tourvieille de Labrouhe, D. (2003). *Del Mildiu Del Girasol: Manejo De Estrategias De Control Para Una Eficacia Durable, 2º Congreso Argentino De Girasol*. Available at: www.asagir.org.ar/acerca-de-2do-congreso-484 63–70 (section sanidad).

Tourvieille de Labrouhe, D., Serre, F., Walser, P., Roche, S., and Vear, F. (2008). Quantitative resistance to downy mildew (*Plasmopara halstedii*) in sunflower (*Helianthus annuus*). *Euphytica* 164, 433–444. doi: 10.1007/s10681-008-9698-1

Tourvieille de Labrouhe, D., Walser, P., Jolivot, D., Roche, S., Serre, F., Leguillon, M., et al. (2012). "Proposal for improvement of sunflower downy mildew race nomenclature," in *Proceedings of the 18th International Sunflower Conference*, Mar-del-Plata, 328–333.

Vear, F., Gentzbittel, L., Philippon, J., Mouzeyar, S., Mestries, E., Roeckel-Drevet, P., et al. (1997). The genetics of resistance to five races of downy mildew (*Plasmopara halstedii*) in sunflower (*Helianthus annuus* L.). *Theor. Appl. Genet.* 95, 584–589. doi: 10.1007/s001220050599
Vear, F., Serieys, H., Petit, A., Serre, F., Boudon, J.-P., Roche, S., et al. (2008a). "Origins of major genes for downy mildew resistance in sunflower," in *Proceedings of the 17th International Sunflower Conference*, Cordoba.

Vear, F., Serre, F., Jouan-Dufournel, I., Bert, P. F., Roche, S., Walser, P., et al. (2008b). Inheritance of quantitative resistance to downy mildew (*Plasmopara halstedii*) in sunflower (*Helianthus annuus* L.). *Euphytica* 164, 561–570. doi: 10.1007/s10681-008-9759-5

Vincourt, P., As-sadi, F., Bordat, A., Langlade, N. B., Gouzy, J., Pouilly, N., et al. (2012). Consensus mapping of major resistance genes and independent QTL for quantitative resistance to sunflower downy mildew. *Theor. Appl. Genet.* 125, 909–920. doi: 10.1007/s00122-012-1882-y

Viranyi, F., Gulya, T. J., and Tourvieille de Labrouhe, D. (2015). Recent changes in the pathogenic variability of *Plasmopara halstedii* (sunflower downy mildew) populations from different continents. *Helia* 38, 149–162. doi: 10.1515/helia-2015-0009

Vleeshouwers, V. G. A. A., Rietman, H., Krenek, P., Champouret, N., Young, C., Oh, S.-K., et al. (2008). Effector genomics accelerates discovery and functional profiling of potato disease resistance and *Phytophthora infestans* avirulence genes. *PLoS One* 3:e2875. doi: 10.1371/journal.pone.0002875

Voorrips, R. E. (2002). MapChart: software for the graphical presentation of linkage maps and QTLs. *J. Hered.* 93, 77–78. doi: 10.1093/jhered/93.1.77

Wieckhorst, S., Bachlava, E., Dussle, C. M., Tang, S., Gao, W., Saski, C., et al. (2010). Fine mapping of the sunflower resistance locus Pl (ARG) introduced from

the wild species *Helianthus argophyllus*. *Theor. Appl. Genet.* 121, 1633–1644. doi: 10.1007/s00122-010-1416-4

Zhang, M., Liu, Z., and Jan, C. C. (2016). Molecular mapping of a rust resistance gene R-14 in cultivated sunflower line PH 3. *Mol. Breed.* 36:32. doi: 10.1007/s11032-016-0456-0

Zhang, Z. W., Ma, G. J., Zhao, J., Markell, S. G., and Qi, L. L. (2017). Discovery and introgression of the wild sunflower-derived novel downy mildew resistance gene Pl (19) in confection sunflower (*Helianthus annuus* L.). *Theor. Appl. Genet.* 130, 29–39. doi: 10.1007/s00122-016-2786-z

Zimmer, D. E., and Kinman, M. L. (1972). Downy mildew resistance in cultivated sunflower and its inheritance. *Crop Sci.* 12, 749–751. doi: 10.2135/cropsci1972.0011183X001200060009x

**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.

Copyright © 2018 Pecrix, Penouilh-Suzette, Muños, Vear and Godiard. 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.