Identification and characterization of Crr1a, a gene for resistance to clubroot disease (Plasmodiophora brassicae Woronin) in Brassica rapa L

BjuA03.BNT1 plays a positive role in resistance to clubroot disease in resynthesized Brassica juncea L

Clubroot has become a major obstacle in rapeseed production. Breeding varieties resistant to clubroot is the most effective method for disease management. However, the clubroot-resistant germplasm of rapeseed remains limited. To tackle this challenge, we synthesized the clubroot-resistant mustard, CT19, via distant hybridization, and subsequently an F2 segregating population was created by intercrossing CT19 with a clubroot-susceptible germplasm CS15. A major-effect clubroot resistance QTL qCRa3-1 on chromosome A03 was identified through QTL scanning. Transcriptome analyses of CT19 and CS15 revealed that the mechanisms conferring resistance to Plasmodiophora brassica likely involved the regulation of flavonoid metabolism, fatty acid metabolism, and sulfur metabolism. By combining the results from transcriptome, QTL mapping, and gene sequencing, a candidate gene BjuA03.BNT1, encoding NLR (nucleotide-binding domain leucine-rich repeat-containing receptors) protein, was obtained. Intriguingly, comparing with CT19, a base T insertion was discovered in the BjuA03.BNT1 gene's coding sequence in CS15, resulting an alteration within the LRR conserved domain. Overexpression of BjuA03.BNT1 from CT19 notably enhanced the resistance to clubroot in Arabidopsis. Our investigations revealed that BjuA03.BNT1 regulated the resistance to clubroot by modulating fatty acid synthesis and the structure of cell wall. These results are highly relevant for molecular breeding to improve clubroot resistance in rapeseed.

Comparative Analysis of Transcriptomes Reveals Pathways and Verifies Candidate Genes for Clubroot Resistance in Brassica oleracea

Clubroot, a soil-borne disease caused by Plasmodiophora brassicae, is one of the most destructive diseases of Brassica oleracea all over the world. However, the mechanism of clubroot resistance remains unclear. In this research, transcriptome sequencing was conducted on root samples from both resistant (R) and susceptible (S) B. oleracea plants infected by P. brassicae. Then the comparative analysis was carried out between the R and S samples at different time points during the infection stages to reveal clubroot resistance related pathways and candidate genes. Compared with 0 days after inoculation, a total of 4991 differential expressed genes were detected from the S pool, while only 2133 were found from the R pool. Gene function enrichment analysis found that the effector-triggered immunity played a major role in the R pool, while the pathogen-associated molecular pattern triggered immune response was stronger in the S pool. Simultaneously, candidate genes were identified through weighted gene co-expression network analysis, with Bol010786 (CNGC13) and Bol017921 (SD2-5) showing potential for conferring resistance to clubroot. The findings of this research provide valuable insights into the molecular mechanisms underlying clubroot resistance and present new avenues for further research aimed at enhancing the clubroot resistance of B. oleracea through breeding.

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.

Multiple transcription factors involved in the response of Chinese cabbage against Plasmodiophora brassicae

Clubroot disease, which is caused by the obligate biotrophic protist Plasmodiophora brassicae, leads to the formation of galls, commonly known as pathogen-induced tumors, on the roots of infected plants. The identification of crucial regulators of host tumor formation is essential to unravel the mechanisms underlying the proliferation and differentiation of P. brassicae within plant cells. To gain insight into this process, transcriptomic analysis was conducted to identify key genes associated with both primary and secondary infection of P. brassicae in Chinese cabbage. Our results demonstrate that the k-means clustering of subclass 1, which exhibited specific trends, was closely linked to the infection process of P. brassicae. Of the 1610 differentially expressed genes (DEGs) annotated in subclass 1, 782 were identified as transcription factors belonging to 49 transcription factor families, including bHLH, B3, NAC, MYB_related, WRKY, bZIP, C2H2, and ERF. In the primary infection, several genes, including the predicted Brassica rapa probable pectate lyase, RPM1-interacting protein 4-like, L-type lectin-domain-containing receptor kinase, G-type lectin S-receptor-like serine, B. rapa photosystem II 22 kDa protein, and MLP-like protein, showed significant upregulation. In the secondary infection stage, 45 of 50 overlapping DEGs were upregulated. These upregulated DEGs included the predicted B. rapa endoglucanase, long-chain acyl-CoA synthetase, WRKY transcription factor, NAC domain-containing protein, cell division control protein, auxin-induced protein, and protein variation in compound-triggered root growth response-like and xyloglucan glycosyltransferases. In both the primary and secondary infection stages, the DEGs were predicted to be Brassica rapa putative disease resistance proteins, L-type lectin domain-containing receptor kinases, ferredoxin-NADP reductases, 1-aminocyclopropane-1-carboxylate synthases, histone deacetylases, UDP-glycosyltransferases, putative glycerol-3-phosphate transporters, and chlorophyll a-binding proteins, which are closely associated with plant defense responses, biosynthetic processes, carbohydrate transport, and photosynthesis. This study revealed the pivotal role of transcription factors in the initiation of infection and establishment of intracellular parasitic relationships during the primary infection stage, as well as the proliferation and differentiation of the pathogen within the host cell during the secondary infection stage.

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.

Comparative transcriptome analysis of canola carrying a single vs stacked resistance genes against clubroot

Pyramiding resistance genes may expand the efficacy and scope of a canola variety against clubroot (Plasmodiophora brassicae), a serious threat to canola production in western Canada. However, the mechanism(s) of multigenic resistance, especially the potential interaction among clubroot resistance (CR) genes, are not well understood. In this study, transcriptome was compared over three canola (Brassica napus L.) inbred/hybrid lines carrying a single CR gene in chromosome A03 (CRaM, Line 16) or A08 (Crr1rutb, Line 20), and both genes (CRaM+Crr1rutb, Line 15) inoculated with a field population (L-G2) of P. brassicae pathotype X, a new variant found in western Canada recently. The line16 was susceptible, while lines 15 and 20 were partially resistant. Functional annotation identified differential expression of genes (DEGs) involved in biosynthetic processes responsive to stress and regulation of cellular process; The Venn diagram showed that the partially resistant lines 15 and 20 shared 1,896 differentially expressed genes relative to the susceptible line 16, and many of these DEGs are involved in defense responses, activation of innate immunity, hormone biosynthesis and programmed cell death. The transcription of genes involved in Pathogen-Associated Molecular Pattern (PAMP)-Triggered and Effector-Triggered Immunity (PTI and ETI) was particularly up-regulated, and the transcription level was higher in line 15 (CRaM + Crr1rutb) than in line 20 (Crr1rutb only) for most of the DEGs. These results indicated that the partial resistance to the pathotype X was likely conferred by the CR gene Crr1rutb for both lines 15 and 20 that functioned via the activation of both PTI and ETI signaling pathways. Additionally, these two CR genes might have synergistic effects against the pathotype X, based on the higher transcription levels of defense-related DEGs expressed by inoculated line 15, highlighting the benefit of gene stacking for improved canola resistance as opposed to a single CR gene alone.

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.

Pathotyping Systems and Pathotypes of Plasmodiophora brassicae-Navigating toward the Optimal Classification

Plasmodiophora brassicae Woronin, an obligate biotrophic soil-borne pathogen, poses a significant threat to cruciferous crops worldwide by causing the devastating disease known as clubroot. Pathogenic variability in P. brassicae populations has been recognized since the 1930s based on its interactions with Brassica species. Over time, numerous sets of differential hosts have been developed and used worldwide to explore the pathogenic variability within P. brassicae populations. These sets encompass a range of systems, including the Williams system, the European Clubroot Differential set (ECD), the Brassica napus set, the Japanese Clubroot Differential Set, the Canadian Clubroot Differential Set (CCS), the Korean Clubroot Differential Set, and the Chinese Sinitic Clubroot Differential set (SCD). However, all existing systems possess both advantages as well as limitations regarding the detection of pathotypes from various Brassica species and their corresponding virulence pattern on Brassica genotypes. This comprehensive review aims to compare the main differential systems utilized in classifying P. brassicae pathotypes worldwide. Their strengths, limitations, and implications are evaluated, thereby enhancing our understanding of pathogenic variability.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

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.

Map-based cloning and functional analysis of a major quantitative trait locus, BolC.Pb9.1, controlling clubroot resistance in a wild Brassica relative (Brassica macrocarpa)

KEY MESSAGE

A causal gene BoUGT76C2, conferring clubroot resistance in wild Brassica oleracea, was identified and functionally characterized. Clubroot is a devastating soil-borne disease caused by the obligate biotrophic pathogen Plasmodiophora brassica (P. brassicae), which poses a great threat to Brassica oleracea (B. oleracea) production. Although several QTLs associated with clubroot resistance (CR) have been mapped in cultivated B. oleracea, none have been cloned in B. oleracea. Previously, we found that the wild B. oleracea B2013 showed high resistance to clubroot. In this study, we constructed populations using B2013 and broccoli line 90196. CR in B2013 is quantitatively inherited, and a major QTL, BolC.Pb9.1, was identified on C09 using QTL-seq and linkage analysis. The BolC.Pb9.1 was finely mapped to a 56 kb genomic region using F2:3 populations. From the target region, the candidate BoUGT76C2 showed nucleotide variations between the parents, and was inducible in response to P. brassicae infection. We generated BoUGT76C2 overexpression lines in the 90196 background, which showed significantly enhanced resistance to P. brassicae compared to the WT line, suggesting that BoUGT76C2 corresponds to the resistance gene BolC.Pb.9.1. This is the first report on the CR gene map-based cloning and functional analysis from wild relatives, which provides a theoretical basis to the understanding of the molecular mechanism of CR, and lays a foundation to improve the CR of cultivated B. oleracea.

A Ca2+ sensor BraCBL1.2 involves in BraCRa-mediated clubroot resistance in Chinese cabbage

Clubroot disease caused by Plasmodiophora brassicae (P. brassicae) severely threatens the cultivation of Cruciferous plants, especially Chinese cabbage. Recently, resistance genes in plants have been reported to encode for a Ca2+-permeable channel in the plasma membrane, which can mediate the cytosolic Ca2+ increase in plant cells upon pathogen attack. However, the downstream Ca2+ sensor and decoder are still unknown. In this study, we identified the virulent and avirulent P. brassicae isolates (Pbs) of two near isogenic lines, CR 3-2 and CS 3-2, with CR 3-2 harboring clubroot resistant gene BraCRa. The transcriptomic analysis was then conducted with CR 3-2 after inoculating with virulent isolate PbE and avirulent isolate Pb4. From the differentially expressed genes of transcriptomic data, we identified a Ca2+-sensor encoding gene, BraCBL1.2, that was highly induced in CR 3-2 during infection by Pb4 but not by PbE. Moreover, GUS histochemical staining and subcellular localization analysis revealed that BraCBL1.2 was specifically expressed in the root hair cells of Arabidopsis and encoded a putative Ca2+ sensor localized in the plasma membrane. We also developed an assay to investigate the BraCRa-mediated hypersensitive response (HR) in tobacco leaves. The results suggest that BraCBL1.2 is involved in the BraCRa-mediated plant ETI immune response against P. brassicae. In addition, we verified that overexpression of BraCBL1.2 enhanced clubroot resistance in Arabidopsis. Collectively, our data identified the involvement of a Ca2+ sensor in BraCRa-mediated clubroot resistance in Chinese cabbage, providing a theoretical basis for further research on the resistance of Chinese cabbage to P. brassicae.

Natural variation in Arabidopsis responses to Plasmodiophora brassicae reveals an essential role for Resistance to Plasmodiophora brasssicae 1 (RPB1)

Despite the identification of clubroot resistance genes in various Brassica crops our understanding of the genetic basis of immunity to Plasmodiophora brassicae infection in the model plant Arabidopsis thaliana remains limited. To address this issue, we performed a screen of 142 natural accessions and identified 11 clubroot-resistant Arabidopsis lines. Genome-wide association analysis identified several genetic loci significantly linked with resistance. Three genes from two of these loci were targeted for deletion by CRISPR/Cas9 mutation in resistant accessions Est-1 and Uod-1. Deletion of Resistance to Plasmodiophora brassicae 1 (RPB1) rendered both lines susceptible to the P. brassicae pathotype P1+. Further analysis of rpb1 knock-out Est-1 and Uod-1 lines showed that the RPB1 protein is required for activation of downstream defence responses, such as the expression of phytoalexin biosynthesis gene CYP71A13. RPB1 has recently been shown to encode a cation channel localised in the endoplasmic reticulum. The clubroot susceptible Arabidopsis accession Col-0 lacks a functional RPB1 gene; when Col-0 is transformed with RPB1 expression driven by its native promoter it is capable of activating RPB1 transcription in response to infection, but this is not sufficient to confer resistance. Transient expression of RPB1 in Nicotiana tabacum induced programmed cell death in leaves. We conclude that RPB1 is a critical component of the defence response to P. brassicae infection in Arabidopsis, acting downstream of pathogen recognition but required for the elaboration of effective resistance.

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.

WeiTsing, a pericycle-expressed ion channel, safeguards the stele to confer clubroot resistance

Plant roots encounter numerous pathogenic microbes that often cause devastating diseases. One such pathogen, Plasmodiophora brassicae (Pb), causes clubroot disease and severe yield losses on cruciferous crops worldwide. Here, we report the isolation and characterization of WeiTsing (WTS), a broad-spectrum clubroot resistance gene from Arabidopsis. WTS is transcriptionally activated in the pericycle upon Pb infection to prevent pathogen colonization in the stele. Brassica napus carrying the WTS transgene displayed strong resistance to Pb. WTS encodes a small protein localized in the endoplasmic reticulum (ER), and its expression in plants induces immune responses. The cryoelectron microscopy (cryo-EM) structure of WTS revealed a previously unknown pentameric architecture with a central pore. Electrophysiology analyses demonstrated that WTS is a calcium-permeable cation-selective channel. Structure-guided mutagenesis indicated that channel activity is strictly required for triggering defenses. The findings uncover an ion channel analogous to resistosomes that triggers immune signaling in the pericycle.

Key Advances in the New Era of Genomics-Assisted Disease Resistance Improvement of Brassica Species

Disease resistance improvement remains a major focus in breeding programs as diseases continue to devastate Brassica production systems due to intensive cultivation and climate change. Genomics has paved the way to understand the complex genomes of Brassicas, which has been pivotal in the dissection of the genetic underpinnings of agronomic traits driving the development of superior cultivars. The new era of genomics-assisted disease resistance breeding has been marked by the development of high-quality genome references, accelerating the identification of disease resistance genes controlling both qualitative (major) gene and quantitative resistance. This facilitates the development of molecular markers for marker assisted selection and enables genome editing approaches for targeted gene manipulation to enhance the genetic value of disease resistance traits. This review summarizes the key advances in the development of genomic resources for Brassica species, focusing on improved genome references, based on long-read sequencing technologies and pangenome assemblies. This is further supported by the advances in pathogen genomics, which have resulted in the discovery of pathogenicity factors, complementing the mining of disease resistance genes in the host. Recognizing the co-evolutionary arms race between the host and pathogen, it is critical to identify novel resistance genes using crop wild relatives and synthetic cultivars or through genetic manipulation via genome-editing to sustain the development of superior cultivars. Integrating these key advances with new breeding techniques and improved phenotyping using advanced data analysis platforms will make disease resistance improvement in Brassica species more efficient and responsive to current and future demands.

Understanding the Proteomes of Plant Development and Stress Responses in Brassica Crops

Brassica crops have great economic value due to their rich nutritional content and are therefore grown worldwide as oilseeds, vegetables, and condiments. Deciphering the molecular mechanisms associated with the advantageous phenotype is the major objective of various Brassica improvement programs. As large technological advancements have been achieved in the past decade, the methods to understand molecular mechanisms underlying the traits of interest have also taken a sharp upturn in plant breeding practices. Proteomics has emerged as one of the preferred choices nowadays along with genomics and other molecular approaches, as proteins are the ultimate effector molecules responsible for phenotypic changes in living systems, and allow plants to resist variable environmental stresses. In the last two decades, rapid progress has been made in the field of proteomics research in Brassica crops, but a comprehensive review that collates the different studies is lacking. This review provides an inclusive summary of different proteomic studies undertaken in Brassica crops for cytoplasmic male sterility, oil content, and proteomics of floral organs and seeds, under different biotic and abiotic stresses including post-translational modifications of proteins. This comprehensive review will help in understanding the role of different proteins in controlling plant phenotypes, and provides information for initiating future studies on Brassica breeding and improvement programs.

Canola with Stacked Genes Shows Moderate Resistance and Resilience against a Field Population of Plasmodiophora brassicae (Clubroot) Pathotype X

Genetic resistance is a cornerstone for managing clubroot (Plasmodiophora brassicae). However, when used repeatedly, a clubroot resistance (CR) gene can be broken rapidly. In this study, canola inbred/hybrid lines carrying one or two CR genes (Rcr1/CRa^M and Crr1^rutb) were assessed against P. brassicae pathotype X by repeated exposure to the same inoculum source under a controlled environment. Lines carrying two CR genes, either Rcr1 + Crr1^rutb or CRa^M + Crr1^rutb, showed partial resistance. Selected lines were inoculated with a field pathotype X population (L-G3) at 5 × 10⁶ resting spores/g soil, and all clubs were returned to the soil they came from six weeks after inoculation. The planting was repeated for five cycles, with diseased roots being returned to the soil after each cycle. The soil inoculum was quantified using qPCR before each planting cycle. All lines with a single CR gene were consistently susceptible, maintaining high soil inoculum levels over time. The lines carrying two CR genes showed much lower clubroot severity, resulting in a 10-fold decline in soil inoculum. These results showed that the CR-gene stacking provided moderate resistance against P. brassicae pathotype X, which may also help reduce the pathogen inoculum buildup in soil.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

The clubroot pathogen Plasmodiophora brassicae: A profile update

BACKGROUND

Plasmodiophora brassicae is the causal agent of clubroot disease of cruciferous plants and one of the biggest threats to the rapeseed (Brassica napus) and brassica vegetable industry worldwide.

DISEASE SYMPTOMS

In the advanced stages of clubroot disease wilting, stunting, yellowing, and redness are visible in the shoots. However, the typical symptoms of the disease are the presence of club-shaped galls in the roots of susceptible hosts that block the absorption of water and nutrients.

HOST RANGE

Members of the family Brassicaceae are the primary host of the pathogen, although some members of the family, such as Bunias orientalis, Coronopus squamatus, and Raphanus sativus, have been identified as being consistently resistant to P. brassicae isolates with variable virulence profile.

TAXONOMY

Class: Phytomyxea; Order: Plasmodiophorales; Family: Plasmodiophoraceae; Genus: Plasmodiophora; Species: Plasmodiophora brassicae (Woronin, 1877).

DISTRIBUTION

Clubroot disease is spread worldwide, with reports from all continents except Antarctica. To date, clubroot disease has been reported in more than 80 countries.

PATHOTYPING

Based on its virulence on different hosts, P. brassicae is classified into pathotypes or races. Five main pathotyping systems have been developed to understand the relationship between P. brassicae and its hosts. Nowadays, the Canadian clubroot differential is extensively used in Canada and has so far identified 36 different pathotypes based on the response of a set of 13 hosts.

EFFECTORS AND RESISTANCE

After the identification and characterization of the clubroot pathogen SABATH-type methyltransferase PbBSMT, several other effectors have been characterized. However, no avirulence gene is known, hindering the functional characterization of the five intercellular nucleotide-binding (NB) site leucine-rich-repeat (LRR) receptors (NLRs) clubroot resistance genes validated to date.

IMPORTANT LINK

Canola Council of Canada is constantly updating information about clubroot and P. brassicae as part of their Canola Encyclopedia: https://www.canolacouncil.org/canola-encyclopedia/diseases/clubroot/.

PHYTOSANITARY CATEGORIZATION

PLADBR: EPPO A2 list; Annex designation 9E.

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.

Fine mapping and candidate gene analysis of CRA3.7 conferring clubroot resistance in Brassica rapa

In this study, we fine-mapped a clubroot resistance gene CRA3.7 in Chinese cabbage and developed its closely linked marker syau-InDel3008 for marker-assisted selection in CR cultivars breeding. Chinese cabbage is an important leafy vegetable rich in many nutrients widely grown in China. Clubroot disease caused by an obligate biotrophic pathogen Plasmodiophora brassicae was rapidly spread and challenged to Chinese cabbage production. A clubroot resistance (CR) gene, CRA3.7, was mapped on chromosome A03 of Brassica rapa. A Chinese cabbage line 'CR510', which harbor homozygous resistance locus CRA3.7 was selected from a BC4F3 family. 'CR510' was crossed with a clubroot susceptible Chinese cabbage inbred line '59-1'. Total 51 recombinant plants were identified from an F₂ population including 3000 individuals. These recombinants were selfed and the clubroot resistance of F_2/3 families was evaluated. Finally, a clubroot resistance gene CRA3.7 was fine-mapped to an interval of approximately 386 kb between marker syau-InDel3024 and syau-InDel3008. According to the reference genome, total 54 genes including five encoding the TIR-NBS-LRR proteins was annotated in the fine-mapped region. Further, nine candidate's gene expression in parental lines at 7, 14 and 21 days after inoculation of P. brassicae were evaluated. Bra019376, Bra019401, Bra019403 and Bra019410 are highly expressed in 'CR510' than '59-1'. Gene sequence of Bra019410 from 'CR510' was cloned and identified different from CRa. Therefore, Bra019376, Bra019401, Bra019403 and Bra019410 are the most likely candidates for CRA3.7. Our research provides a valuable germplasm resource against P. brassicae Pb3 and CRA3.7 closely linked marker for marker-assisted selection in CR cultivars breeding.

Identification of a genomic region containing genes involved in resistance to four pathotypes of Plasmodiophora brassicae in Brassica rapa turnip ECD02

Clubroot, caused by Plasmodiophora brassicae, is an important disease of brassica crops worldwide. Vegetable turnip (Brassica rapa L.) have proven to be a source of clubroot resistance genes effective against many pathotypes of P. brassicae. The F₁ progeny from the cross B. rapa canola ACDC (susceptible, S) × B. rapa turnip ECD02 (resistant, R) were backcrossed with ACDC, then self-pollinated to produce BC₁ S₁ lines. All the F₁ plants were resistant to four pathotypes (3A, 3D, 3H, and 5X) of P. brassicae. Segregation for R and S in BC₁ to each pathotype was 1:1 and resistance reactions were highly correlated. From whole genome sequencing, 192.1 M sequences with 96% template coverage from ECD02, and 478.9 M sequences with 92% coverage from ACDC, were aligned with the reference genome of B. rapa. Genotyping-by-sequencing was performed on the BC₁ population. The number of aligned short reads per plant in the BC₁ ranged from 1.4 to 8.5 M sequences with 4-8% template coverage. We obtained 1,344 high-quality single-nucleotide polymorphism (SNP) loci with a mean missing rate at 0.27% and distributed them on 10 chromosomes. A single co-localized quantitative trait loci (QTL), designated as Rcr9^ECD02 on chromosome A08, conferred resistance to the four pathotypes. The QTL explained 68.9-74.4% of phenotypic variation with the logarithm of the odds values of 24.3 to 31.1. Bulked segregant analysis was performed, and 14 SNP markers linked to the gene were developed using the Kompetitive Allele Specific PCR. Rcr9^ECD02 was mapped into an interval of 2.2 cM, flanked by CF_A08_10664692 and CF_A08_12230973, which spanned 1.51 Mb on the chromosome and included 219 B. rapa genes. Four of these genes (BraA08g012910.3C, BraA08g012920.3C, BraA08g013130.3C, and BraA08g013630.3C) encoded disease resistance proteins.

Predicting Cloned Disease Resistance Gene Homologs (CDRHs) in Radish, Underutilised Oilseeds, and Wild Brassicaceae Species

Brassicaceae crops, including Brassica, Camelina and Raphanus species, are among the most economically important crops globally; however, their production is affected by several diseases. To predict cloned disease resistance (R) gene homologs (CDRHs), we used the protein sequences of 49 cloned R genes against fungal and bacterial diseases in Brassicaceae species. In this study, using 20 Brassicaceae genomes (17 wild and 3 domesticated species), 3172 resistance gene analogs (RGAs) (2062 nucleotide binding-site leucine-rich repeats (NLRs), 497 receptor-like protein kinases (RLKs) and 613 receptor-like proteins (RLPs)) were identified. CDRH clusters were also observed in Arabis alpina, Camelina sativa and Cardamine hirsuta with assigned chromosomes, consisting of 62 homogeneous (38 NLR, 17 RLK and 7 RLP clusters) and 10 heterogeneous RGA clusters. This study highlights the prevalence of CDRHs in the wild relatives of the Brassicaceae family, which may lay the foundation for rapid identification of functional genes and genomics-assisted breeding to develop improved disease-resistant Brassicaceae crop cultivars.

Genome-wide identification of biotin carboxyl carrier subunits of acetyl-CoA carboxylase in Brassica and their role in stress tolerance in oilseed Brassica napus

Background

Biotin carboxyl carrier protein (BCCP) is a subunit of Acetyl CoA-carboxylase (ACCase) which catalyzes the conversion of acetyl-CoA to malonyl-CoA in a committed step during the de novo biosynthesis of fatty acids. Lipids, lipid metabolites, lipid-metabolizing and -modifying enzymes are known to play a role in biotic and abiotic stress tolerance in plants. In this regard, an understanding of the Brassica napus BCCP genes will aid in the improvement of biotic and abiotic stress tolerance in canola.

Results

In this study, we identified 43 BCCP genes in five Brassica species based on published genome data. Among them, Brassica rapa, Brassica oleracea, Brassica nigra, Brassica napus and Brassica juncea had six, seven, seven, 10 and 13 BCCP homologs, respectively. Phylogenetic analysis categorized them into five classes, each with unique conserved domains. The promoter regions of all BCCP genes contained stress-related cis-acting elements as determined by cis-element analysis. We identified four and three duplicated gene pairs (segmental) in B. napus and B. juncea respectively, indicating the role of segmental duplication in the expansion of this gene family. The Ka/Ks ratios of orthologous gene pairs between Arabidopsis thaliana and five Brassica species were mostly less than 1.0, implying that purifying selection, i.e., selective removal of deleterious alleles, played a role during the evolution of Brassica genomes. Analysis of 10 BnaBCCP genes using qRT-PCR showed a different pattern of expression because of exposure of the plants to biotic stresses, such as clubroot and sclerotinia diseases, and abiotic stresses such as drought, low temperature and salinity stresses.

Conclusions

The identification and functional analysis of the Brassica BCCPs demonstrated that some of these genes might play important roles in biotic and abiotic stress responses. Results from this study could lay the foundation for a better understanding of these genes for the improvement of Brassica crops for stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08920-y.

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.

Genetic and molecular analysis reveals that two major loci and their interaction confer clubroot resistance in canola introgressed from rutabaga

Clubroot disease caused by Plasmodiophora brassicae is one of the serious threats to canola (Brassica napus L. subsp. napus) production. The evolution of new pathotypes rendering available resistances ineffective compel the introgression of new resistance into canola and extend our understanding of the genetic and molecular basis of the resistance. In this paper, we report the genetic and molecular basis of clubroot resistance in canola, introgressed from a rutabaga (B. napus L. subsp. rapifera Metzg. 'Polycross'), by using a doubled-haploid (DH) mapping population. Whole-genome resequencing (WGRS)-based bulked segregant analysis followed by genetic mapping and expression analysis of the genes in resistant and susceptible DH lines at 7 and 14 d after inoculation were carried out. Following this approach, two major quantitative trait loci (QTL) located at 14.41-15.44 Mb of A03 and at 9.96-11.09 Mb of A08 chromosomes and their interaction was observed to confer resistance to pathotypes 3H, 3A, and 3D. Analysis of the genes from the two QTL regions suggested that decreased expression of sugar transporter genes (BnaA03g29290D and BnaA03g29310D) may play an important role in resistance conferred by the A03 QTL, while increased expression of the toll/interleukin-1 receptor (TIR)-nucleotide binding (NB)-leucine rich repeat (LRR) (TNL) genes (BnaA08g10100D, BnaA08g09220D, and BnaA08g10540D) could be the major determinant of the resistance conferred by the A08 QTL. Single-nucleotide polymorphism (SNP) allele-specific polymerase chain reaction (PCR)-based markers, which could be detected by agarose gel electrophoresis, were also developed from the two QTL regions for use in breeding including pyramiding of multiple clubroot resistance genes.

Introgression of clubroot resistant gene into L. from based on homoeologous exchange

Clubroot is a soil-borne disease in cabbage (Brassica oleracea L. var. capitata L.) caused by Plasmodiophora brassicae, which poses a great threat to cabbage production. However, clubroot resistance (CR) genes in Brassica rapa could be introduced into the cabbage via breeding to make it clubroot resistant. In this study, CR genes from B. rapa were introduced into the cabbage genome and the mechanism of gene introgression was explored. Two methods were used to create CR materials: (i) The fertility of CR Ogura CMS cabbage germplasms containing CRa was restored by using an Ogura CMS restorer. After cytoplasmic replacement and microspore culture, CRa-positive microspore individuals were obtained. (ii) Distant hybridization was performed between cabbage and B. rapa, which contained three CR genes (CRa, CRb, and Pb8.1). Finally, BC2 individuals containing all three CR genes were obtained. Inoculation results showed that both CRa-positive microspore individuals and BC2 individuals containing three CR genes were resistant to race 4 of P. brassicae. Sequencing results from CRa-positive microspore individuals with specific molecular markers and genome-wide association study (GWAS) showed penetration at the homologous position of the cabbage genome by a 3.42 Mb CRa containing a fragment from B. rapa; indicating homoeologous exchange (HE) as the theoretical basis for the introgression of CR resistance. The successful introduction of CR into the cabbage genome in the present study can provide useful clues for creating introgression lines within other species of interest.

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.

Efficient marker-assisted breeding for clubroot resistance in elite Pol-CMS rapeseed varieties by updating the PbBa8.1 locus

Clubroot disease poses a severe threat to rapeseed (Brassica napus) production worldwide and has recently been spreading across China at an unprecedented pace. Breeding and cultivation of resistant varieties constitute a promising and environment-friendly approach to mitigating this threat. In this study, the clubroot resistance locus PbBa8.1 was successfully transferred into SC4, a shared paternal line of three elite varieties in five generations by marker-assisted backcross breeding. Kompetitive allele specific PCR (KASP) markers of clubroot resistance gene PbBa8.1 and its linked high erucic acid gene (FAE1) were designed and applied for foreground selection, and 1,000 single-nucleotide polymorphisms (SNPs) were selected and used for the background selection. This breeding strategy produced recombinants with the highest recovery ratio of the recurrent parent genome (> 95%) at BC2F2 while breaking the linkage with FAE1 during the selection. An updated version of the paternal line (SC4R) was generated at BC2F3, showing significantly improved clubroot resistance at the seedling stage via artificial inoculation, and was comparable to that of the donor parent. Field trials of the three elite varieties and their updated versions in five environments indicated similar agronomic appearance and final yield. The introduced breeding strategy precisely pyramids the PbBa8.1 and FAE1 loci with the assistance of technical markers in a shorter period and could be applied to other desirable traits for directional improvement in the future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01305-9.

Transferring of clubroot-resistant locus CRd from Chinese cabbage (Brassica rapa) to canola (Brassica napus) through interspecific hybridization

Clubroot, caused by Plasmodiophora brassicae is one of the most severe threats to brassica species in China and worldwide. Breeding for clubroot resistant varieties is one of the best ways to overcome this disease. In this study, we introduced clubroot resistance (CR) gene CRd from Chinese cabbage (85-74) into elite Brassica napus inbred line Zhongshuang 11 through interspecific hybridization and subsequent backcrossing with whole-genome molecular marker-assisted selection (MAS). The resistant test of CRd to P. brassicae isolates was evaluated in the greenhouse as well as in field conditions. Close linkage markers and the whole-chromosome background marker selection approach improved the recovery rate from 78.3% in BC₁ to 100% in BC₃F₁. The improved clubroot-resistant variety, Zhongshuang11R, was successfully selected in the BC₃F₂ generation. The greenhouse and field resistant tests revealed that Zhongshuang11R was resistant to P. brassicae pathotypes. The agronomic characteristics of Zhongshuang11R were similar to those of its recurrent parental line, including oil content, composition of fatty acid, plant height, primary effective branches, grain yield per plant and thousand-seed weight. In addition, the oil quality could satisfy the quality requirements for commercial rapeseed oil. Our results will enrich the resistant resources of canola and will certainly accelerate clubroot resistance breeding programs in B. napus.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

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.

Histopathology of the Plasmodiophora brassicae-Chinese Cabbage Interaction in Hosts Carrying Different Sources of Resistance

Clubroot is a serious soil-borne disease of crucifers caused by the obligate parasite Plasmodiophora brassicae. The genetic basis and histopathology of clubroot resistance in two Chinese cabbage (Brassica rapa ssp. pekinensis) inbred lines Bap055 and Bap246, challenged with pathotype 4 of P. brassicae, was evaluated. The Chinese cabbage cultivar "Juxin" served as a susceptible check. The resistance in Bap055 was found to be controlled by the CRa gene, while resistance in Bap246 fit a model of control by unknown recessive gene. Infection of the roots by P. brassicae was examined by inverted microscopy. Despite their resistance, primary and secondary infection were observed to occur in Bap055 and Bap246. Primary infection was detected at 2 days post-inoculation (DPI) in "Juxin," at 4 DPI in Bap055, and at 6 DPI in Bap246. Infection occurred most quickly on "Juxin," with 60% of the root hairs infected at 10 DPI, followed by Bap055 (31% of the root hairs infected at 12 DPI) and Bap246 (20% of the root hairs infected at 14 DPI). Secondary infection of "Juxin" was first observed at 8 DPI, while in Bap055 and Bap246, secondary infection was first observed at 10 DPI. At 14 DPI, the percentage of cortical infection in "Juxin," Bap055 and Bap246 was 93.3, 20.0, and 11.1%, respectively. Although cortical infection was more widespread in Bap055 than in Bap246, secondary infection in both of these hosts was restricted relative to the susceptible check, and the vascular system remained intact. A large number of binucleate secondary plasmodia were observed in "Juxin" and the vascular system was disrupted at 16 DPI; in Bap055 and Bap246, only a few secondary plasmodia were visible, with no binucleate secondary plasmodia. The defense mechanisms and expression of resistance appears to differ between Chinese cabbage cultivars carrying different sources of resistance.

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.

Development of novel clubroot resistant rapeseed lines (Brassica napus L.) effective against Japanese field isolates by marker assisted selection

Clubroot is an important disease infectible to cruciferous plants and a major threat to rapeseed production in Japan. However, no clubroot resistant rapeseed cultivars have been released. We surveyed pathotype variation of six isolates collected from rapeseed fields and found they were classified as pathotype groups 2 and 4 using Japanese F₁ Chinese cabbage cultivars. We produced the resynthesized clubroot resistant Brassica napus harboring two resistant loci, Crr1 and Crr2, by interspecific crossing and developed resistant rapeseed lines for southern and northern regions by marker-assisted selection and backcrossing. We improved the DNA marker for erucic acid content to remove linkage drag between Crr1 and high erucic acid content and successfully selected lines with clubroot resistance and zero erucic acid for northern regions. A novel line, 'Tohoku No. 106', suitable for southern regions showed stable resistance against all six isolates and high performance in infested fields. We conclude that Crr1 and Crr2 are important genes for CR rapeseed breeding and marker-assisted selection is effective in improving clubroot resistance.

Identification of Micro Ribonucleic Acids and Their Targets in Response to Plasmodiophora brassicae Infection in Brassica napus

Clubroot disease, which is caused by the soil-borne pathogen Plasmodiophora brassicae War (P. brassicae), is one of the oldest and most destructive diseases of Brassica and cruciferous crops in the world. Plant microRNAs [micro ribonucleic acids (miRNAs)] play important regulatory roles in several developmental processes. Although the role of plant miRNAs in plant-microbe interaction has been extensively studied, there are only few reports on the specific functions of miRNAs in response to P. brassicae. This study investigated the roles of miRNAs and their targets during P. brassicae infection in a pair of Brassica napus near-isogenic lines (NILs), namely clubroot-resistant line 409R and clubroot-susceptible line 409S. Small RNA sequencing (sRNA-seq) and degradome-seq were performed on root samples of 409R and 409S with or without P. brassicae inoculation. sRNA-seq identified a total of 48 conserved and 72 novel miRNAs, among which 18 had a significant differential expression in the root of 409R, while only one miRNA was differentially expressed in the root of 409S after P. brassicae inoculation. The degradome-seq analysis identified 938 miRNA target transcripts, which are transcription factors, enzymes, and proteins involved in multiple biological processes and most significantly enriched in the plant hormone signal transduction pathway. Between 409R and 409S, we found eight different degradation pathways in response to P. brassicae infection, such as those related to fatty acids. By combining published transcriptome data, we identified a total of six antagonistic miRNA-target pairs in 409R that are responsive to P. brassicae infection and involved in pathways associated with root development, hypersensitive cell death, and chloroplast metabolic synthesis. Our results reveal that P. brassicae infection leads to great changes in miRNA pool and target transcripts. More interestingly, these changes are different between 409R and 409S. Clarification of the crosstalk between miRNAs and their targets may shed new light on the possible mechanisms underlying the pathogen resistance against P. brassicae.

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.

Transcriptome and Coexpression Network Analyses Reveal Hub Genes in Chinese Cabbage (Brassica rapa L. ssp. pekinensis) During Different Stages of Plasmodiophora brassicae Infection

Clubroot, caused by the soil-borne protist Plasmodiophora brassicae, is one of the most destructive diseases of Chinese cabbage worldwide. However, the clubroot resistance mechanisms remain unclear. In this study, in both clubroot-resistant (DH40R) and clubroot-susceptible (DH199S) Chinese cabbage lines, the primary (root hair infection) and secondary (cortical infection) infection stages started 2 and 5 days after inoculation (dai), respectively. With the extension of the infection time, cortical infection was blocked and complete P. brassica resistance was observed in DH40R, while disease scales of 1, 2, and 3 were observed at 8, 13, and 22 dai in DH199S. Transcriptome analysis at 0, 2, 5, 8, 13, and 22 dai identified 5,750 relative DEGs (rDEGs) between DH40R and DH199S. The results indicated that genes associated with auxin, PR, disease resistance proteins, oxidative stress, and WRKY and MYB transcription factors were involved in clubroot resistance regulation. In addition, weighted gene coexpression network analysis (WGCNA) identified three of the modules whose functions were highly associated with clubroot-resistant, including ten hub genes related to clubroot resistance (ARF2, EDR1, LOX4, NHL3, NHL13, NAC29, two AOP1, EARLI 1, and POD56). These results provide valuable information for better understanding the molecular regulatory mechanism of Chinese cabbage clubroot resistance.

Influence of Genotype on High Glucosinolate Synthesis Lines of Brassica rapa

This study was conducted to investigate doubled haploid (DH) lines produced between high GSL (HGSL) Brassica rapa ssp. trilocularis (yellow sarson) and low GSL (LGSL) B. rapa ssp. chinensis (pak choi) parents. In total, 161 DH lines were generated. GSL content of HGSL DH lines ranged from 44.12 to 57.04 μmol·g⁻¹·dry weight (dw), which is within the level of high GSL B. rapa ssp. trilocularis (47.46 to 59.56 μmol g⁻¹ dw). We resequenced five of the HGSL DH lines and three of the LGSL DH lines. Recombination blocks were formed between the parental and DH lines with 108,328 single-nucleotide polymorphisms in all chromosomes. In the measured GSL, gluconapin occurred as the major substrate in HGSL DH lines. Among the HGSL DH lines, BrYSP_DH005 had glucoraphanin levels approximately 12-fold higher than those of the HGSL mother plant. The hydrolysis capacity of GSL was analyzed in HGSL DH lines with a Korean pak choi cultivar as a control. Bioactive compounds, such as 3-butenyl isothiocyanate, 4-pentenyl isothiocyanate, 2-phenethyl isothiocyanate, and sulforaphane, were present in the HGSL DH lines at 3-fold to 6.3-fold higher levels compared to the commercial cultivar. The selected HGSL DH lines, resequencing data, and SNP identification were utilized for genome-assisted selection to develop elite GSL-enriched cultivars and the industrial production of potential anti-cancerous metabolites such as gluconapin and glucoraphanin.

Discovery of BrATG6 and its potential role in Brassica rapa L. resistance to infection by Plasmodiophora brassicae

Clubroot disease, caused by Plasmodiophora brassicae infection, occurs in cruciferous vegetable crops in many areas of the world, sometimes leading to yield loss. In this study, a differentially expressed protein (0305), was found between control and clubroot-diseased Chinese cabbage (Brassica rapa L.) roots through two-dimensional electrophoresis. Mass spectrometry analysis showed that Bra003466 was highly matched to protein 0305. Because the sequence of Bra003466 had 89% percent identity with ATG6 of Arabidopsis thaliana and other Brassica, the gene was named as BrATG6. However, 790 bp sequences were mismatched with the cDNA sequence of the Bra003466 gene from the Brassica database. In this study, we cloned the cDNA of Bra003466 and found the BrATG6 was highly expressed in roots among all organs. When plants were inoculated with P. brassicae Woronin, the expression of BrATG6 was significantly increased in infected roots of Chinese cabbage. This result was verified by reverse transcription-qPCR and in situ hybridization. Examination of disease resistance showed that, compared with wild type plants, A. thaliana ATG6 deletion mutants were more easily infected by P. brassicae than WT. This shows that BrATG6 may play a potential role in the resistance of B. rapa to P. brassicae infection.

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.