Alien Domains Shaped the Modular Structure of Plant NLR Proteins

The N-terminal domains of NLR immune receptors exhibit structural and functional similarities across divergent plant lineages

Nucleotide-binding domain and leucine-rich repeat (NLR) proteins are a prominent class of intracellular immune receptors in plants. However, our understanding of plant NLR structure and function is limited to the evolutionarily young flowering plant clade. Here, we describe an extended spectrum of NLR diversity across divergent plant lineages and demonstrate the structural and functional similarities of N-terminal domains that trigger immune responses. We show that the broadly distributed coiled-coil (CC) and toll/interleukin-1 receptor (TIR) domain families of nonflowering plants retain immune-related functions through translineage activation of cell death in the angiosperm Nicotiana benthamiana. We further examined a CC subfamily specific to nonflowering lineages and uncovered an essential N-terminal MAEPL motif that is functionally comparable with motifs in resistosome-forming CC-NLRs. Consistent with a conserved role in immunity, the ectopic activation of CCMAEPL in the nonflowering liverwort Marchantia polymorpha led to profound growth inhibition, defense gene activation, and signatures of cell death. Moreover, comparative transcriptomic analyses of CCMAEPL activity delineated a common CC-mediated immune program shared across evolutionarily divergent nonflowering and flowering plants. Collectively, our findings highlight the ancestral nature of NLR-mediated immunity during plant evolution that dates its origin to at least ∼500 million years ago.

A comprehensive review of soybean RNL and TIR domain proteins

Both prokaryotic and eukaryotic organisms use the nucleotide-binding domain/leucine-rich repeat (NBD/LRR)-triggered immunity (NLR-triggered immunity) signaling pathway to defend against pathogens. Plant NLRs are intracellular immune receptors that can bind to effector proteins secreted by pathogens. Dicotyledonous plants express a type of NLR known as TIR domain-containing NLRs (TNLs). TIR domains are enzymes that catalyze the production of small molecules that are essential for immune signaling and lead to plant cell death. The activation of downstream TNL signaling components, such as enhanced disease susceptibility 1 (EDS1), phytoalexin deficient 4 (PAD4), and senescence-associated gene 101 (SAG101), is facilitated by these small molecules. Helper NLRs (hNLRs) and the EDS1-PAD4/SAG101 complex associate after activation, causing the hNLRs to oligomerize, translocate to the plasma membrane (PM), and produce cation-selective channels. According to a recent theory, cations enter cells through pores created by oligomeric hNLRs and trigger cell death. Occasionally, TNLs can self-associate to create higher-order oligomers. Here, we categorized soybean TNLs based on the protein domains that they possess. We believe that TNLs may help soybean plants effectively fight pathogens by acting as a source of genetic resistance. In summary, the purpose of this review is to elucidate the range of TNLs that are expressed in soybean.

NLR immune receptors: structure and function in plant disease resistance

Nucleotide-binding and leucine-rich repeat receptors (NLRs) are a diverse family of intracellular immune receptors that play crucial roles in recognizing and responding to pathogen invasion in plants. This review discusses the overall model of NLR activation and provides an in-depth analysis of the different NLR domains, including N-terminal executioner domains, the nucleotide-binding oligomerization domain (NOD) module, and the leucine-rich repeat (LRR) domain. Understanding the structure-function relationship of these domains is essential for developing effective strategies to improve plant disease resistance and agricultural productivity.

PlantNLRatlas: a comprehensive dataset of full- and partial-length NLR resistance genes across 100 chromosome-level plant genomes

Plants have evolved two layers of protection against biotic stress: PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI). The primary mechanism of ETI involves nucleotide-binding leucine-rich repeat immune receptors (NLRs). Although NLR genes have been studied in several plant species, a comprehensive database of NLRs across a diverse array of species is still lacking. Here, we present a thorough analysis of NLR genes across 100 high-quality plant genomes (PlantNLRatlas). The PlantNLRatlas includes a total of 68,452 NLRs, of which 3,689 are full-length and 64,763 are partial-length NLRs. The majority of NLR groups were phyletically clustered. In addition, the domain sequences were found to be highly conserved within each NLR group. Our PlantNLRatlas dataset is complementary to RefPlantNLR, a collection of NLR genes which have been experimentally confirmed. The PlantNLRatlas should prove helpful for comparative investigations of NLRs across a range of plant groups, including understudied taxa. Finally, the PlantNLRatlas resource is intended to help the field move past a monolithic understanding of NLR structure and function.

NB-LRR Lineage-Specific Equipment Is Sorted Out by Sequence Pattern Adaptation and Domain Segment Shuffling

The nucleotide-binding and leucine-rich repeat (NB-LRR) genes, also known as resistance (R)-genes, play an important role in the activation of immune responses. In recent years, large-scale studies have been performed to highlight the diversification of plant NB-LRR repertories. It is well known that, to provide new functionalities, NB-LRR sequences are subject to duplication, domain fusions and acquisition and other kinds of mutations. Although some mechanisms that govern NB-LRR protein domain adaptations have been uncovered, to retrace the plant-lineage-specific evolution routes of R protein structure, a multi-genome comparative analysis was performed. This study allowed us to define groups of genes sharing homology relationships across different species. It is worth noting that the most populated groups contained well-characterized R proteins. The arsenal profile of such groups was investigated in five botanical families, including important crop species, to underline specific adaptation signatures. In addition, the dissection of 70 NB domains of well-characterized R-genes revealed the NB core motifs from which the three main R protein classes have been diversified. The structural remodeling of domain segments shaped the specific NB-LRR repertoires observed in each plant species. This analysis provided new evolutionary and functional insights on NB protein domain shuffling. Taken together, such findings improved our understanding of the molecular adaptive selection mechanisms occurring at plant R loci.

Multilevel evolution shapes the function of NB-LRR encoding genes in plant innate immunity

A sophisticated innate immune system based on diverse pathogen receptor genes (PRGs) evolved in the history of plant life. To reconstruct the direction and magnitude of evolutionary trajectories of a given gene family, it is critical to detect the ancestral signatures. The rearrangement of functional domains made up the diversification found in PRG repertoires. Structural rearrangement of ancient domains mediated the NB-LRR evolutionary path from an initial set of modular proteins. Events such as domain acquisition, sequence modification and temporary or stable associations are prominent among rapidly evolving innate immune receptors. Over time PRGs are continuously shaped by different forces to find their optimal arrangement along the genome. The immune system is controlled by a robust regulatory system that works at different scales. It is important to understand how the PRG interaction network can be adjusted to meet specific needs. The high plasticity of the innate immune system is based on a sophisticated functional architecture and multi-level control. Due to the complexity of interacting with diverse pathogens, multiple defense lines have been organized into interconnected groups. Genomic architecture, gene expression regulation and functional arrangement of PRGs allow the deployment of an appropriate innate immunity response.

Prediction of NB-LRR resistance genes based on full-length sequence homology

The activation of plant immunity is mediated by resistance (R)-gene receptors, also known as nucleotide-binding leucine-rich repeat (NB-LRR) genes, which in turn trigger the authentic defense response. R-gene identification is a crucial goal for both classic and modern plant breeding strategies for disease resistance. The conventional method identifies NB-LRR genes using a protein motif/domain-based search (PDS) within an automatically predicted gene set of the respective genome assembly. PDS proved to be imprecise since repeat masking prior to automatic genome annotation unwittingly prevented comprehensive NB-LRR gene detection. Furthermore, R-genes have diversified in a species-specific manner, so that NB-LRR gene identification cannot be universally standardized. Here, we present the full-length Homology-based R-gene Prediction (HRP) method for the comprehensive identification and annotation of a genome's R-gene repertoire. Our method has substantially addressed the complex genomic organization of tomato (Solanum lycopersicum) NB-LRR gene loci, proving to be more performant than the well-established RenSeq approach. HRP efficiency was also tested on three differently assembled and annotated Beta sp. genomes. Indeed, HRP identified up to 45% more full-length NB-LRR genes compared to previous approaches. HRP also turned out to be a more refined strategy for R-gene allele mining, testified by the identification of hitherto undiscovered Fom-2 homologs in five Cucurbita sp. genomes. In summary, our high-performance method for full-length NB-LRR gene discovery will propel the identification of novel R-genes towards development of improved cultivars.

Transcriptomic and genomic analysis provides new insights in molecular and genetic processes involved in zucchini ZYMV tolerance

BACKGROUND

Cucurbita pepo is highly susceptible to Zucchini yellow mosaic virus (ZYMV) and the resistance found in several wild species cannot be considered as complete or broad-spectrum resistance. In this study, a source of tolerance introgressed in C. pepo (381e) from C. moschata, in True French (TF) background, was investigated 12 days post-inoculation (DPI) at transcriptomic and genomic levels.

RESULTS

The comparative RNA-sequencing (RNA-Seq) of TF (susceptible to ZYMV) and 381e (tolerant to ZYMV) allowed the evaluation of about 33,000 expressed transcripts and the identification of 146 differentially expressed genes (DEGs) in 381e, mainly involved in photosynthesis, transcription, cytoskeleton organization and callose synthesis. By contrast, the susceptible cultivar TF triggered oxidative processes related to response to biotic stimulus and activated key regulators of plant virus intercellular movement. In addition, the discovery of variants located in transcripts allowed the identification of two chromosome regions rich in Single Nucleotide Polymorphisms (SNPs), putatively introgressed from C. moschata, containing genes exclusively expressed in 381e.

CONCLUSION

381e transcriptome analysis confirmed a global improvement of plant fitness by reducing the virus titer and movement. Furthermore, genes implicated in ZYMV tolerance in C. moschata introgressed regions were detected. Our work provides new insight into the plant virus recovery process and a better understanding of the molecular basis of 381e tolerance.

The APAF1_C/WD40 repeat domain-encoding gene from the sea lettuce Ulva mutabilis sheds light on the evolution of NB-ARC domain-containing proteins in green plants

We advance Ulva's genetic tractability and highlight its value as a model organism by characterizing its APAF1_C/WD40 domain-encoding gene, which belongs to a family that bears homology to R genes. The multicellular chlorophyte alga Ulva mutabilis (Ulvophyceae, Ulvales) is native to coastal ecosystems worldwide and attracts both high socio-economic and scientific interest. To further understand the genetic mechanisms that guide its biology, we present a protocol, based on adapter ligation-mediated PCR, for retrieving flanking sequences in U. mutabilis vector-insertion mutants. In the created insertional library, we identified a null mutant with an insertion in an apoptotic protease activating factor 1 helical domain (APAF1_C)/WD40 repeat domain-encoding gene. Protein domain architecture analysis combined with phylogenetic analysis revealed that this gene is a member of a subfamily that arose early in the evolution of green plants (Viridiplantae) through the acquisition of a gene that also encoded N-terminal nucleotide-binding adaptor shared by APAF-1, certain R-gene products and CED-4 (NB-ARC) and winged helix-like (WH-like) DNA-binding domains. Although phenotypic analysis revealed no mutant phenotype, gene expression levels in control plants correlated to the presence of bacterial symbionts, which U. mutabilis requires for proper morphogenesis. In addition, our analysis led to the discovery of a putative Ulva nucleotide-binding site and leucine-rich repeat (NBS-LRR) Resistance protein (R-protein), and we discuss how the emergence of these R proteins in green plants may be linked to the evolution of the APAF1_C/WD40 protein subfamily.

NLRexpress-A bundle of machine learning motif predictors-Reveals motif stability underlying plant Nod-like receptors diversity

Examination of a collection of over 80,000 Plant Nod-like receptors (NLRs) revealed an overwhelming sequence diversity underlying functional specificity of pathogen detection, signaling and cooperativity. The NLR canonical building blocks-CC/TIR/RPW8, NBS and LRR-contain, however, a number of conserved sequence motifs showing a significant degree of invariance amongst different NLR groups. To identify these motifs we developed NLRexpress-a bundle of 17 machine learning (ML)-based predictors, able to swiftly and precisely detect CC, TIR, NBS, and LRR motifs while minimizing computing time without accuracy losses-aimed as an instrument scalable for screening overall proteomes, transcriptomes or genomes for identifying integral NLRs and discriminating them against incomplete sequences lacking key motifs. These predictors were further used to screen a subset of ∼34,000 regular plant NLR sequences. Motifs were analyzed using unsupervised ML techniques to assess the structural correlations hidden underneath pattern variabilities. Both the NB-ARC switch domain which admittedly is the most conserved region of NLRs and the highly diverse LRR domain with its vastly variable lengths and repeat irregularities-show well-defined relations between motif subclasses, highlighting the importance of structural invariance in shaping NLR sequence diversity. The online NLRexpress webserver can be accessed at https://nlrexpress.biochim.ro.

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.

Large-scale gene gains and losses molded the NLR defense arsenal during the Cucurbita evolution

Genome-wide annotation reveals that the gene birth-death process of the Cucurbita R family is associated with a species-specific diversification of TNL and CNL protein classes. The Cucurbitaceae family includes nearly 1000 plant species known universally as cucurbits. Cucurbita genus includes many economically important worldwide crops vulnerable to more than 200 pathogens. Therefore, the identification of pathogen-recognition genes is of utmost importance for this genus. The major class of plant-resistance (R) genes encodes nucleotide-binding site and leucine-rich repeat (NLR) proteins, and is divided into three sub-classes namely, TIR-NB-LRR (TNL), CC-NB-LRR (CNL) and RPW8-NB-LRR (RNL). Although the characterization of the NLR gene family has been carried out in important Cucurbita species, this information is still linked to the availability of sequenced genomes. In this study, we analyzed 40 de novo transcriptomes and 5 genome assemblies, which were explored to investigate the Cucurbita expressed-NLR (eNLR) and NLR repertoires using an ad hoc gene annotation approach. Over 1850 NLR-encoding genes were identified, finely characterized and compared to 96 well-characterized plant R-genes. The maximum likelihood analyses revealed an unusual diversification of CNL/TNL genes and a strong RNL conservation. Indeed, several gene gain and loss events have shaped the Cucurbita NLR family. Finally, to provide a first validation step Cucurbita, eNLRs were explored by real-time PCR analysis. The NLR repertories of the 12 Cucurbita species presented in this paper will be useful to discover novel R-genes.

Polymorphisms in plants to restrict losses to pathogens: From gene family expansions to complex network evolution

Genetic polymorphisms are the basis of the natural diversity seen in all life on earth, also in plant-pathogen interactions. Initially, studies on plant-pathogen interaction focused on reporting phenotypic variation in resistance properties and on the identification of underlying major genes. Nowadays, the field of plant-pathogen interactions is moving from focusing on families of single dominant genes involved in gene-for-gene interactions to an understanding of the plant immune system in the context of a much more complex signaling network and quantitative resistance. Simultaneously, studies on pathosystems from the wild and genome analyses advanced, revealing tremendous variation in natural plant populations. It is now imperative to place studies on genetic diversity and evolution of plant-pathogen interactions in the appropriate molecular biological, as well as evolutionary, context.

Genotype-Specific Expression and NLR Repertoire Contribute to Phenotypic Resistance Diversity in Plantago lanceolata

High levels of phenotypic variation in resistance appears to be nearly ubiquitous across natural host populations. Molecular processes contributing to this variation in nature are still poorly known, although theory predicts resistance to evolve at specific loci driven by pathogen-imposed selection. Nucleotide-binding leucine-rich repeat (NLR) genes play an important role in pathogen recognition, downstream defense responses and defense signaling. Identifying the natural variation in NLRs has the potential to increase our understanding of how NLR diversity is generated and maintained, and how to manage disease resistance. Here, we sequenced the transcriptomes of five different Plantago lanceolata genotypes when inoculated by the same strain of obligate fungal pathogen Podosphaera plantaginis. A de novo transcriptome assembly of RNA-sequencing data yielded 24,332 gene models with N50 value of 1,329 base pairs and gene space completeness of 66.5%. The gene expression data showed highly varying responses where each plant genotype demonstrated a unique expression profile in response to the pathogen, regardless of the resistance phenotype. Analysis on the conserved NB-ARC domain demonstrated a diverse NLR repertoire in P. lanceolata consistent with the high phenotypic resistance diversity in this species. We find evidence of selection generating diversity at some of the NLR loci. Jointly, our results demonstrate that phenotypic resistance diversity results from a crosstalk between different defense mechanisms. In conclusion, characterizing the architecture of resistance in natural host populations may shed unprecedented light on the potential of evolution to generate variation.

The Tomato Interspecific NB-LRR Gene Arsenal and Its Impact on Breeding Strategies

Tomato (Solanum lycopersicum L.) is a model system for studying the molecular basis of resistance in plants. The investigation of evolutionary dynamics of tomato resistance (R)-loci provides unique opportunities for identifying factors that promote or constrain genome evolution. Nucleotide-binding domain and leucine-rich repeat (NB-LRR) receptors belong to one of the most plastic and diversified families. The vast amount of genomic data available for Solanaceae and wild tomato relatives provides unprecedented insights into the patterns and mechanisms of evolution of NB-LRR genes. Comparative analysis remarked a reshuffling of R-islands on chromosomes and a high degree of adaptive diversification in key R-loci induced by species-specific pathogen pressure. Unveiling NB-LRR natural variation in tomato and in other Solanaceae species offers the opportunity to effectively exploit genetic diversity in genomic-driven breeding programs with the aim of identifying and introducing new resistances in tomato cultivars. Within this motivating context, we reviewed the repertoire of NB-LRR genes available for tomato improvement with a special focus on signatures of adaptive processes. This issue is still relevant and not thoroughly investigated. We believe that the discovery of mechanisms involved in the generation of a gene with new resistance functions will bring great benefits to future breeding strategies.

Inferring RPW8-NLRs's evolution patterns in seed plants: case study in Vitis vinifera

Genomic and transcriptomic studies in plants and, more in deep, in grapevine reveal that the disease-resistance RNL gene family is highly variable. RNLs (RPW8-NLRs) are a phylogenetically distinct class of nucleotide oligomerization domain (NOD)-like receptors (NLRs) identified in plants. Two RNLs, namely, the NRG1 (N Requirement Gene 1) and the ADR1 (Activated Disease Resistance 1), have been characterized; however, little is known about the RNL evolutionary history in higher plants. To trace the diversification of RNL gene subfamily, we scanned the NLR proteins of 73 plant genomes belonging to 29 taxa, revealing a noticeable diversification across species and within the same genus or botanic family together with a conspicuous expansion in important crop species. To explore the RNL variability in Vitis vinifera and gain information with respect to their structure, evolutionary diversification of five grape genomes ('Aglianico', 'Falanghina', 'Sultanina', 'Tannat', and 'Nebbiolo') has been compared to the reference genome ('Pinot Noir'). The number of RNLs ranged from 6 ('Sultanina') to 14 ('Nebbiolo'), in contrast to the 10 'Pinot Noir' RNLs. The phylogenetic study on grapevine RNLs revealed that all collapsed into NRG1-clade, rather than four. To investigate more in depth the means of intraspecific variability of grape RNL copies, a transcriptomic profiling in response to powdery mildew (PM) infection was carried out through qRT-PCRs and public databases interrogation. The RNL expression variability identified in transcriptome data sets supports the hypothesis of a functional expansion/contraction in grapevine varieties. Although no direct correlations between grapevine PM-resistance and RNL expression was identified, our work can provide good candidates for functional studies able to elucidate the putative "helper" role of RNLs in grape immune signalling.