An angiosperm NLR Atlas reveals that NLR gene reduction is associated with ecological specialization and signal transduction component deletion

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Co-expression network analysis and identification of core genes in the interaction between wheat and Puccinia striiformis f. sp. tritici

The epidemic of stripe rust, caused by the pathogen Puccinia striiformis f. sp. tritici (Pst), would reduce wheat (Triticum aestivum) yields seriously. Traditional experimental methods are difficult to discover the interaction between wheat and Pst. Multi-omics data analysis provides a new idea for efficiently mining the interactions between host and pathogen. We used 140 wheat-Pst RNA-Seq data to screen for differentially expressed genes (DEGs) between low susceptibility and high susceptibility samples, and carried out Gene Ontology (GO) enrichment analysis. Based on this, we constructed a gene co-expression network, identified the core genes and interacted gene pairs from the conservative modules. Finally, we checked the distribution of Nucleotide-binding and leucine-rich repeat (NLR) genes in the co-expression network and drew the wheat NLR gene co-expression network. In order to provide accessible information for related researchers, we built a web-based visualization platform to display the data. Based on the analysis, we found that resistance-related genes such as TaPR1, TaWRKY18 and HSP70 were highly expressed in the network. They were likely to be involved in the biological processes of Pst infecting wheat. This study can assist scholars in conducting studies on the pathogenesis and help to advance the investigation of wheat-Pst interaction patterns.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Dipterocarpoidae genomics reveal their demography and adaptations to Asian rainforests

Dipterocarpoideae species form the emergent layer of Asian rainforests. They are the indicator species for Asian rainforest distribution, but they are severely threatened. Here, to understand their adaptation and population decline, we assemble high-quality genomes of seven Dipterocarpoideae species including two autotetraploid species. We estimate the divergence time between Dipterocarpoideae and Malvaceae and within Dipterocarpoideae to be 108.2 (97.8‒118.2) and 88.4 (77.7‒102.9) million years ago, and we identify a whole genome duplication event preceding dipterocarp lineage diversification. We find several genes that showed a signature of selection, likely associated with the adaptation to Asian rainforests. By resequencing of two endangered species, we detect an expansion of effective population size after the last glacial period and a recent sharp decline coinciding with the history of local human activities. Our findings contribute to understanding the diversification and adaptation of dipterocarps and highlight anthropogenic disturbances as a major factor in their endangered status.

Comprehensive analysis of the aldehyde dehydrogenase gene family in Phaseolus vulgaris L. and their response to saline-alkali stress

Background

Aldehyde dehydrogenase (ALDH) scavenges toxic aldehyde molecules by catalyzing the oxidation of aldehydes to carboxylic acids. Although ALDH gene family members in various plants have been extensively studied and were found to regulate plant response to abiotic stress, reports on ALDH genes in the common bean (Phaseolus vulgaris L.) are limited. In this study, we aimed to investigate the effects of neutral (NS) and basic alkaline (AS) stresses on growth, physiological and biochemical indices, and ALDH activity, ALDH gene expression of common bean. In addition, We used bioinformatics techniques to analyze the physical and chemical properties, phylogenetic relationships, gene replication, collinearity, cis-acting elements, gene structure, motifs, and protein structural characteristics of PvALDH family members.

Results

We found that both NS and AS stresses weakened the photosynthetic performance of the leaves, induced oxidative stress, inhibited common bean growth, and enhanced the antioxidative system to scavenge reactive oxygen species. Furthermore, we our findings revealed that ALDH in the common bean actively responds to NS or AS stress by inducing the expression of PvALDH genes. In addition, using the established classification criteria and phylogenetic analysis, 27 PvALDHs were identified in the common bean genome, belonging to 10 ALDH families. The primary expansion mode of PvALDH genes was segmental duplication. Cis-acting elemental analysis showed that PvALDHs were associated with abiotic stress and phytohormonal responses. Gene expression analysis revealed that the PvALDH gene expression was tissue-specific. For instance, PvALDH3F1 and PvALDH3H1 were highly expressed in flower buds and flowers, respectively, whereas PvALDH3H2 and PvALDH2B4 were highly expressed in green mature pods and young pods, respectively. PvALDH22A1 and PvALDH11A2 were highly expressed in leaves and young trifoliates, respectively; PvALDH18B2 and PvALDH18B3 were highly expressed in stems and nodules, respectively; and PvALDH2C2 and PvALDH2C3 were highly expressed in the roots. PvALDHs expression in the roots responded positively to NS-AS stress, and PvALDH2C3, PvALDH5F1, and PvALDH10A1 were significantly (P < 0.05) upregulated in the roots.

Conclusion

These results indicate that AS stress causes higher levels of oxidative damage than NS stress, resulting in weaker photosynthetic performance and more significant inhibition of common bean growth. The influence of PvALDHs potentially modulates abiotic stress response, particularly in the context of saline-alkali stress. These findings establish a basis for future research into the potential roles of ALDHs in the common bean.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

The chromosome-level genome of the submerged plant Cryptocoryne crispatula provides insights into the terrestrial-freshwater transition in Araceae

Plant terrestrialization (i.e. the transition to a terrestrial environment) is a significant evolutionary event that has been intensively studied. While certain plant lineages, particularly in angiosperms, have re-adapted to freshwater habitats after colonizing terrene, however, the molecular mechanism of the terrestrial-freshwater (T-F) transition remains limited. Here, the basal monocot Araceae was selected as the study object to explore the T-F transition adaptation mechanism by comparative genomic analysis. Our findings revealed that the substitution rates significantly increased in the lineage of freshwater Araceae, which may promote their adaptation to the freshwater habitat. Additionally, 20 gene sets across all four freshwater species displayed signs of positive selection contributing to tissue development and defense responses in freshwater plants. Comparative synteny analysis showed that genes specific to submerged plants were enriched in cellular respiration and photosynthesis. In contrast, floating plants were involved in regulating gene expression, suggesting that gene and genome duplications may provide the original material for plants to adapt to the freshwater environment. Our study provides valuable insights into the genomic aspects of the transition from terrestrial to aquatic environments in Araceae, laying the groundwork for future research in the angiosperm.

Telomere-to-telomere and gap-free genome assembly of a susceptible grapevine species (Thompson Seedless) to facilitate grape functional genomics

Grapes are globally recognized as economically significant fruit trees. Among grape varieties, Thompson Seedless holds paramount influence for fresh consumption and for extensive applications in winemaking, drying, and juicing. This variety is one of the most efficient genotypes for grape genetic modification. However, the lack of a high-quality genome has impeded effective breeding efforts. Here, we present the high-quality reference genome of Thompson Seedless with all 19 chromosomes represented as 19 contiguous sequences (N50 = 27.1 Mb) with zero gaps and prediction of all telomeres and centromeres. Compared with the previous assembly (TSv1 version), the new assembly incorporates an additional 31.5 Mb of high-quality sequenced data with annotation of a total of 30 397 protein-coding genes. We also performed a meticulous analysis to identify nucleotide-binding leucine-rich repeat genes (NLRs) in Thompson Seedless and two wild grape varieties renowned for their disease resistance. Our analysis revealed a significant reduction in the number of two types of NLRs, TIR-NB-LRR (TNL) and CC-NB-LRR (CNL), in Thompson Seedless, which may have led to its sensitivity to many fungal diseases, such as powdery mildew, and an increase in the number of a third type, RPW8 (resistance to powdery mildew 8)-NB-LRR (RNL). Subsequently, transcriptome analysis showed significant enrichment of NLRs during powdery mildew infection, emphasizing the pivotal role of these elements in grapevine's defense against powdery mildew. The successful assembly of a high-quality Thompson Seedless reference genome significantly contributes to grape genomics research, providing insight into the importance of seedlessness, disease resistance, and color traits, and these data can be used to facilitate grape molecular breeding efforts.

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.

Immunobiodiversity: Conserved and specific immunity across land plants and beyond

Angiosperms represent most plants that humans cultivate, grow, and eat. However, angiosperms are only one of five major land plant lineages. As a whole lineage, plants also include algal groups. All these clades represent a tremendous genetic diversity that can be investigated to reveal the evolutionary history of any given mechanism. In this review, we describe the current model of the plant immune system, discuss its evolution based on the recent literature, and propose future directions for the field. In angiosperms, plant-microbe interactions have been intensively studied, revealing essential cell surface and intracellular immune receptors, as well as metabolic and hormonal defense pathways. Exploring diversity at the genomic and functional levels demonstrates the conservation of these pathways across land plants, some of which are beyond plants. On basis of the conserved mechanisms, lineage-specific variations have occurred, leading to diversified reservoirs of immune mechanisms. In rare cases, this diversity has been harnessed and successfully transferred to other species by integration of wild immune receptors or engineering of novel forms of receptors for improved resistance to pathogens. We propose that exploring further the diversity of immune mechanisms in the whole plant lineage will reveal completely novel sources of resistance to be deployed in crops.

The plant immune receptor SNC1 monitors helper NLRs targeted by a bacterial effector

Plants deploy intracellular receptors to counteract pathogen effectors that suppress cell-surface-receptor-mediated immunity. To what extent pathogens manipulate intracellular receptor-mediated immunity, and how plants tackle such manipulation, remains unknown. Arabidopsis thaliana encodes three similar ADR1 class helper nucleotide-binding domain leucine-rich repeat receptors (ADR1, ADR1-L1, and ADR1-L2), which are crucial in plant immunity initiated by intracellular receptors. Here, we report that Pseudomonas syringae effector AvrPtoB suppresses ADR1-L1- and ADR1-L2-mediated cell death. ADR1, however, evades such suppression by diversifying into two ubiquitination sites targeted by AvrPtoB. The intracellular sensor SNC1 interacts with and guards the CCR domains of ADR1-L1/L2. Removal of ADR1-L1/L2 or delivery of AvrPtoB activates SNC1, which then signals through ADR1 to trigger immunity. Our work elucidates the long-sought-after function of SNC1 in defense, and also how plants can use dual strategies, sequence diversification, and a multi-layered guard-guardee system, to counteract pathogen's attack on core immunity functions.

Evolution of antimicrobial cysteine-rich peptides in plants

KEY MESSAGE

We analyzed the evolutionary pattern of cysteine-rich peptides (CRPs) to infer the relationship between CRP copy number and plant ecotype, and the origin of bi-domains CRPs. Plants produce cysteine-rich peptides (CRPs) that have long-lasting broad-spectrum antimicrobial activity to protect themselves from various groups of pathogens. We analyzed 240 plant genomes, ranging from algae to eudicots, and discovered that CRPs are widely distributed in plants. Our comparative genomics results revealed that CRP genes have been amplified through both whole genome and local tandem duplication. The copy number of these genes varied significantly across lineages and was associated with the plant ecotype. This may be due to their resistance to changing pathogenic environments. The conserved and lineage-specific CRP families contribute to diverse antimicrobial activities. Furthermore, we investigated the unique bi-domain CRPs that result from unequal crossover events. Our findings provide a unique evolutionary perspective on CRPs and insights into their antimicrobial and symbiosis characteristics.

Plant cell surface immune receptors-Novel insights into function and evolution

Plants use surface resident and intracellular immune receptors to provide robust immunity against microbial infections. The contribution of the two receptor types to plant immunity differs spatially and temporally. The ongoing identification of new plant cell surface immune receptors and their microbial-derived immunogenic ligands reveal a previously unexpected complexity of plant surface sensors involved in the detection of specific microbial species. Comparative analyses of the plant species distribution of cell surface immune receptors indicate that plants harbor larger sets of genus- or species-specific surface receptors in addition to very few widespread pattern sensors. Leucine-rich repeat surface and intracellular immune sensors emerge as two polymorphic receptor classes whose evolutionary trajectories appear to be linked. This is consistent with their functional cooperativity in providing full plant immunity.

New Biochemical Principles for NLR Immunity in Plants

While working for the United States Department of Agriculture on the North Dakota Agricultural College campus in Fargo, North Dakota, in the 1940s and 1950s, Harold H. Flor formulated the genetic principles for coevolving plant host-pathogen interactions that govern disease resistance or susceptibility. His 'gene-for-gene' legacy runs deep in modern plant pathology and continues to inform molecular models of plant immune recognition and signaling. In this review, we discuss recent biochemical insights to plant immunity conferred by nucleotide-binding domain/leucine-rich-repeat (NLR) receptors, which are major gene-for-gene resistance determinants in nature and cultivated crops. Structural and biochemical analyses of pathogen-activated NLR oligomers (resistosomes) reveal how different NLR subtypes converge in various ways on calcium (Ca2+) signaling to promote pathogen immunity and host cell death. Especially striking is the identification of nucleotide-based signals generated enzymatically by plant toll-interleukin 1 receptor (TIR) domain NLRs. These small molecules are part of an emerging family of TIR-produced cyclic and noncyclic nucleotide signals that steer immune and cell-death responses in bacteria, mammals, and plants. A combined genetic, molecular, and biochemical understanding of plant NLR activation and signaling provides exciting new opportunities for combatting diseases in crops. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Contribution of Duplicated Nucleotide-Binding Leucine-Rich Repeat (NLR) Genes to Wheat Disease Resistance

Wheat has a large and diverse repertoire of NLRs involved in disease resistance, with over 1500 NLRs detected in some studies. These NLR genes occur as singletons or clusters containing copies of NLRs from different phylogenetic clades. The number of NLRs and cluster size can differ drastically among ecotypes and cultivars. Primarily, duplication has led to the evolution and diversification of NLR genes. Among the various mechanisms, whole genome duplication (WGD) is the most intense and leading cause, contributing to the complex evolutionary history and abundant gene set of hexaploid wheat. Tandem duplication or recombination is another major mechanism of NLR gene expansion in wheat. The diversity and divergence of duplicate NLR genes are responsible for the broad-spectrum resistance of most plant species with limited R genes. Understanding the mechanisms underlying the rapid evolution and diversification of wheat NLR genes will help improve disease resistance in crops. The present review focuses on the diversity and divergence of duplicate NLR genes and their contribution to wheat disease resistance. Moreover, we provide an overview of disease resistance-associated gene duplication and the underlying strategies in wheat.

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.

Plant immune signaling network mediated by helper NLRs

Plant nucleotide-binding leucine-rich repeat receptors (NLRs) are intracellular immune receptors for pathogen recognition and signaling. They include sensor NLRs (sNLRs) that detect pathogens, and helper NLRs, which transduce downstream immune signals. During immune responses, both membrane-localized pattern recognition receptors (PRRs) and sNLRs rely on helper NLRs for signal transduction. The Arabidopsis helper NLRs, ADR1s and NRG1s, along with their interacting lipase-like protein dimers, are differentially required by sNLRs. Recent structural and biochemical analyses suggest that they assemble into oligomeric resistosomes with lipase-like protein dimers upon perception of small molecules produced from enzymatic activities of upstream TIR-type sNLRs. As a result, ADR1s and NRG1s form membrane calcium channels to induce immune responses and cell death. In contrast, Solanaceous NRC clade helper NLRs transduce signals from many sNLRs and some PRRs. Here, we summarize the recent advances in plant helper NLR research, with a focus on their structural and biochemical mechanisms in immune signaling.

Genome-wide comparative analysis of the valine glutamine motif containing genes in four Ipomoea species

BACKGROUND

Genes with valine glutamine (VQ) motifs play an essential role in plant growth, development, and resistance to biotic and abiotic stresses. However, little information on the VQ genes in sweetpotato and other Ipomoea species is available.

RESULTS

This study identified 55, 58, 50 and 47 VQ genes from sweetpotato (I. batatas), I.triflida, I. triloba and I. nil, respectively. The phylogenetic analysis revealed that the VQ genes formed eight clades (I-VII), and the members in the same group exhibited similar exon-intron structure and conserved motifs distribution. The distribution of the VQ genes among the chromosomes of Ipomoea species was disproportional, with no VQ genes mapped on a few of each species' chromosomes. Duplication analysis suggested that segmental duplication significantly contributes to their expansion in sweetpotato, I.trifida, and I.triloba, while the segmental and tandem duplication contributions were comparable in I.nil. Cis-regulatory elements involved in stress responses, such as W-box, TGACG-motif, CGTCA-motif, ABRE, ARE, MBS, TCA-elements, LTR, and WUN-motif, were detected in the promoter regions of the VQ genes. A total of 30 orthologous groups were detected by syntenic analysis of the VQ genes. Based on the analysis of RNA-seq datasets, it was found that the VQ genes are expressed distinctly among different tissues and hormone or stress treatments. A total of 40 sweetpotato differentially expressed genes (DEGs) refer to biotic (sweetpotato stem nematodes and Ceratocystis fimbriata pathogen infection) or abiotic (cold, salt and drought) stress treatments were detected. Moreover, IbVQ8, IbVQ25 and IbVQ44 responded to the five stress treatments and were selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis.

CONCLUSIONS

Our study may provide new insights into the evolution of VQ genes in the four Ipomoea genomes and contribute to the future molecular breeding of sweetpotatoes.

Genome-wide analysis of NBS-LRR genes revealed contribution of disease resistance from Saccharum spontaneum to modern sugarcane cultivar

INTRODUCTION

During plant evolution, nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) genes have made significant contributions to plant disease resistance. With many high-quality plant genomes sequenced, identification and comprehensive analyses of NBS-LRR genes at whole genome level are of great importance to understand and utilize them.

METHODS

In this study, we identified the NBS-LRR genes of 23 representative species at whole genome level, and researches on NBS-LRR genes of four monocotyledonous grass species, Saccharum spontaneum, Saccharum officinarum, Sorghum bicolor and Miscanthus sinensis, were focused.

RESULTS AND DISCUSSION

We found that whole genome duplication, gene expansion, and allele loss could be factors affecting the number of NBS-LRR genes in the species, and whole genome duplication is likely to be the main cause of the number of NBS-LRR genes in sugarcane. Meanwhile, we also found a progressive trend of positive selection on NBS-LRR genes. These studies further elucidated the evolutionary pattern of NBS-LRR genes in plants. Transcriptome data from multiple sugarcane diseases revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern sugarcane cultivars, and the proportion was significantly higher than the expected. This finding reveals that S. spontaneum has a greater contribution to disease resistance for modern sugarcane cultivars. In addition, we observed allelespecific expression of seven NBS-LRR genes under leaf scald, and 125 NBS-LRR genes responding to multiple diseases were identified. Finally, we built a plant NBS-LRR gene database to facilitate subsequent analysis and use of NBSLRR genes obtained here. In conclusion, this study complemented and completed the research of plant NBS-LRR genes, and discussed how NBS-LRR genes responding to sugarcane diseases, which provided a guide and genetic resources for further research and utilization of NBS-LRR genes.

Genome-wide characterization of aldehyde dehydrogenase gene family members in groundnut (Arachis hypogaea) and the analysis under saline-alkali stress

Groundnut or peanut (Arachis hypogaea) is a legume crop. Its seeds are rich in protein and oil. Aldehyde dehydrogenase (ALDH, EC: 1.2.1.3) is an important enzyme involved in detoxification of aldehyde and cellular reactive oxygen species, as well as in attenuation of lipid peroxidation-meditated cellular toxicity under stress conditions. However, few studies have been identified and analyzed about ALDH members in Arachis hypogaea. In the present study, 71 members of the ALDH superfamily (AhALDH) were identified using the reference genome obtained from the Phytozome database. A systematic analysis of the evolutionary relationship, motif, gene structure, cis-acting elements, collinearity, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and expression patterns was conducted to understand the structure and function of AhALDHs. AhALDHs exhibited tissue-specific expression, and quantitative real-time PCR identified significant differences in the expression levels of AhALDH members under saline-alkali stress. The results revealed that some AhALDHs members could be involved in response to abiotic stress. Our findings on AhALDHs provide insights for further study.

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.

Variation in plant Toll/Interleukin-1 receptor domain protein dependence on ENHANCED DISEASE SUSCEPTIBILITY 1

Toll/Interleukin-1 receptor (TIR) domains are integral to immune systems across all kingdoms. In plants, TIRs are present in nucleotide-binding leucine-rich repeat (NLR) immune receptors, NLR-like, and TIR-only proteins. Although TIR-NLR and TIR signaling in plants require the ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) protein family, TIRs persist in species that have no EDS1 members. To assess whether particular TIR groups evolved with EDS1, we searched for TIR-EDS1 co-occurrence patterns. Using a large-scale phylogenetic analysis of TIR domains from 39 algal and land plant species, we identified 4 TIR families that are shared by several plant orders. One group occurred in TIR-NLRs of eudicots and another in TIR-NLRs across eudicots and magnoliids. Two further groups were more widespread. A conserved TIR-only group co-occurred with EDS1 and members of this group elicit EDS1-dependent cell death. In contrast, a maize (Zea mays) representative of TIR proteins with tetratricopeptide repeats was also present in species without EDS1 and induced EDS1-independent cell death. Our data provide a phylogeny-based plant TIR classification and identify TIRs that appear to have evolved with and are dependent on EDS1, while others have EDS1-independent activity.

Comparative-genomic analysis reveals dynamic NLR gene loss and gain across Apiaceae species

Introduction: Nucleotide-binding leucine-rich repeat (NLR) genes play a crucial role in green plants' responding to various pathogens. Genome-scale evolutionary studies of NLR genes are important for discovering and applying functional NLR genes. However, little is known about the evolution of NLR genes in the Apiaceae family including agricultural and medical plants. Methods: In this study, comparative genomic analysis was performed in four Apiaceae species to trace the dynamic evolutionary patterns of NLR genes during speciation in this family. Results: The results revealed different number of NLR genes in these four Apiaceae species, namely, Angelica sinensis (95), Coriandrum sativum (183), Apium graveolens (153) and Daucus carota (149). Phylogenetic analysis demonstrated that NLR genes in these four species were derived from 183 ancestral NLR lineages and experienced different levels of gene-loss and gain events. The contraction pattern of the ancestral NLR lineages was discovered during the evolution of D. carota, whereas a different pattern of contraction after first expansion of NLR genes was observed for A. sinensis, C. sativum and A. graveolens. Discussion: Taken together, rapid and dynamic gene content variation has shaped evolutionary history of NLR genes in Apiaceae species.

Genome-wide characterization of the NLR gene family in tomato (Solanum lycopersicum) and their relatedness to disease resistance

Nucleotide-binding leucine-rich-repeat receptors (NLR), the largest group of genes associated with plant disease resistance (R), have attracted attention due to their crucial role in protecting plants from pathogens. Genome-wide studies of NLRs have revealed conserved domains in the annotated tomato genome. The 321 NLR genes identified in the tomato genome have been randomly mapped to 12 chromosomes. Phylogenetic analysis and classification of NLRs have revealed that 211 genes share full-length domains categorized into three major clades (CNL, TNL, and RNL); the remaining 110 NLRs share partial domains and are classified in CN, TN, and N according to their motifs and gene structures. The cis-regulatory elements of NLRs exhibit the maximum number of these elements and are involved in response to biotic and abiotic stresses, pathogen recognition, and resistance. Analysis of the phylogenetic relationship between tomato NLRs and orthologs in other species has shown conservation among Solanaceae members and variation with A. thaliana. Synteny and Ka/Ks analyses of Solanum lycopersicum and Solanum tuberosum orthologs have underscored the importance of NLR conservation and diversification from ancestral species millions of years ago. RNA-seq data and qPCR analysis of early and late blight diseases in tomatoes revealed consistent NLR expression patterns, including upregulation in infected compared to control plants (with some exceptions), suggesting the role of NLRs as key regulators in early blight resistance. Moreover, the expression levels of NLRs associated with late blight resistance (Solyc04g007060 [NRC4] and Solyc10g008240 [RIB12]) suggested that they regulate S. lycopersicum resistance to P. infestans. These findings provide important fundamental knowledge for understanding NLR evolution and diversity and will empower the broader characterization of disease resistance genes for pyramiding through speed cloning to develop disease-tolerant varieties.

EDS1 complexes are not required for PRR responses and execute TNL-ETI from the nucleus in Nicotiana benthamiana

Heterodimeric complexes incorporating the lipase-like proteins EDS1 with PAD4 or SAG101 are central hubs in plant innate immunity. EDS1 functions encompass signal relay from TIR domain-containing intracellular NLR-type immune receptors (TNLs) towards RPW8-type helper NLRs (RNLs) and, in Arabidopsis thaliana, bolstering of signaling and resistance mediated by cell-surface pattern recognition receptors (PRRs). Increasing evidence points to the activation of EDS1 complexes by small molecule binding. We used CRISPR/Cas-generated mutant lines and agroinfiltration-based complementation assays to interrogate functions of EDS1 complexes in Nicotiana benthamiana. We did not detect impaired PRR signaling in N. benthamiana lines deficient in EDS1 complexes or RNLs. Intriguingly, in assays monitoring functions of SlEDS1-NbEDS1 complexes in N. benthamiana, mutations within the SlEDS1 catalytic triad could abolish or enhance TNL immunity. Furthermore, nuclear EDS1 accumulation was sufficient for N. benthamiana TNL (Roq1) immunity. Reinforcing PRR signaling in Arabidopsis might be a derived function of the TNL/EDS1 immune sector. Although Solanaceae EDS1 functionally depends on catalytic triad residues in some contexts, our data do not support binding of a TNL-derived small molecule in the triad environment. Whether and how nuclear EDS1 activity connects to membrane pore-forming RNLs remains unknown.

Evolution patterns of NBS genes in the genus Dendrobium and NBS-LRR gene expression in D. officinale by salicylic acid treatment

BACKGROUND

Dendrobium officinale Kimura et Migo, which contains rich polysaccharides, flavonoids and alkaloids, is a Traditional Chinese Medicine (TCM) with important economic benefits, while various pathogens have brought huge losses to its industrialization. NBS gene family is the largest class of plant disease resistance (R) genes, proteins of which are widely distributed in the upstream and downstream of the plant immune systems and are responsible for receiving infection signals and regulating gene expression respectively. It is of great significance for the subsequent disease resistance breeding of D. officinale to identify NBS genes by using the newly published high-quality chromosome-level D. officinale genome.

RESULTS

In this study, a total of 655 NBS genes were uncovered from the genomes of D. officinale, D. nobile, D. chrysotoxum, V. planifolia, A. shenzhenica, P. equestris and A. thaliana. The phylogenetic results of CNL-type protein sequences showed that orchid NBS-LRR genes have significantly degenerated on branches a and b. The Dendrobium NBS gene homology analysis showed that the Dendrobium NBS genes have two obvious characteristics: type changing and NB-ARC domain degeneration. Because the NBS-LRR genes have both NB-ARC and LRR domains, 22 D. officinale NBS-LRR genes were used for subsequent analyses, such as gene structures, conserved motifs, cis-elements and functional annotation analyses. All these results suggested that D. officinale NBS-LRR genes take part in the ETI system, plant hormone signal transduction pathway and Ras signaling pathway. Finally, there were 1,677 DEGs identified from the salicylic acid (SA) treatment transcriptome data of D. officinale. Among them, six NBS-LRR genes (Dof013264, Dof020566, Dof019188, Dof019191, Dof020138 and Dof020707) were significantly up-regulated. However, only Dof020138 was closely related to other pathways from the results of WGCNA, such as pathogen identification pathways, MAPK signaling pathways, plant hormone signal transduction pathways, biosynthetic pathways and energy metabolism pathways.

CONCLUSION

Our results revealed that the NBS gene degenerations are common in the genus Dendrobium, which is the main reason for the diversity of NBS genes, and the NBS-LRR genes generally take part in D. officinale ETI system and signal transduction pathways. In addition, the D. officinale NBS-LRR gene Dof020138, which may have an important breeding value, is indirectly activated by SA in the ETI system.

Activation and Regulation of NLR Immune Receptor Networks

Plants have many types of immune receptors that recognize diverse pathogen molecules and activate the innate immune system. The intracellular immune receptor family of nucleotide-binding domain leucine-rich repeat-containing proteins (NLRs) perceives translocated pathogen effector proteins and executes a robust immune response, including programmed cell death. Many plant NLRs have functionally specialized to sense pathogen effectors (sensor NLRs) or to execute immune signaling (helper NLRs). Sub-functionalized NLRs form a network-type receptor system known as the NLR network. In this review, we highlight the concept of NLR networks, discussing how they are formed, activated and regulated. Two main types of NLR networks have been described in plants: the ACTIVATED DISEASE RESISTANCE 1/N REQUIREMENT GENE 1 network and the NLR-REQUIRED FOR CELL DEATH network. In both networks, multiple helper NLRs function as signaling hubs for sensor NLRs and cell-surface-localized immune receptors. Additionally, the networks are regulated at the transcriptional and posttranscriptional levels, and are also modulated by other host proteins to ensure proper network activation and prevent autoimmunity. Plant pathogens in turn have converged on suppressing NLR networks, thereby facilitating infection and disease. Understanding the NLR immune system at the network level could inform future breeding programs by highlighting the appropriate genetic combinations of immunoreceptors to use while avoiding deleterious autoimmunity and suppression by pathogens.

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.

Connecting the dots between cell surface- and intracellular-triggered immune pathways in plants

Plants can detect microbial molecules via surface-localized pattern-recognition receptors (PRRs) and intracellular immune receptors from the nucleotide-binding, leucine-rich repeat receptor (NLR) family. The corresponding pattern-triggered (PTI) and effector-triggered (ETI) immunity were long considered separate pathways, although they converge on largely similar cellular responses, such as calcium influx and overlapping gene reprogramming. A number of studies recently uncovered genetic and molecular interconnections between PTI and ETI, highlighting the complexity of the plant immune network. Notably, PRR- and NLR-mediated immune responses require and potentiate each other to reach an optimal immune output. How PTI and ETI connect to confer robust immunity in different plant species, including crops will be an exciting future research area.

Concerted expansion and contraction of immune receptor gene repertoires in plant genomes

Recent reports suggest that cell-surface and intracellular immune receptors function synergistically to activate robust defence against pathogens, but whether they co-evolve is unclear. Here we determined the numbers of cell-surface and intracellular immune receptors in 350 species. Surprisingly, the number of receptor genes that are predicted to encode cell-surface and intracellular immune receptors is strongly correlated. We suggest this is consistent with mutual potentiation of immunity initiated by cell-surface and intracellular receptors being reflected in the concerted co-evolution of the size of their repertoires across plant species.

Genome-Wide Characterization of Nitrogenase Reductase (nifH) Genes in the Sweet Potato [Ipomoea batatas (L.) Lam] and Its Wild Ancestors

The sweet potato (Ipomoea batatas (L.) Lam.) is an important and widely grown crop, and the nitrogenase reductase (nifH) gene is the most widely sequenced marker gene used to identify nitrogen-fixing bacteria and archaea. There have been many examples of the isolation of the diazotrophic endophytes in sweet potatoes, and there has been no report on whether sweet potatoes and their wild ancestors harbored nifH genes. In this study, a comprehensive analysis of nifH genes has been conducted on these species by using bioinformatics and molecular biology methods. A total of 20, 19 and 17 nifH genes were identified for the first time in sweet potatoes, I. trifida and I. triloba, respectively. Based on a phylogenetic analysis, all of the nifH genes, except for g10233.t1, itf14g14040.t1 and itb14g15470.t1, were clustered into five independent clades: I, II, III, IV and V. The nifH genes clustered in the same phylogenetic branch showed a more similar distribution of conserved motifs and exons–introns than those of the other ones. All of the identified genes were further mapped on the 15 chromosomes of the sweet potato, I. trifida and I. triloba. No segmental duplication was detected in each genome of three Ipomoea species, and 0, 8 and 7 tandemly duplicated gene pairs were detected in the genome of the sweet potato, I. trifida and I. triloba, respectively. Synteny analysis between the three Ipomoea species revealed that there were 7, 7 and 8 syntenic gene pairs of nifH genes detected between the sweet potato and I. trifida, between the sweet potato and I. triloba and between I. trifida and I. triloba, respectively. All of the duplicated and syntenic nifH genes were subjected to purifying selection inside duplicated genomic elements during speciation, except for the tandemly duplicated gene pair itf11g07340.t2_itf11g07340.t3, which was subjected to positive selection. Different expression profiles were detected in the sweet potato, I. trifida and I. triloba. According to the above results, four nifH genes of the sweet potato (g950, g16683, g27094 and g33987) were selected for quantitative real-time polymerase chain reaction (qRT-PCR) analysis in two sweet potato cultivars (Eshu 15 and Long 9) under nitrogen deficiency (N0) and normal (N1) conditions. All of them were upregulated in the N1 treatment and were consistent with the analysis of the RNA-seq data. We hope that these results will provide new insights into the nifH genes in the sweet potato and its wild ancestors and will contribute to the molecular breeding of sweet potatoes in the future.

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.

The chromosome-level genome of a free-floating aquatic weed Pistia stratiotes provides insights into its rapid invasion

Pistia stratiotes (Araceae), commonly referred to as water lettuce, is one of the most notorious weeds that cause severe damage to the economy and natural ecosystems of infested areas. In order to explore the mechanism of its rapid invasion, here, we assembled a high-quality chromosome-level genome for P. stratiotes based on the Illumina sequencing, PacBio sequencing, and Hi-C scaffolding technology. The assembled genome is 311.87 Mb in size with a contig N50 of 1.08 Mb. The contigs were further anchored on 14 pseudochromosomes with a scaffold N50 of 21.21 Mb. A total of 20,356 protein-coding genes were predicted, of which 79.35% were functionally annotated here. Evolutionary analysis showed that P. stratiotes and Colocasia esculenta were clustered together as sister lineages that diverged approximately 61 Mya. The synteny analyses indicated that two whole-genome duplication (WGD) events occurred within a short period in P. stratiotes. Moreover, comparative genome analysis indicated that the expansion of gene families corresponding to disease resistance might contribute to rapid invasion in P. stratiotes. Also, we analyzed the disease-resistance gene family (NBS-LRR) involved in plant defense. A genome-wide search in P. stratiotes genome identified 85 NBS-LRR genes in this study. In conclusion, our present study provides some new insights into the evolution of the invasive aquatic plant P. stratiotes. Our reference genome will also provide valuable resources for future invasion genomics research programs.

Genome-Wide Analysis of NBS-LRR Genes From an Early-Diverging Angiosperm Euryale ferox

NBS-LRR genes are the largest gene family in plants conferring resistance to pathogens. At present, studies on the evolution of NBS-LRR genes in angiosperms mainly focused on monocots and eudicots, while studies on NBS-LRR genes in the basal angiosperms are limited. Euryale ferox represents an early-diverging angiosperm order, Nymphaeales, and confronts various pathogens during its lifetime, which can cause serious economic losses in terms of yield and quality. In this study, we performed a genome-wide identification and analysis of NBS-LRR genes in E. ferox. All 131 identified NBS-LRR genes could be divided into three subclasses according to different domain combinations, including 18 RNLs, 40 CNLs, and 73 TNLs. The E. ferox NBS-LRR genes are unevenly distributed on 29 chromosomes; 87 genes are clustered at 18 multigene loci, and 44 genes are singletons. Gene duplication analysis revealed that segmental duplications acted as a major mechanism for NBS-LRR gene expansions but not for RNL genes, because 18 RNL genes were scattered over 11 chromosomes without synteny loci, indicating that the expansion of RNL genes could have been caused by ectopic duplications. Ancestral gene reconciliation based on phylogenetic analysis revealed that there were at least 122 ancestral NBS-LRR lineages in the common ancestor of the three Nymphaeaceae species, suggesting that NBS-LRR genes expanded slightly during speciation in E. ferox. Transcriptome analysis showed that the majority of NBS-LRR genes were at a low level of expression without pathogen stimulation. Overall, this study characterized the profile of NBS-LRR genes in E. ferox and should serve as a valuable resource for disease resistance breeding in E. ferox.

Molecular innovations in plant TIR-based immunity signaling

A protein domain (Toll and Interleukin-1 receptor [TIR]-like) with homology to animal TIRs mediates immune signaling in prokaryotes and eukaryotes. Here, we present an overview of TIR evolution and the molecular versatility of TIR domains in different protein architectures for host protection against microbial attack. Plant TIR-based signaling emerges as being central to the potentiation and effectiveness of host defenses triggered by intracellular and cell-surface immune receptors. Equally relevant for plant fitness are mechanisms that limit potent TIR signaling in healthy tissues but maintain preparedness for infection. We propose that seed plants evolved a specialized protein module to selectively translate TIR enzymatic activities to defense outputs, overlaying a more general function of TIRs.

Thirty years of resistance: Zig-zag through the plant immune system

Understanding the plant immune system is crucial for using genetics to protect crops from diseases. Plants resist pathogens via a two-tiered innate immune detection-and-response system. The first plant Resistance (R) gene was cloned in 1992 . Since then, many cell-surface pattern recognition receptors (PRRs) have been identified, and R genes that encode intracellular nucleotide-binding leucine-rich repeat receptors (NLRs) have been cloned. Here, we provide a list of characterized PRRs and NLRs. In addition to immune receptors, many components of immune signaling networks were discovered over the last 30 years. We review the signaling pathways, physiological responses, and molecular regulation of both PRR- and NLR-mediated immunity. Recent studies have reinforced the importance of interactions between the two immune systems. We provide an overview of interactions between PRR- and NLR-mediated immunity, highlighting challenges and perspectives for future research.

The genome of the rice variety LTH provides insight into its universal susceptibility mechanism to worldwide rice blast fungal strains

The widely used rice variety Lijiangxintuanheigu (LTH) shows a universal susceptibility to thousands of Magnaporthe oryzae isolates, the causal agent of devastating rice blast, making LTH an ideal line in resistance (R) gene cloning. However, the underlying genetic mechanism of the universal susceptibility has not been fully revealed because of the lack of a high-quality genome. Here, we took a genomic approach together with experimental assays to investigate LTH’s universal susceptibility to rice blast. Using Nanopore long reads, we assembled a chromosome-level genome. Millions of genomic variants were detected by comparing LTH with 10 other rice varieties, of which large-effect variants could affect plant immunity. Gene family analyses show that the number of R genes and leucine-rich repeat receptor-like protein kinase (LRR-RLK)-encoding genes decrease significantly in LTH. Rice blast resistance genes called Pi genes are either absent or disrupted by genomic variations. Additionally, residual R genes of LTH are likely under weak pathogen selection pressure, and other plant defense-related genes are weakly induced by rice blast. In contrast, the pattern-triggered immunity (PTI) of LTH is normal, as demonstrated by experimental assays. Therefore, we conclude that weak effector-trigger immunity (ETI)-mediated primarily by Pi genes but not PTI results in the universal susceptibility of LTH to rice blast. The attenuated ETI of LTH may be also associated with reduced numbers of R genes and LRR-RLKs, and minimally functional residual defense-related genes. Finally, we demonstrate the use of the LTH genome by rapid cloning of the Pi gene Piak from a resistant variety.

Genome-wide comparative analysis of the nucleotide-binding site-encoding genes in four Ipomoea species

The nucleotide-binding site (NBS)-encoding gene is a major type of resistance (R) gene, and its diverse evolutionary patterns were analyzed in different angiosperm lineages. Until now, no comparative studies have been done on the NBS encoding genes in Ipomoea species. In this study, various numbers of NBS-encoding genes were identified across the whole genome of sweet potato (Ipomoea batatas) (#889), Ipomoea trifida (#554), Ipomoea triloba (#571), and Ipomoea nil (#757). Gene analysis showed that the CN-type and N-type were more common than the other types of NBS-encoding genes. The phylogenetic analysis revealed that the NBS-encoding genes formed three monophyletic clades: CNL, TNL, and RNL, which were distinguished by amino acid motifs. The distribution of the NBS-encoding genes among the chromosomes was non-random and uneven; 83.13, 76.71, 90.37, and 86.39% of the genes occurred in clusters in sweet potato, I. trifida, I. triloba, and I. nil, respectively. The duplication pattern analysis reveals the presence of higher segmentally duplicated genes in sweet potatoes than tandemly duplicated ones. The opposite trend was found for the other three species. A total of 201 NBS-encoding orthologous genes were found to form synteny gene pairs between any two of the four Ipomea species, suggesting that each of the synteny gene pairs was derived from a common ancestor. The gene expression patterns were acquired by analyzing using the published datasets. To explore the candidate resistant genes in sweet potato, transcriptome analysis has been carried out using two resistant (JK20 and JK274) and susceptible cultivars (Tengfei and Santiandao) of sweet potato for stem nematodes and Ceratocystis fimbriata pathogen, respectively. A total of 11 differentially expressed genes (DEGs) were found in Tengfei and JK20 for stem nematodes and 19 DEGs in Santiandao and JK274 for C. fimbriata. Moreover, six DEGs were further selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis. The results may provide new insights into the evolution of NBS-encoding genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Genome-wide Identification and Evolutionary Analysis of NBS-LRR Genes From Secale cereale

Secale cereale is an important crop in the Triticeae tribe of the Poaceae family, and it has unique agronomic characteristics and genome properties. It possesses resistance to many diseases and serves as an important resource for the breeding of other Triticeae crops. We performed a genome-wide study on S. cereale to identify the largest group of plant disease resistance genes (R genes), the nucleotide-binding site-leucine-rich repeat receptor (NBS-LRR) genes. In its genome, 582 NBS-LRR genes were identified, including one from the RNL subclass and 581 from the CNL subclass. The NBS-LRR gene number in the S. cereale genome is greater than that in barley and the diploid wheat genomes. S. cereale chromosome 4 contains the largest number of NBS-LRR genes among the seven chromosomes, which is different from the pattern in barley and the genomes B and D of wheat but similar to that in the genome A of wheat. Further synteny analysis suggests that more NBS-LRR genes on chromosome 4 have been inherited from a common ancestor by S. cereale and the wheat genome A than the wheat genomes B and D. Phylogenetic analysis revealed that at least 740 NBS-LRR lineages are present in the common ancestor of S. cereale, Hordeum vulgare and Triticum urartu. However, most of them have only been inherited by one or two species, with only 65 of them preserved in all three species. The S. cereale genome inherited 382 of these ancestral NBS-LRR lineages, but 120 of them have been lost in both H. vulgare and T. urartu. This study provides the full NBS-LRR profile of the S. cereale genome, which is a resource for S. cereale breeding and indicates that S. cereale can be an important material for the molecular breeding of other Triticeae crops.

Frequent Gene Duplication/Loss Shapes Distinct Evolutionary Patterns of NLR Genes in Arecaceae Species

Nucleotide-binding leucine-rich repeat (NLR) genes play a key role in plant immune responses and have co-evolved with pathogens since the origin of green plants. Comparative genomic studies on the evolution of NLR genes have been carried out in several angiosperm lineages. However, most of these lineages come from the dicot clade. In this study, comparative analysis was performed on NLR genes from five Arecaceae species to trace the dynamic evolutionary pattern of the gene family during species speciation in this monocot lineage. The results showed that NLR genes from the genomes of Elaeis guineensis (262), Phoenix dactylifera (85), Daemonorops jenkinsiana (536), Cocos nucifera (135) and Calamus simplicifolius (399) are highly variable. Frequent domain loss and alien domain integration have occurred to shape the NLR protein structures. Phylogenetic analysis revealed that NLR genes from the five genomes were derived from dozens of ancestral genes. D. jenkinsiana and E. guineensis genomes have experienced “consistent expansion” of the ancestral NLR lineages, whereas a pattern of “first expansion and then contraction” of NLR genes was observed for P. dactylifera, C. nucifera and C. simplicifolius. The results suggest that rapid and dynamic gene content and structure variation have shaped the NLR profiles of Arecaceae species.

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.