Expression and Role of Biosynthetic, Transporter, Receptor, and Responsive Genes for Auxin Signaling during Clubroot Disease Development

Bioinformatics and functional analysis of EDS1 genes in Brassica napus in response to Plasmodiophora brassicae infection

Enhanced Disease Susceptibility 1 (EDS1) is a key regulator of plant-pathogen-associated molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) responses. In the Brassica napus genome, we identified six novel EDS1 genes, among which four were responsive to clubroot infection, a major rapeseed disease resistant to chemical control. Developing resistant cultivars is a potent and economically viable strategy to control clubroot infection. Bioinformatics analysis revealed conserved domains and structural uniformity in Bna-EDS1 homologs. Bna-EDS1 promoters harbored elements associated with diverse phytohormones and stress responses, highlighting their crucial roles in plant defense. A functional analysis was performed with Bna-EDS1 overexpression and RNAi transgenic lines. Bna-EDS1 overexpression boosted resistance to clubroot and upregulated defense-associated genes (PR1, PR2, ICS1, and CBP60), while Bna-EDS1 RNAi increased plant susceptibility, indicating suppression of the defense signaling pathway downstream of NBS-LRRs. RNA-Seq analysis identified key transcripts associated with clubroot resistance, including phenylpropanoid biosynthesis. Activation of SA regulator NPR1, defense signaling markers PR1 and PR2, and upregulation of MYC-TFs suggested that EDS1-mediated clubroot resistance potentially involves the SA pathway. Our findings underscore the pivotal role of Bna-EDS1-dependent mechanisms in resistance of B. napus to clubroot disease, and provide valuable insights for fortifying resistance against Plasmodiophora brassicae infection in rapeseed.

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.

Application of the NanoString nCounter System as an Alternative Method to Investigate Molecular Mechanisms Involved in Host Plant Responses to Plasmodiophora brassicae

Clubroot, caused by the soilborne pathogen Plasmodiophora brassicae, is an important disease of canola (Brassica napus) and other crucifers. The recent application of RNA sequencing (RNA-seq) technologies to study P. brassicae−host interactions has generated large amounts of gene expression data, improving knowledge of the molecular mechanisms of pathogenesis and host resistance. Quantitative PCR (qPCR) analysis has been widely applied to examine the expression of a limited number of genes and to validate the results of RNA-seq studies, but may not be ideal for analyzing larger suites of target genes or increased sample numbers. Moreover, the need for intermediate steps such as cDNA synthesis may introduce variability that could affect the accuracy of the data generated by qPCR. Here, we report the validation of gene expression data from a previous RNA-seq study of clubroot using the NanoString nCounter System, which achieves efficient gene expression quantification in a fast and simple manner. We first confirm the robustness of the NanoString system by comparing the results with those generated by qPCR and RNA-seq and then discuss the importance of some candidate genes for resistance or susceptibility to P. brassicae in the host. The results show that the expression of genes measured using NanoString have a high correlation with the values obtained using the other two technologies, with R > 0.90 and p < 0.01, and the same expression patterns for most genes. The three methods (qPCR, RNA-seq, and NanoString) were also compared in terms of laboratory procedures, time, and cost. We propose that the NanoString nCounter System is a robust, sensitive, highly reproducible, and simple technology for gene expression analysis. NanoString could become a common alternative to qPCR to validate RNA-seq data or to create panels of genes for use as markers of resistance/susceptibility when plants are challenged with different P. brassicae pathotypes.

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.

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.

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

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

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

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