Refining the Life Cycle of Plasmodiophora brassicae

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.

Biological Control of Rapeseed Clubroot (Plasmodiophora brassicae) Using the Endophytic Fungus Didymella macrostoma P2

Didymella macrostoma P2 was isolated from rapeseed (Brassica napus), and it is an endophyte of rapeseed and an antagonist of three rapeseed pathogens, Botrytis cinerea, Leptosphaeria biglobosa, and Sclerotinia sclerotiorum. However, whether P2 has a suppressive effect on infection of rapeseed by the clubroot pathogen Plasmodiophora brassicae remains unknown. This study was conducted to detect production of antimicrobials by P2 and to determine the efficacy of the antimicrobials and P2 pycnidiospores in suppression of rapeseed clubroot. The results showed that cultural filtrates (CFs) of P2 in potato dextrose broth and the substances in pycnidiospore mucilages exuded from P2 pycnidia were inhibitory to P. brassicae. In the indoor experiment, seeds of the susceptible rapeseed cultivar Zhongshuang No. 9 treated with P2 CF and the P2 pycnidiospore suspension (P2 SS, 1 × 107 spores/ml) reduced clubroot severity by 31 to 70% on the 30-day-old seedlings compared with the control (seeds treated with water). P2 was reisolated from the roots of the seedlings in the treatment of P2 SS; the average isolation frequency in the healthy roots (26%) was much higher than that (5%) in the diseased roots. In the field experiment, seeds of another susceptible rapeseed cultivar, Huayouza 50 (HYZ50), treated with P2 CF, P2 CE (chloroform extract of P2 CF, 30 μg/ml), and P2 SS reduced clubroot severity by 29 to 48% on 60-day-old seedlings and by 28 to 59% on adult plants (220 days old) compared with the control treatment. The three P2 treatments on HYZ50 produced significantly (P < 0.05) higher seed yield than the control treatment on this rapeseed cultivar, and they even generated seed yield similar to that produced by the resistant rapeseed cultivar Shengguang 165R in one of the two seasons. These results suggest that D. macrostoma P2 is an effective biocontrol agent against rapeseed clubroot.

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.

Identification and Characterization of High-Molecular-Weight Proteins Secreted by Plasmodiophora brassicae That Suppress Plant Immunity

Plasmodiophora brassicae is an obligate intracellular parasitic protist that causes clubroot disease on cruciferous plants. So far, some low-molecular-weight secreted proteins from P. brassicae have been reported to play an important role in plant immunity regulation, but there are few reports on its high-molecular-weight secreted proteins. In this study, 35 putative high-molecular-weight secreted proteins (>300 amino acids) of P. brassicae (PbHMWSP) genes that are highly expressed during the infection stage were identified using transcriptome analysis and bioinformatics prediction. Then, the secretory activity of 30 putative PbHMWSPs was confirmed using the yeast signal sequence trap system. Furthermore, the genes encoding 24 PbHMWSPs were successfully cloned and their functions in plant immunity were studied. The results showed that ten PbHMWSPs could inhibit flg22-induced reactive oxygen burst, and ten PbHMWSPs significantly inhibited the expression of the SA signaling pathway marker gene PR1a. In addition, nine PbHMWSPs could inhibit the expression of a marker gene of the JA signaling pathway. Therefore, a total of 19 of the 24 tested PbHMWSPs played roles in suppressing the immune response of plants. Of these, it is worth noting that PbHMWSP34 can inhibit the expression of JA, ET, and several SA signaling pathway marker genes. The present study is the first to report the function of the high-molecular-weight secreted proteins of P. brassicae in plant immunity, which will enrich the theory of interaction mechanisms between the pathogens and plants.

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.

Telomere-to-telomere Genome Assembly of the Clubroot Pathogen Plasmodiophora Brassicae

Plasmodiophora brassicae (Woronin, 1877), a biotrophic, obligate parasite, is the causal agent of clubroot disease in brassicas. The clubroot pathogen has been reported in more than 80 countries worldwide, causing economic losses of hundreds of millions every year. Despite its widespread impact, very little is known about the molecular strategies it employs to induce the characteristic clubs in the roots of susceptible hosts during infection, nor about the mechanisms it uses to overcome genetic resistance. Here, we provide the first telomere-to-telomere complete genome of P. brassicae. We generated ∼27 Gb of Illumina, Oxford Nanopore, and PacBio HiFi data from resting spores of strain Pb3A and produced a 25.3 Mb assembly comprising 20 chromosomes, with an N50 of 1.37 Mb. The BUSCO score, the highest reported for any member of the group Rhizaria (Eukaryota: 88.2%), highlights the limitations within the Eukaryota database for members of this lineage. Using available transcriptomic data and protein evidence, we annotated the Pb3A genome, identifying 10,521 protein-coding gene models. This high-quality, complete genome of P. brassicae will serve as a crucial resource for the plant pathology community to advance the much-needed understanding of the evolution of the clubroot pathogen.

The Plasmodiophora brassicae Golgi-localized UPF0016 protein PbGDT1 mediates calcium but not manganese transport in yeast and Nicotiana benthamiana

Manganese and calcium homeostasis and signalling, in eukaryotic organisms, are regulated through membrane located pumps, channels and exchangers, including the Mn2+/Ca2+ uncharacterized protein family 0016 (UPF0016). Here we show that Plasmodiophora brassicae PbGDT1 is a member of the UPF0016 and an ortholog of Saccharomyces cerevisiae Gdt1p (GCR Dependent Translation Factor 1) protein involved in manganese homeostasis as well as the calcium mediated stress response in yeast. PbGDT1 complemented the ScGdt1p and ScPMR1 (Ca2+ ATPase) double null mutant under elevated calcium stress but not under elevated manganese conditions. In both yeast and Nicotiana benthamiana, PbGDT1 localizes to the Golgi apparatus, with additional ER association in N. benthamiana. Expression of PbGDT1 in N. benthamiana, suppresses BAX-triggered cell death, further highlighting the importance of calcium homeostasis in maintaining cell physiology and integrity in a stress environment.

Environmental conditions modulate the effect of epigenetic factors controlling the response of Arabidopsis thaliana to Plasmodiophora brassicae

The resistance of Arabidopsis thaliana to clubroot, a major disease of Brassicaceae caused by the obligate protist Plasmodiophora brassicae, is controlled in part by epigenetic factors. The detection of some of these epigenetic quantitative trait loci (QTLepi) has been shown to depend on experimental conditions. The aim of the present study was to assess whether and how temperature and/or soil water availability influenced both the detection and the extent of the effect of response QTLepi. The epigenetic recombinant inbred line (epiRIL) population, derived from the cross between ddm1-2 and Col-0 (partially resistant and susceptible to clubroot, respectively), was phenotyped for response to P. brassicae under four abiotic conditions including standard conditions, a 5°C temperature increase, drought, and flooding. The abiotic constraints tested had a significant impact on both the leaf growth of the epiRIL population and the outcome of the epiRIL-pathogen interaction. Linkage analysis led to the detection of a total of 31 QTLepi, 18 of which were specific to one abiotic condition and 13 common to at least two environments. EpiRIL showed significant plasticity under epigenetic control, which appeared to be specific to the traits evaluated and to the abiotic conditions. These results highlight that the environment can affect the epigenetic architecture of plant growth and immune responses and advance our understanding of the epigenetic factors underlying plasticity in response to climate change.

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.

Magnesium oxide nanoparticles reduce clubroot by regulating plant defense response and rhizosphere microbial community of tumorous stem mustard (Brassica juncea var. tumida)

Clubroot, caused by Plasmodiophora brassicae, is a major disease that significantly impairs the yield of cruciferous crops and causes significant economic losses across the globe. The prevention of clubroot, especially in tumorous stem mustard (without resistant varieties), are is limited and primarily relies on fungicides. Engineered nanoparticles have opened up new avenues for the management of plant diseases, but there is no report on their application in the prevention of clubroot. The results showed that the control efficacy of 500 mg/L MgO NPs against clubroot was 54.92%. However, when the concentration was increased to 1,500 and 2,500 mg/L, there was no significant change in the control effect. Compared with CK, the average fresh and dry weight of the aerial part of plants treated with MgO NPs increased by 392.83 and 240.81%, respectively. Compared with the F1000 treatment, increases were observed in the content of soil available phosphorus (+16.72%), potassium (+9.82%), exchangeable magnesium (+24.20%), and water-soluble magnesium (+20.64%) in the 1,500 mg/L MgO NPs treatment. The enzyme-linked immune sorbent assay (ELISA) results showed that the application of MgO NPs significantly increased soil peroxidase (POD, +52.69%), alkaline protease (AP, +41.21%), alkaline phosphatase (ALP, +79.26%), urease (+52.69%), and sucrase (+56.88%) activities; And also increased plant L-phenylalanine ammonla-lyase (PAL, +70.49%), polyphenol oxidase (PPO, +36.77%), POD (+38.30%), guaiacol peroxidase (POX, +55.46%) activities and salicylic acid (SA, +59.86%) content. However, soil and plant catalase (CAT, -27.22 and - 19.89%, respectively), and plant super oxidase dismutase (SOD, -36.33%) activities were significantly decreased after the application of MgO NPs. The metagenomic sequencing analysis showed that the MgO NPs treatments significantly improved the α-diversity of the rhizosphere soil microbial community. The relative abundance of beneficial bacteria genera in the rhizosphere soil, including Pseudomonas, Sphingopyxis, Acidovorax, Variovorax, and Bosea, was significantly increased. Soil metabolic functions, such as oxidative phosphorylation (ko00190), carbon fixation pathways in prokaryotes (ko00720), indole alkaloid biosynthesis (ko00901), and biosynthesis of various antibiotics (ko00998) were significantly enriched. These results suggested that MgO NPs might control clubroot by promoting the transformation and utilization of soil nutrients, stimulating plant defense responses, and enriching soil beneficial bacteria.

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.

Testing Effects of Seed Treatments against Clubroot Disease in Various Oilseed Rape Hybrids

Clubroot disease, caused by the protist pathogen Plasmodiophora brassicae, is an emerging threat to cruciferous crops, including oilseed rape (Brassica napus L.). Most of the current commercial cultivars are highly susceptible, and efficient management tools are lacking practical implementation. Over three years and three experimental periods, we studied the effects of isotianil in comparison with Bacillus amyloliquefaciens QST713-HiCFU against clubroot disease under greenhouse experiments. Our results show control effects, which were strongly dependent on seasons, host plant genotype, and clubroot isolates: isotianil and B. amyloliquefaciens QST713-HiCFU reduced disease severity consistently at variable, but field-relevant spore concentrations of clubroot isolates; with seed treatments showing superior effects compared to drench applications. The co-application of isotianil with B. amyloliquefaciens QST713-HiCFU could, in some cases, increase the efficacy. Interestingly, all studied hybrids reacted to treatments, albeit to a somewhat different extent. When tested against a field isolate, the results obtained with the single spore isolate were partially confirmed but with greater variability. Overall, the generally positive effects of isotianil and B. amyloliquefaciens QST713-HiCFU on the reduction of clubroot were repeatedly observed. The inoculation of clubroot disease with different spore counts indicates a dose-response effect for tested products. This study highlights the importance of performing experiments holistically over multiple, consecutive seasons, with various isolates, application types, and different genetic resources of host plants.

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.

Phagocytosis underpins the biotrophic lifestyle of intracellular parasites in the class Phytomyxea (Rhizaria)

Phytomyxea are intracellular biotrophic parasites infecting plants and stramenopiles, including the agriculturally impactful Plasmodiophora brassicae and the brown seaweed pathogen Maullinia ectocarpii. They belong to the clade Rhizaria, where phagotrophy is the main mode of nutrition. Phagocytosis is a complex trait of eukaryotes, well documented for free-living unicellular eukaryotes and specific cellular types of animals. Data on phagocytosis in intracellular, biotrophic parasites are scant. Phagocytosis, where parts of the host cell are consumed at once, is seemingly at odds with intracellular biotrophy. Here we provide evidence that phagotrophy is part of the nutritional strategy of Phytomyxea, using morphological and genetic data (including a novel transcriptome of M. ectocarpii). We document intracellular phagocytosis in P. brassicae and M. ectocarpii by transmission electron microscopy and fluorescent in situ hybridization. Our investigations confirm molecular signatures of phagocytosis in Phytomyxea and hint at a small specialized subset of genes used for intracellular phagocytosis. Microscopic evidence confirms the existence of intracellular phagocytosis, which in Phytomyxea targets primarily host organelles. Phagocytosis seems to coexist with the manipulation of host physiology typical of biotrophic interactions. Our findings resolve long debated questions on the feeding behaviour of Phytomyxea, suggesting an unrecognized role for phagocytosis in biotrophic interactions.

Silicon ameliorates clubroot responses in canola (Brassica napus): A "multi-omics"-based investigation into possible mechanisms

Clubroot disease, caused by Plasmodiophora brassicae Woronin, results in severe yield losses in Brassica crops, including canola. Silicon (Si) mitigates several stresses and enhances plant resistance to phytopathogens. We investigated the effects of Si on clubroot disease symptoms in canola at two concentrations of Si, Si: soil in 1: 100 w/w (Si1.0) and Si: soil in 1:200 w/w (Si0.5) under greenhouse conditions. In addition, the effects of Si on P. brassicae-induced gene expression, endogenous levels of phytohormones and metabolites were studied using "omics" approaches. Si application reduced clubroot symptoms and improved plant growth parameters. Gene expression analysis revealed increased transcript-level responses in Si1.0 compared to Si0.5 plants at 7-, 14-, and 21-days post-inoculation (dpi). Pathogen-induced transcript-level changes were affected by Si treatment, with genes related to antioxidant activity (e.g., POD, CAT), phytohormone biosynthesis and signalling (e.g., PDF1.2, NPR1, JAZ, IPT, TAA), nitrogen metabolism (e.g., NRT, AAT), and secondary metabolism (e.g., PAL, BCAT4) exhibiting differential expression. Endogenous levels of phytohormones (e.g., auxin, cytokinin), a majority of the amino acids and secondary metabolites (e.g., glucosinolates) were increased at 7 dpi, followed by a decrease at 14- and 21-dpi due to Si-treatment. Stress hormones such as abscisic acid (ABA), salicylic acid (SA), and jasmonic acid (JA) also decreased at the later time points in Si0.5, and Si1.0 treated plants. Si appears to improve clubroot symptoms while enhancing plant growth and associated metabolic processes, including nitrogen metabolism and secondary metabolite biosynthesis.

Getting to the root of a club - Understanding developmental manipulation by the clubroot pathogen

Plasmodiophora brassicae Wor., the clubroot pathogen, is the perfect example of an "atypical" plant pathogen. This soil-borne protist and obligate biotrophic parasite infects the roots of cruciferous crops, inducing galls or clubs that lead to wilting, loss of productivity, and plant death. Unlike many other agriculturally relevant pathosystems, research into the molecular mechanisms that underlie clubroot disease and Plasmodiophora-host interactions is limited. After release of the first P. brassicae genome sequence and subsequent availability of transcriptomic data, the clubroot research community have implicated the involvement of phytohormones during the clubroot pathogen's manipulation of host development. Herein we review the main events leading to the formation of root galls and describe how modulation of select phytohormones may be key to modulating development of the plant host to the benefit of the pathogen. Effector-host interactions are at the base of different strategies employed by pathogens to hijack plant cellular processes. This is how we suspect the clubroot pathogen hijacks host plant metabolism and development to induce nutrient-sink roots galls, emphasizing a need to deepen our understanding of this master manipulator.

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.

The plant trans-Golgi network component ECHIDNA regulates defense, cell death, and endoplasmic reticulum stress

The trans-Golgi network (TGN) acts as a central platform for sorting and secreting various cargoes to the cell surface, thus being essential for the full execution of plant immunity. However, the fine-tuned regulation of TGN components in plant defense and stress response has been not fully elucidated. Our study revealed that despite largely compromising penetration resistance, the loss-of-function mutation of the TGN component protein ECHIDNA (ECH) induced enhanced postinvasion resistance to powdery mildew in Arabidopsis thaliana. Genetic and transcriptome analyses and hormone profiling demonstrated that ECH loss resulted in salicylic acid (SA) hyperaccumulation via the ISOCHORISMATE SYNTHASE 1 biosynthesis pathway, thereby constitutively activating SA-dependent innate immunity that was largely responsible for the enhanced postinvasion resistance. Furthermore, the ech mutant displayed accelerated SA-independent spontaneous cell death and constitutive POWDERY MILDEW RESISTANCE 4-mediated callose depositions. In addition, ECH loss led to a chronically prolonged endoplasmic reticulum stress in the ech mutant. These results provide insights into understanding the role of TGN components in the regulation of plant immunity and stress responses.

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.

Comparing the infection biology and gene expression differences of Plasmodiophora brassicae primary and secondary zoospores

Plasmodiophora brassicae (Wor.) is an obligate plant pathogen affecting Brassicae worldwide. To date, there is very little information available on the biology and molecular basis of P. brassicae primary and secondary zoospore infections. To examine their roles, we used microscope to systematically investigate the infection differences of P. brassicae between samples inoculated separately with resting spores and secondary zoospores. The obvious development of P. brassicae asynchrony that is characterized by secondary plasmodium, resting sporangial plasmodium, and resting spores was observed at 12 days in Brassica rapa inoculated with resting spores but not when inoculated with secondary zoospores at the same time. Inoculation with resting spores resulted in much more development of zoosporangia clusters than inoculation with secondary zoospores in non-host Spinacia oleracea. The results indicated that primary zoospore infection played an important role in the subsequent development. To improve our understanding of the infection mechanisms, RNA-seq analysis was performed. Among 18 effectors identified in P. brassicae, 13 effectors were induced in B. rapa seedlings inoculated with resting spores, which suggested that the pathogen and host first contacted, and more effectors were needed. Corresponding to those in B. rapa, the expression levels of most genes involved in the calcium-mediated signaling pathway and PTI pathway were higher in plants inoculated with resting spores than in those inoculated with secondary zoospores. The ETI pathway was suppressed after inoculation with secondary zoospores. The genes induced after inoculation with resting spores were suppressed in B. rapa seedlings inoculated with secondary zoospores, which might be important to allow a fully compatible interaction and contribute to a susceptible reaction in the host at the subsequent infection stage. The primary zoospores undertook an more important interaction with plants.

Ethylene Plays a Dual Role during Infection by Plasmodiophora brassicae of Arabidopsis thaliana

Plasmodiophora brassicae infection leads to hypertrophy of host roots and subsequent formation of galls, causing huge economic losses to agricultural producers of Cruciferae plants. Ethylene (ET) has been reported to play a vital role against necrotrophic pathogens in the classic immunity system. More clues suggested that the defense to pathogens in roots may be different from the acrial. The ET pathway may play a positive role in the infection of P. brassicae, as shown by recent transcriptome profiling. However, the molecular basis of ET remains poorly understood. In this study, we investigated the potential role of ethylene against P. brassicae infection in an ein3/eil1 double-mutant of Arabidopsis thaliana (A. thaliana). After infection, ein3/eil1 (Disease Index/DI: 93) showed more susceptibility compared with wild type (DI: 75). Then, we inoculated A. thaliana Columbia-0 (Col-0) with P. brassicae by 1-aminocyclopropane-1-carboxylic acid (ACC) and pyrazinamide (PZA), respectively. It was found that the symptoms of infected roots with ACC were more serious than those with PZA at 20 dpi (day post infection). However, the DI were almost the same in different treatments at 30 dpi. WRKY75 can be directly regulated by ET and was upregulated at 7 dpi with ACC, as shown by qRT-PCR. The wrky75-c mutant of A. thaliana (DI: 93.75) was more susceptible than the wild type in Arabidopsis. Thus, our work reveals the dual roles of ET in infection of P. brassicae and provides evidence of ET in root defense against pathogens.

What Can We Learn from -Omics Approaches to Understand Clubroot Disease?

Clubroot is one of the most economically significant diseases worldwide. As a result, many investigations focus on both curing the disease and in-depth molecular studies. Although the first transcriptome dataset for the clubroot disease describing the clubroot disease was published in 2006, many different pathogen-host plant combinations have only recently been investigated and published. Articles presenting -omics data and the clubroot pathogen Plasmodiophora brassicae as well as different host plants were analyzed to summarize the findings in the richness of these datasets. Although genome data for the protist have only recently become available, many effector candidates have been identified, but their functional characterization is incomplete. A better understanding of the life cycle is clearly required to comprehend its function. While only a few proteome studies and metabolome analyses were performed, the majority of studies used microarrays and RNAseq approaches to study transcriptomes. Metabolites, comprising chemical groups like hormones were generally studied in a more targeted manner. Furthermore, functional approaches based on such datasets have been carried out employing mutants, transgenic lines, or ecotypes/cultivars of either Arabidopsis thaliana or other economically important host plants of the Brassica family. This has led to new discoveries of potential genes involved in disease development or in (partial) resistance or tolerance to P. brassicae. The overall contribution of individual experimental setups to a larger picture will be discussed in this review.

A Loop-Mediated Isothermal DNA Amplification (LAMP) Assay for Detection of the Clubroot Pathogen Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a serious threat to cruciferous crops around the world. The resting spores of P. brassicae are a primary source of infection and can survive in soil for many years. Detection of resting spores in soil is essential for forecasting clubroot prevalence. Detection of P. brassicae has been relying on plant bioassays or PCR-based methods. The loop-mediated isothermal DNA amplification (LAMP) is a promising approach for microorganism detection with the advantage of high sensitivity, accuracy, and convenience in viewing. In this study, we developed a LAMP assay for detection of P. brassicae in soil, roots, and seeds. This method can detect P. brassicae at a minimal amount of 1 fg of plasmid DNA or 10 resting spores in the soil. Compared with conventional PCR, the LAMP was more sensitive in detection of P. brassicae at the lower levels in soil samples. In conclusion, we elaborated a sensitive, accurate, and easy-to-use LAMP assay to detect P. brassicae, which will facilitate sustainable clubroot management and planning.

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.

Sugar Transporters in Plasmodiophora brassicae: Genome-Wide Identification and Functional Verification

Plasmodiophora brassicae, an obligate intracellular pathogen, can hijack the host’s carbohydrates for survival. When the host plant is infected by P. brassicae, a large amount of soluble sugar accumulates in the roots, especially glucose, which probably facilitates the development of this pathogen. Although a complete glycolytic and tricarboxylic acid cycle (TCA) cycle existed in P. brassicae, very little information about the hexose transport system has been reported. In this study, we screened 17 putative sugar transporters based on information about their typical domains. The structure of these transporters showed a lot of variation compared with that of other organisms, especially the number of transmembrane helices (TMHs). Phylogenetic analysis indicated that these sugar transporters were far from the evolutionary relationship of other organisms and were unique in P. brassicae. The hexose transport activity assay indicated that eight transporters transported glucose or fructose and could restore the growth of yeast strain EBY.VW4000, which was deficient in hexose transport. The expression level of these glucose transporters was significantly upregulated at the late inoculation time when resting spores and galls were developing and a large amount of energy was needed. Our study provides new insights into the mechanism of P. brassicae survival in host cells by hijacking and utilizing the carbohydrates of the host.

Comparative transcriptome analysis reveals distinct responsive biological processes in radish genotypes contrasting for Plasmodiophora brassicae interaction

Plasmodiophora brassicae is a protozoan pathogen that causes clubroot disease, which is one of the most destructive diseases for Brassica crops, including radish. However, little is known about the molecular mechanism of clubroot resistance in radish. In this study, we performed a comparative transcriptome analysis between resistant and susceptible radish inoculated with P. brassicae. More differentially expressed genes (DEGs) were identified at 28 days after inoculation (DAI) compared to 7 DAI in both genotypes. Gene ontology (GO) and KEGG enrichment indicated that stress/defense response, secondary metabolic biosynthesis, hormone metabolic process, and cell periphery are directly involved in the defense response process. Further analysis of the transcriptome revealed that effector-triggered immunity (ETI) plays key roles in the defense response. The plant hormones jasmonic acid (JA), ethylene (ET), and abscisic acid (ABA) related genes are activated in clubroot defense in the resistant line. Auxin (AUX) hormone related genes are activated in the developing galls of susceptible radish. Our study provides a global transcriptional overview for clubroot development for insights into the P. brassicae defense mechanisms in radish.

Starch content changes and metabolism-related gene regulation of Chinese cabbage synergistically induced by Plasmodiophora brassicae infection

Clubroot is one of the major diseases adversely affecting Chinese cabbage (Brassica rapa) yield and quality. To precisely characterize the Plasmodiophora brassicae infection on Chinese cabbage, we developed a dual fluorescent staining method for simultaneously examining the pathogen, cell structures, and starch grains. The number of starch (amylopectin) grains increased in B. rapa roots infected by P. brassicae, especially from 14 to 21 days after inoculation. Therefore, the expression levels of 38 core starch metabolism genes were investigated by quantitative real-time PCR. Most genes related to starch synthesis were up-regulated at seven days after the P. brassicae inoculation, whereas the expression levels of the starch degradation-related genes increased at 14 days after the inoculation. Then genes encoding the core enzymes involved in starch metabolism were investigated by assessing their chromosomal distributions, structures, duplication events, and synteny among Brassica species. Genome comparisons indicated that 38 non-redundant genes belonging to six core gene families related to starch metabolism are highly conserved among Arabidopsis thaliana, B. rapa, Brassica nigra, and Brassica oleracea. Genome sequencing projects have revealed that P. brassicae obtained host nutrients by manipulating plant metabolism. Starch may serve as a carbon source for P. brassicae colonization as indicated by the histological observation and transcriptomic analysis. Results of this study may elucidate the evolution and expression of core starch metabolism genes and provide researchers with novel insights into the pathogenesis of clubroot in B. rapa.

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.

A Novel Target (Oxidation Resistant 2) in Arabidopsis thaliana to Reduce Clubroot Disease Symptoms via the Salicylic Acid Pathway without Growth Penalties

The clubroot disease (Plasmodiophora brassicae) is one of the most damaging diseases worldwide among brassica crops. Its control often relies on resistant cultivars, since the manipulation of the disease hormones, such as salicylic acid (SA) alters plant growth negatively. Alternatively, the SA pathway can be increased by the addition of beneficial microorganisms for biocontrol. However, this potential has not been exhaustively used. In this study, a recently characterized protein Oxidation Resistant 2 (OXR2) from Arabidopsis thaliana is shown to increase the constitutive pathway of SA defense without decreasing plant growth. Plants overexpressing AtOXR2 (OXR2-OE) show strongly reduced clubroot symptoms with improved plant growth performance, in comparison to wild type plants during the course of infection. Consequently, oxr2 mutants are more susceptible to clubroot disease. P. brassicae itself was reduced in these galls as determined by quantitative real-time PCR. Furthermore, we provide evidence for the transcriptional downregulation of the gene encoding a SA-methyltransferase from the pathogen in OXR2-OE plants that could contribute to the phenotype.

Candidate Effectors of Plasmodiophora brassicae Pathotype 5X During Infection of Two Brassica napus Genotypes

Clubroot, caused by Plasmodiophora brassicae, is one of the most important diseases of canola (Brassica napus) in Canada. Disease management relies heavily on planting clubroot resistant (CR) cultivars, but in recent years, new resistance-breaking pathotypes of P. brassicae have emerged. Current efforts against the disease are concentrated in developing host resistance using traditional genetic breeding, omics and molecular biology. However, because of its obligate biotrophic nature, limited resources have been dedicated to investigating molecular mechanisms of pathogenic infection. We previously performed a transcriptomic study with the cultivar resistance-breaking pathotype 5X on two B. napus hosts presenting contrasting resistance/susceptibility, where we evaluated the mechanisms of host response. Since cultivar-pathotype interactions are very specific, and pathotype 5X is one of the most relevant resistance-breaking pathotypes in Canada, in this study, we analyze the expression of genes encoding putative secreted proteins from this pathotype, predicted using a bioinformatics pipeline, protein modeling and orthologous comparisons with effectors from other pathosystems. While host responses were found to differ markedly in our previous study, many common effectors are found in the pathogen while infecting both hosts, and the gene response among biological pathogen replicates seems more consistent in the effectors associated with the susceptible interaction, especially at 21 days after inoculation. The predicted effectors indicate the predominance of proteins with interacting domains (e.g., ankyrin), and genes bearing kinase and NUDIX domains, but also proteins with protective action against reactive oxygen species from the host. Many of these genes confirm previous predictions from other clubroot studies. A benzoic acid/SA methyltransferase (BSMT), which methylates SA to render it inactive, showed high levels of expression in the interactions with both hosts. Interestingly, our data indicate that E3 ubiquitin proteasome elements are also potentially involved in pathogenesis. Finally, a gene with similarity to indole-3-acetaldehyde dehydrogenase is a promising candidate effector because of its involvement in indole acetic acid synthesis, since auxin is one of the major players in clubroot development.

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.

Molecular Pathotyping of Plasmodiophora brassicae-Genomes, Marker Genes, and Obstacles

Here we review the usefulness of the currently available genomic information for the molecular identification of pathotypes. We focused on effector candidates and genes implied to be pathotype specific and tried to connect reported marker genes to Plasmodiophora brassicae genome information. The potentials for practical applications, current obstacles and future perspectives are discussed.

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

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