The architecture of the Plasmodiophora brassicae nuclear and mitochondrial genomes

Decoding the Arsenal: Protist Effectors and Their Impact on Photosynthetic Hosts

Interactions between various microbial pathogens including viruses, bacteria, fungi, oomycetes, and their plant hosts have traditionally been the focus of phytopathology. In recent years, a significant and growing interest in the study of eukaryotic microorganisms not classified among fungi or oomycetes has emerged. Many of these protists establish complex interactions with photosynthetic hosts, and understanding these interactions is crucial in understanding the dynamics of these parasites within traditional and emerging types of farming, including marine aquaculture. Many phytopathogenic protists are biotrophs with complex polyphasic life cycles, which makes them difficult or impossible to culture, a fact reflected in a wide gap in the availability of comprehensive genomic data when compared to fungal and oomycete plant pathogens. Furthermore, our ability to use available genomic resources for these protists is limited by the broad taxonomic distance that these organisms span, which makes comparisons with other genomic datasets difficult. The current rapid progress in genomics and computational tools for the prediction of protein functions and interactions is revolutionizing the landscape in plant pathology. This is also opening novel possibilities, specifically for a deeper understanding of protist effectors. Tools like AlphaFold2 enable structure-based function prediction of effector candidates with divergent protein sequences. In turn, this allows us to ask better biological questions and, coupled with innovative experimental strategies, will lead into a new era of effector research, especially for protists, to expand our knowledge on these elusive pathogens and their interactions with photosynthetic hosts. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Telomere-to-telomere Genome Assembly of the Clubroot Pathogen Plasmodiophora Brassicae

Plasmodiophora brassicae (Woronin, 1877), a biotrophic, obligate parasite, is the causal agent of clubroot disease in brassicas. The clubroot pathogen has been reported in more than 80 countries worldwide, causing economic losses of hundreds of millions every year. Despite its widespread impact, very little is known about the molecular strategies it employs to induce the characteristic clubs in the roots of susceptible hosts during infection, nor about the mechanisms it uses to overcome genetic resistance. Here, we provide the first telomere-to-telomere complete genome of P. brassicae. We generated ∼27 Gb of Illumina, Oxford Nanopore, and PacBio HiFi data from resting spores of strain Pb3A and produced a 25.3 Mb assembly comprising 20 chromosomes, with an N50 of 1.37 Mb. The BUSCO score, the highest reported for any member of the group Rhizaria (Eukaryota: 88.2%), highlights the limitations within the Eukaryota database for members of this lineage. Using available transcriptomic data and protein evidence, we annotated the Pb3A genome, identifying 10,521 protein-coding gene models. This high-quality, complete genome of P. brassicae will serve as a crucial resource for the plant pathology community to advance the much-needed understanding of the evolution of the clubroot pathogen.

Comparative genomics of Ascetosporea gives new insight into the evolutionary basis for animal parasitism in Rhizaria

BACKGROUND

Ascetosporea (Endomyxa, Rhizaria) is a group of unicellular parasites infecting aquatic invertebrates. They are increasingly being recognized as widespread and important in marine environments, causing large annual losses in invertebrate aquaculture. Despite their importance, little molecular data of Ascetosporea exist, with only two genome assemblies published to date. Accordingly, the evolutionary origin of these parasites is unclear, including their phylogenetic position and the genomic adaptations that accompanied the transition from a free-living lifestyle to parasitism. Here, we sequenced and assembled three new ascetosporean genomes, as well as the genome of a closely related amphizoic species, to investigate the phylogeny, origin, and genomic adaptations to parasitism in Ascetosporea.

RESULTS

Using a phylogenomic approach, we confirm the monophyly of Ascetosporea and show that Paramyxida group with Mikrocytida, with Haplosporida being sister to both groups. We report that the genomes of these parasites are relatively small (12-36 Mb) and gene-sparse (~ 2300-5200 genes), while containing surprisingly high amounts of non-coding sequence (~ 70-90% of the genomes). Performing gene-tree aware ancestral reconstruction of gene families, we demonstrate extensive gene losses at the origin of parasitism in Ascetosporea, primarily of metabolic functions, and little gene gain except on terminal branches. Finally, we highlight some functional gene classes that have undergone expansions during evolution of the group.

CONCLUSIONS

We present important new genomic information from a lineage of enigmatic but important parasites of invertebrates and illuminate some of the genomic innovations accompanying the evolutionary transition to parasitism in this lineage. Our results and data provide a genetic basis for the development of control measures against these parasites.

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.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

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.

Protocol: rhPCR and SNaPshot assays to distinguish Plasmodiophora brassicae pathotype clusters

BACKGROUND

Clubroot of canola (Brassica napus), caused by the soilborne pathogen Plasmodiophora brassicae, has become a serious threat to canola production in Canada. The deployment of clubroot-resistant (CR) cultivars is the most commonly used management strategy; however, the widespread cultivation of CR canola has resulted in the emergence of new pathotypes of P. brassicae capable of overcoming resistance. Several host differential sets have been reported for pathotype identification, but such testing is time-consuming, labor-intensive, and based on phenotypic classifications. The development of rapid and objective methods that allow for efficient, cost-effective and convenient pathotyping would enable testing of a much larger number of samples in shorter times. The aim of this study was to develop two pathotyping assays, an RNase H2-dependent PCR (rhPCR) assay and a SNaPshot assay, which could quickly differentiate P. brassicae pathotypes.

RESULTS

Both assays clearly distinguished between pathotype clusters in a collection of 38 single-spore isolates of P. brassicae. Additional isolates pathotyped from clubbed roots and samples from blind testing also were correctly clustered. The rhPCR assay generated clearly differentiating electrophoretic bands without non-specific amplification. The SNaPshot assay was able to detect down to a 10% relative allelic proportion in a 10:90 template mixture with both single-spore isolates and field isolates when evaluated in a relative abundance test.

CONCLUSIONS

This study describes the development of two rapid and sensitive technologies for P. brassicae pathotyping. The high-throughput potential and accuracy of both assays makes them promising as SNP-based pathotype identification tools for clubroot diagnostics. rhPCR is a highly sensitive approach that can be optimized into a quantitative assay, while the main advantages of SNaPshot are its ability to multiplex samples and alleles in a single reaction and the detection of up to four allelic variants per target site.

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.

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.

Metagenomics approach for Polymyxa betae genome assembly enables comparative analysis towards deciphering the intracellular parasitic lifestyle of the plasmodiophorids

Genomic knowledge of the tree of life is biased to specific groups of organisms. For example, only six full genomes are currently available in the rhizaria clade. Here, we have applied metagenomic techniques enabling the assembly of the genome of Polymyxa betae (Rhizaria, Plasmodiophorida) RES F41 isolate from unpurified zoospore holobiont and comparison with the A26-41 isolate. Furthermore, the first P. betae mitochondrial genome was assembled. The two P. betae nuclear genomes were highly similar, each with just ~10.2 k predicted protein coding genes, ~3% of which were unique to each isolate. Extending genomic comparisons revealed a greater overlap with Spongospora subterranea than with Plasmodiophora brassicae, including orthologs of the mammalian cation channel sperm-associated proteins, raising some intriguing questions about zoospore physiology. This work validates our metagenomics pipeline for eukaryote genome assembly from unpurified samples and enriches plasmodiophorid genomics; providing the first full annotation of the P. betae genome.

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.

Genetic Structure of Plasmodiophora brassicae Populations Virulent on Clubroot Resistant Canola (Brassica napus)

Clubroot, caused by Plasmodiophora brassicae Woronin, is a significant threat to the canola (Brassica napus L.) industry in Canada. Clubroot resistance has been overcome in more than 200 fields since 2013, representing one of the biggest challenges to sustainable canola production. The genetic structure of 36 single-spore isolates derived from 12 field isolates of P. brassicae collected before and after the introduction of clubroot resistant (CR) canola cultivars (2005-2014) was evaluated by simple sequence repeat (SSR) marker analysis. Polymorphisms were detected in 32 loci with the identification of 93 distinct alleles. A low level of genetic diversity was found among the single-spore isolates. Haploid linkage disequilibrium and number of migrants suggested that recombination and migration were rare or almost absent in the tested P. brassicae population. A relatively clear relationship was found between the genetic structure and virulence phenotypes of the pathogen as defined on the differential hosts of Somé et al., Williams, and the Canadian Clubroot Differential (CCD) set. Although genetic variability within each pathotype group, as classified on each differential system, was low, significant genetic differentiation was observed among the pathotypes. The highest correlation between genetic structure and virulence was found among matrices produced with genetic data and the hosts of the CCD set, with a threshold index of disease of 50% to distinguish susceptible from resistant reactions. Genetically homogeneous single-spore isolates provided a more complete and clearer picture of the population genetic structure of P. brassicae, and the results suggest some promise for the development of pathotype-specific primers.

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.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

Recent advances in understanding mitochondrial genome diversity

Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.