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.

Black Sigatoka or black leaf streak disease, caused by the Dothideomycete fungus Pseudocercospora fijiensis (previously: Mycosphaerella fijiensis), is the most significant foliar disease of banana worldwide. Due to the lack of effective host resistance, management of this disease requires frequent fungicide applications, which greatly increase the economic and environmental costs to produce banana. Weekly applications in most banana plantations lead to rapid evolution of fungicide-resistant strains within populations causing disease-control failures throughout the world. Given its extremely high economic importance, two strains of P. fijiensis were sequenced and assembled with the aid of a new genetic linkage map. The 74-Mb genome of P. fijiensis is massively expanded by LTR retrotransposons, making it the largest genome within the Dothideomycetes. Melting-curve assays suggest that the genomes of two closely related members of the Sigatoka disease complex, P. eumusae and P. musae, also are expanded. Electrophoretic karyotyping and analyses of molecular markers in P. fijiensis field populations showed chromosome-length polymorphisms and high genetic diversity. Genetic differentiation was also detected using neutral markers, suggesting strong selection with limited gene flow at the studied geographic scale. Frequencies of fungicide resistance in fungicide-treated plantations were much higher than those in untreated wild-type P. fijiensis populations. A homologue of the Cladosporium fulvum Avr4 effector, PfAvr4, was identified in the P. fijiensis genome. Infiltration of the purified PfAVR4 protein into leaves of the resistant banana variety Calcutta 4 resulted in a hypersensitive-like response. This result suggests that Calcutta 4 could carry an unknown resistance gene recognizing PfAVR4. Besides adding to our understanding of the overall Dothideomycete genome structures, the P. fijiensis genome will aid in developing fungicide treatment schedules to combat this pathogen and in improving the efficiency of banana breeding programs.

Fusarium oxysporum (Fo) belongs to a group of soil-borne hyphomycetes that are taxonomically collated in the Fusarium oxysporum Species Complex (FOSC). Hitherto, those infecting bananas were placed in the forma specialis cubense (Foc). Recently, however, these genetically different Foc lineages were recognized as new Fusarium spp. placed in the Fusarium of Banana Complex (FOBC). A member of this complex F. odoratissimum II-5 that uniquely comprises the so-called Tropical Race 4 (TR4), is a major problem sweeping through production zones of Cavendish banana in several regions of the world. Because of this, there is an urgent need for a phenotyping method that allows the screening for resistance to TR4 of large numbers of banana genotypes. Most Fusarium species produce three types of spores: macroconidia, microconidia and the persistent chlamydospores that can contaminate soils for many years. Inoculum production has been an important bottleneck for efficient phenotyping due to the low or variable number of conidia and the elaborate laboratory procedures requiring specific infrastructure. Here, we report a rapid, simple and high-yielding spore production method for nine F. oxysporum formae speciales as well as the biocontrol species Fo47 and Fo618-12. For Fusarium spp. causing Fusarium wilt or Panama disease of banana, we used the protocol for four species comprising the recognized physiological races, including Tropical Race 4 (TR4). We subsequently tested the produced inoculum in comparative inoculation trials on banana plants to evaluate their efficiency. All assays resulted in typical symptoms within 10 weeks; significant differences in final disease ratings were observed, depending on inoculum concentration. Pouring inoculum directly onto banana plants showed the most consistent and reproducible results, as expressed in external wilting, internal discoloration and determined by real-time PCR assays on entire rhizomes. Moreover, this method allows the inoculation of 250 plants per hour by one individual thereby facilitating the phenotyping of large mutant and breeding populations.

Cross-kingdom small RNA (sRNA) silencing has recently emerged as a mechanism facilitating fungal colonization and disease development. Here we characterized RNAi pathways in Zymoseptoria tritici, a major fungal pathogen of wheat, and assessed their contribution to pathogenesis. Computational analysis of fungal sRNA and host mRNA sequencing datasets was used to define the global sRNA populations in Z. tritici and predict their mRNA targets in wheat. 389 in planta-induced sRNA loci were identified. sRNAs generated from some of these loci were predicted to target wheat mRNAs including those potentially involved in pathogen defense. However, molecular approaches failed to validate targeting of selected wheat mRNAs by fungal sRNAs. Mutant strains of Z. tritici carrying deletions of genes encoding key components of RNAi such as Dicer-like (DCL) and Argonaute (AGO) proteins were generated, and virulence bioassays suggested that these are dispensable for full infection of wheat. Nonetheless, our results did suggest the existence of non-canonical DCL-independent pathway(s) for sRNA biogenesis in Z. tritici. dsRNA targeting essential fungal genes applied in vitro or generated from an RNA virus vector in planta in a procedure known as HIGS (Host-Induced Gene Silencing) was ineffective in preventing Z. tritici growth or disease. We also demonstrated that Z. tritici is incapable of dsRNA uptake. Collectively, our data suggest that RNAi approaches for gene function analyses in this fungal species and potentially also as a control measure may not be as effective as has been demonstrated for some other plant pathogenic fungi.

Oomycetes in the class Saprolegniomycetidae of the Eukaryotic kingdom Stramenopila have evolved as severe pathogens of amphibians, crustaceans, fish and insects, resulting in major losses in aquaculture and damage to aquatic ecosystems. We have sequenced the 63 Mb genome of the fresh water fish pathogen, Saprolegnia parasitica. Approximately 1/3 of the assembled genome exhibits loss of heterozygosity, indicating an efficient mechanism for revealing new variation. Comparison of S. parasitica with plant pathogenic oomycetes suggests that during evolution the host cellular environment has driven distinct patterns of gene expansion and loss in the genomes of plant and animal pathogens. S. parasitica possesses one of the largest repertoires of proteases (270) among eukaryotes that are deployed in waves at different points during infection as determined from RNA-Seq data. In contrast, despite being capable of living saprotrophically, parasitism has led to loss of inorganic nitrogen and sulfur assimilation pathways, strikingly similar to losses in obligate plant pathogenic oomycetes and fungi. The large gene families that are hallmarks of plant pathogenic oomycetes such as Phytophthora appear to be lacking in S. parasitica, including those encoding RXLR effectors, Crinkler's, and Necrosis Inducing-Like Proteins (NLP). S. parasitica also has a very large kinome of 543 kinases, 10% of which is induced upon infection. Moreover, S. parasitica encodes several genes typical of animals or animal-pathogens and lacking from other oomycetes, including disintegrins and galactose-binding lectins, whose expression and evolutionary origins implicate horizontal gene transfer in the evolution of animal pathogenesis in S. parasitica.

Sensing external signals and transducing these into intracellular responses requires a molecular signaling system that is crucial for every living organism. Two important eukaryotic signal transduction pathways that are often interlinked are G-protein signaling and phospholipid signaling. Heterotrimeric G-protein subunits activated by G-protein-coupled receptors (GPCRs) are typical stimulators of phospholipid signaling enzymes such as phosphatidylinositol phosphate kinases (PIPKs) or phospholipase C (PLC). However, a direct connection between the two pathways likely exists in oomycetes and slime molds, as they possess a unique class of GPCRs that have a PIPK as an accessory domain. In principle, these so-called GPCR-PIPKs have the capacity of perceiving an external signal (via the GPCR domain) that, via PIPK, directly activates downstream phospholipid signaling. Here we reveal the sporadic occurrence of GPCR-PIPKs in all eukaryotic supergroups, except for plants. Notably, all species having GPCR-PIPKs are unicellular microorganisms that favor aquatic environments. Phylogenetic analysis revealed that GPCR-PIPKs are likely ancestral to eukaryotes and significantly expanded in the last common ancestor of oomycetes. In addition to GPCR-PIPKs, we identified five hitherto-unknown classes of GPCRs with accessory domains, four of which are universal players in signal transduction. Similarly to GPCR-PIPKs, this enables a direct coupling between extracellular sensing and downstream signaling. Overall, our findings point to an ancestral signaling system in eukaryotes where GPCR-mediated sensing is directly linked to downstream responses.IMPORTANCE G-protein-coupled receptors (GPCRs) are central sensors that activate eukaryotic signaling and are the primary targets of human drugs. In this report, we provide evidence for the widespread though limited presence of a novel class of GPCRs in a variety of unicellular eukaryotes. These include free-living organisms and organisms that are pathogenic for plants, animals, and humans. The novel GPCRs have a C-terminal phospholipid kinase domain, pointing to a direct link between sensing external signals via GPCRs and downstream intracellular phospholipid signaling. Genes encoding these receptors were likely present in the last common eukaryotic ancestor and were lost during the evolution of higher eukaryotes. We further describe five other types of GPCRs with a catalytic accessory domain, the so-called GPCR-bigrams, four of which may potentially have a role in signaling. These findings shed new light onto signal transduction in microorganisms and provide evidence for alternative eukaryotic signaling pathways.

Banana is the most popular and most exported fruit and also a major food crop for millions of people around the world. Despite its importance and the presence of serious disease threats, research into this crop is limited. One of those is Panama disease or Fusarium wilt. In the previous century Fusarium wilt wiped out the "Gros Michel" based banana industry in Central America. The epidemic was eventually quenched by planting "Cavendish" bananas. However, 50 years ago the disease recurred, but now on "Cavendish" bananas. Since then the disease has spread across South-East Asia, to the Middle-East and the Indian subcontinent and leaped into Africa. Here, we report the presence of Fusarium oxysporum f.sp. cubense Tropical Race 4 (Foc TR4) in "Cavendish" plantations in Laos, Myanmar, and Vietnam. A combination of classical morphology, DNA sequencing, and phenotyping assays revealed a very close relationship between the Foc TR4 strains in the entire Greater Mekong Subregion (GMS), which is increasingly prone to intensive banana production. Analyses of single-nucleotide polymorphisms enabled us to initiate a phylogeography of Foc TR4 across three geographical areas-GMS, Indian subcontinent, and the Middle East revealing three distinct Foc TR4 sub-lineages. Collectively, our data place these new incursions in a broader agroecological context and underscore the need for awareness campaigns and the implementation of validated quarantine measures to prevent further international dissemination of Foc TR4.

Phytophthora species are oomycete plant pathogens with such major social and economic impact that genome sequences have been determined for Phytophthora infestans, P. sojae and P. ramorum. Pepsin-like aspartic proteinases (APs) are produced in a wide variety of species (from bacteria to humans) and contain conserved motifs and landmark residues. APs fulfil critical roles in infectious organisms and their host cells. Annotation of Phytophthora APs would provide invaluable information for studies into their roles in the physiology of Phytophthora species and interactions with their hosts.

Pseudocercospora fijiensis is the causal agent of the black leaf streak disease (BLSD) of banana. Bananas are important global export commodities and a major staple food. Their susceptibility to BLSD pushes disease management towards excessive fungicide use, largely relying on multisite inhibitors and sterol demethylation inhibitors (DMIs). These fungicides are ubiquitous in plant disease control, targeting the CYP51 enzyme. We examined sensitivity to DMIs in P. fijiensis field isolates collected from various major banana production zones in Colombia, Costa Rica, Dominican Republic, Ecuador, the Philippines, Guadalupe, Martinique and Cameroon and determined the underlying genetic reasons for the observed phenotypes.

Succinate dehydrogenase inhibitor (SDHI) fungicides are widely used for the control of a broad range of fungal diseases. This has been the most rapidly expanding fungicide group in terms of new molecules discovered and introduced for agricultural use over the past fifteen years. A particular pattern of differential sensitivity (resistance) to the stretched heterocycle amide SDHIs (SHA-SDHIs), a subclass of chemically-related SDHIs, was observed in naïve Zymoseptoria tritici populations not previously exposed to these chemicals. Subclass-specific resistance was confirmed at the enzyme level but did not correlate with the genotypes of the succinate dehydrogenase (SDH) encoding genes. Mapping and characterization of the molecular mechanisms responsible for standing SHA-SDHI resistance in natural field isolates identified a gene paralog of SDHC, termed ZtSDHC3, which encodes for an alternative C subunit of succinate dehydrogenase, named alt-SDHC. Using reverse genetics, we showed that alt-SDHC associates with the three other SDH subunits, leading to a fully functional enzyme and that a unique Qp-site residue within the alt-SDHC protein confers SHA-SDHI resistance. Enzymatic assays, computational modelling and docking simulations for the two SQR enzymes (altC-SQR, WT_SQR) enabled us to describe enzyme-inhibitor interactions at an atomistic level and to propose rational explanations for differential potency and resistance across SHA-SDHIs. European Z. tritici populations displayed a presence (20-30%) / absence polymorphism of ZtSDHC3, as well as differences in ZtSDHC3 expression levels and splicing efficiency. These polymorphisms have a strong impact on SHA-SDHI resistance phenotypes. Characterization of the ZtSDHC3 promoter in European Z. tritici populations suggests that transposon insertions are associated with the strongest resistance phenotypes. These results establish that a dispensable paralogous gene determines SHA-SDHIs fungicide resistance in natural populations of Z. tritici. This study paves the way to an increased awareness of the role of fungicidal target paralogs in resistance to fungicides and demonstrates the paramount importance of population genomics in fungicide discovery.

Crop protection anilinopyrimidine (AP) fungicides were introduced more than 20 years ago for the control of a range of diseases caused by ascomycete plant pathogens, and in particular for the control of gray mold caused by Botrytis cinerea. Although early mode of action studies suggested an inhibition of methionine biosynthesis, the molecular target of this class of fungicides was never fully clarified. Despite AP-specific resistance having been described in B. cinerea field isolates and in multiple other targeted species, the underlying resistance mechanisms were unknown. It was therefore expected that the genetic characterization of resistance mechanisms would permit the identification of the molecular target of these fungicides. In order to explore the widest range of possible resistance mechanisms, AP-resistant B. cinerea UV laboratory mutants were generated and the mutations conferring resistance were determined by combining whole-genome sequencing and reverse genetics. Genetic mapping from a cross between a resistant field isolate and a sensitive reference isolate was used in parallel and led to the identification of an additional molecular determinant not found from the characterized UV mutant collection. Together, these two approaches enabled the characterization of an unrivaled diversity of resistance mechanisms. In total, we report the elucidation of resistance-conferring mutations within nine individual genes, two of which are responsible for almost all instances of AP resistance in the field. All identified resistance-conferring genes encode proteins that are involved in mitochondrial processes, suggesting that APs primarily target the mitochondria. The functions of these genes and their possible interactions are discussed in the context of the potential mode of action for this important class of fungicides.

The Dothideomycete Pseudocercospora fijiensis, previously Mycosphaerella fijiensis, is the causal agent of black Sigatoka, one of the most destructive diseases of bananas and plantains. Disease management depends on fungicide applications, with a major contribution from sterol demethylation-inhibitors (DMIs). The continued use of DMIs places considerable selection pressure on natural P. fijiensis populations, enabling the selection of novel genotypes with reduced sensitivity. The hitherto explanatory mechanism for this reduced sensitivity was the presence of non-synonymous point mutations in the target gene Pfcyp51, encoding the sterol 14α-demethylase enzyme. Here, we demonstrate a second mechanism involved in DMI sensitivity of P. fijiensis. We identified a 19-bp element in the wild-type (wt) Pfcyp51 promoter that concatenates in strains with reduced DMI sensitivity. A polymerase chain reaction (PCR) assay identified up to six Pfcyp51 promoter repeats in four field populations of P. fijiensis in Costa Rica. We used transformation experiments to swap the wt promoter of a sensitive field isolate with a promoter from a strain with reduced DMI sensitivity that comprised multiple insertions. Comparative in vivo phenotyping showed a functional and proportional up-regulation of Pfcyp51, which consequently decreased DMI sensitivity. Our data demonstrate that point mutations in the Pfcyp51 coding domain, as well as promoter inserts, contribute to the reduced DMI sensitivity of P. fijiensis. These results provide new insights into the importance of the appropriate use of DMIs and the need for the discovery of new molecules for black Sigatoka management.

CRISPR/Cas has become the state-of-the-art technology for genetic manipulation in diverse organisms, enabling targeted genetic changes to be performed with unprecedented efficiency. Here we report on the first establishment of robust CRISPR/Cas editing in the important necrotrophic plant pathogen Botrytis cinerea based on the introduction of optimized Cas9-sgRNA ribonucleoprotein complexes (RNPs) into protoplasts. Editing yields were further improved by development of a novel strategy that combines RNP delivery with cotransformation of transiently stable vectors containing telomeres, which allowed temporary selection and convenient screening for marker-free editing events. We demonstrate that this approach provides superior editing rates compared to existing CRISPR/Cas-based methods in filamentous fungi, including the model plant pathogen Magnaporthe oryzae. Genome sequencing of edited strains revealed very few additional mutations and no evidence for RNP-mediated off-targeting. The high performance of telomere vector-mediated editing was demonstrated by random mutagenesis of codon 272 of the sdhB gene, a major determinant of resistance to succinate dehydrogenase inhibitor (SDHI) fungicides by in bulk replacement of the codon 272 with codons encoding all 20 amino acids. All exchanges were found at similar frequencies in the absence of selection but SDHI selection allowed the identification of novel amino acid substitutions which conferred differential resistance levels towards different SDHI fungicides. The increased efficiency and easy handling of RNP-based cotransformation is expected to accelerate molecular research in B. cinerea and other fungi.

The plant microbiome represents an enormous untapped resource for discovering novel genes and bioactive compounds. Previously, we isolated Pseudomonas sp. SH-C52 from the rhizosphere of sugar beet plants grown in a soil suppressive to the fungal pathogen Rhizoctonia solani and showed that its antifungal activity is, in part, attributed to the production of the chlorinated 9-amino-acid lipopeptide thanamycin (Mendes et al., 2011). To get more insight into its biosynthetic repertoire, the genome of Pseudomonas sp. SH-C52 was sequenced and subjected to in silico, mutational and functional analyses. The sequencing revealed a genome size of 6.3 Mb and 5579 predicted ORFs. Phylogenetic analysis placed strain SH-C52 within the Pseudomonas corrugata clade. In silico analysis for secondary metabolites revealed a total of six non-ribosomal peptide synthetase (NRPS) gene clusters, including the two previously described NRPS clusters for thanamycin and the 2-amino acid antibacterial lipopeptide brabantamide. Here we show that thanamycin also has activity against an array of other fungi and that brabantamide A exhibits anti-oomycete activity and affects phospholipases of the late blight pathogen Phytophthora infestans. Most notably, mass spectrometry led to the discovery of a third lipopeptide, designated thanapeptin, with a 22-amino-acid peptide moiety. Seven structural variants of thanapeptin were found with varying degrees of activity against P. infestans. Of the remaining four NRPS clusters, one was predicted to encode for yet another and unknown lipopeptide with a predicted peptide moiety of 8-amino acids. Collectively, these results show an enormous metabolic potential for Pseudomonas sp. SH-C52, with at least three structurally diverse lipopeptides, each with a different antimicrobial activity spectrum.

Pythium ultimum is a ubiquitous oomycete plant pathogen responsible for a variety of diseases on a broad range of crop and ornamental species.

In eukaryotes phospholipase D (PLD) is involved in many cellular processes. Currently little is known about PLDs in oomycetes. Here we report that the oomycete plant pathogen Phytophthora infestans has a large repertoire of PLDs divided over six subfamilies: PXPH-PLD, PXTM-PLD, TM-PLD, PLD-likes, and type A and B sPLD-likes. Since the latter have signal peptides we developed a method using metabolically labelled phospholipids to monitor if P. infestans secretes PLD. In extracellular medium of ten P. infestans strains PLD activity was detected as demonstrated by the production of phosphatidic acid and the PLD specific marker phosphatidylalcohol.

The haploid fungus Pseudocercospora fijiensis causes black Sigatoka in banana and is chiefly controlled by extensive fungicide applications, threatening occupational health and the environment. The 14α-Demethylase Inhibitors (DMIs) are important disease control fungicides, but they lose sensitivity in a rather gradual fashion, suggesting an underlying polygenic genetic mechanism. In spite of this, evidence found thus far suggests that P. fijiensis cyp51 gene mutations are the main responsible factor for sensitivity loss in the field. To better understand the mechanisms involved in DMI resistance, in this study we constructed a genetic map using DArTseq markers on two F1 populations generated by crossing two different DMI resistant strains with a sensitive strain. Analysis of the inheritance of DMI resistance in the F1 populations revealed two major and discrete DMI-sensitivity groups. This is an indicative of a single major responsible gene. Using the DMI-sensitivity scorings of both F1 populations and the generation of genetic linkage maps, the sensitivity causal factor was located in a single genetic region. Full agreement was found for genetic markers in either population, underlining the robustness of the approach. The two maps indicated a similar genetic region where the Pfcyp51 gene is found. Sequence analyses of the Pfcyp51 gene of the F1 populations also revealed a matching bimodal distribution with the DMI resistant. Amino acid substitutions in P. fijiensis CYP51 enzyme of the resistant progeny were previously correlated with the loss of DMI sensitivity. In addition, the resistant progeny inherited a Pfcyp51 gene promoter insertion, composed of a repeat element with a palindromic core, also previously correlated with increased gene expression. This genetic approach confirms that Pfcyp51 is the single explanatory gene for reduced sensitivity to DMI fungicides in the analysed P. fijiensis strains. Our study is the first genetic analysis to map the underlying genetic factors for reduced DMI efficacy.

Pathogens deploy a wide range of pathogenicity factors, including a plethora of proteases, to modify host tissue or manipulate host defences. Metalloproteases (MPs) have been implicated in virulence in several animal and plant pathogens. Here we investigated the repertoire of MPs in 46 stramenopile species including 37 oomycetes, 5 diatoms, and 4 brown algae. Screening their complete proteomes using hidden Markov models (HMMs) trained for MP detection resulted in over 4,000 MPs, with most species having between 65 and 100 putative MPs. Classification in clans and families according to the MEROPS database showed a highly diverse MP repertoire in each species. Analyses of domain composition, orthologous groups, distribution, and abundance within the stramenopile lineage revealed a few oomycete-specific MPs and MPs potentially related to lifestyle. In-depth analyses of MPs in the plant pathogen Phytophthora infestans revealed 91 MPs, divided over 21 protein families, including 25 MPs with a predicted signal peptide or signal anchor. Expression profiling showed different patterns of MP gene expression during pre-infection and infection stages. When expressed in leaves of Nicotiana benthamiana, 12 MPs changed the sizes of lesions caused by inoculation with P. infestans; with 9 MPs the lesions were larger, suggesting a positive effect on the virulence of P. infestans, while 3 MPs had a negative effect, resulting in smaller lesions. To the best of our knowledge, this is the first systematic inventory of MPs in oomycetes and the first study pinpointing MPs as potential pathogenicity factors in Phytophthora.

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.