Genotyping-by-sequencing reveals three QTL for clubroot resistance to six pathotypes of Plasmodiophora brassicae in Brassica rapa

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.

Multiomics analysis of a resistant European turnip ECD04 during clubroot infection reveals key hub genes underlying resistance mechanism

The clubroot disease has become a worldwide threat for crucifer crop production, due to its soil-borne nature and difficulty to eradicate completely from contaminated field. In this study we used an elite resistant European fodder turnip ECD04 and investigated its resistance mechanism using transcriptome, sRNA-seq, degradome and gene editing. A total of 1751 DEGs were identified from three time points after infection, among which 7 hub genes including XTH23 for cell wall assembly and two CPK28 genes in PTI pathways. On microRNA, we identified 17 DEMs and predicted 15 miRNA-target pairs (DEM-DEG). We validated two pairs (miR395-APS4 and miR160-ARF) by degradome sequencing. We investigated the miR395-APS4 pair by CRISPR-Cas9 mediated gene editing, the result showed that knocking-out APS4 could lead to elevated clubroot resistance in B. napus. In summary, the data acquired on transcriptional response and microRNA as well as target genes provide future direction especially gene candidates for genetic improvement of clubroot resistance on Brassica species.

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.

A CRISPR/Cas9-based vector system enables the fast breeding of selection-marker-free canola with Rcr1-rendered clubroot resistance

Breeding for disease resistance in major crops is of crucial importance for global food security and sustainability. However, common biotechnologies such as traditional transgenesis or genome editing do not provide an ideal solution, whereas transgenic crops free of selection markers such as cisgenic/intragenic crops might be suitable. In this study, after cloning and functional verification of the Rcr1 gene for resistance to clubroot (Plasmodiophora brassicae), we confirmed that the genes Rcr1, Rcr2, Rcr4, and CRa from Brassica rapa crops and the resistance gene from B. napus oilseed rape cv. 'Mendel' on chromosome A03 were identical in their coding regions. We also determined that Rcr1 has a wide distribution in Brassica breeding materials and renders potent resistance against multiple representative clubroot strains in Canada. We then modified a CRISPR/Cas9-based cisgenic vector system and found that it enabled the fast breeding of selection-marker-free transgenic crops with add-on traits, with selection-marker-free canola (B. napus) germplasms with Rcr1-rendered stable resistance to clubroot disease being successfully developed within 2 years. In the B. napus background, the intragenic vector system was able to remove unwanted residue sequences from the final product with high editing efficiency, and off-target mutations were not detected. Our study demonstrates the potential of applying this breeding strategy to other crops that can be transformed by Agrobacterium. Following the streamlined working procedure, intragenic germplasms can be developed within two generations, which could significantly reduce the breeding time and labor compared to traditional introgression whilst still achieving comparable or even better breeding results.

Identification of Clubroot (Plasmodiophora brassicae) Resistance Loci in Chinese Cabbage (Brassica rapa ssp. pekinensis) with Recessive Character

The soil-borne pathogen Plasmodiophora brassicae is the causal agent of clubroot, a major disease in Chinese cabbage (Brassica rapa ssp. pekinensis). The host's resistance genes often confer immunity to only specific pathotypes and may be rapidly overcome. Identification of novel clubroot resistance (CR) from germplasm sources is necessary. In this study, Bap246 was tested by being crossed with different highly susceptible B. rapa materials and showed recessive resistance to clubroot. An F2 population derived from Bap246 × Bac1344 was used to locate the resistance Quantitative Trait Loci (QTL) by Bulk Segregant Analysis Sequencing (BSA-Seq) and QTL mapping methods. Two QTL on chromosomes A01 (4.67-6.06 Mb) and A08 (10.42-11.43 Mb) were found and named Cr4Ba1.1 and Cr4Ba8.1, respectively. Fifteen and eleven SNP/InDel markers were used to narrow the target regions in the larger F2 population to 4.67-5.17 Mb (A01) and 10.70-10.84 Mb (A08), with 85 and 19 candidate genes, respectively. The phenotypic variation explained (PVE) of the two QTL were 30.97% and 8.65%, respectively. Combined with gene annotation, mutation site analysis, and real-time quantitative polymerase chain reaction (qRT-PCR) analysis, one candidate gene in A08 was identified, namely Bra020861. And an insertion and deletion (InDel) marker (co-segregated) named Crr1-196 was developed based on the gene sequence. Bra013275, Bra013299, Bra013336, Bra013339, Bra013341, and Bra013357 in A01 were the candidate genes that may confer clubroot resistance in Chinese cabbage. The resistance resource and the developed marker will be helpful in Brassica breeding programs.

Identification of QTLs for resistance to 10 pathotypes of Plasmodiophora brassicae in Brassica oleracea cultivar ECD11 through genotyping-by-sequencing

KEY MESSAGE

Two major quantitative trait loci (QTLs) and five minor QTLs for 10 pathotypes were identified on chromosomes C01, C03, C04 and C08 through genotyping-by-sequencing from Brassica oleracea. Clubroot caused by Plasmodiophora brassicae is an important disease in brassica crops. Managing clubroot disease of canola on the Canadian prairie is challenging due to the continuous emergence of new pathotypes. Brassica oleracea is considered a major source of quantitative resistance to clubroot. Genotyping-by-sequencing (GBS) was performed in the parental lines; T010000DH3 (susceptible), ECD11 (resistant) and 124 BC1 plants. A total of 4769 high-quality polymorphic SNP loci were obtained and distributed on 9 chromosomes of B. oleracea. Evaluation of 124 BC1S1 lines for resistance to 10 pathotypes: 3A, 2B, 5C, 3D, 5G, 3H, 8J, 5K, 5L and 3O of P. brassicae, was carried out. Seven QTLs, 5 originating from ECD11 and 2 from T010000DH3, were detected. One major QTL designated as Rcr_C03-1 on C03 contributed 16.0-65.6% of phenotypic variation explained (PVE) for 8 pathotypes: 2B, 5C, 5G, 3H, 8J, 5K, 5L and 3O. Another major QTL designated as Rcr_C08-1 on C08 contributed 8.3 and 23.5% PVE for resistance to 8J and 5K, respectively. Five minor QTLs designated as Rcr_C01-1, Rcr_C03-2, Rcr_C03-3, Rcr_C04-1 and Rcr_C08-2 were detected on chromosomes C01, C03, C04 and C08 that contributed 8.3-23.5% PVE for 5 pathotypes each of 3A, 2B, 3D, 8J and 5K. There were 1, 10 and 4 genes encoding TIR-NBS-LRR/CC-NBS-LRR class disease resistance proteins in the Rcr_C01-1, Rcr_C03-1 and Rcr_C08-1 flanking regions. The syntenic regions of the two major QTLs Rcr_C03-1 and Rcr_C08-1 in the B. rapa genome 'Chiifu' were searched.

Brassica napus Haploid and Double Haploid Production and Its Latest Applications

Rapeseed is one of the most important oil crops in the world. Increasing demand for oil and limited agronomic capabilities of present-day rapeseed result in the need for rapid development of new, superior cultivars. Double haploid (DH) technology is a fast and convenient approach in plant breeding as well as genetic research. Brassica napus is considered a model species for DH production based on microspore embryogenesis; however, the molecular mechanisms underlying microspore reprogramming are still vague. It is known that morphological changes are accompanied by gene and protein expression patterns, alongside carbohydrate and lipid metabolism. Novel, more efficient methods for DH rapeseed production have been reported. This review covers new findings and advances in Brassica napus DH production as well as the latest reports related to agronomically important traits in molecular studies employing the double haploid rapeseed lines.

The soil bacterial community regulates germination of Plasmodiophora brassicae resting spores rather than root exudates

Clubroot, caused by Plasmodiophora brassicae, is a severe soil-borne disease that restricts the production of cruciferous crops worldwide. A better understanding of biotic and abiotic factors regulating germination of P. brassicae resting spores in the soil is significant for developing novel control methods. Previous studies reported that root exudates can trigger P. brassicae resting spore germination, thus enabling a targeted attack of P. brassicae on host plant roots. However, we found that native root exudates collected under sterile conditions from host or non-host plants cannot stimulate the germination of sterile spores, indicating that root exudates may not be direct stimulation factors. Instead, our studies demonstrate that soil bacteria are essential for triggering germination. Through 16s rRNA amplicon sequencing analysis, we found that certain carbon sources and nitrate can reshape the initial microbial community to an inducing community leading to the germination of P. brassicae resting spores. The stimulating communities significantly differed in composition and abundance of bacterial taxa compared to the non-stimulating ones. Several enriched bacterial taxa in stimulating community were significantly correlated with spore germination rates and may be involved as stimulation factors. Based on our findings, a multi-factorial 'pathobiome' model comprising abiotic and biotic factors is proposed to represent the putative plant-microbiome-pathogen interactions associated with breaking spore dormancy of P. brassicae in soil. This study presents novel views on P. brassicae pathogenicity and lays the foundation for novel sustainable control strategies of clubroot.

Advances in Biological Control and Resistance Genes of Brassicaceae Clubroot Disease-The Study Case of China

Clubroot disease is a soil-borne disease caused by Plasmodiophora brassicae. It occurs in cruciferous crops exclusively, and causes serious damage to the economic value of cruciferous crops worldwide. Although different measures have been taken to prevent the spread of clubroot disease, the most fundamental and effective way is to explore and use disease-resistance genes to breed resistant varieties. However, the resistance level of plant hosts is influenced both by environment and pathogen race. In this work, we described clubroot disease in terms of discovery and current distribution, life cycle, and race identification systems; in particular, we summarized recent progress on clubroot control methods and breeding practices for resistant cultivars. With the knowledge of these identified resistance loci and R genes, we discussed feasible strategies for disease-resistance breeding in the future.

Marker-Assisted Pyramiding of CRa and CRd Genes to Improve the Clubroot Resistance of Brassica rapa

Clubroot, caused by Plasmodiophora brassicae, is an economically important soil-borne disease that threatens Brassicaceae crops worldwide. In recent years, the incidence area of Chinese cabbage (Brassica rapa ssp. pekinensis) clubroot disease has increased, which severely affects the yield and quality of Chinese cabbage. The resistance of varieties harboring the single clubroot-resistance (CR) gene is easily broken through by P. brassicae pathotypes. CRa and CRd, genetically identified in B. rapa, are CR genes known to be highly resistant to different P. brassicaea pathotypes. In our study, we perform the gene pyramiding of CRa and CRd in Chinese cabbages through marker-assisted selection (MAS), and develop homozygous pyramided lines. The newly generated pyramided lines exhibit greater resistance to six different pathotypes than that of two parental lines carrying a single CR gene. This study provides new CR-gene-pyramided lines for the development of clubroot-resistant Brassica varieties for future breeding programs.

Early-stage responses to Plasmodiophora brassicae at the transcriptome and metabolome levels in clubroot resistant and susceptible oilseed Brassica napus

Clubroot, a devastating soil-borne root disease, in Brassicaceae is caused by Plasmodiophora brassicae Woronin (P. brassicae W.), an obligate biotrophic protist. Plant growth and development, as well as seed yield of Brassica crops, are severely affected due to this disease. Several reports described the molecular responses of B. napus to P. brassicae; however, information on the early stages of pathogenesis is limited. In this study, we have used transcriptomics and metabolomics to characterize P. brassicae pathogenesis at 1-, 4-, and 7-days post-inoculation (DPI) in clubroot resistant (CR) and susceptible (CS) doubled-haploid (DH) canola lines. When we compared between inoculated and uninoculated groups, a total of 214 and 324 putative genes exhibited differential expression (q-value < 0.05) at one or more time-points in the CR and CS genotypes, respectively. When the inoculated CR and inoculated CS genotypes were compared, 4765 DEGs were differentially expressed (q-value < 0.05) at one or more time-points. Several metabolites related to organic acids (e.g., citrate, pyruvate), amino acids (e.g., proline, aspartate), sugars, and mannitol, were differentially accumulated in roots in response to pathogen infection when the CR and CS genotypes were compared. Several DEGs also corresponded to differentially accumulated metabolites, including pyrroline-5-carboxylate reductase (BnaC04g11450D), citrate synthase (BnaC02g39080D), and pyruvate kinase (BnaC04g23180D) as detected by transcriptome analysis. Our results suggest important roles for these genes in mediating resistance to clubroot disease. To our knowledge, this is the first report of an integrated transcriptome and metabolome analysis aimed at characterizing the molecular basis of resistance to clubroot in canola.

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.

Improvement of Resistance to Clubroot Disease in the Ogura CMS Restorer Line R2163 of Brassica napus

Oilseed rape (Brassica napus) has significant heterosis and Ogura CMS is a major way to use it. Ogura CMS has the advantages of complete and stable male sterility and easy-to-breed maintainers. Therefore, to breed better restorers has become an important goal for this system. Incidentally, clubroot is a soil-borne disease that is difficult to control by fungicidal chemicals, and it has been the main disease of oilseed rape in recent years in China, severely restricting the development of the oilseed rape industry. At present, the most effective method for controlling clubroot disease is to cultivate resistant varieties. One Ogura CMS restorer line (R2163) has shown much better combining ability, but lacks the clubroot disease resistance. This study was carried out to improve R2163 through marker-assisted backcross breeding (MABB). The resistant locus PbBa8.1 was introduced into the restorer R2163, and we then selected R2163R with clubroot disease resistance. Using the new restorer R2163R as the male parent and the sterile lines 116A and Z11A as the female parent, the improved, new resistant hybrids Kenyouza 741R and Huayouza 706R performed well, providing strong resistance and good agronomic traits. This work advances the utilization of heterosis and breeding for clubroot disease resistance in B. napus.

Evolutionary expansion and expression dynamics of cytokinin-catabolizing CKX gene family in the modern amphidiploid mustard (Brassica sp.)

Plant cytokinins (CKs) promote development and physiological processes, drought tolerance, root architecture, and ultimately crop productivity. Biologically active CKs (iP, tZ, and cZ) are precisely maintained in the vegetative and floral tissues through their irreversible degradation by developmentally regulated CK-catabolizing cytokinin oxidase/dehydrogenase (CKX) enzyme. A meta-analysis of CKX proteins was performed through an exhaustive exploration of multiple genome databases of cyanobacteria, bryophyte, monocot and eudicot plants to reveal the intricate evolutionary profiles of CKX enzymes specific to the family Brassicaceae. At least 175 unique paralogous/orthologous CKX sequences were successfully retrieved and phylogenetically clustered into distinct groups. Observations of structural divergences among paralogous sequences compared to their orthologs indicated that the progenitor CKX sequence had been subjected to massive structural modifications, possibly as a result of the evolutionary split between monocots and eudicots. An analysis of dN/dS comparisons of orthologous genes revealed that segmental CKX gene duplications have evolved primarily under purifying selection. Further, 24 CKX genes with conserved signature domain were identified in the amphidiploid Brassica juncea genome (AABB; 2n = 36). Genetic evolution of paralogous and orthologous genes was largely responsible for the expansion of CKX homoeologs in the amphidiploid Brassica genomes. Also, comparative analyses of 1.5 kb-long upstream regulatory regions of BjCKX genes identified various development- and stress-responsive elements. Spatial and temporal expression profiles of CKX genes were primarily attributed to their structural diversity observed in the 5'-regulatory regions along with species evolution. This data suggested that CKX duplicate genes had partitioned their spatial expression (= function) during evolution. These findings illustrated the evolutionary importance of CKX genes during plant development, and also suggested their deployment for future crop improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03294-0.

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.

R gene triplication confers European fodder turnip with improved clubroot resistance

Clubroot is one of the most important diseases for many important cruciferous vegetables and oilseed crops worldwide. Different clubroot resistance (CR) loci have been identified from only limited species in Brassica, making it difficult to compare and utilize these loci. European fodder turnip ECD04 is considered one of the most valuable resources for CR breeding. To explore the genetic and evolutionary basis of CR in ECD04, we sequenced the genome of ECD04 using de novo assembly and identified 978 candidate R genes. Subsequently, the 28 published CR loci were physically mapped to 15 loci in the ECD04 genome, including 62 candidate CR genes. Among them, two CR genes, CRA3.7.1 and CRA8.2.4, were functionally validated. Phylogenetic analysis revealed that CRA3.7.1 and CRA8.2.4 originated from a common ancestor before the whole-genome triplication (WGT) event. In clubroot susceptible Brassica species, CR-gene homologues were affected by transposable element (TE) insertion, resulting in the loss of CR function. It can be concluded that the current functional CR genes in Brassica rapa and non-functional CR genes in other Brassica species were derived from a common ancestral gene before WGT. Finally, a hypothesis for CR gene evolution is proposed for further discussion.

In silico integration of disease resistance QTL, genes and markers with the Brassica juncea physical map

Brassica juncea (AABB), Indian mustard, is a source of disease resistance genes for a wide range of pathogens. The availability of reference genome sequences for B. juncea has made it possible to characterise the genomic structure and distribution of these disease resistance genes. Potentially functional disease resistance genes can be identified by co-localization with genetically mapped disease resistance quantitative trait loci (QTL). Here we identify and characterise disease resistance gene analogs (RGAs), including nucleotide-binding site-leucine-rich repeat (NLR), receptor-like kinase (RLK) and receptor-like protein (RLP) classes, and investigate their association with disease resistance QTL intervals. The molecular genetic marker sequences for four white rust (Albugo candida) disease resistance QTL, six blackleg (Leptosphaeria maculans) disease resistance QTL and BjCHI1, a gene cloned from B. juncea for hypocotyl rot disease, were extracted from previously published studies and used to compare with candidate RGAs. Our results highlight the complications for the identification of functional resistance genes, including the duplicated appearance of genetic markers for several resistance loci, including Ac2(t), AcB1-A4.1, AcB1-A5.1, Rlm6 and PhR2 in both the A and B genomes, due to the presence of homoeologous regions. Furthermore, the white rust loci, Ac2(t) and AcB1-A4.1, mapped to the same position on chromosome A04 and may be different alleles of the same gene. Despite these challenges, a total of nine candidate genomic regions hosting 14 RLPs, 28 NLRs and 115 RLKs were identified. This study facilitates the mapping and cloning of functional resistance genes for applications in crop improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01309-5.

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.

Allelic variation of a clubroot resistance gene (Crr1a) in Japanese cultivars of Chinese cabbage (Brassica rapa L.)

Clubroot resistance (CR) is an important trait in Chinese cabbage breeding worldwide. Although Crr1a, the gene responsible for clubroot-resistance, has been cloned and shown to encode the NLR protein, its allelic variation and molecular function remain unknown. Here, we investigated the sequence variation and function of three Crr1a alleles cloned from six CR F₁ cultivars of Chinese cabbage. Gain-of-function analysis revealed that Crr1a^Kinami90_a isolated from the cv. 'Kinami 90' conferred clubroot resistance as observed for Crr1a^G004 . Because two susceptible alleles commonly lacked 172 amino acids in the C-terminal region, we investigated clubroot resistance in transgenic Arabidopsis harboring the chimeric Crr1a, in which 172 amino acids of the functional alleles were fused to the susceptible alleles. The fusion of the C-terminal region to the susceptible alleles restored resistance, indicating that their susceptibility was caused by the lack of the C-terminus. We developed DNA markers to detect the two functional Crr1a alleles, and demonstrated that the functional Crr1a alleles were frequently found in European fodder turnips, whereas they were rarely introduced into Japanese CR cultivars of Chinese cabbage. These results would contribute to CR breeding via marker-assisted selection and help our understanding of the molecular mechanisms underlying clubroot resistance.

Development of molecular markers based on CRa gene sequencing of different clubroot disease-resistant cultivars of Chinese cabbage

BACKGROUND

CRa is a key gene in Chinese cabbage (Brassica rapa ssp. pekinensis) that confers resistance to Plasmodiophora brassicae. In order to efficiently screen the clubroot resistance (CR) gene CRa in breeding, two functional codominant markers of the CRa gene were developed.

METHODS AND RESULTS

In this study, through comparing the CRa allele sequences in resistant and susceptible cultivars of Chinese cabbage, we found two insertion and deletion of sequence variations in the fourth exon between resistant and susceptible cultivars. Two functional codominant markers for CRa gene were obtained based on the variations, namely, CRaEX04-1 and CRaEX04-3. The lengths of the extended fragment of CRaEX04-1 marker were 321 bp and 186 bp in resistant and susceptible cultivars, respectively. In contrast, those of CRaEX04-3 were 704 bp and 413 bp, respectively. We verified the genetic stability between the developed markers and CRa gene using 57 Chinese cabbage cultivars with known resistance and two genetic populations. The results showed that the marker identification was completely consistent with the known phenotypes in 57 cultivars. The marker identification results followed the 3:1 of Mendel's first law in the F₂ population, and the 1:1 of Mendel's first law in the BC₁.

CONCLUSIONS

CRaEX04-1 and CRaEX04-3 can be used as a practical molecular marker for breeding and germplasm resource creation of clubroot disease-resistant Chinese cabbage.

Marker-Assisted Pyramiding of Genes for Multilocular Ovaries, Self-Compatibility, and Clubroot Resistance in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Molecular marker-assisted gene pyramiding combined with backcrossing has been widely applied for crop variety improvement. Molecular marker identification could be used in the early stage of breeding to achieve the rapid and effective pyramiding of multiple genes. To create high-quality germplasm for Chinese cabbage breeding, multi-gene pyramiding for self-compatibility, multilocular, and clubroot resistance was performed through molecular marker-assisted selection. The results showed that self-compatibility and multilocular traits were controlled by a pair of recessive genes. Two flanking markers, sau_um190 and cun_246a, and marker Teo-1, based on the gene sequence related to multilocular ovaries, were used for multilocular ovary trait selection. Two flanking markers, SCF-6 and SC-12, and marker Sal-SLGI /PK1+PK4, based on the gene sequence, were used for self-compatibility selection. Two flanking markers, TCR74 and TCR79, closely linked to clubroot resistance gene CRb, were used as foreground selection markers. Based on Chinese cabbage genomic information, 111 SSR markers covering 10 chromosomes were applied for background selection. After multiple generations of selection, a multi-gene pyramided line from a BC4F2 population with self-compatibility, multilocular ovaries, and clubroot resistance was obtained with a high genomic background recovery rate. The improved pyramided line is expected to be utilized as a potential material in further breeding programs.

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.

Identification and Fine-Mapping of Clubroot (Plasmodiophora brassicae) Resistant QTL in Brassica rapa

European fodder turnips (Brassica rapa ssp. rapifera) were identified as sources of clubroot resistance (CR) and have been widely used in Brassica resistance breeding. An F2 population derived from a cross between a resistant turnip and a susceptible Chinese cabbage was used to determine the inheritance and locating the resistance Quantitative Trait Loci (QTLs). The parents showed to be very resistant/susceptible to the field isolates (pathotype 4) of clubroot from Henan in China. After inoculation, 27 very resistant or susceptible individuals were selected to construct bulks, respectively. Next-generation-sequencing-based Bulk Segregant Analysis Sequencing (BSA-Seq) was used and located resistance QTL on chromosome A03 (3.3–7.5 Mb) and A08 (0.01–6.5 Mb), named Bcr1 and Bcr2, respectively. Furthermore, an F3 population including 180 families derived from F2 individuals was phenotyped and used to verify and narrow candidate regions. Ten and seven Kompetitive Allele-Specific PCR (KASP) markers narrowed the target regions to 4.3–4.78 Mb (A03) and 0.02–0.79 Mb (A08), respectively. The phenotypic variation explained (PVE) of the two QTLs were 33.3% and 13.3% respectively. The two candidate regions contained 99 and 109 genes. In the A03 candidate region, there were three candidate R genes, namely Bra006630, Bra006631 and Bra006632. In the A08 candidate region, there were two candidate R genes, namely Bra030815 and Bra030846.

Fine Mapping of Clubroot Resistance Loci CRA8.1 and Candidate Gene Analysis in Chinese Cabbage (Brassica rapa L.)

Clubroot is caused by Plasmodiophora brassicae, which threatens Brassicaceae crop production worldwide. In recent years, there has been an outbreak and rapid spread of clubroot in many major cruciferous crop-producing areas of China. In this study, we identified a cabbage material DingWen (DW) with different resistant capabilities from Huashuang5R (H5R) and Huayouza62R of Brassica napus, which are currently used as the main resistant cultivars for clubroot management in China. We used a next-generation sequencing-based bulked segregant analysis approach, combined with genetic mapping to identify clubroot-resistant (CR) genes from F₁ population generated from a cross between the DW (CR) and HZSX (clubroot susceptible). The CR locus of DW (named CRA8.1) was mapped to a region between markers A08-4346 and A08-4853, which contains two different loci CRA8.1a and CRA8.1b after fine mapping. The CRA8.1b loci contain a fragment of 395 kb between markers A08-4624 and A08-4853 on A08 chromosome, and it is responsible for the resistance to PbZj and PbXm isolates. However, together with CRA8.1a, corresponding to a 765-kb region between markers A08-4346 and A08-4624, then it can confer resistance to PbXm ⁺. Finally, through expression analysis between resistant and susceptible materials, two genes encoding TIR-NBS-LRR proteins (BraA08g039211E and BraA08g039212E) and one gene encoding an RLP protein (BraA08g039193E) were identified to be the most likely CR candidates for the peculiar resistance in DW.

A Proteome-Level Investigation Into Plasmodiophora brassicae Resistance in Brassica napus Canola

Clubroot of Brassicaceae, an economically important soil borne disease, is caused by Plasmodiophora brassicae Woronin, an obligate, biotrophic protist. This disease poses a serious threat to canola and related crops in Canada and around the globe causing significant losses. The pathogen is continuously evolving and new pathotypes are emerging, which necessitates the development of novel resistant canola cultivars to manage the disease. Proteins play a crucial role in many biological functions and the identification of differentially abundant proteins (DAP) using proteomics is a suitable approach to understand plant-pathogen interactions to assist in the development of gene specific markers for developing clubroot resistant (CR) cultivars. In this study, P. brassicae pathotype 3 (P3H) was used to challenge CR and clubroot susceptible (CS) canola lines. Root samples were collected at three distinct stages of pathogenesis, 7-, 14-, and 21-days post inoculation (DPI), protein samples were isolated, digested with trypsin and subjected to liquid chromatography with tandem mass spectrometry (LC-MS/MS) analysis. A total of 937 proteins demonstrated a significant (q-value < 0.05) change in abundance in at least in one of the time points when compared between control and inoculated CR-parent, CR-progeny, CS-parent, CS-progeny and 784 proteins were significantly (q < 0.05) changed in abundance in at least in one of the time points when compared between the inoculated- CR and CS root proteomes of parent and progeny across the three time points tested. Functional annotation of differentially abundant proteins (DAPs) revealed several proteins related to calcium dependent signaling pathways. In addition, proteins related to reactive oxygen species (ROS) biochemistry, dehydrins, lignin, thaumatin, and phytohormones were identified. Among the DAPs, 73 putative proteins orthologous to CR proteins and quantitative trait loci (QTL) associated with eight CR loci in different chromosomes including chromosomes A3 and A8 were identified. Proteins including BnaA02T0335400WE, BnaA03T0374600WE, BnaA03T0262200WE, and BnaA03T0464700WE are orthologous to identified CR loci with possible roles in mediating clubroot responses. In conclusion, these results have contributed to an improved understanding of the mechanisms involved in mediating response to P. brassicae in canola at the protein level.

Identification of Candidate Genes for Clubroot-Resistance in Brassica oleracea Using Quantitative Trait Loci-Sequencing

Clubroot caused by Plasmodiophora brassicae is a devastating disease of cabbage (Brassica oleracea). To identify quantitative trait loci (QTLs) for clubroot resistance (CR) in B. oleracea, genomic resequencing was carried out in two sets of extreme pools, group I and group II, which were constructed separately from 110 and 74 F2 cloned lines derived from the cross between clubroot-resistant (R) cabbage "GZ87" (against race 4) and susceptible (S) cabbage "263." Based on the QTL-sequencing (QTL-Seq) analysis of group I and group II, three QTLs (i.e., qCRc7-2, qCRc7-3, and qCRc7-4) were determined on the C07 chromosome. RNA-Seq and qRT-PCR were conducted in the extreme pools of group II before and after inoculation, and two potential candidate genes (i.e., Bol037115 and Bol042270), which exhibiting upregulation after inoculation in the R pool but downregulation in the S pool, were identified from the three QTLs on C07. A functional marker "SWU-OA" was developed from qCRc7-4 on C07, exhibiting ∼95% accuracy in identifying CR in 56 F2 lines. Our study will provide valuable information on resistance genes against P. brassicae and may accelerate the breeding process of B. oleracea with CR.

Fingerprint construction through genotyping by sequencing for applied breeding in Brassica rapa

This study evaluated genotyping by sequencing (GBS) protocol for fingerprinting Brassica rapa and the data derived were more reliable than the re-sequencing data of B. rapa. Of the 10 enzyme solutions used to analyze the numbers of genotypes and single nucleotide polymorphisms (SNPs) in B. rapa, five solutions showed better results, namely: A (HaeIII, 450-500 bp), E (RsaI+HaeIII, 500-550 bp), F (RsaI+HaeIII, 500-600 bp), G (RsaI+HaeIII, 'All' fragment), and J (RsaI+EcoRV-HF®, 'All' fragment). The five enzyme solutions showed less than 40% similarity in different individuals from various samples, and 90% similarity in between two individuals from one sample. The E enzyme solution was most suitable for fingerprinting B. rapa revealing well-distributed SNPs in the whole genome. Of the 82 highly inbred lines and 18 F1 lines of B. rapa sequenced by GBS in E enzyme solution, known parents of 10 F1 lines were verified and male parents were discovered for 8 F1 lines that had only known female parents. This study provided a valuable method for screening parents for F1 lines in B. rapa for applied breeding through efficient evaluation of GBS with varied library construction strategies.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Local Duplication of TIR-NBS-LRR Gene Marks Clubroot Resistance in Brassica napus cv. Tosca

Clubroot, caused by Plasmodiophora brassicae infection, is a disease of growing importance in cruciferous crops, including oilseed rape (Brassica napus). The affected plants exhibit prominent galling of the roots that impairs their capacity for water and nutrient uptake, which leads to growth retardation, wilting, premature ripening, or death. Due to the scarcity of effective means of protection against the pathogen, breeding of resistant varieties remains a crucial component of disease control measures. The key aspect of the breeding process is the identification of genetic factors associated with variable response to the pathogen exposure. Although numerous clubroot resistance loci have been described in Brassica crops, continuous updates on the sources of resistance are necessary. Many of the resistance genes are pathotype-specific, moreover, resistance breakdowns have been reported. In this study, we characterize the clubroot resistance locus in the winter oilseed rape cultivar "Tosca." In a series of greenhouse experiments, we evaluate the disease severity of P. brassicae-challenged "Tosca"-derived population of doubled haploids, which we genotype with Brassica 60 K array and a selection of SSR/SCAR markers. We then construct a genetic map and narrow down the resistance locus to the 0.4 cM fragment on the A03 chromosome, corresponding to the region previously described as Crr3. Using Oxford Nanopore long-read genome resequencing and RNA-seq we review the composition of the locus and describe a duplication of TIR-NBS-LRR gene. Further, we explore the transcriptomic differences of the local genes between the clubroot resistant and susceptible, inoculated and control DH lines. We conclude that the duplicated TNL gene is a promising candidate for the resistance factor. This study provides valuable resources for clubroot resistance breeding programs and lays a foundation for further functional studies on clubroot resistance.

Development of Ogura CMS restorers in Brassica oleracea subspecies via direct RfoB gene transformation

The Ogura CMS Rfo^B restorer developing via Rfo^B gene transformation was utilized to produce specific morphological Ogura CMS restorers and clubroot resistance lines in Brassica oleracea subspecies. Brassica oleracea vegetables including cabbage, cauliflower, kohlrabi, Brussels sprouts and Chinese kale are morphologically very different despite being members of the same species. The Ogura cytoplasmic male sterility (CMS) system is the most stable strategy for the hybrid breeding of these species. However, this limits the utilization of some excellent genes due to the lack of fertile restorer genes in the system. Herein, to efficaciously use Ogura CMS, the Ogura CMS Rfo^B restorer was produced by transforming the modified Rfo^B restorer gene into the Ogura CMS line 'CMS2016' of B. oleracea var. capitata. This gene was shown to recover fertility of natural Ogura CMS lines in B. oleracea subspecies and create transient Ogura CMS Rfo^B restorers such as the clubroot resistance Ogura CMS Rfo^B restorer. Interestingly, clubroot resistant individuals without transgenic elements were screened in the progenies of hybrids between B. oleracea inbred lines and the clubroot resistance Ogura CMS Rfo^B restorer. In addition, 18 different morphological Ogura CMS restorers were developed to specifically recover fertile of Ogura CMS cultivars in B. oleracea subspecies.

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.

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

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

Linkage mapping and QTL analysis of flowering time using ddRAD sequencing with genotype error correction in Brassica napus

BACKGROUND

Brassica napus is an important oilseed crop cultivated worldwide. During domestication and breeding of B. napus, flowering time has been a target of selection because of its substantial impact on yield. Here we use double digest restriction-site associated DNA sequencing (ddRAD) to investigate the genetic basis of flowering in B. napus. An F₂ mapping population was derived from a cross between an early-flowering spring type and a late-flowering winter type.

RESULTS

Flowering time in the mapping population differed by up to 25 days between individuals. High genotype error rates persisted after initial quality controls, as suggested by a genotype discordance of ~ 12% between biological sequencing replicates. After genotype error correction, a linkage map spanning 3981.31 cM and compromising 14,630 single nucleotide polymorphisms (SNPs) was constructed. A quantitative trait locus (QTL) on chromosome C2 was detected, covering eight flowering time genes including FLC.

CONCLUSIONS

These findings demonstrate the effectiveness of the ddRAD approach to sample the B. napus genome. Our results also suggest that ddRAD genotype error rates can be higher than expected in F₂ populations. Quality filtering and genotype correction and imputation can substantially reduce these error rates and allow effective linkage mapping and QTL analysis.

Comparative Transcriptome Analysis of Rutabaga (Brassica napus) Cultivars Indicates Activation of Salicylic Acid and Ethylene-Mediated Defenses in Response to Plasmodiophora brassicae

Clubroot, caused by Plasmodiophora brassicae Woronin, is an important soilborne disease of Brassica napus L. and other crucifers. To improve understanding of the mechanisms of resistance and pathogenesis in the clubroot pathosystem, the rutabaga (B. napus subsp. rapifera Metzg) cultivars ‘Wilhelmsburger’ (resistant) and ‘Laurentian’ (susceptible) were inoculated with P. brassicae pathotype 3A and their transcriptomes were analyzed at 7, 14, and 21 days after inoculation (dai) by RNA sequencing (RNA-seq). Thousands of transcripts with significant changes in expression were identified in each host at each time-point in inoculated vs. non-inoculated plants. Molecular responses at 7 and 14 dai supported clear differences in the clubroot response mechanisms of the two genotypes. Both the resistant and the susceptible cultivars activated receptor-like protein (RLP) genes, resistance (R) genes, and genes involved in salicylic acid (SA) signaling as clubroot defense mechanisms. In addition, genes related to calcium signaling and genes encoding leucine-rich repeat (LRR) receptor kinases, the respiratory burst oxidase homolog (RBOH) protein, and transcription factors such as WRKYs, ethylene responsive factors, and basic leucine zippers (bZIPs), appeared to be upregulated in ‘Wilhelmsburger’ to restrict P. brassicae development. Some of these genes are essential components of molecular defenses, including ethylene (ET) signaling and the oxidative burst. Our study highlights the importance of activation of genes associated with SA- and ET-mediated responses in the resistant cultivar. A set of candidate genes showing contrasting patterns of expression between the resistant and susceptible cultivars was identified and includes potential targets for further study and validation through approaches such as gene editing.

The Use of Genetic and Gene Technologies in Shaping Modern Rapeseed Cultivars (Brassica napus L.)

Since their domestication, Brassica oilseed species have undergone progressive transformation allied with the development of breeding and molecular technologies. The canola (Brassica napus) crop has rapidly expanded globally in the last 30 years with intensive innovations in canola varieties, providing for a wider range of markets apart from the food industry. The breeding efforts of B. napus, the main source of canola oil and canola meal, have been mainly focused on improving seed yield, oil quality, and meal quality along with disease resistance, abiotic stress tolerance, and herbicide resistance. The revolution in genetics and gene technologies, including genetic mapping, molecular markers, genomic tools, and gene technology, especially gene editing tools, has allowed an understanding of the complex genetic makeup and gene functions in the major bioprocesses of the Brassicales, especially Brassica oil crops. Here, we provide an overview on the contributions of these technologies in improving the major traits of B. napus and discuss their potential use to accomplish new improvement targets.

Identification of resistance loci in Chinese and Canadian canola/rapeseed varieties against Leptosphaeria maculans based on genome-wide association studies

Background

The fungal pathogen Leptosphaeria maculans (Lm). causes blackleg disease on canola/rapeseed in many parts of the world. It is important to use resistant cultivars to manage the disease and minimize yield losses. In this study, twenty-two Lm isolates were used to identify resistance genes in a collection of 243 canola/rapeseed (Brassica napus L.) accessions from Canada and China. These Lm isolates carry different compliments of avirulence genes, and the investigation was based on a genome-wide association study (GWAS) and genotype-by-sequencing (GBS).

Results

Using the CROP-SNP pipeline, a total of 81,471 variants, including 78,632 SNPs and 2839 InDels, were identified. The GWAS was performed using TASSEL 5.0 with GLM + Q model. Thirty-two and 13 SNPs were identified from the Canadian and Chinese accessions, respectively, tightly associated with blackleg resistance with P values < 1 × 10^− 4. These SNP loci were distributed on chromosomes A03, A05, A08, A09, C01, C04, C05, and C07, with the majority of them on A08 followed by A09 and A03. The significant SNPs identified on A08 were all located in a 2010-kb region and associated with resistance to 12 of the 22 Lm isolates. Furthermore, 25 resistance gene analogues (RGAs) were identified in these regions, including two nucleotide binding site (NBS) domain proteins, fourteen RLKs, three RLPs and six TM-CCs. These RGAs can be the potential candidate genes for blackleg resistance.

Conclusion

This study provides insights into potentially new genomic regions for discovery of additional blackleg resistance genes. The identified regions associated with blackleg resistance in the germplasm collection may also contribute directly to the development of canola varieties with novel resistance genes against blackleg of canola.

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.

Quantitative Trait Locus Mapping of Clubroot Resistance and Plasmodiophora brassicae Pathotype Banglim-Specific Marker Development in Brassica rapa

Clubroot resistance is an economically important trait in Brassicaceae crops. Although many quantitative trait loci (QTLs) for clubroot resistance have been identified in Brassica, disease-related damage continues to occur owing to differences in host variety and constant pathogen variation. Here, we investigated the inheritance of clubroot resistance in a double haploid population developed by crossing clubroot resistant and susceptible lines "09CR500" and "09CR501", respectively. The resistance of "09CR500" to Plasmodiophora brassicae pathotype "Banglim" was controlled as a single dominant gene, with the segregation of resistance and susceptibility being nearly 1:1. PbBrA08^Banglim was identified as having a logarithm of odds value of 7.9-74.8, and a phenotypic variance of 26.0-97.1% with flanking marker "09CR.11390652" in A08. After aligning QTL regions to the B. rapa reference genome, 11 genes were selected as candidates. PbBrA08^Banglim was located near Crr1, CRs, and Rcr9 loci, but differences were validated by marker analysis, gene structural variations, and gene expression levels, as well as phenotypic responses to the pathotype. Genotyping using the "09CR.11390652" marker accurately distinguished the Banglim-resistance phenotypes in the double haploid population. Thus, the developed marker will be useful in Brassica breeding programs, marker-assisted selection, and gene pyramiding to identify and develop resistant cultivars.

Genetics of Clubroot and Fusarium Wilt Disease Resistance in Brassica Vegetables: The Application of Marker Assisted Breeding for Disease Resistance

The genus Brassica contains important vegetable crops, which serve as a source of oil seed, condiments, and forages. However, their production is hampered by various diseases such as clubroot and Fusarium wilt, especially in Brassica vegetables. Soil-borne diseases are difficult to manage by traditional methods. Host resistance is an important tool for minimizing disease and many types of resistance (R) genes have been identified. More than 20 major clubroot (CR) disease-related loci have been identified in Brassica vegetables and several CR-resistant genes have been isolated by map-based cloning. Fusarium wilt resistant genes in Brassica vegetables have also been isolated. These isolated R genes encode the toll-interleukin-1 receptor/nucleotide-binding site/leucine-rice-repeat (TIR-NBS-LRR) protein. DNA markers that are linked with disease resistance allele have been successfully applied to improve disease resistance through marker-assisted selection (MAS). In this review, we focused on the recent status of identifying clubroot and Fusarium wilt R genes and the feasibility of using MAS for developing disease resistance cultivars in Brassica vegetables.

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.

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.

iTRAQ-based quantitative analysis reveals proteomic changes in Chinese cabbage (Brassica rapa L.) in response to Plasmodiophora brassicae infection

Clubroot disease is one of the major diseases affecting Brassica crops, especially Chinese cabbage (Brassica rapa L. ssp. pekinensis), which is known to be highly susceptible to the disease. In this study, the obligate biotrophic protist Plasmodiophora brassicae Woronin was used to infect the roots of Chinese cabbage seedlings. The disease symptoms were noticeable at 28 and 35 days after inoculation (DAI) in the susceptible (CM) line. Using isobaric tags for relative and absolute quantitation (iTRAQ) analysis, a total of 5,003 proteins of differential abundance were identified in the resistant/susceptible lines, which could be quantitated by dipeptide or polypeptide segments. Gene ontology (GO) analysis indicated that the differentially expressed proteins (DEPs) between the susceptible (CM) and resistant (CCR) lines were associated with the glutathione transferase activity pathway, which could catalyze the combination of glutathione and other electrophilic compounds to protect plants from disease. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the DEPs may be significantly enriched cytokinin signaling or arginine biosynthesis pathways, both of which are responses to stimuli and are plant defense reactions. The cytokinins may facilitate cell division in the shoot, resulting in the hypertrophy and formation of galls and the presentation of typical clubroot symptoms. In this study, the proteomic results provide a new perspective for creating germplasm resistance to P. brassicae, as well as a genetic basis for breeding to improve Chinese cabbage.

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.

Identifying candidate genes for Phytophthora capsici resistance in pepper (Capsicum annuum) via genotyping-by-sequencing-based QTL mapping and genome-wide association study

Phytophthora capsici (Leon.) is a globally prevalent, devastating oomycete pathogen that causes root rot in pepper (Capsicum annuum). Several studies have identified quantitative trait loci (QTL) underlying resistance to P. capsici root rot (PcRR). However, breeding for pepper cultivars resistant to PcRR remains challenging due to the complexity of PcRR resistance. Here, we combined traditional QTL mapping with GWAS to broaden our understanding of PcRR resistance in pepper. Three major-effect loci (5.1, 5.2, and 5.3) conferring broad-spectrum resistance to three isolates of P. capsici were mapped to pepper chromosome P5. In addition, QTLs with epistatic interactions and minor effects specific to isolate and environment were detected on other chromosomes. GWAS detected 117 significant SNPs across the genome associated with PcRR resistance, including SNPs on chromosomes P5, P7, and P11 that colocalized with the QTLs identified here and in previous studies. Clusters of candidate nucleotide-binding site-leucine-rich repeat (NBS-LRR) and receptor-like kinase (RLK) genes were predicted within the QTL and GWAS regions; such genes often function in disease resistance. These candidate genes lay the foundation for the molecular dissection of PcRR resistance. SNP markers associated with QTLs for PcRR resistance will be useful for marker-assisted breeding and genomic selection in pepper breeding.