Transcriptome analysis of newly classified bZIP transcription factors of Brassica rapa in cold stress response

•136 putative bZIP transcription members were identified in Brassica rapa.•LCRs are the significant region for molecular functions of BrbZIP proteins.•75 AtbZIP and 136 BrbZIP proteins were classified into nine groups according to UARR and LCRs.•The expression data of Brassica rapa bZIP genes to cold stress response was displayed.•Six genes containing N-rich domain might be involved in cold stress response.

Plant bZIP transcription factors play crucial roles in biological processes. In this study, 136 putative bZIP transcription members were identified in Brassica rapa. The bZIP family can be divided into nine groups according to the specific amino acid rich domain in B. rapa and Arabidopsis thaliana. To screen the cold stress responsive BrbZIP genes, we evaluated whether the transcription patterns of the BrbZIP genes were enhanced by cold treatment in the inbred lines, Chiifu and Kenshin, by microarray data analysis and qRT-PCR. The expression level of six genes increased significantly in Kenshin, but these genes were unchanged in Chiifu. These findings suggest that the six genes that encoded proteins containing N-rich regions might be involved in cold stress response. The results presented herein provide valuable information regarding the molecular basis of the bZIP transcription factors and their potential function in regulation growth and development, particularly in cold stress response.