Copy number variation and transcriptional polymorphisms of Phytophthora sojae RXLR effector genes Avr1a and Avr3a

Bacterial host adaptation through sequence and structural variations of a single type III effector gene

Molecular mechanisms underlying quantitative variations of pathogenicity remain elusive. Here, we identified the Xanthomonas campestris XopJ6 effector that triggers disease resistance in cauliflower and Arabidopsis thaliana. XopJ6 is a close homolog of the Ralstoniapseudosolanacearum PopP2 YopJ family acetyltransferase. XopJ6 is recognized by the RRS1-R/RPS4 NLR pair that integrates a WRKY decoy domain mimicking effector targets. We identified a XopJ6 natural variant carrying a single residue substitution in XopJ6 WRKY-binding site that disrupts interaction with WRKY proteins. This mutation allows XopJ6 to evade immune perception while retaining some XopJ6 virulence functions. Interestingly, xopJ6 resides in a Tn3-family transposon likely contributing to xopJ6 copy number variation (CNV). Using synthetic biology, we demonstrate that xopJ6 CNV tunes pathogen virulence on Arabidopsis through gene dosage-mediated modulation of xopJ6 expression. Together, our findings highlight how sequence and structural genetic variations restricted at a particular effector gene contribute to bacterial host adaptation.

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.

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.

Fusarium-produced vitamin B6 promotes the evasion of soybean resistance by Phytophthora sojae

Plants can be infected by multiple pathogens concurrently in natural systems. However, pathogen-pathogen interactions have rarely been studied. In addition to the oomycete Phytophthora sojae, fungi such as Fusarium spp. also cause soybean root rot. In a 3-year field investigation, we discovered that P. sojae and Fusarium spp. frequently coexisted in diseased soybean roots. Out of 336 P. sojae-soybean-Fusarium combinations, more than 80% aggravated disease. Different Fusarium species all enhanced P. sojae infection when co-inoculated on soybean. Treatment with Fusarium secreted non-proteinaceous metabolites had an effect equal to the direct pathogen co-inoculation. By screening a Fusarium graminearum mutant library, we identified Fusarium promoting factor of Phytophthora sojae infection 1 (Fpp1), encoding a zinc alcohol dehydrogenase. Fpp1 is functionally conserved in Fusarium and contributes to metabolite-mediated infection promotion, in which vitamin B6 (VB6) produced by Fusarium is key. Transcriptional and functional analyses revealed that Fpp1 regulates two VB6 metabolism genes, and VB6 suppresses expression of soybean disease resistance-related genes. These results reveal that co-infection with Fusarium promotes loss of P. sojae resistance in soybean, information that will inform the sustainable use of disease-resistant crop varieties and provide new strategies to control soybean root rot.

Genome analysis of Phytophthora cactorum strains associated with crown- and leather-rot in strawberry

Phytophthora cactorum has two distinct pathotypes that cause crown rot and leather rot in strawberry (Fragaria × ananassa). Strains of the crown rot pathotype can infect both the rhizome (crown) and fruit tissues, while strains of the leather rot pathotype can only infect the fruits of strawberry. The genome of a highly virulent crown rot strain, a low virulent crown rot strain, and three leather rot strains were sequenced using PacBio high fidelity (HiFi) long read sequencing. The reads were de novo assembled to 66.4-67.6 megabases genomes in 178-204 contigs, with N50 values ranging from 892 to 1,036 kilobases. The total number of predicted complete genes in the five P. cactorum genomes ranged from 17,286 to 17,398. Orthology analysis identified a core secretome of 8,238 genes. Comparative genomic analysis revealed differences in the composition of potential virulence effectors, such as putative RxLR and Crinklers, between the crown rot and the leather rot pathotypes. Insertions, deletions, and amino acid substitutions were detected in genes encoding putative elicitors such as beta elicitin and cellulose-binding domain proteins from the leather rot strains compared to the highly virulent crown rot strain, suggesting a potential mechanism for the crown rot strain to escape host recognition during compatible interaction with strawberry. The results presented here highlight several effectors that may facilitate the tissue-specific colonization of P. cactorum in strawberry.

Pathogenesis mechanisms of phytopathogen effectors

Plants commonly face the threat of invasion by a wide variety of pathogens and have developed sophisticated immune mechanisms to defend against infectious diseases. However, successful pathogens have evolved diverse mechanisms to overcome host immunity and cause diseases. Different cell structures and unique cellular organelles carried by plant cells endow plant-specific defense mechanisms, in addition to the common framework of innate immune system shared by both plants and animals. Effectors serve as crucial virulence weapons employed by phytopathogens to disarm the plant immune system and promote infection. Here we summarized the many diverse strategies by which phytopathogen effectors overcome plant defense and prospected future perspectives. This article is categorized under: Infectious Diseases > Molecular and Cellular Physiology.

Molecular mechanisms of Phytophthora sojae avirulence effectors escaping host recognition

Phytophthora sojae is a well-known destructive oomycete pathogen, which causes soybean stem and root rot and poses a serious threat to global food security. Growing soybean cultivars with the appropriate resistance to P. sojae (Rps) genes are the primary management strategy to reduce losses. In most Phytophthora pathosystems, host resistance protein encoded by a specific R gene in the plant recognizes corresponding RxLR effector protein, encoded by an avirulence gene. This gene-for-gene relationship has been exploited to help breeders and agronomists deploy soybean cultivars. To date, 6 Rps genes have been incorporated into commercial soybean germplasm and trigger plant immunity in response to 8 P. sojae avirulence effectors. The incorporation of Rps genes in the soybean population creates selection pressure in favor of novel pathotypes of P. sojae. The 8 avirulence genes evolved to evade the host immune system, driven by genetic selection pressures. Understanding the evading strategies has important reference value for the prevention and control of Phytophthora stem and root rot. This investigation primarily highlights the research on the strategies of P. sojae avirulence effector evasion of host recognition, looking forward to creating durable resistance genes and thereby enabling successful disease management.

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.

Mutations in PpAGO3 Lead to Enhanced Virulence of Phytophthora parasitica by Activation of 25-26 nt sRNA-Associated Effector Genes

Oomycetes represent a unique group of plant pathogens that are destructive to a wide range of crops and natural ecosystems. Phytophthora species possess active small RNA (sRNA) silencing pathways, but little is known about the biological roles of sRNAs and associated factors in pathogenicity. Here we show that an AGO gene, PpAGO3, plays a major role in the regulation of effector genes hence the pathogenicity of Phytophthora parasitica. PpAGO3 was unique among five predicted AGO genes in P. parasitica, showing strong mycelium stage-specific expression. Using the CRISPR-Cas9 technology, we generated PpAGO3^ΔRGG1-3 mutants that carried a deletion of 1, 2, or 3 copies of the N-terminal RGG motif (QRGGYD) but failed to obtain complete knockout mutants, which suggests its vital role in P. parasitica. These mutants showed increased pathogenicity on both Nicotiana benthamiana and Arabidopsis thaliana plants. Transcriptome and sRNA sequencing of PpAGO3^ΔRGG1 and PpAGO3^ΔRGG3 showed that these mutants were differentially accumulated with 25–26 nt sRNAs associated with 70 predicted cytoplasmic effector genes compared to the wild-type, of which 13 exhibited inverse correlation between gene expression and 25–26 nt sRNA accumulation. Transient overexpression of the upregulated RXLR effector genes, PPTG_01869 and PPTG_15425 identified in the mutants PpAGO3^ΔRGG1 and PpAGO3^ΔRGG3, strongly enhanced N. benthamiana susceptibility to P. parasitica. Our results suggest that PpAGO3 functions together with 25–26 nt sRNAs to confer dynamic expression regulation of effector genes in P. parasitica, thereby contributing to infection and pathogenicity of the pathogen.

ATAC-Seq Reveals the Landscape of Open Chromatin and cis-Regulatory Elements in the Phytophthora sojae Genome

Nucleosome-free open chromatin often harbors transcription factor (TF)-binding sites that are associated with active cis-regulatory elements. However, analysis of open chromatin regions has rarely been applied to oomycete or fungal plant pathogens. In this study, we performed the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to identify open chromatin and cis-regulatory elements in Phytophthora sojae at the mycelial stage. We identified 10,389 peaks representing nucleosome-free regions (NFRs). The peaks were enriched in gene-promoter regions and associated with 40% of P. sojae genes; transcription levels were higher for genes with multiple peaks than genes with a single peak and were higher for genes with a single peak than genes without peak. Chromatin accessibility was positively correlated with gene transcription level. Through motif discovery based on NFR peaks in core promoter regions, 25 candidate cis-regulatory motifs with evidence of TF-binding footprints were identified. These motifs exhibited various preferences for location in the promoter region and associations with the transcription level of their target genes, which included some putative pathogenicity-related genes. As the first study revealing the landscape of open chromatin and the correlation between chromatin accessibility and gene transcription level in oomycetes, the results provide a technical reference and data resources for future studies on the regulatory mechanisms of gene transcription.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evasion of plant immunity by microbial pathogens

Plant pathogenic viruses, bacteria, fungi and oomycetes cause destructive diseases in natural habitats and agricultural settings, thereby threatening plant biodiversity and global food security. The capability of plants to sense and respond to microbial infection determines the outcome of plant-microorganism interactions. Host-adapted microbial pathogens exploit various infection strategies to evade or counter plant immunity and eventually establish a replicative niche. Evasion of plant immunity through dampening host recognition or the subsequent immune signalling and defence execution is a crucial infection strategy used by different microbial pathogens to cause diseases, underpinning a substantial obstacle for efficient deployment of host genetic resistance genes for sustainable disease control. In this Review, we discuss current knowledge of the varied strategies microbial pathogens use to evade the complicated network of plant immunity for successful infection. In addition, we discuss how to exploit this knowledge to engineer crop resistance.

Genomic Variations and Mutational Events Associated with Plant-Pathogen Interactions

Simple Summary

Plants, unlike animals, do not have defender cells or an adaptive immune system. Instead, plants rely on each cell’s innate immunity and systemic signals emitted from infection sites. On the other hand, not all plants, even within the same species, are genetically identical, and their genetic backgrounds determine how well they respond to stress factors. Through evolution, plants have acquired various defense mechanisms that play important roles in the never-ending fight between plants and pathogens. Genetic variation in relation to plant disease resistance can thus be contextualized to provide new insights into these defense mechanisms and evolutionary processes that lead to resistance to pathogens. By focusing on genetic variations and mutational events linked with plant–pathogen interactions, the paper explores how genome compartments facilitate plant and pathogen evolutionary processes.

Abstract

Phytopathologists are actively researching the molecular basis of plant–pathogen interactions. The mechanisms of responses to pathogens have been studied extensively in model crop plant species and natural populations. Today, with the rapid expansion of genomic technologies such as DNA sequencing, transcriptomics, proteomics, and metabolomics, as well as the development of new methods and protocols, data analysis, and bioinformatics, it is now possible to assess the role of genetic variation in plant–microbe interactions and to understand the underlying molecular mechanisms of plant defense and microbe pathogenicity with ever-greater resolution and accuracy. Genetic variation is an important force in evolution that enables organisms to survive in stressful environments. Moreover, understanding the role of genetic variation and mutational events is essential for crop breeders to produce improved cultivars. This review focuses on genetic variations and mutational events associated with plant–pathogen interactions and discusses how these genome compartments enhance plants’ and pathogens’ evolutionary processes.

Epigenetics of wheat-rust interaction: an update

The outcome of different host-pathogen interactions is influenced by both genetic and epigenetic systems, which determine the response of plants to pathogens and vice versa. This review highlights key molecular mechanisms and conceptual advances involved in epigenetic research and the progress made in epigenetics of wheat-rust interactions. Epigenetics implies the heritable changes in the way of gene expression as a consequence of the modification of DNA bases, histone proteins, and/or non-coding-RNA biogenesis without disturbing the underlying nucleotide sequence. The changes occurring between DNA and its surrounding chromatin without altering its DNA sequence and leading to significant changes in the genome of any organism are called epigenetic changes. Epigenetics has already been used successfully to explain the mechanism of human pathogens and in the identification of pathogen-induced modifications within various host plants. Wheat rusts are one of the most vital fungal diseases throughout the major wheat-growing areas of the world. The epigenome in plant pathogens causing diseases such as wheat rusts is mysterious. The investigations of host and pathogen epigenetics in the wheat rusts system can offer a piece of suitable evidence for elucidation of the molecular basis of host-pathogen interaction. Besides, the information on the epigenetic regulation of the genes involved in resistance or pathogenicity will provide better insights into the complex resistance signaling pathways and could provide answers to certain key questions, such as whether epigenetic regulation of certain genes is imparting resistance to host in response of certain pathogen elicitors or not. In the last few years, there has been an upsurge in research on the host as well as pathogen epigenetics and its outcome in plant-pathogen interactions. This review summarizes the progress made in the areas related to the epigenetic control of host-pathogen interaction with particular emphasis on wheat rusts.

Mutations in the Promoter and Coding Regions of Avr3a Cause Gain of Virulence of Phytophthora sojae to Rps3a in Soybean

Phytophthora sojae threatens soybean production worldwide, and the cultivation of soybean cultivars carrying Rps genes is the most effective way to control this pathogen. However, DNA mutations in the Avr genes of P. sojae can escape recognization of the corresponding Rps genes, leading to the loss of soybean resistance. In this study, we investigated sequence polymorphism and transcript level of the Avr3a gene in Chinese isolates of P. sojae. Twenty-four mutations resulting in five unique Avr3a alleles were discovered in the Avr3a coding region from 32 P. sojae isolates. The Avr3a transcripts were detectable in the isolates containing Avr3a(I), Avr3a(II), Avr3a(III), and Avr3a(IV) but not in the isolates containing Avr3a(V). Promoter and 5'-UTR sequence analysis revealed eight unique mutations in the promoter region of Avr3a(V), suggesting that the mutations could result in the loss of Avr3a(V) transcription. Virulence tests indicated the isolates containing Avr3a(II) and Avr3a(IV) were virulent, suggesting that the mutations in the coding regions of Avr3a(II) and Avr3a(IV) caused the gain of virulence to Rps3a. Based on DNA mutations of Avr3a in virulent alleles, two SNP markers and one PCR-based marker were developed successfully for detecting the virulence of P. sojae isolates to Rps3a. These findings provide new insights into escape mechanisms of Avr3a and effective support for accurate pathotype identification of P. sojae using molecular methods.

Genomics accelerated isolation of a new stem rust avirulence gene-wheat resistance gene pair

Stem rust caused by the fungus Puccinia graminis f. sp. tritici (Pgt) is a devastating disease of the global staple crop wheat. Although this disease was largely controlled in the latter half of the twentieth century, new virulent strains of Pgt, such as Ug99, have recently evolved^1,2. These strains have caused notable losses worldwide and their continued spread threatens global wheat production. Breeding for disease resistance provides the most cost-effective control of wheat rust diseases³. A number of rust resistance genes have been characterized in wheat and most encode immune receptors of the nucleotide-binding leucine-rich repeat (NLR) class⁴, which recognize pathogen effector proteins known as avirulence (Avr) proteins⁵. However, only two Avr genes have been identified in Pgt so far, AvrSr35 and AvrSr50 (refs. ^6,7), and none in other cereal rusts^8,9. The Sr27 resistance gene was first identified in a wheat line carrying an introgression of the 3R chromosome from Imperial rye¹⁰. Although not deployed widely in wheat, Sr27 is widespread in the artificial crop species Triticosecale (triticale), which is a wheat-rye hybrid and is a host for Pgt^11,12. Sr27 is effective against Ug99 (ref. ¹³) and other recent Pgt strains^14,15. Here, we identify both the Sr27 gene in wheat and the corresponding AvrSr27 gene in Pgt and show that virulence to Sr27 can arise experimentally and in the field through deletion mutations, copy number variation and expression level polymorphisms at the AvrSr27 locus.

Signaling Pathways and Downstream Effectors of Host Innate Immunity in Plants

Phytopathogens, such as biotrophs, hemibiotrophs and necrotrophs, pose serious stress on the development of their host plants, compromising their yields. Plants are in constant interaction with such phytopathogens and hence are vulnerable to their attack. In order to counter these attacks, plants need to develop immunity against them. Consequently, plants have developed strategies of recognizing and countering pathogenesis through pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). Pathogen perception and surveillance is mediated through receptor proteins that trigger signal transduction, initiated in the cytoplasm or at the plasma membrane (PM) surfaces. Plant hosts possess microbe-associated molecular patterns (P/MAMPs), which trigger a complex set of mechanisms through the pattern recognition receptors (PRRs) and resistance (R) genes. These interactions lead to the stimulation of cytoplasmic kinases by many phosphorylating proteins that may also be transcription factors. Furthermore, phytohormones, such as salicylic acid, jasmonic acid and ethylene, are also effective in triggering defense responses. Closure of stomata, limiting the transfer of nutrients through apoplast and symplastic movements, production of antimicrobial compounds, programmed cell death (PCD) are some of the primary defense-related mechanisms. The current article highlights the molecular processes involved in plant innate immunity (PII) and discusses the most recent and plausible scientific interventions that could be useful in augmenting PII.

Comparative genomics of two inbred lines of the potato cyst nematode Globodera rostochiensis reveals disparate effector family-specific diversification patterns

BACKGROUND

Potato cyst nematodes belong to the most harmful pathogens in potato, and durable management of these parasites largely depends on host-plant resistances. These resistances are pathotype specific. The current Globodera rostochiensis pathotype scheme that defines five pathotypes (Ro1 - Ro5) is both fundamentally and practically of limited value. Hence, resistant potato varieties are used worldwide in a poorly informed manner.

RESULTS

We generated two novel reference genomes of G. rostochiensis inbred lines derived from a Ro1 and a Ro5 population. These genome sequences comprise 173 and 189 scaffolds respectively, marking a ≈ 24-fold reduction in fragmentation as compared to the current reference genome. We provide copy number variations for 19 effector families. Four dorsal gland effector families were investigated in more detail. SPRYSECs, known to be implicated in plant defence suppression, constitute by far the most diversified family studied herein with 60 and 99 variants in Ro1 and Ro5 distributed over 18 and 26 scaffolds. In contrast, CLEs, effectors involved in feeding site induction, show strong physical clustering. The 10 and 16 variants cluster on respectively 2 and 1 scaffolds. Given that pathotypes are defined by their effectoromes, we pinpoint the disparate nature of the contributing effector families in terms of sequence diversification and loss and gain of variants.

CONCLUSIONS

Two novel reference genomes allow for nearly complete inventories of effector diversification and physical organisation within and between pathotypes. Combined with insights we provide on effector family-specific diversification patterns, this constitutes a basis for an effectorome-based virulence scheme for this notorious pathogen.

Expression of several Phytophthora cinnamomi putative RxLRs provides evidence for virulence roles in avocado

Phytophthora cinnamomi is a plant pathogenic oomycete that causes Phytophthora root rot of avocado (PRR). Currently, there is a limited understanding of the molecular interactions underlying this disease. Other Phytophthora species employ an arsenal of effector proteins to manipulate host physiology, of which the RxLR effectors contribute to virulence by interfering with host immune responses. The aim of this study was to identify candidate RxLR effectors in P. cinnamomi that play a role in establishing PRR, and to infer possible functions for these effectors. We identified 61 candidate RxLR genes which were expressed during infection of a susceptible avocado rootstock. Several of these genes were present in multiple copies in the P. cinnamomi genome, suggesting that they may contribute to pathogen fitness. Phylogenetic analysis of the manually predicted RxLR protein sequences revealed 12 P. cinnamomi RxLRs that were related to characterised effectors in other Phytophthora spp., providing clues to their functions in planta. Expression profiles of nine more RxLRs point to possible virulence roles in avocado-highlighting a way forward for studies of this interaction. This study represents the first investigation of the expression of P. cinnamomi RxLR genes during the course of avocado infection, and puts forward a pipeline to pinpoint effector genes with potential as virulence determinants, providing a foundation for the future functional characterization of RxLRs that contribute to P. cinnamomi virulence in avocado.

Comparative Analysis of Host-Associated Variation in Phytophthora cactorum

Phytophthora cactorum is often described as a generalist pathogen, with isolates causing disease in a range of plant species. It is the causative agent of two diseases in the cultivated strawberry, crown rot (CR; causing whole plant collapse) and leather rot (LR; affecting the fruit). In the cultivated apple, P. cactorum causes girdling bark rots on the scion (collar rot) and rootstock (crown rot), as well as necrosis of the fine root system (root rot) and fruit rots. We investigated evidence for host specialisation within P. cactorum through comparative genomic analysis of 18 isolates. Whole genome phylogenetic analysis provided genomic support for discrete lineages within P. cactorum, with well-supported non-recombining clades for strawberry CR and apple infecting isolates specialised to strawberry crowns and apple tissue. Isolates of strawberry CR are genetically similar globally, while there is more diversity in apple-infecting isolates. We sought to identify the genetic basis of host specialisation, demonstrating gain and loss of effector complements within the P. cactorum phylogeny, representing putative determinants of host boundaries. Transcriptomic analysis highlighted that those effectors found to be specific to a single host or expanded in the strawberry lineage are amongst those most highly expressed during infection of strawberry and give a wider insight into the key effectors active during strawberry infection. Many effectors that had homologues in other Phytophthoras that have been characterised as avirulence genes were present but not expressed in our tested isolate. Our results highlight several RxLR-containing effectors that warrant further investigation to determine whether they are indeed virulence factors and host-specificity determinants for strawberry and apple. Furthermore, additional work is required to determine whether these effectors are suitable targets to focus attention on for future resistance breeding efforts.

Biotechnological approaches in management of oomycetes diseases

Plant pathogenic oomycetes cause significant impact on agriculture and, therefore, their management is utmost important. Though conventional methods to combat these pathogens (resistance breeding and use of fungicides) are available but these are limited by the availability of resistant cultivars due to evolution of new pathogenic races, development of resistance in the pathogens against agrochemicals and their potential hazardous effects on the environment and human health. This has fuelled a continual search for novel and alternate strategies for management of phytopathogens. The recent advances in oomycetes genome (Phytophthora infestans, P. ramorum, P. sojae, Pythium ultimum, Albugo candida etc.) would further help in understanding host-pathogen interactions essentially needed for designing effective management strategies. In the present communication the novel and alternate strategies for the management of oomycetes diseases are discussed.

A Cas12a-based gene editing system for Phytophthora infestans reveals monoallelic expression of an elicitor

Phytophthora infestans is a destructive pathogen of potato and a model for investigations of oomycete biology. The successful application of a CRISPR gene editing system to P. infestans is so far unreported. We discovered that it is difficult to express CRISPR/Cas9 but not a catalytically inactive form in transformants, suggesting that the active nuclease is toxic. We were able to achieve editing with CRISPR/Cas12a using vectors in which the nuclease and its guide RNA were expressed from a single transcript. Using the elicitor gene Inf1 as a target, we observed editing of one or both alleles in up to 13% of transformants. Editing was more efficient when guide RNA processing relied on the Cas12a direct repeat instead of ribozyme sequences. INF1 protein was not made when both alleles were edited in the same transformant, but surprisingly also when only one allele was altered. We discovered that the isolate used for editing, 1306, exhibited monoallelic expression of Inf1 due to insertion of a copia-like element in the promoter of one allele. The element exhibits features of active retrotransposons, including a target site duplication, long terminal repeats, and an intact polyprotein reading frame. Editing occurred more often on the transcribed allele, presumably due to differences in chromatin structure. The Cas12a system not only provides a tool for modifying genes in P. infestans, but also for other members of the genus by expanding the number of editable sites. Our work also highlights a natural mechanism that remodels oomycete genomes.

Transcriptional Variability Associated With CRISPR-Mediated Gene Replacements at the Phytophthora sojae Avr1b-1 Locus

Transcriptional plasticity enables oomycetes to rapidly adapt to environmental challenges including emerging host resistance. For example, the soybean pathogen Phytophthora sojae can overcome resistance conferred by the host resistance gene Rps1b through natural silencing of its corresponding effector gene, Avr1b-1. With the Phytophthora CRISPR/Cas9 genome editing system, it is possible to generate site-specific knock-out (KO) and knock-in (KI) mutants and to investigate the biological functions of target genes. In this study, the Avr1b-1 gene was deleted from the P. sojae genome using a homology-directed recombination strategy that replaced Avr1b-1 with a gene encoding the fluorescent protein mCherry. As expected, all selected KO transformants gained virulence on Rps1b plants, while infection of plants lacking Rps1b was not compromised. When a sgRNA-resistant version of Avr1b-1 was reintroduced into the Avr1b-1 locus of an Avr1b KO transformant, KI transformants with a well-transcribed Avr1b-1 gene were unable to infect Rps1b-containing soybeans. However, loss of expression of the incoming Avr1b-1 gene was frequently observed in KI transformants, which resulted in these transformants readily infecting Rps1b soybeans. A similar variability in the expression levels of the incoming gene was observed with AVI- or mCherry-tagged Avr1b-1 constructs. Our results suggest that Avr1b-1 may be unusually susceptible to transcriptional variation.

Silent control: microbial plant pathogens evade host immunity without coding sequence changes

Both animals and plants have evolved a robust immune system to surveil and defeat invading pathogenic microbes. Evasion of host immune surveillance is the key for pathogens to initiate successful infection. To evade the host immunity, plant pathogens evolved a variety of strategies such as masking themselves from host immune recognitions, blocking immune signaling transductions, reprogramming immune responses and adapting to immune microenvironmental changes. Gain of new virulence genes, sequence and structural variations enables plant pathogens to evade host immunity through changes in the genetic code. However, recent discoveries demonstrated that variations at the transcriptional, post-transcriptional, post-translational and glycome level enable pathogens to cope with the host immune system without coding sequence changes. The biochemical modification of pathogen associated molecular patterns and silencing of effector genes emerged as potent ways for pathogens to hide from host recognition. Altered processing in mRNA activities provide pathogens with resilience to microenvironment changes. Importantly, these hiding variants are directly or indirectly modulated by catalytic enzymes or enzymatic complexes and cannot be revealed by classical genomics alone. Unveiling these novel host evasion mechanisms in plant pathogens enables us to better understand the nature of plant disease and pinpoints strategies for rational diseases management in global food protection.

Integrated Analysis of Gene Expression, SNP, InDel, and CNV Identifies Candidate Avirulence Genes in Australian Isolates of the Wheat Leaf Rust Pathogen Puccinia triticina

The leaf rust pathogen, Puccinia triticina (Pt), threatens global wheat production. The deployment of leaf rust (Lr) resistance (R) genes in wheat varieties is often followed by the development of matching virulence in Pt due to presumed changes in avirulence (Avr) genes in Pt. Identifying such Avr genes is a crucial step to understand the mechanisms of wheat-rust interactions. This study is the first to develop and apply an integrated framework of gene expression, single nucleotide polymorphism (SNP), insertion/deletion (InDel), and copy number variation (CNV) analysis in a rust fungus and identify candidate avirulence genes. Using a long-read based de novo genome assembly of an isolate of Pt ('Pt104') as the reference, whole-genome resequencing data of 12 Pt pathotypes derived from three lineages Pt104, Pt53, and Pt76 were analyzed. Candidate avirulence genes were identified by correlating virulence profiles with small variants (SNP and InDel) and CNV, and RNA-seq data of an additional three Pt isolates to validate expression of genes encoding secreted proteins (SPs). Out of the annotated 29,043 genes, 2392 genes were selected as SP genes with detectable expression levels. Small variant comparisons between the isolates identified 27-40 candidates and CNV analysis identified 14-31 candidates for each Avr gene, which when combined, yielded the final 40, 64, and 69 candidates for AvrLr1, AvrLr15, and AvrLr24, respectively. Taken together, our results will facilitate future work on experimental validation and cloning of Avr genes. In addition, the integrated framework of data analysis that we have developed and reported provides a more comprehensive approach for Avr gene mining than is currently available.

Insight Into Function and Subcellular Localization of Plasmopara viticola Putative RxLR Effectors

Grapevine downy mildew, caused by oomycete fungus Plasmopara viticola, is one of the most devastating diseases of grapes across the major production regions of the world. Although many putative effector molecules have been identified from this pathogen, the functions of the majority of these are still unknown. In this study, we analyzed the potential function of 26 P. viticola effectors from the highly virulent strain YL. Using transient expression in leaf cells of the tobacco Nicotiana benthamiana, we found that the majority of the effectors could suppress cell death triggered by BAX and INF1, while seven could induce cell death. The subcellular localization of effectors in N. benthamiana was consistent with their localization in cells of Vitis vinifera. Those effectors that localized to the nucleus (17/26) showed a variety of subnuclear localization. Ten of the effectors localized predominantly to the nucleolus, whereas the remaining seven localized to nucleoplasm. Interestingly, five of the effectors were strongly related in sequence and showed identical subcellular localization, but had different functions in N. benthamiana leaves and expression patterns in grapevine in response to P. viticola. This study highlights the potential functional diversity of P. viticola effectors.

Genomic Investigation of the Strawberry Pathogen Phytophthora fragariae Indicates Pathogenicity Is Associated With Transcriptional Variation in Three Key Races

The oomycete Phytophthora fragariae is a highly destructive pathogen of cultivated strawberry (Fragaria × ananassa), causing the root rotting disease, "red core". The host-pathogen interaction has a well described gene-for-gene resistance relationship, but to date neither candidate avirulence nor resistance genes have been identified. We sequenced a set of American, Canadian, and United Kingdom isolates of known race type, along with three representatives of the closely related pathogen of the raspberry (Rubus idaeus), P. rubi, and found a clear population structure, with a high degree of nucleotide divergence seen between some race types and abundant private variation associated with race types 4 and 5. In contrast, between isolates defined as United Kingdom races 1, 2, and 3 (UK1-2-3) there was no evidence of gene loss or gain; or the presence of insertions/deletions (INDELs) or Single Nucleotide Polymorphisms (SNPs) within or in proximity to putative pathogenicity genes could be found associated with race variation. Transcriptomic analysis of representative UK1-2-3 isolates revealed abundant expression variation in key effector family genes associated with pathogen race; however, further long read sequencing did not reveal any long range polymorphisms to be associated with avirulence to race UK2 or UK3 resistance, suggesting either control in trans or other stable forms of epigenetic modification modulating gene expression. This work reveals the combined power of population resequencing to uncover race structure in pathosystems and in planta transcriptomic analysis to identify candidate avirulence genes. This work has implications for the identification of putative avirulence genes in the absence of associated expression data and points toward the need for detailed molecular characterisation of mechanisms of effector regulation and silencing in oomycete plant pathogens.

Discriminant haplotypes of avirulence genes of Phytophthora sojae lead to a molecular assay to predict phenotypes

The soybean-Phytophthora sojae interaction operates on a gene-for-gene relationship, where the product of a resistance gene (Rps) in the host recognizes that of an avirulence gene (Avr) in the pathogen to generate an incompatible reaction. To exploit this form of resistance, one must match with precision the appropriate Rps gene with the corresponding Avr gene. Currently, this association is evaluated by phenotyping assays that are labour-intensive and often imprecise. To circumvent this limitation, we sought to develop a molecular assay that would reveal the avirulence allele of the seven main Avr genes (Avr1a, Avr1b, Avr1c, Avr1d, Avr1k, Avr3a, and Avr6) in order to diagnose with precision the pathotypes of P. sojae isolates. For this purpose, we analysed the genomic regions of these Avr genes in 31 recently sequenced isolates with different virulence profiles and identified discriminant mutations between avirulence and virulence alleles. Specific primers were designed to generate amplicons of a distinct size, and polymerase chain reaction conditions were optimized in a final assay of two parallel runs. When tested on the 31 isolates of known virulence, the assay accurately revealed all avirulence alleles. The test was further assessed and compared to a phenotyping assay on 25 isolates of unknown virulence. The two assays matched in 97% (170/175) of the interactions studied. Interestingly, the sole cases of discrepancy were obtained with Avr3a, which suggests a possible imperfect interaction with Rps3a. This molecular assay offers a powerful and reliable tool to exploit and study with greater precision soybean resistance against P. sojae.

Editing of an effector gene promoter sequence impacts plant-Phytophthora interaction

Pathogen avirulence (Avr) effectors interplay with corresponding plant resistance (R) proteins and activate robust plant immune responses. Although the expression pattern of Avr genes has been tied to their functions for a long time, it is still not clear how Avr gene expression patterns impact plant-microbe interactions. Here, we selected PsAvr3b, which shows a typical effector gene expression pattern from a soybean root pathogen Phytophthora sojae. To modulate gene expression, we engineered PsAvr3b promoter sequences by in situ substitution with promoter sequences from Actin (constitutive expression), PsXEG1 (early expression), and PsNLP1 (later expression) using the CRISPR/Cas9. PsAvr3b driven by different promoters resulted in distinct expression levels across all the tested infection time points. Importantly, those mutants with low PsAvr3b expression successfully colonized soybean plants carrying the cognate R gene Rps3b. To dissect the difference in plant responses to the PsAvr3b expression level, we conducted RNA-sequencing of different infection samples at 24 h postinfection and found soybean immune genes, including a few previously unknown genes that are associated with resistance. Our study highlights that fine-tuning in Avr gene expression impacts the compatibility of plant disease and provides clues to improve crop resistance in disease control management.

Effector gene silencing mediated by histone methylation underpins host adaptation in an oomycete plant pathogen

The relentless adaptability of pathogen populations is a major obstacle to effective disease control measures. Increasing evidence suggests that gene transcriptional polymorphisms are a strategy deployed by pathogens to evade host immunity. However, the underlying mechanisms of transcriptional plasticity remain largely elusive. Here we found that the soybean root rot pathogen Phytophthora sojae evades the soybean Resistance gene Rps1b through transcriptional polymorphisms in the effector gene Avr1b that occur in the absence of any sequence variation. Elevated levels of histone H3 Lysine27 tri-methylation (H3K27me3) were observed at the Avr1b locus in a naturally occurring Avr1b-silenced strain but not in an Avr1b-expressing strain, suggesting a correlation between this epigenetic modification and silencing of Avr1b. To genetically test this hypothesis, we edited the gene, PsSu(z)12, encoding a core subunit of the H3K27me3 methyltransferase complex by using CRISPR/Cas9, and obtained three deletion mutants. H3K27me3 depletion within the Avr1b genomic region correlated with impaired Avr1b gene silencing in these mutants. Importantly, these mutants lost the ability to evade immune recognition by soybeans carrying Rps1b. These data support a model in which pathogen effector transcriptional polymorphisms are associated with changes in chromatin epigenetic marks, highlighting epigenetic variation as a mechanism of pathogen adaptive plasticity.

Haplotype-Phased Genome Assembly of Virulent Phytophthora ramorum Isolate ND886 Facilitated by Long-Read Sequencing Reveals Effector Polymorphisms and Copy Number Variation

Phytophthora ramorum is a destructive pathogen that causes sudden oak death disease. The genome sequence of P. ramorum isolate Pr102 was previously produced, using Sanger reads, and contained 12 Mb of gaps. However, isolate Pr102 had shown reduced aggressiveness and genome abnormalities. In order to produce an improved genome assembly for P. ramorum, we performed long-read sequencing of highly aggressive P. ramorum isolate CDFA1418886 (abbreviated as ND886). We generated a 60.5-Mb assembly of the ND886 genome using the Pacific Biosciences (PacBio) sequencing platform. The assembly includes 302 primary contigs (60.2 Mb) and nine unplaced contigs (265 kb). Additionally, we found a 'highly repetitive' component from the PacBio unassembled unmapped reads containing tandem repeats that are not part of the 60.5-Mb genome. The overall repeat content in the primary assembly was much higher than the Pr102 Sanger version (48 versus 29%), indicating that the long reads have captured repetitive regions effectively. The 302 primary contigs were phased into 345 haplotype blocks and 222,892 phased variants, of which the longest phased block was 1,513,201 bp with 7,265 phased variants. The improved phased assembly facilitated identification of 21 and 25 Crinkler effectors and 393 and 394 RXLR effector genes from two haplotypes. Of these, 24 and 25 RXLR effectors were newly predicted from haplotypes A and B, respectively. In addition, seven new paralogs of effector Avh207 were found in contig 54, not reported earlier. Comparison of the ND886 assembly with Pr102 V1 assembly suggests that several repeat-rich smaller scaffolds within the Pr102 V1 assembly were possibly misassembled; these regions are fully encompassed now in ND886 contigs. Our analysis further reveals that Pr102 is a heterokaryon with multiple nuclear types in the sequences corresponding to contig 10 of ND886 assembly.

The Nuclear-Localized RxLR Effector PvAvh74 From Plasmopara viticola Induces Cell Death and Immunity Responses in Nicotiana benthamiana

Downy mildew is one of the most serious diseases of grapevine (Vitis spp). The causal agent of grapevine downy mildew, Plasmopara viticola, is an obligate biotrophic oomycete. Although oomycete pathogens such as P. viticola are known to secrete RxLR effectors to manipulate host immunity, there have been few studies of the associated mechanisms by which these may act. Here, we show that a candidate P. viticola RxLR effector, PvAvh74, induces cell death in Nicotiana benthamiana leaves. Using agroinfiltration, we found that nuclear localization, two putative N-glycosylation sites, and 427 amino acids of the PvAvh74 carboxyl terminus were necessary for cell-death-inducing activity. Using virus-induced gene silencing (VIGS), we found that PvAvh74-induced cell death in N. benthamiana requires EDS1, NDR1, SGT1, RAR1, and HSP90, but not BAK1. The MAPK cascade components MEK2, WIPK, and SIPK were also involved in PvAvh74-induced cell death in N. benthamiana. Transient expression of PvAvh74 could suppress Phytophthora capsici colonization of N. benthamiana, which suggests that PvAvh74 elicits plant immune responses. Suppression of P. capsici colonization also was dependent on nuclear localization of PvAvh74. Additionally, PvAvh74-triggered cell death could be suppressed by another effector, PvAvh8, from the same isolate. This work provides a framework to further investigate the interactions of PvAvh74 and other RxLR effectors with host immunity.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Stable predictive markers for Phytophthora sojae avirulence genes that impair infection of soybean uncovered by whole genome sequencing of 31 isolates

BACKGROUND

The interaction between oomycete plant pathogen Phytophthora sojae and soybean is characterized by the presence of avirulence (Avr) genes in P. sojae, which encode for effectors that trigger immune responses and resistance in soybean via corresponding resistance genes (Rps). A recent survey highlighted a rapid diversification of P. sojae Avr genes in soybean fields and the need to deploy new Rps genes. However, the full genetic diversity of P. sojae isolates remains complex and dynamic and is mostly characterized on the basis of phenotypic associations with differential soybean lines.

RESULTS

We sequenced the genomes of 31 isolates of P. sojae, representing a large spectrum of the pathotypes found in soybean fields, and compared all the genetic variations associated with seven Avr genes (1a, 1b, 1c, 1d, 1k, 3a, 6) and how the derived haplotypes matched reported phenotypes in 217 interactions. We discovered new variants, copy number variations and some discrepancies with the virulence of previously described isolates with Avr genes, notably with Avr1b and Avr1c. In addition, genomic signatures revealed 11.5% potentially erroneous phenotypes. When these interactions were re-phenotyped, and the Avr genes re-sequenced over time and analyzed for expression, our results showed that genomic signatures alone accurately predicted 99.5% of the interactions.

CONCLUSIONS

This comprehensive genomic analysis of seven Avr genes of P. sojae in a population of 31 isolates highlights that genomic signatures can be used as accurate predictors of phenotypes for compatibility with Rps genes in soybean. Our findings also show that spontaneous mutations, often speculated as a source of aberrant phenotypes, did not occur within the confines of our experiments and further suggest that epigenesis or gene silencing do not account alone for previous discordance between genotypes and phenotypes. Furthermore, on the basis of newly identified virulence patterns within Avr1c, our results offer an explanation why Rps1c has failed more rapidly in the field than the reported information on virulence pathotypes.

Gene expression polymorphism underpins evasion of host immunity in an asexual lineage of the Irish potato famine pathogen

BACKGROUND

Outbreaks caused by asexual lineages of fungal and oomycete pathogens are a continuing threat to crops, wild animals and natural ecosystems (Fisher MC, Henk DA, Briggs CJ, Brownstein JS, Madoff LC, McCraw SL, Gurr SJ, Nature 484:186-194, 2012; Kupferschmidt K, Science 337:636-638, 2012). However, the mechanisms underlying genome evolution and phenotypic plasticity in asexual eukaryotic microbes remain poorly understood (Seidl MF, Thomma BP, BioEssays 36:335-345, 2014). Ever since the 19th century Irish famine, the oomycete Phytophthora infestans has caused recurrent outbreaks on potato and tomato crops that have been primarily caused by the successive rise and migration of pandemic asexual lineages (Goodwin SB, Cohen BA, Fry WE, Proc Natl Acad Sci USA 91:11591-11595, 1994; Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10:e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Cooke DEL, Cano LM, Raffaele S, Bain RA, Cooke LR, Etherington GJ, Deahl KL, Farrer RA, Gilroy EM, Goss EM, et al. PLoS Pathog 8:e1002940, 2012). However, the dynamics of genome evolution within these clonal lineages have not been determined. The objective of this study was to use a comparative genomics and transcriptomics approach to determine the molecular mechanisms that underpin phenotypic variation within a clonal lineage of P. infestans.

RESULTS

Here, we reveal patterns of genomic and gene expression variation within a P. infestans asexual lineage by comparing strains belonging to the South American EC-1 clone that has dominated Andean populations since the 1990s (Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Delgado RA, Monteros-Altamirano AR, Li Y, Visser RGF, van der Lee TAJ, Vosman B, Plant Pathol 62:1081-1088, 2013; Forbes GA, Escobar XC, Ayala CC, Revelo J, Ordonez ME, Fry BA, Doucett K, Fry WE, Phytopathology 87:375-380, 1997; Oyarzun PJ, Pozo A, Ordonez ME, Doucett K, Forbes GA, Phytopathology 88:265-271, 1998). We detected numerous examples of structural variation, nucleotide polymorphisms and loss of heterozygosity within the EC-1 clone. Remarkably, 17 genes are not expressed in one of the two EC-1 isolates despite apparent absence of sequence polymorphisms. Among these, silencing of an effector gene was associated with evasion of disease resistance conferred by a potato immune receptor.

CONCLUSIONS

Our findings highlight the molecular changes underpinning the exceptional genetic and phenotypic plasticity associated with host adaptation in a pandemic clonal lineage of a eukaryotic plant pathogen. We observed that the asexual P. infestans lineage EC-1 can exhibit phenotypic plasticity in the absence of apparent genetic mutations resulting in virulence on a potato carrying the Rpi-vnt1.1 gene. Such variant alleles may be epialleles that arose through epigenetic changes in the underlying genes.

Expression of the Potato Late Blight Resistance Gene Rpi-phu1 and Phytophthora infestans Effectors in the Compatible and Incompatible Interactions in Potato

This study describes late blight resistance of potato breeding lines resulting from crosses between cultivar 'Sárpo Mira' and Rpi-phu1 gene donors. The progeny is investigated for the presence of Rpi-Smira1 and Rpi-phu1 resistance (R) genes. Interestingly, in detached-leaflet tests, plants with both R genes withstood the infection of the Phytophthora infestans isolate virulent to each gene separately, due to either interaction of these genes or the presence of additional resistance loci. The interaction was studied further in three chosen breeding lines on the transcriptional level. The Rpi-phu1 expression, measured over 5 days, revealed different patterns depending on the outcome of the interaction with P. infestans: it increased in infected plants whereas it remained low and stable when infection was unsuccessful. The expression patterns of P. infestans effectors Avr-vnt1, AvrSmira1, and Avr8, recognized by the Rpi-phu1, Rpi-Smira1, and Rpi-Smira2 genes, respectively, were evaluated in the same experimental setup. This is the first report that the Avr-vnt1 effector expression is not switched off permanently in virulent isolates to avoid recognition by an R protein but can reappear in a postbiotrophic phase and is present constantly when infecting plants without the corresponding R gene. Both a plant and a pathogen can react to the other interacting side by changing the transcript accumulation of R genes or effectors.

Genome sequence of Plasmopara viticola and insight into the pathogenic mechanism

Plasmopara viticola causes downy mildew disease of grapevine which is one of the most devastating diseases of viticulture worldwide. Here we report a 101.3 Mb whole genome sequence of P. viticola isolate 'JL-7-2' obtained by a combination of Illumina and PacBio sequencing technologies. The P. viticola genome contains 17,014 putative protein-coding genes and has ~26% repetitive sequences. A total of 1,301 putative secreted proteins, including 100 putative RXLR effectors and 90 CRN effectors were identified in this genome. In the secretome, 261 potential pathogenicity genes and 95 carbohydrate-active enzymes were predicted. Transcriptional analysis revealed that most of the RXLR effectors, pathogenicity genes and carbohydrate-active enzymes were significantly up-regulated during infection. Comparative genomic analysis revealed that P. viticola evolved independently from the Arabidopsis downy mildew pathogen Hyaloperonospora arabidopsidis. The availability of the P. viticola genome provides a valuable resource not only for comparative genomic analysis and evolutionary studies among oomycetes, but also enhance our knowledge on the mechanism of interactions between this biotrophic pathogen and its host.

Genome sequence of Plasmopara viticola and insight into the pathogenic mechanism

Plasmopara viticola causes downy mildew disease of grapevine which is one of the most devastating diseases of viticulture worldwide. Here we report a 101.3 Mb whole genome sequence of P. viticola isolate 'JL-7-2' obtained by a combination of Illumina and PacBio sequencing technologies. The P. viticola genome contains 17,014 putative protein-coding genes and has ~26% repetitive sequences. A total of 1,301 putative secreted proteins, including 100 putative RXLR effectors and 90 CRN effectors were identified in this genome. In the secretome, 261 potential pathogenicity genes and 95 carbohydrate-active enzymes were predicted. Transcriptional analysis revealed that most of the RXLR effectors, pathogenicity genes and carbohydrate-active enzymes were significantly up-regulated during infection. Comparative genomic analysis revealed that P. viticola evolved independently from the Arabidopsis downy mildew pathogen Hyaloperonospora arabidopsidis. The availability of the P. viticola genome provides a valuable resource not only for comparative genomic analysis and evolutionary studies among oomycetes, but also enhance our knowledge on the mechanism of interactions between this biotrophic pathogen and its host.

A candidate RxLR effector from Plasmopara viticola can elicit immune responses in Nicotiana benthamiana

BACKGROUND

Diverse plant pathogens deliver effectors into plant cells to alter host processes. Oomycete pathogen encodes a large number of putative RxLR effectors which are likely to play a role in manipulating plant defense responses. The secretome of Plasmopara viticola (downy mildew of grapevine) contains at least 162 candidate RxLR effectors discovered in our recent studies, but their roles in infection and pathogenicity remain to be determined. Here, we characterize in depth one of the putative RxLR effectors, PvRxLR16, which has been reported to induce cell death in Nicotiana benthamiana in our previous study.

RESULTS

The nuclear localization, W/Y/L motifs, and a putative N-glycosylation site in C-terminal of PvRxLR16 were essential for cell death-inducing activity. Suppressor of G-two allele of Skp1 (SGT1), heat shock protein 90 (HSP90) and required for Mla12 resistance (RAR1), but not somatic embryogenesis receptor-like kinase (SERK3), were required for the cell death response triggered by PvRxLR16 in N. benthamiana. Some mitogen-activated protein kinases and transcription factors were also involved in the perception of PvRxLR16 by N. benthamiana. PvRxLR16 could also significantly enhance plant resistance to Phytophthora capsici and the nuclear localization was required for this ability. However, some other PvRxLR effectors could suppress defense responses and disease resistance induced by PvRxLR16, suggesting that it may not trigger host cell death or immune responses during physiological infection under natural conditions.

CONCLUSION

These data demonstrate that PvRxLR16 may be recognized by endogenous proteins in nucleus to trigger immune responses in N. benthamiana, which in turn can be suppressed by other PvRxLR effectors.

Copy number variation contributes to cryptic genetic variation in outbreak lineages of Cryptococcus gattii from the North American Pacific Northwest

Background

Copy number variants (CNVs) are a class of structural variants (SVs) and are defined as fragments of DNA that are present at variable copy number in comparison with a reference genome. Recent advances in bioinformatics methodologies and sequencing technologies have enabled the high-resolution quantification of genome-wide CNVs. In pathogenic fungi SVs have been shown to alter gene expression, influence host specificity, and drive fungicide resistance, but little attention has focused specifically on CNVs. Using publicly available sequencing data, we identified 90 isolates across 212 Cryptococcus gattii genomes that belong to the VGII subgroups responsible for the recent deadly outbreaks in the North American Pacific Northwest. We generated CNV profiles for each sample to investigate the prevalence and function of CNV in C. gattii.

Results

We identified eight genetic clusters among publicly available Illumina whole genome sequence data from 212 C. gattii isolates through population structure analysis. Three clusters represent the VGIIa, VGIIb, and VGIIc subgroups from the North American Pacific Northwest. CNV was bioinformatically predicted and affected ~300–400 Kilobases (Kb) of the C. gattii VGII subgroup genomes. Sixty-seven loci, encompassing 58 genes, showed highly divergent patterns of copy number variation between VGII subgroups. Analysis of PFam domains within divergent CN variable genes revealed enrichment of protein domains associated with transport, cell wall organization and external encapsulating structure.

Conclusions

CNVs may contribute to pathological and phenotypic differences observed between the C. gattii VGIIa, VGIIb, and VGIIc subpopulations. Genes overlapping with population differentiated CNVs were enriched for several virulence related functional terms. These results uncover novel candidate genes to examine the genetic and functional underpinnings of C. gattii pathogenicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3044-0) contains supplementary material, which is available to authorized users.

Molecular Soybean-Pathogen Interactions

Soybean hosts a wide variety of pathogens that cause significant yield losses. The importance of soybean as a major oilseed crop has led to research focused on its interactions with pathogens, such as Soybean mosaic virus, Pseudomonas syringae, Phytophthora sojae, Phakopsora pachyrhizi, and Heterodera glycines. Pioneering work on soybean's interactions with these organisms, which represent the five major pathogen groups (viruses, bacteria, oomycetes, fungi, and nematodes), has contributed to our understanding of the molecular mechanisms underlying virulence and immunity. These mechanisms involve conserved and unique features that validate the need for research in both soybean and homologous model systems. In this review, we discuss identification of effectors and their functions as well as resistance gene-mediated recognition and signaling. We also point out areas in which model systems and recent advances in resources and tools have provided opportunities to gain deeper insights into soybean-pathogen interactions.

Studying the Mechanism of Plasmopara viticola RxLR Effectors on Suppressing Plant Immunity

The RxLR effector family, produced by oomycete pathogens, may manipulate host physiological and biochemical events inside host cells. A group of putative RxLR effectors from Plasmopara viticola have been recently identified by RNA-Seq analysis in our lab. However, their roles in pathogenesis are poorly understood. In this study, we attempted to characterize 23 PvRxLR effector candidates identified from a P. viticola isolate “ZJ-1-1.” During host infection stages, expression patterns of the effector genes were varied and could be categorized into four different groups. By using transient expression assays in Nicotiana benthamiana, we found that 17 of these effector candidates fully suppressed programmed cell death elicited by a range of cell death-inducing proteins, including BAX, INF1, PsCRN63, PsojNIP, PvRxLR16 and R3a/Avr3a. We also discovered that all these PvRxLRs could target the plant cell nucleus, except for PvRxLR55 that localized to the membrane. Furthermore, we identified a single effector, PvRxLR28, that showed the highest expression level at 6 hpi. Functional analysis revealed that PvRxLR28 could significantly enhance susceptibilities of grapevine and tobacco to pathogens. These results suggest that most P. viticola effectors tested in this study may act as broad suppressors of cell death to manipulate immunity in plant.

Elucidating the Role of Effectors in Plant-Fungal Interactions: Progress and Challenges

Pathogenic fungi have diverse growth lifestyles that support fungal colonization on plants. Successful colonization and infection for all lifestyles depends upon the ability to modify living host plants to sequester the necessary nutrients required for growth and reproduction. Secretion of virulence determinants referred to as “effectors” is assumed to be the key governing factor that determines host infection and colonization. Effector proteins are capable of suppressing plant defense responses and alter plant physiology to accommodate fungal invaders. This review focuses on effector molecules of biotrophic and hemibiotrophic plant pathogenic fungi, and the mechanism required for the release and uptake of effector molecules by the fungi and plant cells, respectively. We also place emphasis on the discovery of effectors, difficulties associated with predicting the effector repertoire, and fungal genomic features that have helped promote effector diversity leading to fungal evolution. We discuss the role of specific effectors found in biotrophic and hemibiotrophic fungi and examine how CRISPR/Cas9 technology may provide a new avenue for accelerating our ability in the discovery of fungal effector function.

Oomycete interactions with plants: infection strategies and resistance principles

The Oomycota include many economically significant microbial pathogens of crop species. Understanding the mechanisms by which oomycetes infect plants and identifying methods to provide durable resistance are major research goals. Over the last few years, many elicitors that trigger plant immunity have been identified, as well as host genes that mediate susceptibility to oomycete pathogens. The mechanisms behind these processes have subsequently been investigated and many new discoveries made, marking a period of exciting research in the oomycete pathology field. This review provides an introduction to our current knowledge of the pathogenic mechanisms used by oomycetes, including elicitors and effectors, plus an overview of the major principles of host resistance: the established R gene hypothesis and the more recently defined susceptibility (S) gene model. Future directions for development of oomycete-resistant plants are discussed, along with ways that recent discoveries in the field of oomycete-plant interactions are generating novel means of studying how pathogen and symbiont colonizations overlap.

Strain Specific Factors Control Effector Gene Silencing in Phytophthora sojae

The Phytophthora sojae avirulence gene Avr3a encodes an effector that is capable of triggering immunity on soybean plants carrying the resistance gene Rps3a. P. sojae strains that express Avr3a are avirulent to Rps3a plants, while strains that do not are virulent. To study the inheritance of Avr3a expression and virulence towards Rps3a, genetic crosses and self-fertilizations were performed. A cross between P. sojae strains ACR10 X P7076 causes transgenerational gene silencing of Avr3a allele, and this effect is meiotically stable up to the F₅ generation. However, test-crosses of F₁ progeny (ACR10 X P7076) with strain P6497 result in the release of silencing of Avr3a. Expression of Avr3a in the progeny is variable and correlates with the phenotypic penetrance of the avirulence trait. The F₁ progeny from a direct cross of P6497 X ACR10 segregate for inheritance for Avr3a expression, a result that could not be explained by parental imprinting or heterozygosity. Analysis of small RNA arising from the Avr3a gene sequence in the parental strains and hybrid progeny suggests that the presence of small RNA is necessary but not sufficient for gene silencing. Overall, we conclude that inheritance of the Avr3a gene silenced phenotype relies on factors that are variable among P. sojae strains.

Recent Progress in RXLR Effector Research

Some of the most devastating oomycete pathogens deploy effector proteins, with the signature amino acid motif RXLR, that enter plant cells to promote virulence. Research on the function and evolution of RXLR effectors has been very active over the decade that has transpired since their discovery. Comparative genomics indicate that RXLR genes play a major role in virulence for Phytophthora and downy mildew species. Importantly, gene-for-gene resistance against these oomycete lineages is based on recognition of RXLR proteins. Comparative genomics have revealed several mechanisms through which this resistance can be broken, most notably involving epigenetic control of RXLR gene expression. Structural studies have revealed a core fold that is present in the majority of RXLR proteins, providing a foundation for detailed mechanistic understanding of virulence and avirulence functions. Finally, functional studies have demonstrated that suppression of host immunity is a major function for RXLR proteins. Host protein targets are being identified in a variety of plant cell compartments. Some targets comprise hubs that are also manipulated by bacteria and fungi, thereby revealing key points of vulnerability in the plant immune network.

A Recent Expansion of the RXLR Effector Gene Avrblb2 Is Maintained in Global Populations of Phytophthora infestans Indicating Different Contributions to Virulence

The introgression of disease resistance (R) genes encoding immunoreceptors with broad-spectrum recognition into cultivated potato appears to be the most promising approach to achieve sustainable management of late blight caused by the oomycete pathogen Phytophthora infestans. Rpi-blb2 from Solanum bulbocastanum shows great potential for use in agriculture based on preliminary potato disease trials. Rpi-blb2 confers immunity by recognizing the P. infestans avirulence effector protein AVRblb2 after it is translocated inside the plant cell. This effector belongs to the RXLR class of effectors and is under strong positive selection. Structure-function analyses revealed a key polymorphic amino acid (position 69) in AVRblb2 effector that is critical for activation of Rpi-blb2. In this study, we reconstructed the evolutionary history of the Avrblb2 gene family and further characterized its genetic structure in worldwide populations. Our data indicate that Avrblb2 evolved as a single-copy gene in a putative ancestral species of P. infestans and has recently expanded in the Phytophthora spp. that infect solanaceous hosts. As a consequence, at least four variants of AVRblb2 arose in P. infestans. One of these variants, with a Phe residue at position 69, evades recognition by the cognate resistance gene. Surprisingly, all Avrblb2 variants are maintained in pathogen populations. This suggests a potential benefit for the pathogen in preserving duplicated versions of AVRblb2, possibly because the variants may have different contributions to pathogen fitness in a diversified solanaceous host environment.

Adaptive genomic structural variation in the grape powdery mildew pathogen, Erysiphe necator

BACKGROUND

Powdery mildew, caused by the obligate biotrophic fungus Erysiphe necator, is an economically important disease of grapevines worldwide. Large quantities of fungicides are used for its control, accelerating the incidence of fungicide-resistance. Copy number variations (CNVs) are unbalanced changes in the structure of the genome that have been associated with complex traits. In addition to providing the first description of the large and highly repetitive genome of E. necator, this study describes the impact of genomic structural variation on fungicide resistance in Erysiphe necator.

RESULTS

A shotgun approach was applied to sequence and assemble the genome of five E. necator isolates, and RNA-seq and comparative genomics were used to predict and annotate protein-coding genes. Our results show that the E. necator genome is exceptionally large and repetitive and suggest that transposable elements are responsible for genome expansion. Frequent structural variations were found between isolates and included copy number variation in EnCYP51, the target of the commonly used sterol demethylase inhibitor (DMI) fungicides. A panel of 89 additional E. necator isolates collected from diverse vineyard sites was screened for copy number variation in the EnCYP51 gene and for presence/absence of a point mutation (Y136F) known to result in higher fungicide tolerance. We show that an increase in EnCYP51 copy number is significantly more likely to be detected in isolates collected from fungicide-treated vineyards. Increased EnCYP51 copy numbers were detected with the Y136F allele, suggesting that an increase in copy number becomes advantageous only after the fungicide-tolerant allele is acquired. We also show that EnCYP51 copy number influences expression in a gene-dose dependent manner and correlates with fungal growth in the presence of a DMI fungicide.

CONCLUSIONS

Taken together our results show that CNV can be adaptive in the development of resistance to fungicides by providing increasing quantitative protection in a gene-dosage dependent manner. The results of this work not only demonstrate the effectiveness of using genomics to dissect complex traits in organisms with very limited molecular information, but also may have broader implications for understanding genomic dynamics in response to strong selective pressure in other pathogens with similar genome architectures.

Transcriptomic analysis of the phytopathogenic oomycete Phytophthora cactorum provides insights into infection-related effectors

Background

Phytophthora cactorum, a hemibiotrophic oomycete pathogen, can cause destructive diseases on numerous crops worldwide, leading to essential economic losses every year. However, little has been known about its molecular pathogenicity mechanisms. To gain insight into its repertoire of effectors, the P. cactorum transcriptome was investigated using Illumina RNA-seq.

Results

We first demonstrated an in vitro inoculation method that can be used to mimic natural cyst germination on host plants. Over 28 million cDNA reads were obtained for five life cycle stages (mycelium, sporangium, zoospore, cyst and germinating cyst) and de novo assembled into 21,662 unique genes. By comparisons with 11 public databases, 88.99% of the unique genes were annotated, including 15,845 mapped to the gene models of the annotated relative Phytophthora infestans. Using TribeMCL, 5,538 gene families conserved across P. cactorum and other three completely sequenced Phytophthora pathogen species were determined. In silico analyses revealed that 620 P. cactorum effector homologues including 94 RXLR effector candidates matched known or putative virulence genes in other oomycetes. About half of the RXLR effector candidates were predicted to share a conserved structure unit, termed the WY-domain fold. A subset of the effector genes were checked and validated by PCR amplification. Transcriptional experiments indicated that effector genes were differentially expressed during the life cycle and host infection stages of P. cactorum. Ectopic expression in Nicotiana benthamiana revealed that RXLR, elicitin and NLP effectors can trigger plant cell death. These effectors are highly conserved across oomycete species. Single nucleotide polymorphisms for RXLR effectors were detected in a collection of P. cactorum isolates from different countries and hosts.

Conclusions

This study demonstrates the comprehensive sequencing, de novo assembly, and analyses of the transcriptome of P. cactorum life cycle stages. In the absence of genome sequence, transcriptome data is important for infection-related gene discovery in P. cactorum, as demonstrated here for the effector genes. The first look at the transcriptome and effector arsenal of P. cactorum provides valuable data to elucidate the pathogenicity basis of this broad-host-range pathogen.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-980) contains supplementary material, which is available to authorized users.