Genome-wide identification and expression analysis of the bZIP transcription factors, and functional analysis in response to drought and cold stresses in pear (Pyrus breschneideri)

Background Transcription factors (TFs) are involved in many important biological processes, including cell stretching, histological differentiation, metabolic activity, seed storage, gene regulation, and response to abiotic and biotic stresses. Little is known about the functions, evolutionary history, and expression patterns of basic region-leucine zipper TF family genes in pear, despite the release of the genome of Chinese white pears ("Dangshansuli"). Results Overall, 92 bZIP genes were identified in the pear genome (Pyrus breschneideri). Of these, 83 were randomly distributed on all 17 chromosomes except chromosome 4, and the other 9 genes were located on loose scaffolding. The genes were divided into 14 subgroups. Whole-genome duplications, dispersed duplication, and purifying selection for whole-genome duplications are the main reasons for the expansion of the PbrbZIP gene family. The analysis of functional annotation enrichment indicated that most of the functions of PbrbZIP genes were enriched in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways involved in the abiotic stress response. Next, expression analysis and virus-induced gene silencing results indicated that PbrbZIP genes might play critical roles in response to drought and cold stresses, especially for the genes from subgroups A, C, G, I, and S. Conclusions Ninety-two PbrbZIP genes were identified from the pear genome and classified into 14 subgroups. PbrbZIP genes were mainly expanded from whole-genome duplications and dispersed duplications and retained by purifying selection. PbrbZIP genes were induced by cold and drought stresses and played important roles in drought and cold tolerance. These results provided useful information for further increasing the tolerance of pears to stresses and a foundation to study the cold and drought tolerance mechanism of PbrbZIP genes.