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

Divergent molecular pathways govern temperature-dependent wheat stem rust resistance genes

The wheat stem rust pathogen Puccinia graminis f. sp. tritici (Pgt) causes severe crop losses worldwide. Several stem rust resistance (Sr) genes exhibit temperature-dependent immune responses. Sr6-mediated resistance is enhanced at lower temperatures, whereas Sr13 and Sr21 resistances are enhanced at higher temperatures. Here, we clone Sr6 using mutagenesis and resistance gene enrichment and sequencing (MutRenSeq), identifying it to encode a nucleotide-binding leucine-rich repeat (NLR) protein with an integrated BED domain. Sr6 temperature sensitivity is also transferred to wheat plants transformed with the Sr6 gene. Differential gene expression analysis of near-isogenic lines inoculated with Pgt at varying temperatures reveals that genes upregulated in the low-temperature-effective Sr6 response differ from those upregulated in the high-temperature-effective responses associated with Sr13 and Sr21. These findings highlight divergent molecular pathways involved in temperature-sensitive immunity and inform future strategies for deployment and engineering of genetic resistance in response to a changing climate.

The resistance awakens: Diversity at the DNA, RNA, and protein levels informs engineering of plant immune receptors from Arabidopsis to crops

Plants rely on germline-encoded, innate immune receptors to sense pathogens and initiate the defense response. The exponential increase in quality and quantity of genomes, RNA-seq datasets, and protein structures has underscored the incredible biodiversity of plant immunity. Arabidopsis continues to serve as a valuable model and theoretical foundation of our understanding of wild plant diversity of immune receptors, while expansion of study into agricultural crops has also revealed distinct evolutionary trajectories and challenges. Here, we provide the classical context for study of both intracellular nucleotide-binding, leucine-rich repeat receptors and surface-localized pattern recognition receptors at the levels of DNA sequences, transcriptional regulation, and protein structures. We then examine how recent technology has shaped our understanding of immune receptor evolution and informed our ability to efficiently engineer resistance. We summarize current literature and provide an outlook on how researchers take inspiration from natural diversity in bioengineering efforts for disease resistance from Arabidopsis and other model systems to crops.

From Lesions to Lessons: Two Decades of Filamentous Plant Pathogen Genomics

Many filamentous microorganisms, such as fungi and oomycetes, have evolved the ability to colonize plants and cause devastating crop diseases. Coevolutionary conflicts with their hosts have shaped the genomes of these plant pathogens. Over the past 20 years, genomics and genomics-enabled technologies have revealed remarkable diversity in genome size, architecture, and gene regulatory mechanisms. Technical and conceptual advances continue to provide novel insights into evolutionary dynamics, diversification of distinct genomic compartments, and facilitated molecular disease diagnostics. In this review, we discuss how genomics has advanced our understanding of genome organization and plant-pathogen coevolution and provide a perspective on future developments in the field. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Majority of the Highly Variable NLRs in Maize Share Genomic Location and Contain Additional Target-Binding Domains

Nucleotide-binding, leucine-rich repeat (LRR) proteins (NLRs) are a major class of immune receptors in plants. NLRs include both conserved and rapidly evolving members; however, their evolutionary trajectory in crops remains understudied. Availability of crop pan-genomes enables analysis of the recent events in the evolution of this highly complex gene family within domesticated species. Here, we investigated the NLR complement of 26 nested association mapping (NAM) founder lines of maize. We found that maize has just four main subfamilies containing rapidly evolving highly variable NLR (hvNLR) receptors. Curiously, three of these phylogenetically distinct hvNLR lineages are located in adjacent clusters on chromosome 10. Members of the same hvNLR clade show variable expression and methylation across lines and tissues, which is consistent with their rapid evolution. By combining sequence diversity analysis and AlphaFold2 computational structure prediction, we predicted ligand-binding sites in the hvNLRs. We also observed novel insertion domains in the LRR regions of two hvNLR subfamilies that likely contribute to target recognition. To make this analysis accessible, we created NLRCladeFinder, a Google Colaboratory notebook, that accepts any newly identified NLR sequence, places it in the evolutionary context of the maize pan-NLRome, and provides an updated clade alignment, phylogenetic tree, and sequence diversity information for the gene of interest. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Reversible ubiquitination conferred by domain shuffling controls paired NLR immune receptor complex homeostasis in plant immunity

Plant intracellular NLR immune receptors can function individually or in pairs to detect pathogen effectors and activate immune responses. NLR homeostasis has to be tightly regulated to ensure proper defense without triggering autoimmunity. However, in contrast to singleton NLRs, the mechanisms controlling the paired NLRs complex homeostasis are less understood. The paired Arabidopsis RRS1/RPS4 immune receptor complex confers disease resistance through effector recognition mediated by the integrated WRKY domain of RRS1. Here, through proximity labeling, we reveal a ubiquitination-deubiquitination cycle that controls the homeostasis of the RRS1/RPS4 complex. E3 ligase RARE directly binds and ubiquitinates RRS1's WRKY domain to promote its proteasomal degradation, thereby destabilizing RPS4 indirectly and compromising the stability and function of the RRS1/RPS4 complex. Conversely, the deubiquitinating enzymes UBP12/UBP13 deubiquitinate RRS1's WRKY domain, counteracting RARE's effects. Interestingly, the abundance of WRKY transcription factors WRKY70 and WRKY41 is also regulated by RARE and UBP12/UBP13. Phylogenetic analysis suggests this regulation likely transferred from WRKY70/WRKY41 to RRS1 upon WRKY domain integration. Our findings improve our understanding of homeostatic regulation of paired NLR complex and uncover a paradigm whereby domain integration can co-opt preexisting post-translational modification to regulate novel protein functions.

The role of Exo70s in plant defense against pathogens and insect pests and their application for crop breeding

Plant diseases caused by pathogens and pests lead to crop losses, posing a threat to global food security. The secretory pathway is an integral component of plant defense. The exocyst complex regulates the final step of the secretory pathway and is thus essential for secretory defense. In the last decades, several subunits of the exocyst complex have been reported to be involved in plant defense, especially Exo70s. This comprehensive review focuses on the functions of the exocyst Exo70s in plant immunity, particularly in recognizing pathogen and pest signatures. We discussed Exo70's interactions with immune receptors and other immune-related proteins, its symbiotic relationships with microbes, and its role in non-host resistance. Finally, we discussed the future engineering breeding of crops with resistance to pathogens and pests based on our current understanding of Exo70s.

Chromosome-level genome assembly of Pinus massoniana provides insights into conifer adaptive evolution

Pinus massoniana, a conifer of significant economic and ecological value in China, is renowned for its wide adaptability and oleoresin production. We sequenced and assembled the chromosomal-level P. massoniana genome, revealing 80,366 protein-coding genes and significant gene family expansions associated with stress response and plant-pathogen interactions. Long-intron genes, which are predominantly presented in low-copy gene families, are strongly linked to the recent long terminal repeat burst in the Pinus genome. By reanalyzing population transcriptomic data, we identified genetic markers linked to oleoresin synthesis, including those within the CYP450 and TPS gene families. The results suggest that the genes of the resin terpene biosynthesis pathway can be activated in several cell types, and the oleoresin yield may depend on the rate-limiting enzymes. Using a multiomics algorithm, we identified several regulatory factors, including PmMYB4 and PmbZIP2, that interact with TPS and CYP450 genes, potentially playing a role in oleoresin production. This was further validated through molecular genetics analyses. We observed signatures of adaptive evolution in dispersed duplicates and horizontal gene transfer events that have contributed to the species adaptation. This study provides insights for further research into the evolutionary biology of conifers and lays the groundwork for genomic-assisted breeding and sustainable management of Masson pine.

Deciphering the landscape and evolutionary trajectory of NLR immune receptors in Dioscorea alata

Dioscorea alata, a key tuber crop for global food security, is threatened by anthracnose disease caused by Colletotrichum gloeosporioides. However, identification of functional resistance genes against C. gloeosporioides in D. alata is challenging due to low flowering and hybridization efficiency of this plant. Nucleotide-binding leucine-rich repeat (NLR) genes constitute the largest group of plant disease resistance genes, from which functional genes against diverse pathogens across various crops have been cloned. In this study, a comprehensive genome-wide analysis identified 346 NLR genes from D. alata, including one RNL and 345 CNLs. These NLRs were unequally distributed on 20 chromosomes, with chromosome 3 harboring the highest number (78 NLR genes). The majority of NLR genes (91%) were located in multigene clusters, implying that tandem or proximal duplication was the primary driving force for NLR gene expansion in D. alata. Comparative analysis of Dioscoreaceae species revealed high variability and differential expansion patterns of NLR genes. In addition, transcriptome profiling of D. alata post-infection with C. gloeosporioides identified 12 differentially expressed NLR genes. In summary, this study sheds new light on the genetic architecture and evolutionary dynamics of D. alata NLR genes, offering valuable insights for cloning functional genes against C. gloeosporioides.

The pan-NLRome analysis based on 23 genomes reveals the diversity of NLRs in Brassica napus

Brassica napus, a globally significant oilseed crop, exhibits a wide distribution across diverse climatic zones. B. napus is being increasingly susceptible to distinct diseases, such as blackleg, clubroot and sclerotinia stem rot, leading to substantial reductions in yield. Nucleotide-binding site leucine-rich repeat genes (NLRs), the most pivotal family of resistance genes, can be effectively harnessed by identifying and uncovering their diversity to acquire premium disease-resistant gene resources. Here, we collected the genomes of 23 accessions and established the first comprehensive pan-NLRome in B. napus by leveraging multiple genomic resources. We observe significant variation in the number of NLR genes across different B. napus accessions, ranging from 189 to 474. Notably, TNL (TIR-NBS-LRR) genes constitute approximately half of the total count, indicating their predominant presence in B. napus. The number of NLRs in the C subgenome is significantly higher than that in the A subgenome, and chromosome C09 exhibits the highest density of NLR genes with featuring multiple NLR clusters. Domain analysis reveals that the integrated domains significantly enhance the diversity of NLRs, with B3 DNA binding, VQ, and zinc fingers being the most prevalent integrated domains. Pan-genomic analysis reveals that the core type of NLR genes, which is present in most accessions, constitutes approximately 58% of the total NLRs. Furthermore, we conduct a comparative analysis of the diversity of NLR genes across distinct ecotypes, leading to the identification of ecotype-specific NLRs and their integrated domains. In conclusion, our study effectively addresses the limitations of a single reference genome and provides valuable insights into the diversity of NLR genes in B. napus, thereby contributing to disease resistance breeding.

The blast pathogen effector AVR-Pik binds and stabilizes rice heavy metal-associated (HMA) proteins to co-opt their function in immunity

Intracellular nucleotide-binding domain and leucine-rich repeat-containing (NLR) receptors play crucial roles in immunity across multiple domains of life. In plants, a subset of NLRs contain noncanonical integrated domains that are thought to have evolved from host targets of pathogen effectors to serve as pathogen baits. However, the functions of host proteins with similarity to NLR integrated domains and the extent to which they are targeted by pathogen effectors remain largely unknown. Here, we show that the blast fungus effector AVR-Pik binds a subset of related rice proteins containing a heavy metal-associated (HMA) domain, one of the domains that has repeatedly integrated into plant NLR immune receptors. We find that AVR-Pik binding stabilizes the rice small HMA (sHMA) proteins OsHIPP19 and OsHIPP20. Knockout of OsHIPP20 causes enhanced disease resistance towards the blast pathogen, indicating that OsHIPP20 is a susceptibility gene (S-gene). We propose that AVR-Pik has evolved to bind HMA domain proteins and co-opt their function to suppress immunity. Yet this binding carries a trade-off, it triggers immunity in plants carrying NLR receptors with integrated HMA domains.

A single NLR gene confers resistance to leaf and stripe rust in wheat

Nucleotide-binding leucine-rich repeat (NLR) disease resistance genes typically confer resistance against races of a single pathogen. Here, we report that Yr87/Lr85, an NLR gene from Aegilops sharonensis and Aegilops longissima, confers resistance against both P. striiformis tritici (Pst) and Puccinia triticina (Pt) that cause stripe and leaf rust, respectively. Yr87/Lr85 confers resistance against Pst and Pt in wheat introgression as well as transgenic lines. Comparative analysis of Yr87/Lr85 and the cloned Triticeae NLR disease resistance genes shows that Yr87/Lr85 contains two distinct LRR domains and that the gene is only found in Ae. sharonensis and Ae. longissima. Allele mining and phylogenetic analysis indicate multiple events of Yr87/Lr85 gene flow between the two species and presence/absence variation explaining the majority of resistance to wheat leaf rust in both species. The confinement of Yr87/Lr85 to Ae. sharonensis and Ae. longissima and the resistance in wheat against Pst and Pt highlight the potential of these species as valuable sources of disease resistance genes for wheat improvement.

A trade-off between investment in molecular defense repertoires and growth in plants

Given the negative fitness effects that pathogens impose on their hosts, the benefits of resistance should be universal. However, there is marked variation across plant species in the number of nucleotide-binding leucine-rich repeat receptors, which form a cornerstone of defense. The growth-defense trade-off hypothesis predicts costs associated with defense investment to generate variation in these traits. Our analysis comparing features of the intracellular immune-receptor repertoires with trait data of 187 species shows that in wild plants, the size of the molecular defense repertoire correlates negatively with growth. By contrast, we do not find evidence for a growth-defense trade-off in agricultural plants. Our cross-species approach highlights the central role of defense investment in shaping ecological trait variation and its sensitivity to domestication.

Structure-aware annotation of leucine-rich repeat domains

Protein domain annotation is typically done by predictive models such as HMMs trained on sequence motifs. However, sequence-based annotation methods are prone to error, particularly in calling domain boundaries and motifs within them. These methods are limited by a lack of structural information accessible to the model. With the advent of deep learning-based protein structure prediction, existing sequenced-based domain annotation methods can be improved by taking into account the geometry of protein structures. We develop dimensionality reduction methods to annotate repeat units of the Leucine Rich Repeat solenoid domain. The methods are able to correct mistakes made by existing machine learning-based annotation tools and enable the automated detection of hairpin loops and structural anomalies in the solenoid. The methods are applied to 127 predicted structures of LRR-containing intracellular innate immune proteins in the model plant Arabidopsis thaliana and validated against a benchmark dataset of 172 manually-annotated LRR domains.

High-quality chromosome scale genome assemblies of two important Sorghum inbred lines, Tx2783 and RTx436

Sorghum bicolor (L.) Moench is a significant grass crop globally, known for its genetic diversity. High quality genome sequences are needed to capture the diversity. We constructed high-quality, chromosome-level genome assemblies for two vital sorghum inbred lines, Tx2783 and RTx436. Through advanced single-molecule techniques, long-read sequencing and optical maps, we improved average sequence continuity 19-fold and 11-fold higher compared to existing Btx623 v3.0 reference genome and obtained 19 and 18 scaffolds (N50 of 25.6 and 14.4) for Tx2783 and RTx436, respectively. Our gene annotation efforts resulted in 29 612 protein-coding genes for the Tx2783 genome and 29 265 protein-coding genes for the RTx436 genome. Comparative analyses with 26 plant genomes which included 18 sorghum genomes and 8 outgroup species identified around 31 210 protein-coding gene families, with about 13 956 specific to sorghum. Using representative models from gene trees across the 18 sorghum genomes, a total of 72 579 pan-genes were identified, with 14% core, 60% softcore and 26% shell genes. We identified 99 genes in Tx2783 and 107 genes in RTx436 that showed functional enrichment specifically in binding and metabolic processes, as revealed by the GO enrichment Pearson Chi-Square test. We detected 36 potential large inversions in the comparison between the BTx623 Bionano map and the BTx623 v3.1 reference sequence. Strikingly, these inversions were notably absent when comparing Tx2783 or RTx436 with the BTx623 Bionano map. These inversion were mostly in the pericentromeric region which is known to have low complexity regions and harder to assemble and suggests the presence of potential artifacts in the public BTx623 reference assembly. Furthermore, in comparison to Tx2783, RTx436 exhibited 324 883 additional Single Nucleotide Polymorphisms (SNPs) and 16 506 more Insertions/Deletions (INDELs) when using BTx623 as the reference genome. We also characterized approximately 348 nucleotide-binding leucine-rich repeat (NLR) disease resistance genes in the two genomes. These high-quality genomes serve as valuable resources for discovering agronomic traits and structural variation studies.

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 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.