Differentially expressed bZIP transcription factors confer multi-tolerances in Gossypium hirsutum L

Basic leucine zipper (bZIP) transcription factor plays an important role in various biological processes, such as response to biotic and abiotic stresses. In this study we performed a systematic investigation and analysis of bZIP gene family in Gossypium hirsutum to predict their functions in response to different abiotic stresses. A total of 207 bZIP genes were identified from Gossypium hirsutum genome and classified into 13 subfamilies through phylogenetic analysis, which was testified by the analysis of conserved motifs and exon-intron structures. Annotation of GHbZIPs was performed based on well-studied Arabidopsis bZIPs to speculate the gene function. RNA-seq analysis was conducted to identify the co-expressed and differentially expressed bZIPs under cold, heat, salt and PEG treatments. Promoter analysis and interaction network of GHbZIP proteins demonstrated that ABA-activated signaling pathway was pivotal in the regulation of GHbZIPs, and GHbZIPs involved in ER stress were supposed to function through interaction with other GHbZIPs and ABA pathway. Cis-elements in the upstream and downstream of GHbZIPs interaction network were also discussed. These findings provided us with clues about functions of bZIP in Gossypium hirsutum.