Uncovering the dynamic evolution of nucleotide‐binding site‐leucine‐rich repeat (NBS‐LRR) genes in Brassicaceae

Plant genomes harbor dozens to hundreds of nucleotide‐binding site‐leucine‐rich repeat (NBS‐LRR) genes; however, the long‐term evolutionary history of these resistance genes has not been fully understood. This study focuses on five Brassicaceae genomes and the Carica papaya genome to explore changes in NBS‐LRR genes that have taken place in this Rosid II lineage during the past 72 million years. Various numbers of NBS‐LRR genes were identified from Arabidopsis lyrata (198), A. thaliana (165), Brassica rapa (204), Capsella rubella (127), Thellungiella salsuginea (88), and C. papaya (51). In each genome, the identified NBS‐LRR genes were found to be unevenly distributed among chromosomes and most of them were clustered together. Phylogenetic analysis revealed that, before and after Brassicaceae speciation events, both toll/interleukin‐1 receptor‐NBS‐LRR (TNL) genes and non‐toll/interleukin‐1 receptor‐NBS‐LRR (nTNL) genes exhibited a pattern of first expansion and then contraction, suggesting that both subclasses of NBS‐LRR genes were responding to pathogen pressures synchronically. Further, by examining the gain/loss of TNL and nTNL genes at different evolutionary nodes, this study revealed that both events often occurred more drastically in TNL genes. Finally, the phylogeny of nTNL genes suggested that this NBS‐LRR subclass is composed of two separate ancient gene types: RPW8‐NBS‐LRR and Coiled‐coil‐NBS‐LRR.