Soil microbiota influences clubroot disease by modulating Plasmodiophora brassicae and Brassica napus transcriptomes

Distinct prokaryotic and eukaryotic communities and networks in two agricultural fields of central Japan with different histories of maize-cabbage rotation

Crop rotation is an important agricultural practice for homeostatic crop cultivation. Here, we applied high-throughput sequencing of ribosomal RNA gene amplicons to investigate soil biota in two fields of central Japan with different histories of maize-cabbage rotation. We identified 3086 eukaryotic and 17,069 prokaryotic sequence variants (SVs) from soil samples from two fields rotating two crops at three different growth stages. The eukaryotic and prokaryotic communities in the four sample groups of two crops and two fields were clearly distinguished using β-diversity analysis. Redundancy analysis showed the relationships of the communities in the fields to pH and nutrient, humus, and/or water content. The complexity of eukaryotic and prokaryotic networks was apparently higher in the cabbage-cultivated soils than those in the maize-cultivated soils. The node SVs (nSVs) of the networks were mainly derived from two eukaryotic phyla: Ascomycota and Cercozoa, and four prokaryotic phyla: Pseudomonadota, Acidobacteriota, Actinomycetota, and Gemmatimonadota. The networks were complexed by cropping from maize to cabbage, suggesting the formation of a flexible network under crop rotation. Ten out of the 16 eukaryotic nSVs were specifically found in the cabbage-cultivated soils were derived from protists, indicating the potential contribution of protists to the formation of complex eukaryotic networks.

Fructan exohydrolases (FEHs) are upregulated by salicylic acid together with defense-related genes in non-fructan accumulating plants

The identification of several fructan exohydrolases (FEHs, EC 3.2.1.80) in non-fructan accumulating plants raised the question of their roles. FEHs may be defense-related proteins involved in the interactions with fructan-accumulating microorganisms. Since known defense-related proteins are upregulated by defense-related phytohormones, we tested the hypothesis that FEHs of non-fructan accumulating plants are upregulated by salicylic acid (SA), jasmonic acid (JA) and ethylene (ET) using the model plant Arabidopsis thaliana and the agronomically relevant and genetically related species Brassica napus. By sequence homologies with the two known FEH genes of A. thaliana, At6-FEH, and At6&1-FEH, the genes coding for the putative B. napus FEHs, Bn6-FEH and Bn6&1-FEH, were identified. Plants were treated at root level with SA, methyl jasmonate (MeJA) or 1-aminocyclopropane-1-carboxylic acid (ACC). The transcript levels of defense-related and FEH genes were measured after treatments. MeJA and ACC did not upregulate FEHs, while HEL (HEVEIN-LIKE PREPROTEIN) expression was enhanced by both phytohormones. In both species, the expression of AOS, encoding a JA biosynthesis enzyme, was enhanced by MeJA and that of the defensine PDF1.2 and the ET signaling transcription factor ERF1/2 by ACC. In contrast, SA not only increased the expression of genes encoding antimicrobial proteins (PR1 and HEL) and the defense-related transcription factor WRKY70 but also that of FEH genes, in particular 6&1-FEH genes. This result supports the putative role of FEHs as defense-related proteins. Genotypic variability of SA-mediated FEH regulation (transcript level and activities) was observed among five varieties of B. napus, suggesting different susceptibilities toward fructan-accumulating pathogens.

The soil bacterial community regulates germination of Plasmodiophora brassicae resting spores rather than root exudates

Clubroot, caused by Plasmodiophora brassicae, is a severe soil-borne disease that restricts the production of cruciferous crops worldwide. A better understanding of biotic and abiotic factors regulating germination of P. brassicae resting spores in the soil is significant for developing novel control methods. Previous studies reported that root exudates can trigger P. brassicae resting spore germination, thus enabling a targeted attack of P. brassicae on host plant roots. However, we found that native root exudates collected under sterile conditions from host or non-host plants cannot stimulate the germination of sterile spores, indicating that root exudates may not be direct stimulation factors. Instead, our studies demonstrate that soil bacteria are essential for triggering germination. Through 16s rRNA amplicon sequencing analysis, we found that certain carbon sources and nitrate can reshape the initial microbial community to an inducing community leading to the germination of P. brassicae resting spores. The stimulating communities significantly differed in composition and abundance of bacterial taxa compared to the non-stimulating ones. Several enriched bacterial taxa in stimulating community were significantly correlated with spore germination rates and may be involved as stimulation factors. Based on our findings, a multi-factorial 'pathobiome' model comprising abiotic and biotic factors is proposed to represent the putative plant-microbiome-pathogen interactions associated with breaking spore dormancy of P. brassicae in soil. This study presents novel views on P. brassicae pathogenicity and lays the foundation for novel sustainable control strategies of clubroot.

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.

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.

Perspectives for integrated insect pest protection in oilseed rape breeding

In the past, breeding for incorporation of insect pest resistance or tolerance into cultivars for use in integrated pest management schemes in oilseed rape/canola (Brassica napus) production has hardly ever been approached. This has been largely due to the broad availability of insecticides and the complexity of dealing with high-throughput phenotyping of insect performance and plant damage parameters. However, recent changes in the political framework in many countries demand future sustainable crop protection which makes breeding approaches for crop protection as a measure for pest insect control attractive again. At the same time, new camera-based tracking technologies, new knowledge-based genomic technologies and new scientific insights into the ecology of insect-Brassica interactions are becoming available. Here we discuss and prioritise promising breeding strategies and direct and indirect breeding targets, and their time-perspective for future realisation in integrated insect pest protection of oilseed rape. In conclusion, researchers and oilseed rape breeders can nowadays benefit from an array of new technologies which in combination will accelerate the development of improved oilseed rape cultivars with multiple insect pest resistances/tolerances in the near future.

Infection by a eukaryotic gut parasite in wild Daphnia sp. associates with a distinct bacterial community

Host-associated bacterial communities play an important role in host fitness and resistance to diseases. Yet, few studies have investigated tripartite interaction between a host, parasite and host-associated bacterial communities in natural settings. Here, we use 16S rRNA gene amplicon sequencing to compare gut- and body- bacterial communities of wild water fleas belonging to the Daphnia longispina complex, between uninfected hosts and those infected with the common and virulent eukaryotic gut parasite Caullerya mesnili (Family: Ichthyosporea). We report community-level changes in host-associated bacteria with the presence of the parasite infection; namely decreased alpha diversity and increased beta diversity at the site of infection, i.e. host gut (but not host body). We also report decreased abundance of bacterial taxa proposed elsewhere to be beneficial for the host, and an appearance of taxa specifically associated with infected hosts. Our study highlights the host-microbiota-infection link in a natural system and raises questions about the role of host-associated microbiota in natural disease epidemics as well as the functional roles of bacteria specifically associated with infected hosts.

Response of Bacterial Community to the Occurrence of Clubroot Disease in Chinese Cabbage

Clubroot disease is a common soilborne disease caused by Plasmodiophora brassicas Wor. and widely occurs in Chinese cabbage. Soil microorganisms play vital roles in the occurrence and development of plant diseases. The changes in the soil bacterial community could indicate the severity of plant disease and provide the basis for its control. This study focused on the bacterial community of the clubroot disease-infected soil-root system with different severity aiming to reveal the composition and structure of soil bacteria and identified potential biomarker bacteria of the clubroot disease. In the clubroot disease-infected soil, the bacterial community is mainly composed of Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, Bacilli, Thermolrophilia, Bacteroidia, Gemmatimonadetes, Subgroup_6, Deltaproteobacteria, KD4-96, and some other classes, while the major bacterial classes in the infected roots were Oxyphotobacteria, Gammaproteobacteria, Alphaproteobacteria, Actinobacteria, Bacilli, Bacteroidia, Saccharimonadia, Thermoleophilia, Clostridia, Chloroflexia, and some other classes. The severe clubroot disease soil-root system was found to possess a poorer bacterial richness, evenness, and better coverage. Additionally, a significant difference was observed in the structure of the bacterial community between the high-severity (HR) and healthy (LR) soil-root system. Bacillus asahii and Noccaea caerulescens were identified as the differential bacteria between the LR and HR soil and roots, respectively. pH was demonstrated as a vital factor that was significantly associated with the abundance of B. asahii and N. caerulescens. This study provides novel insight into the relationship between soil bacteria and the pathogen of clubroot disease in Chinese cabbage. The identification of resistant species provides candidates for the monitoring and biocontrol of the clubroot disease.

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.

Nitrogen Supply and Host-Plant Genotype Modulate the Transcriptomic Profile of Plasmodiophora brassicae

Nitrogen fertilization can affect the susceptibility of Brassica napus to the telluric pathogen Plasmodiophora brassicae. Our previous works highlighted that the influence of nitrogen can strongly vary regarding plant cultivar/pathogen strain combinations, but the underlying mechanisms are unknown. The present work aims to explore how nitrogen supply can affect the molecular physiology of P. brassicae through its life epidemiological cycle. A time-course transcriptome experiment was conducted to study the interaction, under two conditions of nitrogen supply, between isolate eH and two B. napus genotypes (Yudal and HD-018), harboring (or not harboring) low nitrogen-conditional resistance toward this isolate (respectively). P. brassicae transcriptional patterns were modulated by nitrogen supply, these modulations being dependent on both host-plant genotype and kinetic time. Functional analysis allowed the identification of P. brassicae genes expressed during the secondary phase of infection, which may play a role in the reduction of Yudal disease symptoms in low-nitrogen conditions. Candidate genes included pathogenicity-related genes ("NUDIX," "carboxypeptidase," and "NEP-proteins") and genes associated to obligate biotrophic functions of P. brassicae. This work illustrates the importance of considering pathogen's physiological responses to get a better understanding of the influence of abiotic factors on clubroot resistance/susceptibility.

Control of Streptomyces alfalfae XY25 T Over Clubroot Disease and Its Effect on Rhizosphere Microbial Community in Chinese Cabbage Field Trials

Clubroot caused by Plasmodiophora brassicae is one of the most destructive diseases in cruciferous crops. Streptomyces alfalfae XY25 ^T , a biological control agent, exhibited great ability to relieve clubroot disease, regulate rhizosphere bacterial and fungal communities in Chinese cabbage, and promote its growth in greenhouse. Therefore, field experiments were carried out to investigate the effects of S. alfalfae XY25 ^T on clubroot and rhizosphere microbial community in Chinese cabbage. Results showed that the control efficiency of clubroot by S. alfalfae XY25 ^T was 69.4%. Applying the agent can alleviate soil acidification; increase the contents of soil organic matter, available nitrogen, available phosphorus, and available potassium; and enhance activities of invertase, urease, catalase, and alkaline phosphatase. During Chinese cabbage growth, bacterial diversity decreased first and then increased, and fungal diversity decreased gradually after inoculation with S. alfalfae XY25 ^T . High-throughput sequencing analysis showed that the main bacterial phyla were Proteobacteria, Bacteroidetes, Acidobacteria, and Planctomycetes, and the major fungal phyla were Ascomycota and Basidiomycota in rhizosphere soil. The dominant bacterial genera were Flavobacterium, Candidatus, Pseudomonas, Stenotrophomonas, Sphingomonas, Flavisolibacter, and Gemmatimonbacteria with no significant difference in abundance, and the major fungal genera were Monographella, Aspergillus, Hypocreales, Chytridiaceae, Fusarium, Pleosporales, Agaricales, Mortierella, and Pleosporales. The significant differences were observed among Pleosporales, Basidiomycota, Colletotrichum, two strains attributed to Agaricales, and another two unidentified fungi by using S. alfalfae XY25 ^T . Moreover, quantitative real-time PCR results indicated that P. brassicae content was significantly decreased after the agent inoculation. In conclusion, S. alfalfae XY25 ^T can affect rhizosphere microbial communities; therefore, applying the agent is an effective approach to reduce the damage caused by clubroot.

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.