Genome-wide identification and evolutionary analyses of bZIP transcription factors in wheat and its relatives and expression profiles of anther development related TabZIP genes

Background: Among the largest and most diverse transcription factor families in plants, basic leucine zipper (bZIP) family participate in regulating various processes, including floral induction and development, stress and hormone signaling, photomorphogenesis, seed maturation and germination, and pathogen defense. Although common wheat (Triticum aestivum L.) is one of the most widely cultivated and consumed food crops in the world, there is no comprehensive analysis of bZIPs in wheat, especially those involved in anther development. Previous studies have demonstrated wheat, T. urartu, Ae. tauschii, barley and Brachypodium are evolutionarily close in Gramineae family, however, the real evolutionary relationship still remains mysterious. Results: In this study, 187 bZIP family genes were comprehensively identified from current wheat genome. 98, 96 and 107 members of bZIP family were also identified from the genomes of T. urartu, Ae. tauschii and barley, respectively. Orthology analyses suggested 69.4 % of TubZIPs were orthologous to 68.8 % of AetbZIPs and wheat had many more in-paralogs in the bZIP family than its relatives. It was deduced wheat had a closer phylogenetic relationship with barley and Brachypodium than T. urartu and Ae. tauschii. bZIP proteins in wheat, T. urartu and Ae. tauschii were divided into 14 subgroups based on phylogenetic analyses. Using Affymetrix microarray data, 48 differentially expressed TabZIP genes were identified to be related to anther development from comparison between the male sterility line and the restorer line. Genes with close evolutionary relationship tended to share similar gene structures. 15 of 23 selected TabZIP genes contained LTR elements in their promoter regions. Expression of 21 among these 23 TabZIP genes were obviously responsive to low temperature. These 23 TabZIP genes all exhibited distinct tissue-specific expression pattern. Among them, 11 TabZIP genes were predominantly expressed in anther and most of them showed over-dominance expression mode in the cross combination TY806 x BS366. Conclusions: The genome-wide identification provided an overall insight of bZIP gene family in wheat and its relatives. The evolutionary relationship of wheat and its relatives was proposed based on orthology analyses. Microarray and expression analyses suggested the potential involvement of bZIP genes in anther development and facilitated selection of anther development related gene for further functional characterization.