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.

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

Nucleotide-binding leucine-rich-repeat (NLR) genes comprise the largest family of plant disease-resistance genes. Angiosperm NLR genes are phylogenetically divided into the TNL, CNL, and RNL subclasses. NLR copy numbers and subclass composition vary tremendously across angiosperm genomes. However, the evolutionary associations between genomic NLR content and ecological adaptation, or between NLR content and signal transduction components, are poorly characterized because of limited genome availability. In this study, we established an angiosperm NLR atlas (ANNA, https://biobigdata.nju.edu.cn/ANNA/) that includes NLR genes from over 300 angiosperm genomes. Using ANNA, we revealed that NLR copy numbers differ up to 66-fold among closely related species owing to rapid gene loss and gain. Interestingly, NLR contraction was associated with adaptations to aquatic, parasitic, and carnivorous lifestyles. The convergent NLR reduction in aquatic plants resembles the lack of NLR expansion during the long-term evolution of green algae before the colonization of land. A co-evolutionary pattern between NLR subclasses and plant immune pathway components was also identified, suggesting that immune pathway deficiencies may drive TNL loss. Finally, we identified a conserved TNL lineage that may function independently of the EDS1-SAG101-NRG1 module. Collectively, these findings provide new insights into the evolution of NLR genes in the context of ecological adaptation and genome content variation.

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.

Distinct Patterns of Gene Gain and Loss: Diverse Evolutionary Modes of NBS-Encoding Genes in Three Solanaceae Crop Species

Plant resistance conferred by nucleotide binding site (NBS)-encoding resistance genes plays a key role in the defense against various pathogens throughout the entire plant life cycle. However, comparative analyses for the systematic evaluation and determination of the evolutionary modes of NBS-encoding genes among Solanaceae species are rare. In this study, 447, 255, and 306 NBS-encoding genes were identified from the genomes of potato, tomato, and pepper, respectively. These genes usually clustered as tandem arrays on chromosomes; few existed as singletons. Phylogenetic analysis indicated that three subclasses [TNLs (TIR-NBS-LRR), CNLs (CC-NBS-LRR), and RNLs (RPW8-NBS-LRR)] each formed a monophyletic clade and were distinguished by unique exon/intron structures and amino acid motif sequences. By comparing phylogenetic and systematic relationships, we inferred that the NBS-encoding genes in the present genomes of potato, tomato, and pepper were derived from 150 CNL, 22 TNL, and 4 RNL ancestral genes, and underwent independent gene loss and duplication events after speciation. The NBS-encoding genes therefore exhibit diverse and dynamic evolutionary patterns in the three Solanaceae species, giving rise to the discrepant gene numbers observed today. Potato shows a “consistent expansion” pattern, tomato exhibits a pattern of “first expansion and then contraction,” and pepper presents a “shrinking” pattern. The earlier expansion of CNLs in the common ancestor led to the dominance of this subclass in gene numbers. However, RNLs remained at low copy numbers due to their specific functions. Along the evolutionary process of NBS-encoding genes in Solanaceae, species-specific tandem duplications contributed the most to gene expansions.

Large-Scale Analyses of Angiosperm Nucleotide-Binding Site-Leucine-Rich Repeat Genes Reveal Three Anciently Diverged Classes with Distinct Evolutionary Patterns

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes make up the largest plant disease resistance gene family (R genes), with hundreds of copies occurring in individual angiosperm genomes. However, the expansion history of NBS-LRR genes during angiosperm evolution is largely unknown. By identifying more than 6,000 NBS-LRR genes in 22 representative angiosperms and reconstructing their phylogenies, we present a potential framework of NBS-LRR gene evolution in the angiosperm. Three anciently diverged NBS-LRR classes (TNLs, CNLs, and RNLs) were distinguished with unique exon-intron structures and DNA motif sequences. A total of seven ancient TNL, 14 CNL, and two RNL lineages were discovered in the ancestral angiosperm, from which all current NBS-LRR gene repertoires were evolved. A pattern of gradual expansion during the first 100 million years of evolution of the angiosperm clade was observed for CNLs. TNL numbers remained stable during this period but were eventually deleted in three divergent angiosperm lineages. We inferred that an intense expansion of both TNL and CNL genes started from the Cretaceous-Paleogene boundary. Because dramatic environmental changes and an explosion in fungal diversity occurred during this period, the observed expansions of R genes probably reflect convergent adaptive responses of various angiosperm families. An ancient whole-genome duplication event that occurred in an angiosperm ancestor resulted in two RNL lineages, which were conservatively evolved and acted as scaffold proteins for defense signal transduction. Overall, the reconstructed framework of angiosperm NBS-LRR gene evolution in this study may serve as a fundamental reference for better understanding angiosperm NBS-LRR genes.