Genetics of Clubroot Resistance in Brassica Species

Optimizing Clubroot Management and the Role of Canola Cultivar Mixtures

The sustainable cultivation of canola is under threat from clubroot disease (Plasmodiophora brassicae). The pathogen's resting spores can survive in the soil for extended periods, complicating disease management. Therefore, effective clubroot control requires a combination of tactics that provide multiple layers of protection. Management strategies have focused on pathogen avoidance and reducing disease levels in infested fields. The sanitation of machinery and field equipment remains the most effective method for preventing the pathogen's introduction into non-infested fields. For disease reduction, crop rotation, liming, chemical control, and host resistance are commonly employed, with the use of clubroot-resistant cultivars being the most effective to date. However, resistance breakdown has been observed within four years of the introduction of new cultivars, jeopardizing the long-term effectiveness of this approach. A promising yet underexplored strategy is the use of cultivar mixtures. This approach leverages mechanisms such as the dilution effect, the barrier effect, induced resistance, disruptive selection, and the compensatory effect to control the disease. Cultivar mixtures have the potential to reduce the impact of clubroot on canola production while preserving pathogen population structure, thereby minimizing the likelihood of resistance breakdown. Given its potential, further research into cultivar mixtures as a management strategy for clubroot disease is warranted.

Pyramiding of triple Clubroot resistance loci conferred superior resistance without negative effects on agronomic traits in Brassica napus

Clubroot disease caused by Plasmodiophora brassicae is becoming a serious threat to rapeseed (Brassica napus) production worldwide. Breeding resistant varieties using CR (clubroot resistance) loci is the most promising solution. Using marker-assisted selection and speed-breeding technologies, we generated Brassica napus materials in homozygous or heterozygous states using CRA3.7, CRA08.1, and CRA3.2 loci in the elite parental line of the Zhongshuang11 background. We developed three elite lines with two CR loci in different combinations and one line with three CR loci at the homozygous state. In our study, we used six different clubroot strains (Xinmin, Lincang, Yuxi, Chengdu, Chongqing, and Jixi) which are categorized into three groups based on our screening results. The newly pyramided lines with two or more CR loci displayed better disease resistance than the parental lines carrying single CR loci. There is an obvious gene dosage effect between CR loci and disease resistance levels. For example, pyramided lines with triple CR loci in the homozygous state showed superior resistance for all pathogens tested. Moreover, CR loci in the homozygous state are better on disease resistance than the heterozygous state. More importantly, no negative effect was observed on agronomic traits for the presence of multiple CR loci in the same background. Overall, these data suggest that the pyramiding of triple clubroot resistance loci conferred superior resistance with no negative effects on agronomic traits in Brassica napus.

Resynthesizing Brassica napus with race specific resistance genes and race non-specific QTLs to multiple races of Plasmodiophora brassicae

Clubroot disease in canola (Brassica napus) continues to spread across the Canadian prairies. Growing resistant cultivars is considered the most economical means of controlling the disease. However, sources of resistance to clubroot in B. napus are very limited. In this study, we conducted interspecific crosses using a B. rapa line (T19) carrying race-specific resistance genes and two B. oleracea lines, ECD11 and JL04, carrying race non-specific QTLs. Employing embryo rescue and conventional breeding methods, we successfully resynthesized a total of eight B. napus lines, with four derived from T19 × ECD11 and four from T19 × JL04. Additionally, four semi-resynthesized lines were developed through crosses with a canola line (DH16516). Testing for resistance to eight significant races of Plasmodiophora brassicae was conducted on seven resynthesized lines and four semi-resynthesized lines. All lines exhibited high resistance to the strains. Confirmation of the presence of clubroot resistance genes/QTLs was performed in the resynthesized lines using SNP markers linked to race-specific genes in T19 and race non-specific QTLs in ECD11. The developed B. napus germplasms containing clubroot resistance are highly valuable for the development of canola cultivars resistant to clubroot.

Advances in and Perspectives on Transgenic Technology and CRISPR-Cas9 Gene Editing in Broccoli

Broccoli, a popular international Brassica oleracea crop, is an important export vegetable in China. Broccoli is not only rich in protein, vitamins, and minerals but also has anticancer and antiviral activities. Recently, an Agrobacterium-mediated transformation system has been established and optimized in broccoli, and transgenic transformation and CRISPR-Cas9 gene editing techniques have been applied to improve broccoli quality, postharvest shelf life, glucoraphanin accumulation, and disease and stress resistance, among other factors. The construction and application of genetic transformation technology systems have led to rapid development in broccoli worldwide, which is also good for functional gene identification of some potential traits in broccoli. This review comprehensively summarizes the progress in transgenic technology and CRISPR-Cas9 gene editing for broccoli over the past four decades. Moreover, it explores the potential for future integration of digital and smart technologies into genetic transformation processes, thus demonstrating the promise of even more sophisticated and targeted crop improvements. As the field continues to evolve, these innovations are expected to play a pivotal role in the sustainable production of broccoli and the enhancement of its nutritional and health benefits.

Fine mapping and candidate gene analysis of CRA8.1.6, which confers clubroot resistance in turnip (Brassica rapa ssp. rapa)

Clubroot disease poses a significant threat to Brassica crops, necessitating ongoing updates on resistance gene sources. In F2 segregants of the clubroot-resistant inbred line BrT18-6-4-3 and susceptible DH line Y510, the genetic analysis identified a single dominant gene responsible for clubroot resistance. Through bulk segregant sequencing analysis and kompetitive allele-specific polymerase chain reaction assays, CRA8.1.6 was mapped within 110 kb (12,255-12,365 Mb) between markers L-CR11 and L-CR12 on chromosome A08. We identified B raA08g015220.3.5C as the candidate gene of CRA8.1.6. Upon comparison with the sequence of disease-resistant material BrT18-6-4-3, we found 249 single-nucleotide polymorphisms, seven insertions, six deletions, and a long terminal repeat (LTR) retrotransposon (5,310 bp) at 909 bp of the first intron. However, the LTR retrotransposon was absent in the coding sequence of the susceptible DH line Y510. Given the presence of a non-functional LTR insertion in other materials, it showed that the LTR insertion might not be associated with susceptibility. Sequence alignment analysis revealed that the fourth exon of the susceptible line harbored two deletions and an insertion, resulting in a frameshift mutation at 8,551 bp, leading to translation termination at the leucine-rich repeat domain's C-terminal in susceptible material. Sequence alignment of the CDS revealed a 99.4% similarity to Crr1a, which indicate that CRA8.1.6 is likely an allele of the Crr1a gene. Two functional markers, CRA08-InDel and CRA08-KASP1, have been developed for marker-assisted selection in CR turnip cultivars. Our findings could facilitate the development of clubroot-resistance turnip cultivars through marker-assisted selection.

Proteome- and metabolome-level changes during early stages of clubroot infection in Brassica napus canola

Clubroot is a destructive root disease of canola (Brassica napus L.) caused by Plasmodiophora brassicae Woronin. Despite extensive research into the molecular responses of B. napus to P. brassicae, there is limited information on proteome- and metabolome-level changes in response to the pathogen, especially during the initial stages of infection. In this study, we have investigated the proteome- and metabolome- level changes in the roots of clubroot-resistant (CR) and -susceptible (CS) doubled-haploid (DH) B. napus lines, in response to P. brassicae pathotype 3H at 1-, 4-, and 7-days post-inoculation (DPI). Root proteomes were analyzed using nanoflow liquid chromatography coupled with tandem mass spectrometry (nano LC-MS/MS). Comparisons of pathogen-inoculated and uninoculated root proteomes revealed 2515 and 1556 differentially abundant proteins at one or more time points (1-, 4-, and 7-DPI) in the CR and CS genotypes, respectively. Several proteins related to primary metabolites (e.g., amino acids, fatty acids, and lipids), secondary metabolites (e.g., glucosinolates), and cell wall reinforcement-related proteins [e.g., laccase, peroxidases, and plant invertase/pectin methylesterase inhibitors (PInv/PMEI)] were identified. Eleven nucleotides and nucleoside-related metabolites, and eight fatty acids and sphingolipid-related metabolites were identified in the metabolomics study. To our knowledge, this is the first report of root proteome-level changes and associated alterations in metabolites during the early stages of P. brassicae infection in B. napus.

RNA-Seq Bulked Segregant Analysis of an Exotic B. napus ssp. napobrassica (Rutabaga) F2 Population Reveals Novel QTLs for Breeding Clubroot-Resistant Canola

In this study, a rutabaga (Brassica napus ssp. napobrassica) donor parent FGRA106, which exhibited broad-spectrum resistance to 17 isolates representing 16 pathotypes of Plasmodiophora brassicae, was used in genetic crosses with the susceptible spring-type canola (B. napus ssp. napus) accession FG769. The F2 plants derived from a clubroot-resistant F1 plant were screened against three P. brassicae isolates representing pathotypes 3A, 3D, and 3H. Chi-square (χ2) goodness-of-fit tests indicated that the F2 plants inherited two major clubroot resistance genes from the CR donor FGRA106. The total RNA from plants resistant (R) and susceptible (S) to each pathotype were pooled and subjected to bulked segregant RNA-sequencing (BSR-Seq). The analysis of gene expression profiles identified 431, 67, and 98 differentially expressed genes (DEGs) between the R and S bulks. The variant calling method indicated a total of 12 (7 major + 5 minor) QTLs across seven chromosomes. The seven major QTLs included: BnaA5P3A.CRX1.1, BnaC1P3H.CRX1.2, and BnaC7P3A.CRX1.1 on chromosomes A05, C01, and C07, respectively; and BnaA8P3D.CRX1.1, BnaA8P3D.RCr91.2/BnaA8P3H.RCr91.2, BnaA8P3H.Crr11.3/BnaA8P3D.Crr11.3, and BnaA8P3D.qBrCR381.4 on chromosome A08. A total of 16 of the DEGs were located in the major QTL regions, 13 of which were on chromosome C07. The molecular data suggested that clubroot resistance in FGRA106 may be controlled by major and minor genes on both the A and C genomes, which are deployed in different combinations to confer resistance to the different isolates. This study provides valuable germplasm for the breeding of clubroot-resistant B. napus cultivars in Western Canada.

Effects of Exogenous Ergothioneine on Brassica rapa Clubroot Development Revealed by Transcriptomic Analysis

Clubroot disease is a soil-borne disease caused by Plasmodiophora brassicae that leads to a serious yield reduction in cruciferous plants. In this study, ergothioneine (EGT) was used to culture P. brassicae resting spores, the germination of which was significantly inhibited. Further exogenous application of EGT and P. brassicae inoculation in Chinese cabbage showed that EGT promoted root growth and significantly reduced the incidence rate and disease index. To further explore the mechanism by which EGT improves the resistance of Chinese cabbage to clubroot, a Chinese cabbage inbred line BJN3-2 susceptible to clubroot treated with EGT was inoculated, and a transcriptome analysis was conducted. The transcriptome sequencing analysis showed that the differentially expressed genes induced by EGT were significantly enriched in the phenylpropanoid biosynthetic pathway, and the genes encoding related enzymes involved in lignin synthesis were upregulated. qRT-PCR, peroxidase activity, lignin and flavonoid content determination showed that EGT promoted the lignin and flavonoid synthesis of Chinese cabbage and improved its resistance to clubroot. This study provides a new insight for the comprehensive prevention and control of cruciferous clubroot and for further study of the effects of EGT on clubroot disease.

Mapping of a novel clubroot disease resistance locus in Brassica napus and related functional identification

Clubroot disease, caused by Plasmodiophora brassicae, is a devastating disease that results in substantial yield loss in Brassicaceae crops worldwide. In this study, we identified a clubroot disease resistance (CR) Brassica napus, "Kc84R," which was obtained by mutation breeding. Genetic analysis revealed that the CR trait of "Kc84R" was controlled by a single dominant locus. We used the bulked segregant analysis sequencing (BSA-seq) approach, combined with genetic mapping based on single nucleotide polymorphism (SNP) markers to identify CR loci from the F₂ population derived from crossing CR "Kc84R" and clubroot susceptible "855S." The CR locus was mapped to a region between markers BnSNP14198336 and BnSNP14462201 on the A03 chromosome, and this fragment of 267 kb contained 68 annotated candidate genes. Furthermore, we performed the CR relation screening of candidate genes with the model species Arabidopsis. An ERF family transcriptional activator, BnERF034, was identified to be associated with the CR, and the corresponding Arabidopsis homozygous knockout mutants exhibited more pronounced resistance compared with the wild-type Col-0 and the transgenic lines of BnERF034 in response to P. brassicae infection. Additionally, the expression analysis between resistant and susceptible materials indicated that BnERF034 was identified to be the most likely CR candidate for the resistance in Kc84R. To conclude, this study reveals a novel gene responsible for CR. Further analysis of BnERF034 may reveal the molecular mechanisms underlying the CR of plants and provide a theoretical basis for Brassicaceae resistance breeding.

Germplasm Enhancement and Identification of Loci Conferring Resistance against Plasmodiophora brassicae in Broccoli

In order to breed broccoli and other Brassica materials to be highly resistant to clubroot disease, 41 Brassicaceae varieties were developed and identified between 2020 and 2021. Seven known clubroot genes were used for screening these materials. In addition, the resistant and susceptible broccoli cultivars were designed for observing their differences in the infection process with Plasmodiophora brassicae. The results showed that 90% of total materials had carried more than two clubroot resistance genes: one material carried two disease resistance genes, four materials carried seven genes for clubroot resistance, two materials carried six genes for clubroot resistance, and in total 32% of these materials carried five genes for clubroot resistance. As a result, several new genotypes of Brassicaceae germplasm were firstly created and obtained based on distant hybridization and identification of loci conferring resistance against Plasmodiophora brassicae in this study. We found and revealed that similar infection models of Plasmodiophora brassicae occurred in susceptible and resistant cultivars of broccoli, but differences in infection efficiency of Plasmodiophora brassicae also existed in both materials. For resistant broccoli plants, a small number of conidia formed in the root hair, and only a few spores could enter the cortex without forming sporangia while sporangia could form in susceptible plants. Our study could provide critical Brassica materials for breeding resistant varieties and new insight into understanding the mechanism of plant resistance.

Multi-Omics Approaches to Improve Clubroot Resistance in Brassica with a Special Focus on Brassica oleracea L

Brassica oleracea is an agronomically important species of the Brassicaceae family, including several nutrient-rich vegetables grown and consumed across the continents. But its sustainability is heavily constrained by a range of destructive pathogens, among which, clubroot disease, caused by a biotrophic protist Plasmodiophora brassicae, has caused significant yield and economic losses worldwide, thereby threatening global food security. To counter the pathogen attack, it demands a better understanding of the complex phenomenon of Brassica-P. brassicae pathosystem at the physiological, biochemical, molecular, and cellular levels. In recent years, multiple omics technologies with high-throughput techniques have emerged as successful in elucidating the responses to biotic and abiotic stresses. In Brassica spp., omics technologies such as genomics, transcriptomics, ncRNAomics, proteomics, and metabolomics are well documented, allowing us to gain insights into the dynamic changes that transpired during host-pathogen interactions at a deeper level. So, it is critical that we must review the recent advances in omics approaches and discuss how the current knowledge in multi-omics technologies has been able to breed high-quality clubroot-resistant B. oleracea. This review highlights the recent advances made in utilizing various omics approaches to understand the host resistance mechanisms adopted by Brassica crops in response to the P. brassicae attack. Finally, we have discussed the bottlenecks and the way forward to overcome the persisting knowledge gaps in delivering solutions to breed clubroot-resistant Brassica crops in a holistic, targeted, and precise way.

Identification of Novel Locus RsCr6 Related to Clubroot Resistance in Radish (Raphanus sativus L.)

Clubroot is a devastating disease that causes substantial yield loss worldwide. However, the inheritance and molecular mechanisms of clubroot resistance during pathogen infection in radish remain largely unclear. In this study, we investigated the inheritance of clubroot resistance in the F₂ population derived from crossing clubroot-resistant (CR) and clubroot-susceptible inbred lines "GLX" and "XNQ," respectively. Genetic analysis revealed that a single dominant gene controlled the clubroot resistance of "GLX" with a Mendelian ratio of resistance and susceptibility of nearly 3:1. Bulked segregant analysis combined with whole-genome resequencing (BSA-seq) was performed to detect the target region of RsCr6 on chromosome Rs8. Linkage analysis revealed that the RsCr6 locus was located between two markers, HB321 and HB331, with an interval of approximately 92 kb. Based on the outcomes of transcriptome analysis, in the RsCr6 locus, the R120263140 and R120263070 genes with a possible relation to clubroot resistance were considered candidate genes. In addition, three core breeding materials containing the two reported quantitative trait loci (QTLs) and our novel locus RsCr6 targeting clubroot resistance were obtained using marker-assisted selection (MAS) technology. This study reveals a novel locus responsible for clubroot resistance in radishes. Further analysis of new genes may reveal the molecular mechanisms underlying the clubroot resistance of plants and provide a theoretical basis for radish resistance breeding.

Control Strategies of Clubroot Disease Caused by Plasmodiophora brassicae

The clubroot disease caused by the soil-borne pathogen Plasmodiophora brassicae is one of the most important diseases of cruciferous crops worldwide. As with many plant pathogens, the spread is closely related to the cultivation of suitable host plants. In addition, temperature and water availability are crucial determinants for the occurrence and reproduction of clubroot disease. Current global changes are contributing to the widespread incidence of clubroot disease. On the one hand, global trade and high prices are leading to an increase in the cultivation of the host plant rapeseed worldwide. On the other hand, climate change is improving the living conditions of the pathogen P. brassicae in temperate climates and leading to its increased occurrence. Well-known ways to control efficiently this disease include arable farming strategies: growing host plants in wide crop rotations, liming the contaminated soils, and using resistant host plants. Since chemical control of the clubroot disease is not possible or not ecologically compatible, more and more alternative control options are being investigated. In this review, we address the challenges for its control, with a focus on biological control options.

Identification of Two Major QTLs in Brassica napus Lines With Introgressed Clubroot Resistance From Turnip Cultivar ECD01

Plasmodiophora brassicae causes clubroot disease in brassica crops worldwide. Brassica rapa, a progenitor of Brassica napus (canola), possesses important sources for resistance to clubroot. A doubled haploid (DH) population consisting of 84 DH lines were developed from a Backcross2 (BC₂) plant through an interspecific cross of B. rapa turnip cv. ECD01 (resistant, R) with canola line DH16516 (susceptible, S) and then backcrossed with DH16516 as the recurrent parent. The DH lines and their parental lines were tested for resistance to four major pathotypes (3A, 3D, 3H, and 5X) of P. brassicae identified from canola. The R:S segregation ratio for pathotype 3A was 1:3, and 3:1 for pathotypes 3D, 3H, and 5X. From genotyping by sequencing (GBS), a total of 355.3 M short reads were obtained from the 84 DH lines, ranging from 0.81 to 11.67 M sequences per line. The short reads were aligned into the A-genome of B. napus "Darmor-bzh" version 4.1 with a total of 260 non-redundant single-nucleotide polymorphism (SNP) sites. Two quantitative trait loci (QTLs), Rcr10 ^ECD01 and Rcr9 ^ECD01 , were detected for the pathotypes in chromosomes A03 and A08, respectively. Rcr10 ^ECD01 and Rcr9 ^ECD01 were responsible for resistance to 3A, 3D, and 3H, while only one QTL, Rcr9 ^ECD01 , was responsible for resistance to pathotype 5X. The logarithm of the odds (LOD) values, phenotypic variation explained (PVE), additive (Add) values, and confidence interval (CI) from the estimated QTL position varied with QTL, with a range of 5.2-12.2 for LOD, 16.2-43.3% for PVE, 14.3-25.4 for Add, and 1.5-12.0 cM for CI. The presence of the QTLs on the chromosomes was confirmed through the identification of the percentage of polymorphic variants using bulked-segregant analysis. There was one gene encoding a disease resistance protein and 24 genes encoding proteins with function related to plant defense response in the Rcr10 ^ECD01 target region. In the Rcr9 ^ECD01 region, two genes encoded disease resistance proteins and 10 genes encoded with defense-related function. The target regions for Rcr10 ^ECD01 and Rcr9 ^ECD01 in B. napus were homologous to the 11.0-16.0 Mb interval of chromosome A03 and the 12.0-14.5 Mb interval of A08 in B. rapa "Chiifu" reference genome, respectively.

Multi-Omic Investigation of Low-Nitrogen Conditional Resistance to Clubroot Reveals Brassica napus Genes Involved in Nitrate Assimilation

Nitrogen fertilization has been reported to influence the development of clubroot, a root disease of Brassicaceae species, caused by the obligate protist Plasmodiophora brassicae. Our previous works highlighted that low-nitrogen fertilization induced a strong reduction of clubroot symptoms in some oilseed rape genotypes. To further understand the underlying mechanisms, the response to P. brassicae infection was investigated in two genotypes "Yudal" and HD018 harboring sharply contrasted nitrogen-driven modulation of resistance toward P. brassicae. Targeted hormone and metabolic profiling, as well as RNA-seq analysis, were performed in inoculated and non-inoculated roots at 14 and 27 days post-inoculation, under high and low-nitrogen conditions. Clubroot infection triggered a large increase of SA concentration and an induction of the SA gene markers expression whatever the genotype and nitrogen conditions. Overall, metabolic profiles suggested that N-driven induction of resistance was independent of SA signaling, soluble carbohydrate and amino acid concentrations. Low-nitrogen-driven resistance in "Yudal" was associated with the transcriptional regulation of a small set of genes, among which the induction of NRT2- and NR-encoding genes. Altogether, our results indicate a possible role of nitrate transporters and auxin signaling in the crosstalk between plant nutrition and partial resistance to pathogens.

Using wild relatives and related species to build climate resilience in Brassica crops

Climate change will have major impacts on crop production: not just increasing drought and heat stress, but also increasing insect and disease loads and the chance of extreme weather events and further adverse conditions. Often, wild relatives show increased tolerances to biotic and abiotic stresses, due to reduced stringency of selection for yield and yield-related traits under optimum conditions. One possible strategy to improve resilience in our modern-day crop cultivars is to utilize wild relative germplasm in breeding, and attempt to introgress genetic factors contributing to greater environmental tolerances from these wild relatives into elite crop types. However, this approach can be difficult, as it relies on factors such as ease of hybridization and genetic distance between the source and target, crossover frequencies and distributions in the hybrid, and ability to select for desirable introgressions while minimizing linkage drag. In this review, we outline the possible effects that climate change may have on crop production, introduce the Brassica crop species and their wild relatives, and provide an index of useful traits that are known to be present in each of these species that may be exploitable through interspecific hybridization-based approaches. Subsequently, we outline how introgression breeding works, what factors affect the success of this approach, and how this approach can be optimized so as to increase the chance of recovering the desired introgression lines. Our review provides a working guide to the use of wild relatives and related crop germplasm to improve biotic and abiotic resistances in Brassica crop species.

Using wild relatives and related species to build climate resilience in Brassica crops

Climate change will have major impacts on crop production: not just increasing drought and heat stress, but also increasing insect and disease loads and the chance of extreme weather events and further adverse conditions. Often, wild relatives show increased tolerances to biotic and abiotic stresses, due to reduced stringency of selection for yield and yield-related traits under optimum conditions. One possible strategy to improve resilience in our modern-day crop cultivars is to utilize wild relative germplasm in breeding, and attempt to introgress genetic factors contributing to greater environmental tolerances from these wild relatives into elite crop types. However, this approach can be difficult, as it relies on factors such as ease of hybridization and genetic distance between the source and target, crossover frequencies and distributions in the hybrid, and ability to select for desirable introgressions while minimizing linkage drag. In this review, we outline the possible effects that climate change may have on crop production, introduce the Brassica crop species and their wild relatives, and provide an index of useful traits that are known to be present in each of these species that may be exploitable through interspecific hybridization-based approaches. Subsequently, we outline how introgression breeding works, what factors affect the success of this approach, and how this approach can be optimized so as to increase the chance of recovering the desired introgression lines. Our review provides a working guide to the use of wild relatives and related crop germplasm to improve biotic and abiotic resistances in Brassica crop species.

Identification of resistance loci against new pathotypes of Plasmodiophora brassicae in Brassica napus based on genome-wide association mapping

Genetic resistance is a successful strategy for management of clubroot (Plasmodiophora brassicae) of brassica crops, but resistance can break down quickly. Identification of novel sources of resistance is especially important when new pathotypes arise. In the current study, the reaction of 177 accessions of Brassica napus to four new, virulent pathotypes of P. brassicae was assessed. Each accession was genotyped using genotyping by sequencing to identify and map novel sources of clubroot resistance using mixed linear model (MLM) analysis. The majority of accessions were highly susceptible (70–100 DSI), but a few accessions exhibited strong resistance (0–20 DSI) to pathotypes 5X (21 accessions), 3A (8), 2B (7), and 3D (15), based on the Canadian Clubroot Differential system. In total, 301,753 SNPs were mapped to 19 chromosomes. Population structure analysis indicated that the 177 accessions belong to seven major populations. SNPs were associated with resistance to each pathotype using MLM. In total, 13 important SNP loci were identified, with 9 SNPs mapped to the A-genome and 4 to the C-genome. The SNPs were associated with resistance to pathotypes 5X (2 SNPs), 3A (4), 2B (5) and 3D (6). A Blast search of 1.6 Mb upstream and downstream from each SNP identified 13 disease-resistance genes or domains. The distance between a SNP locus and the nearest resistance gene ranged from 0.04 to 0.74 Mb. The resistant lines and SNP markers identified in this study can be used to breed for resistance to the most prevalent new pathotypes of P. brassicae in Canada.

Clubroot in Brassica: recent advances in genomics, breeding, and disease management

Clubroot disease, caused by Plasmodiophora brassicae, affects Brassica oilseed and vegetable production worldwide. This review is focused on various aspects of clubroot disease and its management, including understanding the pathogen and resistance in the host plants. Advances in genetics, molecular biology techniques, and omics research have helped to identify several major loci, QTL, and genes from the Brassica genomes involved in the control of clubroot resistance. Transcriptomic studies have helped to extend our understanding of the mechanism of infection by the pathogen and the molecular basis of resistance/susceptibility in the host plants. A comprehensive understanding of the clubroot disease and host resistance would allow developing a better strategy by integrating the genetic resistance with cultural practices to manage this disease from a long-term perspective.

Virulence Spectrum of Single-Spore and Field Isolates of Plasmodiophora brassicae Able to Overcome Resistance in Canola (Brassica napus)

Clubroot, caused by Plasmodiophora brassicae Woronin, is an important disease of canola (Brassica napus L.) that is managed mainly by planting clubroot-resistant (CR) cultivars. Field isolates of P. brassicae can be heterogeneous mixtures of various pathotypes, making assessments of the genetics of host-pathogen interactions challenging. Thirty-four single-spore isolates were obtained from nine field isolates of the pathogen collected from CR canola cultivars. The virulence patterns of the single-spore and field isolates were assessed on the 13 host genotypes of the Canadian Clubroot Differential (CCD) set, which includes the differentials of Williams and Somé et al. Indices of disease (IDs) severity of 25, 33, and 50% (±95% confidence interval) were compared as potential thresholds to distinguish between resistant and susceptible reactions, with an ID of 50% giving the most consistent responses for pathotype classification purposes. With this threshold, 13 pathotypes could be distinguished based on the CCD system, 7 on the differentials of Williams, and 3 on the hosts of Somé et al. The highest correlations were observed among virulence matrices generated using the three threshold IDs on the CCD set. Genetically homogeneous single-spore isolates gave a clearer profile of the P. brassicae pathotype structure. Novel pathotypes, not reported in Canada previously, were identified among the isolates. This large collection of single-spore isolates can serve as a reference in screening and breeding for clubroot resistance.

Comparing the Infection Biology of Plasmodiophora brassicae in Clubroot Susceptible and Resistant Hosts and Non-hosts

The potential infection biology of Plasmodiophora brassicae in resistant hosts and non-hosts is still not completely understood. Clubroot resistance assay on European clubroot differentials (ECD) set revealed that ECD10 (Brassica napus) and ECD4 (Brassica rapa) show a complete resistance to the tested P. brassicae isolate in contrast to highly susceptible hosts Westar (B. napus) and ECD5 (B. rapa). Previously, we used fluorescent probe-based confocal microscopy (FCM) to refine the life cycle of P. brassicae and indicate the important time points during its infection in Arabidopsis. Here, we used FCM to systematically investigate the infection of P. brassicae in two resistant host species ECD10 and ECD4 and two non-host crops wheat and barley at each indicated time points, compared with two susceptible hosts Westar and ECD5. We found that P. brassicae can initiate the primary infection phase and produce uninucleate primary plasmodia in both resistant hosts and non-hosts just like susceptible hosts at 2 days post-inoculation (dpi). Importantly, P. brassicae can develop into zoosporangia and secondary zoospores and release the secondary zoospores from the zoosporangia in resistant hosts at 7 dpi, comparable to susceptible hosts. However, during the secondary infection phase, no secondary plasmodium was detected in the cortical cells of both resistant hosts in contrast to massive secondary plasmodia present in the cortex tissue of two susceptible hosts leading to root swelling at 15 dpi. In both non-host crops, only uninucleate primary plasmodia were observed throughout roots at 7 and 15 dpi. Quantitative PCR based on DNA revealed that the biomass of P. brassicae has no significant increase from 2 dpi in non-host plants and from 7 dpi in resistant host plants, compared to the huge biomass increase in susceptible host plants from 2 to 25 dpi. Our study reveals that the primary infection phase in the root epidermis and the secondary infection phase in the cortex tissue are, respectively, blocked in non-hosts and resistant hosts, contributing to understanding of cellular and molecular mechanisms underlying clubroot non-host and host resistance.

Two Clubroot-Resistance Genes, Rcr3 and Rcr9wa, Mapped in Brassica rapa Using Bulk Segregant RNA Sequencing

Genetic resistance is widely used to manage clubroot (Plasmodiophora brassicae) in brassica crops, but new pathotypes have recently been identified on canola (Brassica napus) on the Canadian prairies. Resistance effective against both the most prevalent pathotype (3H, based on the Canadian Clubroot Differential system) and the new pathotypes is needed. BC₁ plants of Brassica rapa from a cross of line 96-6990-2 (clubroot resistance originating from turnip cultivar ‘Waaslander’) and a susceptible doubled-haploid line, ACDC, exhibited a 1:1 segregation for resistance against pathotypes 3H and 5X. A resistance gene designated as Rcr3 was mapped initially based on the percentage of polymorphic variants using bulked segregant RNA sequencing (BSR-Seq) and further mapped using Kompetitive Allele Specific PCR. DNA variants were identified by assembling short reads against a reference genome of B. rapa. Rcr3 was mapped into chromosome A08. It was flanked by single nucleotide polymorphisms (SNP) markers (A90_A08_SNP_M12 and M16) between 10.00 and 10.23 Mb, in an interval of 231.6 Kb. There were 32 genes in the Rcr3 interval. Three genes (Bra020951, Bra020974, and Bra020979) were annotated with disease resistance mechanisms, which are potential candidates for Rcr3. Another resistance gene, designated as Rcr9^wa, for resistance to pathotype 5X was mapped, with the flanking markers (A90_A08_SNP_M28 and M79) between 10.85 and 11.17 Mb using the SNP sites identified through BSR-Seq for Rcr3. There were 44 genes in the Rcr9^wa interval, three of which (Bra020827, Bra020828, Bra020814) were annotated as immune-system-process related genes, which are potential candidates for Rcr9^wa.

QTL Mapping and Inheritance of Clubroot Resistance Genes Derived From Brassica rapa subsp. rapifera (ECD 02) Reveals Resistance Loci and Distorted Segregation Ratios in Two F 2 Populations of Different Crosses

In this study, Brassica rapa subsp. rapifera (ECD 02) which exhibits broad-spectrum resistance to many Canadian Plasmodiophora brassicae isolates was crossed with two clubroot-susceptible B. rapa accessions to produce two F ₂ populations. The F ₂ plants were screened against P. brassicae pathotypes 3H, 5X, and 5G. The Chi-square goodness of fit test showed that the vast majority (≈75%) of the crosses that produced the F ₂ populations showed segregation ratios of 9R:7S, 7R:9S, 13R:3S, 3R:13S, 5R:11S, 11R:5S, and 1R:15S. These were modifications of the 15R:1S ratio expected for the inheritance of two dominant major clubroot resistance (CR) genes from ECD 02. The distorted segregation ratios suggest that the two resistance genes are on different chromosomes and that two genes interact in an epistatic manner to confer resistance. Genotyping was conducted with 144 PCR-based markers in the two F ₂ populations. Linkage and QTL analysis with the polymorphic markers identified two QTLs on chromosome A03 to be associated with resistance to P. brassicae pathotypes 5X and 5G in Popl#1 while only the second QTL on chromosome A03 was associated with resistance to pathotypes 5X and 5G in Popl#2. The QTLs clustered in genomic regions on the A03 chromosome of B. rapa where the CRa/CRb ^Kato gene(s) are mapped. In addition, the Crr1 gene on the A08 chromosome of B. rapa was detected in the two F ₂ populations. Therefore, the phenotypic and molecular data confirm the existence of two CR genes in ECD 02. This is the first study that shows that major dominant genes in Brassica interact in a non-additive manner to confer resistance to different P. brassicae pathotypes. Key Message: This study provides knowledge on the inheritance and type of gene action for clubroot resistance derived from Brassica rapa subsp. rapifera (ECD 02). The results indicated that duplicate recessive and recessive suppression epistatic interactions, digenic additivity and complementary gene action between the CRa/CRb ^Kato gene(s) on the A03 and the Crr1 gene on the A08 chromosome of B. rapa controlled clubroot resistance to P. brassicae pathotypes 3H, 5X and 5G.

Utilization of Ogura CMS germplasm with the clubroot resistance gene by fertility restoration and cytoplasm replacement in Brassica oleracea L

Clubroot disease, a major plant root disease caused by Plasmodiophora brassicae, has become one of the most destructive diseases among cultivated cruciferous vegetables. However, clubroot-resistant Brassica oleracea materials are rare. A few clubroot-resistant cabbage varieties are available on the market, but all are Ogura cytoplasmic male sterile (CMS) types. Therefore, in this study, to reutilize the clubroot-resistant Ogura CMS germplasm of cabbage, a new fertility-restored Ogura CMS material, 16Q2-11, was used as a bridge to transfer the clubroot resistance (CR) gene from the Ogura CMS cytoplasm to the normal cytoplasm by a two-step method (a fertility restoration and cytoplasm replacement method). In the first cross for fertility restoration of Ogura CMS clubroot-resistant cabbage (FRCRC), 16Q2-11 was used as a restorer to cross with Ogura CMS materials containing the CR gene CRb2. Eleven Rfo-positive progenies were generated, of which four contained CRb2: F8-514, F8-620, F8-732 and F8-839. After inoculation with race 4 of P. brassicae, these four CRb2-positive individuals showed resistance. Furthermore, F8-514 and F8-839 were then used as male parents in the second cross of FRCRC to cross with cabbage inbred lines, resulting in the successful introgression of the CRb2 gene into the inbred lines. All offspring produced from this step of cross, which had a normal cytoplasm, showed a high resistance to race 4 of P. brassicae and could be utilized for the breeding of clubroot-resistant cabbage varieties in the future. This is the first time that the Ogura CMS restorer has been used to restore the fertility of Ogura CMS clubroot-resistant cabbages, which could improve germplasm diversity in cabbage and provide a reference method for using CMS germplasm in Brassica crops.

An update on the arsenal: mining resistance genes for disease management of Brassica crops in the genomic era

Brassica species include many economically important crops that provide nutrition and health-promoting substances to humans worldwide. However, as with all crops, their production is constantly threatened by emerging viral, bacterial, and fungal diseases, whose incidence has increased in recent years. Traditional methods of control are often costly, present limited effectiveness, and cause environmental damage; instead, the ideal approach is to mine and utilize the resistance genes of the Brassica crop hosts themselves. Fortunately, the development of genomics, molecular genetics, and biological techniques enables us to rapidly discover and apply resistance (R) genes. Herein, the R genes identified in Brassica crops are summarized, including their mapping and cloning, possible molecular mechanisms, and application in resistance breeding. Future perspectives concerning how to accurately discover additional R gene resources and efficiently utilize these genes in the genomic era are also discussed.

Comparative Transcriptome Analysis Reveals Key Pathways and Hub Genes in Rapeseed During the Early Stage of Plasmodiophora brassicae Infection

Rapeseed (Brassica napus L., AACC, 2n = 38) is one of the most important oil crops around the world. With intensified rapeseed cultivation, the incidence and severity of clubroot infected by Plasmodiophora brassicae Wor. (P. brassicae) has increased very fast, which seriously impedes the development of rapeseed industry. Therefore, it is very important and timely to investigate the mechanisms and genes regulating clubroot resistance (CR) in rapeseed. In this study, comparative transcriptome analysis was carried out on two rapeseed accessions of R- (resistant) and S- (susceptible) line. Three thousand one hundred seventy-one and 714 differentially expressed genes (DEGs) were detected in the R- and S-line compared with the control groups, respectively. The results indicated that the CR difference between the R- and S-line had already shown during the early stage of P. brassicae infection and the change of gene expression pattern of R-line exhibited a more intense defensive response than that of S-line. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of 2,163 relative-DEGs, identified between the R- and S-line, revealed that genes participated in plant hormone signal transduction, fatty acid metabolism, and glucosinolate biosynthesis were involved in regulation of CR. Further, 12 hub genes were identified from all relative-DEGs with the help of weighted gene co-expression network analysis. Haplotype analysis indicated that the natural variations in the coding regions of some hub genes also made contributed to CR. This study not only provides valuable information for CR molecular mechanisms, but also has applied implications for CR breeding in rapeseed.

Integrating long noncoding RNAs and mRNAs expression profiles of response to Plasmodiophora brassicae infection in Pakchoi (Brassica campestris ssp. chinensis Makino)

The biotrophic protist Plasmodiophora brassicae causes serious damage to Brassicaceae crops grown worldwide. However, the molecular mechanism of the Brassica rapa response remains has not been determined. Long noncoding RNA and mRNA expression profiles in response to Plasmodiophora brassicae infection were investigated using RNA-seq on the Chinese cabbage inbred line C22 infected with P. brassicae. Approximately 5,193 mRNAs were significantly differentially expressed, among which 1,345 were upregulated and 3,848 were downregulated. The GO enrichment analysis shows that most of these mRNAs are related to the defense response. Meanwhile, 114 significantly differentially expressed lncRNAs were identified, including 31 upregulated and 83 downregulated. Furthermore, a total of 2,344 interaction relationships were detected between 1,725 mRNAs and 103 lncRNAs with a correlation coefficient greater than 0.8. We also found 15 P. brassicaerelated mRNAs and 16 lncRNA interactions within the correlation network. The functional annotation showed that 15 mRNAs belong to defense response proteins (66.67%), protein phosphorylation (13.33%), root hair cell differentiation (13.33%) and regulation of salicylic acid biosynthetic process (6.67%). KEGG annotation showed that the vast majority of these genes are involved in the biosynthesis of secondary metabolism pathways and plant-pathogen interactions. These results provide a new perspective on lncRNA-mRNA network function and help to elucidate the molecular mechanism of P. brassicae infection.

Resolution of quantitative resistance to clubroot into QTL-specific metabolic modules

Plant disease resistance is often under quantitative genetic control. Thus, in a given interaction, plant cellular responses to infection are influenced by resistance or susceptibility alleles at different loci. In this study, a genetic linkage analysis was used to address the complexity of the metabolic responses of Brassica napus roots to infection by Plasmodiophora brassicae. Metabolome profiling and pathogen quantification in a segregating progeny allowed a comparative mapping of quantitative trait loci (QTLs) involved in resistance and in metabolic adjustments. Distinct metabolic modules were associated with each resistance QTL, suggesting the involvement of different underlying cellular mechanisms. This approach highlighted the possible role of gluconasturtiin and two unknown metabolites in the resistance conferred by two QTLs on chromosomes C03 and C09, respectively. Only two susceptibility biomarkers (glycine and glutathione) were simultaneously linked to the three main resistance QTLs, suggesting the central role of these compounds in the interaction. By contrast, several genotype-specific metabolic responses to infection were genetically unconnected to resistance or susceptibility. Likewise, variations of root sugar profiles, which might have influenced pathogen nutrition, were not found to be related to resistance QTLs. This work illustrates how genetic metabolomics can help to understand plant stress responses and their possible links with disease.

Analysing the genetic architecture of clubroot resistance variation in Brassica napus by associative transcriptomics

Clubroot is a destructive soil-borne pathogen of Brassicaceae that causes significant recurrent reductions in yield of cruciferous crops. Although there is some resistance in oilseed rape (a crop type of the species Brassica napus), the genetic basis of that resistance is poorly understood. In this study, we used an associative transcriptomics approach to elucidate the genetic basis of resistance to clubroot pathotype ECD 17/31/31 across a genetic diversity panel of 245 accessions of B. napus. A single nucleotide polymorphism (SNP) association analysis was performed with 256,397 SNPs distributed across the genome of B. napus and combined with transcript abundance data of 53,889 coding DNA sequence (CDS) gene models. The SNP association analysis identified two major loci (on chromosomes A2 and A3) controlling resistance and seven minor loci. Within these were a total of 86 SNP markers. Altogether, 392 genes were found in these regions. Another 21 genes were implicated as potentially involved in resistance using gene expression marker (GEM) analysis. After GO enrichment analysis and InterPro functional analysis of the identified genes, 82 candidate genes were identified as having roles in clubroot resistance. These results provide useful information for marker-assisted breeding which could lead to acceleration of pyramiding of multiple clubroot resistance genes in new varieties.

Understanding the Resistance Mechanism in Brassica napus to Clubroot Caused by Plasmodiophora brassicae

Exploring the mechanism of plant resistance has become the basis for selection of resistance varieties but reports on revealing resistant mechanism in Brassica napus against Plasmodiophora brassicae are rare. In this study, RNA-seq was conducted in the clubroot-resistant B. napus breeding line ZHE-226 and in the clubroot-susceptible rapeseed cultivar Zhongshuang 11 at 0, 3, 6, 9, and 12 days after inoculation. Strong alteration was detected specifically in ZHE-226 as soon as the root hair infection happened, and significant promotion was found in ZHE-226 on cell division or cell cycle, DNA repair and synthesis, protein synthesis, signaling, antioxidation, and secondary metabolites. Combining results from physiological, biochemical, and histochemical assays, our study highlights an effective signaling in ZHE-226 in response to P. brassicae. This response consists of a fast initiation of receptor kinases by P. brassicae; the possible activation of host intercellular G proteins which might, together with an enhanced Ca²⁺ signaling, promote the production of reactive oxygen species; and programmed cell death in the host. Meanwhile, a strong ability to maintain homeostasis of auxin and cytokinin in ZHE-226 might effectively limit the formation of clubs on host roots. Our study provides initial insights into resistance mechanism in rapeseed to P. brassicae.

Quantitative resistance to clubroot infection mediated by transgenerational epigenetic variation in Arabidopsis

Quantitative disease resistance, often influenced by environmental factors, is thought to be the result of DNA sequence variants segregating at multiple loci. However, heritable differences in DNA methylation, so-called transgenerational epigenetic variants, also could contribute to quantitative traits. Here, we tested this possibility using the well-characterized quantitative resistance of Arabidopsis to clubroot, a Brassica major disease caused by Plasmodiophora brassicae. For that, we used the epigenetic recombinant inbred lines (epiRIL) derived from the cross ddm1-2 × Col-0, which show extensive epigenetic variation but limited DNA sequence variation. Quantitative loci under epigenetic control (QTL^epi ) mapping was carried out on 123 epiRIL infected with P. brassicae and using various disease-related traits. EpiRIL displayed a wide range of continuous phenotypic responses. Twenty QTL^epi were detected across the five chromosomes, with a bona fide epigenetic origin for 16 of them. The effect of five QTL^epi was dependent on temperature conditions. Six QTL^epi co-localized with previously identified clubroot resistance genes and QTL in Arabidopsis. Co-localization of clubroot resistance QTL^epi with previously detected DNA-based QTL reveals a complex model in which a combination of allelic and epiallelic variations interacts with the environment to lead to variation in clubroot quantitative resistance.

Comparative Transcriptome Analysis Reveals Key Pathways and Hub Genes in Rapeseed During the Early Stage of Plasmodiophora brassicae Infection

Rapeseed (Brassica napus L., AACC, 2n = 38) is one of the most important oil crops around the world. With intensified rapeseed cultivation, the incidence and severity of clubroot infected by Plasmodiophora brassicae Wor. (P. brassicae) has increased very fast, which seriously impedes the development of rapeseed industry. Therefore, it is very important and timely to investigate the mechanisms and genes regulating clubroot resistance (CR) in rapeseed. In this study, comparative transcriptome analysis was carried out on two rapeseed accessions of R- (resistant) and S- (susceptible) line. Three thousand one hundred seventy-one and 714 differentially expressed genes (DEGs) were detected in the R- and S-line compared with the control groups, respectively. The results indicated that the CR difference between the R- and S-line had already shown during the early stage of P. brassicae infection and the change of gene expression pattern of R-line exhibited a more intense defensive response than that of S-line. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of 2,163 relative-DEGs, identified between the R- and S-line, revealed that genes participated in plant hormone signal transduction, fatty acid metabolism, and glucosinolate biosynthesis were involved in regulation of CR. Further, 12 hub genes were identified from all relative-DEGs with the help of weighted gene co-expression network analysis. Haplotype analysis indicated that the natural variations in the coding regions of some hub genes also made contributed to CR. This study not only provides valuable information for CR molecular mechanisms, but also has applied implications for CR breeding in rapeseed.

Analysis of genome-wide variants through bulked segregant RNA sequencing reveals a major gene for resistance to Plasmodiophora brassicae in Brassica oleracea

Two cabbage (Brassica oleracea) cultivars 'Tekila' and 'Kilaherb' were identified as resistant to several pathotypes of Plasmodiophora brassicae. In this study, we identified a clubroot resistance gene (Rcr7) in 'Tekila' for resistance to pathotype 3 of P. brassicae from a segregating population derived from 'Tekila' crossed with the susceptible line T010000DH3. Genetic mapping was performed by identifying the percentage of polymorphic variants (PPV), a new method proposed in this study, through bulked segregant RNA sequencing. Chromosome C7 carried the highest PPV (42%) compared to the 30-34% in the remaining chromosomes. A peak with PPV (56-73%) was found within the physical interval 41-44 Mb, which indicated that Rcr7 might be located in this region. Kompetitive Allele-Specific PCR was used to confirm the association of Rcr7 with SNPs in the region. Rcr7 was flanked by two SNP markers and co-segregated with three SNP markers in the segregating population of 465 plants. Seven genes encoding TIR-NBS-LRR disease resistance proteins were identified in the target region, but only two genes, Bo7g108760 and Bo7g109000, were expressed. Resistance to pathotype 5X was also mapped to the same region as Rcr7. B. oleracea lines including 'Kilaherb' were tested with five SNP markers for Rcr7 and for resistance to pathotype 3; 11 of 25 lines were resistant, but 'Kilaherb' was the only line that carried the SNP alleles associated with Rcr7. The presence of Rcr7 in 'Kilaherb' for resistance to both pathotypes 3 and 5X was confirmed through linkage analysis.

Combinations of Independent Dominant Loci Conferring Clubroot Resistance in All Four Turnip Accessions (Brassica rapa) From the European Clubroot Differential Set

Clubroot disease is devastating to Brassica crop production when susceptible cultivars are planted in infected fields. European turnips are the most resistant sources and their resistance genes have been introduced into other crops such oilseed rape (Brassica napus L.), Chinese cabbage and other Brassica vegetables. The European clubroot differential (ECD) set contains four turnip accessions (ECD1-4). These ECD turnips exhibited high levels of resistance to clubroot when they were tested under controlled environmental conditions with Canadian field isolates. Gene mapping of the clubroot resistance genes in ECD1-4 were performed and three independent dominant resistance loci were identified. Two resistance loci were mapped on chromosome A03 and the third on chromosome A08. Each ECD turnip accession contained two of these three resistance loci. Some resistance loci were homozygous in ECD accessions while others showed heterozygosity based on the segregation of clubroot resistance in 20 BC₁ families derived from ECD1 to 4. Molecular markers were developed linked to each clubroot resistance loci for the resistance gene introgression in different germplasm.

Identification of Quantitative Trait Loci for Clubroot Resistance in Brassica oleracea With the Use of Brassica SNP Microarray

Increasing clubroot resistance (CR) of Brassica oleracea by ascertaining the molecular mechanisms has been the key focus in modern B. oleracea breeding. In order to identify the quantitative trait loci (QTLs) associated with CR in B. oleracea, 94 F2 vegetative lines which were developed by tissue culture of selfed seeds from the F1 generation between a clubroot-resistant B. oleracea inbred line and a susceptible line, were identified for disease incidence and six CR-associated traits under a lab inoculation by Plasmodiophora brassicae and were genotyped with the 60K Brassica SNP array. Significant correlations were detected for numbers of fibrous roots and P. brassicae content in roots with disease incidence. Nine linkage groups were constructed from 565 bins which covered around 3,000 SNPs, spanning 1,028 cM of the B. oleracea genome with an average distance of 1.82 cM between adjacent bins. A total of 23 QTLs were identified for disease incidence and the other two correlated traits, individually explaining 6.1-17.8% of the phenotypic variation. Several overlaps were detected among traits, including one three-traits-overlapped locus on linkage group C08 and two important overlapped regions between the two CR-associated traits on C06. The QTLs were compared with known CR loci/genes and the novelty of our QTLs was discussed.

Resynthesis of Brassica juncea for resistance to Plasmodiophora brassicae pathotype 3

The oilseed crop Brassica juncea carries many desirable traits; however, resistance to clubroot disease, caused by Plasmodiophora brassicae, is not available in this species. We are the first to report the clubroot resistant resynthesized B. juncea lines, developed through interspecific crosses between a clubroot resistant B. rapa ssp. rapifera and two susceptible B. nigra lines, and the stability of the resistance in self-pollinated generations. The interspecific nature of the resynthesized B. juncea plants was confirmed by using A- and B-genome specific SSR markers, and flow cytometric analysis of nuclear DNA content. Self-pollinated progeny (S₁ and S₂) of the resynthesized B. juncea plants were evaluated for resistance to P. brassicae pathotype 3. The S₁ and S₂ progenies of one of the resynthesized B. juncea lines were resistant to this pathotype. However, resistance was lost in 6 to 13% plants of the S₂ progenies derived from the second resynthesized B. juncea line; this apparently resulted from the loss of the genomic region carrying resistance due to meiotic anomalies.

Screening of clubroot-resistant varieties and transfer of clubroot resistance genes to Brassica napus using distant hybridization

Clubroot is an economically important disease affecting plants in the family Cruciferae worldwide. In this study, a collection of 50 Cruciferae accessions was screened using Plasmodiophora brassicae pathotype 4 in China. Eight of these demonstrated resistance, including three Chinese cabbages, two cabbages, one radish, one kale, and one Brassica juncea. The three clubroot-resistant Chinese cabbages (1003, 1007 and 1008) were then used to transfer the clubroot resistance genes to B. napus by distant hybridization combined with embryo rescue. Three methods including morphological identification, cytology identification, and molecular marker-assisted selection were used to determine hybrid authenticity, and 0, 2, and 4 false hybrids were identified by these three methods, respectively. In total, 297 true hybrids were identified. Clubroot resistance markers and artificial inoculation were utilized to determine the source of clubroot resistance in the true hybrids. As a result, two simple sequence repeat (SSR) and two intron polymorphic (IP) markers linked to clubroot resistance genes were identified, the clubroot resistance genes of 1007 and 1008 were mapped to A03. At last, 159 clubroot-resistant hybrids were obtained by clubroot resistance markers and artificial inoculation. These intermediate varieties will be used as the 'bridge material' of clubroot resistance for further B. napus breeding.

Genome-wide identification and role of MKK and MPK gene families in clubroot resistance of Brassica rapa

Mitogen-activated protein kinase (MAPK or MPK) cascades play key roles in responses to various biotic stresses, as well as in plant growth and development. However, the responses of MPK and MPK kinase (MKK) in Chinese cabbage (Brassica rapa ssp. pekinensis) to Plasmodiophora brassicae, a causal agent of clubroot disease in Brassica crops, are still not clear. In the present study, a total of 11 B. rapa MKK (BraMKK) and 30 BraMPK genes were identified and unevenly distributed in 6 and 10 chromosomes, respectively. The synteny analysis indicated that these genes experienced whole-genome triplication and segmental and tandem duplication during or after the divergence of B. rapa, accompanied by the loss of three MKK and two MPK orthologs of Arabidopsis. The BraMKK and BraMPK genes were classified into four groups with similar intron/exon structures and conserved motifs in each group. A quantitative PCR analysis showed that the majority of BraMKK and BraMPK genes were natively expressed in roots, hypocotyls, and leaves, whereas 5 BraMKK and 16 BraMPK genes up-regulated in the roots upon P. brassicae infection. Additionally, these 5 BraMKK and 16 BraMPK genes exhibited a significantly different expression pattern between a pair of clubroot-resistant/susceptible near-isogenic lines (NILs). Furthermore, the possible modules of MKK-MPK involved in B. rapa-P. brassicae interaction are also discussed. The present study will provide functional clues for further characterization of the MAPK cascades in B. rapa.

Combinations of Independent Dominant Loci Conferring Clubroot Resistance in All Four Turnip Accessions (Brassica rapa) From the European Clubroot Differential Set

Clubroot disease is devastating to Brassica crop production when susceptible cultivars are planted in infected fields. European turnips are the most resistant sources and their resistance genes have been introduced into other crops such oilseed rape (Brassica napus L.), Chinese cabbage and other Brassica vegetables. The European clubroot differential (ECD) set contains four turnip accessions (ECD1-4). These ECD turnips exhibited high levels of resistance to clubroot when they were tested under controlled environmental conditions with Canadian field isolates. Gene mapping of the clubroot resistance genes in ECD1-4 were performed and three independent dominant resistance loci were identified. Two resistance loci were mapped on chromosome A03 and the third on chromosome A08. Each ECD turnip accession contained two of these three resistance loci. Some resistance loci were homozygous in ECD accessions while others showed heterozygosity based on the segregation of clubroot resistance in 20 BC₁ families derived from ECD1 to 4. Molecular markers were developed linked to each clubroot resistance loci for the resistance gene introgression in different germplasm.

Evaluating Changes in Cell-Wall Components Associated with Clubroot Resistance Using Fourier Transform Infrared Spectroscopy and RT-PCR

Clubroot disease is a serious threat to canola production in western Canada and many parts of the world. Rcr1 is a clubroot resistance (CR) gene identified recently and its molecular mechanisms in mediating CR have been studied using several omics approaches. The current study aimed to characterize the biochemical changes in the cell wall of canola roots connecting to key molecular mechanisms of this CR gene identified in prior studies using Fourier transform infrared (FTIR) spectroscopy. The expression of nine genes involved in phenylpropanoid metabolism was also studied using qPCR. Between susceptible (S) and resistance (R) samples, the most notable biochemical changes were related to an increased biosynthesis of lignin and phenolics. These results were supported by the transcription data on higher expression of BrPAL1. The up-regulation of PAL is indicative of an inducible defence response conferred by Rcr1; the activation of this basal defence gene via the phenylpropanoid pathway may contribute to clubroot resistance conferred by Rcr1. The data indicate that several cell-wall components, including lignin and pectin, may play a role in defence responses against clubroot. Principal components analysis of FTIR data separated non-inoculated samples from inoculated samples, but not so much between inoculated S and inoculated R samples. It is also shown that FTIR spectroscopy can be a useful tool in studying plant-pathogen interaction at cellular levels.

Clubroot resistance QTL are modulated by nitrogen input in Brassica napus

Nitrogen levels can modulate the effectiveness of clubroot resistance in an isolate- and host-specific manner. While the same QTL were detected under high and low nitrogen, their effects were altered. Clubroot, caused by Plasmodiophora brassicae, is one of the most damaging diseases of oilseed rape and is known to be affected by nitrogen fertilization. However, the genetic factors involved in clubroot resistance have not been characterized under nitrogen-limiting conditions. This study aimed to assess the variability of clubroot resistance under different nitrogen levels and to characterize the impact of nitrogen supply on genetic resistance factors. Linkage analyses and a genome-wide association study were conducted to detect QTL for clubroot resistance and evaluate their sensitivity to nitrogen. The clubroot response of a set of 92 diverse oilseed rape accessions and 108 lines derived from a cross between 'Darmor-bzh' (resistant) and 'Yudal' (susceptible) was studied in the greenhouse under high- and low-nitrogen conditions, following inoculation with the P. brassicae isolates eH and K92-16. Resistance to each isolate was controlled by a major QTL and a few small-effects QTL. While the same QTL were detected under both high and low nitrogen, their effects were altered. Clubroot resistance to isolate eH, but not K92-16, was greater under a low-N supply versus a high-N supply. New sources of resistance were found among the oilseed rape accessions under both low and high-N conditions. The results are discussed relative to the literature and from a crop improvement perspective.

Comparative Transcriptome Analysis between Broccoli (Brassica oleracea var. italica) and Wild Cabbage (Brassica macrocarpa Guss.) in Response to Plasmodiophora brassicae during Different Infection Stages

Clubroot, one of the most devastating diseases to the Brassicaceae family, is caused by the obligate biotrophic pathogen Plasmodiophora brassicae. However, studies of the molecular basis of disease resistance are still poor especially in quantitative resistance. In the present paper, two previously identified genotypes, a clubroot-resistant genotype (wild cabbage, B2013) and a clubroot-susceptible genotype (broccoli, 90196) were inoculated by P. brassicae for 0 (T0), 7 (T7), and 14 (T14) day after inoculation (DAI). Gene expression pattern analysis suggested that response changes in transcript level of two genotypes under P. brassicae infection were mainly activated at the primary stage (T7). Based on the results of DEGs functional enrichments from two infection stages, genes associated with cell wall biosynthesis, glucosinolate biosynthesis, and plant hormone signal transduction showed down-regulated at T14 compared to T7, indicating that defense responses to P. brassicae were induced earlier, and related pathways were repressed at T14. In addition, the genes related to NBS-LRR proteins, SA signal transduction, cell wall and phytoalexins biosynthesis, chitinase, Ca²⁺ signals and RBOH proteins were mainly up-regulated in B2013 by comparing those of 90196, indicating the pathways of response defense to clubroot were activated in the resistant genotype. This is the first report about comparative transcriptome analysis for broccoli and its wild relative during the different stages of P. brassicae infection and the results should be useful for molecular assisted screening and breeding of clubroot-resistant genotypes.

Pathotype Classification of Plasmodiophora brassicae Isolates Using Clubroot-Resistant Cultivars of Chinese Cabbage

Clubroot disease caused by Plasmodiophora brassicae is one of the most serious diseases in Brassica crops worldwide. In this study, the pathotypes of 12 Korean P. brassicae field isolates were determined using various Chinese cabbage including 22 commercial cultivars from Korea, China, and Japan, and 15 inbred lines. All P. brassicae isolates exhibited the typical clubroot disease on non-clubroot resistant cultivar, indicating that the isolates were highly pathogenic. According to the reactions on the Williams' hosts, the 12 field isolates were initially classified into five races. However, when these isolates were inoculated onto clubroot-resistant (CR) cultivars of Chinese cabbage, several isolates led to different disease responses even though the isolates have been assigned to the same race by the Williams' host responses. Based on the pathogenicity results, the 12 field isolates were reclassified into four different groups: pathotype 1 (GN1, GN2, GS, JS, and HS), 2 (DJ and KS), 3 (HN1, PC, and YC), and 4 (HN2 and SS). In addition, the CR cultivars from Korea, China, and Japan exhibited distinguishable disease responses to the P. brassicae isolates, suggesting that the 22 cultivars used in this study, including the non-CR cultivars, are classified into four different host groups based on their disease resistance. Combining these findings, the four differential hosts of Chinese cabbage and four pathotype groups of P. brassicae might provide an efficient screening system for resistant cultivars and a new foundation of breeding strategies for CR Chinese cabbage.

A Genome-Wide Association Study Reveals New Loci for Resistance to Clubroot Disease in Brassica napus

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. However, the yield and quality of rapeseed were largely decreased by clubroot (Plasmodiophora brassicae Woronin). Therefore, it is of great importance for screening more resistant germplasms or genes and improving the resistance to P. brassicae in rapeseed breeding. In this study, a massive resistant identification for a natural global population was conducted in two environments with race/pathotype 4 of P. brassicae which was the most predominant in China, and a wide range of phenotypic variation was found in the population. In addition, a genome-wide association study of 472 accessions for clubroot resistance (CR) was performed with 60K Brassica Infinium SNP arrays for the first time. In total, nine QTLs were detected, seven of which were novel through integrative analysis. Furthermore, additive effects in genetic control of CR in rapeseed among the above loci were found. By bioinformatic analyses, the candidate genes of these loci were predicted, which indicated that TIR-NBS gene family might play an important role in CR. It is believable that the results presented in our study could provide valuable information for understanding the genetic mechanism and molecular regulation of CR.

A Genome-Wide Association Study Reveals New Loci for Resistance to Clubroot Disease in Brassica napus

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. However, the yield and quality of rapeseed were largely decreased by clubroot (Plasmodiophora brassicae Woronin). Therefore, it is of great importance for screening more resistant germplasms or genes and improving the resistance to P. brassicae in rapeseed breeding. In this study, a massive resistant identification for a natural global population was conducted in two environments with race/pathotype 4 of P. brassicae which was the most predominant in China, and a wide range of phenotypic variation was found in the population. In addition, a genome-wide association study of 472 accessions for clubroot resistance (CR) was performed with 60K Brassica Infinium SNP arrays for the first time. In total, nine QTLs were detected, seven of which were novel through integrative analysis. Furthermore, additive effects in genetic control of CR in rapeseed among the above loci were found. By bioinformatic analyses, the candidate genes of these loci were predicted, which indicated that TIR-NBS gene family might play an important role in CR. It is believable that the results presented in our study could provide valuable information for understanding the genetic mechanism and molecular regulation of CR.

Shotgun Label-free Proteomic Analysis of Clubroot (Plasmodiophora brassicae) Resistance Conferred by the Gene Rcr1 in Brassica rapa

Clubroot, caused by the plasmodiophorid pathogen Plasmodiophora brassicae, is one of the most serious diseases on Brassica crops worldwide and a major threat to canola production in western Canada. Host resistance is the key strategy for clubroot management on canola. Several clubroot resistance (CR) genes have been identified, but the mechanisms associated with these CR genes are poorly understood. In the current study, a label-free shotgun proteomic approach was used to profile and compare the proteomes of Brassica rapa carrying and not carrying the CR gene Rcr1 in response to P. brassicae infection. A total of 527 differentially accumulated proteins (DAPs) were identified between the resistant (with Rcr1) and susceptible (without Rcr1) samples, and functional annotation of these DAPs indicates that the perception of P. brassicae and activation of defense responses are triggered via an unique signaling pathway distinct from common modes of recognition receptors reported with many other plant-pathogen interactions; this pathway appears to act in a calcium-independent manner through a not-well-defined cascade of mitogen-activated protein kinases and may require the ubiquitin-26S proteasome found to be related to abiotic stresses, especially the cold-stress tolerance in other studies. Both up-regulation of defense-related and down-regulation of pathogenicity-related metabolism was observed in plants carrying Rcr1, and these functions may all contribute to the CR mediated by Rcr1. These results, combined with those of transcriptomic analysis reported earlier, improved our understanding of molecular mechanisms associated with Rcr1 and CR at large, and identified candidate metabolites or pathways related to specific resistance mechanisms. Deploying CR genes with different modes of action may help improve the durability of CR.

Genetics and molecular mapping of resistance to Plasmodiophora brassicae pathotypes 2, 3, 5, 6, and 8 in rutabaga (Brassica napus var. napobrassica)

Clubroot disease, caused by Plasmodiophora brassicae, is a threat to the production of Brassica crops including oilseed B. napus. In Canada, several pathotypes of this pathogen, such as pathotypes 2, 3, 5, 6, and 8, were identified, and resistance to these pathotypes was found in a rutabaga (B. napus var. napobrassica) genotype. In this paper, we report the genetic basis and molecular mapping of this resistance by use of F₂, backcross (BC₁), and doubled haploid (DH) populations generated from crossing of this rutabaga line to a susceptible spring B. napus canola line. The F₁, F₂, and BC₁ populations were evaluated for resistance to pathotype 3, and the DH population was evaluated for resistance to pathotypes 2, 3, 5, 6, and 8. A 3:1 segregation in F₂ and a 1:1 segregation in BC₁ were found for resistance to pathotype 3, and a 1:1 segregation was found in the DH population for resistance to all pathotypes. Molecular mapping by using the DH population identified a genomic region on chromosome A8 carrying resistance to all five pathotypes. This suggests that a single gene or a cluster of genes, located in this genomic region, is involved in the control of resistance to these pathotypes.

Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L.)

Clubroot is a devastating disease caused by Plasmodiophora brassicae and results in severe losses of yield and quality in Brassica crops. Many clubroot resistance genes and markers are available in Brassica rapa but less is known in Brassica oleracea. Here, we applied the genotyping-by-sequencing (GBS) technique to construct a high-resolution genetic map and identify clubroot resistance (CR) genes. A total of 43,821 SNPs were identified using GBS data for two parental lines, one resistant and one susceptible lines to clubroot, and 18,187 of them showed >5× coverage in the GBS data. Among those, 4,103 were credibly genotyped for all 78 F₂ individual plants. These markers were clustered into nine linkage groups spanning 879.9 cM with an average interval of 1.15 cM. Quantitative trait loci (QTLs) survey based on three rounds of clubroot resistance tests using F_2:3 progenies revealed two and single major QTLs for Race 2 and Race 9 of P. brassicae, respectively. The QTLs show similar locations to the previously reported CR loci for Race 4 in B. oleracea but are in different positions from any of the CR loci found in B. rapa. We utilized two reference genome sequences in this study. The high-resolution genetic map developed herein allowed us to reposition 37 and 2 misanchored scaffolds in the 02–12 and TO1000DH genome sequences, respectively. Our data also support additional positioning of two unanchored 3.3 Mb scaffolds into the 02–12 genome sequence.

Resistance to Plasmodiophora brassicae in Brassica rapa and Brassica juncea genotypes From China

Clubroot disease, caused by Plasmodiophora brassicae Woronin, has become a major problem in cruciferous crops worldwide. Chinese cabbage (Brassica rapa), pak choi (B. rapa), and mustard (B. juncea) are important vegetable crops in China. Development of clubroot-resistant cultivars of these crops is urgently needed. In this study, 71 B. rapa and B. juncea genotypes from China, including cultivars and inbred lines, were evaluated for resistance to three P. brassicae pathotypes. A significant interaction was observed between the P. brassicae pathotypes and the Brassica genotypes. Pathotype 3, as defined on the differentials of Williams, exhibited the weakest virulence on all plant material. By contrast, pathotypes 5 and 6 were both highly pathogenic on most of the tested genotypes. In all, 10 of the 14 Chinese cabbage cultivars were resistant to all three pathotypes, while 4 were resistant only to a specific pathotype. Seven of eight progenies obtained from the selfing of Chinese cabbage cultivars were resistant to pathotype 3 but most were susceptible to pathotypes 5 and 6. Most inbred lines of Chinese cabbage and all inbred lines of pak choi and mustard were susceptible to all three pathotypes but their susceptibility was lower to pathotype 3 than to pathotypes 5 and 6.