RGAugury: A pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants

Citation

Li, P., Quan, X., Jia, G., Xiao, J., Cloutier, S., You, F.M. (2016). RGAugury: A pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants. BMC Genomics, [online] 17(1), http://dx.doi.org/10.1186/s12864-016-3197-x

Plain language summary

Each plant cell contains a genome: a linear string of DNA base pairs which make each plant unique. Resistance gene analogs (RGAs) are potential R-genes in plants and can be predicted using a bioinformatics approach. Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. RGAs have conserved domains and motifs that play specific roles in pathogens' resistance. A user-friendly and efficient program, named RGAugury, has been developed for large-scale genome-wide RGA predictions of sequenced plant genomes. RGAugury was evaluated using the well-annotated Arabidopsis genome. RGAugury was also successfully applied to predict RGAs for 50 additional sequenced plant genomes. A user-friendly web interface was implemented to ease command line operations, facilitate visualization and simplify result management for multiple datasets.

Abstract

Background: Resistance gene analogs (RGAs), such as NBS-encoding proteins, receptor-like protein kinases (RLKs) and receptor-like proteins (RLPs), are potential R-genes that contain specific conserved domains and motifs. Thus, RGAs can be predicted based on their conserved structural features using bioinformatics tools. Computer programs have been developed for the identification of individual domains and motifs from the protein sequences of RGAs but none offer a systematic assessment of the different types of RGAs. A user-friendly and efficient pipeline is needed for large-scale genome-wide RGA predictions of the growing number of sequenced plant genomes. Results: An integrative pipeline, named RGAugury, was developed to automate RGA prediction. The pipeline first identifies RGA-related protein domains and motifs, namely nucleotide binding site (NB-ARC), leucine rich repeat (LRR), transmembrane (TM), serine/threonine and tyrosine kinase (STTK), lysin motif (LysM), coiled-coil (CC) and Toll/Interleukin-1 receptor (TIR). RGA candidates are identified and classified into four major families based on the presence of combinations of these RGA domains and motifs: NBS-encoding, TM-CC, and membrane associated RLP and RLK. All time-consuming analyses of the pipeline are paralleled to improve performance. The pipeline was evaluated using the well-annotated Arabidopsis genome. A total of 98.5, 85.2, and 100% of the reported NBS-encoding genes, membrane associated RLPs and RLKs were validated, respectively. The pipeline was also successfully applied to predict RGAs for 50 sequenced plant genomes. A user-friendly web interface was implemented to ease command line operations, facilitate visualization and simplify result management for multiple datasets. Conclusions: RGAugury is an efficiently integrative bioinformatics tool for large scale genome-wide identification of RGAs. It is freely available at Bitbucket: https://bitbucket.org/yaanlpc/rgaugury.

Publication date

2016-11-02