Diversity, abundance, and domain architecture of plant NLR proteins in Fabaceae

The resistance (R) gene family in plants is a vital component of the plant defense system, enabling host resistance against pathogens through interactions with pathogen effector proteins. These R genes often encode nucleotide-binding (NB-ARC or N) and leucine-rich-repeat (LRR or L) domains, collectively forming the NLR protein family. The NLR proteins have been widely explored in crops from Poaceae and Brassicaceae, but limited studies are available for crops in other families, including Fabaceae. To address this gap, we conducted a comprehensive genome-wide analysis of putative NLR proteins in nine Fabaceae crops, including Glycine max, Lupinus angustifolius, Medicago truncatula, Pisum sativum, Phaseolus vulgaris, Trifolium pratense, Vigna angularis, Vigna radiata, and Vigna unguiculata. Our study revealed a substantial variation in the number of NLR proteins, independent of genome size. Notably, the NB-ARC domain exhibited a preferential co-occurrence with a specific LRR domain (IPR001611) in Fabaceae. Furthermore, through protein signature analysis, we identified both species-specific and shared domains across the nine crops. By classifying the identified proteins into seven distinct classes (N, L, CN, TN, NL, CNL, and TNL), we observed species-specific clustering within the CN, TN, and CNL classes, reflecting the diversification of species within Fabaceae. This genome-wide study enhances our understanding of the NLR protein repertoire and comprehensive protein signatures in nine Fabaceae species and provides valuable insights into plant defense mechanisms.