Genome sequence of the necrotrophic plant pathogen Pythium ultimum reveals original pathogenicity mechanisms and effector repertoire

Sensitivity of Globisporangium ultimum to the fungicide metalaxyl is enhanced by the infection with a toti-like mycovirus

In recent years, numerous oomycete mycoviruses have been discovered; however, very few studies have focused on their effects on the host oomycete phenotype. In this study, we investigated the impact of toti-like Pythium ultimum RNA virus 2 (PuRV2) infection on the phytopathogenic soil-borne oomycete Globisporangium ultimum, which serves as a model species for Globisporangium and Pythium, specifically the UOP226 isolate in Japan. We generated a PuRV2-free isogenic line through hyphal tip isolation using high-temperature culture and subsequently compared the phenotypic characteristics and gene expression profiles of UOP226 and the PuRV2-free isogenic line. Our findings revealed that the metalaxyl sensitivity of UOP226 was greater than that of the PuRV2-free isogenic line, whereas the mycelial growth rate and colony morphology remained unchanged in the absence of the fungicide. Furthermore, transcriptome analyses using RNA-seq revealed significant downregulation of ABC-type transporter genes, which are involved in fungicide sensitivity, in UOP226. Our results suggest that PuRV2 infection influences the ecology of G. ultimum in agricultural ecosystems where metalaxyl is applied.

Duplication and neofunctionalization of a horizontally transferred xyloglucanase as a facet of the Red Queen coevolutionary dynamic

Oomycete protists share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a distant region of the tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls, a barrier to pathogen invasion and a rich source of carbohydrates. Using a combination of phylogenomics and functional assays, we investigate the diversification of a horizontally transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 xyloglucanase paralogs retained in P. sojae. Using heterologous expression in yeast, we show consistent evidence that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional variants analyzed subtend a phylogenetic node close to the fungi-to-oomycete transfer, suggesting the horizontally transferred gene was a bona fide xyloglucanase. Expression of three xyloglucanase paralogs in Nicotiana benthamiana triggers high-reactive oxygen species (ROS) generation, while others inhibit ROS responses to bacterial immunogens, demonstrating that the paralogs differentially stimulate pattern-triggered immunity. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze the production of variant breakdown profiles, suggesting that secretion of variant xyloglucanases increases efficiency of xyloglucan breakdown as well as diversifying the damage-associated molecular patterns released. We suggest that this pattern of neofunctionalization and the variant host responses represent an aspect of the Red Queen host-pathogen coevolutionary dynamic.

Complete telomere-to-telomere genomes uncover virulence evolution conferred by chromosome fusion in oomycete plant pathogens

Variations in chromosome number are occasionally observed among oomycetes, a group that includes many plant pathogens, but the emergence of such variations and their effects on genome and virulence evolution remain ambiguous. We generated complete telomere-to-telomere genome assemblies for Phytophthora sojae, Globisporangium ultimum, Pythium oligandrum, and G. spinosum. Reconstructing the karyotype of the most recent common ancestor in Peronosporales revealed that frequent chromosome fusion and fission drove changes in chromosome number. Centromeres enriched with Copia-like transposons may contribute to chromosome fusion and fission events. Chromosome fusion facilitated the emergence of pathogenicity genes and their adaptive evolution. Effectors tended to duplicate in the sub-telomere regions of fused chromosomes, which exhibited evolutionary features distinct to the non-fused chromosomes. By integrating ancestral genomic dynamics and structural predictions, we have identified secreted Ankyrin repeat-containing proteins (ANKs) as a novel class of effectors in P. sojae. Phylogenetic analysis and experiments further revealed that ANK is a specifically expanded effector family in oomycetes. These results revealed chromosome dynamics in oomycete plant pathogens, and provided novel insights into karyotype and effector evolution.

Advances in understanding grapevine downy mildew: From pathogen infection to disease management

Plasmopara viticola is geographically widespread in grapevine-growing regions. Grapevine downy mildew disease, caused by this biotrophic pathogen, leads to considerable yield losses in viticulture annually. Because of the great significance of grapevine production and wine quality, research on this disease has been widely performed since its emergence in the 19th century. Here, we review and discuss recent understanding of this pathogen from multiple aspects, including its infection cycle, disease symptoms, genome decoding, effector biology, and management and control strategies. We highlight the identification and characterization of effector proteins with their biological roles in host-pathogen interaction, with a focus on sustainable control methods against P. viticola, especially the use of biocontrol agents and environmentally friendly compounds.

Plant Growth Promotion and Stress Tolerance Enhancement through Inoculation with Bacillus proteolyticus OSUB18

The isolation of B. proteolyticus OSUB18 from switchgrass unveiled its significant potential in both the enhancement of plant growth and the suppression of plant diseases in our previous study. The elucidation of the related mechanisms governing this intricate plant-microbe interaction involved the utilization of the model plant Arabidopsis thaliana. In our comprehensive study on Arabidopsis, OSUB18 treatment was found to significantly alter root architecture and enhance plant growth under various abiotic stresses. An RNA-seq analysis revealed that OSUB18 modified gene expression, notably upregulating the genes involved in glucosinolate biosynthesis and plant defense, while downregulating those related to flavonoid biosynthesis and wound response. Importantly, OSUB18 also induces systemic resistance in Arabidopsis against a spectrum of bacterial and fungal pathogens and exhibits antagonistic effects on phytopathogenic bacteria, fungi, and oomycetes, highlighting its potential as a beneficial agent in plant stress management and pathogen resistance. Overall, our findings substantiate that OSUB18 exerts a stimulatory influence on plant growth and health, potentially attributed to the remodeling of root architecture, defense signaling, and the comprehensive mitigation of various biotic and abiotic stresses.

Resolving the Globisporangium ultimum (Pythium ultimum) species complex

The Globisporangium ultimum (formerly Pythium ultimum) species complex was previously composed of two morphological varieties: var. ultimum and var. sporangiiferum. Prior attempts to resolve this morphology-based species complex using molecular techniques have been inconclusive or conflicting. The increased availability of sequenced genomes and isolates identified as G. ultimum var. ultimum and var. sporangiiferum has allowed us to examine these relationships at a higher resolution and with a broader scope than previously possible. Using comparative genomics, we identified highly variable gene regions and designed primers for four new protein-coding genes for phylogenetics. These were then used alongside three known markers to generate a nuclear multigene genealogy of the species complex. From a collection of 163 isolates belonging to the target taxa, a subset of 29 was chosen to be included in this study (verified with nuclear rDNA internal transcribed spacer 1 [ITS1] and mitochondrial cytochrome c oxidase subunit 1 [cox1] sequences). Seventeen isolates of var. ultimum were selected to be representative of variations in genotype, morphology, and geographic collection location. The 12 isolates of var. sporangiiferum included all available specimens identified either morphologically (in previous studies) or through sequence similarity with ITS1 and cox1. Based on the fulfillment of reciprocal monophyly and observed genealogical concordance under the genealogical concordance phylogenetic species recognition, we determined that the Globisporangium ultimum species complex is composed of four genetically distinct species: Globisporangium ultimum, Globisporangium sporangiiferum, Globisporangium solveigiae, and Globisporangium bothae.

Cell-type-specific responses to fungal infection in plants revealed by single-cell transcriptomics

Pathogen infection is a dynamic process. Here, we employ single-cell transcriptomics to investigate plant response heterogeneity. By generating an Arabidopsis thaliana leaf atlas encompassing 95,040 cells during infection by a fungal pathogen, Colletotrichum higginsianum, we unveil cell-type-specific gene expression, notably an enrichment of intracellular immune receptors in vasculature cells. Trajectory inference identifies cells that had different interactions with the invading fungus. This analysis divulges transcriptional reprogramming of abscisic acid signaling specifically occurring in guard cells, which is consistent with a stomatal closure dependent on direct contact with the fungus. Furthermore, we investigate the transcriptional plasticity of genes involved in glucosinolate biosynthesis in cells at the fungal infection sites, emphasizing the contribution of the epidermis-expressed MYB122 to disease resistance. This work underscores spatially dynamic, cell-type-specific plant responses to a fungal pathogen and provides a valuable resource that supports in-depth investigations of plant-pathogen interactions.

Duplication and neofunctionalization of a horizontally-transferred xyloglucanase as a facet of the red queen co-evolutionary dynamic

Oomycetes are heterotrophic protists that share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a separate and distant region of the eukaryotic tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls. This is a key trait in the pathology of many oomycetes, as the plant cell wall represents a primary barrier to pathogen invasion and a rich source of carbohydrates. Many of the HGT gene families identified have undergone multiple rounds of duplication. Using a combination of phylogenomic analysis and functional assays, we investigate the diversification of a horizontally-transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 genes retained in P. sojae among a complex pattern of gene duplications and losses. Using a phenotype assay, based on heterologous expression in yeast, we show that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional xyloglucanase variants analysed subtend an ancestral node close to the fungi-oomycetes gene transfer, suggesting the horizontally-transferred gene was a bona fide xyloglucanase. Expression of xyloglucanase paralogs in Nicotiana benthamiana triggers distinct patterns of reactive oxygen species (ROS) generation, demonstrating that enzyme variants differentially stimulate pattern-triggered immunity in plants. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze production of variant breakdown profiles, suggesting that secretion of multiple xyloglucanase variants increases efficiency of xyloglucan breakdown, as well as potentially diversifying the range of Damage-Associated Molecular Patterns (DAMPs) released during pathogen attack. We suggest that such patterns of protein neofunctionalization, and variant host responses, represent an aspect of the Red Queen host-pathogen co-evolutionary dynamic.

Significance Statement

The oomycetes are a diverse group of eukaryotic microbes that include some of the most devastating pathogens of plants. Oomycetes perceive, invade, and colonize plants in similar ways to fungi, in part because they acquired the genes to attack and feed on plants from fungi. These genes are predicted to be useful to oomycete plant pathogens because they have undergone multiple rounds of gene duplication. One key enzyme for attacking plant cell wall structures is called xyloglucanase. Xyloglucanase in the oomycetes has undergone multiple rounds of gene duplication, leading to variants including an enzyme with a C-terminal extension that increases activity. Some xyloglucanase variants trigger unique patterns of reactive oxygen species (ROS) in planta, and generate different profiles of cell wall breakdown products - such outcomes could act to mystify and increase the workload of the plant immune system, allowing successful pathogens to proliferate.

Comparative genomics reveals the presence of simple sequence repeats in genes related to virulence in plant pathogenic Pythium ultimum and Pythium vexans

In this study, we evaluated the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) in the complete genome and transcriptomic sequences of the plant pathogenic species of Pythium to acquire a better knowledge of their genome structure and evolution. Among the species, P. ultimum had the highest RA and RD of SSRs in the genomic sequences, whereas P. vexans had the highest RA and RD in the transcriptomic sequences. The genomic and transcriptomic sequences of P. aphanidermatum showed the lowest RA and RD of SSRs. Trinucleotide SSRs were the most prevalent class in both genomic and transcriptomic sequences, while dinucleotide SSRs were the least prevalent. The G + C content of the transcriptomic sequences was found to be positively correlated with the number (r = 0.601) and RA (r = 0.710) of SSRs. A motif conservation study revealed the highest number of unique motifs in P. vexans (9.9%). Overall, a low conservation of motifs was observed among the species (25.9%). A gene enrichment study revealed P. vexans and P. ultimum carry SSRs in their genes that are directly connected to virulence, whereas the remaining two species, P. aphanidermatum and P. arrhenomanes, harbour SSRs in genes involved in transcription, translation, and ATP binding. In an effort to enhance the genomic resources, a total of 11,002 primers from the transcribed regions were designed for the pathogenic Pythium species. Furthermore, the unique motifs identified in this work could be employed as molecular probes for species identification.

Comparative Genomic Analysis Reveals Gene Content Diversity, Phylogenomic Contour, Putative Virulence Determinants, and Potential Diagnostic Markers within Pythium insidiosum Traits

Pythium insidiosum has successfully evolved into a human/animal filamentous pathogen, causing pythiosis, a life-threatening disease, worldwide. The specific rDNA-based genotype of P. insidiosum (clade I, II, or III) is associated with the different hosts and disease prevalence. Genome evolution of P. insidiosum can be driven by point mutations, pass vertically to the offspring, and diverge into distinct lineages, leading to different virulence, including the ability to be unrecognized by the host. We conducted comprehensive genomic comparisons of 10 P. insidiosum strains and 5 related Pythium species using our online "Gene Table" software to investigate the pathogen's evolutionary history and pathogenicity. In total, 245,378 genes were found in all 15 genomes and grouped into 45,801 homologous gene clusters. Gene contents among P. insidiosum strains varied by as much as 23%. Our results showed a strong agreement between the phylogenetic analysis of 166 core genes (88,017 bp) identified across all genomes and the hierarchical clustering analysis of gene presence/absence profiles, suggesting divergence of P. insidiosum into two groups, clade I/II and clade III strains, and the subsequent segregation of clade I and clade II. A stringent gene content comparison using the Pythium Gene Table provided 3263 core genes exclusively presented in all P. insidiosum strains but no other Pythium species, which could involve host-specific pathogenesis and serve as biomarkers for diagnostic purposes. More studies focusing on characterizing the biological function of the core genes (including the just-identified putative virulence genes encoding hemagglutinin/adhesin and reticulocyte-binding protein) are needed to explore the biology and pathogenicity of this pathogen.

Cell polarity in the protist-to-animal transition

A signature feature of the animal kingdom is the presence of epithelia: sheets of polarized cells that both insulate the organism from its environment and mediate interactions with it. Epithelial cells display a marked apico-basal polarity, which is highly conserved across the animal kingdom, both in terms of morphology and of molecular regulators. How did this architecture first evolve? Although the last eukaryotic common ancestor almost certainly possessed a simple form of apico-basal polarity (marked by the presence of one or several flagella at a single cellular pole), comparative genomics and evolutionary cell biology reveal that the polarity regulators of animal epithelial cells have a surprisingly complex and stepwise evolutionary history. Here, we retrace their evolutionary assembly. We suggest that the "polarity network" that polarized animal epithelial cells evolved by integration of initially independent cellular modules that evolved at distinct steps of our evolutionary ancestry. The first module dates back to the last common ancestor of animals and amoebozoans and involved Par1, extracellular matrix proteins, and the integrin-mediated adhesion complex. Other regulators, such as Cdc42, Dlg, Par6 and cadherins evolved in ancient unicellular opisthokonts, and might have first been involved in F-actin remodeling and filopodial dynamics. Finally, the bulk of "polarity proteins" as well as specialized adhesion complexes evolved in the metazoan stem-line, in concert with the newly evolved intercellular junctional belts. Thus, the polarized architecture of epithelia can be understood as a palimpsest of components of distinct histories and ancestral functions, which have become tightly integrated in animal tissues.

Complete Genome Sequence of Pythium oligandrum, Isolated from Rhizosphere Soils of Chinese Angelica sinensis

Most of the Pythium species are pathogenic to a wide range of economically important crops and, sometimes, can even cause diseases in animals and humans. An exception is that the soil-inhabiting P. oligandrum is an effective biocontrol agent against a diverse suite of pathogens and promotes plant growth. In this work, we sequenced the whole genome of P. oligandrum PO-1, isolated from rhizosphere soils of Chinese Angelica sinensis, using a combination of long-read single-molecule real-time sequencing technology (Pacific Biosciences [PacBio]) and Illumina sequencing. The 2.5-Gb and 5.2-Gb bases were generated respectively. The sequencing depths were 93× with PacBio and 145× with Illumina sequencing. With the PacBio sequencing results further corrected by Illumina sequencing, the genome was assembled into 71 scaffolds with a total size of 39.10 Mb (N₅₀ = 1.45 Mb; L₅₀ = 9)and the longest scaffold is 3.49 Mb. Genome annotation identifies 15,632 protein-coding genes and 0.47 Mb of transposable elements. Our genomic assembly and annotation have been greatly improved compared with the already released three genomes of P. oligandrum. This genomic data will provide valuable information to understand the mechanism underlying its biocontrol potentials and will also facilitate the dissection of genome evolution and environmental adaptation within the genus Pythium. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2022.

SWATH-MS based quantitative proteomics analysis reveals novel proteins involved in PAMP triggered immunity against potato late blight pathogen Phytophthora infestans

Potato is the most important non-grain food in the world, while late blight caused by Phytophthora infestans seriously threatens the production of potato. Since pathogen-associated molecular patterns (PAMPs) are relatively conserved, PAMP-triggered immunity (PTI) can provide durable resistance to late blight for potato. However, knowledge of the regulatory mechanisms of PTI against oomycete pathogens at protein levels remains limited due to the small number of identified proteins. In the present work, changes in the proteome profile of Nicotiana benthamiana leaves upon P. infestans PAMP induction were examined using the SWATH-MS (sequential windowed acquisition of all theoretical mass spectra) approach, which provides quantification of protein abundances and large-scale identification of PTI-related proteins. A total of 4401 proteins have been identified, of which 1429 proteins were differentially expressed at least at one time point of 8, 12, 24 and 48 h after PAMP induction, compared with the expression at 0 h when immediately after PAMP induction. They were further analyzed by expression clustering and gene ontology (GO) enrichment analysis. Through functional verification, six novel DEPs of 19 candidates were proved to be involved in PTI responses, including mitochondrial phosphate carrier protein (MPT) 3, vesicle-associated membrane protein (VAMP) 714, lysophospholipase (LysoPL) 2, ascorbate peroxidase (APX) 1, heat shock 70 kDa protein (HSP) 2 and peptidyl-prolyl cis-trans isomerase FKBP (FKBP) 15-1. Taken together, the time course approach and the resulting large-scale proteomic analyses have enlarged our understanding of PTI mechanisms and provided a valuable resource for the discovery of complex protein networks involved in the resistance response of potato to late blight.

Identification of a C2H2 Transcription Factor (PsCZF3) Associated with RxLR Effectors and Carbohydrate-Active Enzymes in Phytophthora sojae Based on WGCNA

Phytophthora sojae is a destructive soybean pathogen that orchestrates various secreted proteins (effectors) to modulate plant immunity and facilitate infection. Although a number of effectors have been identified and functionally studied in P. sojae, the way these molecules are regulated is marginally known. In this study, we performed a weighted gene correlation network analysis (WGCNA) based on digital RNA-seq, which enabled the identification of a transcription factor (PsCZF3) in P. sojae. This transcription factor is a C₂H₂-type zinc finger protein that regulates the transcription of 35 RxLR effectors during the early infection stage. Phylogenetic analysis revealed that PsCZF3 is a highly conserved protein across oomycetes, suggesting that this regulation mechanism may broadly exist in oomycete species. In addition, by building a subnetwork of PsCZF3 and correlated genes, we also found that PsCZF3 contributed to the transcriptional regulation of carbohydrate-active enzymes. Our findings suggest that the activation of PsCZF3 facilitates P. sojae infection by up-regulating RxLR effectors and carbohydrate-active enzymes.

Genome of Pythium myriotylum Uncovers an Extensive Arsenal of Virulence-Related Genes among the Broad-Host-Range Necrotrophic Pythium Plant Pathogens

The Pythium (Peronosporales, Oomycota) genus includes devastating plant pathogens that cause widespread diseases and severe crop losses. Here, we have uncovered a far greater arsenal of virulence factor-related genes in the necrotrophic Pythium myriotylum than in other Pythium plant pathogens. The genome of a plant-virulent P. myriotylum strain (~70 Mb and 19,878 genes) isolated from a diseased rhizome of ginger (Zingiber officinale) encodes the largest repertoire of putative effectors, proteases, and plant cell wall-degrading enzymes (PCWDEs) among the studied species. P. myriotylum has twice as many predicted secreted proteins than any other Pythium plant pathogen. Arrays of tandem duplications appear to be a key factor of the enrichment of the virulence factor-related genes in P. myriotylum. The transcriptomic analysis performed on two P. myriotylum isolates infecting ginger leaves showed that proteases were a major part of the upregulated genes along with PCWDEs, Nep1-like proteins (NLPs), and elicitin-like proteins. A subset of P. myriotylum NLPs were analyzed and found to have necrosis-inducing ability from agroinfiltration of tobacco (Nicotiana benthamiana) leaves. One of the heterologously produced infection-upregulated putative cutinases found in a tandem array showed esterase activity with preferences for longer-chain-length substrates and neutral to alkaline pH levels. Our results allow the development of science-based targets for the management of P. myriotylum-caused disease, as insights from the genome and transcriptome show that gene expansion of virulence factor-related genes play a bigger role in the plant parasitism of Pythium spp. than previously thought. IMPORTANCE Pythium species are oomycetes, an evolutionarily distinct group of filamentous fungus-like stramenopiles. The Pythium genus includes several pathogens of important crop species, e.g., the spice ginger. Analysis of our genome from the plant pathogen Pythium myriotylum uncovered a far larger arsenal of virulence factor-related genes than found in other Pythium plant pathogens, and these genes contribute to the infection of the plant host. The increase in the number of virulence factor-related genes appears to have occurred through the mechanism of tandem gene duplication events. Genes from particular virulence factor-related categories that were increased in number and switched on during infection of ginger leaves had their activities tested. These genes have toxic activities toward plant cells or activities to hydrolyze polymeric components of the plant. The research suggests targets to better manage diseases caused by P. myriotylum and prompts renewed attention to the genomics of Pythium plant pathogens.

Evolution of pathogenicity in obligate fungal pathogens and allied genera

Obligate fungal pathogens (ascomycetes and basidiomycetes) and oomycetes are known to cause diseases in cereal crop plants. They feed on living cells and most of them have learned to bypass the host immune machinery. This paper discusses some of the factors that are associated with pathogenicity drawing examples from ascomycetes, basidiomycetes and oomycetes, with respect to their manifestation in crop plants. The comparisons have revealed a striking similarity in the three groups suggesting convergent pathways that have arisen from three lineages independently leading to an obligate lifestyle. This review has been written with the intent, that new information on adaptation strategies of biotrophs, modifications in pathogenicity strategies and population dynamics will improve current strategies for breeding with stable resistance.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

Unraveling Plant Cell Death during Phytophthora Infection

Oomycetes form a distinct phylogenetic lineage of fungus-like eukaryotic microorganisms, of which several hundred organisms are considered among the most devastating plant pathogens—especially members of the genus Phytophthora. Phytophthora spp. have a large repertoire of effectors that aid in eliciting a susceptible response in host plants. What is of increasing interest is the involvement of Phytophthora effectors in regulating programed cell death (PCD)—in particular, the hypersensitive response. There have been numerous functional characterization studies, which demonstrate Phytophthora effectors either inducing or suppressing host cell death, which may play a crucial role in Phytophthora’s ability to regulate their hemi-biotrophic lifestyle. Despite several advances in techniques used to identify and characterize Phytophthora effectors, knowledge is still lacking for some important species, including Phytophthora cinnamomi. This review discusses what the term PCD means and the gap in knowledge between pathogenic and developmental forms of PCD in plants. We also discuss the role cell death plays in the virulence of Phytophthora spp. and the effectors that have so far been identified as playing a role in cell death manipulation. Finally, we touch on the different techniques available to study effector functions, such as cell death induction/suppression.

Profile of the in silico secretome of the palm dieback pathogen, Fusarium oxysporum f. sp. albedinis, a fungus that puts natural oases at risk

Understanding biotic changes that occur alongside climate change constitute a research priority of global significance. Here, we address a plant pathogen that poses a serious threat to life on natural oases, where climate change is already taking a toll and severely impacting human subsistence. Fusarium oxysporum f. sp. albedinis is a pathogen that causes dieback disease on date palms, a tree that provides several critical ecosystem services in natural oases; and consequently, of major importance in this vulnerable habitat. Here, we assess the current state of global pathogen spread, we annotate the genome of a sequenced pathogen strain isolated from the native range and we analyse its in silico secretome. The palm dieback pathogen secretes a large arsenal of effector candidates including a variety of toxins, a distinguished profile of secreted in xylem proteins (SIX) as well as an expanded protein family with an N-terminal conserved motif [SG]PC[KR]P that could be involved in interactions with host membranes. Using agrobiodiversity as a strategy to decrease pathogen infectivity, while providing short term resilient solutions, seems to be widely overcome by the pathogen. Hence, the urgent need for future mechanistic research on the palm dieback disease and a better understanding of pathogen genetic diversity.

Whole genome sequencing and phylogenomic analysis show support for the splitting of genus Pythium

The genus Pythium (nom. cons.) sensu lato (s.l.) is composed of many important species of plant pathogens. Early molecular phylogenetic studies suggested paraphyly of Pythium, which led to a formal proposal by Uzuhashi and colleagues in 2010 to split the genus into Pythium sensu stricto (s.s.), Elongisporangium, Globisporangium, Ovatisporangium (= Phytopythium), and Pilasporangium using morphological characters and phylogenies of the mt cytochrome c oxidase subunit 2 (cox2) and D1-D2 domains of nuc 28S rDNA. Although the split was fairly justified by the delineating morphological characters, there were weaknesses in the molecular analyses, which created reluctance in the scientific community to adopt these new genera for the description of new species. In this study, this issue was addressed using phylogenomics. Whole genomes of 109 strains of Pythium and close relatives were sequenced, assembled, and annotated. These data were combined with 10 genomes sequenced in previous studies. Phylogenomic analyses were performed with 148 single-copy genes represented in at least 90% of the taxa in the data set. The results showed support for the division of Pythium s.l. The status of alternative generic names that have been used for species of Pythium in the past (e.g., Artotrogus, Cystosiphon, Eupythium, Nematosporangium, Rheosporangium, Sphaerosporangium) was investigated. Based on our molecular analyses and review of the Pythium generic concepts, we urge the scientific community to adopt the generic names Pythium, Elongisporangium, Globisporangium, and their concepts as proposed by Uzuhashi and colleagues in 2010 in their work going forward. In order to consolidate the taxonomy of these genera, some of the recently described Pythium spp. are transferred to Elongisporangium and Globisporangium.

Comparative analyses of complete Peronosporaceae (Oomycota) mitogenome sequences - insights into structural evolution and phylogeny

Members of the Peronosporaceae (Oomycota, Chromista), which currently consists of 25 genera and approximately 1000 recognised species, are responsible for disease on a wide range of plant hosts. Molecular phylogenetic analyses over the last two decades have improved our understanding of evolutionary relationships within Peronosporaceae. To date, 16 numbered and three named clades have been recognised; it is clear from these studies that the current taxonomy does not reflect evolutionary relationships. Whole organelle genome sequences are an increasingly important source of phylogenetic information, and in this study we present comparative and phylogenetic analyses of mitogenome sequences from 15 of the 19 currently recognized clades of Peronosporaceae, including 44 newly assembled sequences. Our analyses suggest strong conservation of mitogenome size and gene content across Peronosporaceae but, as previous studies have suggested, limited conservation of synteny. Specifically, we identified 28 distinct syntenies amongst the 71 examined isolates. Moreover, 19 of the isolates contained inverted or direct repeats, suggesting repeated sequences may be more common than previously thought. In terms of phylogenetic relationships, our analyses of 34 concatenated mitochondrial gene sequences resulted in a topology that was broadly consistent with previous studies. However, unlike previous studies concatenated mitochondrial sequences provided strong support for higher level relationships within the family.

A Comprehensive Assessment of the Secretome Responsible for Host Adaptation of the Legume Root Pathogen Aphanomyces euteiches

The soil-borne oomycete pathogen Aphanomyces euteiches causes devastating root rot diseases in legumes such as pea and alfalfa. The different pathotypes of A. euteiches have been shown to exhibit differential quantitative virulence, but the molecular basis of host adaptation has not yet been clarified. Here, we re-sequenced a pea field reference strain of A. euteiches ATCC201684 with PacBio long-reads and took advantage of the technology to generate the mitochondrial genome. We identified that the secretome of A. euteiches is characterized by a large portfolio of secreted proteases and carbohydrate-active enzymes (CAZymes). We performed Illumina sequencing of four strains of A. euteiches with contrasted specificity to pea or alfalfa and found in different geographical areas. Comparative analysis showed that the core secretome is largely represented by CAZymes and proteases. The specific secretome is mainly composed of a large set of small, secreted proteins (SSP) without any predicted functional domain, suggesting that the legume preference of the pathogen is probably associated with unknown functions. This study forms the basis for further investigations into the mechanisms of interaction of A. euteiches with legumes.

"Core" RxLR effectors in phytopathogenic oomycetes: A promising way to breeding for durable resistance in plants?

Phytopathogenic oomycetes are known to successfully infect their hosts due to their ability to secrete effector proteins. Of interest to many researchers are effectors with the N-terminal RxLR motif (Arginine-any amino acid-Leucine-Arginine). Owing to advances in genome sequencing, we can now comprehend the high level of diversity among oomycete effectors, and similarly, their conservation within and among species referred to here as "core" RxLR effectors (CREs). Currently, there is a considerable number of CREs that have been identified in oomycetes. Functional characterization of these CREs propose their virulence role with the potential of targeting central cellular processes that are conserved across diverse plant species. We reason that effectors that are highly conserved and recognized by the host, could be harnessed in engineering plants for durable as well as broad-spectrum resistance.

Specific interaction of an RNA-binding protein with the 3'-UTR of its target mRNA is critical to oomycete sexual reproduction

Sexual reproduction is an essential stage of the oomycete life cycle. However, the functions of critical regulators in this biological process remain unclear due to a lack of genome editing technologies and functional genomic studies in oomycetes. The notorious oomycete pathogen Pythium ultimum is responsible for a variety of diseases in a broad range of plant species. In this study, we revealed the mechanism through which PuM90, a stage-specific Puf family RNA-binding protein, regulates oospore formation in P. ultimum. We developed the first CRISPR/Cas9 system-mediated gene knockout and in situ complementation methods for Pythium. PuM90-knockout mutants were significantly defective in oospore formation, with empty oogonia or oospores larger in size with thinner oospore walls compared with the wild type. A tripartite recognition motif (TRM) in the Puf domain of PuM90 could specifically bind to a UGUACAUA motif in the mRNA 3′ untranslated region (UTR) of PuFLP, which encodes a flavodoxin-like protein, and thereby repress PuFLP mRNA level to facilitate oospore formation. Phenotypes similar to PuM90-knockout mutants were observed with overexpression of PuFLP, mutation of key amino acids in the TRM of PuM90, or mutation of the 3′-UTR binding site in PuFLP. The results demonstrated that a specific interaction of the RNA-binding protein PuM90 with the 3′-UTR of PuFLP mRNA at the post-transcriptional regulation level is critical for the sexual reproduction of P. ultimum.

Oomycetes are a class of eukaryotic microorganisms with life cycles and growth habits similar to filamentous fungi, but are not true fungi. Although sexual reproduction, which produce oospores, is an essential stage of life cycle, the functions of critical regulators in this biological process remain unclear due to a lack of genome editing technologies and functional genomic studies in oomycetes. In this study, we developed the first CRISPR/Cas9 system-mediated gene knockout and in situ complementation methods for Pythium ultimum, a notorious oomycete pathogen that is responsible for a variety of diseases in a broad range of plant species. We further identified the Puf family RNA-binding protein PuM90 and the flavodoxin-like protein PuFLP as major functional factors involved in P. ultimum oospore formation. We proposed a new model that PuM90 acts as a stage-specific post-transcriptional regulator by specifically binding to the 3′-UTR of PuFLP and then repressing PuFLP mRNA level. This study describes new technologies and data that will help to elucidate sexual reproduction and post-transcriptional regulation in oomycetes.

Novel Fusari- and Toti-like Viruses, with Probable Different Origins, in the Plant Pathogenic Oomycete Globisporangiumultimum

To further classify the oomycete viruses that have been discovered in recent years, we investigated virus infection in the plant-parasitic oomycete Globisporangium ultimum in Japan. Double-stranded RNA detection, high-throughput sequencing, and RT-PCR revealed that the G. ultimum isolate UOP226 contained two viruses related to fusarivirus and totivirus, named Pythium ultimum RNA virus 1 (PuRV1) and Pythium ultimum RNA virus 2 (PuRV2), respectively. Phylogenetic analysis of the deduced amino acid sequence of the RNA-dependent RNA polymerase (RdRp) showed that fusari-like PuRV1 belonged to a different phylogenetic group than Plasmopara viticola lesion-associated fusari virus (PvlaFV) 1-3 from oomycete Plasmopara viticola. Codon usage bias of the PuRV1 RdRp gene was more similar to those of fungi than Globisporangium and Phytophthora, suggesting that the PuRV1 ancestor horizontally transmitted to G. ultimum ancestor from fungi. Phylogenetic analysis of the deduced amino acid sequence of the RdRp of toti-like PuRV2 showed a monophyletic group with the other toti-like oomycete viruses from Globisporangium, Phytophthora, and Pl. viticola. However, the nucleotide sequences of toti-like oomycete viruses were not so homologous, suggesting the possibility of convergent evolution of toti-like oomycete viruses.

Insights into the Host Specificity of a New Oomycete Root Pathogen, Pythium brassicum P1: Whole Genome Sequencing and Comparative Analysis Reveals Contracted Regulation of Metabolism, Protein Families, and Distinct Pathogenicity Repertoire

Pythium brassicum P1 Stanghellini, Mohammadi, Förster, and Adaskaveg is an oomycete root pathogen that has recently been characterized. It only attacks plant species belonging to Brassicaceae family, causing root necrosis, stunting, and yield loss. Since P. brassicum P1 is limited in its host range, this prompted us to sequence its whole genome and compare it to those of broad host range Pythium spp. such as P. aphanidermatum and P. ultimum var. ultimum. A genomic DNA library was constructed with a total of 374 million reads. The sequencing data were assembled using SOAPdenovo2, yielding a total genome size of 50.3 Mb contained in 5434 scaffolds, N50 of 30.2 Kb, 61.2% G+C content, and 13,232 putative protein-coding genes. Pythium brassicum P1 had 175 species-specific gene families, which is slightly below the normal average. Like P. ultimum, P. brassicum P1 genome did not encode any classical RxLR effectors or cutinases, suggesting a significant difference in virulence mechanisms compared to other oomycetes. Pythium brassicum P1 had a much smaller proportions of the YxSL sequence motif in both secreted and non-secreted proteins, relative to other Pythium species. Similarly, P. brassicum P1 had the fewest Crinkler (CRN) effectors of all the Pythium species. There were 633 proteins predicted to be secreted in the P. brassicum P1 genome, which is, again, slightly below average among Pythium genomes. Pythium brassicum P1 had only one cadherin gene with calcium ion-binding LDRE and DxND motifs, compared to Pythium ultimum having four copies. Pythium brassicum P1 had a reduced number of proteins falling under carbohydrate binding module and hydrolytic enzymes. Pythium brassicum P1 had a reduced complement of cellulase and pectinase genes in contrast to P. ultimum and was deficient in xylan degrading enzymes. The contraction in ABC transporter families in P. brassicum P1 is suggested to be the result of a lack of diversity in nutrient uptake and therefore host range.

Genome Analysis of Phytophthora nicotianae JM01 Provides Insights into Its Pathogenicity Mechanisms

Phytophthora nicotianae is a widely distributed plant pathogen that can cause serious disease and cause significant economic losses to various crops, including tomatoes, tobacco, onions, and strawberries. To understand its pathogenic mechanisms and explore strategies for controlling diseases caused by this pathogen, we sequenced and analyzed the whole genome of Ph. nicotianae JM01. The Ph. nicotianae JM01 genome was assembled using a combination of approaches including shotgun sequencing, single-molecule sequencing, and the Hi-C technique. The assembled Ph. nicotianae JM01 genome is about 95.32 Mb, with contig and scaffold N50 54.23 kb and 113.15 kb, respectively. The average GC content of the whole-genome is about 49.02%, encoding 23,275 genes. In addition, we identified 19.15% of interspersed elements and 0.95% of tandem elements in the whole genome. A genome-wide phylogenetic tree indicated that Phytophthora diverged from Pythium approximately 156.32 Ma. Meanwhile, we found that 252 and 285 gene families showed expansion and contraction in Phytophthora when compared to gene families in Pythium. To determine the pathogenic mechanisms Ph. nicotianae JM01, we analyzed a suite of proteins involved in plant-pathogen interactions. The results revealed that gene duplication contributed to the expansion of Cell Wall Degrading Enzymes (CWDEs) such as glycoside hydrolases, and effectors such as Arg-Xaa-Leu-Arg (RXLR) effectors. In addition, transient expression was performed on Nicotiana benthamiana by infiltrating with Agrobacterium tumefaciens cells containing a cysteine-rich (SCR) protein. The results indicated that SCR can cause symptoms of hypersensitive response. Moreover, we also conducted comparative genome analysis among four Ph. nicotianae genomes. The completion of the Ph. nicotianae JM01 genome can not only help us understand its genomic characteristics, but also help us discover genes involved in infection and then help us understand its pathogenic mechanisms.

A Small Cysteine-Rich Phytotoxic Protein of Phytophthora capsici Functions as Both Plant Defense Elicitor and Virulence Factor

Small cysteine-rich (SCR) proteins, including fungal avirulence proteins, play important roles in pathogen-plant interactions. SCR protein-encoding genes have been discovered in the genomes of Phytophthora pathogens but their functions during pathogenesis remain obscure. Here, we report the characterization of one Phytophthora capsici SCR protein (namely, SCR82) with similarity to Phytophthora cactorum phytotoxic protein PcF. The scr82 gene has 10 allelic sequences in the P. capsici population. Homologs of SCR82 were not identified in fungi or other organisms but in Phytophthora relative species. Initially, scr82 was weakly expressed during the mycelium, sporangium, and zoospore stages but quickly upregulated when the infection initiated. Both ectopic expression of SCR82 and recombinant yeast-expressed protein (rSCR82) caused cell death on tomato leaves. Upon treatment, rSCR82 induced plant defense responses, including the induction of defense gene expression, reactive oxygen species burst, and callose deposition. Knockout of scr82 in P. capsici by CRISPR/Cas9 severely impaired its virulence on host plants and significantly reduced its resistance against oxidative stress. Inversely, its overexpression increased the pathogen's virulence and tolerance to oxidative stress. Our results collectively demonstrate that SCR82 functions as both an important virulence factor and plant defense elicitor, which is conserved across Phytophthora spp.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Predicting Lifestyle from Positive Selection Data and Genome Properties in Oomycetes

As evidenced in parasitism, host and niche shifts are a source of genomic and phenotypic diversification. Exemplary is a reduction in the core metabolism as parasites adapt to a particular host, while the accessory genome often maintains a high degree of diversification. However, selective pressures acting on the genome of organisms that have undergone recent lifestyle or host changes have not been fully investigated. Here, we developed a comparative genomics approach to study underlying adaptive trends in oomycetes, a eukaryotic phylum with a wide and diverse range of economically important plant and animal parasitic lifestyles. Our analysis reveals converging evolution on biological processes for oomycetes that have similar lifestyles. Moreover, we find that certain functions, in particular carbohydrate metabolism, transport, and signaling, are important for host and environmental adaptation in oomycetes. Given the high correlation between lifestyle and genome properties in our oomycete dataset, together with the known convergent evolution of fungal and oomycete genomes, we developed a model that predicts plant pathogenic lifestyles with high accuracy based on functional annotations. These insights into how selective pressures correlate with lifestyle may be crucial to better understand host/lifestyle shifts and their impact on the genome.

Diversity and Function of Appressoria

Endophytic, saprobic, and pathogenic fungi have evolved elaborate strategies to obtain nutrients from plants. Among the diverse plant-fungi interactions, the most crucial event is the attachment and penetration of the plant surface. Appressoria, specialized infection structures, have been evolved to facilitate this purpose. In this review, we describe the diversity of these appressoria and classify them into two main groups: single-celled appressoria (proto-appressoria, hyaline appressoria, melanized (dark) appressoria) and compound appressoria. The ultrastructure of appressoria, their initiation, their formation, and their function in fungi are discussed. We reviewed the molecular mechanisms regulating the formation and function of appressoria, their strategies to evade host defenses, and the related genomics and transcriptomics. The current review provides a foundation for comprehensive studies regarding their evolution and diversity in different fungal groups.

Comparative analyses of saprotrophy in Salisapilia sapeloensis and diverse plant pathogenic oomycetes reveal lifestyle-specific gene expression

Oomycetes include many devastating plant pathogens. Across oomycete diversity, plant-infecting lineages are interspersed by non-pathogenic ones. Unfortunately, our understanding of the evolution of lifestyle switches is hampered by a scarcity of data on the molecular biology of saprotrophic oomycetes, ecologically important primary colonizers of dead tissue that can serve as informative reference points for understanding the evolution of pathogens. Here, we established Salisapilia sapeloensis as a tractable system for the study of saprotrophic oomycetes. We generated multiple transcriptomes from S. sapeloensis and compared them with (i) 22 oomycete genomes and (ii) the transcriptomes of eight pathogenic oomycetes grown under 13 conditions. We obtained a global perspective on gene expression signatures of oomycete lifestyles. Our data reveal that oomycete saprotrophs and pathogens use similar molecular mechanisms for colonization but exhibit distinct expression patterns. We identify a S. sapeloensis-specific array and expression of carbohydrate-active enzymes and putative regulatory differences, highlighted by distinct expression levels of transcription factors. Salisapilia sapeloensis expresses only a small repertoire of candidates for virulence-associated genes. Our analyses suggest lifestyle-specific gene regulatory signatures and that, in addition to variation in gene content, shifts in gene regulatory networks underpin the evolution of oomycete lifestyles.

Pythium Damping-Off and Root Rot of Capsicum annuum L.: Impacts, Diagnosis, and Management

Capsicum annuum L. is a significant horticulture crop known for its pungent varieties and used as a spice. The pungent character in the plant, known as capsaicinoid, has been discovered to have various health benefits. However, its production has been affected due to various exogenous stresses, including diseases caused by a soil-borne pathogen, Pythium spp. predominantly affecting the Capsicum plant in younger stages and causing damping-off, this pathogen can incite root rot in later plant growth stages. Due to the involvement of multiple Pythium spp. and their capability to disperse through various routes, their detection and diagnosis have become crucial. However, the quest for a point-of-care technology is still far from over. The use of an integrated approach with cultural and biological techniques for the management of Pythium spp. can be the best and most sustainable alternative to the traditionally used and hazardous chemical approach. The lack of race-specific resistance genes against Pythium spp. can be compensated with the candidate quantitative trait loci (QTL) genes in C. annuum L. This review will focus on the epidemiological factors playing a major role in disease spread, the currently available diagnostics in species identification, and the management strategies with a special emphasis on Pythium spp. causing damping-off and root rot in different cultivars of C. annuum L.

A Conserved Oomycete CRN Effector Targets Tomato TCP14-2 to Enhance Virulence

Phytophthora spp. secrete vast arrays of effector molecules during infection to aid in host colonization. The crinkling and necrosis (CRN) protein family forms an extensive repertoire of candidate effectors that accumulate in the host nucleus to perturb processes required for immunity. Here, we show that CRN12_997 from Phytophthora capsici binds a TCP transcription factor, SlTCP14-2, to inhibit its immunity-associated activity against Phytophthora spp. Coimmunoprecipitation and bimolecular fluorescence complementation studies confirm a specific CRN12_997-SlTCP14-2 interaction in vivo. Coexpression of CRN12_997 specifically counteracts the TCP14-enhanced immunity phenotype, suggesting that CRN mediated perturbation of SlTCP14-2 function. We show that SlTCP14-2 associates with nuclear chromatin and that CRN12_997 diminishes SlTCP14-2 DNA binding. Collectively, our data support a model in which SlTCP14-2 associates with chromatin to enhance immunity. The interaction between CRN12_997 and SlTCP14-2 reduces DNA binding of the immune regulator. We propose that the modulation of SlTCP14-2 chromatin affinity, caused by CRN12-997, enhances susceptibility to P. capsici.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Horizontal Gene Transfer and Tandem Duplication Shape the Unique CAZyme Complement of the Mycoparasitic Oomycetes Pythium oligandrum and Pythium periplocum

Crop protection strategies that are effective but that reduce our reliance on chemical pesticides are urgently needed to meet the UN sustainable development goals for global food security. Mycoparasitic oomycetes such as Pythium oligandrum and Pythium periplocum, have potential for the biological control of plant diseases that threaten crops and have attracted much attention due to their abilities to antagonize plant pathogens and modulate plant immunity. Studies of the molecular and genetic determinants of mycoparasitism in these species have been less well developed than those of their fungal counterparts. Carbohydrate-active enzymes (CAZymes) from P. oligandrum and P. periplocum are predicted to be important components of mycoparasitism, being involved in the degradation of the cell wall of their oomycete and fungal prey species. To explore the evolution of CAZymes of these species we performed an in silico identification and comparison of the full CAZyme complement (CAZyome) of the two mycoparasitic Pythium species (P. oligandrum and P. periplocum), with seven other Pythium species, and four Phytophthora species. Twenty CAZy gene families involved in the degradation of cellulose, hemicellulose, glucan, and chitin were expanded in, or unique to, mycoparasitic Pythium species and several of these genes were expressed during mycoparasitic interactions with either oomycete or fungal prey, as revealed by RNA sequencing and quantitative qRT-PCR. Genes from three of the cellulose and chitin degrading CAZy families (namely AA9, GH5_14, and GH19) were expanded via tandem duplication and predominantly located in gene sparse regions of the genome, suggesting these enzymes are putative pathogenicity factors able to undergo rapid evolution. In addition, five of the CAZy gene families were likely to have been obtained from other microbes by horizontal gene transfer events. The mycoparasitic species are able to utilize complex carbohydrates present in fungal cell walls, namely chitin and N-acetylglucosamine for growth, in contrast to their phytopathogenic counterparts. Nonetheless, a preference for the utilization of simple sugars for growth appears to be a common trait within the oomycete lineage.

What are the Top 10 Unanswered Questions in Molecular Plant-Microbe Interactions?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.The past few decades have seen major discoveries in the field of molecular plant-microbe interactions. As the result of technological and intellectual advances, we are now able to answer questions at a level of mechanistic detail that we could not have imagined possible 20 years ago. The MPMI Editorial Board felt it was time to take stock and reassess. What big questions remain unanswered? We knew that to identify the fundamental, overarching questions that drive our research, we needed to do this as a community. To reach a diverse audience of people with different backgrounds and perspectives, working in different areas of plant-microbe interactions, we queried the more than 1,400 participants at the 2019 International Congress on Molecular Plant-Microbe Interactions meeting in Glasgow. This group effort resulted in a list of ten, broad-reaching, fundamental questions that influence and inform our research. Here, we introduce these Top 10 unanswered questions, giving context and a brief description of the issues. Each of these questions will be the subject of a detailed review in the coming months. We hope that this process of reflecting on what is known and unknown and identifying the themes that underlie our research will provide a framework to use going forward, giving newcomers a sense of the mystery of the big questions and inspiring new avenues and novel insights.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Effector prediction and characterization in the oomycete pathogen Bremia lactucae reveal host-recognized WY domain proteins that lack the canonical RXLR motif

Pathogens that infect plants and animals use a diverse arsenal of effector proteins to suppress the host immune system and promote infection. Identification of effectors in pathogen genomes is foundational to understanding mechanisms of pathogenesis, for monitoring field pathogen populations, and for breeding disease resistance. We identified candidate effectors from the lettuce downy mildew pathogen Bremia lactucae by searching the predicted proteome for the WY domain, a structural fold found in effectors that has been implicated in immune suppression as well as effector recognition by host resistance proteins. We predicted 55 WY domain containing proteins in the genome of B. lactucae and found substantial variation in both sequence and domain architecture. These candidate effectors exhibit several characteristics of pathogen effectors, including an N-terminal signal peptide, lineage specificity, and expression during infection. Unexpectedly, only a minority of B. lactucae WY effectors contain the canonical N-terminal RXLR motif, which is a conserved feature in the majority of cytoplasmic effectors reported in Phytophthora spp. Functional analysis of 21 effectors containing WY domains revealed 11 that elicited cell death on wild accessions and domesticated lettuce lines containing resistance genes, indicative of recognition of these effectors by the host immune system. Only two of the 11 recognized effectors contained the canonical RXLR motif, suggesting that there has been an evolutionary divergence in sequence motifs between genera; this has major consequences for robust effector prediction in oomycete pathogens.

The Effector Repertoire of the Hop Downy Mildew Pathogen Pseudoperonospora humuli

Pseudoperonospora humuli is an obligate biotrophic oomycete that causes downy mildew (DM), one of the most destructive diseases of cultivated hop that can lead to 100% crop loss in susceptible cultivars. We used the published genome of P. humuli to predict the secretome and effectorome and analyze the transcriptome variation among diverse isolates and during infection of hop leaves. Mining the predicted coding genes of the sequenced isolate OR502AA of P. humuli revealed a secretome of 1,250 genes. We identified 296 RXLR and RXLR-like effector-encoding genes in the secretome. Among the predicted RXLRs, there were several WY-motif-containing effectors that lacked canonical RXLR domains. Transcriptome analysis of sporangia from 12 different isolates collected from various hop cultivars revealed 754 secreted proteins and 201 RXLR effectors that showed transcript evidence across all isolates with reads per kilobase million (RPKM) values > 0. RNA-seq analysis of OR502AA-infected hop leaf samples at different time points after infection revealed highly expressed effectors that may play a relevant role in pathogenicity. Quantitative RT-PCR analysis confirmed the differential expression of selected effectors. We identified a set of P. humuli core effectors that showed transcript evidence in all tested isolates and elevated expression during infection. These effectors are ideal candidates for functional analysis and effector-assisted breeding to develop DM resistant hop cultivars.

A Global Survey of Carbohydrate Esterase Families 1 and 10 in Oomycetes

Carbohydrate-active enzymes (CAZymes) are a cornerstone in the phytopathogenicity of filamentous microbes. CAZymes are required for every step of a successful infection cycle-from penetration, to nutrient acquisition (during colonization), to exit and dispersal. Yet, CAZymes are not a unique feature of filamentous pathogens. They are found across eukaryotic genomes and including, for example, saprotrophic relatives of major pathogens. Comparative genomics and functional analyses revealed that CAZyme content is shaped by a multitude of factors, including utilized substrate, lifestyle, and host preference. Yet, family size alone says little about usage. Indeed, in a previous study, we found that genes putatively coding for the CAZyme families of carbohydrate esterase (CE)1 and CE10, while not specifically enriched in number, were suggested to have lifestyle-specific gene expression patterns. Here, we used comparative genomics and a clustering approach to understand how the repertoire of the CE1- and CE10-encoding gene families is shaped across oomycete evolution. These data are combined with comparative transcriptomic analyses across homologous clusters within the gene families. We find that CE1 and CE10 have been reduced in number in biotrophic oomycetes independent of the phylogenetic relationship of the biotrophs to each other. The reduction in CE1 is different from that observed for CE10: While in CE10 specific clusters of homologous sequences show convergent reduction, CE1 reduction is caused by species-specific losses. Comparative transcriptomics revealed that some clusters of CE1 or CE10 sequences have a higher expression than others, independent of the species composition within them. Further, we find that CE1- and CE10-encoding genes are mainly induced in plant pathogens and that some homologous genes show lifestyle-specific gene expression levels during infection, with hemibiotrophs showing the highest expression levels.

Prediction and Characterization of RXLR Effectors in Pythium Species

RXLR effectors, a class of secreted proteins that are transferred into host cells to manipulate host immunity, have been reported to widely exist in oomycetes, including those from genera Phytophthora, Hyaloperonospora, Albugo, and Saprolegnia. However, in Pythium species, no RXLR effector has yet been characterized, and the origin and evolution of such virulent effectors are still unknown. Here, we developed a modified regular expression method for de novo identification of RXLRs and characterized 359 putative RXLR effectors in nine Pythium species. Phylogenetic analysis revealed that all oomycetous RXLRs formed a single superfamily, suggesting that they might have a common ancestor. RXLR effectors from Pythium and Phytophthora species exhibited similar sequence features, protein structures, and genome locations. In particular, there were significantly more RXLR proteins in the mosquito biological control agent P. guiyangense than in the other eight Pythium species, and P. guiyangense RXLRs might be the result of gene duplication and genome rearrangement events, as indicated by synteny analysis. Expression pattern analysis of RXLR-encoding genes in the plant pathogen P. ultimum detected transcripts of the majority of the predicted RXLR genes, with some RXLR effectors induced in infection stages and one RXLR showing necrosis-inducing activity. Furthermore, all predicted RXLR genes were cloned from two biocontrol agents, P. oligandrum and P. periplocum, and three of the RXLR genes were found to induce a defense response in Nicotiana benthamiana. Taken together, our findings represent the first evidence of RXLR effectors in Pythium species, providing valuable information on their evolutionary patterns and the mechanisms of their interactions with diverse hosts.

Fungi, bacteria and oomycota opportunistically isolated from the seagrass, Zostera marina

Fungi in the marine environment are often neglected as a research topic, despite that fungi having critical roles on land as decomposers, pathogens or endophytes. Here we used culture-dependent methods to survey the fungi associated with the seagrass, Zostera marina, also obtaining bacteria and oomycete isolates in the process. A total of 108 fungi, 40 bacteria and 2 oomycetes were isolated. These isolates were then taxonomically identified using a combination of molecular and phylogenetic methods. The majority of the fungal isolates were classified as belonging to the classes Eurotiomycetes, Dothideomycetes, and Sordariomycetes. Most fungal isolates were habitat generalists like Penicillium sp. and Cladosporium sp., but we also cultured a diverse set of rare taxa including possible habitat specialists like Colletotrichum sp. which may preferentially associate with Z. marina leaf tissue. Although the bulk of bacterial isolates were identified as being from known ubiquitous marine lineages, we also obtained several Actinomycetes isolates and a Phyllobacterium sp. We identified two oomycetes, another understudied group of marine microbial eukaryotes, as Halophytophthora sp. which may be opportunistic pathogens or saprophytes of Z. marina. Overall, this study generates a culture collection of fungi which adds to knowledge of Z. marina associated fungi and highlights a need for more investigation into the functional and evolutionary roles of microbial eukaryotes associated with seagrasses.

A terminal α3-galactose modification regulates an E3 ubiquitin ligase subunit in Toxoplasma gondii

Skp1, a subunit of E3 Skp1/Cullin-1/F-box protein ubiquitin ligases, is modified by a prolyl hydroxylase that mediates O2 regulation of the social amoeba Dictyostelium and the parasite Toxoplasma gondii The full effect of hydroxylation requires modification of the hydroxyproline by a pentasaccharide that, in Dictyostelium, influences Skp1 structure to favor assembly of Skp1/F-box protein subcomplexes. In Toxoplasma, the presence of a contrasting penultimate sugar assembled by a different glycosyltransferase enables testing of the conformational control model. To define the final sugar and its linkage, here we identified the glycosyltransferase that completes the glycan and found that it is closely related to glycogenin, an enzyme that may prime glycogen synthesis in yeast and animals. However, the Toxoplasma enzyme catalyzes formation of a Galα1,3Glcα linkage rather than the Glcα1,4Glcα linkage formed by glycogenin. Kinetic and crystallographic experiments showed that the glycosyltransferase Gat1 is specific for Skp1 in Toxoplasma and also in another protist, the crop pathogen Pythium ultimum The fifth sugar is important for glycan function as indicated by the slow-growth phenotype of gat1Δ parasites. Computational analyses indicated that, despite the sequence difference, the Toxoplasma glycan still assumes an ordered conformation that controls Skp1 structure and revealed the importance of nonpolar packing interactions of the fifth sugar. The substitution of glycosyltransferases in Toxoplasma and Pythium by an unrelated bifunctional enzyme that assembles a distinct but structurally compatible glycan in Dictyostelium is a remarkable case of convergent evolution, which emphasizes the importance of the terminal α-galactose and establishes the phylogenetic breadth of Skp1 glycoregulation.

Mitochondrial genome sequence of Phytophthora sansomeana and comparative analysis of Phytophthora mitochondrial genomes

Phytophthora sansomeana infects soybean and causes root rot. It was recently separated from the species complex P. megasperma sensu lato. In this study, we sequenced and annotated its complete mitochondrial genome and compared it to that of nine other Phytophthora species. The genome was assembled into a circular molecule of 39,618 bp with a 22.03% G+C content. Forty-two protein coding genes, 25 tRNA genes and two rRNA genes were annotated in this genome. The protein coding genes include 14 genes in the respiratory complexes, four ATP synthase genes, 16 ribosomal proteins genes, a tatC translocase gene, six conserved ORFs and a unique orf402. The tRNA genes encode tRNAs for 19 amino acids. Comparison among mitochondrial genomes of 10 Phytophthora species revealed three inversions, each covering multiple genes. These genomes were conserved in gene content with few exceptions. A 3' truncated atp9 gene was found in P. nicotianae. All 10 Phytophthora species, as well as other oomycetes and stramenopiles, lacked tRNA genes for threonine in their mitochondria. Phylogenomic analysis using the mitochondrial genomes supported or enhanced previous findings of the phylogeny of Phytophthora spp.

The Genome of Peronospora belbahrii Reveals High Heterozygosity, a Low Number of Canonical Effectors, and TC-Rich Promoters

Along with Plasmopara destructor, Peronosopora belbahrii has arguably been the economically most important newly emerging downy mildew pathogen of the past two decades. Originating from Africa, it has started devastating basil production throughout the world, most likely due to the distribution of infested seed material. Here, we present the genome of this pathogen and results from comparisons of its genomic features to other oomycetes. The assembly of the nuclear genome was around 35.4 Mbp in length, with an N₅₀ scaffold length of around 248 kbp and an L₅₀ scaffold count of 46. The circular mitochondrial genome consisted of around 40.1 kbp. From the repeat-masked genome, 9,049 protein-coding genes were predicted, out of which 335 were predicted to have extracellular functions, representing the smallest secretome so far found in peronosporalean oomycetes. About 16% of the genome consists of repetitive sequences, and, based on simple sequence repeat regions, we provide a set of microsatellites that could be used for population genetic studies of P. belbahrii. P. belbahrii has undergone a high degree of convergent evolution with other obligate parasitic pathogen groups, reflecting its obligate biotrophic lifestyle. Features of its secretome, signaling networks, and promoters are presented, and some patterns are hypothesized to reflect the high degree of host specificity in Peronospora species. In addition, we suggest the presence of additional virulence factors apart from classical effector classes that are promising candidates for future functional studies.

Recent advances in oomycete genomics

The oomycetes are a class of ubiquitous, filamentous microorganisms that include some of the biggest threats to global food security and natural ecosystems. Within the oomycete class are highly diverse species that infect a broad range of animals and plants. Some of the most destructive plant pathogens are oomycetes, such as Phytophthora infestans, the agent of potato late blight and the cause of the Irish famine. Recent years have seen a dramatic increase in the number of sequenced oomycete genomes. Here we review the latest developments in oomycete genomics and some of the important insights that have been gained. Coupled with proteomic and transcriptomic analyses, oomycete genome sequences have revealed tremendous insights into oomycete biology, evolution, genome organization, mechanisms of infection, and metabolism. We also present an updated phylogeny of the oomycete class using a phylogenomic approach based on the 65 oomycete genomes that are currently available.

Skp1 isoforms are differentially modified by a dual function prolyl 4-hydroxylase/N-acety lglucosaminyltransferase in a plant pathogen

Skp1 is hydroxylated by an O2-dependent prolyl hydroxylase (PhyA) that contributes to O2-sensing in the social amoeba Dictyostelium and the mammalian pathogen Toxoplasma gondii. HO-Skp1 is subject to glycosylation and the resulting pentasaccharide affects Skp1 conformation in a way that influences association of Skp1 with F-box proteins, and potentially the assembly of E3(SCF) ubiquitin ligase complexes that mediate the polyubiquitination of target proteins that are degraded in the 26S-proteasome. To investigate the conservation and specificity of these modifications, we analyzed proteins from the oomycete Pythium ultimum, an important crop plant pathogen. Putative coding sequences for Pythium's predicted PhyA and first glycosyltransferase in the predicted five-enzyme pathway, a GlcNAc-transferase (Gnt1), predict a bifunctional enzyme (Phgt) that, when expressed in Dictyostelium, rescued a knockout of phyA but not gnt1. Though recombinant Phgt was also unable to glycosylate Dictyostelium HO-Skp1, it could hydrolyze UDP-GlcNAc and modify a synthetic hydroxypeptide from Dictyostelium Skp1. Pythium encodes two highly similar Skp1 isoforms, but only Skp1A was efficiently hydroxylated and glycosylated in vitro. While kinetic analysis revealed no evidence for processive processing of Skp1, the physical linkage of the two activities implies dedication to Skp1 in vivo. These findings indicate a widespread occurrence of the Skp1 modification pathway across protist phylogeny, suggest that both Gnt1 and PhyA are specific for Skp1 and indicate that the second Skp1 provides a bypass mechanism for O2-regulation in Pythium and other protists that conserve this gene.

Promising Perspectives for Detection, Identification, and Quantification of Plant Pathogenic Fungi and Oomycetes through Targeting Mitochondrial DNA

Fungi and oomycetes encompass many pathogens affecting crops worldwide. Their effective control requires screening pathogens across the local and international trade networks along with the monitoring of pathogen inocula in the field. Fundamentals to all of these concerns are their efficient detection, identification, and quantification. The use of molecular markers showed the best promise in the field of plant pathogen diagnostics. However, despite the unquestionable benefits of DNA-based methods, two significant limitations are associated with their use. The first limitation concerns the insufficient level of sensitivity due to the very low and uneven distribution of pathogens in plant material. The second limitation pertains to the inability of widely used diagnostic assays to detect cryptic species. Targeting mtDNA appears to provide a solution to these challenges. Its high copy number in microbial cells makes mtDNA an attractive target for developing highly sensitive assays. In addition, previous studies on different pathogen taxa indicated that mitogenome sequence variation could improve cryptic species delimitation accuracy. This review sheds light on the potential application of mtDNA for pathogen diagnostics. This paper covers a brief description of qPCR and DNA barcoding as two major strategies enabling the diagnostics of plant pathogenic fungi and oomycetes. Both strategies are discussed along with the potential use of mtDNA, including their strengths and weaknesses.

Whole-genome sequencing of Puccinia striiformis f. sp. tritici mutant isolates identifies avirulence gene candidates

Background

The stripe rust pathogen, Puccinia striiformis f. sp. tritici (Pst), threats world wheat production. Resistance to Pst is often overcome by pathogen virulence changes, but the mechanisms of variation are not clearly understood. To determine the role of mutation in Pst virulence changes, in previous studies 30 mutant isolates were developed from a least virulent isolate using ethyl methanesulfonate (EMS) mutagenesis and phenotyped for virulence changes. The progenitor isolate was sequenced, assembled and annotated for establishing a high-quality reference genome. In the present study, the 30 mutant isolates were sequenced and compared to the wide-type isolate to determine the genomic variation and identify candidates for avirulence (Avr) genes.

Results

The sequence reads of the 30 mutant isolates were mapped to the wild-type reference genome to identify genomic changes. After selecting EMS preferred mutations, 264,630 and 118,913 single nucleotide polymorphism (SNP) sites and 89,078 and 72,513 Indels (Insertion/deletion) were detected among the 30 mutant isolates compared to the primary scaffolds and haplotigs of the wild-type isolate, respectively. Deleterious variants including SNPs and Indels occurred in 1866 genes. Genome wide association analysis identified 754 genes associated with avirulence phenotypes. A total of 62 genes were found significantly associated to 16 avirulence genes after selection through six criteria for putative effectors and degree of association, including 48 genes encoding secreted proteins (SPs) and 14 non-SP genes but with high levels of association (P ≤ 0.001) to avirulence phenotypes. Eight of the SP genes were identified as avirulence-associated effectors with high-confidence as they met five or six criteria used to determine effectors.

Conclusions

Genome sequence comparison of the mutant isolates with the progenitor isolate unraveled a large number of mutation sites along the genome and identified high-confidence effector genes as candidates for avirulence genes in Pst. Since the avirulence gene candidates were identified from associated SNPs and Indels caused by artificial mutagenesis, these avirulence gene candidates are valuable resources for elucidating the mechanisms of the pathogen pathogenicity, and will be studied to determine their functions in the interactions between the wheat host and the Pst pathogen.