CRISPR-Cas is a recent and powerful addition to the molecular toolbox which allows programmable genome editing. It has been used to modify genes in a wide variety of organisms, but only two alga to date. Here we present a methodology to edit the genome of Thalassiosira pseudonana, a model centric diatom with both ecological significance and high biotechnological potential, using CRISPR-Cas.

Nuclei of arbuscular endomycorrhizal fungi have been described as highly diverse due to their asexual nature and absence of a single cell stage with only one nucleus. This has raised fundamental questions concerning speciation, selection and transmission of the genetic make-up to next generations. Although this concept has become textbook knowledge, it is only based on studying a few loci, including 45S rDNA. To provide a more comprehensive insight into the genetic makeup of arbuscular endomycorrhizal fungi, we applied de novo genome sequencing of individual nuclei of Rhizophagus irregularis. This revealed a surprisingly low level of polymorphism between nuclei. In contrast, within a nucleus, the 45S rDNA repeat unit turned out to be highly diverged. This finding demystifies a long-lasting hypothesis on the complex genetic makeup of arbuscular endomycorrhizal fungi. Subsequent genome assembly resulted in the first draft reference genome sequence of an arbuscular endomycorrhizal fungus. Its length is 141 Mbps, representing over 27,000 protein-coding gene models. We used the genomic sequence to reinvestigate the phylogenetic relationships of Rhizophagus irregularis with other fungal phyla. This unambiguously demonstrated that Glomeromycota are more closely related to Mucoromycotina than to its postulated sister Dikarya.

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.

In 2013, in response to an epidemic of ash dieback disease in England the previous year, we launched a Facebook-based game called Fraxinus to enable non-scientists to contribute to genomics studies of the pathogen that causes the disease and the ash trees that are devastated by it. Over a period of 51 weeks players were able to match computational alignments of genetic sequences in 78% of cases, and to improve them in 15% of cases. We also found that most players were only transiently interested in the game, and that the majority of the work done was performed by a small group of dedicated players. Based on our experiences we have built a linear model for the length of time that contributors are likely to donate to a crowd-sourced citizen science project. This model could serve a guide for the design and implementation of future crowd-sourced citizen science initiatives.

The oomycete pathogen Phytophthora infestans causes potato late blight, and as a potato and tomato specialist pathogen, is seemingly poorly adapted to infect plants outside the Solanaceae. Here, we report the unexpected finding that P. infestans can infect Arabidopsis thaliana when another oomycete pathogen, Albugo laibachii, has colonized the host plant. The behaviour and speed of P. infestans infection in Arabidopsis pre-infected with A. laibachii resemble P. infestans infection of susceptible potato plants. Transcriptional profiling of P. infestans genes during infection revealed a significant overlap in the sets of secreted-protein genes that are induced in P. infestans upon colonization of potato and susceptible Arabidopsis, suggesting major similarities in P. infestans gene expression dynamics on the two plant species. Furthermore, we found haustoria of A. laibachii and P. infestans within the same Arabidopsis cells. This Arabidopsis-A. laibachii-P. infestans tripartite interaction opens up various possibilities to dissect the molecular mechanisms of P. infestans infection and the processes occurring in co-infected Arabidopsis cells.

Rust fungal pathogens of wheat (Triticum spp.) affect crop yields worldwide. The molecular mechanisms underlying the virulence of these pathogens remain elusive, due to the limited availability of suitable molecular genetic research tools. Notably, the inability to perform high-throughput analyses of candidate virulence proteins (also known as effectors) impairs progress. We previously established a pipeline for the fast-forward screens of rust fungal candidate effectors in the model plant Nicotiana benthamiana. This pipeline involves selecting candidate effectors in silico and performing cell biology and protein-protein interaction assays in planta to gain insight into the putative functions of candidate effectors. In this study, we used this pipeline to identify and characterize sixteen candidate effectors from the wheat yellow rust fungal pathogen Puccinia striiformis f sp tritici. Nine candidate effectors targeted a specific plant subcellular compartment or protein complex, providing valuable information on their putative functions in plant cells. One candidate effector, PST02549, accumulated in processing bodies (P-bodies), protein complexes involved in mRNA decapping, degradation, and storage. PST02549 also associates with the P-body-resident ENHANCER OF mRNA DECAPPING PROTEIN 4 (EDC4) from N. benthamiana and wheat. We propose that P-bodies are a novel plant cell compartment targeted by pathogen effectors.

Animal and plant eukaryotic pathogens, such as the human malaria parasite Plasmodium falciparum and the potato late blight agent Phytophthora infestans, are widely divergent eukaryotic microbes. Yet they both produce secretory virulence and pathogenic proteins that alter host cell functions. In P. falciparum, export of parasite proteins to the host erythrocyte is mediated by leader sequences shown to contain a host-targeting (HT) motif centered on an RxLx (E, D, or Q) core: this motif appears to signify a major pathogenic export pathway with hundreds of putative effectors. Here we show that a secretory protein of P. infestans, which is perceived by plant disease resistance proteins and induces hypersensitive plant cell death, contains a leader sequence that is equivalent to the Plasmodium HT-leader in its ability to export fusion of green fluorescent protein (GFP) from the P. falciparum parasite to the host erythrocyte. This export is dependent on an RxLR sequence conserved in P. infestans leaders, as well as in leaders of all ten secretory oomycete proteins shown to function inside plant cells. The RxLR motif is also detected in hundreds of secretory proteins of P. infestans, Phytophthora sojae, and Phytophthora ramorum and has high value in predicting host-targeted leaders. A consensus motif further reveals E/D residues enriched within approximately 25 amino acids downstream of the RxLR, which are also needed for export. Together the data suggest that in these plant pathogenic oomycetes, a consensus HT motif may reside in an extended sequence of approximately 25-30 amino acids, rather than in a short linear sequence. Evidence is presented that although the consensus is much shorter in P. falciparum, information sufficient for vacuolar export is contained in a region of approximately 30 amino acids, which includes sequences flanking the HT core. Finally, positional conservation between Phytophthora RxLR and P. falciparum RxLx (E, D, Q) is consistent with the idea that the context of their presentation is constrained. These studies provide the first evidence to our knowledge that eukaryotic microbes share equivalent pathogenic HT signals and thus conserved mechanisms to access host cells across plant and animal kingdoms that may present unique targets for prophylaxis across divergent pathogens.

Unconventional integrated domains in plant intracellular immune receptors of the nucleotide-binding leucine-rich repeat (NLRs) type can directly bind translocated effector proteins from pathogens and thereby initiate an immune response. The rice (Oryza sativa) immune receptor pairs Pik-1/Pik-2 and RGA5/RGA4 both use integrated heavy metal-associated (HMA) domains to bind the effectors AVR-Pik and AVR-Pia, respectively, from the rice blast fungal pathogen Magnaporthe oryzae These effectors both belong to the MAX effector family and share a core structural fold, despite being divergent in sequence. How integrated domains in NLRs maintain specificity of effector recognition, even of structurally similar effectors, has implications for understanding plant immune receptor evolution and function. Here, using plant cell death and pathogenicity assays and protein-protein interaction analyses, we show that the rice NLR pair Pikp-1/Pikp-2 triggers an immune response leading to partial disease resistance toward the "mis-matched" effector AVR-Pia in planta and that the Pikp-HMA domain binds AVR-Pia in vitro We observed that the HMA domain from another Pik-1 allele, Pikm, cannot bind AVR-Pia, and it does not trigger a plant response. The crystal structure of Pikp-HMA bound to AVR-Pia at 1.9 Å resolution revealed a binding interface different from those formed with AVR-Pik effectors, suggesting plasticity in integrated domain-effector interactions. The results of our work indicate that a single NLR immune receptor can bait multiple pathogen effectors via an integrated domain, insights that may enable engineering plant immune receptors with extended disease resistance profiles.

Plant NLR immune receptors are multidomain proteins that can function as specialized sensor/helper pairs. Paired NLR immune receptors are generally thought to function via negative regulation, where one NLR represses the activity of the second and detection of pathogen effectors relieves this repression to initiate immunity. However, whether this mechanism is common to all NLR pairs is not known. Here, we show that the rice NLR pair Pikp-1/Pikp-2, which confers resistance to strains of the blast pathogen Magnaporthe oryzae (syn. Pyricularia oryzae) expressing the AVR-PikD effector, functions via receptor cooperation, with effector-triggered activation requiring both NLRs to trigger the immune response. To investigate the mechanism of Pikp-1/Pikp-2 activation, we expressed truncated variants of these proteins, and made mutations in previously identified NLR sequence motifs. We found that any domain truncation, in either Pikp-1 or Pikp-2, prevented cell death in the presence of AVR-PikD, revealing that all domains are required for activity. Further, expression of individual Pikp-1 or Pikp-2 domains did not result in cell death. Mutations in the conserved P-loop and MHD sequence motifs in both Pikp-1 and Pikp-2 prevented cell death activation, demonstrating that these motifs are required for the function of the two partner NLRs. Finally, we showed that Pikp-1 and Pikp-2 associate to form homo- and hetero-complexes in planta in the absence of AVR-PikD; on co-expression the effector binds to Pikp-1 generating a tri-partite complex. Taken together, we provide evidence that Pikp-1 and Pikp-2 form a fine-tuned system that is activated by AVR-PikD via receptor cooperation rather than negative regulation.

Accelerated gene evolution is a hallmark of pathogen adaptation and specialization following host-jumps. However, the molecular processes associated with adaptive evolution between host-specific lineages of a multihost plant pathogen remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), host specialization on different grass hosts is generally associated with dynamic patterns of gain and loss of virulence effector genes that tend to define the distinct genetic lineages of this pathogen. Here, we unravelled the biochemical and structural basis of adaptive evolution of APikL2, an exceptionally conserved paralog of the well-studied rice-lineage specific effector AVR-Pik. Whereas AVR-Pik and other members of the six-gene AVR-Pik family show specific patterns of presence/absence polymorphisms between grass-specific lineages of M. oryzae, APikL2 stands out by being ubiquitously present in all blast fungus lineages from 13 different host species. Using biochemical, biophysical and structural biology methods, we show that a single aspartate to asparagine polymorphism expands the binding spectrum of APikL2 to host proteins of the heavy-metal associated (HMA) domain family. This mutation maps to one of the APikL2-HMA binding interfaces and contributes to an altered hydrogen-bonding network. By combining phylogenetic ancestral reconstruction with an analysis of the structural consequences of allelic diversification, we revealed a common mechanism of effector specialization in the AVR-Pik/APikL2 family that involves two major HMA-binding interfaces. Together, our findings provide a detailed molecular evolution and structural biology framework for diversification and adaptation of a fungal pathogen effector family following host-jumps.

In plants, nucleotide-binding domain and leucine-rich repeat (NLR)-containing proteins can form receptor networks to confer hypersensitive cell death and innate immunity. One class of NLRs, known as NLR required for cell death (NRCs), are central nodes in a complex network that protects against multiple pathogens and comprises up to half of the NLRome of solanaceous plants. Given the prevalence of this NLR network, we hypothesised that pathogens convergently evolved to secrete effectors that target NRC activities. To test this, we screened a library of 165 bacterial, oomycete, nematode, and aphid effectors for their capacity to suppress the cell death response triggered by the NRC-dependent disease resistance proteins Prf and Rpi-blb2. Among 5 of the identified suppressors, 1 cyst nematode protein and 1 oomycete protein suppress the activity of autoimmune mutants of NRC2 and NRC3, but not NRC4, indicating that they specifically counteract a subset of NRC proteins independently of their sensor NLR partners. Whereas the cyst nematode effector SPRYSEC15 binds the nucleotide-binding domain of NRC2 and NRC3, the oomycete effector AVRcap1b suppresses the response of these NRCs via the membrane trafficking-associated protein NbTOL9a (Target of Myb 1-like protein 9a). We conclude that plant pathogens have evolved to counteract central nodes of the NRC immune receptor network through different mechanisms. Coevolution with pathogen effectors may have driven NRC diversification into functionally redundant nodes in a massively expanded NLR network.

Nucleotide-binding domain leucine-rich repeat (NLR) immune receptors are important components of plant and metazoan innate immunity that can function as individual units or as pairs or networks. Upon activation, NLRs form multiprotein complexes termed resistosomes or inflammasomes. Although metazoan paired NLRs, such as NAIP/NLRC4, form hetero-complexes upon activation, the molecular mechanisms underpinning activation of plant paired NLRs, especially whether they associate in resistosome hetero-complexes, is unknown. In asterid plant species, the NLR required for cell death (NRC) immune receptor network is composed of multiple resistance protein sensors and downstream helpers that confer immunity against diverse plant pathogens. Here, we show that pathogen effector-activation of the NLR proteins Rx (confers virus resistance), and Bs2 (confers bacterial resistance) leads to oligomerization of their helper NLR, NRC2. Activated Rx does not oligomerize or enter into a stable complex with the NRC2 oligomer and remains cytoplasmic. In contrast, activated NRC2 oligomers accumulate in membrane-associated puncta. We propose an activation-and-release model for NLRs in the NRC immune receptor network. This points to a distinct activation model compared with mammalian paired NLRs.

In eukaryotes, targeted protein degradation (TPD) typically depends on a series of interactions among ubiquitin ligases that transfer ubiquitin molecules to substrates leading to degradation by the 26S proteasome. We previously identified that the bacterial effector protein SAP05 mediates ubiquitin-independent TPD. SAP05 forms a ternary complex via interactions with the von Willebrand Factor Type A (vWA) domain of the proteasomal ubiquitin receptor Rpn10 and the zinc-finger (ZnF) domains of the SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) and GATA BINDING FACTOR (GATA) transcription factors (TFs). This leads to direct TPD of the TFs by the 26S proteasome. Here, we report the crystal structures of the SAP05-Rpn10vWA complex at 2.17 Å resolution and of the SAP05-SPL5ZnF complex at 2.20 Å resolution. Structural analyses revealed that SAP05 displays a remarkable bimodular architecture with two distinct nonoverlapping surfaces, a "loop surface" with three protruding loops that form electrostatic interactions with ZnF, and a "sheet surface" featuring two β-sheets, loops, and α-helices that establish polar interactions with vWA. SAP05 binding to ZnF TFs involves single amino acids responsible for multiple contacts, while SAP05 binding to vWA is more stable due to the necessity of multiple mutations to break the interaction. In addition, positioning of the SAP05 complex on the 26S proteasome points to a mechanism of protein degradation. Collectively, our findings demonstrate how a small bacterial bimodular protein can bypass the canonical ubiquitin-proteasome proteolysis pathway, enabling ubiquitin-independent TPD in eukaryotic cells. This knowledge holds significant potential for the creation of TPD technologies.

Emerging and re-emerging pathogens imperil public health and global food security. Responding to these threats requires improved surveillance and diagnostic systems. Despite their potential, genomic tools have not been readily applied to emerging or re-emerging plant pathogens such as the wheat yellow (stripe) rust pathogen Puccinia striiformis f. sp. tritici (PST). This is due largely to the obligate parasitic nature of PST, as culturing PST isolates for DNA extraction remains slow and tedious.

The filamentous oomycete plant pathogen Phytophthora infestans causes late blight, an economically important disease, on members of the nightshade family (Solanaceae), such as the crop plants potato and tomato. The related plant Nicotiana benthamiana is a model system to study plant-pathogen interactions, and the susceptibility of N. benthamiana to Phytophthora species varies from susceptible to resistant. Little is known about the extent to which plant basal immunity, mediated by membrane receptors that recognise conserved pathogen-associated molecular patterns (PAMPs), contributes to P. infestans resistance.

Rust fungi are obligate biotrophic pathogens that cause considerable damage on crop plants. Puccinia graminis f. sp. tritici, the causal agent of wheat stem rust, and Melampsora larici-populina, the poplar leaf rust pathogen, have strong deleterious impacts on wheat and poplar wood production, respectively. Filamentous pathogens such as rust fungi secrete molecules called disease effectors that act as modulators of host cell physiology and can suppress or trigger host immunity. Current knowledge on effectors from other filamentous plant pathogens can be exploited for the characterisation of effectors in the genome of recently sequenced rust fungi. We designed a comprehensive in silico analysis pipeline to identify the putative effector repertoire from the genome of two plant pathogenic rust fungi. The pipeline is based on the observation that known effector proteins from filamentous pathogens have at least one of the following properties: (i) contain a secretion signal, (ii) are encoded by in planta induced genes, (iii) have similarity to haustorial proteins, (iv) are small and cysteine rich, (v) contain a known effector motif or a nuclear localization signal, (vi) are encoded by genes with long intergenic regions, (vii) contain internal repeats, and (viii) do not contain PFAM domains, except those associated with pathogenicity. We used Markov clustering and hierarchical clustering to classify protein families of rust pathogens and rank them according to their likelihood of being effectors. Using this approach, we identified eight families of candidate effectors that we consider of high value for functional characterization. This study revealed a diverse set of candidate effectors, including families of haustorial expressed secreted proteins and small cysteine-rich proteins. This comprehensive classification of candidate effectors from these devastating rust pathogens is an initial step towards probing plant germplasm for novel resistance components.

Rust fungi cause serious yield reductions on crops, including wheat, barley, soybean, coffee, and represent real threats to global food security. Of these fungi, the flax rust pathogen Melampsora lini has been developed most extensively over the past 80 years as a model to understand the molecular mechanisms that underpin pathogenesis. During infection, M. lini secretes virulence effectors to promote disease. The number of these effectors, their function and their degree of conservation across rust fungal species is unknown. To assess this, we sequenced and assembled de novo the genome of M. lini isolate CH5 into 21,130 scaffolds spanning 189 Mbp (scaffold N50 of 31 kbp). Global analysis of the DNA sequence revealed that repetitive elements, primarily retrotransposons, make up at least 45% of the genome. Using ab initio predictions, transcriptome data and homology searches, we identified 16,271 putative protein-coding genes. An analysis pipeline was then implemented to predict the effector complement of M. lini and compare it to that of the poplar rust, wheat stem rust and wheat stripe rust pathogens to identify conserved and species-specific effector candidates. Previous knowledge of four cloned M. lini avirulence effector proteins and two basidiomycete effectors was used to optimize parameters of the effector prediction pipeline. Markov clustering based on sequence similarity was performed to group effector candidates from all four rust pathogens. Clusters containing at least one member from M. lini were further analyzed and prioritized based on features including expression in isolated haustoria and infected leaf tissue and conservation across rust species. Herein, we describe 200 of 940 clusters that ranked highest on our priority list, representing 725 flax rust candidate effectors. Our findings on this important model rust species provide insight into how effectors of rust fungi are conserved across species and how they may act to promote infection on their hosts.

Understanding the mechanisms and timescales of plant pathogen outbreaks requires a detailed genome-scale analysis of their population history. The fungus Magnaporthe (Syn. Pyricularia) oryzae-the causal agent of blast disease of cereals- is among the most destructive plant pathogens to world agriculture and a major threat to the production of rice, wheat, and other cereals. Although M. oryzae is a multihost pathogen that infects more than 50 species of cereals and grasses, all rice-infecting isolates belong to a single genetically defined lineage. Here, we combined the two largest genomic datasets to reconstruct the genetic history of the rice-infecting lineage of M. oryzae based on 131 isolates from 21 countries.

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

Cooperation between receptors from the nucleotide-binding, leucine-rich repeats (NLR) superfamily is important for intracellular activation of immune responses. NLRs can function in pairs that, upon pathogen recognition, trigger hypersensitive cell death and stop pathogen invasion. Natural selection drives specialization of host immune receptors towards an optimal response, whilst keeping a tight regulation of immunity in the absence of pathogens. However, the molecular basis of co-adaptation and specialization between paired NLRs remains largely unknown. Here, we describe functional specialization in alleles of the rice NLR pair Pik that confers resistance to strains of the blast fungus Magnaporthe oryzae harbouring AVR-Pik effectors. We revealed that matching pairs of allelic Pik NLRs mount effective immune responses, whereas mismatched pairs lead to autoimmune phenotypes, a hallmark of hybrid necrosis in both natural and domesticated plant populations. We further showed that allelic specialization is largely underpinned by a single amino acid polymorphism that determines preferential association between matching pairs of Pik NLRs. These results provide a framework for how functionally linked immune receptors undergo co-adaptation to provide an effective and regulated immune response against pathogens. Understanding the molecular constraints that shape paired NLR evolution has implications beyond plant immunity given that hybrid necrosis can drive reproductive isolation.