Genome-wide identification of the BURP domain-containing genes in Phaseolus vulgaris

Plant-specific BURP domain-containing proteins have an essential role in the plant's development and stress responses. Although BURP domain-containing proteins have been identified in several plant species, genome-wide analysis of the BURP gene family has not been investigated in the common bean. In the present study, we identified 11 BURP family members in the common bean (Phaseolus vulgaris) genome with a comprehensive in silico analysis. Pairwise alignment and phylogenetic analyses grouped PvBURP members into four subfamilies [RD-22 like (3), PG1 beta-like (4), BNM2-like (3), and USP-like (1)] according to their amino acid motifs, protein domains and intron-exon structure. The physical and biochemical characteristics of amino acids, motif and intron-exon structure, and cis-regulatory elements of BURPs members were determined. Promoter regions of BURP members included stress, light, and hormone response-related cis-elements. Therefore, expression profiles of PvBURP genes were identified with in silico tools and qRT-PCR analyses under stress (salt and drought) and hormone treatment (ABA, IAA) in the current study. While significant activity changes were not observed in BURP genes in RNA-seq data sets related to salt stress, it was determined that some BURP genes were expressed differently in those with drought stress. We identified 12 different miRNA, including miRNA395, miRNA156, miRNA169, miRNA171, miRNA319, and miRNA390, targeting the nine PvBURP genes using two different in silico tools based on perfect or near-perfect complementarity to their targets. Here we present the first study to identify and characterize the BURP genes in common bean using whole-genome analysis, and the findings may serve as a reference for future functional research in common bean.