Genome-Wide Characterization and Expression Analysis of Gene Family Under Abiotic Stress in .

gene family is one of the largest transcription factor families. It plays an important role in plant growth, metabolic, and environmental response. However, complete genome-wide investigation of gene family in remains unexplained. In this study, 66 putative genes in the genome of were identified. And their evolutionary classification, physicochemical properties, conserved domain, functional differentiation, and the expression level under different stress conditions were further analyzed. All the members were clustered into 13 subfamilies (A-K, M, and S). A total of 10 conserved motifs were found in GubZIP proteins. Members from the same subfamily shared highly similar gene structures and conserved domains. Tandem duplication events acted as a major driving force for the evolution of gene family in . Cis-acting elements and protein-protein interaction networks showed that in one subfamily are involved in multiple functions, while some from different subfamilies may share the same functional category. The miRNA network targeting GubZIPs showed that the regulation at the transcriptional level may affect protein-protein interaction networks. We suspected that domain-mediated interactions may categorize a protein family into subfamilies in G. uralensis. Furthermore, the tissue-specific gene expression patterns of GubZIPs were analyzed using the public RNA-seq data. Moreover, gene expression level of 66 bZIP family members under abiotic stress treatments was quantified by using qRT-PCR. The results of this study may serve as potential candidates for functional characterization in the future.