Evolution of the number of LRRs in plant disease resistance genes

The largest group of plant resistance (R) genes contain the regions that encode the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains (NBS-LRR genes). To gain new resistance, amino acid substitutions and changes in number of the LRRs that recognize the presence of pathogens are considered important. In this study, we focus on the evolution of the number of LRRs and analyze the genome data of five plant species, Arabidopsis thaliana, Oryza sativa, Medicago truncatula, Lotus japonicus and Populus trichocarpa. We first categorized the NBS-LRR genes in each species into groups and subgroups based on the phylogenetic relationships of their NBS domain sequences. Then we estimated the evolutionary rate of the number of LRRs relative to the synonymous divergence in the NBS domain sequences by a maximum likelihood method assuming the single stepwise mutation model. The estimates ranged from 4.5 to 600 and differed between groups in the same species or between species. This indicated different roles played by different groups of the NBS-LRR genes within a species or the effects of various life history characteristics, such as generation time, of the species. We also tested the fit of the model to the data using the variance of number of LRRs in each subgroup. In some subgroups in some plants (16 out of 174 subgroups), the results of simulation using the estimated rates significantly deviated from the observed data. Those subgroups may have undergone different modes of selection from the other subgroups.