Transcriptome Arofile of Brassica rapa L. Reveals the Involvement of Jasmonic Acid, Ethylene, and Brassinosteroid Signaling Pathways in Clubroot Resistance

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.

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.

Transcriptional Response of Wolfberry to Infestation with the Endophytic Fusarium nematophilum Strain NQ8GII4

Fusarium nematophilum NQ8GII4 is an endophytic fungus isolated from the root of healthy wolfberry (Lycium barbarum). Previous studies have reported that NQ8GII4 could dwell in wolfberry roots and enhance the defense responses in wolfberry against root rot, which is caused by F. oxysporum. To further elucidate the molecular mechanism of wolfberry disease resistance induced by NQ8GII4, in the present study, we adopted RNA sequencing analysis to profile the transcriptome of wolfberry response to NQ8GII4 infestation over a time course of 3 and 7 days postinoculation. Gene ontology enrichment analysis revealed that differentially expressed genes (DEGs) were enriched in biological regulation, response to stimulus, signaling, detoxification, immune system process, transporter activity, electron carrier activity, transcription factor activity, nucleic acid binding transcription factor, and antioxidant activity. Through Kyoto Encyclopedia of Genes and Genomes analysis, it was found that many of these DEGs were enriched in pathways related to plant-pathogen interactions, hormone signal transduction, and the phenylpropanoid biosynthesis pathway in wolfberry. This result suggested that innate immunity, phytohormone signaling, and numerous phenylpropanoid compounds comprise a complex defense network in wolfberry. Chloroplast 50S ribosomal proteins were consistently located at the core position of the response in wolfberry following infestation with NQ8GII4 analyzed by the protein-protein interaction network. This study elucidated the molecular mechanism underlying the interaction between NQ8GII4 and wolfberry, clarified the wolfberry immune response network to endophytic fungi infestation, identified candidate resistance genes in wolfberry, and provided a fundamental date for subsequent work.

Metabonomic profiling of clubroot-susceptible and clubroot-resistant radish and the assessment of disease-resistant metabolites

Plasmodiophora brassicae causes a serious threat to cruciferous plants including radish (Raphanus sativus L.). Knowledge on the pathogenic regularity and molecular mechanism of P. brassicae and radish is limited, especially on the metabolism level. In the present study, clubroot-susceptible and clubroot-resistant cultivars were inoculated with P. brassicae Race 4, root hairs initial infection of resting spores (10⁷ CFU/mL) at 24 h post-inoculation and root galls symptom arising at cortex splitting stage were identified on both cultivars. Root samples of cortex splitting stage of two cultivars were collected and used for untargeted metabonomic analysis. We demonstrated changes in metabolite regulation and pathways during the cortex splitting stage of diseased roots between clubroot-susceptible and clubroot-resistant cultivars using untargeted metabonomic analysis. We identified a larger number of differentially regulated metabolites and heavier metabolite profile changes in the susceptible cultivar than in the resistant counterpart. The metabolites that were differentially regulated in both cultivars were mostly lipids and lipid-like molecules. Significantly regulated metabolites and pathways according to the P value and variable important in projection score were identified. Moreover, four compounds, including ethyl α-D-thioglucopyranoside, imipenem, ginsenoside Rg1, and 6-gingerol, were selected, and their anti-P. brassicae ability and effects on seedling growth were verified on the susceptible cultivar. Except for ethyl α-D-thioglucopyranoside, the remaining could inhibit clubroot development of varing degree. The use of 5 mg/L ginsenoside Rg1 + 5 mg/L 6-gingerol resulted in the lowest disease incidence and disease index among all treatments and enhanced seedling growth. The regulation of pathways or metabolites of carbapenem and ginsenoside was further explored. The results provide a preliminary understanding of the interaction between radish and P. brassicae at the metabolism level, as well as the development of measures for preventing clubroot.

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.

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.

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.

Identification and Characterization of Circular RNAs in Brassica rapa in Response to Plasmodiophora brassicae

Plasmodiophora brassicae is a soil-borne pathogen that attacks the roots of cruciferous plants and causes clubroot disease. CircRNAs are noncoding RNAs, widely existing in plant and animal species. Although knowledge of circRNAs has been updated continuously and rapidly, information about circRNAs in the regulation of clubroot disease resistance is extremely limited in Brassica rapa. Here, Chinese cabbage (BJN 222) containing clubroot resistance genes (CRa) against P. brassicae Pb4 was susceptible to PbE. To investigate the mechanism of cicRNAs responsible for clubroot disease resistance in B. rapa, circRNA-seq was performed with roots of 'BJN 222' at 0, 8, and 23 days post-inoculated (dpi) with Pb4 and PbE. A total of 231 differentially expressed circRNAs were identified between the groups. Based on the differentially expressed circRNAs, the circRNA-miRNA-mRNA network was constructed using the target genes directly or indirectly related to plant resistance. Upregulated novel_circ_000495 suppressed the expression of miR5656-y, leading to the upregulation of Bra026508, which might cause plant resistance. Our results provide new insights into clubroot resistance mechanisms and lay a foundation for further studies exploring complex gene regulation networks in B. rapa.

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.

Identification of Arabidopsis Phospholipase A Mutants With Increased Susceptibility to Plasmodiophora brassicae

Clubroot, caused by the obligate parasite Plasmodiophora brassicae, is one of the most devastating diseases of canola (Brassica napus) in Canada. The identification of novel genes that contribute to clubroot resistance is important for the sustainable management of clubroot, as these genes may be used in the development of resistant canola cultivars. Phospholipase As (PLAs) play important roles in plant defense signaling and stress tolerance, and thus are attractive targets for crop breeding. However, since canola is an allopolyploid and has multiple copies of each PLA gene, it is time-consuming to test the functions of PLAs directly in this crop. In contrast, the model plant Arabidopsis thaliana has a simpler genetic background and only one copy of each PLA. Therefore, it would be reasonable and faster to validate the potential utility of PLA genes in Arabidopsis first. In this study, we identified seven homozygous atpla knockout/knockdown mutants of Arabidopsis, and tested their performance following inoculation with P. brassicae. Four mutants (pla₁-iiα, pla₁-iγ3, pla₁-iii, ppla-iiiβ, ppla-iiiδ) developed more severe clubroot than the wild-type, suggesting increased susceptibility to P. brassicae. The homologs of these Arabidopsis PLAs (AtPLAs) in B. napus (BnPLAs) were identified through Blast searches and phylogenic analysis. Expression of the BnPLAs was subsequently examined in transcriptomic datasets generated from canola infected by P. brassicae, and promising candidates for further characterization identified.

Endomembrane-Targeting Plasmodiophora brassicae Effectors Modulate PAMP Triggered Immune Responses in Plants

Plasmodiophora brassicae is a devastating obligate, intracellular, biotrophic pathogen that causes clubroot disease in crucifer plants. Disease progression is regulated by effector proteins secreted by P. brassicae. Twelve P. brassicae putative effectors (PbPEs), expressed at various stages of disease development [0, 2, 5, 7, 14, 21, and 28 days post inoculation (DPI)] in Arabidopsis and localizing to the plant endomembrane system, were studied for their roles in pathogenesis. Of the 12 PbPEs, seven showed an inhibitory effect on programmed cell death (PCD) as triggered by the PCD inducers, PiINF1 (Phytophthora infestans Infestin 1) and PiNPP1 (P. infestans necrosis causing protein). Showing the strongest level of PCD suppression, PbPE15, a member of the 2-oxoglutarate (2OG) and Fe (II)-dependent oxygenase superfamily and with gene expression during later stages of infection, appears to have a role in tumorigenesis as well as defense signaling in plants. PbPE13 produced an enhanced PiINF1-induced PCD response. Transient expression, in Nicotiana benthamiana leaves of these PbPEs minus the signal peptide (SP) (^Δsp PbPEGFPs), showed localization to the endomembrane system, targeting the endoplasmic reticulum (ER), Golgi bodies and nucleo-cytoplasm, suggesting roles in manipulating plant cell secretion and vesicle trafficking. ^Δsp PbPE13GFP localized to plasma membrane (PM) lipid rafts with an association to plasmodesmata, suggesting a role at the cell-to-cell communication junction. Membrane relocalization of ^Δsp PbPE13GFP, triggered by flagellin N-terminus of Pseudomonas aeruginosa (flg22 - known to elicit a PAMP triggered immune response in plants), supports its involvement in raft-mediated immune signaling. This study is an important step in deciphering P. brassicae effector roles in the disruption of plant immunity to clubroot disease.

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

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

Root Transcriptome and Metabolome Profiling Reveal Key Phytohormone-Related Genes and Pathways Involved Clubroot Resistance in Brassica rapa L

Plasmodiophora brassicae, an obligate biotrophic pathogen-causing clubroot disease, can seriously affect Brassica crops worldwide, especially Chinese cabbage. Understanding the transcriptome and metabolome profiling changes during the infection of P. brassicae will provide key insights in understanding the defense mechanism in Brassica crops. In this study, we estimated the phytohormones using targeted metabolome assays and transcriptomic changes using RNA sequencing (RNA-seq) in the roots of resistant (BrT24) and susceptible (Y510-9) plants at 0, 3, 9, and 20 days after inoculation (DAI) with P. brassicae. Differentially expressed genes (DEGs) in resistant vs. susceptible lines across different time points were identified. The weighted gene co-expression network analysis of the DEGs revealed six pathways including "Plant-pathogen interaction" and "Plant hormone signal transduction" and 15 hub genes including pathogenic type III effector avirulence factor gene (RIN4) and auxin-responsive protein (IAA16) to be involved in plants immune response. Inhibition of Indoleacetic acid, cytokinin, jasmonate acid, and salicylic acid contents and changes in related gene expression in R-line may play important roles in regulation of clubroot resistance (CR). Based on the combined metabolome profiling and hormone-related transcriptomic responses, we propose a general model of hormone-mediated defense mechanism. This study definitely enhances our current understanding and paves the way for improving CR in Brassica rapa.

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

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

CRb and PbBa8.1 Synergically Increases Resistant Genes Expression upon Infection of Plasmodiophora brassicae in Brassica napus

PbBa8.1 and CRb are two clubroot-resistant genes that are important for canola breeding in China. Previously, we combined these resistant genes and developed a pyramid-based, homozygous recurrent inbred line (618R), the results of which showed strong resistance to Plasmodiophora brassicae field isolates; however, the genetic mechanisms of resistance were unclear. In the present work, we conducted comparative RNA sequencing (RNA-Seq) analysis between 618R and its parental lines (305R and 409R) in order to uncover the transcriptomic response of the superior defense mechanisms of 618R and to determine how these two different resistant genes coordinate with each other. Here, we elucidated that the number and expression of differentially expressed genes (DEGs) in 618R are significantly higher than in the parental lines, and PbBa8.1 shares more DEGs and plays a dominant role in the pyramided line. The common DEGs among the lines largely exhibit non-additive expression patterns and enrichment in resistance pathways. Among the enriched pathways, plant-pathogen interaction, plant hormone signaling transduction, and secondary metabolites are the key observation. However, the expressions of the salicylic acid (SA) signaling pathway and reactive oxygen species (ROS) appear to be crucial regulatory components in defense response. Our findings provide comprehensive transcriptomic insight into understanding the interactions of resistance gene pyramids in single lines and can facilitate the breeding of improved resistance in Brassica napus.