Molecular evolution of methylesterase family genes and the BnMES34 is a positive regulator of Plasmodiophora brassicae stress response in Arabidopsis.

Methylesterases (MES) are involved in hydrolysis of carboxylic esters, which have substantial roles in plant metabolic activities and defense mechanisms. This study aimed to comprehensively investigate Brassica napus BnMESs and characterize their role in response to Plasmodiophora brassicae stress. Forty-four BnMES members were identified and categorized into three groups based on their phylogenetic relationships and structural similarities. Through functional predictions in the promoter regions and analysis of RNA-Seq data, BnMES emerged as pivotal in growth, development, and stress responses to B. napus, particularly BnMES34, was strongly induced in response to P. brassicae infection. Gene Ontology analyses highlighted BnMES34's role in regulation of plant disease resistance responses. Furthermore, overexpression of BnMES34 in A. thaliana exhibited milder clubroot symptoms, and reduced disease indices, suggesting positive regulatory role of BnMES34 in plant's response to P. brassicae stress. Molecular docking and enzyme activity verification indicated that BnMES34 has the ability to generate salicylic acid via methyl salicylate, and further experimentally validated in vivo. This discovery indicates that the overexpression of BnMES34 in Arabidopsis confers resistance against clubroot disease. Overall, our research suggests that BnMES34 has a beneficial regulatory role in enhancing stress resistance to P. brassicae in B. napus.

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.

Expression and Role of Response Regulating, Biosynthetic and Degrading Genes for Cytokinin Signaling during Clubroot Disease Development.

The obligate biotroph Plasmodiophora brassicae causes clubroot disease in oilseeds and vegetables of the Brassicaceae family, and cytokinins play a vital role in clubroot formation. In this study, we examined the expression patterns of 17 cytokinin-related genes involved in the biosynthesis, signaling, and degradation in Chinese cabbage inoculated with the Korean pathotype group 4 isolate of P. brassicae, Seosan. This isolate produced the most severe clubroot symptoms in Chinese cabbage cultivar "Bullam-3-ho" compared to three other Korean geographical isolates investigated. BrIPT1, a cytokinin biosynthesis gene, was induced on Day 1 and Day 28 in infected root tissues and the upregulation of this biosynthetic gene coincided with the higher expression of the response regulators BrRR1, on both Days and BrRR6 on Day 1 and 3. BrRR3 and 4 genes were also induced during gall enlargement on Day 35 in leaf tissues. The BrRR4 gene, which positively interact with phytochrome B, was consistently induced in leaf tissues on Day 1, 3, and 14 in the inoculated plants. The cytokinin degrading gene BrCKX3-6 were induced on Day 14, before gall initiation. BrCKX2,3,6 were induced until Day 28 and their expression was downregulated on Day 35. This insight improves our current understanding of the role of cytokinin signaling genes in clubroot disease development.

Transcriptome Analysis Identifies Plasmodiophora brassicae Secondary Infection Effector Candidates.

Plasmodiophora brassicae (Wor.) is an obligate intracellular plant pathogen affecting Brassicas worldwide. Identification of effector proteins is key to understanding the interaction between P. brassicae and its susceptible host plants. To date, there is very little information available on putative effector proteins secreted by P. brassicae during a secondary infection of susceptible host plants, resulting in root gall production. A bioinformatics pipeline approach to RNA-Seq data from Arabidopsis thaliana (L.) Heynh. root tissues at 17, 20, and 24 d postinoculation (dpi) identified 32 small secreted P. brassicae proteins (SSPbPs) that were highly expressed over this secondary infection time frame. Functional signal peptides were confirmed for 31 of the SSPbPs, supporting the accuracy of the pipeline designed to identify secreted proteins. Expression profiles at 0, 2, 5, 7, 14, 21, and 28 dpi verified the involvement of some of the SSPbPs in secondary infection. For seven of the SSPbPs, a functional domain was identified using Blast2GO and 3D structure analysis and domain functionality was confirmed for SSPbP22, a kinase localized to the cytoplasm and nucleus.

Genome-wide identification of genes encoding putative secreted E3 ubiquitin ligases and functional characterization of PbRING1 in the biotrophic protist Plasmodiophora brassicae.

The E3 ubiquitin ligases are key regulators of protein ubiquitination, which have been shown to be involved in a variety of cellular responses to both biotic and abiotic stresses in eukaryotes. However, the E3 ubiquitin ligase homologues in the soil-borne plant pathogen Plasmodiophora brassicae, the causal agent of clubroot disease of crucifer crops worldwide, remain largely unknown. In this study, we characterized secreted E3 ubiquitin ligases, a group of proteins known to be involved in virulence in many pathogens, in a plasmodiophorid P. brassicae. Genome-wide search in the P. brassicae genome retrieved 139 putative E3 ubiquitin ligases, comprising of 115 RING, 15 HECT, 1 HECT-like, and 8 U-box E3 ubiquitin ligases. Among these E3 ubiquitin ligases, 11 RING, 1 U-box, and 3 HECT were found to harbor signal peptide. Based on published RNA-seq data (Schwelm et al. in Sci Rep 5:11153, 2015), we found that these genes were differentially expressed in distinct life stages including germinating spores, maturing spores, and plasmodia. We characterized one potential secreted E3 ubiquitin ligase, PbRING1 (PBRA_000499). Yeast invertase assay showed that PbRING1 harbors a functional N-terminal signal peptide. PbRING1 also harbors a really interested new gene (RING) domain at its C terminus, which was found to display the E3 ligase activity in vitro. Collectively, this study provides a comprehensive insight into the reservoir of putative secreted E3 ligases in P. brassicae.

Computational analysis of the Plasmodiophora brassicae genome: mitochondrial sequence description and metabolic pathway database design.

Plasmodiophora brassicae is an obligate biotrophic pathogenic protist responsible for clubroot, a root gall disease of Brassicaceae species. In addition to the reference genome of the P. brassicae European e3 isolate and the draft genomes of Canadian or Chinese isolates, we present the genome of eH, a second European isolate. Refinement of the annotation of the eH genome led to the identification of the mitochondrial genome sequence, which was found to be bigger than that of Spongospora subterranea, another plant parasitic Plasmodiophorid phylogenetically related to P. brassicae. New pathways were also predicted, such as those for the synthesis of spermidine, a polyamine up-regulated in clubbed regions of roots. A P. brassicae pathway genome database was created to facilitate the functional study of metabolic pathways in transcriptomics approaches. These available tools can help in our understanding of the regulation of P. brassicae metabolism during infection and in response to diverse constraints.

Mimicking the Host Regulation of Salicylic Acid: A Virulence Strategy by the Clubroot Pathogen Plasmodiophora brassicae.

The plant hormone salicylic acid (SA) plays a critical role in defense against biotrophic pathogens such as Plasmodiophora brassicae, which is an obligate pathogen of crucifer species and the causal agent of clubroot disease of canola (Brassica napus). P. brassicae encodes a protein, predicted to be secreted, with very limited homology to benzoic acid (BA)/SA-methyltransferase, designated PbBSMT. PbBSMT has a SA- and an indole-3-acetic acid-binding domain, which are also present in Arabidopsis thaliana BSMT1 (AtBSMT1) and, like AtBSMT1, has been shown to methylate BA and SA. In support of the hypothesis that P. brassicae uses PbBSMT to overcome SA-mediated defenses by converting SA into inactive methyl salicylate (MeSA), here, we show that PbBSMT suppresses local defense and provide evidence that PbBSMT is much more effective than AtBSMT1 at suppressing the levels of SA and its associated effects. Basal SA levels in Arabidopsis plants that constitutively overexpress PbBSMT compared with those in Arabidopsis wild-type Col-0 (WT) were reduced approximately 80% versus only a 50% reduction in plants overexpressing AtBSMT1. PbBSMT-overexpressing plants were more susceptible to P. brassicae than WT plants; they also were partially compromised in nonhost resistance to Albugo candida. In contrast, AtBSMT1-overexpressing plants were not more susceptible than WT to either P. brassicae or A. candida. Furthermore, transgenic Arabidopsis and tobacco plants overexpressing PbBSMT exhibited increased susceptibility to virulent Pseudomonas syringae pv. tomato DC3000 (DC3000) and virulent Pseudomonas syringae pv. tabaci, respectively. Gene-mediated resistance to DC3000/AvrRpt2 and tobacco mosaic virus (TMV) was also compromised in Arabidopsis and Nicotiana tabacum 'Xanthi-nc' plants overexpressing PbBSMT, respectively. Transient expression of PbBSMT or AtBSMT1 in lower leaves of N. tabacum Xanthi-nc resulted in systemic acquired resistance (SAR)-like enhanced resistance to TMV in the distal systemic leaves. Chimeric grafting experiments revealed that, similar to SAR, the development of a PbBSMT-mediated SAR-like phenotype was also dependent on the MeSA esterase activity of NtSABP2 in the systemic leaves. Collectively, these results strongly suggest that PbBSMT is a novel effector, which is secreted by P. brassicae into its host plant to deplete pathogen-induced SA accumulation.

The immunophilin repertoire of Plasmodiophora brassicae and functional analysis of PbCYP3 cyclophilin.

Plasmodiophora brassicae is a soil-borne pathogen that belongs to Rhizaria, an almost unexplored eukaryotic organism group. This pathogen requires a living host for growth and multiplication, which makes molecular analysis further complicated. To broaden our understanding of a plasmodiophorid such as P. brassicae, we here chose to study immunophilins, a group of proteins known to have various cellular functions, including involvement in plant defense and pathogen virulence. Searches in the P. brassicae genome resulted in 20 putative immunophilins comprising of 11 cyclophilins (CYPs), 7 FK506-binding proteins (FKBPs) and 2 parvulin-like proteins. RNAseq data showed that immunophilins were differentially regulated in enriched life stages such as germinating spores, maturing spores, and plasmodia, and infected Brassica hosts (B. rapa, B. napus and B. oleracea). PbCYP3 was highly induced in all studied life stages and during infection of all three Brassica hosts, and hence was selected for further analysis. PbCYP3 was heterologously expressed in Magnaporthe oryzae gene-inactivated ΔCyp1 strain. The new strain ΔCyp1+ overexpressing PbCYP3 showed increased virulence on rice compared to the ΔCyp1 strain. These results suggest that the predicted immunophilins and particularly PbCYP3 are activated during plant infection. M. oryzae is a well-studied fungal pathogen and could be a valuable tool for future functional studies of P. brassicae genes, particularly elucidating their role during various infection phases.

Integrated omics study of lipid droplets from Plasmodiophora brassicae.

Plasmodiophora brassicae causes clubroot disease in cruciferous. In this report, lipid droplets were observed in the resting spores of P. brassicae. 295 lipid droplet-associated proteins were identified and categorized into nine groups. Transcriptome analysis of these proteins during three different zoosporic stages revealed differences in gene expression pattern. GO enrichment analysis revealed that these proteins associated with lipid droplets were mainly linked to biosynthesis and metabolism. GC-MS analysis revealed that lipid droplets contain seven types of free fatty acids: saturated fatty acids C16:0 and C18:0, and unsaturated fatty acids C18:1Δ9, C18:1Δ11, C18:2, C20:4 and C20:5. P. brassicae accumulated a large amount of triacylglycerols (TAGs). We systematically analyzed the putative proteins involved in TAG biosynthesis and its metabolic pathway. KEGG pathway analysis defined 3390 genes, including 167 genes involved in lipid metabolism. Transcriptome analysis revealed that 162 candidate enzymes involved in lipid metabolism were differential expressed. Our omics studies are the first to investigate the lipid droplet organelles in P. brassicae, providing a reference resource to study protist lipid droplets.