Publications

Plant pathogens compromise crop yields. Plants have evolved robust innate immunity that depends in part on intracellular Nucleotide-binding, Leucine Rich-Repeat (NLR) immune receptors that activate defense responses upon detection of pathogen-derived effectors. Most sensor NLRs that detect effectors require the activity of helper NLRs, but how helper NLRs support sensor NLR function is poorly understood. Many Solanaceae NLRs require the NRC (NLR-Required for Cell death) class of helper NLRs. We show here that Rpi-amr3, a sensor NLR from Solanum americanum, detects AVRamr3 from the potato late blight pathogen, Phytophthora infestans, and activates oligomerization of the helper NLR NRC2 into a high-molecular weight resistosome. The NRC2 resistosome also forms upon recognition of P. infestans effector AVRamr1 by another sensor NLR, Rpi-amr1. The ATP-binding motif of Rpi-amr3 is required for NRC2 resistosome formation, but not for interaction with the cognate effector. The NRC2 resistosome can be activated by AVRamr3 homologs from other Phytophthora species. Mechanistic understanding of NRC resistosome formation will underpin engineering crops with durable disease resistance.

The NRC immune receptor network has evolved in asterid plants from a pair of linked genes into a genetically dispersed and phylogenetically structured network of sensor and helper NLR (nucleotide-binding domain and leucine-rich repeat-containing) proteins. In some species, such as the model plant Nicotiana benthamiana and other Solanaceae, the NRC network forms up to half of the NLRome, and NRCs are scattered throughout the genome in gene clusters of varying complexities. Here, we describe NRCX, an atypical, but essential member of the NRC family that lacks canonical features of these NLR helper proteins, such as a functional N-terminal MADA motif and the capacity to trigger autoimmunity. In contrast to other NRCs, systemic gene silencing of NRCX markedly impairs plant growth resulting in a dwarf phenotype. Remarkably, dwarfism of NRCX silenced plants is partially dependent on NRCX paralogs NRC2 and NRC3, but not NRC4. Despite its negative impact on plant growth when silenced systemically, transient RNA interference of NRCX in mature N. benthamiana leaves doesn’t result in visible cell death phenotypes. However, alteration of NRCX expression modulates the hypersensitive response mediated by NRC2 and NRC3 in a manner consistent with a negative role for NRCX in the NRC network. We conclude that NRCX is an atypical member of the NRC network that has evolved to contribute to the homeostasis of this genetically unlinked NLR network.

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.

Plant pathogens cause recurrent epidemics that threaten crop yield and global food security. Efforts to retool the plant immune system have been limited to modifying natural components and can be nullified by the emergence of new pathogen races. Therefore, there is a need to develop made-to-order synthetic plant immune receptors with resistance tailored to the pathogen genotypes present in the field. Here we show that plant immune receptors can be used as scaffolds for VHH nanobody fusions that bind fluorescent proteins (FPs). The receptor-nanobody fusions signal in the presence of the corresponding FP and confer resistance against plant viruses expressing FPs. Given that nanobodies can be raised against virtually any molecule, immune receptor-nanobody fusions have the potential to generate resistance against all major plant pathogens and pests.

In order to infect a new host species, the pathogen must evolve to enhance infection and transmission in the novel environment. Although we often think of evolution as a process of accumulation, it is also a process of loss. Here, we document an example of regressive evolution in the Irish potato famine pathogen (Phytophthora infestans) lineage, providing evidence that a key sequence motif in the effector PexRD54 has degenerated following a host jump. We began by looking at PexRD54 and PexRD54-like sequences from across Phytophthora species. We found that PexRD54 emerged in the common ancestor of Phytophthora clade 1b and 1c species, and further sequence analysis showed that a key functional motif, the C-terminal ATG8-interacting motif (AIM), was also acquired at this point in the lineage. A closer analysis showed that the P. mirabilis PexRD54 (PmPexRD54) AIM appeared unusual, the otherwise-conserved central residue mutated from a glutamate to a lysine. We aimed to determine whether this PmPexRD54 AIM polymorphism represented an adaptation to the Mirabilis jalapa host environment. We began by characterizing the M. jalapa ATG8 family, finding that they have a unique evolutionary history compared to previously characterized ATG8s. Then, using co-immunoprecipitation and isothermal titration calorimetry assays, we showed that both full-length PmPexRD54 and the PmPexRD54 AIM peptide bind very weakly to the M. jalapa ATG8s. Through a combination of binding assays and structural modelling, we showed that the identity of the residue at the position of the PmPexRD54 AIM polymorphism can underpin high-affinity binding to plant ATG8s. Finally, we conclude that the functionality of the PexRD54 AIM was lost in the P. mirabilis lineage, perhaps owing to as-yet-unknown pressure on this effector in the new host environment.

Objectives Accurate single-nucleotide polymorphisms (SNP) calls are crucial for robust evolutionary and population genetic inferences in genomic analyses. Such inferences can reveal the time-scales and processes associated with the emergence and spread of pandemic plant pathogens, such as the rice blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae). However, the specificity and sensitivity of SNP calls depend on the filtering parameters applied to the data. Here, we used a benchmarking approach to evaluate the impact of SNP calling on different population genetic analyses of the rice blast fungus, namely genetic clustering, topology of phylogenetic reconstructions and estimation of evolutionary rates. Results To benchmark SNP calling parameters, we generated a gold standard set of validated SNPs by sequencing nine M. oryzae genomes with both Illumina short-reads and Oxford Nanopore Technologies (ONT). We used the gold standard set of SNPs to identify the SNP calling parameter configuration that maximizes sensitivity and specificity. We found that the choice of parameter configurations can substantially change the number of ascertained SNPs, preferentially affecting SNPs segregating at low population frequency. However, SNP calling parameter configurations did not significantly affect the clustering of isolates in clonal lineages, the monophyly of each clonal lineage, and the estimation of evolutionary rates. We leverage the evolutionary rates obtained from each SNP calling parameter configuration to generate divergence time estimates that take into account the uncertainty associated with both the estimation of evolutionary rates and SNP calling. Our analysis indicates that M. oryzae clonal lineage expansions took place ~300 years ago.