A rapid molecular diagnostic tool to discriminate alleles of avirulence genes and haplotypes of Phytophthora sojae using high-resolution melting analysis.

Effectors encoded by avirulence genes (Avr) interact with the Phytophthora sojae resistance gene (Rps) products to generate incompatible interactions. The virulence profile of P. sojae is rapidly evolving as a result of the large-scale deployment of Rps genes in soybean. For a successful exploitation of Rps genes, it is recommended that soybean growers use cultivars containing the Rps genes corresponding to Avr genes present in P. sojae populations present in their fields. Determination of the virulence profile of P. sojae isolates is critical for the selection of soybean cultivars. High-resolution melting curve (HRM) analysis is a powerful tool, first applied in medicine, for detecting mutations with potential applications in different biological fields. Here, we report the development of an HRM protocol, as an original approach to discriminate effectors, to differentiate P. sojae haplotypes for six Avr genes. An HRM assay was performed on 24 P. sojae isolates with different haplotypes collected from soybean fields across Canada. The results clearly confirmed that the HRM assay discriminated different virulence genotypes. Moreover, the HRM assay was able to differentiate multiple haplotypes representing small allelic variations. HRM-based prediction was validated by phenotyping assays. This HRM assay provides a unique, cost-effective and efficient tool to predict virulence pathotypes associated with six different Avr (1b, 1c, 1d, 1k, 3a and 6) genes from P. sojae, which can be applied in the deployment of appropriate Rps genes in soybean fields.

A unique effector secreted by Pseudozyma flocculosa mediates its biocontrol activity.

BACKGROUND: Pseudozyma flocculosa is a highly efficient biocontrol agent (BCA) of powdery mildews whose mode of action remains elusive. It is known to secrete unique effectors during its interaction with powdery mildews but effectors have never been shown to be part of the arsenal of a BCA. Here, we characterize the role of the effector Pf2826 released by Pseudozyma flocculosa during its tripartite interaction with barley and the pathogen fungus Blumeria graminis f. sp. hordei. RESULTS: We utilized CRISPR-Cas9-based genome editing and confirmed that secreted P. flocculosa effector Pf2826 is required for full biocontrol activity. We monitored the localization of the effector Pf2826 with C-terminal mCherry tag and found it localized around the haustoria and on powdery mildew spores. His-tagged Pf2826 recombinant protein was expressed, purified, and used as bait in a pull-down assay from total proteins extracted during the tripartite interaction. Potential interactors were identified by LC-MS/MS analysis after removing unspecific interactions found in the negative controls. A two-way yeast two-hybrid assay validated that Pf2826 interacted with barley pathogenesis-related (PR) proteins HvPR1a and chitinase and with an effector protein from powdery mildew. CONCLUSIONS: In contrast to the usual modes of action of competition, parasitism, and antibiosis ascribed to BCAs, this study shows that effector pf2826 plays a vital role in the biocontrol activity of P. flocculosa by interacting with plant PR proteins and a powdery mildew effector, altering the host-pathogen interaction.

RXLR effector gene Avr3a from Phytophthora sojae is recognized by Rps8 in soybean.

The use of resistance genes in elite soybean cultivars is one of the most widely used methods to manage Phytophthora sojae. This method relies on effector-triggered immunity, where a Resistant to P. sojae (Rps) gene product from the plant recognizes a specific effector from the pathogen, encoded by an avirulence (Avr) gene. Many Avr genes from P. sojae have been identified in the last decade, allowing a better exploitation of this type of resistance. The objective of the present study was to identify the Avr gene triggering immunity derived from the soybean resistance gene Rps8. The analysis of a segregating F2 progeny coupled with a genotyping-by-sequencing approach led to the identification of a putative Avr8 locus. The investigation of this locus using whole-genome sequencing data from 31 isolates of P. sojae identified Avr3a as the likely candidate for Avr8. Long-read sequencing also revealed that P. sojae isolates can carry up to five copies of the Avr3a gene, compared to the four previously reported. Haplotype and transcriptional analyses showed that amino acid changes and absence of Avr3a transcripts from P. sojae isolates caused changes in virulence towards Rps8. Functional analyses using CRISPR/Cas9 knockout and constitutive expression demonstrated that Rps8 interacted with Avr3a. We also showed that a specific allele of Avr3a is recognized by Rps3a but not Rps8. While Rps3a and Rps8 have been previously described as closely linked, this is the first report of a clear distinction hitherto undefined between these two resistance genes.

Mapping of partial resistance to Phytophthora sojae in soybean PIs using whole-genome sequencing reveals a major QTL.

In the last decade, more than 70 quantitative trait loci (QTL) related to soybean [Glycine max (L.) Merr.] partial resistance (PR) against Phytophthora sojae have been identified by genome-wide association studies (GWAS). However, most of them have either a minor effect on the resistance level or are specific to a single phenotypic variable or one isolate, thereby limiting their use in breeding programs. In this study, we have used an analytical approach combining (a) the phenotypic characterization of a diverse panel of 357 soybean accessions for resistance to P. sojae captured through a single variable, corrected dry weight; (b) a new hydroponic assay allowing the inoculation of a combination of P. sojae isolates covering the spectrum of commercially relevant Rps genes; and (c) exhaustive genotyping through whole-genome resequencing (WGS). This led to the identification of a novel P. sojae resistance QTL with a relatively major effect compared with the previously reported QTL. The QTL interval, spanning ∼500 kb on chromosome (Chr) 15, does not colocalize with previously reported QTL for P. sojae resistance. Plants carrying the favorable allele at this QTL were 60% more resistant. Eight genes were found to reside in the linkage disequilibrium (LD) block containing the peak single-nucleotide polymorphism (SNP) including Glyma.15G217100, which encodes a major latex protein (MLP)-like protein, with a functional annotation related to pathogen resistance. Expression analysis of Glyma.15G217100 indicated that it was nearly eight times more highly expressed in a group of plant introductions (PIs) carrying the resistant (R) allele compared with those carrying the susceptible (S) allele within a short period after inoculation. These results offer new and valuable options to develop improved soybean cultivars with broad resistance to P. sojae through marker-assisted selection.

Progress and Challenges in Elucidating the Functional Role of Effectors in the Soybean-Phytophthora sojae Interaction.

Phytophthora sojae, the agent responsible for stem and root rot, is one of the most damaging plant pathogens of soybean. To establish a compatible-interaction, P. sojae secretes a wide array of effector proteins into the host cell. These effectors have been shown to act either in the apoplastic area or the cytoplasm of the cell to manipulate the host cellular processes in favor of the development of the pathogen. Deciphering effector-plant interactions is important for understanding the role of P. sojae effectors in disease progression and developing approaches to prevent infection. Here, we review the subcellular localization, the host proteins, and the processes associated with P. sojae effectors. We also discuss the emerging topic of effectors in the context of effector-resistance genes interaction, as well as model systems and recent developments in resources and techniques that may provide a better understanding of the soybean-P. sojae interaction.

A reassessment of flocculosin-mediated biocontrol activity of Pseudozyma flocculosa through CRISPR/Cas9 gene editing.

Pseudozyma flocculosa is an epiphytic yeast with powerful antagonistic activity against powdery mildews. This activity has been associated with the production of a rare antifungal glycolipid, flocculosin. In spite of the discovery of a specific gene cluster for flocculosin synthesis, attempts to ascribe a functional role to the molecule have been hampered by the inability to efficiently transform P. flocculosa. In this study, two different approaches, target gene replacement by homologous recombination (HR) and CRISPR-Cas9 based genome-editing, were utilized to decipher the role of flocculosin in the biocontrol activity of P.flocculosa. It was possible to alter the production of flocculosin through edition of fat1 by HR, but such mutants displayed abnormal phenotypes and the inability to produce sporidia. Sequencing analyses revealed that transformation by HR led to multiple insertions in the genome explaining the pleiotrophic effects of the approach. On the other hand, CRISPR-Cas9 transformation yielded one mutant that was altered specifically in the proper synthesis of flocculosin. Notwithstanding the loss of flocculosin production, such mutant was phenotypically similar to the wild-type, and when tested for its biocontrol activity against powdery mildew, displayed the same efficacy. These results offer strong evidence that flocculosin-mediated antibiosis is not responsible for the mode of action of P. flocculosa and highlight the potential of CRISPR-Cas9 for functional studies of otherwise difficult-to-transform fungi such as P. flocculosa.

Silicon influences the localization and expression of Phytophthora sojae effectors in interaction with soybean.

In plant-pathogen interactions, expression and localization of effectors in the aqueous apoplastic region play a crucial role in the establishment or suppression of pathogen development. Silicon (Si) has been shown to protect plants in several host-pathogen interactions, but its mode of action remains a source of debate. Its deposition in the apoplastic area of plant cells suggests that it might interfere with receptor-effector recognition. In this study, soybean plants treated or not with Si were inoculated with Phytophthora sojae and differences in the ensuing infection process were assessed through different microscopy techniques, transcript analysis of effector and defense genes, and effector (Avr6) localization through immunolocalization and fluorescence labeling. In plants grown without Si, the results showed the rapid (4 d post-inoculation) host recognition by P. sojae through the development of haustorium-like bodies, followed by expression and release of effectors into the apoplastic region. In contrast, Si treatment resulted in limited pathogen development, and significantly lower expression and presence of Avr6 in the apoplastic region. Based on immunolocalization and quantification of Avr6 through fluorescence labeling, our results suggest that the presence of Si in the apoplast interferes with host recognition and/or limits receptor-effector interactions, which leads to an incompatible interaction.

Rhamnose synthase activity is required for pathogenicity of the vascular wilt fungus Verticillium dahliae.

The initial interaction of a pathogenic fungus with its host is complex and involves numerous metabolic pathways and regulatory proteins. Considerable attention has been devoted to proteins that play a crucial role in these interactions, with an emphasis on so-called effector molecules that are secreted by the invading microbe to establish the symbiosis. However, the contribution of other types of molecules, such as glycans, is less well appreciated. Here, we present a random genetic screen that enabled us to identify 58 novel candidate genes that are involved in the pathogenic potential of the fungal pathogen Verticillium dahliae, which causes vascular wilt diseases in over 200 dicotyledonous plant species, including economically important crops. One of the candidate genes that was identified concerns a putative biosynthetic gene involved in nucleotide sugar precursor formation, as it encodes a putative nucleotide-rhamnose synthase/epimerase-reductase (NRS/ER). This enzyme has homology to bacterial enzymes involved in the biosynthesis of the nucleotide sugar deoxy-thymidine diphosphate (dTDP)-rhamnose, a precursor of L-rhamnose, which has been shown to be required for virulence in several human pathogenic bacteria. Rhamnose is known to be a minor cell wall glycan in fungi and has therefore not been suspected as a crucial molecule in fungal-host interactions. Nevertheless, our study shows that deletion of the VdNRS/ER gene from the V. dahliae genome results in complete loss of pathogenicity on tomato and Nicotiana benthamiana plants, whereas vegetative growth and sporulation are not affected. We demonstrate that VdNRS/ER is a functional enzyme in the biosynthesis of uridine diphosphate (UDP)-rhamnose, and further analysis has revealed that VdNRS/ER deletion strains are impaired in the colonization of tomato roots. Collectively, our results demonstrate that rhamnose, although only a minor cell wall component, is essential for the pathogenicity of V. dahliae.

Verticillium dahliae Sge1 differentially regulates expression of candidate effector genes.

The ascomycete fungus Verticillium dahliae causes vascular wilt diseases in hundreds of dicotyledonous plant species. However, thus far, only few V. dahliae effectors have been identified, and regulators of pathogenicity remain unknown. In this study, we investigated the role of the V. dahliae homolog of Sge1, a transcriptional regulator that was previously implicated in pathogenicity and effector gene expression in Fusarium oxysporum. We show that V. dahliae Sge1 (VdSge1) is required for radial growth and production of asexual conidiospores, because VdSge1 deletion strains display reduced radial growth and reduced conidia production. Furthermore, we show that VdSge1 deletion strains have lost pathogenicity on tomato. Remarkably, VdSge1 is not required for induction of Ave1, the recently identified V. dahliae effector that activates resistance mediated by the Ve1 immune receptor in tomato. Further assessment of the role of VdSge1 in the induction of the nine most highly in-planta-induced genes that encode putative effectors revealed differential activity. Although the expression of one putative effector gene in addition to Ave1 was not affected by VdSge1 deletion, VdSge1 appeared to be required for the expression of six putative effector genes, whereas two of the putative effectors genes were found to be negatively regulated by VdSge1. In conclusion, our data suggest that VdSge1 differentially regulates V. dahliae effector gene expression.

Evidence for functional diversification within a fungal NEP1-like protein family.

In this study, we functionally analyzed the gene family encoding necrosis- and ethylene-inducing-like proteins (NLP) of the vascular wilt pathogen Verticillium dahliae. We show that the composition of the NLP gene family varies little among V. dahliae isolates. The cytotoxic activity of NLP family members of a tomato-pathogenic V. dahliae strain was determined, demonstrating that only two of the seven NLP induced plant cell death. The genes encoding these cytotoxic NLP were found to be induced in V. dahliae upon colonization of tomato. Interestingly, targeted deletion of either of the two genes in V. dahliae significantly compromised virulence on tomato as well as on Arabidopsis plants, whereas deletion of only one of the two genes affected virulence on Nicotiana benthamiana. This could be attributed to differential induction of the two NLP genes in V. dahliae upon N. benthamiana colonization, revealing that the in planta induction of NLP genes varies between plant hosts. Intriguingly, one of the NLP genes appears to also affect vegetative growth and conidiospore production, because the corresponding deletion strain produced significantly fewer conidiospores and developed extensive aerial mycelium. In conclusion, we demonstrate that the expanded V. dahliae NLP family shows functional diversification, revealing not only differential cytotoxicity between family members but also that the cytotoxic NLP play a role in vegetative growth and asexual reproduction in addition to their contribution to virulence.

Tomato immune receptor Ve1 recognizes effector of multiple fungal pathogens uncovered by genome and RNA sequencing.

Fungal plant pathogens secrete effector molecules to establish disease on their hosts, and plants in turn use immune receptors to try to intercept these effectors. The tomato immune receptor Ve1 governs resistance to race 1 strains of the soil-borne vascular wilt fungi Verticillium dahliae and Verticillium albo-atrum, but the corresponding Verticillium effector remained unknown thus far. By high-throughput population genome sequencing, a single 50-Kb sequence stretch was identified that only occurs in race 1 strains, and subsequent transcriptome sequencing of Verticillium-infected Nicotiana benthamiana plants revealed only a single highly expressed ORF in this region, designated Ave1 (for Avirulence on Ve1 tomato). Functional analyses confirmed that Ave1 activates Ve1-mediated resistance and demonstrated that Ave1 markedly contributes to fungal virulence, not only on tomato but also on Arabidopsis. Interestingly, Ave1 is homologous to a widespread family of plant natriuretic peptides. Besides plants, homologous proteins were only found in the bacterial plant pathogen Xanthomonas axonopodis and the plant pathogenic fungi Colletotrichum higginsianum, Cercospora beticola, and Fusarium oxysporum f. sp. lycopersici. The distribution of Ave1 homologs, coincident with the presence of Ave1 within a flexible genomic region, strongly suggests that Verticillium acquired Ave1 from plants through horizontal gene transfer. Remarkably, by transient expression we show that also the Ave1 homologs from F. oxysporum and C. beticola can activate Ve1-mediated resistance. In line with this observation, Ve1 was found to mediate resistance toward F. oxysporum in tomato, showing that this immune receptor is involved in resistance against multiple fungal pathogens.

Random insertional mutagenesis in fungal genomes to identify virulence factors.

Agrobacterium tumefaciens-mediated transformation (ATMT) has become an important tool for functional genomics in fungi. ATMT-based approaches such as random insertional mutagenesis and targeted knockout are widely used for gene functional analysis in plant-pathogen interactions. Here, we describe a protocol for the identification of pathogenicity and virulence genes through random insertional mutagenesis using the fungal wilt pathogen Verticillium dahliae as an example for the protocol.

Comparative genomics yields insights into niche adaptation of plant vascular wilt pathogens.

The vascular wilt fungi Verticillium dahliae and V. albo-atrum infect over 200 plant species, causing billions of dollars in annual crop losses. The characteristic wilt symptoms are a result of colonization and proliferation of the pathogens in the xylem vessels, which undergo fluctuations in osmolarity. To gain insights into the mechanisms that confer the organisms' pathogenicity and enable them to proliferate in the unique ecological niche of the plant vascular system, we sequenced the genomes of V. dahliae and V. albo-atrum and compared them to each other, and to the genome of Fusarium oxysporum, another fungal wilt pathogen. Our analyses identified a set of proteins that are shared among all three wilt pathogens, and present in few other fungal species. One of these is a homolog of a bacterial glucosyltransferase that synthesizes virulence-related osmoregulated periplasmic glucans in bacteria. Pathogenicity tests of the corresponding V. dahliae glucosyltransferase gene deletion mutants indicate that the gene is required for full virulence in the Australian tobacco species Nicotiana benthamiana. Compared to other fungi, the two sequenced Verticillium genomes encode more pectin-degrading enzymes and other carbohydrate-active enzymes, suggesting an extraordinary capacity to degrade plant pectin barricades. The high level of synteny between the two Verticillium assemblies highlighted four flexible genomic islands in V. dahliae that are enriched for transposable elements, and contain duplicated genes and genes that are important in signaling/transcriptional regulation and iron/lipid metabolism. Coupled with an enhanced capacity to degrade plant materials, these genomic islands may contribute to the expanded genetic diversity and virulence of V. dahliae, the primary causal agent of Verticillium wilts. Significantly, our study reveals insights into the genetic mechanisms of niche adaptation of fungal wilt pathogens, advances our understanding of the evolution and development of their pathogenesis, and sheds light on potential avenues for the development of novel disease management strategies to combat destructive wilt diseases.

Generation of a 3D indexed Petunia insertion database for reverse genetics.

BLAST searchable databases containing insertion flanking sequences have revolutionized reverse genetics in plant research. The development of such databases has so far been limited to a small number of model species and normally requires extensive labour input. Here we describe a highly efficient and widely applicable method that we adapted to identify unique transposon-flanking genomic sequences in Petunia. The procedure is based on a multi-dimensional pooling strategy for the collection of DNA samples; up to thousands of different templates are amplified from each of the DNA pools separately, and knowledge of their source is safeguarded by the use of pool-specific (sample) identification tags in one of the amplification primers. All products are combined into a single sample that is subsequently used as a template for unidirectional pyrosequencing. Computational analysis of the clustered sequence output allows automatic assignment of sequences to individual DNA sources. We have amplified and analysed transposon-flanking sequences from a Petunia transposon insertion library of 1000 individuals. Using 30 DNA isolations, 70 PCR reactions and two GS20 sequencing runs, we were able to allocate around 10 000 transposon flanking sequences to specific plants in the library. These sequences have been organized in a database that can be BLAST-searched for insertions into genes of interest. As a proof of concept, we have performed an in silico screen for insertions into members of the NAM/NAC transcription factor family. All in silico-predicted transposon insertions into members of this family could be confirmed in planta.