Abstract Background The basic leucine zipper (bZIP) transcription factor is one of the most abundant and conserved gene families in eukaryotes. In addition to participating in plant development and growth, bZIP transcription factors play crucial roles in various abiotic stress responses and anthocyanin accumulation. Up to now, analysis of bZIP gene family members in pomegranate (Punica granatum) has not been reported. Three published pomegranate genome sequences provide valuable resources for further gene function analysis. Results Using bioinformatics analysis, 65 PgbZIPs were identified and analyzed from the ‘Taishanhong’ pomegranate genome. We divided them into 13 groups (A, B, C, D, E, F, G, H, I, J, K, M, and S) according to the phylogenetic relationship with those of Arabidopsis, each containing a different number of genes. The regularity of exon/intron number and distribution was consistent with the classification of groups in the evolutionary tree. Transcriptome analysis of different tissues showed that members of the PgbZIP gene family were differentially expressed in different developmental stages and tissues of pomegranate. Among them, we selected PgbZIP16 and PgbZIP34 as candidate genes which affect anthocyanin accumulation. The full-length CDS region of PgbZIP16 and PgbZIP34 were cloned from pomegranate petals by homologous cloning technique, encoding 170 and 174 amino acids, which were 510 bp and 522 bp, respectively. Subcellular localization assays suggested that both PgbZIP16 and PgbZIP34 were nucleus-localized. Real-time quantitative PCR (qPCR) was used to explore the expression of PgbZIP16 and PgbZIP34 in the petals of three kinds of ornamental pomegranates at the full flowering stage. The results demonstrated that the expression of PgbZIP16 in red petals was 5.83 times of that in white petals, while PgbZIP34 was 3.9 times. The results of transient expression in tobacco showed that consistent trends were observed in anthocyanin concentration and expression levels of related genes, which both increased and then decreased. Both PgbZIP16 and PgbZIP34 could promote anthocyanin accumulation in tobacco leaves. We obtained transgenic strains overexpressing PgbZIP16, and the histochemical staining for GUS activity showed that overexpressed PgbZIP16 seedlings were expressed in the stem. Transgenic experiments indicated that overexpression of PgbZIP16 significantly upregulated UF3GT, ANS and DFR genes in Arabidopsis and enhanced anthocyanin accumulation. Conclusions The whole genome identification, gene structure, phylogeny, gene cloning, subcellular location and functional verification of the pomegranate bZIP gene family provide a theoretical foundation for the functional study of the PgbZIP gene family and candidate genes for anthocyanin biosynthesis.