# Clubroot Symptoms and Resting Spore Production in a Doubled Haploid Population of Oilseed Rape (Brassica napus) Are Controlled by Four Main QTLs

Andrea Botero-Ramírez, Anne Laperche, Solenn Guichard, Mélanie Jubault, Antoine Gravot, Stephen Strelkov, Maria M. Manzanares-Dauleux

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

## To cite this version:

Andrea Botero-Ramírez, Anne Laperche, Solenn Guichard, Mélanie Jubault, Antoine Gravot, et al. Clubroot Symptoms and Resting Spore Production in a Doubled Haploid Population of Oilseed Rape (Brassica napus) Are Controlled by Four Main QTLs. Frontiers in Plant Science, 2020, 11, 10.3389/fpls.2020.604527. [hal-03137395](https://hal.inrae.fr/hal-03137395)

HAL Id: hal-03137395

https://hal.inrae.fr/hal-03137395

Submitted on 11 Feb 2021

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

# Clubroot Symptoms and Resting Spore Production in a Doubled Haploid Population of Oilseed Rape (*Brassica napus*) Are Controlled by Four Main QTLs

Andrea Botero-Ramírez<fnref n="1" />, Anne Laperche², Solenn Guichard², Mélanie Jubault², Antoine Gravot², Stephen E. Strelkov<fnref n="1" />* and Maria J. Manzanares-Dauleux²*

**OPEN ACCESS**

**Edited by:**

Harsh Raman,
New South Wales Department
of Primary Industries, Australia

**Reviewed by:**

Xiaoming Wu,
Chinese Academy of Agricultural
Sciences, China
Mulatu Geleta,
Swedish University of Agricultural
Sciences, Sweden

**\*Correspondence:**

Stephen E. Strelkov
strelkov@ualberta.ca
Maria J. Manzanares-Dauleux
maria.manzanares-dauleux@inrae.fr

**Specialty section:**

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

**Received:** 09 September 2020

**Accepted:** 20 November 2020

**Published:** 17 December 2020

**Citation:**

Botero-Ramírez A, Laperche A,
Guichard S, Jubault M, Gravot A,
Strelkov SE and
Manzanares-Dauleux MJ (2020)
Clubroot Symptoms and Resting
Spore Production in a Doubled
Haploid Population of Oilseed Rape
(Brassica napus) Are Controlled by
Four Main QTLs.
Front. Plant Sci. 11:604527.
doi: 10.3389/fpls.2020.604527

Clubroot, caused by *Plasmodiophora brassicae* Woronin, is one of the most important diseases of oilseed rape (*Brassica napus* L.). The rapid erosion of monogenic resistance in clubroot-resistant (CR) varieties underscores the need to diversify resistance sources controlling disease severity and traits related to pathogen fitness, such as resting spore production. The genetic control of disease index (DI) and resting spores per plant (RSP) was evaluated in a doubled haploid (DH) population consisting of 114 winter oilseed rape lines, obtained from the cross ‘Aviso’ × ‘Montego,’ inoculated with *P. brassicae* isolate “eH.” Linkage analysis allowed the identification of three quantitative trait loci (QTLs) controlling DI (PbBn_di_A02, PbBn_di_A04, and PbBn_di_C03). A significant decrease in DI was observed when combining effects of the three resistance alleles at these QTLs. Only one QTL, PbBn_rsp_C03, was found to control RSP, reducing resting spore production by 40%. PbBn_rsp_C03 partially overlapped with PbBn_di_C03 in a nucleotide-binding leucine-rich repeat (NLR) gene-containing region. Consideration of both DI and RSP in breeding for clubroot resistance is recommended for the long-term management of this disease.

**Keywords:** oilseed rape, *Plasmodiophora brassicae*, clubroot, linkage analysis, quantitative resistance, *Rhizaria*

## INTRODUCTION

Clubroot, caused by the obligate parasite *Plasmodiophora brassicae* Woronin, is one of the most important diseases of cruciferous crops worldwide, causing significant yield and quality losses in oilseed rape (canola; *Brassica napus* L.) and other Brassicas (Dixon, 2009). A soilborne disease, clubroot is associated with the formation of large galls on the roots of susceptible hosts, which impede water and nutrient uptake. The life cycle of *P. brassicae* consists of three main stages: (i) survival in the soil, (ii) root hair infection, and (iii) cortical infection. The pathogen survives in

<footnote n="1">¹ Department of Agricultural, Food and Nutritional Sciences, Faculty of Agricultural, Life and Environmental Sciences, University of Alberta, Edmonton, AB, Canada, ² Institut de Génétique, Environnement et Protection des Plantes, INRAE, Institut Agro, Université de Rennes 1, Le Rheu, France</footnote>

<footnote>**Abbreviations:** CR, clubroot-resistant; dai, days after inoculation; DH, doubled haploid; DI, disease index; RSP, resting spores per plant.</footnote>
the soil as long-lived resting spores (Kageyama and Asano, 2009). Under favorable conditions, the resting spores germinate to produce primary zoospores, which infect the host root hairs. The presence of plant root exudates can enhance resting spore germination (Macfarlane, 1970; Friberg et al., 2005; Rashid et al., 2013). Primary plasmodia develop within the infected root hairs, eventually giving rise to secondary zoospores. The secondary zoospores penetrate the cortical tissue and develop into intracellular secondary plasmodia, which cleave to produce a new generation of resting spores (Kageyama and Asano, 2009). As the galls decompose, these resting spores are released back into the soil, where they serve as inoculum for future infections. As many as $1 \times 10^{10}$ resting spores per g of gall tissue may be produced on a susceptible *B. napus* host (Hwang et al., 2013).

The persistence of *P. brassicae* resting spores in the soil (Wallenhammar, 1996) makes clubroot a particularly difficult disease to manage. Various strategies have been proposed, including crop rotation, biological control, liming of the soil, and manipulation of the seeding date (Hwang et al., 2014). The most effective strategy for clubroot management, however, is the deployment of CR host cultivars (Rahman et al., 2014). The resistance in most CR oilseed rape cultivars is monogenic and derived from the European winter oilseed rape 'Mendel' (Diederichsen et al., 2014; Fredua-Agyeman et al., 2018). 'Mendel' was one of the first CR winter oilseed rape cultivars with acceptable agronomic performance released in Europe (Diederichsen et al., 2003, 2009). The clubroot resistance in 'Mendel' is based on one dominant and race- or pathotype-specific gene (Diederichsen et al., 2006). While genetic resistance is the most effective and convenient method to manage clubroot, the selection pressure it imposes on *P. brassicae* populations can cause rapid shifts in virulence. These shifts may result in a loss of resistance, as has already been documented in oilseed rape/canola in Europe (Diederichsen et al., 2014; Orgeur et al., 2016) and Canada (Strelkov et al., 2016, 2018).

The diversification and introduction of novel resistance sources can reduce the risk of resistance loss, as can the development of cultivars with polygenic resistance and the implementation of strategies such as gene pyramiding and the planting of multilines or cultivar mixtures (Parlevliet and Zadoks, 1977; Pink and Puddephat, 1999). For example, a combination of major resistance genes and quantitative trait loci (QTLs) could result in CR cultivars with resistance that is more durable. The development of cultivars with more diverse resistance, however, requires a deeper understanding of the genetic control of clubroot resistance (Manzanares-Dauleux et al., 2000).

Most genetic studies of clubroot resistance have focused on understanding the resistance harbored in *Brassica oleracea* and *Brassica rapa*, the ancestral parents of *B. napus*. In *B. rapa*, multiple major genes controlling clubroot resistance have been identified, including *Crr1a*, *Crr1b* (*Crr1a* and *Crr1b* were initially identified as a single locus, *Crr1*), *Crr2*, *Crr3*, *Crr5A*, *CRA*, *CRb*, *CRd*, *CRs*, *Rcr1*, *Rcr2*, *Rcr3*, *Rcr5*, and *Rcr9*<sup>WA</sup> (Matsumoto et al., 1998; Suwabe et al., 2003; Hirai et al., 2004; Saito et al., 2006; Hatakeyama et al., 2013, 2017; Chu et al., 2014; Huang et al., 2017, 2019; Yu et al., 2017; Nguyen et al., 2018; Pang et al., 2018; Laila et al., 2019; Karim et al., 2020). In addition, two QTLs were

reported to control resistance to a Korean isolate of *P. brassicae* classified as per Williams (1966) pathotype 2 (Choi et al., 2020), and two QTLs were found to control resistance to a Chinese isolate of pathotype 7 (Zhu et al., 2019). Another QTL (*Rcr4*) was found to control resistance to Canadian isolates of the pathogen representing pathotypes 2, 3, 5, 6, and 8, classified as per Williams (1966), and two QTLs (*Rcr8* and *Rcr9*) controlled resistance to an isolate of pathotype 5X, as defined on the Canadian Clubroot Differential set (Yu et al., 2017; Strelkov et al., 2018). Three loci (*Crr4*, *CRc*, and *CRk*) controlled resistance to non-pathotyped isolates of *P. brassicae* (Suwabe et al., 2006; Sakamoto et al., 2008). In *B. oleracea*, clubroot resistance has been found to be quantitative and is controlled mainly by QTLs with both major and minor effects (Rocherieux et al., 2004; Neik et al., 2017).

Genetic analyses of the control of clubroot resistance in *B. napus* have identified the presence of a major resistance gene (*Cra*) (Zhang et al., 2016) and nearly 30 different QTLs in various plant populations harboring resistance to multiple *P. brassicae* strains (Manzanares-Dauleux et al., 2000; Werner et al., 2008; Li et al., 2016; Laperche et al., 2017; Aigu et al., 2018; Hejna et al., 2019; Wagner et al., 2019). A series of studies by our group using a DH population derived from the varieties 'Darmor-bzh' and 'Yudal' allowed the identification of a major QTL controlling resistance to *P. brassicae* isolate Pb137-522, while two QTLs were found to govern resistance to isolate K92-16 (Manzanares-Dauleux et al., 2000). In addition, when this DH population was tested against the pathogen isolate "eH," two major and one minor QTLs were found to control resting spore production (Aigu et al., 2018), one moderate and three minor QTLs controlled the pathogen-plant genomic DNA ratio (Wagner et al., 2019), and nine QTLs regulated clubroot severity, measured as a DI (Laperche et al., 2017; Aigu et al., 2018; Wagner et al., 2019). Moreover, other groups have identified different QTLs controlling clubroot resistance in one DH population challenged with seven different isolates of *P. brassicae* (Werner et al., 2008), in an associative transcriptome analysis of 245 accessions inoculated with the European Clubroot Differential (ECD) pathotype 17/31/31 (Hejna et al., 2019), and in a genome wide association study of 472 host accessions inoculated with Williams' pathotype 4 (Li et al., 2016).

Previous studies have included very diverse genetic materials (spring, old winter lines), whose use in breeding programs may be difficult given potential issues such as linkage drag (Zamir, 2001; Yousef and Juvik, 2002; Lecomte et al., 2004). The challenges associated with the introduction of polygenic partial resistance from non-elite materials into elite oilseed rape genetic backgrounds may be one reason for the limited introduction of QTLs into new CR cultivars. Therefore, the identification of loci conferring partial clubroot resistance in recent cultivars with good agronomic performance may be of interest, as breeders could use these more readily.

Sustainable clubroot management requires both short- and long-term approaches. In the short term, the focus should be on minimizing the impact of the pathogen on the host and therefore on crop production; in the longer term, the aim should be on reducing inoculum levels. Resistant cultivars are a convenient and highly effective disease management tool in the
short term, but may not prove durable over the long term. The loss of effective pathogen control may reflect the strong selection pressure imposed on the pathogen, particularly since plants with low disease severity do not necessarily produce low amounts of resting spores (Siemens et al., 2002). A recent evaluation of *P. brassicae* resting spore dynamics in response to the cropping of CR *B. napus* cultivars indicated increases in soil inoculum loads (Ernst et al., 2019).

Plant disease epidemics are highly influenced by pathogen virulence (measured as disease severity) and fitness (measured as the reproductive rate of the pathogen), both of which are generally assumed to be positively correlated (Sacristán and García-Arenal, 2008); clubroot is no exception. However, studies with respect to the relationship between pathogen virulence and reproductive rate have produced contradictory results in some pathosystems (Fox and Williams, 1984; Kover and Schaal, 2002; Montarry et al., 2006). Indeed, there are reports where pathogen inoculum production and virulence are not correlated, or even negatively correlated (Robert et al., 2002; Sacristán et al., 2005; Montarry et al., 2010; Aigu et al., 2018).

Given the assumption of a high correlation between disease severity and pathogen fitness, most genetic studies of clubroot resistance have focused only on disease severity (Manzanares-Dauleux et al., 2000; Werner et al., 2008; Laperche et al., 2017). However, while there is a relationship between *P. brassicae* resting spore production and root gall size (disease severity), these variables are not necessarily correlated, and resting spore production can be affected by host resistance and environmental factors (Murakami et al., 2004; Gravot et al., 2016; Aigu et al., 2018). Therefore, both traits should be considered for genetic analyses and in breeding programs. The new CR ideotypes should combine low resting spore production and low DI. Selection for both traits will facilitate improved disease management over the short and long term. In the short term, such an approach will minimize the direct impact of the pathogen on disease severity and hence on crop production, while in the long term, it will reduce the number of resting spores released into the soil, thereby limiting inoculum build-up and the potential for future epidemics.

This research had two objectives. First, it aimed to identify the QTLs involved in the control of resting spore production per plant (RSP) and clubroot symptoms (DI) in a segregating DH population from the cross of 'Aviso' × 'Montego,' two modern winter oilseed rape cultivars partially resistant to clubroot with good agronomic characteristics. Second, it aimed to identify some recombinant lines of potential interest for breeding efforts, i.e., carrying favorable alleles at multiple QTLs to decrease DI and limit resting spore production.

# MATERIALS AND METHODS

## Pathogen Material and Inoculum Preparation

The *P. brassicae* selection isolate “eH” (Fähling et al., 2003), classified as pathotype P₁ on the system of Somé et al. (1996), was used for all the experiments. To prepare inoculum, resting

spores were extracted from frozen root galls of the universally susceptible Chinese cabbage (*B. rapa* subs. *pekinensis*) cv. Granaat (ECD 05; Buczacki et al., 1975) inoculated with the isolate. Briefly, the galls were thawed at room temperature and then homogenized at maximum speed in a home blender. The resulting homogenate was filtered sequentially, first through cheesecloth and then through 56 and 100-µm diameter pore stainless steel sieves (Retsch, Haan, Germany). The resting spore concentration of the filtered suspension was estimated by counting in a Malassez cell and adjusted to a final concentration of 1 × 10⁷ resting spores ml⁻¹ with sterile distilled water.

## Greenhouse Experiments and Disease Assessment

A 2-year experiment was conducted under greenhouse conditions in 2015 and 2016. The experimental design consisted of completely randomized blocks nested within the 2 years; in total, 114 genotypes (treatments) with four replicates (blocks) were established, with the experimental unit comprising six plants per genotype. Seeds of each genotype were sown in 4-cm-diameter pots (one seed per pot) filled with “Falienor 922016F3” potting mix (Falienor, Vivy, France), which consists of 65% Irish peat, 20% black peat, 15% perlite, and 2% clay (pH = 6.2). The greenhouse was maintained at temperatures between 19 and 22°C under a 16/8 day/night cycle. Plants were fertilized with “Liquoplant FD 134 hiver” nutrient solution (Plantin, Courthéson, France) once or twice a week by sub-irrigation. Inoculations were conducted 7 days after sowing, by applying 1 ml of the *P. brassicae* resting spore suspension to the base of each seedling.

Disease assessment was conducted 54 dai on a 0-3 scale following Manzanares-Dauleux et al. (2000), where 0 = no visible galling, 1 = very light galling usually confined to lateral roots, 2 = moderate galling on lateral roots and the taproot, 2 + = severe galling on all roots but some roots remain healthy, and 3 = one large gall with no remaining healthy roots. The individual severity ratings were then used to calculate a DI using Eq. 1:

$$DI = \frac{[(0 \times n_0) + (25 \times n_1) + (50 \times n_2) + (75 \times n_{2+}) + (100 \times n_3)] \times 100}{N}$$

Where $n_0$, $n_1$, $n_2$, $n_{2+}$, and $n_3$ represent the number of plants in each severity class and N is the total number of plants evaluated. To confirm the pathotype designation of the isolate “eH” as P₁, the isolate was inoculated on the differential hosts of Somé et al. (1996) [*B. napus* cv. Nevin (ECD 06), *B. napus* cv. Wilhelmsburger (ECD 10), and *B. napus* cv. Brutor]. The Chinese cabbage ECD 05 was also included as a susceptible check in all of the experiments.

After disease assessment, all roots in the experimental unit were pooled and stored at -20°C until processing. Resting spores in the pooled root samples were quantified by flow cytometry following Aigu et al. (2018). Briefly, the roots were thawed at room temperature and weighted. Each sample was homogenized in 100 ml of distilled water in a home blender, and the resulting suspension was filtered as described for the inoculum preparation. The resting spore suspensions were diluted in a
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Distribution of the clubroot disease index (DI) for the doubled haploid progeny from a cross of the oilseed rape cultivars 'Aviso' × 'Montego' following inoculation with *Plasmodiophora brassicae* isolate "eH." The parents 'Aviso' and 'Montego' are highlighted.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Distribution of number of *Plasmodiophora brassicae* resting spores produced per plant for the doubled haploid progeny from a cross of the oilseed rape cultivars 'Aviso' × 'Montego' following inoculation with the pathogen isolate "eH." The parents 'Aviso' and 'Montego' are highlighted.</figcaption></figure>

1:20 ratio with distilled water. Resting spores were quantified in a CyFlow flow cytometer (Sysmex partec, Görlitz, Germany) equipped with a 20-mW blue laser (488 nm) and a Forward-Scatter (FCS) detector to determine particle size (Aigu et al., 2018).

A standard curve was generated from a suspension of resting spores, produced as described above for inoculum preparation,

with the spore concentration estimated by counting in a Malassez cell. This solution was then diluted as needed to generate nine spore suspensions with concentrations ranging from 1 × 10⁴ to 6 × 10⁵ resting spores·ml⁻¹. Particle counting was completed for each of the dilution points with the flow cytometer and an injection speed of 20. The measured resting spore concentrations were adjusted with a regression Eq. (2) obtained from the
<figure><img src="image_6.png" /><figcaption>FIGURE 3 | Relationship between disease index (DI) and number of resting spores per plant (RSP) for the doubled-haploid progeny from a cross of the oilseed rape cultivars ‘Aviso’ × ‘Montego.’ The Spearman’s rank correlation coefficient was 0.65 (p-value < 0.001). Numbers over the highlighted points indicate the code of each line with atypical behavior. Green points represent lines with a DI < 20 and an RSP < 5.9 × 10⁷ resting spores plant⁻¹; blue points represent lines with a DI between 20 and 40 and an RSP < 5.9 × 10⁷ resting spores plant⁻¹; orange points represent lines with DI < 20 and an RSP > 8 × 10⁷ resting spores plant⁻¹; yellow points represent lines with DI between 30 and 50 and an RSP > 9 × 10⁷ resting spores plant⁻¹.</figcaption></figure>

standard curve. The obtained regression curve had an R² of 95%. Spore levels were expressed as RSP by dividing the total number of resting spores in the suspension by the number of plants in the pooled sample.

$$RS = 2 \times 10^{-6}x^2 + 2.2283x + 41716 \quad (1)$$

Where RS is the adjusted resting spore count in the suspension and x is the particle count obtained with the flow cytometer.

## Plant Material and Genetic Map

A population of 114 DH lines obtained from the cross of the winter oilseed rape cultivars ‘Aviso’ × ‘Montego’ was chosen for linkage analysis because previous experiments indicated that both parents are partially resistant to clubroot, and therefore the presence of new QTLs controlling the disease was suspected. The parental lines were used as controls in all experiments. The genetic map used for the analysis was described by Delourme et al. (2013). Additional genotyping was conducted with the 60K infinitum array (Clarke et al., 2016), thus leading to an updated map. Given the population size, recombination was not always possible, resulting in a high number of markers at the same genetic position. Since this is not desirable for QTL mapping, we only kept a single marker or unique loci to represent each cluster. The genetic map covers 1947 cM (892 cM for the A genome and 1055 cM for the C genome) at a density of 1.18 markers each cM; it comprises 2301 SNPs representing 831 unique loci. The linkage groups (LGs) with the highest proportion of loci distortion were A02, A03, A09, C01, and C09.

## Statistical Analyses

Statistical analyses were performed with R (R Core Team, 2019). Spearman’s rank correlation coefficient was estimated to evaluate

the correlation between RSP and DI. The mixed linear model presented in Eq. 3 was estimated to analyze RSP and DI. The nlme (Pinheiro et al., 2019) and lsmeans (Lenth, 2016) packages were used to test both random and fixed effects:

$$Y_{ijk} = \mu + G_i + Y_j + B_{k(j)} + GY_{ij} + \varepsilon_{ijk} \quad (2)$$

Where $Y_{ijk}$ is either RSP or DI in each genotype ($i^{th}$) observed in the $j^{th}$ year and in the $k^{th}$ block; $\mu$ is the population mean; $G_i$ is the effect of the $i^{th}$ genotype; $Y_j$ is the effect of the $j^{th}$ year; $B_{k(j)}$ is the effect of the $k^{th}$ block nested with the $j^{th}$ year; $GY_{ik}$ is the effect of the interaction between the $i^{th}$ genotype and the $j^{th}$ year; and $\varepsilon_{ijk}$ are the residuals. $B_{k(j)}$ was considered as random.

Broad sense heritability as defined by Holland et al. (2003) (Eq. 4) was estimated using the model in Eq. 3.

$$H^2 = \frac{\sigma^2 G}{\sigma^2 G + \frac{\sigma^2 Gy}{y} + \frac{\sigma^2 e}{y * b}} \quad (3)$$

Where $H^2$ is the broad sense heritability, $\sigma^2 G$ is the variance of the genotype effect, $\sigma^2 Gy$ is the variance of the genotype × year effect, y is the number of years, b is the number of blocks, and $\sigma^2 e$ is the residual variance.

Adjusted means were calculated using the lsmeans package (Lenth, 2016) for each genotype across years and replications and used as phenotypic values for QTL analyses.

The reactions of the progeny were grouped based on genotypes which possessed the identified QTLs for DI and RSP, and the phenotypic responses were compared using Duncan’s new multiple range test at $P \le 0.05$.
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Genome scan of LOD scores for disease index (DI). (A) The LOD scores estimated by simple interval mapping (SIM). (B) The LOD scores for the chromosomes where QTLs were identified (A02 = chromosome 2, A04 = chromosome 4, C03 = chromosome 13); the LOD scores estimated by multiple QTL mapping are indicated in black, while the LOD scores estimated by SIM are indicated in blue. The red line in both graphs represents the LOD threshold determined by 1000 permutations.</figcaption></figure>

## Linkage Analysis

Linkage analysis was conducted with the R/qlt package (Broman et al., 2003). At first, a simple interval mapping (SIM) was performed to get a rough estimate of the QTLs controlling each trait (RSP and DI). The LOD threshold for these analyses was 3.17, estimated by a 1000 permutation test ($\alpha$ = 0.05). Afterward, multiple QTL mapping was conducted. Manual selection of the QTL model was completed with the functions *addqtl*, *addintqtl*, and *fitqtl*, allowing for QTL-pairwise interaction using the multiple imputation regression method. The QTLs were added one by one and at each step, and two models were compared: one with the already validated QTLs, and the same model including the tested QTL and the corresponding QTL by QTL interactions. Only QTLs whose effect was significant ($\alpha$ = 0.05) according to the ANOVA table were retained in the model. LOD values and $R^2$ values for each QTL were also obtained by the *fitqtl* function. The confidence intervals of the QTLs were estimated with a LOD drop of one unit.

For each QTL, the genes present in the confidence interval were gathered using the reference genome v4 of 'Darmor-bzh' (Chalhoub et al., 2014) by using the physical anchorage of the confidence interval flanking SNP markers.

## RESULTS

### Characterization of the Phenotypic Response

The frequency distribution for DI (Figure 1) and RSP (Figure 2) indicated a continuous distribution, suggesting polygenic control of both variables. In both cases, the parents 'Aviso' and 'Montego' showed intermediate clubroot responses. Mean DI was 36.8 for 'Aviso' and 46.5 for 'Montego'; the RSP was $7.2 \times 10^7$ resting spores plant$^{-1}$ for 'Aviso' and $7.6 \times 10^7$ resting spores plant$^{-1}$ for 'Montego.' RSP was significantly correlated with DI ($p$ value < 0.001) (Figure 3), resulting in a Spearman's coefficient
<table><caption>**TABLE 1** | QTLs controlling clubroot disease index (DI) or number of *Plasmodiophora brassicaceae* resting spores per plant (RSP), and their position in the physical map, identified by multiple QTL mapping in a doubled haploid population obtained from a cross of the oilseed rape cultivars 'Aviso' × 'Montego'.</caption><thead><tr><th>QTL name</th><th>Trait</th><th>Chromosome</th><th>Position (cM)</th><th>CI (cM)</th><th>Markers CI</th><th>LOD</th><th>CI (cM)</th><th>R<sup>2</sup> (%)</th><th>Favorable allele source</th><th>Additive effect</th><th>Position in physical map (Mb)</th><th>Number of genes</th><th>Unannotated or proteins with unknown function</th><th>Genes related with plant resistance</th></tr></thead><tbody><tr><td>PbBn_di_A02</td><td>DI</td><td>A02</td><td>28.9</td><td>28.0-34.0</td><td>BS008863-BS009106</td><td>4.7</td><td>28.0-34.0</td><td>4.5</td><td>Montego</td><td>-6.93</td><td>4.84-5.50</td><td>101</td><td>23</td><td>5</td></tr><tr><td>PbBn_di_A04</td><td></td><td>A04</td><td>57.6</td><td>57.6-59.4</td><td>BS006202-BS006447</td><td>15.2</td><td>57.6-59.4</td><td>18.1</td><td>Montego</td><td>-7.737</td><td>1.34</td><td></td><td></td><td></td></tr><tr><td>PbBn_di_C03</td><td></td><td>C03</td><td>30.8</td><td>28.1-32.6</td><td>BS007532-Bn-C3-p5080881</td><td>30.3</td><td>28.1-32.6</td><td>51</td><td>Aviso</td><td>13.193</td><td>4.09-4.88</td><td>147</td><td>29</td><td>6</td></tr><tr><td>PbBn_rsp_C03</td><td>RSP</td><td>C03</td><td>24.5</td><td>19.1-29.0</td><td>BS012716-Bn-C3-p4469843</td><td>5.1</td><td>19.1-29.0</td><td>18.3</td><td>Aviso</td><td>8.5 × 10<sup>5</sup></td><td>2.89-4.50</td><td>346</td><td>72</td><td>10</td></tr></tbody></table>

CI, confidence interval.

of 0.65. The heritability for both variables was high (85.2% for DI and 84.4% for RSP).

Among the 114 recombinant progeny lines tested, 18% had a DI < 20; 53% had intermediate levels of disease (DI between 20 and 50), while the remaining 29% had a DI between 50.7 and 98. RSP ranged from 5.6 × 10⁷ to 1.4 × 10⁸, with about 54.2% of the lines exhibiting a higher RSP than both parents. Only 30% of the lines produced fewer resting spores than the parental lines (Figure 2).

Lines 88 and 189 developed mild symptoms of clubroot (DI of 7.3 and 10.5, respectively) but had significant resting spore production (RSP of 5.6 × 10⁷ and 5.9 × 10⁷ resting spores plant⁻¹, respectively) (Figure 3). Similarly, although lines 27, 64, 52, and 127 had a DI < 20, the RSP was fairly high, ranging from 8 × 10⁷ in lines 52 and 127 to 9 × 10⁷ resting spores plant⁻¹ in lines 27 and 64. Some lines with intermediate DI also showed high resting spore production, including line 83 (DI 30.5 and 1.2 × 10⁹ resting spores plant⁻¹), line 148 (DI 44.8 and 1.2 × 10⁹ resting spores plant⁻¹), line 120 (DI 36.3 and 9.2 × 10⁸ resting spores plant⁻¹), line 74 (DI 37.3 and 9.7 × 10⁸ resting spores plant⁻¹), line 168 (DI 44.8 and 9.3 × 10⁸ resting spores plant⁻¹), line 48 (DI 46.8 and 9.7 × 10⁸ resting spores plant⁻¹), and lines 109 and 39 (DI 48.9 and 9.2 × 10⁸ resting spores plant⁻¹). In contrast, some lines with an intermediate DI had a lower RSP, including lines 76 (DI 20.8 and 5.8 × 10⁷ resting spores plant⁻¹) and 164 (DI 32.8 and 5.7 × 10⁷ resting spores plant⁻¹) (Figure 3).

## QTLs Controlling DI and Number of Resting Spores per Plant

Quantitative trait loci were assigned names consisting of three parts separated by underscores. The first part of each name includes the initials Pb and Bn to indicate *P. brassicaceae* and *B. napus*, respectively. The second part indicates the trait controlled by the QTL, DI, or RSP, in lowercase letters. The third part indicates the chromosome on which the QTL is located. SIM indicated the presence of two main QTLs controlling DI on chromosomes A04 and C03 (Figure 4), PbBn_di_A04, and PbBn_di_C03. Further analysis allowed the identification of an additional QTL on chromosome A02, PbBn_di_A02, and the final fitted model accounted for 78.4% of the total variation. Most of the variance in DI was controlled by PbBn_di_C03 (51.0%) followed by PbBn_di_A04 (18.1%). A minor effect was detected for PbBn_di_A02, which explained 4.5% of the total variation. Since the population under investigation was a DH, only homozygous lines were tested, and thus only additive effects were detectable, leading to a strong relationship between the percentage of the variance accounted for by the QTLs and their additive effect. It was observed that the QTL with the strongest effect (PbBn_di_C03) also had the highest additive effect (Table 1). No epistasis was found among the QTLs.

Only one QTL controlling RSP, PbBn_rsp_C03 on the chromosome C03, was detected either by SIM or by fitting a multiple QTL model (Figure 5). The fitted model accounted for 18.3% of the total variation for that variable. PbBn_di_C03 and PbBn_rsp_C03 overlapped. On PbBn_di_C03 and PbBn_rsp_C03, the 'Aviso' allele contributed to phenotypes
with lower DI and RSP, respectively, while on PbBn_di_A04 and PbBn_di_A02, the 'Montego' alleles were more favorable for reducing DI (Table 1). The lines that did not have any favorable allele for disease reduction presented the highest DI; an intermediate DI was observed only when the favorable allele of PbBn_di_C03 was present, and the highest reductions in DI occurred whenever the PbBn_di_C03, PbBn_di_A04, and PbBn_di_A02 favorable alleles were present together (Figure 6).

It was confirmed that RSP was not controlled by the QTLs PbBn_di_A02 or PbBn_di_A04, since the presence of their favorable alleles did not cause a difference in this value (Figure 7).

## Physical Anchorage of the Confidence Intervals on the *B. napus* Genome

The markers names, positions, and start and end-points for all QTLs are indicated in Table 1.

Only three of the four identified QTLs could be mapped on the *B. napus* genome: PbBn_di_A02, PbBn_di_C03, and PbBn_rsp_C03. The QTL PbBn_di_A04 could not be mapped because the region where it was positioned corresponded to a chimeric region on the reference genome (v4). PbBn_di_A02 covered 661.48 kb in a region encoding 101 genes, while PbBn_di_C03 covered 788.420 kb in total, encoding 147 genes. Four genes related to plant resistance and defense were found in PbBn_di_A02 QTL. Two were annotated as hypersensitive-induced response protein (HIR), two as a cyclic nucleotide gated channel, and the fourth as a Mlo-like protein 6. In PbBn_di_C03, six genes related to disease resistance or defense were identified, including five TIR-NBS-LRR class disease resistance proteins and one LRR protein kinase-like protein.

The only QTL related to RSP, PbBn_rsp_C03, covered 1614.43 kb. In this QTL, 346 genes were encountered, of which 10 were related to disease resistance or defense, including five TIR-NBS-LRR class disease resistance proteins, three leucine-rich repeat (LRR) protein kinase-like proteins, one defensin-like protein 203, and one WRKY transcription factor 18. PbBn_rsp_C03 QTL overlapped with PbBn_di_C03 in the region from 28.10 to 29.01 cM (corresponding to a 0.39 Mb physical region). In the common region between PbBn_di_C03 and PbBn_rsp_C03 (4.09–4.50 Mb), three TIR-NBS-LRR class disease resistance proteins were found (Table 2).

## DISCUSSION

Research on clubroot resistance has focused mainly on disease severity expressed as a DI. Resting spore production in host genotypes has not been examined to the same extent, although this is an important measurement of pathogen fitness, affecting inoculum build up in infested fields. In the 'Aviso' × 'Montego' population, the highest observed RSP ($1.4 \times 10^8$ resting spores plant$^{-1}$) was 2.5 times greater than the lowest ($5.6 \times 10^7$ resting spores plant$^{-1}$). Such variation may be important from an epidemiological perspective, since doubling the number of resting spores released into the soil could result in faster and more significant inoculum increases for future crops.

Our results indicate that RSP is not as strongly correlated with DI (Spearman's coefficient of 0.65) as might be expected based on the trade-off hypothesis, where within-host multiplication, within-host transmission, and virulence of pathogens are positively correlated traits (Frank, 1996); and thus, higher disease levels would result in higher within host reproduction rates (Sacristán and García-Arenal, 2008). Nonetheless, these results are consistent with previous reports on the clubroot pathosystem. In a study with *Arabidopsis thaliana*, Siemens et al. (2002) observed that correlation between DI and spore number per root weight was between 0.7 and 0.9, and that one of the main factors determining those correlation values was host resistance. Similarly, Murakami et al. (2004) concluded that resting spore production is host-specific and, therefore, clubroot severity cannot account for resting production on its own. This was especially true for intermediate disease severities; for example, Chinese cabbage and broccoli plants with intermediate levels of clubroot produced approximately $1 \times 10^9$ resting spores plant$^{-1}$, while cabbage plants with the same amount of disease produced only $1 \times 10^8$ resting spores plant$^{-1}$ (Murakami et al., 2004). More recently, Aigu et al. (2018) reported that *P. brassicae* was able to produce high numbers of resting spores in some *B. napus* genotypes with mild symptoms of clubroot, and highlighted the partially resistant genotype 'Darmor-bzh' ($2.7 \times 10^8$ resting spores plant$^{-1}$ and DI = 30).

Linkage analysis showed that the genetic control of RSP and DI are related, since the QTLs PbBn_di_C03 and PbBn_rsp_C03 co-localized, suggesting that the C03 locus controls both traits. The importance of this genomic region is highlighted by the fact that PbBn_di_C03 explains the highest proportion of the variation for DI. The other loci, PbBn_di_A04 and PbBn_di_A02, only control DI and not RSP. We have reported this type of genetic architecture previously, where all QTLs controlling RSP also control DI, but not all DI QTLs control RSP (Laperche et al., 2017; Aigu et al., 2018). Co-localization of two QTLs with an intermediate effect on DI (QTL controlling < 20% of the variation) and a high effect on RSP (QTL controlling > 50% of the variation) has also been reported (Laperche et al., 2017; Aigu et al., 2018). Collectively, these studies suggest that co-localization of DI and RSP QTLs does not depend on the QTL effect (percentage of the variance of the trait that the QTL explains).

In a recent study of the genomic regions controlling DI and RSP in response to *P. brassicae* isolate "eH" in another *B. napus* DH population ('Darmor-bzh' × 'Yudal'), five QTLs controlling DI were located on chromosomes A05, A07, C02, C03, and C09 (Laperche et al., 2017). Two of these QTLs, on chromosomes C02 and C09, were also found to control RSP (Aigu et al., 2018). In the current study, we also identified QTLs on chromosome C03, but ours were in the regions 4.09–4.88 Mb (PbBn_di_C03) and 2.89–4.5 Mb (PbBn_rsp_C03), while the one found by Laperche et al. (2017) occurred in the region 4.6–5.0 Mb. This suggests that the same genomic region may be controlling clubroot resistance in both populations; however, further analysis is required for confirmation. The percentage of the variation explained by PbBn_di_C03 in our study and the QTL on chromosome C03 reported by Laperche et al. (2017)
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Genome scan of LOD scores for the number of resting spores per plant. (A) The LOD scores estimated by simple interval mapping (SIM). (B) The LOD scores for the chromosome C03 where one QTL was identified (C03 = chromosome 13); the LOD scores estimated by multiple QTL mapping are indicated in black, while the LOD scores estimated by SIM are indicated in blue. The red line represents the LOD threshold determined by 1000 permutations.</figcaption></figure>

in the 'Darmor-bzh' × 'Yudal' population was very different. The QTLs we detected on C03 controlled 51.4% of DI and 18.3% of the variation in RSP, while Laperche et al. (2017) found that the QTL on chromosome C03 controlled 7.75% of DI. No QTL controlling RSP was found on chromosome C03 in the 'Darmor-bzh' × 'Yudal' population (Aigu et al., 2018). Similarly, Werner et al. (2008) found QTLs in the LGs N13 and N02 (chromosomes C03, and A02, respectively) in a DH population of *B. napus* obtained from a cross of the DH line 263/11 and the oilseed rape cultivar 'Express' challenged with the *P. brassicae* isolates "01:60," "01:07," and "k." QTLs on chromosome C03 were identified when plants were inoculated with isolates "01:60" and "01:07" and explained 28.6 and 11.7% of the variance in DI, respectively. A QTL on chromosome A02 was also found after inoculation with the isolate "k," explaining 17.6% of the variance in DI. In addition, QTLs at syntenic positions have been identified in *B. oleracea* on chromosome C03 (Nagaoka et al., 2010; Lee et al., 2016; Li et al., 2016) and in *B. rapa* on chromosome A02 (Yu et al., 2017), suggesting that

the C03 genomic region merits further investigation. Indeed, the QTLs on chromosome C03 seem to be involved in the genetic control of response to different isolates, as has been observed in multiple genetic backgrounds harboring different effects (minor to major) depending on both the plant genotype and the *P. brassicae* isolate.

Most genetic analyses of the resistance harbored by oilseed rape against clubroot disease have been performed with very diverse genetic material (spring, old winter lines), and few resistance sources appear to be lines with good agronomic value. In the current study, the parents of the DH population were the cultivars 'Aviso' and 'Montego,' with the former released by Danisco seeds (Holeby, Denmark) in 2000 and the latter released by Limagrain (Saint-Beauzire, France) in 2002. 'Aviso' is a variety with good agronomic performance that shows resistance to other diseases such as blackleg (*Leptosphaeria maculans*) (Stonard et al., 2007; Jestin et al., 2015). These cultivars may represent suitable donors of clubroot resistance in oilseed rape breeding programs, considering that with their good agronomic value, less intensive
<figure><img src="image_9.png" /></figure>

**FIGURE 6** | Boxplot showing clubroot disease index (DI) in the doubled-haploid progeny from a cross of the oilseed rape cultivars 'Aviso' × 'Montego.' The reactions of the progeny are grouped based on genotypes which possess none of the favorable QTL alleles ("None"); the favorable allele of the QTL PbBn_di_C03 ("PbBn_di_C03"); and the favorable alleles of the QTLs PbBn_di_C03, PbBn_di_A04, and PbBn_di_A02 ("PbBn_di_C03 + PbBn_di_A04 + PbBn_di_A02"). Boxes with the same letter within the graph do not differ according to Duncan's new multiple range test at *p* > 0.05.

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

**FIGURE 7** | Boxplot showing *Plasmodiophora brassicae* resting spores per plant in the doubled-haploid progeny from a cross of the oilseed rape cultivars 'Aviso' × 'Montego.' The progeny are grouped based on genotypes which possess the favorable alleles of the QTLs PbBn_rsp_C03, PbBn_di_A04, and PbBn_di_A02 ("PbBn_rsp_C03 + PbBn_di_A04 + PbBn_di_A02"); the favorable alleles of the QTLs PbBn_rsp_C03 and PbBn_di_C03 ("PbBn_rsp_C03 + PbBn_di_C03"); the favorable alleles of the QTLs PbBn_rsp_C03 + PbBn_di_C03 + PbBn_di_A04 + PbBn_di_A02 ("PbBn_rsp_C03 + PbBn_di_C03 + PbBn_di_A04 + PbBn_di_A02"); or none of the favorable alleles ("None"). Boxes with the same letter within the graph do not differ according to Duncan's new multiple range test at *p* > 0.05.
<table><caption>**TABLE 2** | Annotation of the genes related to resistance and defense for the QTLs controlling clubroot disease index (DI) or number of *Plasmodiophora brassicae* resting spores per plant (RSP), in a doubled haploid population obtained from a cross of the oilseed rape cultivars 'Aviso' × 'Montego.'</caption><thead><tr><th>QTL name</th><th>Trait</th><th>Chromosome</th><th>Position in physical map (Mb)</th><th>Gene name</th><th>Position (bp)</th><th>Gene annotation</th><th>Homoeologous gene name</th><th>Arabidopsis thaliana orthologs</th><th>Brassica rapa ortologs</th><th>Brassica oleracea ortologs</th></tr></thead><tbody><tr><td rowspan="5">PbBn_di_A02</td><td rowspan="5">DI</td><td rowspan="5">A02</td><td rowspan="5">4.84-5.50</td><td>BnaA02g09860D</td><td>4945996-4948245</td><td>Hypersensitive-induced response protein (HIR) 2</td><td>BnaC02g13810D</td><td>AT5G54100</td><td>Bra022696</td><td>Bo2g041840</td></tr><tr><td>BnaA02g09870D</td><td>4950703-4954201</td><td>Hypersensitive-induced response protein (HIR) 2</td><td>BnaC02g13820D</td><td>AT5G54095</td><td>Bra022695</td><td>Bo2g042150</td></tr><tr><td>BnaA02g09790D</td><td>4902667-4907022</td><td>Cyclic nucleotide gated channel</td><td>BnaC02g44680D</td><td>AT5G54250</td><td>Bra022702</td><td>Bo2g042020</td></tr><tr><td>BnaA02g10440D</td><td>5360323-5362868</td><td>Cyclic nucleotide gated channel</td><td>BnaC02g14560D</td><td>AT5G53130</td><td>Bra022632</td><td>.</td></tr><tr><td>BnaA02g10080D</td><td>5102140-5105226</td><td>Mlo-like protein 6</td><td>BnaC02g14090D</td><td>AT5G53760</td><td>Bra022673</td><td>Bo2g044430</td></tr><tr><td rowspan="3">PbBn_di_C03</td><td rowspan="3"></td><td rowspan="3">C03</td><td rowspan="3">4.50-4.88</td><td>BnaC03g09410D</td><td>4510238-4510891</td><td>TIR-NBS-LRR class disease resistance proteins</td><td>.</td><td>.</td><td>Bra006487</td><td>Bo3g013450</td></tr><tr><td>BnaC03g09420D</td><td>4511004-4512486</td><td>TIR-NBS-LRR class disease resistance proteins</td><td>.</td><td>.</td><td>.</td><td>.</td></tr><tr><td>BnaC03g10100D</td><td>4867978-4869203</td><td>Leucine-rich repeat protein kinase-like protein</td><td>.</td><td>AT5G20480</td><td>Bra006560</td><td>Bo3g015190</td></tr><tr><td rowspan="8">PbBn_rsp_C03</td><td rowspan="8">RSP</td><td rowspan="8">C03</td><td rowspan="8">2.89-4.09</td><td>BnaC03g08650D</td><td>4082610-4087660</td><td>TIR-NBS-LRR class disease resistance proteins</td><td>.</td><td>AT5G17880</td><td>.</td><td>Bo3g012820</td></tr><tr><td>BnaC03g08660D</td><td>4089063-4094327</td><td>TIR-NBS-LRR class disease resistance proteins</td><td>.</td><td>AT5G17890</td><td>.</td><td>Bo3g012830</td></tr><tr><td>BnaC03g06310D</td><td>3057621-3060722</td><td>Leucine-rich repeat protein kinase-like proteins</td><td>BnaA03g04770D</td><td>AT5G14210</td><td>Bra006243</td><td>Bo3g009670</td></tr><tr><td>BnaC03g07510D</td><td>3539026-3542690</td><td>Leucine-rich repeat protein kinase-like proteins</td><td>BnaA03g05790D</td><td>AT5G16000</td><td>Bra006335</td><td>Bo3g010820</td></tr><tr><td>BnaC03g07920D</td><td>3723953-3729620</td><td>Leucine-rich repeat protein kinase-like proteins</td><td>.</td><td>.</td><td>Bra006366</td><td>Bo3g012140</td></tr><tr><td>BnaC03g08510D</td><td>4041438-4041833</td><td>Defensin-like protein 203</td><td>.</td><td>.</td><td>.</td><td>Bo3g012680</td></tr><tr><td>BnaC03g06770D</td><td>3236365-3238839</td><td>WRKY transcription factor 18</td><td>BnaA03g05230D</td><td>AT5G15130</td><td>Bra006283</td><td>Bo3g010100</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th>QTL name</th><th>Trait</th><th>Chromosome</th><th>Position in physical map (Mb)</th><th>Gene name</th><th>Position (bp)</th><th>Gene annotation</th><th>Homoeologous gene name</th><th>Arabidopsis thaliana orthologs</th><th>Brassica rapa ortologs</th><th>Brassica oleracea ortologs</th></tr></thead><tbody><tr><td>Common region between PbBn_di_C03 and PbBn_rsp_C03</td><td></td><td>C03</td><td>4.09-4.50</td><td>BnaC03g08900D</td><td>4222525-4224460</td><td>TIR-NBS-LRR class disease resistance proteins</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td>BnaC03g08920D</td><td>4227171-4232582</td><td>TIR-NBS-LRR class disease resistance proteins</td><td>BnaAnng27760D</td><td>AT5G18370</td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td>BnaC03g09010D</td><td>4264274-4265687</td><td>TIR-NBS-LRR class disease resistance proteins</td><td>BnaA03g07110D</td><td></td><td>Bra006458</td><td>B03g013130</td></tr></tbody></table>

backcrossing with elite lines would be required, facilitating the transfer of polygenic traits into new cultivars.

The potential utility of the 'Aviso' × 'Montego' population in breeding programs is underscored by the fact that the resting spore production was lower relative to values reported in previous studies. The 'Aviso' × 'Montego' DH population produced between 5.6 × 10⁷ and 1.45 × 10⁸ resting spores plant⁻¹. In contrast, in the 'Darmor-bzh' × 'Yudal' DH population, the RSP values ranged from 1 × 10⁸ and 5 × 10⁹ (Aigu et al., 2018). Hence, the 'Aviso' × 'Montego' DH population used in this research appears promising for the selection of parental lines with low DI and RSP (lines 189 and 88) for the transfer of the identified QTLs into new cultivars. While it has been reported that the simultaneous transfer of even four unlinked QTLs is possible (Hospital and Charcosset, 1997; Lecomte et al., 2004; Steele et al., 2006), the use of at least three markers for each QTL is recommended to make them useful in marker-assisted selection (Hospital and Charcosset, 1997). The transfer of the four QTLs identified in the current study is therefore possible, although additional refinement and validation of the markers is required for marker-assisted selection.

The search we conducted for the genes underlying the QTLs in the 'Darmor-bzh' reference genome allowed the identification of a set of genes, which might be involved in the response of the 'Aviso' × 'Montego' population to *P. brassicae*, on chromosomes A02 and C03. Several of the putative gene products belonged to protein families involved in plant defense and resistance, including HIRs, cyclic nucleotide gated channel, Mlo-like protein 6, TIR-NBS-LRR (TNL) class disease resistance proteins, LRR protein kinase-like, defensin-like protein 203, and WRKY transcription factor 18 (Rushton et al., 1995; Song et al., 1995; Buchanan and Gay, 1996; García-Olmedo et al., 1998; Clough et al., 2000; Eulgem et al., 2000; Dangl and Jones, 2001; Kanzaki et al., 2002; Choi et al., 2010; Zhou et al., 2010; Duan et al., 2013). The presence of a cluster of seven TNL-encoding genes is of special interest. The two clubroot-resistance loci (*Cra* and *Crr1*) that have been molecularly resolved so far (Ueno et al., 2012; Hatakeyama et al., 2013) both encode for TNLs. In addition, the fine mapping of the *Rcr1* locus (Yu et al., 2016) and GWAS have highlighted the recurrent presence of nucleotide-binding LRR (NLR)-gene clusters in clubroot-resistance loci, supporting the importance of this protein family in driving clubroot resistance (Stotz et al., 2018).

The 'Aviso' × 'Montego' population characterized in this study holds promise for the development of new CR oilseed rape cultivars, since the parental lines were recent varieties with good agronomic traits, and some of the tested lines presented very low DI accompanied by low RSP. While much of the emphasis in clubroot resistance breeding traditionally has been on selection based on low disease severity, resting spore production in plant material should also be considered. While breeding based on DI is critical for producing cultivars that perform well in *P. brassicae*-infested fields, the incorporation of lower RSP as a selection criterion will enable more sustainable clubroot management, by selecting plant genotypes where pathogen multiplication (fitness) is reduced.
DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this article will be
made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS

AB-R conducted some of the laboratory analyses, performed the
statistical and linkage analyses, and wrote the manuscript. AL,
AG, and MM-D contributed to development of the research
concept and design of the study and reviewed multiple versions of
the manuscript. AL and MJ designed and directed the execution
of the experiments. SG and MJ helped to carry out the clubroot
tests. SS provided project guidance and extensively edited the
manuscript. All authors contributed to manuscript revision, and
read and approved the submitted version

REFERENCES

Aigu, Y., Laperche, A., Mendes, J., Lariagon, C., Guichard, S., Gravot, A., et al. (2018). Nitrogen supply exerts a major/minor switch between two QTLs controlling *Plasmodiophora brassicae* spore content in rapeseed. *Plant Pathol.* 67, 1574-1581. doi: 10.1111/ppa.12867

Broman, K. W., Wu, H., Sen, S., and Churchill, G. A. (2003). R/qlt: QTL mapping in experimental crosses. *Bioinformatics* 19, 889–890. doi: 10.1093/bioinformatics/btg112

Buchanan, S. G. S. C., and Gay, N. J. (1996). Structural and functional diversity in the leucine-rich repeat family of proteins. *Prog. Biophys. Mol. Biol.* 65, 1–44. doi: 10.1016/S0079-6107(96)00003-X

Buczacki, S. T., Toxopeus, H., Mattusch, P., Johnston, T. D., Dixon, G. R., and Hobolth, L. A. (1975). Study of physiologic specialization in *Plasmodiophora brassicae*: proposals for attempted rationalization through an international approach. *Trans. Br. Mycol. Soc.* 65, 295–303. doi: 10.1016/S0007-1536(75) 80013-1

Chalhoub, B., Denoeud, F., Liu, S., Parkin, I. A. P., Tang, H., Wang, X., et al. (2014). Early allopolyploid evolution in the post-Neolithic *Brassica napus* oilseed genome. *Science* 345, 950–953. doi: 10.1126/science.1253435

Choi, H. W., Kim, Y. J., and Hwang, B. K. (2010). The hypersensitive induced reaction and leucine-rich repeat proteins regulate plant cell death associated with disease and plant immunity. *Mol. Plant Microb. Interact.* 24, 68–78. doi: 10.1094/MPMI-02-10-0030

Choi, S. R., Oh, S. H., Chhapekar, S. S., Dhandapani, V., Lee, C. Y., Rameneni, J. J., et al. (2020). Quantitative trait locus mapping of clubroot resistance and *Plasmodiophora brassicae* pathotype banglim-specific marker development in *Brassica rapa*. *Int. J. Mol. Sci.* 21:4157. doi: 10.3390/ijms21114157

Chu, M., Song, T., Falk, K. C., Zhang, X., Liu, X., Chang, A., et al. (2014). Fine mapping of *Rcr1* and analyses of its effect on transcriptome patterns during infection by *Plasmodiophora brassicae*. *BMC Genom.* 15:1166. doi: 10.1186/1471-2164-15-1166

Clarke, W. E., Higgins, E. E., Plieske, J., Wieseke, R., Sidebottom, C., Khedikar, Y., et al. (2016). A high-density SNP genotyping array for *Brassica napus* and its ancestral diploid species based on optimised selection of single-locus markers in the allotetraploid genome. *Theor. Appl. Genet.* 129, 1887–1899. doi: 10.1007/s00122-016-2746-7

Clough, S. J., Fengler, K. A., Yu, I. C., Lippok, B., Smith, R. K., and Bent, A. F. (2000). The *Arabidopsis dnd1* "defense, no death" gene encodes a mutated cyclic nucleotide-gated ion channel. *Proc. Natl. Acad. Sci. U.S.A.* 97, 9323–9328. doi: 10.1073/pnas.150005697

Dangl, J. L., and Jones, J. D. G. (2001). Plant pathogens and integrated defence responses to infection. *Nature* 411, 826–833.

FUNDING

The mapping population, the genotyping data, and the genetic
map used in this research were developed in the framework of
the national collaborative project entitled RAPOSDYN (ANR-
11-BTBR-0004) funded by the program “Investments for the
Future.” AB-R was awarded two mobility research training grants
under the ERASMUS + Program, one in 2018 and one in 2019.
The University of Alberta provided in-kind support.

ACKNOWLEDGMENTS

We thank the M2 students of Institut Agro/AGROCAMPUS
OUEST who participated in the phenotyping and carried out
the first runs of the analysis. We are grateful to Regine
Delourme for helpful discussions, especially considering the
QTL nomenclature.

Delourme, R., Falentin, C., Fomeju, B. F., Boillot, M., Lassalle, G., André, I., et al. (2013). High-density SNP-based genetic map development and linkage disequilibrium assessment in *Brassica napus* L. *BMC Genom.* 14:120. doi: 10.1186/1471-2164-14-120

Diederichsen, E., Beckmann, J., Schondelmeier, J., and Dreyer, F. (2006). Genetics of clubroot resistance in *Brassica napus* "Mendel". *Acta Hortic.* 706, 307–311.

Diederichsen, E., Deppe, U., and Sacristan, M. D. (2003). ".) Characterization of clubroot resistance in recent winter oilseed rape material," in *Proceedings of the 11th International Rapeseed Congress*, Copenhagen.

Diederichsen, E., Frauen, M., Linders, E. G. A., Hatakeyama, K., and Hirai, M. (2009). Status and perspectives of clubroot resistance breeding in crucifer crops. *J. Plant Growth Regul.* 28, 265–281. doi: 10.1007/s00344-009-9100-0

Diederichsen, E., Frauen, M., and Ludwig-Müller, J. (2014). Clubroot disease management challenges from a German perspective. *Can. J. Plant Pathol.* 36, 85–98. doi: 10.1080/07060661.2013.861871

Dixon, G. R. (2009). The Occurrence and economic impact of *Plasmodiophora brassicae* and clubroot disease. *J. Plant Growth Regul.* 28, 194–202. doi: 10.1007/s00344-009-9090-y

Duan, Y., Guo, J., Shi, X., Guan, X., Liu, F., Bai, P., et al. (2013). Wheat hypersensitive-induced reaction genes *TaHIR1* and *TaHIR3* are involved in response to stripe rust fungus infection and abiotic stresses. *Plant Cell Rep.* 32, 273–283. doi: 10.1007/s00299-012-1361-6

Ernst, T. W., Kher, S., Stanton, D., Rennie, D. C., Hwang, S. F., and Strelkov, S. E. (2019). *Plasmodiophora brassicae* resting spore dynamics in clubroot resistant canola (*Brassica napus*) cropping systems. *Plant Pathol.* 68, 399–408. doi: 10.1111/ppa.12949

Eulgem, T., Rushton, P. J., Robatzek, S., and Somssich, I. E. (2000). The WRKY superfamily of plant transcription factors. *Trends Plant Sci.* 5, 199–206. doi: 10.1016/S1360-1385(00)01600-9

Fähling, M., Graf, H., and Siemens, J. (2003). Pathotype separation of *Plasmodiophora brassicae* by the host plant. *J. Phytopathol.* 151, 425–430. doi: 10.1046/j.1439-0434.2003.00744.x

Fox, D. T., and Williams, P. H. (1984). Correlation of spore production by *Albugo Candida* on *Brassica campestris* and a visual white rust rating scale. *Can. J. Plant Pathol.* 6, 175–178. doi: 10.1080/07060668409501580

Frank, S. A. (1996). Models of parasite virulence. *Q. Rev. Biol.* 71, 37–78. doi: 10.1086/419267

Fredua-Agyeman, R., Hwang, S.-F., Strelkov, S. E., Zhou, Q., and Feindel, D. (2018). Potential loss of clubroot resistance genes from donor parent *Brassica rapa* subsp. *rapifera* (ECD 04) during doubled haploid production. *Plant Pathol.* 67, 892–901. doi: 10.1111/ppa.12816

Friberg, H., Lagerlöf, J., and Rämert, B. (2005). Germination of *Plasmodiophora brassicae* resting spores stimulated by a non-host plant. *Eur. J. Plant Pathol.* 113, 275–281. doi: 10.1007/s10658-005-2797-0
García-Olmedo, F., Molina, A., Alamillo, J. M., and Rodríguez-Palenzuélá, P. (1998). Plant defense peptides. *Biopolymers* 47, 479–491.

Gravot, A., Lemarié, S., Richard, G., Lime, T., Lariagon, C., and Manzanares-Dauleux, M. J. (2016). Flooding affects the development of *Plasmodiophora brassicae* in *Arabidopsis* roots during the secondary phase of infection. *Plant Pathol.* 65, 1153–1160. doi: 10.1111/ppa.12487

Hatakeyama, K., Niwa, T., Kato, T., Ohara, T., Kakizaki, T., and Matsumoto, S. (2017). The tandem repeated organization of NB-LRR genes in the clubroot-resistant *Crb* locus in *Brassica rapa* L. *Mol. Genet. Genom.* 292, 397–405. doi: 10.1007/s00438-016-1281-1

Hatakeyama, K., Suwabe, K., Tomita, R. N., Kato, T., Nunome, T., Fukuoka, H., et al. (2013). Identification and characterization of *Crr1a*, a gene for resistance to clubroot disease (*Plasmodiophora brassicae* Woronin) in *Brassica rapa* L. *PLoS One* 8:e054745. doi: 10.1371/journal.pone.0054745

Hejna, O., Havlickova, L., He, Z., Bancroft, I., and Curn, V. (2019). Analysing the genetic architecture of clubroot resistance variation in *Brassica napus* by associative transcriptomics. *Mol. Breed.* 39:112. doi: 10.1007/s11032-019-1021-4

Hirai, M., Harada, T., Kubo, N., Tsukada, M., Suwabe, K., and Matsumoto, S. (2004). A novel locus for clubroot resistance in *Brassica rapa* and its linkage markers. *Theor. Appl. Genet.* 108, 639–643. doi: 10.1007/s00122-003-1475-x

Holland, J. B., Nyquist, W. E., and Cervantes-Martínez, C. T. (2003). "Estimating and interpreting heritability for plant breeding: an update," in *Plant Breeding Reviews*, ed. J. Janick, (New York, NY: John Wiley & Sons Ltd), 9–112. doi: 10.1002/9780470650202.ch2

Hospital, F., and Charcosset, A. (1997). Marker-assisted introgression of quantitative trait loci. *Genetics* 147, 1469–1485.

Huang, Z., Peng, G., Gossen, B. D., and Yu, F. (2019). Fine mapping of a clubroot resistance gene from turnip using SNP markers identified from bulked segregant RNA-Seq. *Mol. Breed.* 39:131. doi: 10.1007/s11032-019-1038-8

Huang, Z., Peng, G., Liu, X., Deora, A., Falk, K. C., Gossen, B. D., et al. (2017). Fine mapping of a clubroot resistance gene in chinese cabbage using SNP markers identified from bulked segregant RNA sequencing. *Front. Plant Sci.* 8:1448. doi: 10.3389/fpls.2017.01448

Hwang, S. F., Ahmed, H. U., Zhou, Q., Rashid, A., Strelkov, S. E., Gossen, B. D., et al. (2013). Effect of susceptible and resistant canola plants on *Plasmodiophora brassicae* resting spore populations in the soil. *Plant Pathol.* 62, 404–412. doi: 10.1111/j.1365-3059.2012.02636.x

Hwang, S.-F., Howard, R. J., Strelkov, S. E., Gossen, B. D., and Peng, G. (2014). Management of clubroot (*Plasmodiophora brassicae*) on canola (*Brassica napus*) in western Canada. *Can. J. Plant Pathol.* 36, 49–65. doi: 10.1080/07060661.2013.863806

Jestin, C., Bardol, N., Lodé, M., Duffé, P., Domin, C., Vallée, P., et al. (2015). Connected populations for detecting quantitative resistance factors to phoma stem canker in oilseed rape (*Brassica napus* L.). *Mol. Breed.* 35:167. doi: 10.1007/s11032-015-0356-8

Kageyama, K., and Asano, T. (2009). Life cycle of *Plasmodiophora brassicae*. *J. Plant Growth Regul.* 28, 203–211. doi: 10.1007/s00344-009-9101-z

Kanzaki, H., Nirasawa, S., Saitoh, H., Ito, M., Nishihara, M., Terauchi, R., et al. (2002). Overexpression of the wasabi defensin gene confers enhanced resistance to blast fungus (*Magnaporthe grisea*) in transgenic rice. *Theor. Appl. Genet.* 105, 809–814. doi: 10.1007/s00122-001-0817-9

Karim, M., Dakouri, A., Zhang, Y., Chen, Q., Peng, G., Strelkov, S. E., et al. (2020). Two Clubroot-resistance genes, Rcr3 and Rcr9wa, mapped in *Brassica rapa* using bulk segregant RNA sequencing. *Int. J. Mol. Sci.* 21:5033. doi: 10.3390/ijms21145033

Kover, P. X., and Schaal, B. A. (2002). Genetic variation for disease resistance and tolerance among *Arabidopsis thaliana* accessions. *Proc. Natl. Acad. Sci. U.S.A.* 99, 11270–11274. doi: 10.1073/pnas.102288999

Laila, R., Park, J.-I., Robin, A. H. K., Natarajan, S., Vijayakumar, H., Shirasawa, K., et al. (2019). Mapping of a novel clubroot resistance QTL using ddRAD-seq in Chinese cabbage (*Brassica rapa* L.). *BMC Plant Biol.* 19:13. doi: 10.1186/s12870-018-1615-8

Laperche, A., Aigu, Y., Jubault, M., Ollier, M., Guichard, S., Glory, P., et al. (2017). Clubroot resistance QTL are modulated by nitrogen input in *Brassica napus*. *Theor. Appl. Genet.* 130, 669–684. doi: 10.1007/s00122-016-2842-8

Lecomte, L., Duffé, P., Buret, M., Servin, B., Hospital, F., and Causse, M. (2004). Marker-assisted introgression of five QTLs controlling fruit quality

traits into three tomato lines revealed interactions between QTLs and genetic backgrounds. *Theor. Appl. Genet.* 109, 658–668. doi: 10.1007/s00122-004-1674-0

Lee, J., Izzah, N. K., Choi, B.-S., Joh, H. J., Lee, S.-C., Perumal, S., et al. (2016). Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (*Brassica oleracea* L.). *DNA Res.* 23, 29–41. doi: 10.1093/dnares/dsv034

Lenth, R. V. (2016). Least-squares means: the R Package lsmeans. *J. Stat. Softw.* 69, 1–33. doi: 10.18637/jss.v069.i01

Li, L., Luo, Y., Chen, B., Xu, K., Zhang, F., Li, H., et al. (2016). A genome-wide association study reveals new loci for resistance to clubroot disease in *Brassica napus*. *Front. Plant Sci.* 7:1483. doi: 10.3389/fpls.2016.01483

Macfarlane, I. (1970). Germination of resting spores of *Plasmodiophora brassicae*. *Trans. Br. Mycol. Soc.* 55, 97–112. doi: 10.1016/S0007-1536(70)80100-0

Manzanares-Dauleux, M. J., Delourme, R., Baron, F., and Thomas, G. (2000). Mapping of one major gene and of QTLs involved in resistance to clubroot in *Brassica napus*. *Theor. Appl. Genet.* 101, 885–891. doi: 10.1007/s001220051557

Matsumoto, E., Yasui, C., Ohi, M., and Tsukada, M. (1998). Linkage analysis of RFLP markers for clubroot resistance and pigmentation in Chinese cabbage (*Brassica rapa* ssp. *pekinensis*). *Euphytica* 104:79. doi: 10.1023/A:1018370418201

Montarry, J., Corbiere, R., Lesueur, S., Glais, I., and Andrivon, D. (2006). Does selection by resistant hosts trigger local adaptation in plant-pathogen systems? *J. Evol. Biol.* 19, 522–531. doi: 10.1111/j.1420-9101.2005.01005.x

Montarry, J., Hamelin, F., Glais, I., Corbière, R., and Andrivon, D. (2010). Fitness costs associated with unnecessary virulence factors and life history traits: evolutionary insights from the potato late blight pathogen *Phytophthora infestans*. *BMC Evol. Biol.* 10:283. doi: 10.1186/1471-2148-10-283

Murakami, H., Tsuchima, S., Akimoto, T., Kuroyanagi, Y., and Shishido, Y. (2004). Quantitative studies on the relationship between plowing into soil of clubbed roots of preceding crops caused by *Plasmodiophora brassicae* and disease severity in succeeding crops. *J. Soil Sci. Plant Nutr.* 50, 1307–1311. doi: 10.1080/00380768.2004.10408609

Nagaoka, T., Doullah, M. A. U., Matsumoto, S., Kawasaki, S., Ishikawa, T., Hori, H., et al. (2010). Identification of QTLs that control clubroot resistance in *Brassica oleracea* and comparative analysis of clubroot resistance genes between *B. rapa* and *B. oleracea*. *Theor. Appl. Genet.* 120, 1335–1346. doi: 10.1007/s00122-010-1259-z

Neik, T. X., Barbetti, M. J., and Batley, J. (2017). Current status and challenges in identifying disease resistance genes in *Brassica napus*. *Front. Plant Sci.* 8:1788. doi: 10.3389/fpls.2017.01788

Nguyen, M. L., Monkhos, G. F., Komakhin, R. A., and Monakhos, S. G. (2018). The new clubroot resistance locus is located on chromosome A05 in Chinese cabbage (*Brassica rapa* L.). *Russ. J. Genet.* 54, 296–304. doi: 10.1134/S1022795418030080

Orgeur, G., Jestin, C., Delaunay, A., Lebourg, D., Bagot, P., Corbel, A., et al. (2016). Caractérisation des pathotypes de Hernie des crucifères en France et mise au point d'un test pour l'évaluation de la résistance des variétés de colza. *Innov. Agron.* 50, 145–155.

Pang, W., Fu, P., Li, X., Zhan, Z., Yu, S., and Piao, Z. (2018). Identification and mapping of the clubroot resistance gene crd in chinese cabbage (*Brassica rapa* ssp. *pekinensis*). *Front. Plant Sci* 9:653. doi: 10.3389/fpls.2018.00653

Parlevliet, J. E., and Zadoks, J. C. (1977). The integrated concept of disease resistance: a new view including horizontal and vertical resistance in plants. *Euphytica* 26, 5–21. doi: 10.1007/BF00032062

Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D., and R Core Team, (2019). *nlme: Linear and Nonlinear Mixed Effects Models*. Available online at: https://CRAN.R-project.org/package=nlme (accessed July 27, 2019).

Pink, D., and Puddephat, I. (1999). Deployment of disease resistance genes by plant transformation - a 'mix and match' approach. *Trends Plant Sci.* 4, 71–75. doi: 10.1016/S1360-1385(98)01372-7

R Core Team, (2019). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Rahman, H., Peng, G., Yu, F., Falk, K. C., Kulkarni, M., and Selvaraj, G. (2014). Genetics and breeding for clubroot resistance in Canadian spring canola (*Brassica napus* L.). *Can. J. of Plant Pathol.* 36, 122–134. doi: 10.1080/07060661. 2013.862571
Rashid, A., Ahmed, H. U., Xiao, Q., Hwang, S. F., and Strelkov, S. E. (2013). Effects of root exudates and pH on *Plasmodiophora brassicae* resting spore germination and infection of canola (*Brassica napus* L.) root hairs. Crop Prot. 48, 16–23. doi: 10.1016/j.cropro.2012.11.025

Robert, C., Bancal, M.-O., and Lannou, C. (2002). Wheat leaf rust uredospore production and carbon and nitrogen export in relation to lesion size and density. *Phytopathology* 92, 762–768. doi: 10.1094/PHYTO.2002.92.7.762

Rocherieux, J., Glory, P., Giboulot, A., Boury, S., Barbeyron, G., Thomas, G., et al. (2004). Isolate-specific and broad-spectrum QTLs are involved in the control of clubroot in *Brassica oleracea*. *Theor. Appl. Genet.* 108, 1555–1563. doi: 10.1007/s00122-003-1580-x

Rushton, P. J., Macdonald, H., Huttly, A. K., Lazarus, C. M., and Hooley, R. (1995). Members of a new family of DNA-binding proteins bind to a conserved cis-element in the promoters of α-Amy2 genes. *Plant Molec. Biol.* 29, 691–702. doi: 10.1007/BF00041160

Sacristán, S., Fraile, A., Malpica, J. M., and García-Arenal, F. (2005). An analysis of host adaptation and its relationship with virulence in cucumber mosaic virus. *Phytopathology* 95, 827–833. doi: 10.1094/PHYTO-95-0827

Sacristán, S., and García-Arenal, F. (2008). The evolution of virulence and pathogenicity in plant pathogen populations. *Mol. Plant Pathol.* 9, 369–384. doi: 10.1111/j.1364-3703.2007.00460.x

Saito, M., Kubo, N., Matsumoto, S., Suwabe, K., Tsukada, M., and Hirai, M. (2006). Fine mapping of the clubroot resistance gene, *Crr3*, in *Brassica rapa*. *Theor. Appl. Genet.* 114, 81–91. doi: 10.1007/s00122-006-0412-1

Sakamoto, K., Saito, A., Hayashida, N., Taguchi, G., and Matsumoto, E. (2008). Mapping of isolate-specific QTLs for clubroot resistance in Chinese cabbage (*Brassica rapa* L. ssp. *pekinensis*). *Theor. Appl. Genet.* 117, 759–767. doi: 10.1007/s00122-008-0817-0

Siemens, J., Nagel, M., Ludwig-Müller, J., and Sacristán, M. D. (2002). The interaction of *Plasmodiophora brassicae* and *Arabidopsis thaliana*: parameters for disease quantification and screening of mutant lines. *J. Phytopathol.* 150, 592–605. doi: 10.1046/j.1439-0434.2002.00818.x

Somé, A., Manzanares, M. J., Laurens, F., Baron, F., Thomas, G., and Rouxel, F. (1996). Variation for virulence on *Brassica napus* L. amongst *Plasmodiophora brassicae* collections from France and derived single-spore isolates. *Plant Pathol.* 3, 432–439.

Song, W.-Y., Wang, G.-L., Chen, L.-L., Kim, H.-S., Pi, L.-Y., Holsten, T., et al. (1995). A receptor kinase-like protein encoded by the rice disease resistance gene, *Xa21*. *Science* 270, 1804–1806. doi: 10.1126/science.270.5243.1804

Steele, K. A., Price, A. H., Shashidhar, H. E., and Witcombe, J. R. (2006). Marker-assisted selection to introgress rice QTLs controlling root traits into an Indian upland rice variety. *Theor. Appl. Genet.* 112, 208–221. doi: 10.1007/s00122-005-0110-4

Stonard, J., Downes, K., Pirie, E., Fitt, B., and Evans, N. (2007). “Development of phoma stem canker (*Leptosphaeria maculans*) and light leaf spot (*Pyrenopeziza brassicae*) on current and historical oilseed rape cultivars in 2003/04, 2004/05 and 2005/06 UK growing seasons,” in *Proceedings the 12th International Rapeseed Congress, Sustainable Development in Cruciferous Oilseed Crops Production*, Wuhan.

Stotz, H. U., Harvey, P. J., Haddadi, P., Mashanova, A., Kukol, A., Larkan, N. J., et al. (2018). Genomic evidence for genes encoding leucine-rich repeat receptors linked to resistance against the eukaryotic extra- and intracellular *Brassica napus* pathogens *Leptosphaeria maculans* and *Plasmodiophora brassicae*. *PLoS One* 13:e0198201. doi: 10.1371/journal.pone.0198201

Strelkov, S. E., Hwang, S.-F., Manolii, V. P., Cao, T., and Feindel, D. (2016). Emergence of new virulence phenotypes of *Plasmodiophora brassicae* on canola (*Brassica napus*) in Alberta, Canada. *Eur. J. Plant Pathol.* 145, 517–529. doi: 10.1007/s10658-016-0888-8

Strelkov, S. E., Hwang, S.-F., Manolii, V. P., Cao, T., Fredua-Agyeman, R., Harding, M. W., et al. (2018). Virulence and pathotype classification of *Plasmodiophora brassicae* populations collected from clubroot resistant canola (*Brassica napus*)

in Canada. *Can. J. Plant Pathol.* 40, 284–298. doi: 10.1080/07060661.2018.1459851

Suwabe, K., Tsukazaki, H., Iketani, H., Hatakeyama, K., Fujimura, M., Nunome, T., et al. (2003). Identification of two loci for resistance to clubroot (*Plasmodiophora brassicae* Woronin) in *Brassica rapa* L. *Theor. Appl. Genet.* 107, 997–1002. doi: 10.1007/s00122-003-1309-x

Suwabe, K., Tsukazaki, H., Iketani, H., Hatakeyama, K., Kondo, M., Fujimura, M., et al. (2006). Simple sequence repeat-based comparative genomics between *Brassica rapa* and *Arabidopsis thaliana*: the genetic origin of clubroot resistance. *Genetics* 173, 309–319. doi: 10.1534/genetics.104.038968

Ueno, H., Matsumoto, E., Aruga, D., Kitagawa, S., Matsumura, H., and Hayashida, N. (2012). Molecular characterization of the *Cra* gene conferring clubroot resistance in *Brassica rapa*. *Plant Mol. Biol.* 80, 621–629. doi: 10.1007/s11103-012-9971-5

Wagner, G., Laperche, A., Lariagon, C., Marnet, N., Renault, D., Guitton, Y., et al. (2019). Resolution of quantitative resistance to clubroot into QTL-specific metabolic modules. *J. Exp. Bot.* 70, 5375–5390. doi: 10.1093/jxb/erz265

Wallenhammar, A.-C. (1996). Prevalence of *Plasmodiophora brassicae* in a spring oilseed rape growing area in central Sweden and factors influencing soil infestation levels. *Plant Pathol.* 45, 710–719. doi: 10.1046/j.1365-3059.1996.01017-x

Werner, S., Diederichsen, E., Frauen, M., Schondelmaier, J., and Jung, C. (2008). Genetic mapping of clubroot resistance genes in oilseed rape. *Theor. Appl. Genet.* 116:363. doi: 10.1007/s00122-007-0674-2

Williams, P. H. (1966). A system for the determination of races of *Plasmodiophora brassicae* that infect cabbage and rutabaga. *Phytopathology* 56, 624–626.

Yousef, G. G., and Juvik, J. A. (2002). Enhancement of seedling emergence in sweet corn by marker-assisted backcrossing of beneficial QTL. *Crop Sci.* 42, 96–104. doi: 10.2135/cropsci2002.9600

Yu, F., Zhang, X., Huang, Z., Chu, M., Song, T., Falk, K. C., et al. (2016). Identification of genome-wide variants and discovery of variants associated with *Brassica rapa* clubroot resistance gene *Rcr1* through bulked segregant RNA sequencing. *PLoS One* 11:e0153218. doi: 10.1371/journal.pone.0153218

Yu, F., Zhang, X., Peng, G., Falk, K. C., Strelkov, S. E., and Gossen, B. D. (2017). Genotyping-by-sequencing reveals three QTL for clubroot resistance to six pathotypes of *Plasmodiophora brassicae* in *Brassica rapa*. *Sci. Rep.* 7:4156. doi: 10.1038/s41598-017-04903-2

Zamir, D. (2001). Improving plant breeding with exotic genetic libraries. *Nat. Rev. Genet.* 2, 983–989. doi: 10.1038/35103590

Zhang, H., Feng, J., Hwang, S.-F., Strelkov, S. E., Falak, I., Huang, X., et al. (2016). Mapping of clubroot (*Plasmodiophora brassicae*) resistance in canola (*Brassica napus*). *Plant Pathol.* 65, 435–440. doi: 10.1111/ppa.12422

Zhou, L., Cheung, M.-Y., Li, M.-W., Fu, Y., Sun, Z., Sun, S.-M., et al. (2010). Rice hypersensitive induced reaction protein 1 (*OsHIR1*) associates with plasma membrane and triggers hypersensitive cell death. *BMC Plant Biol.* 10:290. doi: 10.1186/1471-2229-10-290

Zhu, H., Zhai, W., Li, X., and Zhu, Y. (2019). Two QTLs controlling clubroot resistance identified from bulked segregant sequencing in pakchoi (*Brassica campestris* ssp. *chinensis* Makino). *Sci. Rep.* 9:9228. doi: 10.1038/s41598-019-44724-z

**Conflict of Interest:** 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 © 2020 Botero-Ramírez, Laperche, Guichard, Jubault, Gravot, Strelkov and Manzanares-Dauleux. 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.