Dominant integration locus drives continuous diversification of plant immune receptors with exogenous domain fusions

A haplotype-resolved chromosome-level assembly and annotation of European hazelnut (C. avellana cv. Jefferson) provides insight into mechanisms of eastern filbert blight resistance

European hazelnut (Corylus avellana L.) is an important tree nut crop. Hazelnut production in North America is currently limited in scalability due to Anisogramma anomala, a fungal pathogen that causes Eastern Filbert Blight (EFB) disease in hazelnut. Successful deployment of EFB resistant cultivars has been limited to the state of Oregon, where the breeding program at Oregon State University (OSU) has released cultivars with a dominant allele at a single resistance locus identified by classical breeding, linkage mapping, and molecular markers. C. avellana cultivar "Jefferson" is resistant to the predominant EFB biotype in Oregon and has been selected by the OSU breeding program as a model for hazelnut genetic and genomic research. Here, we present a near complete, haplotype-resolved chromosome-level hazelnut genome assembly for "Jefferson". This new assembly is a significant improvement over a previously published genome draft. Analysis of genomic regions linked to EFB resistance and self-incompatibility confirmed haplotype splitting and identified new gene candidates that are essential for downstream molecular marker development, thereby facilitating breeding efforts.

Zinc-finger (ZiF) fold secreted effectors form a functionally diverse family across lineages of the blast fungus Magnaporthe oryzae

Filamentous plant pathogens deliver effector proteins into host cells to suppress host defence responses and manipulate metabolic processes to support colonization. Understanding the evolution and molecular function of these effectors provides knowledge about pathogenesis and can suggest novel strategies to reduce damage caused by pathogens. However, effector proteins are highly variable, share weak sequence similarity and, although they can be grouped according to their structure, only a few structurally conserved effector families have been functionally characterized to date. Here, we demonstrate that Zinc-finger fold (ZiF) secreted proteins form a functionally diverse effector family in the blast fungus Magnaporthe oryzae. This family relies on the Zinc-finger motif for protein stability and is ubiquitously present in blast fungus lineages infecting 13 different host species, forming different effector tribes. Homologs of the canonical ZiF effector, AVR-Pii, from rice infecting isolates are present in multiple M. oryzae lineages. Wheat infecting strains of the fungus also possess an AVR-Pii like allele that binds host Exo70 proteins and activates the immune receptor Pii. Furthermore, ZiF tribes may vary in the proteins they bind to, indicating functional diversification and an intricate effector/host interactome. Altogether, we uncovered a new effector family with a common protein fold that has functionally diversified in lineages of M. oryzae. This work expands our understanding of the diversity of M. oryzae effectors, the molecular basis of plant pathogenesis and may ultimately facilitate the development of new sources for pathogen resistance.

The unconventional resistance protein PTR recognizes the Magnaporthe oryzae effector AVR-Pita in an allele-specific manner

Blast disease caused by the fungus Magnaporthe oryzae is one of the most devastating rice diseases. Disease resistance genes such as Pi-ta or Pi-ta2 are critical in protecting rice production from blast. Published work reports that Pi-ta codes for a nucleotide-binding and leucine-rich repeat domain protein (NLR) that recognizes the fungal protease-like effector AVR-Pita by direct binding. However, this model was challenged by the recent discovery that Pi-ta2 resistance, which also relies on AVR-Pita detection, is conferred by the unconventional resistance gene Ptr, which codes for a membrane protein with a cytoplasmic armadillo repeat domain. Here, using NLR Pi-ta and Ptr RNAi knockdown and CRISPR/Cas9 knockout mutant rice lines, we found that AVR-Pita recognition relies solely on Ptr and that the NLR Pi-ta has no role in it, indicating that it is not the Pi-ta resistance gene. Different alleles of Ptr confer different recognition specificities. The A allele of Ptr (PtrA) detects all natural sequence variants of the effector and confers Pi-ta2 resistance, while the B allele of Ptr (PtrB) recognizes a restricted set of AVR-Pita alleles and, thereby, confers Pi-ta resistance. Analysis of the natural diversity in AVR-Pita and of mutant and transgenic strains identified one specific polymorphism in the effector sequence that controls escape from PtrB-mediated resistance. Taken together, our work establishes that the M. oryzae effector AVR-Pita is detected in an allele-specific manner by the unconventional rice resistance protein Ptr and that the NLR Pi-ta has no function in Pi-ta resistance and the recognition of AVR-Pita.

Barley MLA3 recognizes the host-specificity effector Pwl2 from Magnaporthe oryzae

Plant nucleotide-binding leucine-rich repeat (NLRs) immune receptors directly or indirectly recognize pathogen-secreted effector molecules to initiate plant defense. Recognition of multiple pathogens by a single NLR is rare and usually occurs via monitoring for changes to host proteins; few characterized NLRs have been shown to recognize multiple effectors. The barley (Hordeum vulgare) NLR gene Mildew locus a (Mla) has undergone functional diversification, and the proteins encoded by different Mla alleles recognize host-adapted isolates of barley powdery mildew (Blumeria graminis f. sp. hordei [Bgh]). Here, we show that Mla3 also confers resistance to the rice blast fungus Magnaporthe oryzae in a dosage-dependent manner. Using a forward genetic screen, we discovered that the recognized effector from M. oryzae is Pathogenicity toward Weeping Lovegrass 2 (Pwl2), a host range determinant factor that prevents M. oryzae from infecting weeping lovegrass (Eragrostis curvula). Mla3 has therefore convergently evolved the capacity to recognize effectors from diverse pathogens.

The exocyst complex is an evolutionary battleground in plant-microbe interactions

Exocytosis is a conserved trafficking pathway that transports secretory vesicles to the extracellular space, replenishes the plasma membrane and is essential for establishing cell polarity. Its spatiotemporal regulation is mediated by an evolutionary conserved octameric tethering complex, the exocyst. In plants, certain subunits of this complex have diversified and acquired multiple functions, including a central role in defense against pathogens and pests. Here, I review the latest evidence suggesting the dramatic expansion and functional diversification of the exocyst subunit Exo70 is likely driven by a coevolutionary arms race, in which Exo70 proteins are repeatedly targeted by effectors from multiple pathogens and, in turn, are monitored by plant immune receptors for pathogen perception.

NLR receptors in plant immunity: making sense of the alphabet soup

Plants coordinately use cell-surface and intracellular immune receptors to perceive pathogens and mount an immune response. Intracellular events of pathogen recognition are largely mediated by immune receptors of the nucleotide binding and leucine rich-repeat (NLR) classes. Upon pathogen perception, NLRs trigger a potent broad-spectrum immune reaction, usually accompanied by a form of programmed cell death termed the hypersensitive response. Some plant NLRs act as multifunctional singleton receptors which combine pathogen detection and immune signaling. However, NLRs can also function in higher order pairs and networks of functionally specialized interconnected receptors. In this article, we cover the basic aspects of plant NLR biology with an emphasis on NLR networks. We highlight some of the recent advances in NLR structure, function, and activation and discuss emerging topics such as modulator NLRs, pathogen suppression of NLRs, and NLR bioengineering. Multi-disciplinary approaches are required to disentangle how these NLR immune receptor pairs and networks function and evolve. Answering these questions holds the potential to deepen our understanding of the plant immune system and unlock a new era of disease resistance breeding.

Asymmetric Evolution of Protein Domains in the Leucine-Rich Repeat Receptor-Like Kinase Family of Plant Signaling Proteins

The coding sequences of developmental genes are expected to be deeply conserved, with cis-regulatory change driving the modulation of gene function. In contrast, proteins with roles in defense are expected to evolve rapidly, in molecular arms races with pathogens. However, some gene families include both developmental and defense genes. In these families, does the tempo and mode of evolution differ between genes with divergent functions, despite shared ancestry and structure? The leucine-rich repeat receptor-like kinase (LRR-RLKs) protein family includes members with roles in plant development and defense, thus providing an ideal system for answering this question. LRR-RLKs are receptors that traverse plasma membranes. LRR domains bind extracellular ligands; RLK domains initiate intracellular signaling cascades in response to ligand binding. In LRR-RLKs with roles in defense, LRR domains evolve faster than RLK domains. To determine whether this asymmetry extends to LRR-RLKs that function primarily in development, we assessed evolutionary rates and tested for selection acting on 11 subfamilies of LRR-RLKs, using deeply sampled protein trees. To assess functional evolution, we performed heterologous complementation assays in Arabidopsis thaliana (Arabidopsis). We found that the LRR domains of all tested LRR-RLK proteins evolved faster than their cognate RLK domains. All tested subfamilies of LRR-RLKs had strikingly similar patterns of molecular evolution, despite divergent functions. Heterologous transformation experiments revealed that multiple mechanisms likely contribute to the evolution of LRR-RLK function, including escape from adaptive conflict. Our results indicate specific and distinct evolutionary pressures acting on LRR versus RLK domains, despite diverse organismal roles for LRR-RLK proteins.

An NLR integrated domain toolkit to identify plant pathogen effector targets

Plant immune receptors, known as NOD-like receptors (NLRs), possess unique integrated decoy domains that enable plants to attract pathogen effectors and initiate a specific immune response. The present study aimed to create a library of these integrated domains (IDs) and screen them with pathogen effectors to identify targets for effector virulence and NLR-effector interactions. This works compiles IDs found in NLRs from seven different plant species and produced a library of 78 plasmid clones containing a total of 104 IDs, representing 43 distinct InterPro domains. A yeast two-hybrid assay was conducted, followed by an in planta interaction test, using 32 conserved effectors from Ralstonia pseudosolanacearum type III. Through these screenings, three interactions involving different IDs (kinase, DUF3542, WRKY) were discovered interacting with two unrelated type III effectors (RipAE and PopP2). Of particular interest was the interaction between PopP2 and ID#85, an atypical WRKY domain integrated into a soybean NLR gene (GmNLR-ID#85). Using a Förster resonance energy transfer-fluorescence lifetime imaging microscopy technique to detect protein-protein interactions in living plant cells, PopP2 was demonstrated to physically associate with ID#85 in the nucleus. However, unlike the known WRKY-containing Arabidopsis RRS1-R NLR receptor, GmNLR-ID#85 could not be acetylated by PopP2 and failed to activate RPS4-dependent immunity when introduced into the RRS1-R immune receptor. The generated library of 78 plasmid clones, encompassing these screenable IDs, is publicly available through Addgene. This resource is expected to be valuable for the scientific community with respect to discovering targets for effectors and potentially engineering plant immune receptors.

Structural mechanism of heavy metal-associated integrated domain engineering of paired nucleotide-binding and leucine-rich repeat proteins in rice

Plant nucleotide-binding and leucine-rich repeat (NLR) proteins are immune sensors that detect pathogen effectors and initiate a strong immune response. In many cases, single NLR proteins are sufficient for both effector recognition and signaling activation. These proteins possess a conserved architecture, including a C-terminal leucine-rich repeat (LRR) domain, a central nucleotide-binding (NB) domain, and a variable N-terminal domain. Nevertheless, many paired NLRs linked in a head-to-head configuration have now been identified. The ones carrying integrated domains (IDs) can recognize pathogen effector proteins by various modes; these are known as sensor NLR (sNLR) proteins. Structural and biochemical studies have provided insights into the molecular basis of heavy metal-associated IDs (HMA IDs) from paired NLRs in rice and revealed the co-evolution between pathogens and hosts by combining naturally occurring favorable interactions across diverse interfaces. Focusing on structural and molecular models, here we highlight advances in structure-guided engineering to expand and enhance the response profile of paired NLR-HMA IDs in rice to variants of the rice blast pathogen MAX-effectors (Magnaporthe oryzae AVRs and ToxB-like). These results demonstrate that the HMA IDs-based design of rice materials with broad and enhanced resistance profiles possesses great application potential but also face considerable challenges.

The integrated LIM-peptidase domain of the CSA1-CHS3/DAR4 paired immune receptor detects changes in DA1 peptidase inhibitors in Arabidopsis

White blister rust, caused by the oomycete Albugo candida, is a widespread disease of Brassica crops. The Brassica relative Arabidopsis thaliana uses the paired immune receptor complex CSA1-CHS3/DAR4 to resist Albugo infection. The CHS3/DAR4 sensor NLR, which functions together with its partner, the helper NLR CSA1, carries an integrated domain (ID) with homology to DA1 peptidases. Using domain swaps with several DA1 homologs, we show that the LIM-peptidase domain of the family member CHS3/DAR4 functions as an integrated decoy for the family member DAR3, which interacts with and inhibits the peptidase activities of the three closely related peptidases DA1, DAR1, and DAR2. Albugo infection rapidly lowers DAR3 levels and activates DA1 peptidase activity, thereby promoting endoreduplication of host tissues to support pathogen growth. We propose that the paired immune receptor CSA1-CHS3/DAR4 detects the actions of a putative Albugo effector that reduces DAR3 levels, resulting in defense activation.

Molecular Mechanisms of the Co-Evolution of Wheat and Rust Pathogens

Wheat (Triticum spp.) is a cereal crop domesticated >8000 years ago and the second-most-consumed food crop nowadays. Ever since mankind has written records, cereal rust diseases have been a painful awareness in antiquity documented in the Old Testament (about 750 B.C.). The pathogen causing the wheat stem rust disease is among the first identified plant pathogens in the 1700s, suggesting that wheat and rust pathogens have co-existed for thousands of years. With advanced molecular technologies, wheat and rust genomes have been sequenced, and interactions between the host and the rust pathogens have been extensively studied at molecular levels. In this review, we summarized the research at the molecular level and organized the findings based on the pathogenesis steps of germination, penetration, haustorial formation, and colonization of the rusts to present the molecular mechanisms of the co-evolution of wheat and rust pathogens.

Engineering effector-triggered immunity in rice: Obstacles and perspectives

Improving rice immunity is one of the most effective approaches to reduce yield loss by biotic factors, with the aim of increasing rice production by 2050 amidst limited natural resources. Triggering a fast and strong immune response to pathogens, effector-triggered immunity (ETI) has intrigued scientists to intensively study and utilize the mechanisms for engineering highly resistant plants. The conservation of ETI components and mechanisms across species enables the use of ETI components to generate broad-spectrum resistance in plants. Numerous efforts have been made to introduce new resistance (R) genes, widen the effector recognition spectrum and generate on-demand R genes. Although engineering ETI across plant species is still associated with multiple challenges, previous attempts have provided an enhanced understanding of ETI mechanisms. Here, we provide a survey of recent reports in the engineering of rice R genes. In addition, we suggest a framework for future studies of R gene-effector interactions, including genome-scale investigations in both rice and pathogens, followed by structural studies of R proteins and effectors, and potential strategies to use important ETI components to improve rice immunity.

Assassination tango: an NLR/NLR-ID immune receptors pair of rapeseed co-operates inside the nucleus to activate cell death

Plant immunity largely relies on intracellular nucleotide-binding domain leucine-rich repeat (NLR) immune receptors. Some plant NLRs carry integrated domains (IDs) that mimic authentic pathogen effector targets. We report here the identification of a genetically linked NLR-ID/NLR pair: BnRPR1 and BnRPR2 in Brassica napus. The NLR-ID carries two ID fusions and the mode of action of the pair conforms to the proposed "integrated sensor/decoy" model. The two NLRs interact and the heterocomplex localizes in the plant-cell nucleus and nucleolus. However, the BnRPRs pair does not operate through a negative regulation as it was previously reported for other NLR-IDs. Cell death is induced only upon co-expression of the two proteins and is dependent on the helper genes, EDS1 and NRG1. The nuclear localization of both proteins seems to be essential for cell death activation, while the IDs of BnRPR1 are dispensable for this purpose. In summary, we describe a new pair of NLR-IDs with interesting features in relation to its regulation and the cell death activation.

An NBS-LRR protein in the Rpp1 locus negates the dominance of Rpp1-mediated resistance against Phakopsora pachyrhizi in soybean

The soybean Rpp1 locus confers resistance to Phakopsora pachyrhizi, causal agent of rust, and resistance is usually dominant over susceptibility. However, dominance of Rpp1-mediated resistance is lost when a resistant genotype (Rpp1 or Rpp1b) is crossed with susceptible line TMG06_0011, and the mechanism of this dominant susceptibility (DS) is unknown. Sequencing the Rpp1 region reveals that the TMG06_0011 Rpp1 locus has a single nucleotide-binding site leucine-rich repeat (NBS-LRR) gene (DS-R), whereas resistant PI 594760B (Rpp1b) is similar to PI 200492 (Rpp1) and has three NBS-LRR resistance gene candidates. Evidence that DS-R is the cause of DS was reflected in virus-induced gene silencing of DS-R in Rpp1b/DS-R or Rpp1/DS-R heterozygous plants with resistance partially restored. In heterozygous Rpp1b/DS-R plants, expression of Rpp1b candidate genes was not significantly altered, indicating no effect of DS-R on transcription. Physical interaction of the DS-R protein with candidate Rpp1b resistance proteins was supported by yeast two-hybrid studies and in silico modeling. Thus, we conclude that suppression of resistance most likely does not occur at the transcript level, but instead probably at the protein level, possibly with Rpp1 function inhibited by binding to the DS-R protein. The DS-R gene was found in other soybean lines, with an estimated allele frequency of 6% in a diverse population, and also found in wild soybean (Glycine soja). The identification of a dominant susceptible NBS-LRR gene provides insight into the behavior of NBS-LRR proteins and serves as a reminder to breeders that the dominance of an R gene can be influenced by a susceptibility allele.

Characterization of defense responses against bacterial pathogens in duckweeds lacking EDS1

ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) mediates the induction of defense responses against pathogens in most angiosperms. However, it has recently been shown that a few species have lost EDS1. It is unknown how defense against disease unfolds and evolves in the absence of EDS1. We utilize duckweeds; a collection of aquatic species that lack EDS1, to investigate this question. We established duckweed-Pseudomonas pathosystems and used growth curves and microscopy to characterize pathogen-induced responses. Through comparative genomics and transcriptomics, we show that the copy number of infection-associated genes and the infection-induced transcriptional responses of duckweeds differ from other model species. Pathogen defense in duckweeds has evolved along different trajectories than in other plants, including genomic and transcriptional reprogramming. Specifically, the miAMP1 domain-containing proteins, which are absent in Arabidopsis, showed pathogen responsive upregulation in duckweeds. Despite such divergence between Arabidopsis and duckweed species, we found conservation of upregulation of certain genes and the role of hormones in response to disease. Our work highlights the importance of expanding the pool of model species to study defense responses that have evolved in the plant kingdom independent of EDS1.

A new NLR disease resistance gene Xa47 confers durable and broad-spectrum resistance to bacterial blight in rice

Bacterial blight (BB) induced by Xanthomonas oryzae pv. oryzae (Xoo) is a devastating bacterial disease in rice. The use of disease resistance (R) genes is the most efficient method to control BB. Members of the nucleotide-binding domain and leucine-rich repeat containing protein (NLR) family have significant roles in plant defense. In this study, Xa47, a new bacterial blight R gene encoding a typical NLR, was isolated from G252 rice material, and XA47 was localized in the nucleus and cytoplasm. Among 180 rice materials tested, Xa47 was discovered in certain BB-resistant materials. Compared with the wild-type G252, the knockout mutants of Xa47 was more susceptible to Xoo. By contrast, overexpression of Xa47 in the susceptible rice material JG30 increased BB resistance. The findings indicate that Xa47 positively regulates the Xoo stress response. Consequently, Xa47 may have application potential in the genetic improvement of plant disease resistance. The molecular mechanism of Xa47 regulation merits additional examination.

A blast fungus zinc-finger fold effector binds to a hydrophobic pocket in host Exo70 proteins to modulate immune recognition in rice

Exocytosis plays an important role in plant-microbe interactions, in both pathogenesis and symbiosis. Exo70 proteins are integral components of the exocyst, an octameric complex that mediates tethering of vesicles to membranes in eukaryotes. Although plant Exo70s are known to be targeted by pathogen effectors, the underpinning molecular mechanisms and the impact of this interaction on infection are poorly understood. Here, we show the molecular basis of the association between the effector AVR-Pii of the blast fungus Maganaporthe oryzae and rice Exo70 alleles OsExo70F2 and OsExo70F3, which is sensed by the immune receptor pair Pii via an integrated RIN4/NOI domain. The crystal structure of AVR-Pii in complex with OsExo70F2 reveals that the effector binds to a conserved hydrophobic pocket in Exo70, defining an effector/target binding interface. Structure-guided and random mutagenesis validates the importance of AVR-Pii residues at the Exo70 binding interface to sustain protein association and disease resistance in rice when challenged with fungal strains expressing effector mutants. Furthermore, the structure of AVR-Pii defines a zinc-finger effector fold (ZiF) distinct from the MAX (Magnaporthe Avrs and ToxB-like) fold previously described for a majority of characterized M. oryzae effectors. Our data suggest that blast fungus ZiF effectors bind a conserved Exo70 interface to manipulate plant exocytosis and that these effectors are also baited by plant immune receptors, pointing to new opportunities for engineering disease resistance.

NLR we there yet? Nucleocytoplasmic coordination of NLR-mediated immunity

Plant intracellular nucleotide-binding leucine-rich repeat immune receptors (NLRs) perceive the activity of pathogen-secreted effector molecules that, when undetected, promote colonisation of hosts. Signalling from activated NLRs converges with and potentiates downstream responses from activated pattern recognition receptors (PRRs) that sense microbial signatures at the cell surface. Efficient signalling of both receptor branches relies on the host cell nucleus as an integration point for transcriptional reprogramming, and on the macromolecular transport processes that mediate the communication between cytoplasm and nucleoplasm. Studies on nuclear pore complexes (NPCs), the nucleoporin proteins (NUPs) that compose NPCs, and nuclear transport machinery constituents that control nucleocytoplasmic transport, have revealed that they play important roles in regulating plant immune responses. Here, we discuss the contributions of nucleoporins and nuclear transport receptor (NTR)-mediated signal transduction in plant immunity with an emphasis on NLR immune signalling across the nuclear compartment boundary and within the nucleus. We also highlight and discuss cytoplasmic and nuclear functions of NLRs and their signalling partners and further consider the potential implications of NLR activation and resistosome formation in both cellular compartments for mediating plant pathogen resistance and programmed host cell death.

Insight into the structure and molecular mode of action of plant paired NLR immune receptors

The specific recognition of pathogen effectors by intracellular nucleotide-binding domain and leucine-rich repeat receptors (NLRs) is an important component of plant immunity. NLRs have a conserved modular architecture and can be subdivided according to their signaling domain that is mostly a coiled-coil (CC) or a Toll/Interleukin1 receptor (TIR) domain into CNLs and TNLs. Single NLR proteins are often sufficient for both effector recognition and immune activation. However, sometimes, they act in pairs, where two different NLRs are required for disease resistance. Functional studies have revealed that in these cases one NLR of the pair acts as a sensor (sNLR) and one as a helper (hNLR). The genes corresponding to such resistance protein pairs with one-to-one functional co-dependence are clustered, generally with a head-to-head orientation and shared promoter sequences. sNLRs in such functional NLR pairs have additional, non-canonical and highly diverse domains integrated in their conserved modular architecture, which are thought to act as decoys to trap effectors. Recent structure-function studies on the Arabidopsis thaliana TNL pair RRS1/RPS4 and on the rice CNL pairs RGA4/RGA5 and Pik-1/Pik-2 are unraveling how such protein pairs function together. Focusing on these model NLR pairs and other recent examples, this review highlights the distinctive features of NLR pairs and their various fascinating mode of action in pathogen effector perception. We also discuss how these findings on NLR pairs pave the way toward improved plant disease resistance.

Show me your ID: NLR immune receptors with integrated domains in plants

Nucleotide-binding and leucine-rich repeat receptors (NLRs) are intracellular plant immune receptors that recognize pathogen effectors secreted into the plant cell. Canonical NLRs typically contain three conserved domains including a central nucleotide binding (NB-ARC) domain, C-terminal leucine-rich repeats (LRRs) and an N-terminal domain. A subfamily of plant NLRs contain additional noncanonical domain(s) that have potentially evolved from the integration of the effector targets in the canonical NLR structure. These NLRs with extra domains are thus referred to as NLRs with integrated domains (NLR-IDs). Here, we first summarize our current understanding of NLR-ID activation upon effector binding, focusing on the NLR pairs Pik-1/Pik-2, RGA4/RGA5, and RRS1/RPS4. We speculate on their potential oligomerization into resistosomes as it was recently shown for certain canonical plant NLRs. Furthermore, we discuss how our growing understanding of the mode of action of NLR-ID continuously informs engineering approaches to design new resistance specificities in the context of rapidly evolving pathogens.

The activity of the RGA5 sensor NLR from rice requires binding of its integrated HMA domain to effectors but not HMA domain self-interaction

The rice nucleotide-binding (NB) and leucine-rich repeat (LRR) domain immune receptors (NLRs) RGA4 and RGA5 form a helper NLR/sensor NLR (hNLR/sNLR) pair that specifically recognizes the effectors AVR-Pia and AVR1-CO39 from the blast fungus Magnaporthe oryzae. While RGA4 contains only canonical NLR domains, RGA5 has an additional unconventional heavy metal-associated (HMA) domain integrated after its LRR domain. This RGA5HMA domain binds the effectors and is crucial for their recognition. Investigation of the three-dimensional structure of the AVR1-CO39/RGA5HMA complex by X-ray crystallography identified a candidate surface for effector binding in the HMA domain and showed that the HMA domain self-interacts in the absence of effector through the same surface. Here, we investigated the relevance of this HMA homodimerization for RGA5 function and the role of the RGA5HMA effector-binding and self-interaction surface in effector recognition. By analysing structure-informed point mutations in the RGA5HMA -binding surface in protein interaction studies and in Nicotiana benthamiana cell death assays, we found that HMA self-interaction does not contribute to RGA5 function. However, the effector-binding surface of RGA5HMA identified by X-ray crystallography is crucial for both in vitro and in vivo effector binding as well as effector recognition. These results support the current hypothesis that noncanonical integrated domains of NLRs act primarily as effector traps and deepen our understanding of the sNLRs' function within NLR pairs.

Current understanding of atypical resistance against fungal pathogens in wheat

Pathogens and pests are a major challenge to global food security. Around one hundred different pests and pathogens challenge wheat, one of the most important food crops in the world. Traditional worldwide use of a few key resistance genes in wheat cultivars has necessitated a diversification of the toolbox of resistance genes in wheat varieties over the coming decades to meet the global production demands. Recent advances in gene discovery and functional characterization of genetic resistance mechanisms in wheat reveal great diversity in the types and effectiveness of the underlying resistance genes. This article summarizes the recent developments in the discovery of non-traditional "atypical" resistance genes in wheat against diverse fungal pathogens.

Genome-wide identification and functional exploration of the legume lectin genes in Brassica napus and their roles in Sclerotinia disease resistance

As one of the largest classes of lectins, legume lectins have a variety of desirable features such as antibacterial and insecticidal activities as well as anti-abiotic stress ability. The Sclerotinia disease (SD) caused by the soil-borne fungus Sclerotinia sclerotiorum is a devastating disease affecting most oil crops such as Brassica napus. Here, we identified 130 legume lectin (LegLu) genes in B. napus, which could be phylogenetically classified into seven clusters. The BnLegLu gene family has been significantly expanded since the whole-genome duplication (WGD) or segmental duplication. Gene structure and conserved motif analysis suggested that the BnLegLu genes were well conserved in each cluster. Moreover, relative to those genes only containing the legume lectin domain in cluster VI-VII, the genes in cluster I-V harbored a transmembrane domain and a kinase domain linked to the legume lectin domain in the C terminus. The expression of most BnLegLu genes was relatively low in various tissues. Thirty-five BnLegLu genes were responsive to abiotic stress, and 40 BnLegLu genes were strongly induced by S. sclerotiorum, with a most significant up-regulation of 715-fold, indicating their functional roles in SD resistance. Four BnLegLu genes were located in the candidate regions of genome-wide association analysis (GWAS) results which resulted from a worldwide rapeseed population consisting of 324 accessions associated with SD. Among them, the positive role of BnLegLus-16 in SD resistance was validated by transient expression in tobacco leaves. This study provides important information on BnLegLu genes, particularly about their roles in SD resistance, which may help targeted functional research and genetic improvement in the breeding of B. napus.

A lineage-specific Exo70 is required for receptor kinase-mediated immunity in barley

In the evolution of land plants, the plant immune system has experienced expansion in immune receptor and signaling pathways. Lineage-specific expansions have been observed in diverse gene families that are potentially involved in immunity but lack causal association. Here, we show that Rps8-mediated resistance in barley to the pathogen Puccinia striiformis f. sp. tritici (wheat stripe rust) is conferred by a genetic module: Pur1 and Exo70FX12, which are together necessary and sufficient. Pur1 encodes a leucine-rich repeat receptor kinase and is the ortholog of rice Xa21, and Exo70FX12 belongs to the Poales-specific Exo70FX clade. The Exo70FX clade emerged after the divergence of the Bromeliaceae and Poaceae and comprises from 2 to 75 members in sequenced grasses. These results demonstrate the requirement of a lineage-specific Exo70FX12 in Pur1-mediated immunity and suggest that the Exo70FX clade may have evolved a specialized role in receptor kinase signaling.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

A genetically linked pair of NLR immune receptors shows contrasting patterns of evolution

Throughout their evolution, plant nucleotide-binding leucine-rich-repeat receptors (NLRs) have acquired widely divergent unconventional integrated domains that enhance their ability to detect pathogen effectors. However, the functional dynamics that drive the evolution of NLRs with integrated domains (NLR-IDs) remain poorly understood. Here, we reconstructed the evolutionary history of an NLR locus prone to unconventional domain integration and experimentally tested hypotheses about the evolution of NLR-IDs. We show that the rice (Oryza sativa) NLR Pias recognizes the effector AVR-Pias of the blast fungal pathogen Magnaporthe oryzae. Pias consists of a functionally specialized NLR pair, the helper Pias-1 and the sensor Pias-2, that is allelic to the previously characterized Pia pair of NLRs: the helper RGA4 and the sensor RGA5. Remarkably, Pias-2 carries a C-terminal DUF761 domain at a similar position to the heavy metal-associated (HMA) domain of RGA5. Phylogenomic analysis showed that Pias-2/RGA5 sensor NLRs have undergone recurrent genomic recombination within the genus Oryza, resulting in up to six sequence-divergent domain integrations. Allelic NLRs with divergent functions have been maintained transspecies in different Oryza lineages to detect sequence-divergent pathogen effectors. By contrast, Pias-1 has retained its NLR helper activity throughout evolution and is capable of functioning together with the divergent sensor-NLR RGA5 to respond to AVR-Pia. These results suggest that opposite selective forces have driven the evolution of paired NLRs: highly dynamic domain integration events maintained by balancing selection for sensor NLRs, in sharp contrast to purifying selection and functional conservation of immune signaling for helper NLRs.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant–pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

Identifying plant genes shaping microbiota composition in the barley rhizosphere

A prerequisite to exploiting soil microbes for sustainable crop production is the identification of the plant genes shaping microbiota composition in the rhizosphere, the interface between roots and soil. Here, we use metagenomics information as an external quantitative phenotype to map the host genetic determinants of the rhizosphere microbiota in wild and domesticated genotypes of barley, the fourth most cultivated cereal globally. We identify a small number of loci with a major effect on the composition of rhizosphere communities. One of those, designated the QRMC-3HS, emerges as a major determinant of microbiota composition. We subject soil-grown sibling lines harbouring contrasting alleles at QRMC-3HS and hosting contrasting microbiotas to comparative root RNA-seq profiling. This allows us to identify three primary candidate genes, including a Nucleotide-Binding-Leucine-Rich-Repeat (NLR) gene in a region of structural variation of the barley genome. Our results provide insights into the footprint of crop improvement on the plant’s capacity of shaping rhizosphere microbes.

A prerequisite to exploiting soil microbes for sustainable crop production is the identification of the plant genes shaping microbiota composition in the rhizosphere. Here, the authors report QTLs and the associated candidate genes underlying rhizosphere microbiome composition in barley.

A Comprehensive Study of the WRKY Transcription Factor Family in Strawberry

WRKY transcription factors play critical roles in plant growth and development or stress responses. Using up-to-date genomic data, a total of 64 and 257 WRKY genes have been identified in the diploid woodland strawberry, Fragaria vesca, and the more complex allo-octoploid commercial strawberry, Fragaria × ananassa cv. Camarosa, respectively. The completeness of the new genomes and annotations has enabled us to perform a more detailed evolutionary and functional study of the strawberry WRKY family members, particularly in the case of the cultivated hybrid, in which homoeologous and paralogous FaWRKY genes have been characterized. Analysis of the available expression profiles has revealed that many strawberry WRKY genes show preferential or tissue-specific expression. Furthermore, significant differential expression of several FaWRKY genes has been clearly detected in fruit receptacles and achenes during the ripening process and pathogen challenged, supporting a precise functional role of these strawberry genes in such processes. Further, an extensive analysis of predicted development, stress and hormone-responsive cis-acting elements in the strawberry WRKY family is shown. Our results provide a deeper and more comprehensive knowledge of the WRKY gene family in strawberry.

Evolutionary footprint of plant immunity

There are pieces of evidence from genomic footprints and fossil records indicating that plants have co-evolved with microbes after terrestrialization for more than 407 million years. Therefore, to truly comprehend plant evolution, we need to understand the co-evolutionary process and history between plants and microbes. Recent developments in genomes and transcriptomes of a vast number of plant species as well as microbes have greatly expanded our knowledge of the evolution of the plant immune system. In this review, we summarize recent advances in the co-evolution between plants and microbes with emphasis on the plant side and point out future research needed for understanding plant-microbial co-evolution. Knowledge of the evolution and variation of the plant immune system will better equip us on designing crops with boosted performance in agricultural fields.

Phylogenomic classification and synteny network analyses deciphered the evolutionary landscape of aldo-keto reductase (AKR) gene superfamily in the plant kingdom

Aldo-keto reductase-domain (PF00248) containing proteins (AKRs) are NAD(P)(H)-dependent oxidoreductases of a multigene superfamily that mediate versatile functions in plants ranging from detoxification, metal chelation, potassium ion efflux to specialized metabolism. To uncover the complete repertoire of AKR gene superfamily in plants, a systematic kingdom-wide identification, phylogeny reconstruction, classification and synteny network clustering analyses were performed in this study using 74 diverse plant genomes. Plant AKRs were omnipresent, legitimately classified into 4 groups (based on phylogeny) and 14 subgroups (based on the ≥ 60% of protein sequence identity). Species composition of AKR subgroups highlights their distinct emergence during plant evolution. Loss of AKR subgroups among plants was apparent and that various lineage-, order/family- and species-specific losses were observed. The subgroups IA, IVB and IVF were flourished and diversified well during plant evolution, likely related to the complexity of plant's specialized metabolism and environmental adaptation. About 65% of AKRs were in genomic synteny regions across the plant kingdom and the AKRs relevant to important functions (e.g. vitamin B₆ metabolism) were in profoundly conserved angiosperm-wide synteny communities. This study underscores the evolutionary landscape of plant AKRs and provides a comprehensive resource to facilitate the functional characterization of them.

New recognition specificity in a plant immune receptor by molecular engineering of its integrated domain

Plant nucleotide-binding and leucine-rich repeat domain proteins (NLRs) are immune sensors that recognize pathogen effectors. Here, we show that molecular engineering of the integrated decoy domain (ID) of an NLR can extend its recognition spectrum to a new effector. We relied for this on detailed knowledge on the recognition of the Magnaporthe oryzae effectors AVR-PikD, AVR-Pia, and AVR1-CO39 by, respectively, the rice NLRs Pikp-1 and RGA5. Both receptors detect their effectors through physical binding to their HMA (Heavy Metal-Associated) IDs. By introducing into RGA5_HMA the AVR-PikD binding residues of Pikp-1_HMA, we create a high-affinity binding surface for this effector. RGA5 variants carrying this engineered binding surface perceive the new ligand, AVR-PikD, and still recognize AVR-Pia and AVR1-CO39 in the model plant N. benthamiana. However, they do not confer extended disease resistance specificity against M. oryzae in transgenic rice plants. Altogether, our study provides a proof of concept for the design of new effector recognition specificities in NLRs through molecular engineering of IDs.

A Systematic Phylogenomic Classification of the Multidrug and Toxic Compound Extrusion Transporter Gene Family in Plants

Multidrug and toxic compound extrusion (MATE) transporters comprise a multigene family that mediates multiple functions in plants through the efflux of diverse substrates including organic molecules, specialized metabolites, hormones, and xenobiotics. MATE classification based on genome-wide studies remains ambiguous, likely due to a lack of large-scale phylogenomic studies and/or reference sequence datasets. To resolve this, we established a phylogeny of the plant MATE gene family using a comprehensive kingdom-wide phylogenomic analysis of 74 diverse plant species. We identified more than 4,000 MATEs, which were classified into 14 subgroups based on a systematic bioinformatics pipeline using USEARCH, blast+ and synteny network tools. Our classification was performed using a four-step process, whereby MATEs sharing ≥ 60% protein sequence identity with a ≤ 1E-05 threshold at different sequence lengths (either full-length, ≥ 60% length, or ≥ 150 amino acids) or retaining in the similar synteny blocks were assigned to the same subgroup. In this way, we assigned subgroups to 95.8% of the identified MATEs, which we substantiated using synteny network clustering analysis. The subgroups were clustered under four major phylogenetic groups and named according to their clockwise appearance within each group. We then generated a reference sequence dataset, the usefulness of which was demonstrated in the classification of MATEs in additional species not included in the original analysis. Approximately 74% of the plant MATEs exhibited synteny relationships with angiosperm-wide or lineage-, order/family-, and species-specific conservation. Most subgroups evolved independently, and their distinct evolutionary trends were likely associated with the development of functional novelties or the maintenance of conserved functions. Together with the systematic classification and synteny network profiling analyses, we identified all the major evolutionary events experienced by the MATE gene family in plants. We believe that our findings and the reference dataset provide a valuable resource to guide future functional studies aiming to explore the key roles of MATEs in different aspects of plant physiology. Our classification framework can also be readily extendable to other (super) families.

Genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq

Background

Nucleotide-binding and leucine-rich repeat (NLR) genes have attracted wide attention due to their crucial role in protecting plants from pathogens. SMRT-RenSeq, combining PacBio sequencing after resistance gene enrichment sequencing (RenSeq), is a powerful method for selectively capturing and sequencing full-length NLRs. Haynaldia villosa, a wild grass species with a proven potential for wheat improvement, confers resistance to multiple diseases. So, genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq can facilitate disease resistance genes exploration.

Results

In this study, SMRT-RenSeq was performed to identify the genome-wide NLR complement of H. villosa. In total, 1320 NLRs were annotated in 1169 contigs, including 772 complete NLRs. All the complete NLRs were phylogenetically analyzed and 11 main clades with special characteristics were derived. NLRs could be captured with high efficiency when aligned with cloned R genes, and cluster expansion in some specific gene loci was observed. The physical location of NLRs to individual chromosomes in H. villosa showed a perfect homoeologous relationship with group 1, 2, 3, 5 and 6 of other Triticeae species, however, NLRs physically located on 4VL were largely in silico predicted to be located on the homoeologous group 7. Fifteen types of integrated domains (IDs) were integrated in 52 NLRs, and Kelch and B3 NLR-IDs were found to have expanded in H. villosa, while DUF948, NAM-associated and PRT_C were detected as unique integrated domains implying the new emergence of NLR-IDs after H. villosa diverged from other species.

Conclusion

SMRT-RenSeq is a powerful tool to identify NLR genes from wild species using the baits of the evolutionary related species with reference sequences. The availability of the NLRs from H. villosa provide a valuable library for R gene mining and transfer of disease resistance into wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08334-w.

A truncated CC-NB-ARC gene TaRPP13L1-3D positively regulates powdery mildew resistance in wheat via the RanGAP-WPP complex-mediated nucleocytoplasmic shuttle

A wheat RPP13-like isoform interacting with WPP1 contributes to quantitative and/or basal resistance to powdery mildew (Blumeria graminis f. sp. tritici) by restricting the development of Bgt conidia. Plant disease resistance (R) genes confer an ability to resist infection by pathogens expressing specific avirulence genes. Recognition of Peronospora parasitica 13-like (RPP13-like) genes belong to the nucleotide-binding site and leucine-rich repeat (NBS-LRR) superfamily and play important roles in resistance to various plant diseases. Previously, we detected a TaRPP13-like gene located on chromosome 3D (TaRPP13L1-3D) in the TaSpl1 resided region, which is strongly induced by the cell death phenotype (Zhang et al. 2021). Here, we investigated the expression and functional role of TaRPP13L1-3D in wheat responding to fungal stress. TaRPP13L1-3D encoded a typical NB-ARC structure characterized by Rx-N and P-loop NTPase domains. TaRPP13L1-3D transcripts were strongly upregulated in wheat by powdery mildew (Blumeria graminis f. sp. tritici; Bgt) and stripe rust (Puccinia striiformis f. sp. tritici; Pst) infection although opposing expression patterns were observed in response to wheat-Bgt in incompatible and compatible backgrounds. Overexpression of TaRPP13L1-3D enhanced disease resistance to Bgt, accompanied by upregulation of the defense-related marker genes encoding phytoalexin-deficient4 (PAD4), thaumatin-like protein (TLP) and chitinase 8-like protein (Chi8L), while silencing of TaRPP13L1-3D disrupted the resistance to Bgt infection. Subcellular localization studies showed that TaRPP13L1-3D is located in both the plasma membrane and nucleus, while yeast-two-hybrid (Y2H) assays indicated that TaRPP13L1-3D interacts with WPP domain-containing protein 1 (TaWPP1). This indicates that TaRPP13L1-3D shuttles between the nucleus and cytoplasm membrane via a mechanism that is mediated by the RanGAP-WPP complex in nuclear pores. This insight into TaRPP13L1-3D will be useful in dissecting the mechanism of fungal resistance in wheat, and understanding the interaction between R gene expression and pathogen defense.

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

Functional diversification gave rise to allelic specialization in a rice NLR immune receptor pair

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.

RefPlantNLR is a comprehensive collection of experimentally validated plant disease resistance proteins from the NLR family

Reference datasets are critical in computational biology. They help define canonical biological features and are essential for benchmarking studies. Here, we describe a comprehensive reference dataset of experimentally validated plant nucleotide-binding leucine-rich repeat (NLR) immune receptors. RefPlantNLR consists of 481 NLRs from 31 genera belonging to 11 orders of flowering plants. This reference dataset has several applications. We used RefPlantNLR to determine the canonical features of functionally validated plant NLRs and to benchmark 5 NLR annotation tools. This revealed that although NLR annotation tools tend to retrieve the majority of NLRs, they frequently produce domain architectures that are inconsistent with the RefPlantNLR annotation. Guided by this analysis, we developed a new pipeline, NLRtracker, which extracts and annotates NLRs from protein or transcript files based on the core features found in the RefPlantNLR dataset. The RefPlantNLR dataset should also prove useful for guiding comparative analyses of NLRs across the wide spectrum of plant diversity and identifying understudied taxa. We hope that the RefPlantNLR resource will contribute to moving the field beyond a uniform view of NLR structure and function.

De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes

We report de novo genome assemblies, transcriptomes, annotations, and methylomes for the 26 inbreds that serve as the founders for the maize nested association mapping population. The number of pan-genes in these diverse genomes exceeds 103,000, with approximately a third found across all genotypes. The results demonstrate that the ancient tetraploid character of maize continues to degrade by fractionation to the present day. Excellent contiguity over repeat arrays and complete annotation of centromeres revealed additional variation in major cytological landmarks. We show that combining structural variation with single-nucleotide polymorphisms can improve the power of quantitative mapping studies. We also document variation at the level of DNA methylation and demonstrate that unmethylated regions are enriched for cis-regulatory elements that contribute to phenotypic variation.

The ETS-ETI cycle: evolutionary processes and metapopulation dynamics driving the diversification of pathogen effectors and host immune factors

The natural diversity of pathogen effectors and host immune components represents a snapshot of the underlying evolutionary processes driving the host-pathogen arms race. In plants, this arms race is manifested by an ongoing cycle of disease and resistance driven by pathogenic effectors that promote disease (effector-triggered susceptibility; ETS) and plant resistance proteins that recognize effector activity to trigger immunity (effector-triggered immunity; ETI). Here we discuss how this ongoing ETS-ETI cycle has shaped the natural diversity of both plant resistance proteins and pathogen effectors. We focus on the evolutionary forces that drive the diversification of the molecules that determine the outcome of plant-pathogen interactions and introduce the concept of metapopulation dynamics (i.e., the introduction of genetic variation from conspecific organisms in different populations) as an alternative mechanism that can introduce and maintain diversity in both host and pathogen populations.

Two NLR immune receptors acquired high-affinity binding to a fungal effector through convergent evolution of their integrated domain

A subset of plant NLR immune receptors carry unconventional integrated domains in addition to their canonical domain architecture. One example is rice Pik-1 that comprises an integrated heavy metal-associated (HMA) domain. Here, we reconstructed the evolutionary history of Pik-1 and its NLR partner, Pik-2, and tested hypotheses about adaptive evolution of the HMA domain. Phylogenetic analyses revealed that the HMA domain integrated into Pik-1 before Oryzinae speciation over 15 million years ago and has been under diversifying selection. Ancestral sequence reconstruction coupled with functional studies showed that two Pik-1 allelic variants independently evolved from a weakly binding ancestral state to high-affinity binding of the blast fungus effector AVR-PikD. We conclude that for most of its evolutionary history the Pik-1 HMA domain did not sense AVR-PikD, and that different Pik-1 receptors have recently evolved through distinct biochemical paths to produce similar phenotypic outcomes. These findings highlight the dynamic nature of the evolutionary mechanisms underpinning NLR adaptation to plant pathogens.

Two NLR immune receptors acquired high-affinity binding to a fungal effector through convergent evolution of their integrated domain

A subset of plant NLR immune receptors carry unconventional integrated domains in addition to their canonical domain architecture. One example is rice Pik-1 that comprises an integrated heavy metal-associated (HMA) domain. Here, we reconstructed the evolutionary history of Pik-1 and its NLR partner, Pik-2, and tested hypotheses about adaptive evolution of the HMA domain. Phylogenetic analyses revealed that the HMA domain integrated into Pik-1 before Oryzinae speciation over 15 million years ago and has been under diversifying selection. Ancestral sequence reconstruction coupled with functional studies showed that two Pik-1 allelic variants independently evolved from a weakly binding ancestral state to high-affinity binding of the blast fungus effector AVR-PikD. We conclude that for most of its evolutionary history the Pik-1 HMA domain did not sense AVR-PikD, and that different Pik-1 receptors have recently evolved through distinct biochemical paths to produce similar phenotypic outcomes. These findings highlight the dynamic nature of the evolutionary mechanisms underpinning NLR adaptation to plant pathogens.

Allelic variants of the NLR protein Rpi-chc1 differentially recognize members of the Phytophthora infestans PexRD12/31 effector superfamily through the leucine-rich repeat domain

Phytophthora infestans is a pathogenic oomycete that causes the infamous potato late blight disease. Resistance (R) genes from diverse Solanum species encode intracellular receptors that trigger effective defense responses upon the recognition of cognate RXLR avirulence (Avr) effector proteins. To deploy these R genes in a durable fashion in agriculture, we need to understand the mechanism of effector recognition and the way the pathogen evades recognition. In this study, we cloned 16 allelic variants of the Rpi-chc1 gene from Solanum chacoense and other Solanum species, and identified the cognate P. infestans RXLR effectors. These tools were used to study effector recognition and co-evolution. Functional and non-functional alleles of Rpi-chc1 encode coiled-coil nucleotide-binding leucine-rich repeat (CNL) proteins, being the first described representatives of the CNL16 family. These alleles have distinct patterns of RXLR effector recognition. While Rpi-chc1.1 recognized multiple PexRD12 (Avrchc1.1) proteins, Rpi-chc1.2 recognized multiple PexRD31 (Avrchc1.2) proteins, both belonging to the PexRD12/31 effector superfamily. Domain swaps between Rpi-chc1.1 and Rpi-chc1.2 revealed that overlapping subdomains in the leucine-rich repeat (LRR) domain are responsible for the difference in effector recognition. This study showed that Rpi-chc1.1 and Rpi-chc1.2 evolved to recognize distinct members of the same PexRD12/31 effector family via the LRR domain. The biased distribution of polymorphisms suggests that exchange of LRRs during host-pathogen co-evolution can lead to novel recognition specificities. These insights will guide future strategies to breed durable resistant varieties.

A Comparative Overview of the Intracellular Guardians of Plants and Animals: NLRs in Innate Immunity and Beyond

Nucleotide-binding domain leucine-rich repeat receptors (NLRs) play important roles in the innate immune systems of both plants and animals. Recent breakthroughs in NLR biochemistry and biophysics have revolutionized our understanding of how NLR proteins function in plant immunity. In this review, we summarize the latest findings in plant NLR biology and draw direct comparisons to NLRs of animals. We discuss different mechanisms by which NLRs recognize their ligands in plants and animals. The discovery of plant NLR resistosomes that assemble in a comparable way to animal inflammasomes reinforces the striking similarities between the formation of plant and animal NLR complexes. Furthermore, we discuss the mechanisms by which plant NLRs mediate immune responses and draw comparisons to similar mechanisms identified in animals. Finally, we summarize the current knowledge of the complex genetic architecture formed by NLRs in plants and animals and the roles of NLRs beyond pathogen detection.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

Analysis of intraspecies diversity reveals a subset of highly variable plant immune receptors and predicts their binding sites

The evolution of recognition specificities by the immune system depends on the generation of receptor diversity and on connecting the binding of new antigens with the initiation of downstream signaling. In plant immunity, the innate Nucleotide-Binding Leucine-Rich Repeat (NLR) receptor family enables antigen binding and immune signaling. In this study, we surveyed the NLR complements of 62 ecotypes of Arabidopsis thaliana and 54 lines of Brachypodium distachyon and identified a limited number of NLR subfamilies that show high allelic diversity. We show that the predicted specificity-determining residues cluster on the surfaces of Leucine-Rich Repeat domains, but the locations of the clusters vary among NLR subfamilies. By comparing NLR phylogeny, allelic diversity, and known functions of the Arabidopsis NLRs, we formulate a hypothesis for the emergence of direct and indirect pathogen-sensing receptors and of the autoimmune NLRs. These findings reveal the recurring patterns of evolution of innate immunity and can inform NLR engineering efforts.

The allelic rice immune receptor Pikh confers extended resistance to strains of the blast fungus through a single polymorphism in the effector binding interface

Arms race co-evolution drives rapid adaptive changes in pathogens and in the immune systems of their hosts. Plant intracellular NLR immune receptors detect effectors delivered by pathogens to promote susceptibility, activating an immune response that halts colonization. As a consequence, pathogen effectors evolve to escape immune recognition and are highly variable. In turn, NLR receptors are one of the most diverse protein families in plants, and this variability underpins differential recognition of effector variants. The molecular mechanisms underlying natural variation in effector recognition by NLRs are starting to be elucidated. The rice NLR pair Pik-1/Pik-2 recognizes AVR-Pik effectors from the blast fungus Magnaporthe oryzae, triggering immune responses that limit rice blast infection. Allelic variation in a heavy metal associated (HMA) domain integrated in the receptor Pik-1 confers differential binding to AVR-Pik variants, determining resistance specificity. Previous mechanistic studies uncovered how a Pik allele, Pikm, has extended recognition to effector variants through a specialized HMA/AVR-Pik binding interface. Here, we reveal the mechanistic basis of extended recognition specificity conferred by another Pik allele, Pikh. A single residue in Pikh-HMA increases binding to AVR-Pik variants, leading to an extended effector response in planta. The crystal structure of Pikh-HMA in complex with an AVR-Pik variant confirmed that Pikh and Pikm use a similar molecular mechanism to extend their pathogen recognition profile. This study shows how different NLR receptor alleles functionally converge to extend recognition specificity to pathogen effectors.

Diversity, structure and function of the coiled-coil domains of plant NLR immune receptors

Plant nucleotide-binding, leucine-rich repeat receptors (NLRs) perceive pathogen avirulence effectors and activate defense responses. Nucleotide-binding, leucine-rich repeat receptors are classified into coiled-coil (CC)-containing and Toll/interleukin-1 receptor (TIR)-containing NLRs. Recent advances suggest that NLR CC domains often function in signaling activation, especially for induction of cell death. In this review, we outline our current understanding of NLR CC domains, including their diversity/classification and structure, their roles in cell death induction, disease resistance, and interaction with other proteins. Furthermore, we provide possible directions for future work.

Recognition and defence of plant-infecting fungal pathogens

Attempted infections of plants with fungi result in diverse outcomes ranging from symptom-less resistance to severe disease and even death of infected plants. The deleterious effect on crop yield have led to intense focus on the cellular and molecular mechanisms that explain the difference between resistance and susceptibility. This research has uncovered plant resistance or susceptibility genes that explain either dominant or recessive inheritance of plant resistance with many of them coding for receptors that recognize pathogen invasion. Approaches based on cell biology and phytochemistry have contributed to identifying factors that halt an invading fungal pathogen from further invasion into or between plant cells. Plant chemical defence compounds, antifungal proteins and structural reinforcement of cell walls appear to slow down fungal growth or even prevent fungal penetration in resistant plants. Additionally, the hypersensitive response, in which a few cells undergo a strong local immune reaction, including programmed cell death at the site of infection, stops in particular biotrophic fungi from spreading into surrounding tissue. In this review, we give a general overview of plant recognition and defence of fungal parasites tracing back to the early 20th century with a special focus on Triticeae and on the progress that was made in the last 30 years.

Comparative Genomics and Functional Studies of Wheat BED-NLR Loci

Nucleotide-binding leucine-rich-repeat (LRR) receptors (NLRs) with non-canonical integrated domains (NLR-IDs) are widespread in plant genomes. Zinc-finger BED (named after the Drosophila proteins Boundary Element-Associated Factor and DNA Replication-related Element binding Factor, named BED hereafter) are among the most frequently found IDs. Five BED-NLRs conferring resistance against bacterial and fungal pathogens have been characterized. However, it is unknown whether BED-NLRs function in a manner similar to other NLR-IDs. Here, we used chromosome-level assemblies of wheat to explore the Yr7 and Yr5a genomic regions and show that, unlike known NLR-ID loci, there is no evidence for a NLR-partner in their vicinity. Using neighbor-network analyses, we observed that BED domains from BED-NLRs share more similarities with BED domains from single-BED proteins and from BED-containing proteins harboring domains that are conserved in transposases. We identified a nuclear localization signal (NLS) in Yr7, Yr5, and the other characterized BED-NLRs. We thus propose that this is a feature of BED-NLRs that confer resistance to plant pathogens. We show that the NLS was functional in truncated versions of the Yr7 protein when expressed in N. benthamiana. We did not observe cell-death upon the overexpression of Yr7 full-length, truncated, and 'MHD' variants in N. benthamiana. This suggests that either this system is not suitable to study BED-NLR signaling or that BED-NLRs require additional components to trigger cell death. These results define novel future directions to further understand the role of BED domains in BED-NLR mediated resistance.