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.