Demystifying biotrophs: FISHing for mRNAs to decipher plant and algal pathogen-host interaction at the single cell level

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.

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.

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.

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 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.

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.

Visualization of Gene Reciprocity among Lactic Acid Bacteria in Yogurt by RNase H-Assisted Rolling Circle Amplification-Fluorescence In Situ Hybridization

Recently, we developed an in situ mRNA detection method termed RNase H-assisted rolling circle amplification-fluorescence in situ hybridization (RHa-RCA-FISH), which can detect even short mRNA in a bacterial cell. However, because this FISH method is sensitive to the sample condition, it is necessary to find a suitable cell permeabilization and collection protocol. Here, we demonstrate its further applicability for detecting intrinsic mRNA expression using lactic acid bacteria (LAB) as a model consortium. Our results show that this method can visualize functional gene expression in LAB cells and can be used for monitoring the temporal transition of gene expression. In addition, we also confirmed that data obtained from bulk analyses such as RNA-seq or microarray do not always correspond to gene expression in individual cells. RHa-RCA-FISH will be a powerful tool to compensate for insufficient data from metatranscriptome analyses while clarifying the carriers of function in microbial consortia. By extending this technique to capture spatiotemporal microbial gene expression at the single-cell level, it will be able to characterize microbial interactions in phytoplankton-bacteria interactions.

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.