Next Article in Journal
Role of the β-Catenin/REG Iα Axis in the Proliferation of Sessile Serrated Adenoma/Polyps Associated with Fusobacterium nucleatum
Next Article in Special Issue
Plasmodiophora brassicae Inoculum Density and Spatial Patterns at the Field Level and Relation to Soil Characteristics
Previous Article in Journal
Seroprevalence and Epidemiology of Toxoplasma gondii in Animals in the Qinghai-Tibetan Plateau Area, China
Previous Article in Special Issue
Molecular Pathotyping of Plasmodiophora brassicae—Genomes, Marker Genes, and Obstacles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Soil-Borne Inoculum of Plasmodiophora brassicae Measured by qPCR on Disease Severity of Clubroot-Resistant Cultivars of Winter Oilseed Rape (Brassica napus L.)

1
Rural Economy and Agricultural Society, HS Konsult AB, Gamla vägen 5G, SE-702 22 Örebro, Sweden
2
Rural Economy and Agricultural Society, HS Konsult AB, P.O. Box 412, SE-751 06 Uppsala, Sweden
3
Rural Economy and Agricultural Society, HS Konsult AB, Brunnby Gård, SE-725 97 Västerås, Sweden
4
RISE, Research Institutes of Sweden AB, Box 187, SE-532 32 Skara, Sweden
*
Author to whom correspondence should be addressed.
Pathogens 2021, 10(4), 433; https://doi.org/10.3390/pathogens10040433
Submission received: 17 March 2021 / Revised: 31 March 2021 / Accepted: 4 April 2021 / Published: 6 April 2021

Abstract

:
Use of resistant cultivars is considered the most effective tool in managing clubroot. Three clubroot-resistant commercial winter oilseed rape (OSR) cultivars and a susceptible ‘Cultivar mix’ were evaluated for disease severity index (DSI) and yield performance in field soils, selected for varying abundance of natural inoculum of Plasmodiophora brassicae. Seven field trials were carried out during 2017–2019 in winter OSR crops, and comparative bioassays were performed in a growth chamber. Substantial variation in clubroot infection between years was observed in the field trials. For Cultivar mix, a negative correlation (y = −252.3ln(x) + 58,897.6) was found between inoculum density and seed yield in five trials, whereas no correlation was found for the resistant cultivars. In bioassays, Cultivar mix exhibited a significantly high correlation between DSIb and number of gene copies g−1 soil (R2 = 0.72). For resistant cvs., Mentor and Alister, correlation was R2 = 0.45 and 0.58, respectively, indicating that resistance was under pressure. In field trials, DSIf of the resistant cultivars was lower (<27). The recommendation is thus to use clubroot-resistant cultivars of OSR as part of Integrated Pest Management in situations where abundance of P. brassicae DNA exceeds 1300 gene copies g−1 soil.

1. Introduction

Rapeseed is an important cornerstone in Swedish cropping, with production of rapeseed oil and rapeseed meal providing high-value food and feed constituents with good nutritional composition [1]. Brassica oilseed crops have been important in Swedish farming for the past 80 years [2] and winter oilseed rape (OSR) (Brassica napus) is currently a profitable crop, with sharply increased acreage in recent years following considerable damage by insect pests in spring OSR. In the main Swedish production areas, the acreage of winter OSR has increased from 90 to 94% of total oilseed rape acreage since 2015. However, increasing outbreaks of plant diseases that survive in the soil as resting spores or resting bodies, such as clubroot, Sclerotinia stem rot, and Verticillium wilt, have also been reported recently in Sweden and globally [3].
Clubroot is the most serious soil-borne disease of brassica oilseed and vegetable crops world-wide, causing appreciable yield losses [4]. The pathogen is an agricultural and biological challenge [5]. The causal agent, the obligate biotroph protist Plasmodiophora brassicae, is a member of the eukaryotic kingdom of Rhiazara, in the novel clade Phytorhiza [6]. The 25.5 Mb genome sequence of P. brassicae single spore isolate e3 was presented for the first time a few years ago [7]. Internationally, the spread of P. brassicae has increased rapidly in spring OSR (canola) in Canada [8,9], North Dakota [10], Finland [11] and Latin-America [12], and in winter OSR in the UK [13], Germany [14], Poland, [15,16], Czech Republic [17], China [18], and other countries.
In recent years, disease outbreaks in winter OSR have been frequently observed in commercial crops in Sweden, particularly in south and central parts of the country. These outbreaks, together with results from soil analysis in Scania (Sweden) [19] and in fields throughout the winter OSR area, show that the spread of P. brassicae has increased sharply. In DNA analyses by the Swedish Seed and Oilseed Growers’ Association on 210 farm samples collected in 2013–2015, P. brassicae DNA was detected in 49% of the fields sampled, but with large variation between fields and regions [20]. Recurring infections are now being recorded in areas in central Sweden where the prevalence of clubroot was high in the 1990s [2].
Control of clubroot is particularly difficult due to the long persistence of the pathogen in the soil [21], which impedes disease control by means of crop rotation. Breeding brassica crops for resistance has been the ultimate objective for many years, but progress is limited [22]. In Sweden, breeding for clubroot resistance in brassica oilseed crops started at Svalöf AB in the 1960s [23] and yielded one resistant cultivar, ‘Tosca’. However, high yield penalty resulted in its withdrawal when the clubroot-resistant cultivar cv. Mendel, with low yield penalty, was released in 2006 by NPZ in Germany. The resistance is based on the dominant gene and was initially effective against most field populations of P. brassicae [14]. However, as a consequence of frequent cropping of resistant cultivars in infested soil, varietal resistance is under pressure in some areas of the UK and has broken down in some cases [13]. Pathotypes capable of causing high levels of disease in resistant cultivars were reported in 2013 in the state of Alberta, Canada [24], with isolates found in these fields showing virulence to all resistant cultivars tested. Confirmed cases of virulent isolates infecting cv. Mendel are reported regularly, with most of the virulent isolates originating from North-east Germany [14]. A study in 2015 showed widespread occurrence in Germany of isolates of P. brassicae that can infect the resistant cultivar cv. Mendel [25].
Managing P. brassicae is one of the major challenges in rapeseed production and active measures must be taken, as the resting spores persist in the soil for up to 17 years [21], with an estimated half-life of 3.6 years. Spread of infected soil by machinery, wind [26,27] and water [28] has been demonstrated. The classic disease symptoms consist of an enlarged root that contains tens of millions of spores per gram of root. The resting spores are robust, with a well-developed mechanism involving various spore walls that provide protection against degradation by soil microorganisms. At high soil moisture and favorable temperature, the resting spores germinate, and zoospores are produced. These infect the root hairs and an intracellular primary phase takes place, followed by release of secondary zoospores [29]. The second phase of infection takes place in rhizodermal cells of the root and the pathogen subsequently induces local hypertrophy caused by mis-regulated plant-derived hormonal pathways [30]. The main part of the life cycle takes place well protected inside the root.
A DNA-based qPCR assay to analyze the presence of P. brassicae DNA in soil samples has been developed within BioSoM, a collaborative thematic project run by the Swedish University of Agricultural Sciences and the agricultural industry to devise biological mapping methods for soil-borne pathogens [31]. Using this assay, the abundance of P. brassicae DNA in soil samples can be determined at two commercial laboratories in Sweden, Eurofins Food and Testing Sweden AB (www.eurofins.se, accessed on 15 January 2021) and Intertek Scanbi Diagnostics AB (www.scanbidiagnostics.com, accessed on 15 January 2021). This fast and reliable qPCR technique has provided insights into the spread of the pathogen across Sweden, with soil analyses from around 30 variety trials in winter OSR and spring OSR showing contamination in 15% of fields tested [32]. Alarming spread of the disease has been demonstrated, with P. brassicae found in 60% of 45 fields on 18 farms in South-west Scania sampled in 2013 [19].
Breeding for resistant cultivars has expanded world-wide as clubroot disease has rapidly spread in oilseed rape production areas world-wide. Since 2012, several resistant winter OSR cultivars, developed based on the cv. Mendel resistance [33], have been available to Swedish OSR growers. However, these resistant cultivars are only partially resistant, designated “varietal resistance” in UK [13], as they are infected to some extent by even low levels of pathogen inoculum [34]. In Swedish winter OSR field trials evaluating clubroot-resistant cultivars harvested in 2014, resistant cultivars showed high yields but also substantial disease incidence, with 16–37% of resistant plants categorized as infected [35]. The highest disease incidence was found at a field trial site with very high abundance of P. brassicae DNA, corresponding to 10 million spores per g of soil, which implies considerable propagation of the pathogen in soil following a clubroot-resistant crop. There is a risk of overcoming this resistance [36], as has already occurred in fields with short OSR rotations in Canada.
Results from earlier Swedish studies on the agronomic performance of resistant breeding lines of spring oilseed turnip rape (Brassica rapa) are currently used for interpreting the results of soil DNA analysis and for providing guidance to growers. These studies showed that, at extremely high levels of inoculum in the soil, severe infection can occur in resistant breeding lines [34].
In the present study, performance of the resistance trait in commercial clubroot-resistant cultivars of winter OSR was evaluated in field soils with a natural inoculum of P. brassicae, to support choice of cultivar with respect to abundance of soil-borne inoculum. The aim was to develop an approach supported by DNA technology for integrated pest management (IPM) of clubroot in winter OSR crops. This was done by assessing infection levels and determining disease severity and yield of resistant and susceptible cultivars of winter OSR grown at field sites with differing abundance of P. brassicae DNA, and by performing comparative studies in a growth chamber in a controlled environment, in order to identify the prerequisites for maintaining sustainable production of winter OSR in fields where P. brassicae occurs. The hypothesis was that the selection pressure for new virulent pathotypes is considerably higher if resistant cultivars are grown in soils with confirmed abundance of P. brassicae DNA.

2. Results

2.1. Field Trials 2017–2018

2.1.1. Influence of Soil-Borne Inoculum of P. brassicae on Agronomic Performance of Winter OSR Cultivars

Four field trials were carried out in 2017. Two fields, located at Simrishamn and Tomelilla in southern Sweden (Scania), were selected for high abundance of P. brassicae DNA, with 1,100,000 and 2,500,000 gene copies g−1 soil, respectively. A field located at Kumla in central Sweden was selected for moderate abundance of disease inoculum, 15,000 gene copies per g−1 soil, and a field at Hallsberg in central Sweden for low levels of disease inoculum, 2500 gene copies g soil−1 (see Material and Methods). The three clubroot-resistant cultivars tested performed best at the two sites with the highest inoculum levels (Simrishamn, Tomelilla), with significantly higher seed yield (kg ha−1) than a susceptible Cultivar mix (p < 0.001) (Table S1). At the field trial site at Kumla with moderate inoculum level, statistically significant differences in seed yield were only obtained between Cultivar mix and clubroot-resistant cv. Alister (p = 0.017). A yield penalty was observed for all three resistant cultivars at the field site with the lowest inoculum level (Hallsberg), where all resistant cultivars produced significantly (p < 0.001) lower seed yield than the susceptible Cultivar mix (Table S1). Seed yield of Cultivar mix at all sites was correlated with the amount of gene copies of P. brassicae per g−1 soil, with seed yield decreasing rapidly when the number of gene copies g−1 soil increased above 4000 (y = −252.3ln(x) + 5897.6; R2 = 0.59) (Figure 1a). No significant correlations were observed for the resistant cultivars (Figure 1b). In terms of oil yield, all three clubroot-resistant cultivars had significantly higher oil yield (kg DM ha−1) than Cultivar mix at Tomelilla (p < 0.001), while at Kumla only cv. Alister differed significantly from Cultivar mix (p = 0.011) (Table S1).
Thousand-seed weight at Kumla was on average 5.7 g for Cultivar mix, 6.4 g for cv. Mentor, 5.7 g for cv. Alister, and 6.5 g for cv. Archimedes. Chlorophyll content was on average 11.2 ppm for Cultivar mix, 13.3 ppm for cv. Mentor, 11.3 ppm for cv. Alister and 14.4 ppm for cv. Archimedes. There was a significant difference (p = 0.014) in plant density between autumn count and spring count at the Tomelilla site, with 34% of Cultivar mix plants dying during winter. The winter-kill reduction in number of Cultivar mix plants was 25% at Simrishamn, 22% at Kumla, and 11% at Hallsberg (Table S1), but there were no significant differences in crop stand between the cultivars in either autumn or spring.

2.1.2. Visual Assessment of Clubroot Disease Severity and Disease Incidence 2017

Disease severity index (DSI) and disease incidence of clubroot varied between the cultivars and were also dependent on location. There was a significant (p < 0.001) effect of resistant cultivars on disease incidence and disease severity at the Simrishamn, Tomelilla and Kumla sites, where disease severity index (DSIf) of the resistant cultivars was 73–91% lower than that of the susceptible Cultivar mix (Table 1). There was a clear correlation between abundance of soil inoculum as indicated by real-time qPCR (gene copies g−1 soil) and DSIf for the susceptible Cultivar mix (y = 5.1678ln(x) − 2.2907; R2 = 0.42) (Figure 2a). Disease severity index (DSIf) rapidly increased as the number of gene copies g−1 soil increased above approximately 4000 gene copies g−1 soil. There was no correlation between gene copies g−1 soil and DSIf (R2 = 0.019) for the resistant cultivars and, even though the roots were infected at a soil inoculum level of 1499 gene copies per g−1 soil, DSIf values were below 30 (Figure 2b).
As regards disease incidence, 100% of plants were infected at the Kumla site, although the average number of gene copies determined was 55,600 g−1 soil (Table 1). Conducive conditions for disease infection and development prevailed at this field site, sown on 9 August (Figure S1). At the Simrishamn and Tomelilla sites, sowing was delayed due to rain, until 25 August and 26 August, respectively. The reduction in DSIf for resistant cultivars was statistically significant (p < 0.001) compared with Cultivar mix (S) and ranged from 77% to 90% at Simrishamn and from 85% to 94% at the Tomelilla site. Disease incidence was high at these sites, with 68% and 63%, respectively, of Cultivar mix plants infected. There were no statistically significant differences in DSIf between the resistant cultivars and Cultivar mix at the Hallsberg site.

2.2. Field Trials 2018–2019

2.2.1. Influence of Clubroot on Agronomic Performance of Winter OSR Cultivars

In 2018–2019, three field trials were carried out in fields at Simrishamn and Tomelilla with very high abundance of P. brassicae DNA (600,000 and 370,000 gene copies g−1 soil, respectively, according to pre-testing) and in a field at Kumla with moderate abundance (50,000 gene copies per g−1 soil) (Table S2). Significantly (p < 0.001) lower yield was observed for cv. Archimedes than for the other cultivars at the Kumla site. At the Simrishamn site, significantly higher yield was obtained for Cultivar mix and cv. Mentor than for cv. Alister (R) (p < 0.001). There was no significant difference in yield between cultivars at the Tomelilla site. Plant counts to estimate winter damage showed no change in plant density between autumn count and spring count at the Simrishamn and Tomelilla sites, whereas there was a significant difference (p < 0.001) in spring count at the Kumla field trial site for cv. Archimedes, with a reduction in plant numbers of 56%. The plant reduction for the other cultivars was on average 17%.

2.2.2. Visual Assessment of Clubroot Severity and Clubroot Incidence

Assessment of disease severity showed no significant difference between cultivars in DSIf in autumn, which ranged from 22.5 for Cultivar mix to 3.3 for cv. Mentor. The assessments indicated that cv. Archimedes suffered from winter damage caused by abiotic factors (Table 2). There was no difference in clubroot severity at the Simrishamn and Tomelilla sites, as DSIf was extremely low (3.3) for Cultivar mix at Simrishamn and no disease was seen at Tomelilla despite the high level of disease inoculum (Table 2). All field trials were assessed for disease a second time, in the stubble after harvest in July 2019, when only 2% and 1% of cv. Mentor plants assessed were infected at Simrishamn and Tomelilla, respectively. Non-conducive conditions prevailed at those sites in 2018 (see Figure S1), as precipitation was low in August and September, while precipitation at the Kumla site after sowing was 39 mm (25–31 August) (Figure S1). All cultivars at Kumla were infected, but no significant difference in DSIf was seen in the autumn assessment, whereas the assessment after harvest showed a statistically significant difference in infection (p = 0.016) between Cultivar mix (DSIf 27.5) and cv. Alister (DSIf 1.7) and cv. Archimedes (DSIf 5.3). Average values from both assessments showed significantly (p < 0.001) lower DSIf values for the resistant cultivars, 76–85% lower than for the susceptible Cultivar mix (Table 2). Regarding disease incidence, on average 34% of the susceptible plants were infected, while there was significantly lower disease incidence (p < 0.004) in the resistant cultivars.
A high correlation between disease severity (DSIf) and seed yield was found for the susceptible Cultivar mix (y = −19.683x + 4123.4; R2 = 0.66) (Figure 3a). At lower levels of disease severity index (<20), the yield reduction was about 390 kg ha−1, while yield was reduced by 2155 kg ha−1 (55%) at DSIf close to 100. No correlation was found for the cultivars defined as resistant (Figure 3b). Likewise, there was a significant negative correlation between disease incidence and yield for Cultivar mix (y = −19.013x + 4300.1, R2 = 0.66), while no correlation was found for cv. Mentor, cv. Alister, or cv. Archimedes.

2.3. Soil Bioassays and Real-Time qPCR Analysis

All cultivars in the field trials were also grown plot-wise in bioassays in a growth chamber under controlled conditions. Real-time qPCR analysis was performed on soil from each plot in the field trials (n = 16) (Table 3 and Table 4). In 2017, there was a statistically significant difference (p < 0.001) between the field sites in average number (n = 4) of gene copies per g−1 soil. The number was significantly highest at the Tomelilla site (1,187,100), followed by Simrishamn (419,200), whereas no significant difference was seen between the Kumla and Hallsberg sites (Table 3). No significant differences were found between cultivars (plots), irrespective of field site or number of gene copies per g−1 soil. In 2018, a significant difference (p < 0.001) between field sites emerged, with the average number (n = 4) of gene copies per g−1 soil being significantly highest at the Simrishamn site (442,300) (Table 4). There was no significant difference in pathogen abundance between Tomelilla and Kumla soil (38,000 and 3,700 gene copies per g−1 soil, respectively). Within the field sites, there was a statistically significant difference in pathogen abundance between samples from plots at Simrishamn (p = 0.015) and Tomelilla (p < 0.001).
The results of the bioassays in 2017 are displayed in Table 3. Under high disease pressure and conducive conditions, considerably higher disease severity index (DSIb) and disease incidence (DIb) were seen compared to DSIf for soil samples from Simrishamn and Tomelilla, where DSIb increased by 48% and 36%, respectively, for Cultivar mix and by 270% and 287%, respectively, for the resistant cultivars. In the latter, DSI exceeded 30, which is considered the limit at which resistance is overcome [36]. In soil from the Kumla site, disease severity in the bioassay was 50% lower than that estimated in the field. For Hallsberg samples, disease severity was at a very low level for all cultivars in the bioassay and no statistically significant difference was seen.
The relationship between DSIb and abundance of gene copies in soil for all cultivars tested is shown in Figure 4. A significantly high correlation between DSIb and gene copies was found for Cultivar mix, (R2 = 0.7165) (Figure 4a), and for cv. Mentor (R2 = 0.4477) (Figure 4b) and cv. Alister (R2 = 0.5771) (Figure 4c). The correlation for cv. Archimedes was low (R2 = 0.23) (Figure 4d). Considering these cultivars individually at a particular DSIb level, e.g., 20, using the equations in Figure 4a, it was found that it only required 7600 gene copies for Cultivar mix to reach DSIb 20, while it required 1,480,000 gene copies for cv. Mentor, 77,000 gene copies for cv. Alister and approximately 960,000 gene copies for cv. Archimedes. This corresponded to 340,000, 177,000, and 2,200,000 spores per g−1 soil, respectively, in the most conducive conditions.
There was a statistically significant correlation between DSIb and number of gene copies for Cultivar mix (p = 0.035), cv. Mentor (p < 0.001) and cv. Alister (p = 0.012) in soil from the Simrishamn site, with a higher number of gene copies in the soil increasing the DSI values. For cv. Archimedes there was no correlation (Table 3). In soil from Kumla, positive correlations were found for cv. Mentor (p = 0.011) and cv. Archimedes (p = 0.017). There was no positive correlation between DSIb and number of gene copies in soil for any cultivar at any of the sites studied in 2019 (Table 4).

3. Discussion

Clubroot is a serious threat to OSR production world-wide. The increasing proportion of clubroot-infected arable land in Sweden has serious consequences for oilseed rape growers [2]. OSR crops are valuable constituents of crop rotations and increasing market prices have resulted in production expanding in recent years, often through increased OSR frequency within rotations in Sweden and other countries [9,13,37,38]. Clubroot-resistant winter OSR cultivars have been available for several years, but the market share is moderate. One constraint is the yield penalty of clubroot-resistant cultivars [35], which was 6–10% relative to susceptible Cultivar mix in the present study when grown in fields with low or no inoculum (Table S1). In this study, the agronomic performance of clubroot-resistant winter OSR cultivars was assessed in the field on soils with different abundance of P. brassicae DNA and in comparative bioassay tests in a growth chamber under conducive conditions for the pathogen.
The severity of clubroot disease varied between the field sites (Table 1), which were selected to represent different levels of disease inoculum determined by real-time qPCR assay [19]. There was also substantial variation in clubroot infection between years (Table 1 and Table 2).
High temperatures in the 2018 season (e.g., 30–37 °C at Kumla) also affected evaporation, and therefore yield differences due to clubroot was probably strongly masked (Table S1), as experienced by others [3]. Yield in the clubroot-resistant cultivars at Kumla in 2018 was only 10% higher than for the susceptible Cultivar mix, despite severe clubroot infection (Table 1), compared with 35% at Simrishamn and 55% at Tomelilla. Winter OSR sowing date (9 August at Kumla, 25–26 August at Simrishamn and Tomelilla), is likely to have played an important role in the results, as the robustness of surviving plants in spring is crucial for the yield outcome (A. Gunnarsson, Svensk Raps, personal communication 2010). For the susceptible Cultivar mix (cvs. Avatar, Dariot, Explicit, Exstorm) in 2017–2018, a negative correlation (y= −252.3ln(x) + 58,897.6) was found between inoculum abundance and seed yield (Figure 1a). Yield losses caused by clubroot are often not observed in winter OSR [39], but the results in Figure 2 clearly show that at an inoculum abundance of 5000 gene copies g−1 soil (corresponding to 12,000 spores g−1 soil), estimated yield loss was 8%, while at an abundance of 850,000 gene copies g−1 soil (corresponding to 2 million spores g−1 soil), estimated yield loss was 44%. For the resistant cultivars, no correlation was found between inoculum abundance and seed yield. Thus, differences in seed yield were due rather to the agronomic performance of the cultivars and the prevailing biotic and abiotic conditions at each site (Figure 1b). The yield reductions observed were in line with those in a field trial on spring OSR (B. napus) in central Sweden, where a yield loss of 50% was found at the highest level of disease incidence [39].
There was a clear correlation between amount of soil inoculum (gene copies g−1 soil) as indicated by real-time qPCR and DSIf for the susceptible Cultivar mix (Figure 2a). DSIf rapidly increased with increasing pathogen abundance above approximately 4000 gene copies g−1 soil (corresponding to 9000 spores g−1 soil). There was no such correlation for the resistant cultivars (R2 = 0.0192), even though these cultivars were infected at a higher rate, particularly at inoculum densities >70,000 gene copies g−1 soil (corresponding to 160,000 spores g−1 soil). The DSIf recorded was close to 30 (Figure 2b), indicating that cultivar resistance was under pressure [40] and confirming our initial hypothesis. The clubroot-resistant cultivars showed partial, rather than complete, resistance, as found in earlier studies on summer oilseed turnip rape [34]. The results also support recent findings in a greenhouse set-up that increasing inoculum load applied to the resistant cv. Mendel increases the infection rate [33]. Canadian studies on clubroot-resistant OSR cultivars have found increased inoculum loads (quantified by qPCR) after cultivation and a peak in P. brassicae DNA in the following year [41].
Disease development of clubroot is strongly influenced by environmental factors, with soil water content long being recognized as a major abiotic factor contributing to the development of clubroot [42]. Disease symptoms can develop after 10–18 h in excessively wet soil [43]. When soil moisture content exceeds 50% soil water-holding capacity, the disease develops very quickly [4]. Disease assessments at the different field sites in this study clearly showed the influence of soil moisture on infection (Table 1 and Table 2). Despite very high abundance of P. brassicae DNA at the Simrishamn and Tomelilla sites in southern Sweden (600,000 and 370,000 gene copies g−1 soil, respectively, corresponding to 1.2 million and 740,000 spores g−1 soil), disease symptoms were absent in autumn 2018.
The extremely low level of clubroot infection at the sites in southern Sweden in 2018 was associated with exceptionally low soil moisture content after sowing and throughout September (Figure S1). Significantly lower seed yield at Simrishamn was recorded for cv. Alister and cv. Archimedes in 2019, but this was not related to clubroot disease as no infection occurred (Table 2). It was probably due to other soil-borne diseases such as Verticillium wilt, which was prevalent but not assessed. The chlorophyll content was higher at that site (p < 0.001), revealing uneven maturity (Table S2). Seed yield at the Kumla site in 2019 was significantly lower for cv. Archimedes than for all other cultivars. This difference was likely related to winter damage due to vigorous growth in autumn, as plant density in spring was reduced by 56%. Uneven maturity was reflected in significantly higher chlorophyll content (p = 0.007).
Temperature is an abiotic factor of great importance for successful clubroot infection. Early studies found that the conducive temperature range for infection is 12–27 °C [44], with galls developing only above 9 °C [44]. More recent studies show that air temperature has a significant effect on clubroot incidence in infected soil, with 40-fold lower risk of clubroot infection at 12 °C compared with 22 °C [11], and with higher clubroot infection at 25 °C compared with 20 °C, particularly at lower pH levels [45].
Winter OSR sown in early August often experiences air temperatures within the conducive range for clubroot infection. In 2017, the field trial site at Kumla was sown when the average air temperature was 20 °C, followed by plentiful precipitation (Figure S1). This resulted in all plants of the susceptible Cultivar mix being infected, with a corresponding DSIf of 92.3.
Recent studies on the survival of P. brassicae spores in different temperature regimes show that spores buried in the soil at 4 °C and 20 °C show better survival than spores on the soil surface at −20 °C and 30 °C [46], also ultraviolet (UV) light has been shown to have a negative effect on spore viability.
In the present study, bare fallow preceded the winter OSR crops at the Kumla site in both years, thus exposing the soil surface to UV radiation and absorbing heat. This heat increment may have been detrimental for P. brassicae resting spores close to the soil surface, as bioassays and field qPCR measurements showed low levels of P. brassicae inoculum (Table 4) relative to number of infected plants in Kumla plots in late autumn (Table 2).
It is evident that, through development of robust resting spores, P. brassicae is well-equipped to survive through many seasonal cycles [47] and in hostile environments. The pH at the field sites ranged from 6.1 to 6.9 in topsoil and from 6.8 to 7.4 at 60–90 cm depth (Table S1), confirming previous observations [48]. However, there was no evident influence of pH on abundance of P. brassicae inoculum in soil. Soil water content appeared to be the most important factor for germination, motility of the resting spores and further steps in the infection cycle. During dry seasons, infection was impeded despite a spore load of about 1.4 million spores g soil−1 (Table 2). A second assessment after harvest in 2019 (Table 2) confirmed that roots were most vulnerable to infection at the seedling stage [47], as clubroot disease incidence was similar as in late autumn 2018. The infection potential is latent and can lead to a false belief among growers that clubroot is no longer a problem, so soil testing for clubroot infection potential in field soil is crucial for determining choice of cultivar. Use of resistant cultivars is a major tool in IPM [49], but based on findings in the present study these cultivars must be used in environments where infection pressure is low to moderate, as otherwise the resistance is put under pressure (Figure 4b–d). Previous recommendations for Sweden state that an abundance of >325,000 gene copies g−1 soil poses a high risk of considerable multiplication of spores in clubroot-resistant cultivars [19]. The results in this study indicate that this threshold should be reduced to 100,000 gene copies g−1 soil to avoid a high risk of pathogen multiplication in conducive conditions.
Populations of a new virulence phenotype of the pathogen capable of overcoming resistance in clubroot-resistant cultivars of OSR have been identified in areas where crop rotations are short [13,24]. Recent studies report severe damage from a virulent isolate on clubroot-resistant cv. Mendel at low spore loads [33]. In bioassays in the present study on soil samples from sites where conditions for pathogen infection were optimal, the DSIb values revealed that resistance was under pressure at high inoculum abundance. Thus, virulent isolates infecting resistant cultivars are likely to be prevalent in the Swedish growing districts represented. The physiological specialization of P. brassicae at these field trial sites requires urgent investigation.
Under field conditions, DSIf was lower than DSIb, but close to the level where resistance of the cultivars tested was under pressure. The influence of high temperatures (<30 °C) on spore survival in the soil and indications of an impact of UV radiation mean that soil sampling to predict the inoculum potential of P. brassicae should be undertaken in cool season conditions. Low levels of disease did not seem to affect yield significantly in susceptible cultivars (Figure 3a). Resistant cultivars may achieve high yield despite high inoculum levels, resulting in multiplication of inoculum and putting resistance under pressure. It is therefore critically important for growers to be proactive and use IPM clubroot control tools such as soil tests based on DNA technology and clubroot-resistant cultivars. This approach will help OSR growers in Sweden and elsewhere to overcome the challenges and maintain OSR cultivation in an environment of global warming and climate change.

4. Materials and Methods

4.1. Field Trial Assessment of Resistant Cultivars of Winter OSR

Four field trials were established in August 2017 and three field trials in August 2018 in winter OSR fields on commercial farms. Sites were selected based on the results of soil testing for prevalence of P. brassicae DNA, which was performed at a commercial laboratory, Eurofins Sweden Testing Agri AB, Kristianstad, Sweden (www.eurofins.se, accessed on 14 January 2021). The field sites were in regions with a high risk of clubroot infection and represented a soil content of P. brassicae DNA defined as low (<1300 gene copies per g−1 soil), moderate (1300–50,000 gene copies per g−1 soil), high (50,000–325,000 gene copies g−1 soil), and very high (>325,000 gene copies g−1 soil), as described in Wallenhammar et al. [19]. Information on the field sites is displayed in Table S3. Crops were managed by the farmers according to common practices in the region, with application of nitrogen, other nutrients, herbicides, and insecticides. A dose of Fe-III-phosphate (Sluxx®®HP, 5 kg ha−1) was applied after sowing to prevent snail damage. Fungicide was applied by the field trial managers at growth stage BBCH 65 to ensure that Sclerotinia stem rot did not reduce yield. Growth regulators were not used. The trials were performed in a randomized block design, with four replicates, 24–30 m2 plots. The clubroot-resistant winter OSR cultivars cv. Mendel, cv. Alister and cv. Archimedes and a control mixture of four susceptible cultivars (Cultivar mix) were sown annually at a rate of 50 germinable seeds m2 and a row spacing of 12 cm. The susceptible Cultivar mix was similar to that used in the Swedish official cultivar trials. In 2017–2018, the cultivars were cv. Avatar, cv. Dariot, cv. Explicit and cv. Exstorm, while in 2018–2019 the cultivars were cv. Avatar, cv. Explicit, cv. Mercedes and cv. Harnas. The seeds were treated with Crusier®® (Thiametoxam 15 mL kg−1) in an industrial batch seed treater (Satec Concept®® ML2000, Agritema, Kiev, Ukraine) at the Department of Seed Technology, Rural Economy and Agricultural Society, Bjärred, Sweden. The seed was kindly provided by the seed companies. The plots were sown at the latest by mid-August to ensure that soil temperatures favored infection. The sowing dates in 2017 were 9 August at Kumla, 15 August at Hallsberg, 25 August at Simrishamn and 26 August at Tomelilla. The dates in 2018 were 15 August at Kumla, 17 August at Simrishamn and 21 August at Tomelilla.
Total harvested area for each plot was 18–22 m2. The plots were harvested with a plot combine harvester, and stage of maturity, chlorophyll content, yield, oil content and thousand seed weight were determined. Seed purity, water content and oil content were used for calculation of yield (9% wc), oil yield, and oil content (NIT-method ISO 12099) at 9% water content (ISO 12099). Plant density was graded in late autumn and in April the following year. Stage of maturity was graded prior to harvest. Thousand-seed weight and chlorophyll content (spectrometric method ISO 10519:1997) were only measured at Kumla in 2018, but at all field sites in 2019. Meteorological data (temperature and rainfall) for each field site were obtained from Lantmet (www.ffe.slu.se/lm/LMHome.cfm?LMSUB=1, accessed on 16 January 2021). Meteorological conditions during sowing and germination of winter OSR varied between the southern and central regions of Sweden where the field trials were carried out in 2017 and 2018, particularly regarding precipitation (Figure S1).

4.2. Soil Sampling and Preparation

4.2.1. Soil Sampling for Selecting Trial Fields

Soil samples to 20 cm depth were collected at 40 randomly selected points along the arms of a “W” sampling pattern, using a soil auger with diameter 22 mm and volume 76 mL (Galko Svets and Lantbruk, Heberg, Sweden). The cores were mixed in a bucket to a composite sample, and handled according to instructions from the laboratory (https://www.svenskraps.se/se/pdf/klumprotsjuka_2020-05-28_Intertek_Analys--Klumprotsjuka_brochure_april%202020.pdf, accessed on 15 January 2021)

4.2.2. Soil Sampling for Determination of Physiochemical Parameters

For determination of physiochemical parameters across the experimental area, 20 soil cores were taken in the top layer (0–30 cm) and 10 cores each in the 30–60 cm and 60–90 cm layers (Table S3).

4.2.3. Soil Sampling for Bioassays

Immediately prior to sowing, 25–30 soil cores were taken randomly from the top 20 cm of the soil profile and pooled plot-wise for bioassay. The cores were mixed in a bucket into a composite sample of approximately 2 L for the plot and transferred to plastic bags. The soil samples were kept in cold storage until arrival at the laboratory, where the soil was poured into foil tins, clods were macerated by hand and stones were removed prior to drying at room temperature to 5% water content. The dried soil was then stored at 10 °C until the bioassay was carried out.
For the bioassay, the composite sample was poured into a plastic bag closed with a cable tie and placed in an outer plastic bag closed with a bag sealer. The bag was then run in a cement mixer (MEEC 63 l, 220 W, 27.5 rpm) for 5 min and 500 g of the soil were weighed into a 1.0 L plastic container (HPPE-plastic), to which three steel balls (diameter 12 mm) and three screw nuts (14 mm) were added. The container was closed with a lid and shaken in a paint mixing apparatus (Svenska Skandex AB) for 45 s at 1000 rpm. Then 50 g of each mixed sample were transferred to a 50 mL Falcon tube, from which subsamples were used for DNA analysis. The rest of the soil was poured back into the composite sample.

4.2.4. Determination of Soil Chemical Properties

Measurements of selected physiochemical parameters in the soils were performed by Euro Eurofins Agro Testing Sweden AB., Kristianstad, Sweden. Soil pH was determined potentiometrically in water (1.0:2.5, wt/wt).

4.3. Soil Bioassay

All cultivars used in the field trials were tested in bioassays. The composite soil sample was divided into four pots (Göttinger, 9 × 9 × 9.5, volume 0.5 L). The bioassays were carried out in a growth chamber in a completely randomized design replicated four times and maintained for a six-week period, to ensure optimal infection according to Wallenhammar [31]. Fifteen seeds per pot of each resistant cultivar or the susceptible Cultivar mix were sown and thinned out to 10 seedlings at true leaf stage. Infected plants were scored according to the following: 0 = no galls; 1 = enlarged lateral roots; 2 = enlarged taproot; 3 = enlarged napiform taproot; 4 = enlarged napiform taproot, lateral roots healthy; 5 = enlarged napiform taproot, lateral roots infected. Disease severity index (DSI) was calculated as
DSI =     Class   no   ×   No .   of   plants   in   each   class ( Total   no .   of   plants )   ×   ( No .   classes     1 )

4.4. Soil DNA Extraction and Real-Time PCR Analysis

DNA was extracted from 350 mg soil samples (two replicates) using the FastDNA®® spin kit for soil (MP Biomedicals, Solon, OH, USA) according to the manufacturer’s instructions and the modifications reported by Wallenhammar et al. [31]. DNA extracts were eluted in 100 µL DES (DNase/Pyrogen-free water) and further purified once using Wizard DNA Clean-Up System (Promega Corporation, Fitchburg, WI, USA) and three times using Illustra MicroSpin S-300 HR Columns (Cytiva, Marlborough, MA, USA), according to the manufacturers’ instructions.
The DNA samples were analyzed by quantitative real-time PCR (qPCR) using the 7300 Real Time PCR System (Applied Biosystems) and quantified as described by Wallenhammar et al. [31]. DNA from each of the two duplicates was analyzed twice in each qPCR reaction. DNA inhibition was assessed by adding 2 µL of P. brassicae standard (150 gene copies) to an extra reaction for one of the duplicate DNA samples and an extra reaction for the negative water control. qPCR analysis performed at the DNA laboratory, Dept. of Soil and Environment, SLU, Skara.

4.5. Field Assessment of Disease

An evaluation of disease was performed on 25 plants per plot, which were uprooted in mid-November, placed in plastic bags and brought to the laboratory, where the roots were carefully rinsed in running tap water and assessed for infection. The roots were scored according to the following: 0 = no galls; 1 = slight galls on lateral roots; 2 = moderate galls (<50% of the root system galled); 3 = severe galls (>50% of the root system galled). Field disease severity index (FDSI) was calculated according to Wallenhammar et al. [34], using the equation presented in Section 4.3. In 2019, an additional disease evaluation was undertaken in the stubble after harvest, on 25 plants randomly uprooted and handled as above.

4.6. Statistics

The results were processed using ANOVA tests followed by Tukey’s HSD-test (p < 0.05) in JMP 15.2 (SAS Inst. Inc. Cary, NC, USA). Interactions between the logarithm of number of gene copies g−1 soil and yield, disease severity index and disease incidence for the cultivars in the field trials and bioassays were calculated using regression analysis in JMP 15.2. All data were used in the regression analysis for the bioassays.

5. Conclusions

Managing clubroot once it establishes in OSR fields is very difficult due to the persistence of the pathogen in soil. Commercial cultivars with clubroot resistance are available from most seed companies and provide high levels of clubroot control if grown in fields where the pathogen inoculum is low, as determined by soil analysis. The response of resistant cultivars to disease development constitutes the basis for current guidance to growers on the best long-term control strategy for clubroot. New resistant cultivars exhibit higher yield potential, in the same range as susceptible cultivars, and are thus an efficient tool for growers. Soil testing based on DNA technology needs to be used to a greater extent to avoid outbreaks in fields cropped with susceptible cultivars, and consistent loss of valuable seed yield.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/pathogens10040433/s1, Figure S1: Meteorological information at the field trial sites, Table S1: Yield parameters at harvest 2018 for the four field trial sites, Table S2: Yield parameters at harvest 2019 for the three field trial sites, Table S3: Location of field sites and soil parameters.

Author Contributions

Conceptualization, A.-C.W. and A.J.; methodology, A.-C.W., A.J., and Z.S.O.; validation, A.-C.W., A.J., and Z.S.O.; investigation, A.-C.W. and A.J.; resources, A.-C.W.; data curation, E.E., Z.S.O., A-C.W., and A.J.; writing—original draft preparation, A.-C.W, E.E., and Z.S.O.; writing—review and editing, A.-C.W, E.E., Z.S.O., and A.J.; project administration, A.-C.W.; funding acquisition, A.-C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Swedish Farmers’ Foundation for Agricultural Research (grant number O-16-20-765), the Foundation for Swedish Oil Plant research and Carl Robert Prytz Donation fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors thank the field experimental units at the Rural Economy and Agricultural Society (Skåne and Örebro branches) for technical support and assistance. We thank Bo Gertsson, Lantmännen Lantbruk, Mats Olsson Sörensson, Syngenta, and Ola Sixtensson, Scandinavian Seed, for providing seeds of resistant cultivars and Albin Gunnarsson, Swedish Seed and Oilseed Growers, for providing ‘Cultivar mix’. We thank Tomas Johansson, RISE, for assisting with the bioassays and the commercial farms for hosting the field trials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ackman, R.D. Canola Fatty Acids- An Ideal Mixture for Health, Nutrition and Food Use. In Canola and Rapeseed: Production, Chemistry, Nutrition and Processing Technology, 1st ed.; Shahidi, F., Ed.; Van Nostrand Reinhold: New York, NY, USA, 1990; pp. 81–98. [Google Scholar]
  2. Wallenhammar, A.-C.; Almquist, C.; Schwelm, A.; Roos, J.; Marcez-Schmidt, K.; Jonsson, A.; Dixelius, C. Clubroot, a persitent threat to Swedish oil seed rape production. Can. J. Plant Pathol. 2014, 36, 135–141. [Google Scholar] [CrossRef]
  3. Zheng, X.; Koopman, B.; Ulber, B.; von Tiderman, A. A Global Survey on Diseases and Pests in Oilseed Rape- Current Challenges and Innovative Strategies of Control. Front. Agron. 2020, 2, 15–21. [Google Scholar] [CrossRef]
  4. Dixon, G.R. The occurrence and economic impact of Plasmodiophora brassicae and clubroot disease. J. Plant Growth Regul. 2009, 28, 194–202. [Google Scholar] [CrossRef]
  5. Dixon, G.R. Clubroot (Plasmodiophora brassicae Woronin)- an agricultural and biological challenge world-wide. Can. J. Plant Pathol. 2014, 36, 5–18. [Google Scholar] [CrossRef]
  6. Sierra, R.; Canas-Duarte, S.J.; Burki, F.; Schwelm, A.; Fogelqvist, J.; Dixelius, C.; González-Garcia, L.N.; Gile, G.H.; Slamovits, C.; Restrepo, S.; et al. Evolutionary origins of Rhiazaran parasites. Mol. Biol. Evol. 2015, 33, 980–983. [Google Scholar] [CrossRef] [Green Version]
  7. Schwelm, A.; Fogelqvist, J.; Knaust, A.; Julke, S.; Lilja, T.; Bonilla-Rosso, G.; Karlsson, M.; Shevchenko, A.; Dhanapandi, V.; Ryun Choi, S.; et al. The Plasmodiophora brassicae genome reveals insightis its life cycle and ancestry of chitin synthases. Nat. Sci. Rep. 2015, 5, 1–12. [Google Scholar]
  8. Strelkov, S.E.; Hwang, S.-F. Clubroot in the Canadian Canola crop: 10 years into the outbreak. Can. J. Plant Pathol. 2014, 36, 27–36. [Google Scholar] [CrossRef]
  9. Stelkov, S.E.; Hwang, S.-F.; Maniolii, V.P.; Turnbull, G.; Freuda-Agyeman, R.; Hollman, K.; Kaus, S. Characterisation of clubroot (Plasmodiophora brassicae) from canola (Brassica napus) in the Peace Country of Alberta, Canada. Can. J. Plant Pathol. 2020, 1–7. [Google Scholar] [CrossRef]
  10. Chapara, V.; Kalwar, N.; Lubenow, L.; Chirumamilla, A. Prevalence of Clubroot on Canola in North Dakota. J. Agron. Agric. Sci. 2019, 2, 8. [Google Scholar] [CrossRef]
  11. Rastas, M.; Latvala, S.; Hannukkala, A. Occurrence of Plasmodiophora brassicae in Finnish turnip rape and oilseed rape fields. Agric. Food Sci. 2012, 21, 141–158. [Google Scholar] [CrossRef] [Green Version]
  12. Botero, A.; García, C.; Gossen, B.D.; Strelkov, S.E.; Todd, C.D.; Bonham-Smith, P.C.; Pérez-López, E. Clubroot disease in Latin America: Distribution and management strategies. Plant Pathol. 2019, 68, 827–833. [Google Scholar] [CrossRef] [Green Version]
  13. McGrann, G.M.D.; Gladders, P.; Smith, J.A.; Burnett, F. Control of clubroot (Plasmodiophora brassicae) in oilseed rape using varietal resistance and soil amendments. Field Crop. Res. 2016, 186, 146–156. [Google Scholar] [CrossRef]
  14. Diedichsen, E.; Frauen, M.; Ludwig-Muller, J. Clubroot management challenges from a German perspective. Can. J. Plant Pathol. 2014, 36, 85–98. [Google Scholar] [CrossRef]
  15. Jedryczka, M.; Kasprzyk, I.; Korbas, M.; Jajor, E.; Kaczmarek, J. Infestation of Polish agricultural soils by Plasmodiophora brassicae along the Polish—Ukrainian border. J. Plant Prot. Res. 2016, 54, 23–241. [Google Scholar] [CrossRef] [Green Version]
  16. Czubbatka-Biénkowska, A.; Kaczmarek, J.; Marzec-Schmidt, M.; Nieróbca, A.; Czajka, A.; Jedryczka, M. Country- wide qPCR based assessment of Plasmodiophora brassicae spread in agricultural soils and recommendations for the cultivation of Brassicaceae crops in Poland. Pathogens 2020, 9, 1070. [Google Scholar] [CrossRef]
  17. Řičařová, V.; Kaczmarek, J.; Strelkov, S.E.; Kazda, J.; Lueders, W.; Rysanek, P.; Manolii, V.; Jedryczka, M. Pathotypes of Plasmodiophora brassicae causing damage to oilseed rape in the Czech Republic and Poland. Eur. J. Plant Pathol. 2016, 145, 559–572. [Google Scholar] [CrossRef] [Green Version]
  18. Chai, A.L.; Xie, X.W.; Shi, Y.X.; Li, B.J. Research status of clubroot (Plasmodiophora brassicae) on cruciferous crops in China. Can. J. Plant Pathol. 2014, 36, 1421–1453. [Google Scholar] [CrossRef]
  19. Wallenhammar, A.-C.; Hansson, F.; Gunnarsson, A.; Jonsson, A. Quantification of Plasmodiophora brassicae Using a DNA-Based Soil Test Facilitates Sustainable Oilseed Rape Production. Plants 2016, 5, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kihlstrand, A. Biologisk markkartering håller koll på klumprotsjuka. Svensk Frötidning 2016, 2, 24. Available online: https://www.svenskraps.se/kunskap/pdf/01854.pdf (accessed on 15 January 2021). (In Swedish).
  21. Wallenhammar, A.-C. Prevalence of Plasmodiophora brassicae in a spring oilseed rape growing area in central Sweden and factors influencing soil infestation levels. Plant Pathol. 1996, 45, 710–719. [Google Scholar] [CrossRef]
  22. Crute, I.R.; Phelps, K.; Barnes, A.; Buczacki, S.T.; Crisp, P. The relationship between genotypes of three Brassica species and collections of Plasmodiophora brassicae Woron. Plant Pathol. 1983, 32, 405–420. [Google Scholar] [CrossRef]
  23. Jönsson, R. Klumprotsjuka. Litteraturöversikt och resultat av resistensarbeten i oljeväxter vid Svalöf 1962–1968. Sveriges Utsädesförenings Tidskr. 1969, 79, 28–56. (In Swedish) [Google Scholar]
  24. Strelkov, S.E.; Hwang, S.-F.; Manolii, V.P.; Cao, T.; Feindel, D. Emergence of new virulence phenotypes pf Plasmodiophora brassicae on canola (Brassica napus). Eur. J. Plant Pathol. 2016, 145, 517–529. [Google Scholar] [CrossRef]
  25. Zamani-Noor, N. Variation in pathotypes and virulence of Plasmodiophora brassicae populations in Germany. Plant Pathol. 2017, 66, 316–324. [Google Scholar] [CrossRef]
  26. Cao, T.; Manolii, V.P.; Hwang, S.F.; Howard, R.J.; Strelkov, S.E. Virulence and spread of Plasmodiophora brassicae [clubroot] in Alberta, Canada. Can. J. Plant Pathol. 2009, 31, 321–329. [Google Scholar] [CrossRef]
  27. Rennie, D.C.; Holtz, M.D.; Turkington, T.K.; Leboldous, J.M.; Hwang, S.-F.; Howard, R.J.; Strelkov, S.E. Movement of Plasmodiophora brassicae resting spores in windblown dust. Can. J. Plant Pathol. 2015, 37, 188–196. [Google Scholar] [CrossRef]
  28. Faggian, R.; Bulman, S.R.; Lawrie, A.C.; Porter, I.J. Specific polymerase chain reaction primers for the detection of Plasmodiophora brassicae in soil and water. Phytopathology 1999, 89, 392–397. [Google Scholar] [CrossRef] [Green Version]
  29. Kageyama, K.; Asano, T. Life cycle of Plasmodiophora brassicae. J. Plant Growth Regul. 2009, 28, 203–211. [Google Scholar] [CrossRef]
  30. Ludvig-Müller, J. Glucosinolates and the clubroot disease: Defense compounds or auxin precursors? Phytochem. Rev. 2009, 8, 135–148. [Google Scholar] [CrossRef]
  31. Wallenhammar, A.-C.; Almquist, C.; Söderström, M.; Jonsson, A. In-field distribution of Plasmodiophoraa brasssice measured using quantitative real-time PCR. Plant Pathol. 2012, 61, 16–28. [Google Scholar] [CrossRef]
  32. Jonsson, A.; Wallenhammar, A.-C.; Dixelius, C. Long-term soil data sets reveals shifts in population of Plasmodiophora brassica. In Proceedings of the 6th International Symposium on Brassica and 18th Crucifer Genetics Workshop, Catania, Italy, 12–16 November 2012. [Google Scholar]
  33. Zamani-Noor, N.; Krohne, I.; Koopmann, B. Greenhouse Evaluation of Clubroot Resistant-Brassica napus cv. Mendel and Its Efficacy Concerning Virulence and Soil Inoculum Levels of Plasmodiophora brassicae. Pathogens 2021, 10, 151. [Google Scholar] [CrossRef] [PubMed]
  34. Wallenhammar, A.-C.; Johnsson, L.; Gerhardson, B. Agronomic performance of partly clubroot-resistant spring oilseed turnip rape. J. Phytopathol. 2000, 148, 495–499. [Google Scholar] [CrossRef]
  35. Gunnarsson, A. Höstrapssorter med klumprotresistens. Field Trial Report of central Sweden. Rural Economy and Agricultural Society Skaraborg: Falköping 2015, 155–156. (In Swedish) [Google Scholar]
  36. Diedichsen, E.; Wagenblatt, B.; Schallehn, V. Production of pure genotype isolates of Plasmodiophora brassicae Wor.—Comparison of inoculations with root hairs containing single sporangiosori or with single resting spores. Eur. J. Plant Pathol. 2016, 145, 621–627. [Google Scholar] [CrossRef]
  37. Hwang, S.-F.; Strelkov, S.E.; Feng, J.; Gossen, B.D.; Howard, R.J. Plasmodiophora brassicae: A review of an emerging a pathogen of the Canadian canola (Brassica napus) crop. Mol. Plant Pathol. 2012, 13, 105–113. [Google Scholar] [CrossRef]
  38. Peng, G.; Lahlali, R.; Hwang, S.-F.; Pageau, D.; Hynes, R.K.; McDonald, M.R.; Gossen, B.D.; Strelkov, S.E. Crop rotation, cultivar resistance, and fungicides/ biofungicides for managing clubroot (Plasmodiophora brassice) on canola. Can. J. Plant Pathol. 2014, 36, 99–112. [Google Scholar] [CrossRef]
  39. Wallenhammar, A.-C. Observations from yield loss from Plasmodiophora brassicae infections in spring oilseed rape. J. Plant Dis. Prot. 1998, 105, 1–7. [Google Scholar]
  40. Liu, L.; Qin, L.; Cheng, X.; Zhang, Y.; Xu, L.; Liu, F.; Tong, C.; Huang, J.; Liu, S.; Wei, Y. Comparing the Infection Biology of Plasmodiophora brassicae in Clubroot Susceptible and Resistant Hosts and Non-hosts. Front Microbiol. 2020, 11, 507036. [Google Scholar] [CrossRef] [PubMed]
  41. Ernst, T.W.; Kher, S.; Stanton, D.; Rennie, D.C.; Hwang, S.F.; Strelkow, S.E. Plasmodiophora brassicae resting spore dynamics in clubroot resistant canola (Brassica napus) cropping systems. Plant Pathol. 2018, 68, 399–408. [Google Scholar] [CrossRef]
  42. Dobson, R.; Gabrielson, R.L.; Baker, A.S. Soil Water Matric Potential Requirements for Root-Hair and Cortical infection of Chinese Cabbage by Plasmodiophora brassicae. Phytopathology 1982, 12, 1598–1600. [Google Scholar] [CrossRef]
  43. Wellman, F.L. Clubroot of crucifers. Tech. Bull. U. S. Dept. Agric. 1930, 181, 1–31. [Google Scholar]
  44. Monteith, J. Relation of soil temperature and soil moisture to infection by Plasmodiophora brassicae. J. Agric. Res. 1924, 28, 549–561. [Google Scholar]
  45. Gossen, B.D.; Adhikari, K.K.C.; McDonald, M.R. Effects on temperature of infection and subsequent development of clubroot under controlled conditions. Plant Pathol. 2012, 61, 593–599. [Google Scholar] [CrossRef]
  46. Zhar, K.; Sarkes, A.; Yang, Y.; Zhou, Q.; Feindel, D.; Harding, M.W. Plasmodiophora brassicae in its environment- effects of temperature and light on resting spore survival in soil. Phytopathology 2021. [CrossRef]
  47. Dixon, G.R. Plasmodiophora brassicae in its environment. J. Plant Growth Regul. 2009, 28, 212–228. [Google Scholar] [CrossRef]
  48. Mattsson, L. Skördevariationer Inom ENSKILDA fält. Storlek och Tänkbara Orsaker; Report 196; Department of Soil Sciences, Division of Soil Fertility, Swedish University of Agricultural Sciences: Uppsala, Sweden, 1995. (In Swedish) [Google Scholar]
  49. Strelkov, S.E.; Hwang, S.-F.; Maniolii, V.P.; Cao, T.; Freuda-Agyeman, R.; Harding, M.W.; Peng, G.; Gossen, B.D.; McDonald, M.R.; Fiendel, D. Virulence and pathotype classification of Plasmodiophora brassicae populations collected from clubroot resistant canola (Brassica napus) in Canada. Can. J. Plant Pathol. 2018, 40, 284–298. [Google Scholar] [CrossRef]
Figure 1. Relationship between gene copies g−1 soil (determined by real-time qPCR) and harvested yield (kg, 9% water content) at the four trial sites (Tomelilla, Simrishamn, Kumla and Hallsberg) in 2018 of: (a) the susceptible Cultivar mix and (b) the clubroot-resistant cultivars cv. Mentor (blue), cv. Alister (yellow) and cv. Archimedes (red). Note: log scale on x-axis.
Figure 1. Relationship between gene copies g−1 soil (determined by real-time qPCR) and harvested yield (kg, 9% water content) at the four trial sites (Tomelilla, Simrishamn, Kumla and Hallsberg) in 2018 of: (a) the susceptible Cultivar mix and (b) the clubroot-resistant cultivars cv. Mentor (blue), cv. Alister (yellow) and cv. Archimedes (red). Note: log scale on x-axis.
Pathogens 10 00433 g001
Figure 2. Relationship between gene copies g−1 soil (determined by real-time qPCR) and disease severity index in field trials (DSIf) for (a) the susceptible Cultivar mix and (b) the clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes, at all four field trial sites in 2017 and at Kumla in 2018. Note: log scale on x-axis.
Figure 2. Relationship between gene copies g−1 soil (determined by real-time qPCR) and disease severity index in field trials (DSIf) for (a) the susceptible Cultivar mix and (b) the clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes, at all four field trial sites in 2017 and at Kumla in 2018. Note: log scale on x-axis.
Pathogens 10 00433 g002
Figure 3. Relationship between disease severity index and seed yield for (a) the susceptible Cultivar mix and (b) resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes sown at four different field sites in 2017 and at Kumla in 2018.
Figure 3. Relationship between disease severity index and seed yield for (a) the susceptible Cultivar mix and (b) resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes sown at four different field sites in 2017 and at Kumla in 2018.
Pathogens 10 00433 g003
Figure 4. Relationship between number of gene copies g−1 soil (determined by real-time qPCR) and disease severity index in bioassay (DSIb) in soil from plots with (a) susceptible Cultivar mix and clubroot-resistant cultivars (b) cv. Mentor, (c) cv. Alister and (d) cv. Archimedes grown in 2017. Note: log scale on x-axis.
Figure 4. Relationship between number of gene copies g−1 soil (determined by real-time qPCR) and disease severity index in bioassay (DSIb) in soil from plots with (a) susceptible Cultivar mix and clubroot-resistant cultivars (b) cv. Mentor, (c) cv. Alister and (d) cv. Archimedes grown in 2017. Note: log scale on x-axis.
Pathogens 10 00433 g004
Table 1. Disease severity index (DSIf) and disease incidence (DIf, percentage of diseased plants) of clubroot visually assessed on 15 December 2017 on roots originating from each plot in field trials harvested in 2018 with clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes and susceptible Cultivar mix. Sites were selected based on field abundance (gene copies g−1 soil) of Plasmodiophora brassicae DNA (Simrishamn 1,100,000; Tomelilla 2,500,000; Kumla 15,000; Hallsberg 2500).
Table 1. Disease severity index (DSIf) and disease incidence (DIf, percentage of diseased plants) of clubroot visually assessed on 15 December 2017 on roots originating from each plot in field trials harvested in 2018 with clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes and susceptible Cultivar mix. Sites were selected based on field abundance (gene copies g−1 soil) of Plasmodiophora brassicae DNA (Simrishamn 1,100,000; Tomelilla 2,500,000; Kumla 15,000; Hallsberg 2500).
Treatment and LocationDSIfDIf (%)
Simrishamn
Cultivar mix52.1b *67.9b
cv. Mentor12.0c15.6cd
cv. Alister8.5c10.8cd
cv. Archimedes5.3c6.0cd
Tomelilla
Cultivar mix55.4b62.5b
cv. Mentor8.3c11.2cd
cv. Alister8.0c9.2cd
cv. Archimedes2.4c4.0cd
Kumla
Cultivar mix92.3a100.0a
cv. Mentor15.5c15.1cd
cv. Alister19.5c16.9c
cv. Archimedes8.3c6.4cd
Hallsberg
Cultivar mix1.8c2.7cd
cv. Mentor0c0d
cv. Alister0c0d
cv. Archimedes0.7c2.0cd
p-value<0.001<0.001
Coefficient of variation56.147.4
* Different letters within each column indicate significant differences according to Tukey’s HSD-test (p < 0.05).
Table 2. Disease severity index (DSIf) and disease incidence (DIf, percentage of diseased plants) of clubroot visually assessed in early December 2018 and July 2019 on winter oilseed rape roots from each plot in field trials 2018/2019 at Simrishamn, Tomelilla, and Kumla for the clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes and the susceptible Cultivar mix. Sites were selected based on field abundance (gene copies g−1 soil) of Plasmodiophora brassicae DNA (Simrishamn 600,000; Tomelilla 370,000; Kumla 50,000).
Table 2. Disease severity index (DSIf) and disease incidence (DIf, percentage of diseased plants) of clubroot visually assessed in early December 2018 and July 2019 on winter oilseed rape roots from each plot in field trials 2018/2019 at Simrishamn, Tomelilla, and Kumla for the clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes and the susceptible Cultivar mix. Sites were selected based on field abundance (gene copies g−1 soil) of Plasmodiophora brassicae DNA (Simrishamn 600,000; Tomelilla 370,000; Kumla 50,000).
DSIf (0–100)DIf (%)
CultivarAutumn 2018July
2019
Both
years
Autumn 2018July
2019
Both
years
Simrishamn
Cultivar mix3.30.01.75.00.02.5
cv. Mentor0.00.70.30.02.01.0
cv. Alister0.00.00.00.00.00.0
cv. Archimedes0.00.00.00.00.00.0
p-valuensnsnsnsnsns
Coeff. of var. 400231475400231418
Tomelilla
Cultivar mix0.00.00.00.00.00.0
cv. Mentor0.01.00.50.01.01.5
cv. Alister1.90.00.95.60.02.8
cv. Archimedes0.00.00.00.00.00.0
p-valuensnsnsnsnsns
Coeff. of var.231400311231400311
Kumla
Cultivar mix22.527.5a *25.0a25.043.2a34.1a
cv. Mentor3.36.9ab5.9b7.510.9ab9.2b
cv. Alister5.81.7b5.1b7.58.0ab4.8b
cv. Archimedes6.75.3b3.8b6.72.0b7.4b
p-valuens0.0160.001ns0.0250.004
Coeff. of var.12899106131109115
Coefficient of variation. * Different letters indicate significant differences according to Tukey’s HSD-test (p < 0.05).
Table 3. Disease severity index (DSIb) and disease incidence (DIb, percentage of diseased plants) of clubroot in bioassays performed in a growth chamber with soil from field trials 2017/2018 with the susceptible Cultivar mix and clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes, and number of gene copies (logarithm at calculation) of Plasmodiophora brassicae. Values are mean for soil samples from each plot (n = 4). Sites were selected based on field abundance (gene copies g−1 soil) of Plasmodiophora brassicae DNA (Simrishamn 1,100,000; Tomelilla 2,500,000; Kumla 15,000; Hallsberg 2500).
Table 3. Disease severity index (DSIb) and disease incidence (DIb, percentage of diseased plants) of clubroot in bioassays performed in a growth chamber with soil from field trials 2017/2018 with the susceptible Cultivar mix and clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes, and number of gene copies (logarithm at calculation) of Plasmodiophora brassicae. Values are mean for soil samples from each plot (n = 4). Sites were selected based on field abundance (gene copies g−1 soil) of Plasmodiophora brassicae DNA (Simrishamn 1,100,000; Tomelilla 2,500,000; Kumla 15,000; Hallsberg 2500).
TreatmentDSIb
(0–100)
DSIb
(%)
Gene Copies
(No. g soil−1)
DSIb × Gene
Copies
DIb × Gene
Copies
Gene copies P. brassicae, average:
Simrishamn 419,200b *
Tomelilla 1,187,100a
Kumla 55,600c
Hallsberg 2800c
p-value <0.001
Coeff. of var. 119
Field-wise:
Simrishamn
Cultivar mix77.5a *94.5a562,300 p = 0.035ns
cv. Mentor21.9bcd30.5bcd435,700 p < 0.001p < 0.001
cv. Alister38.6b49.2b204,100 p = 0.012p = 0.027
cv. Archimedes34.4b46.6b474,700 nsns
p-value0.0200.020ns
Coeff. of var.141134103
Tomelilla
Cultivar mix75.2a86.4a784,300 nsns
cv. Mentor31.6b40.2bc1,129,800 nsns
cv. Alister32.8b42.9b1,603,200 nsns
cv. Archimedes7.0de16.5cde1,231,000 nsns
p-value<0.001nsns
Coeff. of var.716672
Kumla
Cultivar mix29.9bc42.0bc68,300 nsns
cv. Mentor9.0cde14.8de49,500 p = 0.011p = 0.009
cv. Alister7.5cde9.5de38,800 nsns
cv. Archimedes4.4de8.2de65,800 p = 0.017p = 0.009
p-valuensnsns
Coeff. of var.32332376
Hallsberg
Cultivar mix1.3de2.8e3800 nsns
cv. Mentor0.5e1.2e1700 nsns
cv. Alister0.8e1.8e2800 nsns
cv. Archimedes0.0e0.0e3020 nsns
p-valuensnsns
Coeff. of var.32332377
Coefficient of variation. * Different letters indicate significant differences according to Tukey’s HSD-test (p < 0.05).
Table 4. Disease severity index (DSIb) and disease incidence (DIb, percentage of diseased plants) of clubroot in bioassays performed in a growth chamber with soil from the field trials 2018/2019, linked to the number of gene copies (logarithm at calculation) of Plasmodiophora brassicae for susceptible Cultivar mix and clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes. Values are mean for all samples from each plot (n = 4). Sites were selected based on field abundance (gene copies g−1 soil) of P. brassicae DNA (Simrishamn 600,000; Tomelilla 370,000; Kumla 50,000).
Table 4. Disease severity index (DSIb) and disease incidence (DIb, percentage of diseased plants) of clubroot in bioassays performed in a growth chamber with soil from the field trials 2018/2019, linked to the number of gene copies (logarithm at calculation) of Plasmodiophora brassicae for susceptible Cultivar mix and clubroot-resistant cultivars cv. Mentor, cv. Alister, and cv. Archimedes. Values are mean for all samples from each plot (n = 4). Sites were selected based on field abundance (gene copies g−1 soil) of P. brassicae DNA (Simrishamn 600,000; Tomelilla 370,000; Kumla 50,000).
TreatmentDSIb
(0–100)
DIb
(%)
Gene Copies
(No. g soil−1)
DSIb × Gene
Copies
DIb × Gene
Copies
Gene copies P. brassicae, average
Simrishamn 442,300a *
Tomelilla 38,000b
Kumla 3700b
P-value <0.001
Coeff. of var. 113
Field wise
Simrishamn
Cultivar mix15.8a19.4a633,900ansns
cv. Mentor4.6a5.6b471,500abnsns
cv. Alister16.2a19.2a299,900bnsns
cv. Archimedes4.0a5.8b364,300abnsns
p-value0.0200.0200.015
Coeff. of var.1411344.0
Tomelilla
Cultivar mix54.6a61.8 48,400ansns
cv. Mentor22.9b25.7 42,400abnsns
cv. Alister29.8b33.7 27,300cnsns
cv. Archimedes14.6b16.9 35,200bcnsns
p-value<0.001ns<0.001
Coeff. of var.71.266.22.9
Kumla
Cultivar mix2.5 2.5a3400 nsns
cv. Mentor0.7 0.7ab3500 nsns
cv. Alister0.0 0.0b4000 nsns
cv. Archimedes1.7 1.7b3600 nsns
p-valuensnsns
Coeff. of var.32332375
Coefficient of variation. * Different letters indicate significant differences according to Tukey’s HSD-test (p < 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wallenhammar, A.-C.; Omer, Z.S.; Edin, E.; Jonsson, A. Influence of Soil-Borne Inoculum of Plasmodiophora brassicae Measured by qPCR on Disease Severity of Clubroot-Resistant Cultivars of Winter Oilseed Rape (Brassica napus L.). Pathogens 2021, 10, 433. https://doi.org/10.3390/pathogens10040433

AMA Style

Wallenhammar A-C, Omer ZS, Edin E, Jonsson A. Influence of Soil-Borne Inoculum of Plasmodiophora brassicae Measured by qPCR on Disease Severity of Clubroot-Resistant Cultivars of Winter Oilseed Rape (Brassica napus L.). Pathogens. 2021; 10(4):433. https://doi.org/10.3390/pathogens10040433

Chicago/Turabian Style

Wallenhammar, Ann-Charlotte, Zahra Saad Omer, Eva Edin, and Anders Jonsson. 2021. "Influence of Soil-Borne Inoculum of Plasmodiophora brassicae Measured by qPCR on Disease Severity of Clubroot-Resistant Cultivars of Winter Oilseed Rape (Brassica napus L.)" Pathogens 10, no. 4: 433. https://doi.org/10.3390/pathogens10040433

APA Style

Wallenhammar, A. -C., Omer, Z. S., Edin, E., & Jonsson, A. (2021). Influence of Soil-Borne Inoculum of Plasmodiophora brassicae Measured by qPCR on Disease Severity of Clubroot-Resistant Cultivars of Winter Oilseed Rape (Brassica napus L.). Pathogens, 10(4), 433. https://doi.org/10.3390/pathogens10040433

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop