Background: Disease resistance (R) genes from different Rosaceae species have been identified by map-based cloning for resistance breeding. However, there are few reports describing the pattern of R-gene evolution in Rosaceae species because several Rosaceae genome sequences have only recently become available.
Results: Since most disease resistance genes encode NBS-LRR proteins, we performed a systematic genome-wide survey of NBS-LRR genes between five Rosaceae species, namely Fragaria vesca (strawberry), Malus × domestica (apple), Pyrus bretschneideri (pear), Prunus persica (peach) and Prunus mume (mei) which contained 144, 748, 469, 354 and 352 NBS-LRR genes, respectively. A high proportion of multi-genes and similar Ks peaks (Ks = 0.1- 0.2) of gene families in the four woody genomes were detected. A total of 385 species-specific duplicate clades were observed in the phylogenetic tree constructed using all 2067 NBS-LRR genes. High percentages of NBS-LRR genes derived from species-specific duplication were found among the five genomes (61.81% in strawberry, 66.04% in apple, 48.61% in pear, 37.01% in peach and 40.05% in mei). Furthermore, the Ks and Ka/Ks values of TIR-NBS-LRR genes (TNLs) were significantly greater than those of non-TIR-NBS-LRR genes (non-TNLs), and most of the NBS-LRRs had Ka/Ks ratios less than 1, suggesting that they were evolving under a subfunctionalization model driven by purifying selection.
Conclusions: Our results indicate that recent duplications played an important role in the evolution of NBS-LRR genes in the four woody perennial Rosaceae species. Based on the phylogenetic tree produced, it could be inferred that species-specific duplication has mainly contributed to the expansion of NBS-LRR genes in the five Rosaceae species. In addition, the Ks and Ka/Ks ratios suggest that the rapidly evolved TNLs have different evolutionary patterns to adapt to different pathogens compared with non-TNL resistant genes.