Tracing the origin and evolutionary history of plant nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes

Copy number variation and population-specific immune genes in the model vertebrate zebrafish

Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations, we identified a total of 1513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.

Heterotrophic bacteria trigger transcriptome remodelling in the photosynthetic picoeukaryote Micromonas commoda

Marine biogeochemical cycles are built on interactions between surface ocean microbes, particularly those connecting phytoplankton primary producers to heterotrophic bacteria. Details of these associations are not well understood, especially in the case of direct influences of bacteria on phytoplankton physiology. Here we catalogue how the presence of three marine bacteria (Ruegeria pomeroyi DSS-3, Stenotrophomonas sp. SKA14 and Polaribacter dokdonensis MED152) individually and uniquely impact gene expression of the picoeukaryotic alga Micromonas commoda RCC 299. We find a dramatic transcriptomic remodelling by M. commoda after 8 h in co-culture, followed by an increase in cell numbers by 56 h compared with the axenic cultures. Some aspects of the algal transcriptomic response are conserved across all three bacterial co-cultures, including an unexpected reduction in relative expression of photosynthesis and carbon fixation pathways. Expression differences restricted to a single bacterium are also observed, with the Flavobacteriia P. dokdonensis uniquely eliciting changes in relative expression of algal genes involved in biotin biosynthesis and the acquisition and assimilation of nitrogen. This study reveals that M. commoda has rapid and extensive responses to heterotrophic bacteria in ways that are generalizable, as well as in a taxon specific manner, with implications for the diversity of phytoplankton-bacteria interactions ongoing in the surface ocean.

Two adjacent NLR genes conferring quantitative resistance to clubroot disease in Arabidopsis are regulated by a stably inherited epiallelic variation

Clubroot caused by the protist Plasmodiophora brassicae is a major disease affecting cultivated Brassicaceae. Using a combination of quantitative trait locus (QTL) fine mapping, CRISPR-Cas9 validation, and extensive analyses of DNA sequence and methylation patterns, we revealed that the two adjacent neighboring NLR (nucleotide-binding and leucine-rich repeat) genes AT5G47260 and AT5G47280 cooperate in controlling broad-spectrum quantitative partial resistance to the root pathogen P. brassicae in Arabidopsis and that they are epigenetically regulated. The variation in DNA methylation is not associated with any nucleotide variation or any transposable element presence/absence variants and is stably inherited. Variations in DNA methylation at the Pb-At5.2 QTL are widespread across Arabidopsis accessions and correlate negatively with variations in expression of the two genes. Our study demonstrates that natural, stable, and transgenerationally inherited epigenetic variations can play an important role in shaping resistance to plant pathogens by modulating the expression of immune receptors.

Review: Plant microRNAs in pathogen defense: A panacea or a piece of the puzzle?

Plant microRNAs (miRNAs) control key agronomic traits that are associated with their conserved role(s) in development. However, despite a multitude of studies, the utility of miRNAs in plant-pathogen resistance remains less certain. Reviewing the literature identifies three general classes of miRNAs regarding plant pathogen defense. Firstly, a number of evolutionary dynamic 22 nucleotide miRNA families that repress large numbers of plant immunity genes, either directly, or through triggering the biogenesis of secondary siRNAs. However, understanding of their role in defense and of their manipulation to enhance pathogen resistance are still lacking. Secondly, highly conserved miRNAs that indirectly impact disease resistance through their targets that are primarily regulating development or hormone signaling. Any alteration of these miRNAs usually results in pleiotropic impacts, which may alter disease resistance in some plant species, and against some pathogens. Thirdly, are the comparatively diverse and evolutionary dynamic set of non-conserved miRNAs, some of which contribute to pathogen resistance, but whose narrow evolutionary presence will likely restrict their utility. Therefore, reflecting the diverse and evolving nature of plant-pathogen interactions, a complex interplay of plant miRNAs with pathogen responses exists. Any miRNA-based solution for pathogen resistance will likely be highly specific, rather than a general panacea.

Genomic evolution and insights into agronomic trait innovations of Sesamum species

Sesame is an ancient oilseed crop with high oil content and quality. However, the evolutionary history and genetic mechanisms of its valuable agronomic traits remain unclear. Here, we report chromosome-scale genomes of cultivated sesame (Sesamum indicum L.) and six wild Sesamum species, representing all three karyotypes within this genus. Karyotyping and genome-based phylogenic analysis revealed the evolutionary route of Sesamum species from n = 13 to n = 16 and revealed that allotetraploidization occurred in the wild species Sesamum radiatum. Early divergence of the Sesamum genus (48.5-19.7 million years ago) during the Tertiary period and its ancient phylogenic position within eudicots were observed. Pan-genome analysis revealed 9164 core gene families in the 7 Sesamum species. These families are significantly enriched in various metabolic pathways, including fatty acid (FA) metabolism and FA biosynthesis. Structural variations in SiPT1 and SiDT1 within the phosphatidyl ethanolamine-binding protein gene family lead to the genomic evolution of plant-architecture and inflorescence-development phenotypes in Sesamum. A genome-wide association study (GWAS) of an interspecific population and genome comparisons revealed a long terminal repeat insertion and a sequence deletion in DIR genes of wild Sesamum angustifolium and cultivated sesame, respectively; both variations independently cause high susceptibility to Fusarium wilt disease. A GWAS of 560 sesame accessions combined with an overexpression study confirmed that the NAC1 and PPO genes play an important role in upregulating oil content of sesame. Our study provides high-quality genomic resources for cultivated and wild Sesamum species and insights that can improve molecular breeding strategies for sesame and other oilseed crops.

Plant NLR immunity activation and execution: a biochemical perspective

Plants deploy cell-surface and intracellular receptors to detect pathogen attack and trigger innate immune responses. Inside host cells, families of nucleotide-binding/leucine-rich repeat (NLR) proteins serve as pathogen sensors or downstream mediators of immune defence outputs and cell death, which prevent disease. Established genetic underpinnings of NLR-mediated immunity revealed various strategies plants adopt to combat rapidly evolving microbial pathogens. The molecular mechanisms of NLR activation and signal transmission to components controlling immunity execution were less clear. Here, we review recent protein structural and biochemical insights to plant NLR sensor and signalling functions. When put together, the data show how different NLR families, whether sensors or signal transducers, converge on nucleotide-based second messengers and cellular calcium to confer immunity. Although pathogen-activated NLRs in plants engage plant-specific machineries to promote defence, comparisons with mammalian NLR immune receptor counterparts highlight some shared working principles for NLR immunity across kingdoms.

Regulation of plant immunity via small RNA-mediated control of NLR expression

Plants use different receptors to detect potential pathogens: membrane-anchored pattern recognition receptors (PRRs) activated upon perception of pathogen-associated molecular patterns (PAMPs) that elicit pattern-triggered immunity (PTI); and intracellular nucleotide-binding leucine-rich repeat proteins (NLRs) activated by detection of pathogen-derived effectors, activating effector-triggered immunity (ETI). The interconnections between PTI and ETI responses have been increasingly reported. Elevated NLR levels may cause autoimmunity, with symptoms ranging from fitness cost to developmental arrest, sometimes combined with run-away cell death, making accurate control of NLR dosage key for plant survival. Small RNA-mediated gene regulation has emerged as a major mechanism of control of NLR dosage. Twenty-two nucleotide miRNAs with the unique ability to trigger secondary siRNA production from target transcripts are particularly prevalent in NLR regulation. They enhance repression of the primary NLR target, but also bring about repression of NLRs only complementary to secondary siRNAs. We summarize current knowledge on miRNAs and siRNAs in the regulation of NLR expression with an emphasis on 22 nt miRNAs and propose that miRNA and siRNA regulation of NLR levels provides additional links between PTI and NLR defense pathways to increase plant responsiveness against a broad spectrum of pathogens and control an efficient deployment of defenses.

Identification, characterization, and validation of NBS-encoding genes in grass pea

Grass pea is a promising crop with the potential to provide food and fodder, but its genomics has not been adequately explored. Identifying genes for desirable traits, such as drought tolerance and disease resistance, is critical for improving the plant. Grass pea currently lacks known R-genes, including the nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family, which plays a key role in protecting the plant from biotic and abiotic stresses. In our study, we used the recently published grass pea genome and available transcriptomic data to identify 274 NBS-LRR genes. The evolutionary relationships between the classified genes on the reported plants and LsNBS revealed that 124 genes have TNL domains, while 150 genes have CNL domains. All genes contained exons, ranging from 1 to 7. Ten conserved motifs with lengths ranging from 16 to 30 amino acids were identified. We found TIR-domain-containing genes in 132 LsNBSs, with 63 TIR-1 and 69 TIR-2, and RX-CCLike in 84 LsNBSs. We also identified several popular motifs, including P-loop, Uup, kinase-GTPase, ABC, ChvD, CDC6, Rnase_H, Smc, CDC48, and SpoVK. According to the gene enrichment analysis, the identified genes undergo several biological processes such as plant defense, innate immunity, hydrolase activity, and DNA binding. In the upstream regions, 103 transcription factors were identified that govern the transcription of nearby genes affecting the plant excretion of salicylic acid, methyl jasmonate, ethylene, and abscisic acid. According to RNA-Seq expression analysis, 85% of the encoded genes have high expression levels. Nine LsNBS genes were selected for qPCR under salt stress conditions. The majority of the genes showed upregulation at 50 and 200 μM NaCl. However, LsNBS-D18, LsNBS-D204, and LsNBS-D180 showed reduced or drastic downregulation compared to their respective expression levels, providing further insights into the potential functions of LsNBSs under salt stress conditions. They provide valuable insights into the potential functions of LsNBSs under salt stress conditions. Our findings also shed light on the evolution and classification of NBS-LRR genes in legumes, highlighting the potential of grass pea. Further research could focus on the functional analysis of these genes, and their potential use in breeding programs to improve the salinity, drought, and disease resistance of this important crop.

The evolution of plant NLR immune receptors and downstream signal components

Along with the emergence of green plants on this planet one billion years ago, the nucleotide binding site leucine-rich repeat (NLR) gene family originated and diverged into at least three subclasses. Two of them, with either characterized N-terminal toll/interleukin-1 receptor (TIR) or coiled-coil (CC) domain, serve as major types of immune receptor of effector-triggered immunity (ETI) in plants, whereas the one having a N-terminal Resistance to powdery mildew8 (RPW8) domain, functions as signal transfer component to them. In this review, we briefly summarized the history of identification of diverse NLR subclasses across Viridiplantae lineages during the establishment of NLR category, and highlighted recent advances on the evolution of NLR genes and several key downstream signal components under the background of ecological adaption.

Comprehensive Analysis and Functional Verification of the Pinus massoniana NBS-LRR Gene Family Involved in the Resistance to Bursaphelenchus xylophilus

Pinus massoniana Lamb. is a crucial timber and resin conifer in China, but its plantation industry is threatened by outbreaks of pine wilt disease (PWD) caused by Bursaphelenchus xylophilus (pinewood nematode; PWN). However, as of yet, there is no comprehensive analysis of NBS-LRR genes in P. massoniana involved in its defense against PWN. In this study, 507 NBS genes were identified in the transcriptome of resistant and susceptible P. masoniana inoculated with the PWN. The phylogenetic analysis and expression profiles of resistant and susceptible P. massoniana revealed that the up-regulated PmNBS-LRR97 gene was involved in conferring resistance to PWN. The results of real-time quantitative PCR (qRT-PCR) showed that PmNBS-LRR97 was significantly up-regulated after PWN infection, especially in the stems. Subcellular localization indicated that PmNBS-LRR97 located to the cell membrane. PmNBS-LRR97 significantly activated the expression of reactive oxygen species (ROS)-related genes in P. massoniana. In addition, the overexpression of PmNBS-LRR97 was capable of promoting the production of ROS, aiding in plant growth and development. In summary, PmNBS-LRR97 participates in the defense response to PWN and plays an active role in conferring resistance in P. massoniana. This finding provides new insight into the regulatory mechanism of the R gene in P. massoniana.

Genome-wide identification and analysis of the NLR gene family in Medicago ruthenica

Medicago ruthenica, important forage in the legume family, possesses high nutritional value and carries abundant tolerance genes. This study used whole-genome data of M. ruthenica to perform a genome-wide analysis of the nucleotide-binding site-leucine-rich repeat receptor (NLR) gene family, which is the largest family of plant disease resistance genes (R genes). A total of 338 NLR genes were identified in the M. ruthenica genome, including 160 typical genes that contained 80 coiled-coil (CC)-NBS-LRR (CNL) genes, 76 toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes, four resistance to powdery mildew 8 (RPW8)-NBS-LRR (RNL) subclass genes, and 178 atypical NLR genes encoding proteins without at least one important domain. Among its eight chromosomes, M. ruthenica chromosomes 3 and 8 contained most of the NLR genes. More than 40% of all NLR genes were located on these two chromosomes, mainly in multigene clusters. The NLR proteins of M. ruthenica had six highly conserved motifs: P-loop, GLPL, RNBS-D, kinase-2, RNBS-C, and MHDV. Phylogenetic analysis revealed that the NLR genes of M. ruthenica formed three deeply separated clades according to the N-terminal domain of the proteins encoded by these genes. Gene duplication and syntenic analysis suggested four gene duplication types in the NLR genes of M. ruthenica, namely, tandem, proximal, dispersed, and segmental duplicates, which involved 189, 49, 59, and 41 genes, respectively. A total of 41 segmental duplication genes formed 23 NLR gene pairs located on syntenic chromosomal blocks mainly between chromosomes 6 and 7. In addition, syntenic analysis between M. truncatula and M. ruthenica revealed 193 gene pairs located on syntenic chromosomal blocks of the two species. The expression analysis of M. ruthenica NLR genes showed that 303 (89.6%) of the NLR genes were expressed in different varieties. Overall, this study described the full NLR profile of the M. ruthenica genome to provide an important resource for mining disease-resistant genes and disease-resistant breeding.

Predicting Cloned Disease Resistance Gene Homologs (CDRHs) in Radish, Underutilised Oilseeds, and Wild Brassicaceae Species

Brassicaceae crops, including Brassica, Camelina and Raphanus species, are among the most economically important crops globally; however, their production is affected by several diseases. To predict cloned disease resistance (R) gene homologs (CDRHs), we used the protein sequences of 49 cloned R genes against fungal and bacterial diseases in Brassicaceae species. In this study, using 20 Brassicaceae genomes (17 wild and 3 domesticated species), 3172 resistance gene analogs (RGAs) (2062 nucleotide binding-site leucine-rich repeats (NLRs), 497 receptor-like protein kinases (RLKs) and 613 receptor-like proteins (RLPs)) were identified. CDRH clusters were also observed in Arabis alpina, Camelina sativa and Cardamine hirsuta with assigned chromosomes, consisting of 62 homogeneous (38 NLR, 17 RLK and 7 RLP clusters) and 10 heterogeneous RGA clusters. This study highlights the prevalence of CDRHs in the wild relatives of the Brassicaceae family, which may lay the foundation for rapid identification of functional genes and genomics-assisted breeding to develop improved disease-resistant Brassicaceae crop cultivars.

An Overview of PRR- and NLR-Mediated Immunities: Conserved Signaling Components across the Plant Kingdom That Communicate Both Pathways

Cell-surface-localized pattern recognition receptors (PRRs) and intracellular nucleotide-binding domain and leucine-rich repeat receptors (NLRs) are plant immune proteins that trigger an orchestrated downstream signaling in response to molecules of microbial origin or host plant origin. Historically, PRRs have been associated with pattern-triggered immunity (PTI), whereas NLRs have been involved with effector-triggered immunity (ETI). However, recent studies reveal that such binary distinction is far from being applicable to the real world. Although the perception of plant pathogens and the final mounting response are achieved by different means, central hubs involved in signaling are shared between PTI and ETI, blurring the zig-zag model of plant immunity. In this review, we not only summarize our current understanding of PRR- and NLR-mediated immunities in plants, but also highlight those signaling components that are evolutionarily conserved across the plant kingdom. Altogether, we attempt to offer an overview of how plants mediate and integrate the induction of the defense responses that comprise PTI and ETI, emphasizing the need for more evolutionary molecular plant-microbe interactions (EvoMPMI) studies that will pave the way to a better understanding of the emergence of the core molecular machinery involved in the so-called evolutionary arms race between plants and microbes.

Innate Immunity Mechanisms in Marine Multicellular Organisms

The innate immune system provides an adequate response to stress factors and pathogens through pattern recognition receptors (PRRs), located on the surface of cell membranes and in the cytoplasm. Generally, the structures of PRRs are formed by several domains that are evolutionarily conserved, with a fairly high degree of homology in representatives of different species. The orthologs of TLRs, NLRs, RLRs and CLRs are widely represented, not only in marine chordates, but also in invertebrates. Study of the interactions of the most ancient marine multicellular organisms with microorganisms gives us an idea of the evolution of molecular mechanisms of protection against pathogens and reveals new functions of already known proteins in ensuring the body’s homeostasis. The review discusses innate immunity mechanisms of protection of marine invertebrate organisms against infections, using the examples of ancient multicellular hydroids, tunicates, echinoderms, and marine worms in the context of searching for analogies with vertebrate innate immunity. Due to the fact that mucous membranes first arose in marine invertebrates that have existed for several hundred million years, study of their innate immune system is both of fundamental importance in terms of understanding molecular mechanisms of host defense, and of practical application, including the search of new antimicrobial agents for subsequent use in medicine, veterinary and biotechnology.

The Rm1 and Rm2 Resistance Genes to Green Peach Aphid (Myzus persicae) Encode the Same TNL Proteins in Peach (Prunus persica L.)

The green peach aphid (GPA), Myzus persicae, is an important pest of the peach crop. Three major dominant resistance genes have already been detected, Rm1 in the Weeping Flower Peach (WFP) clone, Rm2 in the Rubira clone, and Rm3 in the Fen Shouxing clone. In this study, after NGS resequencing of WFP and Rubira, we found that their genomic sequences in the Rm1 and Rm2 region were similar but very different from that of the susceptible reference peach Lovell. We constructed a BAC library for the GPA-resistant WFP and screened four BAC clones to sequence the target region. The new sequence was 61.7 Kb longer than Lovell and was annotated with four different TIR_NBS_LRR genes. Among them, the TNL1 gene was very overexpressed in WFP leaves 24 h after GPA infestation. This gene was also present and expressed in the Rubira clone and had the same sequence as the candidate Rm3 gene, supporting the hypothesis that the three genes share the same origin. In addition, we identified a second TNL, TNL2, located at 35.4 Kb from TNL1 and slightly overexpressed after GPA infestation. Kasp and size molecular markers were designed for use in marker-assisted selection and were validated in a peach segregating population.

Reference genome assemblies reveal the origin and evolution of allohexaploid oat

Common oat (Avena sativa) is an important cereal crop serving as a valuable source of forage and human food. Although reference genomes of many important crops have been generated, such work in oat has lagged behind, primarily owing to its large, repeat-rich polyploid genome. Here, using Oxford Nanopore ultralong sequencing and Hi-C technologies, we have generated a reference-quality genome assembly of hulless common oat, comprising 21 pseudomolecules with a total length of 10.76 Gb and contig N50 of 75.27 Mb. We also produced genome assemblies for diploid and tetraploid Avena ancestors, which enabled the identification of oat subgenomes and provided insights into oat chromosomal evolution. The origin of hexaploid oat is inferred from whole-genome sequencing, chloroplast genomes and transcriptome assemblies of different Avena species. These findings and the high-quality reference genomes presented here will facilitate the full use of crop genetic resources to accelerate oat improvement.

A reference-quality genome assembly of hexaploid oat variety ‘Sanfensan’ and genome assemblies of its diploid and tetraploid Avena ancestors provide insights into the evolutionary history of allohexaploid oat.

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.

Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)

The most predominant type of resistance (R) genes contain nucleotide-binding sites and leucine-rich repeat (NBS-LRR) domains, characterization of which is helpful for plant resistance improvement. However, the NBS genes of Ipomoea trifida (H.B.K.) remain insufficient to date. In this study, a genome-wide analysis of the NBS-encoding gene in I. trifida (H.B.K.) was carried out. A total of 442 NBS encoding genes were identified, amounting to 1.37% of the total genes of I. trifida (H.B.K.). Based on the analysis of the domains, the identified ItfNBS genes were further classified into seven groups: CNL, NL, CN, N, TNL, TN, and RNL. Phylogenetic analysis showed that the I. trifida NBS genes clustered into three independent clades: RNL, TNL, and CNL. Chromosome location analysis revealed that the distribution of ItfNBS genes in chromosomes was uneven, with a number ranging from 3 to 45. Multiple stress-related regulatory elements were detected in the promoters of the NBS-encoding genes, and their expression profiles were obtained. The qRT-PCR analysis revealed that IbNBS10, IbNBS20, IbNBS258, and IbNBS88 responded to stem nematode infection. These results provide critical proof for further characterization and analysis of NBS-encoding genes with important functions.

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.

Structural basis of NLR activation and innate immune signalling in plants

Animals and plants have NLRs (nucleotide-binding leucine-rich repeat receptors) that recognize the presence of pathogens and initiate innate immune responses. In plants, there are three types of NLRs distinguished by their N-terminal domain: the CC (coiled-coil) domain NLRs, the TIR (Toll/interleukin-1 receptor) domain NLRs and the RPW8 (resistance to powdery mildew 8)-like coiled-coil domain NLRs. CC-NLRs (CNLs) and TIR-NLRs (TNLs) generally act as sensors of effectors secreted by pathogens, while RPW8-NLRs (RNLs) signal downstream of many sensor NLRs and are called helper NLRs. Recent studies have revealed three dimensional structures of a CNL (ZAR1) including its inactive, intermediate and active oligomeric state, as well as TNLs (RPP1 and ROQ1) in their active oligomeric states. Furthermore, accumulating evidence suggests that members of the family of lipase-like EDS1 (enhanced disease susceptibility 1) proteins, which are uniquely found in seed plants, play a key role in providing a link between sensor NLRs and helper NLRs during innate immune responses. Here, we summarize the implications of the plant NLR structures that provide insights into distinct mechanisms of action by the different sensor NLRs and discuss plant NLR-mediated innate immune signalling pathways involving the EDS1 family proteins and RNLs.

Triplication is the main evolutionary driving force of NLP transcription factor family in Chinese cabbage and related species

The NODULE-INCEPTION-like protein (NLP) is a plant-specific transcription factor (TF) family that plays an important role in both signal transduction and nitrate assimilation. However, the NLP gene family in Chinese cabbage (Brassica rapa) has yet to be studied. Here we identified 17, 16, and 32 NLP genes in Chinese cabbage, Brassica oleracea, and Brassica napus, respectively. We found that duplication of those NLP genes almost always originated from genome-wide duplication events. Further analysis (using Arabidopsis as a reference) revealed that the NLP family in Chinese cabbage and B. oleracea was characterized by direct expansion caused by whole-genome duplication. By contrast, indirect expansion characterized B. napus, which arose from hybridization and fusion of the two species. In addition, phylogenetic and homology analyses showed that the Brassica NLP gene family has been highly conserved in evolution. Finally, we also identified optimal codons for four studied species. Altogether, through comparative genome analysis methods, we presented compelling evidence that triplication is the main driving force for the NLP TF family's evolution in Chinese cabbage and related Brassica plants, a process evidently highly conserved. This work will help in better understanding the impact of genome-wide duplication on gene families of plants.

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

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

NB-LRR genes: characteristics in three Solanum species and transcriptional response to Ralstonia solanacearum in tomato

NB-LRR genes in the three Solanum species showed specific constitution characteristics and evolved multiple clusters and duplicates. Some genes could respond to biotic stresses such as tomato bacterial wilt. Nucleotide-binding and leucine-rich repeat (NB-LRR, NLR) is a largest resistance gene family in plants, which plays a key role in response to biotic stresses. In this study, NB-LRR genes in cultivated tomato Solanum lycopersicum (Sl) and its wild relatives S. pennellii (Spe) and S. pimpinellifolium (Spi) were analyzed using bioinformatics approaches. In total, 238, 202 and 217 NB-LRR genes of 8 different types were found in Sl, Spe and Spi, respectively. The three species showed similar genomic characteristics. The NB-LRR genes were mainly distributed on chromosomes 4, 5 and 11 and located at the distal zones, forming multiple clusters and tandem duplicates. A large number of homologs appeared through gene expansion, with most Ka/Ks values being less than 1, indicating that purifying selection had occurred in evolution. These genes were mainly expressed in root and could respond to different biotic stresses. RT-qPCR analysis revealed that SlNLR genes could respond to tomato bacterial wilt, with SlNLR1 probably involved in the resistance response, whereas others being the opposite. The transcription factors (TFs) and interaction proteins that regulate target genes were mainly Dof, NAC and MYB families and kinases. The results provide a basis for the isolation and application of related genes in plant disease resistance breeding.

Expression of an Antiviral Gene GmRUN1 from Soybean Is Regulated via Intron-Mediated Enhancement (IME)

Most of R (resistance) genes encode the protein containing NBS-LRR (nucleotide binding site and leucine-rich repeat) domains. Here, N. benthamiana plants were used for transient expression assays at 3-4 weeks of age. We identified a TNL (TIR-NBS-LRR) encoding gene GmRUN1 that was resistant to both soybean mosaic virus (SMV) and tobacco mosaic virus (TMV). Truncation analysis indicated the importance of all three canonical domains for GmRUN1-mediated antiviral activity. Promoter-GUS analysis showed that GmRUN1 expression is inducible by both salicylic acid (SA) and a transcription factor GmDREB3 via the cis-elements as-1 and ERE (ethylene response element), which are present in its promoter region. Interestingly, GmRUN1 gDNA (genomic DNA) shows higher viral resistance than its cDNA (complementary DNA), indicating the existence of intron-mediated enhancement (IME) for GmRUN1 regulation. We provided evidence that intron2 of GmRUN1 increased the mRNA level of native gene GmRUN1, a soybean antiviral gene SRC7 and also a reporter gene Luciferase, indicating the general transcriptional enhancement of intron2 in different genes. In summary, we identified an antiviral TNL type soybean gene GmRUN1, expression of which was regulated at different layers. The investigation of GmRUN1 gene regulatory network would help to explore the mechanism underlying soybean-SMV interactions.

Plant NLR diversity: the known unknowns of pan-NLRomes

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

Identification and characterization of melon circular RNAs involved in powdery mildew responses through comparative transcriptome analysis

Circular RNAs (circRNAs) are a class of newly discovered non-coding RNAs that are typically derived from a genome's exonic, intronic, and intergenic regions. Recent studies of circRNAs in animals and plants have shown that circRNAs are vital in response to various abiotic and biotic stresses. Powdery mildew disease (PM) is a serious fungal disease threatening the melon industry. We performed whole transcriptome sequencing using the leaves of a PM-resistant (M1) and a PM-susceptible (B29) melon to identify circRNAs and determine their molecular functions. A total of 303 circRNAs were identified and >50% circRNAs were derived from exonic regions. Expression levels were significantly altered in 17 and 23 circRNAs after PM infections in B29 and M1, respectively. Melon circRNAs may participate in the response to biotic stimuli, oxidation reduction, metabolic processes, and the regulation of gene expression based on the functional annotation of circRNA parental genes. Furthermore, 27 circRNAs were predicted to be potential targets or 'sponges' for 18 microRNAs (miRNAs). Our results are the first to identify and characterize circRNA functions in melon and may contribute to a better understanding of the role and regulatory mechanisms of circRNAs in resisting PM.

Molecular evidence of the avocado defense response to Fusarium kuroshium infection: a deep transcriptome analysis using RNA-Seq

Fusarium kuroshium is a novel member of the Ambrosia Fusarium Clade (AFC) that has been recognized as one of the symbionts of the invasive Kuroshio shot hole borer, an Asian ambrosia beetle. This complex is considered the causal agent of Fusarium dieback, a disease that has severely threatened natural forests, landscape trees, and avocado orchards in the last 8 years. Despite the interest in this species, the molecular responses of both the host and F. kuroshium during the infection process and disease establishment remain unknown. In this work, we established an in vitro pathosystem using Hass avocado stems inoculated with F. kuroshium to investigate differential gene expression at 1, 4, 7 and 14 days post-inoculation. RNA-seq technology allowed us to obtain data from both the plant and the fungus, and the sequences obtained from both organisms were analyzed independently. The pathosystem established was able to mimic Fusarium dieback symptoms, such as carbohydrate exudation, necrosis, and vascular tissue discoloration. The results provide interesting evidence regarding the genes that may play roles in the avocado defense response to Fusarium dieback disease. The avocado data set comprised a coding sequence collection of 51,379 UniGenes, from which 2,403 (4.67%) were identified as differentially expressed. The global expression analysis showed that F. kuroshium responsive UniGenes can be clustered into six groups according to their expression profiles. The biologically relevant functional categories that were identified included photosynthesis as well as responses to stress, hormones, abscisic acid, and water deprivation. Additionally, processes such as oxidation-reduction, organization and biogenesis of the cell wall and polysaccharide metabolism were detected. Moreover, we identified orthologues of nucleotide-binding leucine-rich receptors, and their possible action mode was analyzed. In F. kuroshium, we identified 57 differentially expressed genes. Interestingly, the alcohol metabolic process biological category had the highest number of upregulated genes, and the enzyme group in this category may play an important role in the mechanisms of secondary metabolite detoxification. Hydrolytic enzymes, such as endoglucanases and a pectate lyase, were also identified, as well as some proteases. In conclusion, our research was conducted mainly to explain how the vascular tissue of a recognized host of the ambrosia complex responds during F. kuroshium infection since Fusarium dieback is an ambrosia beetle-vectored disease and many variables facilitate its establishment.

Fine Mapping of the Ph-2 Gene Conferring Resistance to Late Blight (Phytophthora infestans) in Tomato

Late blight is a devastating tomato disease. Breeding new varieties with multiple resistance (R) genes is highly effective for preventing late blight. The Ph-2 gene mediates resistance to Phytophthora infestans race T₁ in tomato. In this study, we used an F₂ population derived from a cross between Solanum lycopersicum Moboline (resistant) and LA3988 (susceptible) cultivars for the fine mapping of Ph-2. Two flanking markers, CAPS-1 and CC-Ase, mapped Ph-2 to a 141-kb genomic region containing 21 projected genes, 5 of which were identified as putative R genes. The Solyc10g085460 coding sequence varied significantly between the parents. The markers developed and candidate genes identified in this study shall be useful for the molecular breeding of tomato exhibiting increased late blight resistance and for the cloning of the Ph-2 gene.

Cross Kingdom Immunity: The Role of Immune Receptors and Downstream Signaling in Animal and Plant Cell Death

Both plants and animals are endowed with sophisticated innate immune systems to combat microbial attack. In these multicellular eukaryotes, innate immunity implies the presence of cell surface receptors and intracellular receptors able to detect danger signal referred as damage-associated molecular patterns (DAMPs) and pathogen-associated molecular patterns (PAMPs). Membrane-associated pattern recognition receptors (PRRs), such as Toll-like receptors (TLRs), C-type lectin receptors (CLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are employed by these organisms for sensing different invasion patterns before triggering antimicrobial defenses that can be associated with a form of regulated cell death. Intracellularly, animals nucleotide-binding and oligomerization domain (NOD)-like receptors or plants nucleotide-binding domain (NBD)-containing leucine rich repeats (NLRs) immune receptors likely detect effectors injected into the host cell by the pathogen to hijack the immune signaling cascade. Interestingly, during the co-evolution between the hosts and their invaders, key cross-kingdom cell death-signaling macromolecular NLR-complexes have been selected, such as the inflammasome in mammals and the recently discovered resistosome in plants. In both cases, a regulated cell death located at the site of infection constitutes a very effective mean for blocking the pathogen spread and protecting the whole organism from invasion. This review aims to describe the immune mechanisms in animals and plants, mainly focusing on cell death signaling pathways, in order to highlight recent advances that could be used on one side or the other to identify the missing signaling elements between the perception of the invasion pattern by immune receptors, the induction of defenses or the transmission of danger signals to other cells. Although knowledge of plant immunity is less advanced, these organisms have certain advantages allowing easier identification of signaling events, regulators and executors of cell death, which could then be exploited directly for crop protection purposes or by analogy for medical research.

Integrity of the Post-LRR Domain Is Required for TIR-NB-LRR Function

Plants trigger appropriate defense responses, notably, through intracellular nucleotide-binding (NB) and leucine-rich repeat (LRR)-containing receptors (NLRs) that detect secreted pathogen effector proteins. In NLR resistance genes, the toll/interleukin-1 receptor (TIR)-NB-LRR proteins (TNLs) are an important subfamily, out of which approximately half the members carry a post-LRR (PL) domain of unknown role. We first investigated the requirement of the PL domain for TNL-mediated immune response by mutating the most conserved amino acids across PL domains of Arabidopsis thaliana TNLs. We identified several amino acids in the PL domain of RPS4, required for its ability to trigger a hypersensitive response to AvrRps4 in a Nicotiana tabacum transient assay. Mutating the corresponding amino acids within the PL domain of the tobacco TNL gene N also affected its function. Consequently, our results indicate that the integrity of the PL domain at conserved positions is crucial for at least two unrelated TNLs. We then tested the PL domain specificity for function by swapping PL domains between the paralogs RPS4 and RPS4B. Our results suggest that the PL domain is involved in their TNL pair specificity, 'off state' stability, and NLR complex activation. Considering genetically paired Arabidopsis TNLs, we finally compared the PL and TIR domains of their sensor and executor sequences, respectively. While TIR and PL domains from executors present complete motifs, sensors showed a lack of conservation with degenerated motifs. We here provide a contribution to the functional analysis of the PL domain in order to decipher its role for TNL function.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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

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

Genetics of autoimmunity in plants: an evolutionary genetics perspective

Autoimmunity in plants has been found in numerous hybrids as a form of hybrid necrosis and mutant panels. Uncontrolled cell death is a main cellular outcome of autoimmunity, which negatively impacts growth. Its occurrence highlights the vulnerable nature of the plant immune system. Genetic investigation of autoimmunity in hybrid plants revealed that extreme variation in the immune receptor repertoire is a major contributor, reflecting an evolutionary conundrum that plants face in nature. In this review, we discuss natural variation in the plant immune system and its contribution to fitness. The value of autoimmunity genetics lies in its ability to identify combinations of a natural immune receptor and its partner that are predisposed to triggering autoimmunity. The network of immune components for autoimmunity becomes instrumental in revealing mechanistic details of how immune receptors recognize cellular invasion and activate signaling. The list of autoimmunity-risk variants also allows us to infer evolutionary processes contributing to their maintenance in the natural population. Our approach to autoimmunity, which integrates mechanistic understanding and evolutionary genetics, has the potential to serve as a prognosis tool to optimize immunity in crops.

Genome-wide identification and characterization of NBS-encoding genes in Raphanus sativus L. and their roles related to Fusarium oxysporum resistance

BACKGROUND

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes are important for plant development and disease resistance. Although genome-wide studies of NBS-encoding genes have been performed in several species, the evolution, structure, expression, and function of these genes remain unknown in radish (Raphanus sativus L.). A recently released draft R. sativus L. reference genome has facilitated the genome-wide identification and characterization of NBS-encoding genes in radish.

RESULTS

A total of 225 NBS-encoding genes were identified in the radish genome based on the essential NB-ARC domain through HMM search and Pfam database, with 202 mapped onto nine chromosomes and the remaining 23 localized on different scaffolds. According to a gene structure analysis, we identified 99 NBS-LRR-type genes and 126 partial NBS-encoding genes. Additionally, 80 and 19 genes respectively encoded an N-terminal Toll/interleukin-like domain and a coiled-coil domain. Furthermore, 72% of the 202 NBS-encoding genes were grouped in 48 clusters distributed in 24 crucifer blocks on chromosomes. The U block on chromosomes R02, R04, and R08 had the most NBS-encoding genes (48), followed by the R (24), D (23), E (23), and F (17) blocks. These clusters were mostly homogeneous, containing NBS-encoding genes derived from a recent common ancestor. Tandem (15 events) and segmental (20 events) duplications were revealed in the NBS family. Comparative evolutionary analyses of orthologous genes among Arabidopsis thaliana, Brassica rapa, and Brassica oleracea reflected the importance of the NBS-LRR gene family during evolution. Moreover, examinations of cis-elements identified 70 major elements involved in responses to methyl jasmonate, abscisic acid, auxin, and salicylic acid. According to RNA-seq expression analyses, 75 NBS-encoding genes contributed to the resistance of radish to Fusarium wilt. A quantitative real-time PCR analysis revealed that RsTNL03 (Rs093020) and RsTNL09 (Rs042580) expression positively regulates radish resistance to Fusarium oxysporum, in contrast to the negative regulatory role for RsTNL06 (Rs053740).

CONCLUSIONS

The NBS-encoding gene structures, tandem and segmental duplications, synteny, and expression profiles in radish were elucidated for the first time and compared with those of other Brassicaceae family members (A. thaliana, B. oleracea, and B. rapa) to clarify the evolution of the NBS gene family. These results may be useful for functionally characterizing NBS-encoding genes in radish.

Evolution of NLR Resistance Genes in Magnoliids: Dramatic Expansions of CNLs and Multiple Losses of TNLs

Magnoliids are the third-largest group of angiosperms and occupy a critical position in angiosperm evolution. In the past years, due to the lack of sequenced genomes, the disease resistance gene (R gene) profile of magnoliids remains poorly understood. By the genome-wide identification of 1,832 NLR genes from seven magnoliid genomes, we built a framework for the evolution of magnoliid R genes. TNL genes were completely absent from five magnoliids, presumably due to immune pathway deficiencies. A total of 74 ancestral R genes (70 CNLs, 3 TNLs, and 1 RNL) were recovered in a common ancestor of magnoliids, from which all current NLR gene repertoires were derived. Tandem duplication served as the major drive for NLR genes expansion in seven magnoliid genomes, as most surveyed angiosperms. Due to recent rapid expansions, most magnoliids exhibited "a first expansion followed by a slight contraction and a further stronger expansion" evolutionary pattern, while both Litsea cubeba and Persea americana showed a two-times-repeated pattern of "expansion followed by contraction." The transcriptome analysis of seven different tissues of Saururus chinensis revealed a low expression of most NLR genes, with some R genes displaying a relatively higher expression in roots and fruits. Overall, our study sheds light on the evolution of NLR genes in magnoliids, compensates for insufficiency in major angiosperm lineages, and provides an important reference for a better understanding of angiosperm NLR genes.

Advancement of research on plant NLRs evolution, biochemical activity, structural association, and engineering

In this review, we have included evolution of plant intracellular immune receptors, oligomeric complex formation, enzymatic action, engineering, and mechanisms of immune inspection for appropriate defense outcomes. NLR (Nucleotide binding oligomerization domain containing leucine-rich repeat) proteins are the intracellular immune receptors that recognize pathogen-derived virulence factors to confer effector-triggered immunity (ETI). Activation of plant defense by the NLRs are often conveyed through N-terminal Toll-like/ IL-1 receptor (TIR) or non-TIR (coiled-coils or CC) domains. Homodimerization or self-association property of CC/ TIR domains of plant NLRs contribute to their auto-activity and induction of in planta ectopic cell death. High resolution crystal structures of Arabidopsis thaliana RPS4^TIR, L6^TIR, SNC1^TIR, RPP1^TIR and Muscadinia rotundifolia RPV1^TIR showed that interaction is mediated through one or two distinct interfaces i.e., αA and αE helices comprise AE interface and αD and αE helices were found to form DE interface. By contrast, conserved helical regions were determined for CC domains of plant NLRs. Evolutionary history of NLRs diversification has shown that paired forms were originated from NLR singletons. Plant TIRs executed NAD⁺ hydrolysis activity for cell death promotion. Plant NLRs were found to form large oligomeric complexes as observed in animal inflammasomes. We have also discussed different protein engineering methods includes domain shuffling, and decoy modification that increase effector recognition spectrum of plant NLRs. In summary, our review highlights structural basis of perception of the virulence factors by NLRs or NLR pairs to design novel classes of plant immune receptors.

Distinct Evolutionary Patterns of NBS-Encoding Genes in Three Soapberry Family (Sapindaceae) Species

Nucleotide-binding site (NBS)-type disease resistance genes (R genes) play key roles in plant immune responses and have co-evolved with pathogens over the course of plant lifecycles. Comparative genomic studies tracing the dynamic evolution of NBS-encoding genes have been conducted using many important plant lineages. However, studies on Sapindaceae species have not been performed. In this study, a discrepant number of NBS-encoding genes were identified in the genomes of Xanthoceras sorbifolium (180), Dinnocarpus longan (568), and Acer yangbiense (252). These genes were unevenly distributed and usually clustered as tandem arrays on chromosomes, with few existed as singletons. The phylogenetic analysis revealed that NBS-encoding genes formed three monophyletic clades, RPW8-NBS-LRR (RNL), TIR-NBS-LRR (TNL), and CC-NBS-LRR (CNL), which were distinguished by amino acid motifs. The NBS-encoding genes of the X. sorbifolium, D. longan, and A. yangbiense genomes were derived from 181 ancestral genes (three RNL, 23 TNL, and 155 CNL), which exhibited dynamic and distinct evolutionary patterns due to independent gene duplication/loss events. Specifically, X. sorbifolium exhibited a “first expansion and then contraction” evolutionary pattern, while A. yangbiense and D. longan exhibited a “first expansion followed by contraction and further expansion” evolutionary pattern. However, further expansion in D. longan was stronger than in A. yangbiense after divergence, suggesting that D. longan gained more genes in response to various pathogens. Additionally, the ancient and recent expansion of CNL genes generated the dominance of this subclass in terms of gene numbers, while the low copy number status of RNL genes was attributed to their conserved functions.

Taro Genome Assembly and Linkage Map Reveal QTLs for Resistance to Taro Leaf Blight

Taro (Colocasia esculenta) is a food staple widely cultivated in the humid tropics of Asia, Africa, Pacific and the Caribbean. One of the greatest threats to taro production is Taro Leaf Blight caused by the oomycete pathogen Phytophthora colocasiae Here we describe a de novo taro genome assembly and use it to analyze sequence data from a Taro Leaf Blight resistant mapping population. The genome was assembled from linked-read sequences (10x Genomics; ∼60x coverage) and gap-filled and scaffolded with contigs assembled from Oxford Nanopore Technology long-reads and linkage map results. The haploid assembly was 2.45 Gb total, with a maximum contig length of 38 Mb and scaffold N50 of 317,420 bp. A comparison of family-level (Araceae) genome features reveals the repeat content of taro to be 82%, >3.5x greater than in great duckweed (Spirodela polyrhiza), 23%. Both genomes recovered a similar percent of Benchmarking Universal Single-copy Orthologs, 80% and 84%, based on a 3,236 gene database for monocot plants. A greater number of nucleotide-binding leucine-rich repeat disease resistance genes were present in genomes of taro than the duckweed, ∼391 vs. ∼70 (∼182 and ∼46 complete). The mapping population data revealed 16 major linkage groups with 520 markers, and 10 quantitative trait loci (QTL) significantly associated with Taro Leaf Blight disease resistance. The genome sequence of taro enhances our understanding of resistance to TLB, and provides markers that may accelerate breeding programs. This genome project may provide a template for developing genomic resources in other understudied plant species.

Distinctive Growth and Transcriptional Changes of the Diatom Seminavis robusta in Response to Quorum Sensing Related Compounds

In aquatic habitats, diatoms are frequently found in association with Proteobacteria, many members of which employ cell-to-cell communication via N-acyl homoserine lactones (AHLs). It has been suggested that diatoms could distinguish between beneficial and algicidal bacteria in their surroundings by sensing AHLs. Although some microalgae can interfere with AHL signaling, e.g., by releasing AHL mimics or degrading them, molecular responses to AHLs in microalgae are still unclear. Therefore, we tested the effects of short-chained AHLs, i.e., N-hexanoyl homoserine lactone (C6-HSL), N-3-hydroxyhexanoyl homoserine lactone (OH-C6-HSL), and N-3-oxohexanoyl homoserine lactone (oxo-C6-HSL) and long-chained AHLs, i.e., N-tetradecanoyl homoserine lactone (C14-HSL), N-3-hydroxytetradecanoyl homoserine lactone (OH-C14-HSL), and N-3-oxotetradecanoyl homoserine lactone (oxo-C14-HSL), on growth of the benthic diatom Seminavis robusta. All tested short-chained AHLs did not affect diatom growth, while long-chained AHLs promoted (C14-HSL) or inhibited (OH-C14-HSL and oxo-C14-HSL) growth. To investigate the physiological effects of these long-chained AHLs in more detail, an RNA-seq experiment was performed during which S. robusta was treated with the growth-promoting C14-HSL and the growth-inhibiting oxo-C14-HSL. One tetramic acid was also tested (TA14), a structural rearrangement product of oxo-C14-HSL, which also induced growth inhibition in S. robusta. After 3 days of treatment, analysis revealed that 3,410 genes were differentially expressed in response to at least one of the compounds. In the treatment with the growth-promoting C14-HSL many genes involved in intracellular signaling were upregulated. On the other hand, exposure to growth-inhibiting oxo-C14-HSL and TA14 triggered a switch in lipid metabolism towards increased fatty acid degradation. In addition, oxo-C14-HSL led to downregulation of cell cycle genes, which is in agreement with the stagnation of cell growth in this treatment. Combined, our results indicate that bacterial signaling molecules with high structural similarity induce contrasting physiological responses in S. robusta.

Enzymatic Functions for Toll/Interleukin-1 Receptor Domain Proteins in the Plant Immune System

Rationally engineered improvements to crop plants will be needed to keep pace with increasing demands placed on agricultural systems by population growth and climate change. Engineering of plant immune systems provides an opportunity to increase yields by limiting losses to pathogens. Intracellular immune receptors are commonly used as agricultural disease resistance traits. Despite their importance, how intracellular immune receptors confer disease resistance is still unknown. One major class of immune receptors in dicots contains a Toll/Interleukin-1 Receptor (TIR) domain. The mechanisms of TIR-containing proteins during plant immunity have remained elusive. The TIR domain is an ancient module found in archaeal, bacterial and eukaryotic proteins. In animals, TIR domains serve a structural role by generating innate immune signaling complexes. The unusual animal TIR-protein, SARM1, was recently discovered to function instead as an enzyme that depletes cellular NAD⁺ (nicotinamide adenine dinucleotide) to trigger axonal cell death. Two recent reports have found that plant TIR proteins also have the ability to cleave NAD⁺. This presents a new paradigm from which to consider how plant TIR immune receptors function. Here, we will review recent reports of the structure and function of TIR-domain containing proteins. Intriguingly, it appears that TIR proteins in all kingdoms may use similar enzymatic mechanisms in a variety of cell death and immune pathways. We will also discuss TIR structure–function hypotheses in light of the recent publication of the ZAR1 resistosome structure. Finally, we will explore the evolutionary context of plant TIR-containing proteins and their downstream signaling components across phylogenies and the functional implications of these findings.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

•

Species-wide NLR diversity is high but not unlimited

•

A large fraction of NLR diversity is recovered with 40–50 accessions

•

Presence/absence variation in NLRs is widespread, resulting in a mosaic population

•

A high diversity of NLR-integrated domains favor known virulence targets

Alien Domains Shaped the Modular Structure of Plant NLR Proteins

Plant innate immunity mostly relies on nucleotide-binding (NB) and leucine-rich repeat (LRR) intracellular receptors to detect pathogen-derived molecules and to induce defense responses. A multitaxa reconstruction of NB-domain associations allowed us to identify the first NB–LRR arrangement in the Chlorophyta division of the Viridiplantae. Our analysis points out that the basic NOD-like receptor (NLR) unit emerged in Chlorophytes by horizontal transfer and its diversification started from Toll/interleukin receptor–NB–LRR members. The operon-based genomic structure of Chromochloris zofingiensis NLR copies suggests a functional origin of NLR clusters. Moreover, the transmembrane signatures of NLR proteins in the unicellular alga C. zofingiensis support the hypothesis that the NLR-based immunity system of plants derives from a cell-surface surveillance system. Taken together, our findings suggest that NLRs originated in unicellular algae and may have a common origin with cell-surface LRR receptors.

Plant and Animal Innate Immunity Complexes: Fighting Different Enemies with Similar Weapons

Both animals and plants express intracellular innate immunity receptors known as NLR (NOD-like receptors or nucleotide-binding domain and leucine-rich repeat receptors, respectively). For various mammalian systems, the specific formation of macromolecular structures, such as inflammasomes by activated NLR receptors, has been extensively reported. However, for plant organisms, the formation of such structures was an open scientific question for many years. This year, the first plant 'resistosome' structure was reported, revealing significant structural similarities to mammalian apoptosome and inflammasome structures. In this review, we summarize the key components comprising the mammalian apoptosome/inflammasome structures and the newly discovered plant resistosome, highlighting their commonalities and differences.

A functional investigation of the suppression of CpG and UpA dinucleotide frequencies in plant RNA virus genomes

Frequencies of CpG and UpA dinucleotides in most plant RNA virus genomes show degrees of suppression comparable to those of vertebrate RNA viruses. While pathways that target CpG and UpAs in HIV-1 and echovirus 7 genomes and restrict their replication have been partly characterised, whether an analogous process drives dinucleotide underrepresentation in plant viruses remains undetermined. We examined replication phenotypes of compositionally modified mutants of potato virus Y (PVY) in which CpG or UpA frequencies were maximised in non-structural genes (including helicase and polymerase encoding domains) while retaining protein coding. PYV mutants with increased CpG dinucleotide frequencies showed a dose-dependent reduction in systemic spread and pathogenicity and up to 1000-fold attenuated replication kinetics in distal sites on agroinfiltration of tobacco plants (Nicotiana benthamiana). Even more extraordinarily, comparably modified UpA-high mutants displayed no pathology and over a million-fold reduction in replication. Tobacco plants with knockdown of RDP6 displayed similar attenuation of CpG- and UpA-high mutants suggesting that restriction occurred independently of the plant siRNA antiviral responses. Despite the evolutionary gulf between plant and vertebrate genomes and encoded antiviral strategies, these findings point towards the existence of novel virus restriction pathways in plants functionally analogous to innate defence components in vertebrate cells.

Evolutionary balance between LRR domain loss and young NBS-LRR genes production governs disease resistance in Arachis hypogaea cv. Tifrunner

BACKGROUND

Cultivated peanut (Arachis hypogaea L.) is an important oil and protein crop, but it has low disease resistance; therefore, it is important to reveal the number, sequence features, function, and evolution of genes that confer resistance. Nucleotide-binding site-leucine-rich repeats (NBS-LRRs) are resistance genes that are involved in response to various pathogens.

RESULTS

We identified 713 full-length NBS-LRRs in A. hypogaea cv. Tifrunner. Genetic exchange events occurred on NBS-LRRs in A. hypogaea cv. Tifrunner, which were detected in the same subgenomes and also found in different subgenomes. Relaxed selection acted on NBS-LRR proteins and LRR domains in A. hypogaea cv. Tifrunner. Using quantitative trait loci (QTL), we found that NBS-LRRs were involved in response to late leaf spot, tomato spotted wilt virus, and bacterial wilt in A. duranensis (2 NBS-LRRs), A. ipaensis (39 NBS-LRRs), and A. hypogaea cv. Tifrunner (113 NBS-LRRs). In A. hypogaea cv. Tifrunner, 113 NBS-LRRs were classified as 75 young and 38 old NBS-LRRs, indicating that young NBS-LRRs were involved in response to disease after tetraploidization. However, compared to A. duranensis and A. ipaensis, fewer LRR domains were found in A. hypogaea cv. Tifrunner NBS-LRR proteins, partly explaining the lower disease resistance of the cultivated peanut.

CONCLUSIONS

Although relaxed selection acted on NBS-LRR proteins and LRR domains, LRR domains were preferentially lost in A. hypogaea cv. Tifrunner compared to A. duranensis and A. ipaensis. The QTL results suggested that young NBS-LRRs were important for resistance against diseases in A. hypogaea cv. Tifrunner. Our results provid insight into the greater susceptibility of A. hypogaea cv. Tifrunner to disease compared to A. duranensis and A. ipaensis.

Draft genome sequence of Solanum aethiopicum provides insights into disease resistance, drought tolerance, and the evolution of the genome

BACKGROUND

The African eggplant (Solanum aethiopicum) is a nutritious traditional vegetable used in many African countries, including Uganda and Nigeria. It is thought to have been domesticated in Africa from its wild relative, Solanum anguivi. S. aethiopicum has been routinely used as a source of disease resistance genes for several Solanaceae crops, including Solanum melongena. A lack of genomic resources has meant that breeding of S. aethiopicum has lagged behind other vegetable crops.

RESULTS

We assembled a 1.02-Gb draft genome of S. aethiopicum, which contained predominantly repetitive sequences (78.9%). We annotated 37,681 gene models, including 34,906 protein-coding genes. Expansion of disease resistance genes was observed via 2 rounds of amplification of long terminal repeat retrotransposons, which may have occurred ∼1.25 and 3.5 million years ago, respectively. By resequencing 65 S. aethiopicum and S. anguivi genotypes, 18,614,838 single-nucleotide polymorphisms were identified, of which 34,171 were located within disease resistance genes. Analysis of domestication and demographic history revealed active selection for genes involved in drought tolerance in both "Gilo" and "Shum" groups. A pan-genome of S. aethiopicum was assembled, containing 51,351 protein-coding genes; 7,069 of these genes were missing from the reference genome.

CONCLUSIONS

The genome sequence of S. aethiopicum enhances our understanding of its biotic and abiotic resistance. The single-nucleotide polymorphisms identified are immediately available for use by breeders. The information provided here will accelerate selection and breeding of the African eggplant, as well as other crops within the Solanaceae family.

Large-scale identification and functional analysis of NLR genes in blast resistance in the Tetep rice genome sequence

Tetep is a rice cultivar known for broad-spectrum resistance to blast, a devastating fungal disease. The molecular basis for its broad-spectrum resistance is still poorly understood. Is it because Tetep has many more NLR genes than other cultivars? Or does Tetep possess multiple major NLR genes that can individually confer broad-spectrum resistance to blast? Moreover, are there many interacting NLR pairs in the Tetep genome? We sequenced its genome, obtained a high-quality assembly, and annotated 455 nucleotide-binding site leucine-rich repeat (NLR) genes. We cloned and tested 219 NLR genes as transgenes in 2 susceptible cultivars using 5 to 12 diversified pathogen strains; in many cases, fewer than 12 strains were successfully cultured for testing. Ninety cloned NLRs showed resistance to 1 or more pathogen strains and each strain was recognized by multiple NLRs. However, few NLRs showed resistance to >6 strains, so multiple NLRs are apparently required for Tetep's broad-spectrum resistance to blast. This was further supported by the pedigree analyses, which suggested a correlation between resistance and the number of Tetep-derived NLRs. In developing a method to identify NLR pairs each of which functions as a unit, we found that >20% of the NLRs in the Tetep and 3 other rice genomes are paired. Finally, we designed an extensive set of molecular markers for rapidly introducing clustered and paired NLRs in the Tetep genome for breeding new resistant cultivars. This study increased our understanding of the genetic basis of broad-spectrum blast resistance in rice.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

Stronger When Together: Clustering of Plant NLR Disease resistance Genes

Gene clustering is rare in eukaryotes. However, nucleotide-binding leucine-rich repeat (NLR)-encoding disease resistance (R) genes show consistent clustering in plant genomes. These arrangements are likely to provide coregulatory benefits, as suggested by growing evidence that the gene products of both paired and larger clusters of NLRs act together in triggering immunity. Head-to-head gene pairs where one of the encoded NLRs includes an integrated decoy domain appear to behave differently than clusters evolved from closely related typical NLRs. These patterns may help to explain the broad resistance that most plants have despite their finite number of R genes. By taking into consideration the relationship between genomic arrangement and function, we can improve our understanding of and ability to predict plant immune detection.

Ligand-triggered allosteric ADP release primes a plant NLR complex

Pathogen recognition by nucleotide-binding (NB), leucine-rich repeat (LRR) receptors (NLRs) plays roles in plant immunity. The Xanthomonas campestris pv. campestris effector AvrAC uridylylates the Arabidopsis PBL2 kinase, and the latter (PBL2^UMP) acts as a ligand to activate the NLR ZAR1 precomplexed with the RKS1 pseudokinase. Here we report the cryo-electron microscopy structures of ZAR1-RKS1 and ZAR1-RKS1-PBL2^UMP in an inactive and intermediate state, respectively. The ZAR1^LRR domain, compared with animal NLR^LRR domains, is differently positioned to sequester ZAR1 in an inactive state. Recognition of PBL2^UMP is exclusively through RKS1, which interacts with ZAR1^LRR PBL2^UMP binding stabilizes the RKS1 activation segment, which sterically blocks ZAR1 adenosine diphosphate (ADP) binding. This engenders a more flexible NB domain without conformational changes in the other ZAR1 domains. Our study provides a structural template for understanding plant NLRs.

Animal NLRs continue to inform plant NLR structure and function

Plant NLRs share many of the structural hallmarks of their animal counterparts. At a functional level, the central nucleotide-binding pocket appears to have binding and hydrolysis activities, similar to that of animal NLRs. The TIR domains of plant NLRs have been shown to self-associate, and there is emerging evidence that full-length plant NLRs may do so as well. It is therefore tempting to speculate that plant NLRs may form higher-order complexes similar to those of the mammalian inflammasome. Here we review the available knowledge on structure-function relationships in plant NLRs, focusing on how the information available on animal NLRs informs the mechanism of plant NLR function, and highlight the evidence that innate immunity signalling pathways in multicellular organisms often require the formation of higher-order protein complexes.