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

Draft Genome Sequences of the Oomycete Pilasporangium apinafurcum Strains JCM 30513 and JCM 30514, Formerly Classified as Pythium apinafurcum

Pilasporangium apinafurcum, formerly classified as Pythium apinafurcum, is a unique oomycete that infects plants asymptomatically. Here, we present the draft genome sequences of two variants of P. apinafurcum, JCM 30513 and JCM 30514, isolated from uncultivated field soil in Wakayama Prefecture, Japan.

Complete Mitochondrial Genome of Phytophthora nicotianae and Identification of Molecular Markers for the Oomycetes

Phytophthora nicotianae is one of the most destructive plant pathogens affecting a variety of plants, causing black shank of tobacco, among several other devastating diseases. Herein, we assembled the mitochondrial genome of P. nicotianae and analyzed its gene content and genome structure, performed comparative mitochondrial genomics analysis, and assessed phylogenetic relationships among oomycetes species. The circular mitogenome is 37,561 bp long, with 38 protein-coding genes, 25 transfer RNA (tRNA) genes, and 2 ribosomal RNA genes (rrnl and rrns). The mitochondrial genome showed a biased A/T usage versus G/C. The overall gene content and size of the P. nicotianae mitogenome are identical to those of other published Phytophthora mitogenomes. Interestingly, collinearity analysis using an existing ∼10 k inversion region (including 11 genes and 8 tRNAs) revealed that Phytophthora andina, Phytophthora infestans, Phytophthora mirabilis, Phytophthora ipomoeae, and Phytophthora phaseoli differed from Phytophthora nicotianae, Phytophthora sojae, Phytophthora ramorum, and Phytophthora polonica. Moreover, inverted repeat regions were found to be absent among species of the Peronosporales when compared with species from the Pythiales and Saprolegniales. A phylogenomic investigation based on 29 protein-coding genes demonstrated that Phytophthora is monophyletic, and placed P. nicotianae close to the clade including P. mirabilis, P. ipomoeae, P. andina, P. infestans, and P. phaseoli. Furthermore, we discovered six new candidate DNA molecular markers (rpl6, atp8, nad11, rps2, rps3, and rps4) based on these mitogenomes that would be suitable for species identification in the oomycetes, which have the same identification level as the whole mitogenome and ribosomal DNA sequences. These new molecular markers can not only provide a quick preview of the species without mitogenome information, but will also help to gain better understanding of the oomycetes pathogens and developing treatment or monitoring strategies.

Intrinsic disorder is a common structural characteristic of RxLR effectors in oomycete pathogens

Intrinsic disorder is common in nature and has been studied to play important biological roles in bacterial effectors. However, disorder in oomycete RxLR effectors has not been investigated previously and the roles are unknown. Our results of PONDR VL-XT disorder analysis showed that predicted oomycete RxLR effectors were significantly more disordered than other effectors and secretome. The distribution of disorder content presented preference that RxLR-dEER regions were enriched in disordered residues, suggesting potential role of disorder in effector translocation. In contrast, the disorder content was depleted in the C-terminal regions, especially for W/Y/L motifs. We also found that around 42 % of putative RxLR proteins were predicted to contain at least one α-helix-forming molecular recognition feature (α-MoRF), and most α-MoRFs were located in the C-terminal regions. Furthermore, both of the disorder mutants of PsAvh18 and PcAvh207 lost the cell death-inducing activity, indicating the potential important role of disordered structure in RxLR effector function. Overall, these results demonstrate that intrinsic disorder is a common characteristic of oomycete RxLR proteins, and we postulate that such structure feature may be important for effector translocation or function. This study extends our understanding of RxLR effectors in protein structures, and opens up new directions to explore novel mechanisms of oomycete RxLR effectors.

Finding a partner in the ocean: molecular and evolutionary bases of the response to sexual cues in a planktonic diatom

Microalgae play a major role as primary producers in aquatic ecosystems. Cell signalling regulates their interactions with the environment and other organisms, yet this process in phytoplankton is poorly defined. Using the marine planktonic diatom Pseudo-nitzschia multistriata, we investigated the cell response to cues released during sexual reproduction, an event that demands strong regulatory mechanisms and impacts on population dynamics. We sequenced the genome of P. multistriata and performed phylogenomic and transcriptomic analyses, which allowed the definition of gene gains and losses, horizontal gene transfers, conservation and evolutionary rate of sex-related genes. We also identified a small number of conserved noncoding elements. Sexual reproduction impacted on cell cycle progression and induced an asymmetric response of the opposite mating types. G protein-coupled receptors and cyclic guanosine monophosphate (cGMP) are implicated in the response to sexual cues, which overall entails a modulation of cell cycle, meiosis-related and nutrient transporter genes, suggesting a fine control of nutrient uptake even under nutrient-replete conditions. The controllable life cycle and the genome sequence of P. multistriata allow the reconstruction of changes occurring in diatoms in a key phase of their life cycle, providing hints on the evolution and putative function of their genes and empowering studies on sexual reproduction.

Cloning and functional analysis of succinate dehydrogenase gene PsSDHA in Phytophthora sojae

Succinate dehydrogenase (SDH) is one of the key enzymes of the tricarboxylic acid cycle (TCA cycle) and a proven target of fungicides for true fungi. To explore the roles of the SDHA gene in Phytophthora sojae, we first cloned PsSDHA to construct the PsSDHA silenced expression vector pHAM34-PsSDHA, and then utilized PEG to mediate the P. sojae protoplast transformation experiment. Through transformation screening, we obtained the silenced mutants A1 and A3, which have significant suppressive effect. Further study showed that the hyphae of the silenced mutant strains were shorter and more bifurcated; the growth of the silenced mutants was clearly inhibited in 10% V8 agar medium containing sodium chloride (NaCl), hydrogen peroxide (H2O2) or Congo Red, respectively. The pathogenicity of the silenced mutants was significantly reduced compared with the wild-type strain and the mock. The results could help us better to understand the position and function of SDH in P. sojae and provide a proven target of fungicides for the oomycete.

Effectors of Filamentous Plant Pathogens: Commonalities amid Diversity

Fungi and oomycetes are filamentous microorganisms that include a diversity of highly developed pathogens of plants. These are sophisticated modulators of plant processes that secrete an arsenal of effector proteins to target multiple host cell compartments and enable parasitic infection. Genome sequencing revealed complex catalogues of effectors of filamentous pathogens, with some species harboring hundreds of effector genes. Although a large fraction of these effector genes encode secreted proteins with weak or no sequence similarity to known proteins, structural studies have revealed unexpected similarities amid the diversity. This article reviews progress in our understanding of effector structure and function in light of these new insights. We conclude that there is emerging evidence for multiple pathways of evolution of effectors of filamentous plant pathogens but that some families have probably expanded from a common ancestor by duplication and diversification. Conserved folds, such as the oomycete WY and the fungal MAX domains, are not predictive of the precise function of the effectors but serve as a chassis to support protein structural integrity while providing enough plasticity for the effectors to bind different host proteins and evolve unrelated activities inside host cells. Further effector evolution and diversification arise via short linear motifs, domain integration and duplications, and oligomerization.

Metabolic Diversity and Novelties in the Oomycetes

The eukaryotic microbes called oomycetes include many important saprophytes and pathogens, with the latter exhibiting necrotrophy, biotrophy, or obligate biotrophy. Understanding oomycete metabolism is fundamental to understanding these lifestyles. Genome mining and biochemical studies have shown that oomycetes, which belong to the kingdom Stramenopila, secrete suites of carbohydrate- and protein-degrading enzymes adapted to their environmental niches and produce unusual lipids and energy storage compounds. Despite having limited secondary metabolism, many oomycetes make chemicals for communicating within their species or with their hosts. Horizontal and endosymbiotic gene transfer events have diversified oomycete metabolism, resulting in biochemical pathways that often depart from standard textbook descriptions by amalgamating enzymes from multiple sources. Gene fusions and duplications have further shaped the composition and expression of the enzymes. Current research is helping us learn how oomycetes interact with host and environment, understand eukaryotic diversity and evolution, and identify targets for drugs and crop protection chemicals.

Network and role analysis of autophagy in Phytophthora sojae

Autophagy is an evolutionarily conserved mechanism in eukaryotes with roles in development and the virulence of plant fungal pathogens. However, few reports on autophagy in oomycete species have been published. Here, we identified 26 autophagy-related genes (ATGs) belonging to 20 different groups in Phytophthora sojae using a genome-wide survey, and core ATGs in oomycetes were used to construct a preliminary autophagy pathway model. Expression profile analysis revealed that these ATGs are broadly expressed and that the majority of them significantly increase during infection stages, suggesting a central role for autophagy in virulence. Autophagy in P. sojae was detected using a GFP-PsAtg8 fusion protein and the fluorescent dye MDC during rapamycin and starvation treatment. In addition, autophagy was significantly induced during sporangium formation and cyst germination. Silencing PsAtg6a in P. sojae significantly reduced sporulation and pathogenicity. Furthermore, a PsAtg6a-silenced strain showed haustorial formation defects. These results suggested that autophagy might play essential roles in both the development and infection mechanism of P. sojae.

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.

Phylogenomic Reconstruction of the Oomycete Phylogeny Derived from 37 Genomes

The oomycetes are a class of eukaryotes and include ecologically significant animal and plant pathogens. Single-gene and multigene phylogenetic studies of individual oomycete genera and of members of the larger classes have resulted in conflicting conclusions concerning interspecies relationships among these species, particularly for the Phytophthora genus. The onset of next-generation sequencing techniques now means that a wealth of oomycete genomic data is available. For the first time, we have used genome-scale phylogenetic methods to resolve oomycete phylogenetic relationships. We used supertree methods to generate single-gene and multigene species phylogenies. Overall, our supertree analyses utilized phylogenetic data from 8,355 oomycete gene families. We have also complemented our analyses with superalignment phylogenies derived from 131 single-copy ubiquitous gene families. Our results show that a genome-scale approach to oomycete phylogeny resolves oomycete classes and clades. Our analysis represents an important first step in large-scale phylogenomic analysis of the oomycetes.

The oomycetes are a class of microscopic, filamentous eukaryotes within the Stramenopiles-Alveolata-Rhizaria (SAR) supergroup which includes ecologically significant animal and plant pathogens, most infamously the causative agent of potato blight Phytophthora infestans. Single-gene and concatenated phylogenetic studies both of individual oomycete genera and of members of the larger class have resulted in conflicting conclusions concerning species phylogenies within the oomycetes, particularly for the large Phytophthora genus. Genome-scale phylogenetic studies have successfully resolved many eukaryotic relationships by using supertree methods, which combine large numbers of potentially disparate trees to determine evolutionary relationships that cannot be inferred from individual phylogenies alone. With a sufficient amount of genomic data now available, we have undertaken the first whole-genome phylogenetic analysis of the oomycetes using data from 37 oomycete species and 6 SAR species. In our analysis, we used established supertree methods to generate phylogenies from 8,355 homologous oomycete and SAR gene families and have complemented those analyses with both phylogenomic network and concatenated supermatrix analyses. Our results show that a genome-scale approach to oomycete phylogeny resolves oomycete classes and individual clades within the problematic Phytophthora genus. Support for the resolution of the inferred relationships between individual Phytophthora clades varies depending on the methodology used. Our analysis represents an important first step in large-scale phylogenomic analysis of the oomycetes.

IMPORTANCE The oomycetes are a class of eukaryotes and include ecologically significant animal and plant pathogens. Single-gene and multigene phylogenetic studies of individual oomycete genera and of members of the larger classes have resulted in conflicting conclusions concerning interspecies relationships among these species, particularly for the Phytophthora genus. The onset of next-generation sequencing techniques now means that a wealth of oomycete genomic data is available. For the first time, we have used genome-scale phylogenetic methods to resolve oomycete phylogenetic relationships. We used supertree methods to generate single-gene and multigene species phylogenies. Overall, our supertree analyses utilized phylogenetic data from 8,355 oomycete gene families. We have also complemented our analyses with superalignment phylogenies derived from 131 single-copy ubiquitous gene families. Our results show that a genome-scale approach to oomycete phylogeny resolves oomycete classes and clades. Our analysis represents an important first step in large-scale phylogenomic analysis of the oomycetes.

Evolution of the Sterol Biosynthetic Pathway of Pythium insidiosum and Related Oomycetes Contributes to Antifungal Drug Resistance

Pythiosis is a life-threatening infectious disease caused by the oomycete Pythium insidiosum Direct exposure to Py. insidiosum zoospores can initiate infections of the eye, limb, gastrointestinal tract, or skin/subcutaneous tissue. Treatments for pythiosis have mostly relied on surgery. Antifungal drugs are generally ineffective against Py. insidiosum However, one patient with an invasive Py. insidiosum infection recovered completely following treatment with terbinafine and itraconazole. Additionally, the drug target sterol biosynthetic enzymes have been identified in the oomycete Aphanomyces euteiches It remains an open question whether Py. insidiosum is susceptible to the antifungal drugs and harbors any of the known drug target enzymes. Here, we determined the in vitro susceptibilities of terbinafine and itraconazole against 30 isolates of Py. insidiosum We also analyzed endogenous sterols and searched for genes encoding the sterol biosynthetic enzymes in the genomes of Py. insidiosum and related oomycetes. The susceptibility assay showed that the growth of each of the Py. insidiosum isolates was inhibited by the antifungal agents, but only at difficult-to-achieve concentrations, which explains the clinical resistance of the drugs in the treatment of pythiosis patients. Genome searches of Py. insidiosum and related oomycetes demonstrated that these organisms contained an incomplete set of sterol biosynthetic enzymes. Gas chromatographic mass spectrometry did not detect any sterol end products in Py. insidiosum In conclusion, Py. insidiosum possesses an incomplete sterol biosynthetic pathway. Resistance to antifungal drugs targeting enzymes in the ergosterol biosynthetic pathway in Py. insidiosum was due to modifications or losses of some of the genes encoding the drug target enzymes.

Evolution and diversity of cadherins and catenins

Cadherin genes encode a superfamily of conserved transmembrane proteins that share an adhesive ectodomain composed of tandem cadherin repeats. More than 100 human cadherin superfamily members have been identified, which can be classified into three families: major cadherins, protocadherins and cadherin-related proteins. These superfamily members are involved in diverse fundamental cellular processes including cell-cell adhesion, morphogenesis, cell recognition and signaling. Epithelial cadherin (E-cadherin) is the founding cadherin family member. Its cytoplasmic tail interacts with the armadillo catenins, p120 and β-catenin. Further, α-catenin links the cadherin/armadillo catenin complex to the actin filament network. Even genomes of ancestral metazoan species such as cnidarians and placozoans encode a limited number of distinct cadherins and catenins, emphasizing the conservation and functional importance of these gene families. Moreover, a large expansion of the cadherin and catenin families coincides with the emergence of vertebrates and reflects a major functional diversification in higher metazoans. Here, we revisit and review the functions, phylogenetic classifications and co-evolution of the cadherin and catenin protein families.

PAMPs, PRRs, effectors and R-genes associated with citrus-pathogen interactions

BACKGROUND

Recent application of molecular-based technologies has considerably advanced our understanding of complex processes in plant-pathogen interactions and their key components such as PAMPs, PRRs, effectors and R-genes. To develop novel control strategies for disease prevention in citrus, it is essential to expand and consolidate our knowledge of the molecular interaction of citrus plants with their pathogens.

SCOPE

This review provides an overview of our understanding of citrus plant immunity, focusing on the molecular mechanisms involved in the interactions with viruses, bacteria, fungi, oomycetes and vectors related to the following diseases: tristeza, psorosis, citrus variegated chlorosis, citrus canker, huanglongbing, brown spot, post-bloom, anthracnose, gummosis and citrus root rot.

Identification of effector-like proteins in Trichoderma spp. and role of a hydrophobin in the plant-fungus interaction and mycoparasitism

BACKGROUND

Trichoderma spp. can establish beneficial interactions with plants by promoting plant growth and defense systems, as well as, antagonizing fungal phytopathogens in mycoparasitic interactions. Such interactions depend on signal exchange between both participants and can be mediated by effector proteins that alter the host cell structure and function, allowing the establishment of the relationship. The main purpose of this work was to identify, using computational methods, candidates of effector proteins from T. virens, T. atroviride and T. reesei, validate the expression of some of the genes during a beneficial interaction and mycoparasitism and to define the biological function for one of them.

RESULTS

We defined a catalogue of putative effector proteins from T. virens, T. atroviride and T. reesei. We further validated the expression of 16 genes encoding putative effector proteins from T. virens and T. atroviride during the interaction with the plant Arabidopsis thaliana, and with two anastomosis groups of the phytopathogenic fungus Rhizoctonia solani. We found genes which transcript levels are modified in response to the presence of both plant fungi, as well as genes that respond only to either a plant or a fungal host. Further, we show that overexpression of the gene tvhydii1, a Class II hydrophobin family member, enhances the antagonistic activity of T. virens against R. solani AG2. Further, deletion of tvhydii1 results in reduced colonization of plant roots, while its overexpression increases it.

CONCLUSIONS

Our results show that Trichoderma is able to respond in different ways to the presence of a plant or a fungal host, and it can even distinguish between different strains of fungi of a given species. The putative effector proteins identified here may play roles in preventing perception of the fungus by its hosts, favoring host colonization or protecting it from the host's defense response. Finally, the novel effector protein TVHYDII1 plays a role in plant root colonization by T, virens, and participates in its antagonistic activity against R. solani.

The Comparison of Expressed Candidate Secreted Proteins from Two Arbuscular Mycorrhizal Fungi Unravels Common and Specific Molecular Tools to Invade Different Host Plants

Arbuscular mycorrhizal fungi (AMF), belonging to the fungal phylum Glomeromycota, form mutualistic symbioses with roots of almost 80% of land plants. The release of genomic data from the ubiquitous AMF Rhizophagus irregularis revealed that this species possesses a large set of putative secreted proteins (RiSPs) that could be of major importance for establishing the symbiosis. In the present study, we aimed to identify SPs involved in the establishment of AM symbiosis based on comparative gene expression analyses. We first curated the secretome of the R. irregularis DAOM 197198 strain based on two available genomic assemblies. Then we analyzed the expression patterns of the putative RiSPs obtained from the fungus in symbiotic association with three phylogenetically distant host plants-a monocot, a dicot and a liverwort-in comparison with non-symbiotic stages. We found that 33 out of 84 RiSPs induced in planta were commonly up-regulated in these three hosts. Most of these common RiSPs are small proteins of unknown function that may represent putative host non-specific effector proteins. We further investigated the expressed secretome of Gigaspora rosea, an AM fungal species phylogenetically distant from R. irregularis. G. rosea also presents original symbiotic features, a narrower host spectrum and a restrictive geographic distribution compared to R. irregularis. Interestingly, when analyzing up-regulated G. rosea SPs (GrSPs) in different hosts, a higher ratio of host-specific GrSPs was found compared to RiSPs. Such difference of expression patterns may mirror the restrained host spectrum of G. rosea compared to R. irregularis. Finally, we identified a set of conserved SPs, commonly up-regulated by both fungi in all hosts tested, that could correspond to common keys of AMF to colonize host plants. Our data thus highlight the specificities of two distant AM fungi and help in understanding their conserved and specific strategies to invade different hosts.

Emerging oomycete threats to plants and animals

Oomycetes, or water moulds, are fungal-like organisms phylogenetically related to algae. They cause devastating diseases in both plants and animals. Here, we describe seven oomycete species that are emerging or re-emerging threats to agriculture, horticulture, aquaculture and natural ecosystems. They include the plant pathogens Phytophthora infestans, Phytophthora palmivora, Phytophthora ramorum, Plasmopara obducens, and the animal pathogens Aphanomyces invadans, Saprolegnia parasitica and Halioticida noduliformans For each species, we describe its pathology, importance and impact, discuss why it is an emerging threat and briefly review current research activities.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Diversity, evolution and expression profiles of histone acetyltransferases and deacetylases in oomycetes

BACKGROUND

Oomycetes are a group of fungus-like eukaryotes with diverse microorganisms living in marine, freshwater and terrestrial environments. Many of them are important pathogens of plants and animals, causing severe economic losses. Based on previous study, gene expression in eukaryotic cells is regulated by epigenetic mechanisms such as DNA methylation and histone modification. However, little is known about epigenetic mechanisms of oomycetes.

RESULTS

In this study, we investigated the candidate genes in regulating histone acetylation in oomycetes genomes through bioinformatics approaches and identified a group of diverse histone acetyltransferases (HATs) and histone deacetylases (HDACs), along with three putative novel HATs. Phylogenetic analyses suggested that most of these oomycetes HATs and HDACs derived from distinct evolutionary ancestors. Phylogenetic based analysis revealed the complex and distinct patterns of duplications and losses of HATs and HDACs in oomycetes. Moreover, gene expression analysis unveiled the specific expression patterns of the 33 HATs and 11 HDACs of Phytophthora infestans during the stages of development, infection and stress response.

CONCLUSIONS

In this study, we reveal the structure, diversity and the phylogeny of HATs and HDACs of oomycetes. By analyzing the expression data, we provide an overview of the specific biological stages of these genes involved. Our datasets provide useful inputs to help explore the epigenetic mechanisms and the relationship between genomes and phenotypes of oomycetes.

Mitochondrial genomes and comparative genomics of Aphanomyces astaci and Aphanomyces invadans

The genus Aphanomyces (Saprolegniales, Oomycetes) includes species with a variety of ecologies from saprotrophs to plant and animal parasites. Two important species in this genus are A. astaci, the cause of crayfish plague and its close relative, A. invadans, which causes the epizootic ulcerative syndrome on fish. In this study, we have assembled and annotated the mitochondrial (mt) genomes of A. astaci and A. invadans from the whole genome shotgun sequence reads (PRJNA187372; PRJNA258292, respectively). The assembly was generated from A. astaci Pc-genotype strain APO3 and A. invadans strain NJM9701. The sizes of the mtDNAs were 49,489 bp and 49,061 bp for A. astaci and A. invadans, respectively. The species shared similar genetic content and organization encoding 35 proteins, two ribosomal RNAs, three putative open reading frames and 33 transfer RNAs of 19 amino acids for peptide synthesis. Both species also had a large inverted repeat region (LIR) of approximately 12 kb, the LIR contained large and small ribosomal RNAs and eight protein coding genes. These annotated mt genomes serve as a valuable genetic backbone for further development of diagnostic methods and phylogenetic and migration studies of the animal parasitic species of Aphanomyces.

Resurgence of Less-Studied Smut Fungi as Models of Phytopathogenesis in the Omics Age

The smut fungi form a large, diverse, and nonmonophyletic group of plant pathogens that have long served as both important pests of human agriculture and, also, as fertile organisms of scientific investigation. As modern techniques of molecular genetic analysis became available, many previously studied species that proved refractive to these techniques fell by the wayside and were neglected. Now, as the advent of rapid and affordable next-generation sequencing provides genomic and transcriptomic resources for even these "forgotten" fungi, several species are making a comeback and retaking prominent places in phytopathogenic research. In this review, we highlight several of these smut fungi, with special emphasis on Microbotryum lychnidis-dioicae, an anther smut whose molecular genetic tools have finally begun to catch up with its historical importance in classical genetics and now provide mechanistic insights for ecological studies, evolution of host-pathogen interaction, and investigations of emerging infectious disease.

Characterization of the Verticillium dahliae Exoproteome Involves in Pathogenicity from Cotton-Containing Medium

Verticillium wilt, caused by the Verticillium dahliae phytopathogen, is a devastating disease affecting many economically important crops. Previous studies have shown that the exoproteome of V. dahliae plays a significant role in this pathogenic process, but the components and mechanisms that underlie this remain unclear. In this study, the exoproteome of V. dahliae was induced in a cotton-containing C’zapek-Dox (CCD) medium and quantified using the high-throughput isobaric tag technique for relative and absolute quantification (iTRAQ). Results showed that the abundance of 271 secreted proteins was affected by the CCD medium, of which 172 contain typical signal peptides generally produced by the Golgi/endoplasmic reticulum (ER). These enhanced abundance proteins were predominantly enriched in carbohydrate hydrolases; 126 were classified as carbohydrate-active (CAZymes) and almost all were significantly up-regulated in the CCD medium. Results showed that CAZymes proteins 30 and 22 participate in pectin and cellulose degradation pathways, corresponding with the transcription levels of several genes encoded plant cell wall degradation enzyme activated significantly during cotton infection. In addition, targeted deletion of two pectin lyase genes (VdPL3.1 and VdPL3.3) impaired wilt virulence to cotton. This study demonstrates that the V. dahliae exoproteome plays a crucial role in the development of symptoms of wilting and necrosis, predominantly via the pathogenic mechanisms of plant cell wall degradation as part of host plant infection.

Comparative Genomic Analysis among Four Representative Isolates of Phytophthora sojae Reveals Genes under Evolutionary Selection

Comparative genomic analysis is useful for identifying genes affected by evolutionary selection and for studying adaptive variation in gene functions. In Phytophthora sojae, a model oomycete plant pathogen, the related study is lacking. We compared sequence data among four isolates of P. sojae, which represent its four major genotypes. These isolates exhibited >99.688%, >99.864%, and >98.981% sequence identities at genome, gene, and non-gene regions, respectively. One hundred and fifty-three positive selection and 139 negative selection candidate genes were identified. Between the two categories of genes, the positive selection genes were flanked by larger intergenic regions, poorly annotated in function, and less conserved; they had relatively lower transcription levels but many genes had increased transcripts during infection. Genes coding for predicted secreted proteins, particularly effectors, were overrepresented in positive selection. Several RxLR effector genes were identified as positive selection genes, exhibiting much stronger positive selection levels. In addition, candidate genes with presence/absence polymorphism were analyzed. This study provides a landscape of genomic variation among four representative P. sojae isolates and characterized several evolutionary selection-affected gene candidates. The results suggest a relatively covert two-speed genome evolution pattern in P. sojae and will provide clues for identification of new virulence factors in the oomycete plant pathogens.

Splice Sites Seldom Slide: Intron Evolution in Oomycetes

We examine exon junctions near apparent amino acid insertions and deletions in alignments of orthologous protein-coding genes. In 1,917 ortholog families across nine oomycete genomes, 10–20% of introns are near an alignment gap, indicating at first sight that splice-site displacements are frequent. We designed a robust algorithmic procedure for the delineation of intron-containing homologous regions, and combined it with a parsimony-based reconstruction of intron loss, gain, and splice-site shift events on a phylogeny. The reconstruction implies that 12% of introns underwent an acceptor-site shift, and 10% underwent a donor-site shift. In order to offset gene annotation problems, we amended the procedure with the reannotation of intron boundaries using alignment evidence. The corresponding reconstruction involves much fewer intron gain and splice-site shift events. The frequency of acceptor- and donor-side shifts drops to 4% and 3%, respectively, which are not much different from what one would expect by random codon insertions and deletions. In other words, gaps near exon junctions are mostly artifacts of gene annotation rather than evidence of sliding intron boundaries. Our study underscores the importance of using well-supported gene structure annotations in comparative studies. When transcription evidence is not available, we propose a robust ancestral reconstruction procedure that corrects misannotated intron boundaries using sequence alignments. The results corroborate the view that boundary shifts and complete intron sliding are only accidental in eukaryotic genome evolution and have a negligible impact on protein diversity.

Nonagonal cadherins: A new protein family found within the Stramenopiles

Cadherins, a group of molecules typically associated with planar cell polarity and Wnt signalling, have been little reported outside of the animal kingdom. Here, we identify a new family of cadherins in the Stramenopiles, termed Nonagonal after their 9 transmembrane passes, which contrast to the one or seven passes found in other known cadherin families. Manual curation and experimental validation reveal two subclasses of nonagonal cadherins, depending on the number of uninterrupted extracellular cadherin (EC) modules presented. Firstly, shorter mono-exonic, unimodular, protein models, with 3 to 12 EC domains occur as duplicate paralogs in the saprotrophic Labyrinthulomycetes Aurantiochytrium limanicum and Schizochytrium aggregatum, the gastrointestinal Blastocystis hominis (Blastocystae) and as a single copy gene in the autotrophic Pelagophyte Aureococcus anophagefferens. Larger, single copy, multi-exonal, tri-modular protein models, with up to 72 EC domain in total, are found in the Oomycete genera Albugo, Phytophthora, Pythium and Eurychasma. No homolog was found in the closely related autotrophic Phaeophyceae (brown algae) or Bacillariophyceae (diatoms), nor in several genera of plant and animal pathogenic oomycetes (Aphanomyces, Saprolegnia and Hyaloperonospora). This potential absence was further investigated by synteny analysis of the genome regions flanking the cadherin gene models, which are found to be highly variable. Novel to this new cadherin family is the presence of intercalated laminin and putative carbohydrate binding in tri-modular oomycete cadherins and at the N-terminus of thraustochytrid proteins. As we were unable to detect any homologs of proteins involved in signalling pathways where other cadherin families are involved, we present a conceptual hypothesis on the function of nonagonal cadherin based around the presence of putative carbohydrate binding domains.

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.

Sequencing of the Litchi Downy Blight Pathogen Reveals It Is a Phytophthora Species With Downy Mildew-Like Characteristics

On the basis of its downy mildew-like morphology, the litchi downy blight pathogen was previously named Peronophythora litchii. Recently, however, it was proposed to transfer this pathogen to Phytophthora clade 4. To better characterize this unusual oomycete species and important fruit pathogen, we obtained the genome sequence of Phytophthora litchii and compared it to those from other oomycete species. P. litchii has a small genome with tightly spaced genes. On the basis of a multilocus phylogenetic analysis, the placement of P. litchii in the genus Phytophthora is strongly supported. Effector proteins predicted included 245 RxLR, 30 necrosis-and-ethylene-inducing protein-like, and 14 crinkler proteins. The typical motifs, phylogenies, and activities of these effectors were typical for a Phytophthora species. However, like the genome features of the analyzed downy mildews, P. litchii exhibited a streamlined genome with a relatively small number of genes in both core and species-specific protein families. The low GC content and slight codon preferences of P. litchii sequences were similar to those of the analyzed downy mildews and a subset of Phytophthora species. Taken together, these observations suggest that P. litchii is a Phytophthora pathogen that is in the process of acquiring downy mildew-like genomic and morphological features. Thus P. litchii may provide a novel model for investigating morphological development and genomic adaptation in oomycete pathogens.

Small homologous blocks in phytophthora genomes do not point to an ancient whole-genome duplication

Genomes of the plant-pathogenic genus Phytophthora are characterized by small duplicated blocks consisting of two consecutive genes (2HOM blocks) and by an elevated abundance of similarly aged gene duplicates. Both properties, in particular the presence of 2HOM blocks, have been attributed to a whole-genome duplication (WGD) at the last common ancestor of Phytophthora. However, large intraspecies synteny—compelling evidence for a WGD—has not been detected. Here, we revisited the WGD hypothesis by deducing the age of 2HOM blocks. Two independent timing methods reveal that the majority of 2HOM blocks arose after divergence of the Phytophthora lineages. In addition, a large proportion of the 2HOM block copies colocalize on the same scaffold. Therefore, the presence of 2HOM blocks does not support a WGD at the last common ancestor of Phytophthora. Thus, genome evolution of Phytophthora is likely driven by alternative mechanisms, such as bursts of transposon activity.

SCR96, a small cysteine-rich secretory protein of Phytophthora cactorum, can trigger cell death in the Solanaceae and is important for pathogenicity and oxidative stress tolerance

Peptides and small molecules produced by both the plant pathogen Phytophthora and host plants in the apoplastic space mediate the relationship between the interplaying organisms. Various Phytophthora apoplastic effectors, including small cysteine-rich (SCR) secretory proteins, have been identified, but their roles during interaction remain to be determined. Here, we identified an SCR effector encoded by scr96, one of three novel genes encoding SCR proteins in P. cactorum with similarity to the P. cactorum phytotoxic protein PcF. Together with the other two genes, scr96 was transcriptionally induced throughout the developmental and infection stages of the pathogen. These genes triggered plant cell death (PCD) in the Solanaceae, including Nicotiana benthamiana and tomato. The scr96 gene did not show single nucleotide polymorphisms in a collection of P. cactorum isolates from different countries and host plants, suggesting that its role is essential and non-redundant during infection. Homologues of SCR96 were identified only in oomycetes, but not in fungi and other organisms. A stable protoplast transformation protocol was adapted for P. cactorum using green fluorescent protein as a marker. The silencing of scr96 in P. cactorum caused gene-silenced transformants to lose their pathogenicity on host plants and these transformants were significantly more sensitive to oxidative stress. Transient expression of scr96 partially recovered the virulence of gene-silenced transformants on plants. Overall, our results indicate that the P. cactorum scr96 gene encodes an important virulence factor that not only causes PCD in host plants, but is also important for pathogenicity and oxidative stress tolerance.

Draft Genome Sequence of Biocontrol Agent Pythium oligandrum Strain Po37, an Oomycota

The oomycotaPythium oligandrumPo37 is used as a biocontrol agent of plant diseases. Here, we present the first draft of theP. oligandrumPo37 genome sequence, which comprises 725 scaffolds with a total length of 35.9 Mb and 11,695 predicted protein-coding genes.

CRN13 candidate effectors from plant and animal eukaryotic pathogens are DNA-binding proteins which trigger host DNA damage response

To successfully colonize their host, pathogens produce effectors that can interfere with host cellular processes. Here we investigated the function of CRN13 candidate effectors produced by plant pathogenic oomycetes and detected in the genome of the amphibian pathogenic chytrid fungus Batrachochytrium dendrobatidis (BdCRN13). When expressed in Nicotiana, AeCRN13, from the legume root pathogen Aphanomyces euteiches, increases the susceptibility of the leaves to the oomycete Phytophthora capsici. When transiently expressed in amphibians or plant cells, AeCRN13 and BdCRN13 localize to the cell nuclei, triggering aberrant cell development and eventually causing cell death. Using Förster resonance energy transfer experiments in plant cells, we showed that both CRN13s interact with nuclear DNA and trigger plant DNA damage response (DDR). Mutating key amino acid residues in a predicted HNH-like endonuclease motif abolished the interaction of AeCRN13 with DNA, the induction of DDR and the enhancement of Nicotiana susceptibility to P. capsici. Finally, H2AX phosphorylation, a marker of DNA damage, and enhanced expression of genes involved in the DDR were observed in A. euteiches-infected Medicago truncatula roots. These results show that CRN13 from plant and animal eukaryotic pathogens promotes host susceptibility by targeting nuclear DNA and inducing DDR.

Genome analysis of the sugar beet pathogen Rhizoctonia solani AG2-2IIIB revealed high numbers in secreted proteins and cell wall degrading enzymes

Background

Sugar beet (Beta vulgaris) is a crop cultivated for its high content in sugar, but it is vulnerable to many soil-borne pathogens. One of them is the basidiomycete Rhizoctonia solani. This fungal species has a compatibility system regulating hyphal fusions (anastomosis). Consequently, R. solani species are categorized in anastomosis groups (AGs). AG2-2IIIB isolates are most aggressive on sugar beet. In the present study, we report on the draft genome of R. solani AG2-2IIIB using the Illumina technology. Genome analysis, interpretation and comparative genomics of five sequenced R. solani isolates were carried out.

Results

The draft genome of R. solani AG2-2IIIB has an estimated size of 56.02 Mb. In addition, two normalized EST libraries were sequenced. In total 20,790 of 21,980 AG2-2IIIB isotigs (transcript isoforms) were mapped on the genome with more than 95 % sequence identity. The genome of R. solani AG2-2IIIB was predicted to harbor 11,897 genes and 4908 were found to be isolate-specific. R. solani AG2-2IIIB was predicted to contain 1142 putatively secreted proteins and 473 of them were found to be unique for this isolate. The R. solani AG2-2IIIB genome encodes a high number of carbohydrate active enzymes. The highest numbers were observed for the polysaccharide lyases family 1 (PL-1), glycoside hydrolase family 43 (GH-43) and carbohydrate estarase family 12 (CE-12). Transcription analysis of selected genes representing different enzyme clades revealed a mixed pattern of up- and down-regulation six days after infection on sugar beets featuring variable levels of resistance compared to mycelia of the fungus grown in vitro.

Conclusions

The established R. solani AG2-2IIIB genome and EST sequences provide important information on the gene content, gene structure and transcriptional activity for this sugar beet pathogen. The enriched genomic platform provides an important platform to enhance our understanding of R. solani biology.

Electronic supplementary material

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

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.

Whole genome sequence of the emerging oomycete pathogen Pythium insidiosum strain CDC-B5653 isolated from an infected human in the USA

Pythium insidiosum ATCC 200269 strain CDC-B5653, an isolate from necrotizing lesions on the mouth and eye of a 2-year-old boy in Memphis, Tennessee, USA, was sequenced using a combination of Illumina MiSeq (300 bp paired-end, 14 millions reads) and PacBio (10  Kb fragment library, 356,001 reads). The sequencing data were assembled using SPAdes version 3.1.0, yielding a total genome size of 45.6 Mb contained in 8992 contigs, N50 of 13 Kb, 57% G + C content, and 17,867 putative protein-coding genes. This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession JRHR00000000.

Genomes and virulence difference between two physiological races of Phytophthora nicotianae

BACKGROUND

Black shank is a severe plant disease caused by the soil-borne pathogen Phytophthora nicotianae. Two physiological races of P. nicotianae, races 0 and 1, are predominantly observed in cultivated tobacco fields around the world. Race 0 has been reported to be more aggressive, having a shorter incubation period, and causing worse root rot symptoms, while race 1 causes more severe necrosis. The molecular mechanisms underlying the difference in virulence between race 0 and 1 remain elusive.

FINDINGS

We assembled and annotated the genomes of P. nicotianae races 0 and 1, which were obtained by a combination of PacBio single-molecular real-time sequencing and second-generation sequencing (both HiSeq and MiSeq platforms). Gene family analysis revealed a highly expanded ATP-binding cassette transporter gene family in P. nicotianae. Specifically, more RxLR effector genes were found in the genome of race 0 than in that of race 1. In addition, RxLR effector genes were found to be mainly distributed in gene-sparse, repeat-rich regions of the P. nicotianae genome.

CONCLUSIONS

These results provide not only high quality reference genomes of P. nicotianae, but also insights into the infection mechanisms of P. nicotianae and its co-evolution with the host plant. They also reveal insights into the difference in virulence between the two physiological races.

Diverse Evolutionary Trajectories for Small RNA Biogenesis Genes in the Oomycete Genus Phytophthora

Gene regulation by small RNA pathways is ubiquitous among eukaryotes, but little is known about small RNA pathways in the Stramenopile kingdom. Phytophthora, a genus of filamentous oomycetes, contains many devastating plant pathogens, causing multibillion-dollar damage to crops, ornamental plants, and natural environments. The genomes of several oomycetes including Phytophthora species such as the soybean pathogen P. sojae, have been sequenced, allowing evolutionary analysis of small RNA-processing enzymes. This study examined the evolutionary origins of the oomycete small RNA-related genes Dicer-like (DCL), and RNA-dependent RNA polymerase (RDR) through broad phylogenetic analyses of the key domains. Two Dicer gene homologs, DCL1 and DCL2, and one RDR homolog were cloned and analyzed from P. sojae. Gene expression analysis revealed only minor changes in transcript levels among different life stages. Oomycete DCL1 homologs clustered with animal and plant Dicer homologs in evolutionary trees, whereas oomycete DCL2 homologs clustered basally to the tree along with Drosha homologs. Phylogenetic analysis of the RDR homologs confirmed a previous study that suggested the last common eukaryote ancestor possessed three RDR homologs, which were selectively retained or lost in later lineages. Our analysis clarifies the position of some Unikont and Chromalveolate RDR lineages within the tree, including oomycete homologs. Finally, we analyzed alterations in the domain structure of oomycete Dicer and RDR homologs, specifically focusing on the proposed domain transfer of the DEAD-box helicase domain from Dicer to RDR. Implications of the oomycete domain structure are discussed, and possible roles of the two oomycete Dicer homologs are proposed.

Overexpression of a Phytophthora Cytoplasmic CRN Effector Confers Resistance to Disease, Salinity and Drought in Nicotiana benthamiana

The Crinkler (CRN) effector family is produced by oomycete pathogens and may manipulate host physiological and biochemical events inside host cells. Here, PsCRN161 was identified from Phytophthora sojae based on its broad and strong cell death suppression activities. The effector protein contains two predicted nuclear localization signals and localized to nuclei of plant cells, indicating that it may target plant nuclei to modify host cell physiology and function. The chimeric gene GFP:PsCRN161 driven by the Cauliflower mosaic virus (CaMV) 35S promoter was introduced into Nicotiana benthamiana. The four independent PsCRN161-transgenic lines exhibited increased resistance to two oomycete pathogens (P. parasitica and P. capsici) and showed enhanced tolerance to salinity and drought stresses. Digital gene expression profiling analysis showed that defense-related genes, including ABC transporters, Cyt P450 and receptor-like kinases (RLKs), were significantly up-regulated in PsCRN161-transgenic plants compared with GFP (green fluorescent protein) lines, implying that PsCRN161 expression may protect plants from biotic and abiotic stresses by up-regulation of many defense-related genes. The results reveal previously unknown functions of the oomycete effectors, suggesting that the pathogen effectors could be directly used as functional genes for plant molecular breeding for enhancement of tolerance to biotic and abiotic stresses.

Effector discovery in the fungal wheat pathogen Zymoseptoria tritici

Fungal plant pathogens, such as Zymoseptoria tritici (formerly known as Mycosphaerella graminicola), secrete repertoires of effectors to facilitate infection or trigger host defence mechanisms. The discovery and functional characterization of effectors provides valuable knowledge that can contribute to the design of new and effective disease management strategies. Here, we combined bioinformatics approaches with expression profiling during pathogenesis to identify candidate effectors of Z. tritici. In addition, a genetic approach was conducted to map quantitative trait loci (QTLs) carrying putative effectors, enabling the validation of both complementary strategies for effector discovery. In planta expression profiling revealed that candidate effectors were up-regulated in successive waves corresponding to consecutive stages of pathogenesis, contrary to candidates identified by QTL mapping that were, overall, expressed at low levels. Functional analyses of two top candidate effectors (SSP15 and SSP18) showed their dispensability for Z. tritici pathogenesis. These analyses reveal that generally adopted criteria, such as protein size, cysteine residues and expression during pathogenesis, may preclude an unbiased effector discovery. Indeed, genetic mapping of genomic regions involved in specificity render alternative effector candidates that do not match the aforementioned criteria, but should nevertheless be considered as promising new leads for effectors that are crucial for the Z. tritici-wheat pathosystem.

Genome Sequence and Architecture of the Tobacco Downy Mildew Pathogen Peronospora tabacina

Peronospora tabacina is an obligate biotrophic oomycete that causes blue mold or downy mildew on tobacco (Nicotiana tabacum). It is an economically important disease occurring frequently in tobacco-growing regions worldwide. We sequenced and characterized the genomes of two P. tabacina isolates and mined them for pathogenicity-related proteins and effector-encoding genes. De novo assembly of the genomes using Illumina reads resulted in 4,016 (63.1 Mb, N50 = 79 kb) and 3,245 (55.3 Mb, N50 = 61 kb) scaffolds for isolates 968-J2 and 968-S26, respectively, with an estimated genome size of 68 Mb. The mitochondrial genome has a similar size (approximately 43 kb) and structure to those of other oomycetes, plus several minor unique features. Repetitive elements, primarily retrotransposons, make up approximately 24% of the nuclear genome. Approximately 18,000 protein-coding gene models were predicted. Mining the secretome revealed approximately 120 candidate RxLR, six CRN (candidate effectors that elicit crinkling and necrosis), and 61 WY domain-containing proteins. Candidate RxLR effectors were shown to be predominantly undergoing diversifying selection, with approximately 57% located in variable gene-sparse regions of the genome. Aligning the P. tabacina genome to Hyaloperonospora arabidopsidis and Phytophthora spp. revealed a high level of synteny. Blocks of synteny show gene inversions and instances of expansion in intergenic regions. Extensive rearrangements of the gene-rich genomic regions do not appear to have occurred during the evolution of these highly variable pathogens. These assemblies provide the basis for studies of virulence in this and other downy mildew pathogens.

PhytoPath: an integrative resource for plant pathogen genomics

PhytoPath (www.phytopathdb.org) is a resource for genomic and phenotypic data from plant pathogen species, that integrates phenotypic data for genes from PHI-base, an expertly curated catalog of genes with experimentally verified pathogenicity, with the Ensembl tools for data visualization and analysis. The resource is focused on fungi, protists (oomycetes) and bacterial plant pathogens that have genomes that have been sequenced and annotated. Genes with associated PHI-base data can be easily identified across all plant pathogen species using a BioMart-based query tool and visualized in their genomic context on the Ensembl genome browser. The PhytoPath resource contains data for 135 genomic sequences from 87 plant pathogen species, and 1364 genes curated for their role in pathogenicity and as targets for chemical intervention. Support for community annotation of gene models is provided using the WebApollo online gene editor, and we are working with interested communities to improve reference annotation for selected species.

Genome analyses of the sunflower pathogen Plasmopara halstedii provide insights into effector evolution in downy mildews and Phytophthora

BACKGROUND

Downy mildews are the most speciose group of oomycetes and affect crops of great economic importance. So far, there is only a single deeply-sequenced downy mildew genome available, from Hyaloperonospora arabidopsidis. Further genomic resources for downy mildews are required to study their evolution, including pathogenicity effector proteins, such as RxLR effectors. Plasmopara halstedii is a devastating pathogen of sunflower and a potential pathosystem model to study downy mildews, as several Avr-genes and R-genes have been predicted and unlike Arabidopsis downy mildew, large quantities of almost contamination-free material can be obtained easily.

RESULTS

Here a high-quality draft genome of Plasmopara halstedii is reported and analysed with respect to various aspects, including genome organisation, secondary metabolism, effector proteins and comparative genomics with other sequenced oomycetes. Interestingly, the present analyses revealed further variation of the RxLR motif, suggesting an important role of the conservation of the dEER-motif. Orthology analyses revealed the conservation of 28 RxLR-like core effectors among Phytophthora species. Only six putative RxLR-like effectors were shared by the two sequenced downy mildews, highlighting the fast and largely independent evolution of two of the three major downy mildew lineages. This is seemingly supported by phylogenomic results, in which downy mildews did not appear to be monophyletic.

CONCLUSIONS

The genome resource will be useful for developing markers for monitoring the pathogen population and might provide the basis for new approaches to fight Phytophthora and downy mildew pathogens by targeting core pathogenicity effectors.

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.

Multiple Avirulence Loci and Allele-Specific Effector Recognition Control the Pm3 Race-Specific Resistance of Wheat to Powdery Mildew

In cereals, several mildew resistance genes occur as large allelic series; for example, in wheat (Triticum aestivum and Triticum turgidum), 17 functional Pm3 alleles confer agronomically important race-specific resistance to powdery mildew (Blumeria graminis). The molecular basis of race specificity has been characterized in wheat, but little is known about the corresponding avirulence genes in powdery mildew. Here, we dissected the genetics of avirulence for six Pm3 alleles and found that three major Avr loci affect avirulence, with a common locus_1 involved in all AvrPm3-Pm3 interactions. We cloned the effector gene AvrPm3(a2/f2) from locus_2, which is recognized by the Pm3a and Pm3f alleles. Induction of a Pm3 allele-dependent hypersensitive response in transient assays in Nicotiana benthamiana and in wheat demonstrated specificity. Gene expression analysis of Bcg1 (encoded by locus_1) and AvrPm3 (a2/f2) revealed significant differences between isolates, indicating that in addition to protein polymorphisms, expression levels play a role in avirulence. We propose a model for race specificity involving three components: an allele-specific avirulence effector, a resistance gene allele, and a pathogen-encoded suppressor of avirulence. Thus, whereas a genetically simple allelic series controls specificity in the plant host, recognition on the pathogen side is more complex, allowing flexible evolutionary responses and adaptation to resistance genes.

Comparative Genomics Including the Early-Diverging Smut Fungus Ceraceosorus bombacis Reveals Signatures of Parallel Evolution within Plant and Animal Pathogens of Fungi and Oomycetes

Ceraceosorus bombacis is an early-diverging lineage of smut fungi and a pathogen of cotton trees (Bombax ceiba). To study the evolutionary genomics of smut fungi in comparison with other fungal and oomycete pathogens, the genome of C. bombacis was sequenced and comparative genomic analyses were performed. The genome of 26.09 Mb encodes for 8,024 proteins, of which 576 are putative-secreted effector proteins (PSEPs). Orthology analysis revealed 30 ortholog PSEPs among six Ustilaginomycotina genomes, the largest groups of which are lytic enzymes, such as aspartic peptidase and glycoside hydrolase. Positive selection analyses revealed the highest percentage of positively selected PSEPs in C. bombacis compared with other Ustilaginomycotina genomes. Metabolic pathway analyses revealed the absence of genes encoding for nitrite and nitrate reductase in the genome of the human skin pathogen Malassezia globosa, but these enzymes are present in the sequenced plant pathogens in smut fungi. Interestingly, these genes are also absent in cultivable oomycete animal pathogens, while nitrate reductase has been lost in cultivable oomycete plant pathogens. Similar patterns were also observed for obligate biotrophic and hemi-biotrophic fungal and oomycete pathogens. Furthermore, it was found that both fungal and oomycete animal pathogen genomes are lacking cutinases and pectinesterases. Overall, these findings highlight the parallel evolution of certain genomic traits, revealing potential common evolutionary trajectories among fungal and oomycete pathogens, shaping the pathogen genomes according to their lifestyle.

Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity

How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.

RNA-Seq Analysis of the Expression of Genes Encoding Cell Wall Degrading Enzymes during Infection of Lupin (Lupinus angustifolius) by Phytophthora parasitica

RNA-Seq analysis has shown that over 60% (12,962) of the predicted transcripts in the Phytophthora parasitica genome are expressed during the first 60 h of lupin root infection. The infection transcriptomes included 278 of the 431 genes encoding P. parasitica cell wall degrading enzymes. The transcriptome data provide strong evidence of global transcriptional cascades of genes whose encoded proteins target the main categories of plant cell wall components. A major cohort of pectinases is predominantly expressed early but as infection progresses, the transcriptome becomes increasingly dominated by transcripts encoding cellulases, hemicellulases, β-1,3-glucanases and glycoproteins. The most highly expressed P. parasitica carbohydrate active enzyme gene contains two CBM1 cellulose binding modules and no catalytic domains. The top 200 differentially expressed genes include β-1,4-glucosidases, β-1,4-glucanases, β-1,4-galactanases, a β-1,3-glucanase, an α-1,4-polygalacturonase, a pectin deacetylase and a pectin methylesterase. Detailed analysis of gene expression profiles provides clues as to the order in which linkages within the complex carbohydrates may come under attack. The gene expression profiles suggest that (i) demethylation of pectic homogalacturonan occurs before its deacetylation; (ii) cleavage of the backbone of pectic rhamnogalacturonan I precedes digestion of its side chains; (iii) early attack on cellulose microfibrils by non-catalytic cellulose-binding proteins and enzymes with auxiliary activities may facilitate subsequent attack by glycosyl hydrolases and enzymes containing CBM1 cellulose-binding modules; (iv) terminal hemicellulose backbone residues are targeted after extensive internal backbone cleavage has occurred; and (v) the carbohydrate chains on glycoproteins are degraded late in infection. A notable feature of the P. parasitica infection transcriptome is the high level of transcription of genes encoding enzymes that degrade β-1,3-glucanases during middle and late stages of infection. The results suggest that high levels of β-1,3-glucanases may effectively degrade callose as it is produced by the plant during the defence response.