Printer Friendly

Genome-Wide Identification of bZIP Family Genes Involved in Drought and Heat Stresses in Strawberry (Fragaria vesca).

1. Introduction

The transcription factors containing bZIP domain [1,2] constitute one of the largest gene families in plants. The common feature, bZIP domain with length ranging from 60 to 80 amino acids (aa), includes two unique structures, a highly conserved DNA-binding basic and hinge region and a relatively diversified leucine zipper region [3]. The basic and hinge region consists of 18 aa residues with an invariant motif (N-[X.sub.7]-R/K-[X.sub.9]). In contrast, the downstream leucine zipper region is made up of heptad leucine repeats, a repeated pattern of several amino acids, or other hydrophobic amino acids, such as Ile, Phe, and Met [4]. The leucine zipper region is located exactly 9 amino acids downstream from the C-terminal of the basic region. It has been proved that the leucine zipper region is more favorable to constitute an amphipathic-helix with two helical turns within each heptad, because of its special amino acid composition [1,2].

To date, numerous bZIP members have been investigated and characterized by various eukaryotic genomes. Previous studies have identified 17 bZIP genes in Saccharomyces cerevisiae [5], 89 in rice [2], 56 in human [6], 131 in soybean [7], 125 in maize [4], 75 in Arabidopsis [8], 92 in Sorghum [9], 49 in Ricinus communis L. [10], and 55 in grapevine [11]. Considerable evidence showed that plant bZIP proteins play crucial roles in various aspects of biological processes, including organ and tissue developments [12,13], embryogenesis [14], seed maturation [15], cell elongation [16,17], and nitrogen/carbon and energy metabolism [18-20]. In the recent years, increasing evidence has also indicated that bZIP proteins take part in the regulation of plants' response to abiotic and biotic stresses, including phytohormone abscisic acid (ABA) signaling [21,22], osmosis [19,23], drought [24,25], high salinity [26,27], cold stresses [7,28], pathogen defense [29,30], and light irradiation [31,32].

Although the bZIP transcription factor family in strawberry (F. vesca), apple, and peach has been identified [33], only information on the evolutionary characterization of bZIP genes in Rosaceae was provided. Limited evidence is available regarding basic annotation of structural and functional features in the bZIP family among the economically important strawberry. The cultivated strawberry (F. x ananassa) is a model plant for nonclimacteric fruit species. F. x ananassa derived from 4 diploid ancestors has a relatively complex octaploid genome harboring 56 chromosomes (2n = 8x = 56). In contrast, our experiment system F. vesca, the diploid woodland strawberry, has been sequenced, which offers substantial advantages for molecular and physiological studies [34]. In this research, we performed comprehensive studies on expansion pattern, structural annotation, DNA-binding site specificity, dimerization properties, and expression levels in different stages and abiotic treatments to shed light on bZIP genes in F. vesca.

2. Materials and Methods

2.1. Identification of bZIP Genes in Woodland Strawberry. To generate a comprehensive list of F. vesca bZIP genes, F. vesca v1.1 genome was downloaded from Phytozome (https:// Firstly, a HMM profile of the bZIP domain (PF00170) was downloaded from Pfam ( Secondly, HMMER analysis was performed to search all F. vesca protein sequences and extracted the corresponding coding DNA sequence (CDS) and amino acid sequences. Thirdly, all of the putative bZIP domains were manually verified (E value < 1.0) within Pfam and SMART ( to confirm their presence and completeness of bZIP domain. Molecular weight, instability index, isoelectric points, aliphatic index, and grand average of hydropathicity (GRAVY) of candidate FvbZIP proteins were detected in ExPASy ProtParam server (

2.2. Phylogenetic Analysis and Structures of FvbZIP Gene. Phylogenetic analysis was used for classifying FvbZIP genes into ten clades [33]. The amino acid sequences of bZIP domains located in bZIP proteins from strawberry, A. thaliana (AtbZIP) [8], and rice (OsbZIP) [2] were used to generate a phylogenetic tree through ClustalW alignment and the unrooted Neighbor-joining (N-J) method using MEGA 6.0 [35]. N-J analysis with pairwise deletion and the Jones-Taylor-Thornton (JTT) model was performed using the 1000 replicate bootstraps to support for each node in phylogenetic tree. The exon-intron organizations with intron phase in FvbZIPs were performed by alignment within CDS sequence and corresponding genomic sequence and presented by GSDS 2.0 (

2.3. Conserved Motif Identification of FvbZIP Proteins. The motif distribution of FvbZIP proteins were extracted by MEME (, and parameters are set as follows: the largest number of discovered and conserved motifs set 12 and motif width was from 30 to 50 aa. The remaining parameter settings were kept default.

2.4. Dimerization Properties of FvbZIP Proteins. To describe the foundation of the dimerization stability and speculate the dimerization specificity of 50 FvbZIP members, we named special amino acid residue positions as g, a, b, c, d, e, and f to put in order of each heptad on the basis of standard nomenclature for the seven heptad repeats [36] (see Supplementary Figure S4 available online at The first heptad was manually arranged from four amino acid residues (corresponding amino acid located at g position in the first heptad) before appearance of the first leucine in the bZIP domain. Both boundaries of N-terminal and C-terminal within the FvbZIP leucine zipper region have been the same as applied to bZIP proteins in Arabidopsis [37].

2.5. Plant Materials and Treatments. F. vesca (National Clonal Germplasm Repository accession number PI664444) were cultivated in growth chamber with 22[degrees]C in a 12/12 h dark/light photoperiod. Strawberry plants were transferred into 42[degrees]C chamber for 0, 1, 6, 12, and 24 h posttreatment (hpt) for heat stress treatment. Drought treatment was simulated by depleting water after fully watering plants and sampled at 0 (water content on normalized soil is 83.96%), 48 (79.18%), 96 (63.41%), and 144 (36.78%) hpt. The strawberry plants were rewatered followed 144 h (74.78%) of drought treatment and sampled again 48 h later. All of the experiments were performed three times.

2.6. Quantitative RT-PCR Analysis. Transcriptional changes of FvbZIP genes responding to drought and heat treatments were determined by quantitative RT-PCR analysis on Stratagene Mx3000P Real-Time PCR system using SYBR[R] Premix Ex Taq[TM] (TaKaRa, Japan). The PCR amplification conditions were as follows: 95[degrees]C 10 min, 40 cycles of 15 sec at 95[degrees]C, and 1min at 60[degrees]C; at the end, the melting curve analysis was executed for verifying the specificity of the primer with the following stage: 95[degrees]C for 15 sec, 60[degrees]C for 1 min, and 95[degrees]C for 15 sec. The relative gene expression was determined using an interspacer 18S-26S strawberry gene as an endogenous control gene [38]. The relative expression levels of the target genes were calculated by the [2.sup.-[DELTA][DELTA]Ct] method [39]. Each specific primer pair was listed in Supplementary Table S1. Three independent biological replications were performed at each treatment point.

RPKM (reads per kilobase of exon per million mapped reads) of FvbZIP genes in different tissues and early stage fruit development were directly downloaded from available database ( (Table S5), which reflects mRNA levels. The expression levels (RPKM) for FvbZIP genes in yellow and red fruits of strawberry were from Tianbao Yang's lab in USDA-ARS. To compare the differential expression of bZIP genes (Table S5), the expression level of each gene in the different stages was read and normalized using multiexperiment viewer (MeV) software [40].

3. Results and Discussion

3.1. Identification and Characterization of Strawberry bZIP Family. A total of 50 FvbZIP genes were identified in strawberry genome in Phytozome. The characteristic parameters of all predicted FvbZIP proteins are listed in Table 1, including chromosome location, protein length, molecular weight, theoretical pI, instability index, aliphatic index, and grand average of hydropathicity (GRAVY). The length of amino acids of FvbZIPs is from 114 (mrna26148) to 1106 (mrna03778) aa with an average of 402.54 aa. Compared to those in Arabidopsis [8] (harboring 75 members with an average 321 aa in length) and rice [2] (harboring 89 members with an average 311 aa in length), strawberry has 50 bZIP members with longer average length than these two plants. After divergence of monocots and eudicots, lower frequency of bZIP candidates' evolution occurred in eudicots than monocots [9], which may explain the fact that strawberry and Arabidopsis have fewer bZIP members than rice. The predicted molecular weights of FvbZIP proteins range from 13.51 (mrna26148) to 122.71 (mrna03778) kDa, with pI from 4.73 (mrna29159) to 9.56 (mrna11666) (Table 1).

3.2. Expansion Pattern of bZIP Family in Strawberry, Arabidopsis, and Rice. To survey the extent of species-specific expansion of the bZIP genes in strawberry, Arabidopsis, and rice, we performed a joint phylogenetic analysis of strawberry, Arabidopsis, and rice bZIP proteins. As shown in Arabidopsis and rice analysis [2,8], the joint phylogenetic tree of three species grouped all bZIPs in ten distinct clades (A, B, C, D, E, F, G, H, I, and S) (Supplemental Figure S1A). Six bZIP genes in strawberry, Arabidopsis, and rice could not be conclusively mapped to any clade and was renamed U. We investigated some specific nodes which led to strawberry-, Arabidopsis-, and rice-specific classes within clades (red circle on node in Supplemental Figure S1A). As these nodes located in specific classes among strawberry, Arabidopsis, and rice represent the divergence point in the evolution and indicate the most recent common ancestors (MRCAs) prior to division, we found that some bZIP genes possibly have been present in the MRCAs of strawberry, Arabidopsis, and rice, but one or two species-specific bZIPs were lost in strawberry, Arabidopsis, or rice. Six classes contained only strawberry bZIP genes (Supplemental Figure S1A, red arrows), fifteen classes contained only Arabidopsis bZIP genes (Supplemental Figure S1A, green arrows), and thirty classes only contained rice bZIP genes (Supplemental Figure S1A, yellow arrows), suggesting that gene loss event has occurred in those classes in other species. The interspecies classes revealed that a parallel evolution of bZIP genes happened in three plant species and the orthologous bZIP proteins were observed to play a similar role [41]. The number of classes illustrated a fact that at least 85 ancestral bZIP genes existed before the strawberry, Arabidopsis, and rice split. After splitting from the MRCA, strawberry has lost the largest number of bZIP genes, followed by Arabidopsis and rice (Supplemental Figure S1B). A wide range interspecies comparison of strawberry, Arabidopsis, and rice revealed the species specific and nonspecific proteins.

3.3. Gene Structure of FvbZIP Genes. In order to observe the structure of the FvbZIP genes, we performed an analysis of the number and distribution of exon-intron, which are recorded to play key roles in gene family evolution and further promote to comprehend the emergence and evolution of a given gene [42,43]. As shown in Figure 1, 10 (20%) of total FvbZIP genes exhibit intronless, which occurs exclusively in group S, accounting for 88.9%. About 20% of peach, 19.1% of rice, 15.3% of maize, and 20.4% of castor bean bZIP genes were predicted to be intronless [2,4,10,44]. Among FvbZIP genes containing introns, the number of introns in the open reading frame (ORF) varied from one (mrna29159) to 20 (mrna03778), indicating a great diversity in the FvbZIP members (Figure 1). According to number, position, and splicing phase of the introns located in FvbZIP genes, the considerably conserved basic and hinge regions within the bZIP domain were divided into five patterns (a, b, c, d, and e) (Figure 2; Supplemental Figure S2). A previous study has showed that patterns of intron positions and splicing phases located at the bZIP domain regions were widely regarded as a promoter to understand homology and evolution of bZIP genes [2]. Based on intron location relative to the codon, introns can be classified into three patterns. Introns that do not interrupt the codons, between the first and second bases of codon and between the second and third bases of the codon, are termed phase 0 (P0), phase 1 (P1), and phase 2 (P2), respectively [45].

Patterns a (including 14 FvbZIPs) and e (including 14 FvbZIPs) were the most widespread (Figure 2). Pattern a had an intron in P0 within the hinge region behind the position -6 (Q). Both patterns b and c had one intron each followed the position -22 (R) located at the basic region, but the phase was different. Pattern b had an intron in P0, whereas pattern c had an intron in P2. Pattern d with two introns (both located at P0), one followed position -26 (K) of the basic region and the other one located at the hinge region at position -5, was interrupted by Lys and Ala. Pattern e had no intron distribution neither in the basic nor hinge region. Ten (71.2%) FvbZIPs within pattern g were intronless, while the other 4 had introns located at the outside of their basic and hinge regions (Figures 1 and 2). In particular, the position and phase were observed to be highly conserved when there was an intron interrupting in the hinge region. However, positions of introns located in the basic region were variably distributed in P0 or P2 (Figure 2; Supplemental Figure S2). The similar intron distribution and splicing patterns were observed in the bZIP domain regions of rice, castor bean, and maize, which declares that the intron pattern within the basic and hinge regions of bZIP genes were highly conserved in plants [2,4,10].

Though the lengths of intron were variable (Figure 1), most members clustered in the same clade classified above (Supplemental Figure S1) exhibited the same or similar intron number and splicing phase. For instance, same splicing phases and gene structures were observed in four FvbZIP members within clade A and eight members within clade S; five members in clade I also showed a similar splicing phase and gene structure. The intron pattern a was shared in the whole members in clade A, clade S consisted of the most intronless members showed the intron pattern e, members in clade C showed intron pattern c, and clade D (except for mrna31321) shared pattern d. In conclusion, the conserved splicing phases and gene structures in each clade were in agreement with the classification of FvbZIPs in phylogenetic analysis.

3.4. Conserved Motifs of FvbZIP Proteins. Besides the bZIP domain, some of the other multiple conserved motifs had been detected in Arabidopsis, rice, maize, caster bean, and peach [2,4,8,10,44]. As shown in Figure 3, a total of 12 conserved motifs were identified in FvbZIP genes. The consensus sequence and the amino acid width of those conserved motifs are given in Supplemental Table S2. It can be observed that some motifs are shared by several clades, such as motif 9 which was located in clades C and S, motif 4 shared by clades E and I, and motif 12 which was present in clades A and D. Furthermore, majority of conserved motifs appear in a specific clade, for instance, motifs 2, 3, 5, 6, and 11 in clade D; motifs 7 and 8 in clade A; and motif 10 in clade G (Figure 3). It is a proved hypothesis that the group-specific conserved motifs promote to decide specific functions of FvbZIP genes in each clade and foundation for the functional divergences between members from different clades [46-48]. The conserved motif analysis was further proved to be consistent with phylogenetic relationship and classification of FvbZIP genes.

Notably, motifs 2, 9, and 10 were identified as DOG1 (delay of germination), bZIP C, and MFMR (multifunctional mosaic region) according to Pfam databases, respectively. The previous study has found DOG1-corresponding gene to be related to seed dormancy in a form of quantitative trait locus [49]. It has been suggested that N-terminal half of MFMR motif is rather rich in proline residues and has been termed the PRD [50], and some of these motifs may play a key role in regulating protein-protein interactions [51]. Motif 10 with a Pro-rich domain was supposed to be involved in mediating protein-protein interactions [51]. Motifs 7, 8, and 12 represented potential casein kinase II (CKII) phosphorylation sites (S/TxxD/E) [52], presented as T[FL]DE, [TQS][LM][GC][GDE], and [ST]AE[EA]; TLGE and TLE[DE]; and TVDE, respectively. These CKII phosphorylation sites were considered to be involved in ABA-mediate responses [26,53]. A part of these three motifs also contains another phosphorylation site (R/KxxS/T) associated with Ca2+-dependent protein kinase, presented as R[EQ][ANST]S and GK[DNP][FL][GS], RQ[PQ]T, and RQG[SG] and [AR][TAP]LS in motifs 7, 8, and 12, respectively. Motif4 with a glutamine-rich domain has been proved to play an important role in transcriptional activation [54]. We can speculate a similar function of bZIP proteins due to conserved motifs shared by Arabidopsis, rice, castor bean, and maize, while the roles of other conserved motifs found in FvbZIP proteins are not yet clear.

3.5. DNA-Binding Site Specificity of FvbZIP Proteins. DNA-binding specificity owed by bZIP proteins is determined by certain key amino acids located at the most conserved basic and hinge region of the bZIP domain, and these two regions have a direct interaction with DNA cis-elements [55,56]. The annotation and classification of amino acids located at the basic and hinge region were characterized by diagnostic sequence structure in Arabidopsis and rice [2,8], which indicates that some highly conserved amino acids are present in each clade (Supplementary Figure S3). The first leucine of the leucine zipper region was designed as +1, and C-terminal amino acid of the hinge region was numbered -1 [55]. It is proved that some functional replacements on invariant sites (-18 and -10) by other amino acids will lead to new DNA-binding specificities [57]. In the strawberry bZIP family, these replacements infrequently happened and just occurred in one group. mrna30252 in clade U with a hydrophobic Ile (Isoleucine) residue at position -10 replaced arginine/lysine, indicating that perhaps they are not able to bind DNA or might hold a unique DNA-binding specificity. Furthermore, we can speculate the structural characteristics and DNA-binding specificity of FvbZIP proteins in a clade model, as listed in Supplementary Table S3. Expression level of a particular gene is sometimes determined by binding affinity to specific sites of promoters of genes, which is resulted of variations in their target sites in promoters of different genes [2]. Computational discovery and further experimental verification of DNA binding ability would promote to identify which genes are selectively activated by various FvbZIP genes.

3.6. Dimerization Properties of FvbZIP Proteins. The function of bZIP proteins requires dimerization between parallel coiled-coil constructions [58]. Amino acids present at positions a, d, e, and g near the leucine zipper interface play an important role in regulating oligomerization of leucine zipper domain as well as specificity and stability of dimerization [4]. Thus, we performed detailed characterization of the amino acids located at positions a, d, e, and g to predict their functions. The length of the FvbZIP leucine zipper is variable to range from two to nine heptads. About 18% of amino acids located at position a were Asns (asparagine) in 41 FvbZIP proteins (Figure 4(a)), which is similar to both 16% in humans and H. sapiens, but lesser than observed in Arabidopsis (22%) [37], rice (23%) [2], maize (22%) [4], and castor bean (26%) [10]. The highest frequency of Asn at position a was displayed in the second heptad (accounting for 44%) followed by the fifth heptad (accounting for 36%) (Figure 4(b)), which is also observed earlier for rice and Arabidopsis [2,37]. Since Asns can form a more stable N-N interaction at the a [left and right arrow] a' position than the other a position amino acids, high frequency of Asn at the a position promote to create homodimerizing Leu zippers among the bZIP family [59]. This result raises the possibility that a majority of FvbZIPs prefer to homodimerization. It is observed that minority of charged amino acids (R, K, E, and D) displayed in position a, which may contribute to form hetero-dimerization.

Amino acids located at positions a and d are generally hydrophobic residues and pack on the surface of "knobs and holes" pattern to establish a hydrophobic core area that is important to homo-/hetero-dimerization [60]. The frequency of the conserved Leu responsible for dimer stability at d position, was found in 68% of 50 FvbZIP proteins. The abundance of other hydrophobic amino acid (including I, V, and M) at d position was comparable in FvbZIPs (10%, Figure 4(a)), Arabidopsis (19%), and rice (17%). The Leu at position d is a type of highly conserved and stable aliphatic amino acid residue [60], which is critical for keeping the dimerization stability. compared with the AtbZIP proteins (the proportion of Leu at d position accounts for 56%) [37] and RcbZIP proteins (56%) [10], FvbZIPs showed more abundant Leu located at d position, which suggests that strawberry genes probably hold a stronger stability dimerization on Leu zipper than Arabidopsis.

It is observed that the abundance of charged amino acids displayed at e and g positions, such as acidic amino acids (including E and D) and basic amino acids (including R and K), accounts for 50% and 50%, respectively (Figure 4(a)). Compared to Arabidopsis (charged amino acids at e and g positions were, respectively, accounting for 41% and 53%), the frequencies of charged amino acids located at e position was higher but was lower in strawberry at g position. The charged amino acids located at positions g and e were replaced by alanine can result in creating tetramers instead of dimmers [61]. A greater number of charged amino acids (containing E, D, R, and K) located at e and g positions were found to flank the dimerization interface, which are thought to interact electrostatically to form salt bridges mediating the dimerization and determining the dimerization specificity as well as stability [62].

As dimerization properties of FvbZIP proteins were governed by different contributions of charged amino acids located at e and g positions, we calculated the presence of attractive and repulsive g [left and right arrow] e' pairs in each heptad of strawberry Leu zippers (Supplementary Figure S4) and displayed corresponding frequency histogram in Figure 4(c). Four groups (basic repulsive, acidic repulsive, +/- attractive, and -/+ attractive) were found in complete g [left and right arrow] e' pairs according to both the g and following e position amino acids are charged, while incomplete g [left and right arrow] e' pairs were defined as only one charged position in g or e [10]. The highest frequency of the complete g [left and right arrow] e' pairs was observed to locate at the first heptad for 46%, majority of which was found to be the basic repulsive g [left and right arrow] e' pair. In following three heptads, the frequencies of the complete g [left and right arrow] e' pairs showed a sharp decline. The frequency of complete g e pairs appeared to be increased in the fifth heptad, especially on the -/+ attractive gene pairs. Moreover, attractive g [left and right arrow] e' (including +/- attractive and -/+ attractive) pairs were found most frequently within the second, fifth, and ninth heptad. It is observed that the attractive g [left and right arrow] e' pairs had no tend to form homo- or hetero-dimerization, while the repulsive g [left and right arrow] e' pairs preferred to make a conformation hetero dimerization [62]. Furthermore, biophysical measurements also reveal that repulsive g [left and right arrow] e' pairs present a more primary role in driving dimerization specificity than attractive g [left and right arrow] e' pairs [4]. According to the observation of the dimerization properties in our study, the 50 FvbZIP proteins could be classified into 20 types (Supplementary Figure S4 and Table S4).

3.7. Expression of FvbZIP Genes in Different Organs. Many bZIP proteins have been shown to be involved in the processes of plant growth and development [63]. To explore the tissue-specific expression patterns of genomic scale FvbZIP genes, the transcriptome data were analyzed to determine those genes expressed in different organs and tissues (Figure 5). FvbZIP genes were showed to rank from highest (mrna21832 in style and yellow fruit) to lowest (mrna31321 in embryo) according to their differential expression across main tissues. The FvbZIP genes clustered in the same clade showed different tissue-specific expression patterns. mrna15193, mrna21832, mrna14942, and mrna02284 within clade S showed particularly high expression in most organisms, while the others of clade S presented a relatively low expression in all tissues except for mrna18282 in ghost (refers to the entire seed with its embryo removed) and ovule. The duplicated genes in clades (predicted in the Plant Genome Duplication Database) showed different transcript abundance in specific tissues, such as mrna32022 and mrna18928 in flowered and perianth (Figure 5), which is also supported by the duplicated bZIP genes in tomato [64]. The divergences in expression profiles between paralogs revealed that some of them may acquire new functions after duplication in the evolutionary process.

As an activator of downstream gene expression, members of clade A bZIP genes share function as important regulatory roles in the ABA signaling pathway, abiotic stresses, and tissues development in Arabidopsis [8]. MdbZIP26 was abundantly expressed in roots, stems, and apical buds [65]. OsABF2 showed expression in different tissues of rice and induction of drought, salinity, cold, and ABA. In fact, clade A bZIPs often have played key roles not only in seed development but also in fruit maturation [66]. We further analyzed the expression profiles of the 8 clade A FvbZIP genes in various organs and fruits using transcriptome data. Firstly, the members show different tissue-specific expression. For example, mrna00393 and mrna09110 are moderately expressed in most tissues, while mrna30280 is more enriched only in flower (Figure 5(a)). Secondly, the respective members of the duplicated FvbZIPs pair (mrna14556 and mrna08566, predicted in the Plant Genome Duplication Database) have different pattern specificities. The duplicated pair has a similar expression pattern in each tissue but a different pattern in red and yellow fruits (Figures 5(a) and 5(b)).

It is evident that although highly conserved on amino acids sequence, those duplicated genes are differently regulated [67]. Thirdly, the same gene show different expression profiles in different organs. mrna11837 is expressed at an extremely high level in seedling but very low in majority of other tissues (Figure 5(a)). Overall, the expression of clade A FvbZIP genes has different organ specificities, including a functional diversification coupled with the expansion of this gene family in strawberry. At least one clade A bZIP gene (mrna09110) is expressed in most investigated strawberry organs. It is consistent with previous studies in Arabidopsis and apple [8,65].

3.8. Expression Patterns of FvbZIP Genes in Leaves under Drought and Heat Treatments. We next performed an expression investigation of the clade A FvbZIP genes following exposure to drought and heat. A range of expression patterns was observed and some of the genes clearly upregulated, while others were downregulated (Figures 6(a) and 6(b)). As shown in Figure 6(a), for example, mran00393, mrna08566, mrna30280, and mrna11837 showed a significant upregulation after the plants were treated by drought for 2 days, while all genes in clade A were significantly downregulated at day 4 (Figure 6(a)). These four genes showed a higher expression level than the previous point after rewatering (Figure 6(a)). When exposed to drought, the expression of mrna14556 and mrna28250 decreased gradually (Figure 6(a)). Heat treatment caused mrna30280 transcription levels to increase in all points, and the mrna30280 transcription level peaked 12 h after the heat treatment was started (Figure 6(a)). Similarly, mrna00393 showed an obvious increase in expression at 12 h (Figure 6(a)). mrna14556 and mrna28250 kept low expression during heat treatment (Figure 6(a)).

It is worth noting that the expression of mrna11837 in most organs in Figure 5(a) and mrna30280 in fruit in Figure 5(b) were weak under normal growth conditions but both higher following drought and heat. We conclude that these two genes (mrna11837 and mrna30280) have a different pathway to respond to drought and heat treatments. Drought and heat responses in plants may be a separate pathway, which is consistent with the finding in grape [11]. However, mrna14556 expression was not detected in any organ, even in response to drought or heat. A similar observation was made in grape, apple, and sorghum with some bZIP genes [65,68,69], and we speculate that those genes may be degenerated into a nonfunctional gene, or it is involved in other processes that were not investigated here, such as other abiotic and biotic stresses [68,69], or they may be expressed at other organs and time points after treatment [65]. This can be supported by the fact that the bZIP ABF genes, ABF1 and ABF4, are induced by cold as well as ABF2 and ABF3 expression are induced by high salt in vegetative tissues [26,53]. The constitutive overexpression of ABF3 in rice and Arabidopsis consequently made the transgenic plants more tolerant to stress conditions [70,71]. Moreover, OsbZIP23 and OsbZIP72 play an important role in training drought tolerant rice and resistance to drought treatment by activating ABA signaling [21,22]. The genes upregulated by stresses was suppressed in mutants of OsbABF1 under ABA treatment, which indicates that OsbABF1 is related to abiotic stress responses through ABA pathway in rice. We concluded that the clade A FvbZIP genes showed a range of response patterns when exposed to drought and heat stresses, revealing that they may respond to stress signaling, which is also observed in Arabidopsis, apple, and grape [8,65,68].

4. Conclusion

We performed genome-wide analysis of the strawberry bZIP transcription factor family and conducted a detailed investigation of their evolutionary relationship, structural feature, and organization. In addition, we characterized the expression of some clade A genes following drought and heat treatment. These results may prove useful in developing strategies for the further improvement of stress tolerance in strawberry.

aa:      Amino acids
Ala:     Alanine
Asn:     Asparagine
bZIP:    Basic leucine zipper
DOG:     Delay of germination
Ile:     Isoleucine
Leu:     Leucine
MFMR:    Multifunctional mosaic region
MRCAs:   Most recent common ancestors.

Conflicts of Interest

The authors declare no conflict of interest.


This research was financially supported in part by the Priority Academic of Jiangsu Province, the Ministry of Education-Nanjing Agricultural University Project no. KYT201004, and a scholarship from the China Scholarship Council (CSC, no. 201506850034) to Xiao-Long Wang. It is also supported in part by the Tennessee Agricultural Experiment Station.


[1] W. Iandschultz, P. Johnson, and S. McKnight, "The leucine zipper: A hypothetical structure common to a new class of DNA binding protein," Science, vol. 240, no. 4860, pp. 1759-1764, 1988.

[2] A. Nijhawan, M. Jain, A. K. Tyagi, and J. P. Khurana, "Genomic survey and gene expression analysis of the basic leucine zipper transcription factor family in rice," Plant Physiology, vol. 146, no. 2, pp. 333-350, 2008.

[3] H. C. Hurst, "Transcription factors. 1: bZIP proteins," Protein Profile, vol. 1, no. 2, pp. 123-168, 1993.

[4] K. Wei, J. Chen, Y. Wang et al., "Genome-wide analysis of bZIP-encoding genes in maize," DNA Research, vol. 19, no. 6, pp. 463-476, 2012.

[5] J. Fassler, D. Landsman, A. Acharya, J. R. Moll, M. Bonovich, and C. Vinson, "B-ZIP proteins encoded by the drosophila genome: evaluation of potential dimerization partners," Genome Research, vol. 12, no. 8, pp. 1190-1200, 2002.

[6] C. Vinson, M. Myakishev, A. Acharya, A. A. Mir, J. R. Moll, and M. Bonovich, "Classification of human B-ZIP proteins based on dimerization properties," Molecular and Cellular Biology, vol. 22, no. 18, pp. 6321-6335, 2002.

[7] Y. Liao, H.-F. Zou, W. Wei et al., "Soybean GmbZIP44, GmbZIP62 and GmbZIP78 genes function as negative regulator of ABA signaling and confer salt and freezing tolerance in transgenic Arabidopsis," Planta, vol. 228, no. 2, pp. 225-240, 2008.

[8] M. Jakoby, B. Weisshaar, W. Droge-Laser et al., "bZIP transcription factors in Arabidopsis," Trends in Plant Science, vol. 7, no. 3, pp. 106-111, 2002.

[9] J. Vanitha and S. Ramachandran, "Genome-wide expansion and expression divergence of the basic leucine zipper transcription factors in higher plants with an emphasis on sorghum," Journal of Integrative Plant Biology, vol. 53,no. 3, p. 004, 2011.

[10] Z. Jin, W. Xu, and A. Liu, "Genomic surveys and expression analysis of bZIP gene family in castor bean (Ricinus communis L.)," Planta, vol. 239, no. 2, pp. 299-312, 2014.

[11] J. Liu, N. Chen, F. Chen et al., "Genome-wide analysis and expression profile of the bZIP transcription factor gene family in grapevine (Vitis vinifera)," BMC Genomics, vol. 15, no. 1, p. 1, 2014.

[12] M. Abe, Y. Kobayashi, S. Yamamoto et al., "FD, a bZIP protein mediating signals from the floral pathway integrator FT at the shoot apex," Science, vol. 309, no. 5737, pp. 1052-1056, 2005.

[13] A. B. Silveira, L. Gauer, J. P. Tomaz, P. R. Cardoso, S. Carmello-Guerreiro, and M. Vincentz, "The Arabidopsis AtbZIP9 protein fused to the VP16 transcriptional activation domain alters leaf and vascular development," Plant Science, vol. 172, no. 6, pp. 1148-1156, 2007.

[14] Y. Guan, H. Ren, H. Xie, Z. Ma, and F. Chen, "Identification and characterization of bZIP-type transcription factors involved in carrot (Daucus carota L.) somatic embryogenesis," The Plant Journal, vol. 60, no. 2, pp. 207-217, 2009.

[15] P. Lara, L. Onate-Sanchez, Z. Abraham et al., "Synergistic activation of seed storage protein gene expression in Arabidopsis by ABI3 and two bZIPs related to OPAQUE2," Journal of Biological Chemistry, vol. 278, no. 23, pp. 21003-21011, 2003.

[16] Y. Yin, Q. Zhu, S. Dai, C. Lamb, and R. N. Beachy, "RF2a, a bZIP transcriptional activator of the phloem-specific rice tungro bacilliform virus promoter, functions in vascular development," The EMBO Journal, vol. 16, no. 17, pp. 5247-5259, 1997.

[17] J. Fukazawa, T. Sakai, S. Ishida, I. Yamaguchi, Y. Kamiya, and Y. Takahashi, "Repression of shoot growth, a bZIP transcriptional activator, regulates cell elongation by controlling the level of gibberellins," The Plant Cell, vol. 12, no. 6, pp. 901-915, 2000.

[18] P. Ciceri, F. Locatelli, A. Genga, A. Viotti, and R. J. Schmidt, "The activity of the maize Opaque2 transcriptional activator is regulated diurnally," Plant Physiology, vol. 121, no. 4, pp. 1321-1327, 1999.

[19] F. Weltmeier, A. Ehlert, C. S. Mayer et al., "Combinatorial control of Arabidopsis proline dehydrogenase transcription by specific heterodimerisation of bZIP transcription factors," The EMBO Journal, vol. 25, no. 13, pp. 3133-3143, 2006.

[20] E. Baena-Gonzalez, F. Rolland, J. M. Thevelein, and J. Sheen, "A central integrator of transcription networks in plant stress and energy signalling," Nature, vol. 448, no. 7156, pp. 938-942, 2007.

[21] Y. Xiang, N. Tang, H. Du, H. Ye, and L. Xiong, "Characterization of OsbZIP23 as a key player of the basic leucine zipper transcription factor family for conferring abscisic acid sensitivity and salinity and drought tolerance in rice," Plant Physiology, vol. 148, no. 4, pp. 1938-1952, 2008.

[22] G. Lu, C. Gao, X. Zheng, and B. Han, "Identification of OsbZIP72 as a positive regulator of ABA response and drought tolerance in rice," Planta, vol. 229, no. 3, pp. 605-615, 2009.

[23] A. Balazs, I. Pocsi, Z. Hamari et al., "AtfA bZIP-type transcription factor regulates oxidative and osmotic stress responses in Aspergillus nidulans," Molecular Genetics and Genomics, vol. 283, no. 3, pp. 289-303, 2010.

[24] S. Ying, D.-F. Zhang, J. Fu et al., "Cloning and characterization of a maize bZIP transcription factor, ZmbZIP72, confers drought and salt tolerance in transgenic Arabidopsis," Planta, vol. 235, no. 2, pp. 253-266, 2012.

[25] C. Liu, Y. Wu, and X. Wang, "bZIP transcription factor OsbZIP52/RISBZ5: a potential negative regulator of cold and drought stress response in rice," Planta, vol. 235, no. 6, pp. 1157-1169, 2012.

[26] Y. Uno, T. Furihata, H. Abe, R. Yoshida, K. Shinozaki, and K. Yamaguchi-Shinozaki, "Arabidopsis basic leucine zipper transcription factors involved in an abscisic acid-dependent signal transduction pathway under drought and high-salinity conditions," Proceedings of the National Academy of Sciences, vol. 97, no. 21, pp. 11632-11637, 2000.

[27] C. Huang, J. Zhou, Y. Jie et al., "A ramie (Boehmeria nivea) bZIP transcription factor BnbZIP3 positively regulates drought, salinity and heavy metal tolerance," Molecular Breeding, vol. 36, no. 8, p. 120, 2016.

[28] H. Shimizu, K. Sato, T. Berberich et al., "LIP19, a basic region leucine zipper protein, is a Fos-like molecular switch in the cold signaling of rice plants," Plant and Cell Physiology, vol. 46, no. 10, pp. 1623-1634, 2005.

[29] C. Thurow, A. Schiermeyer, S. Krawczyk, T. Butterbrodt, K. Nickolov, and C. Gatz, "Tobacco bZIP transcription factor TGA2. 2 and related factor TGA2. 1 have distinct roles in plant defense responses and plant development," The Plant Journal, vol. 44, no. 1, pp. 100-113, 2005.

[30] H. Kaminaka, C. Nake, P. Epple et al., "bZIP10-LSD1 antagonism modulates basal defense and cell death in Arabidopsis following infection," The EMBO Journal, vol. 25, no. 18, pp. 4400-4411,2006.

[31] R. Ulm, A. Baumann, A. Oravecz et al., "Genome-wide analysis of gene expression reveals function of the bZIP transcription factor HY5 in the UV-B response of Arabidopsis," Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 5, pp. 1397-1402, 2004.

[32] M. T. Osterlund, C. S. Hardtke, N. Wei, and X. W. Deng, "Targeted destabilization of HY5 during light-regulated development of Arabidopsis," Nature, vol. 405, no. 6785, pp. 462-466, 2000.

[33] X.-L. Wang, Y. Zhong, Z.-M. Cheng, and J.-S. Xiong, "Divergence of the bZIP Gene family in strawberry, peach, and apple suggests multiple modes of gene evolution after duplication," International Journal of Genomics, vol. 2015, Article ID 536943, p. 11, 2015.

[34] V. Shulaev, D. J. Sargent, R. N. Crowhurst et al., "The genome of woodland strawberry (Fragaria vesca)," Nature Genetics, vol. 43, no. 2, pp. 109-116, 2011.

[35] K. Tamura, G. Stecher, D. Peterson, A. Filipski, and S. Kumar, "MEGA6: molecular evolutionary genetics analysis version 6.0," Molecular Biology and Evolution, vol. 30, no. 12, pp. 2725-2729, 2013.

[36] A. McLachlan and M. Stewart, "Tropomyosin coiled-coil interactions: evidence for an unstaggered structure," Journal of Molecular Biology, vol. 98, no. 2, pp. 293-304, 1975.

[37] C. D. Deppmann, A. Acharya, V. Rishi et al., "Dimerization specificity of all 67 B-ZIP motifs in Arabidopsis thaliana: a comparison to Homo sapiens B-ZIP motifs," Nucleic Acids Research, vol. 32, no. 11, pp. 3435-3445, 2004.

[38] W. Wei, Y. Hu, Y.-T. Han, K. Zhang, F.-L. Zhao, and J.-Y. Feng, "The WRKY transcription factors in the diploid woodland strawberry Fragaria vesca: identification and expression analysis under biotic and abiotic stresses," Plant Physiology and Biochemistry, vol. 105, pp. 129-144, 2016.

[39] K. J. Livak and T. D. Schmittgen, "Analysis of relative gene expression data using real-time quantitative PCR and the 2-[DELTA][DELTA]CT method," Methods, vol. 25, no. 4, pp. 402-408, 2001.

[40] A. I. Saeed, N. K. Bhagabati, J. C. Braisted et al., "[9] TM4 microarray software suite," Methods in Enzymology, vol. 411, pp. 134-193, 2006.

[41] K. F. Wei, L. J. Wu, J. Chen, C. Yf, and X. Dx, "Structural evolution and functional diversification analyses of argonaute protein," Journal of Cellular Biochemistry, vol. 113, no. 8, pp. 2576-2585, 2012b.

[42] S. Schwartz, E. Meshorer, and G. Ast, "Chromatin organization marks exon-intron structure," Nature Structural & Molecular Biology, vol. 16, no. 9, pp. 990-995, 2009.

[43] M. J. Betts, R. Guigo, P. Agarwal, and R. B. Russell, "Exon structure conservation despite low sequence similarity: a relic of dramatic events in evolution?" The EMBO Journal, vol. 20, no. 19, pp. 5354-5360, 2001.

[44] M.-Y. Sun, X.-L. Fu, Q.-P. Tan et al., "Analysis of basic leucine zipper genes and their expression during bud dormancy in peach (Prunus persica)," Plant Physiology and Biochemistry, vol. 104, pp. 54-70, 2016.

[45] H. Lin, W. Zhu, J. C. Silva, X. Gu, and C. R. Buell, "Intron gain and loss in segmentally duplicated genes in rice," Genome Biology, vol. 7, no. 5, p. 1, 2006.

[46] K.-D. Scharf, H. Heider, I. Hohfeld, R. Lyck, E. Schmidt, and L. Nover, "The tomato Hsf system: HsfA2 needs interaction with HsfA1 for efficient nuclear import and may be localized in cytoplasmic heat stress granules," Molecular and Cellular Biology, vol. 18, no. 4, pp. 2240-2251, 1998.

[47] P. Doring, E. Treuter, C. Kistner, R. Lyck, A. Chen, and L. Nover, "The role of AHA motifs in the activator function of tomato heat stress transcription factors HsfA1 and HsfA2," The Plant Cell, vol. 12, no. 2, pp. 265-278, 2000.

[48] D. Heerklotz, P. Doring, F. Bonzelius, S. Winkelhaus, and L. Nover, "The balance of nuclear import and export determines the intracellular distribution and function of tomato heat stress transcription factor HsfA2," Molecular and Cellular Biology, vol. 21, no. 5, pp. 1759-1768, 2001.

[49] L. Bentsink, J. Jowett, C. J. Hanhart, and M. Koornneef, "Cloning of DOG1, a quantitative trait locus controlling seed dormancy in Arabidopsis," Proceedings of the National Academy of Sciences, vol. 103, no. 45, pp. 17042-17047, 2006.

[50] Y. Siberil, P. Doireau, and P. Gantet, "Plant bZIP G-box binding factors," European Journal of Biochemistry, vol. 268, no. 22, pp. 5655-5666, 2001.

[51] I. Meier and W. Gruissem, "Novel conserved sequence motifs in plant G-box binding proteins and implications for interactive domains," Nucleic Acids Research, vol. 22, no. 3, pp. 470-478, 1994.

[52] L. A. Pinna, "Casein kinase 2: an 'eminence grise' in cellular regulation?" Biochimica et Biophysica Acta (BBA)-Molecular Cell Research, vol. 1054, no. 3, pp. 267-284, 1990.

[53] C. H-i, H. J-h, H. J-o, K. J-y, and S. Y. Kim, "ABFs, a family of ABA-responsive element binding factors," Journal of Biological Chemistry, vol. 275, no. 3, pp. 1723-1730, 2000.

[54] A. Emili, J. Greenblatt, and C. J. Ingles, "Species-specific interaction of the glutamine-rich activation domains of Sp1 with the TATA box-binding protein," Molecular and Cellular Biology, vol. 14, no. 3, pp. 1582-1593, 1994.

[55] M. Suckow, B. von Wilcken-Bergmann, and B. Muller-Hill, "Identification of three residues in the basic regions of the bZIP proteins GCN4, C/EBP and TAF-1 that are involved in specific DNA binding," The EMBO Journal, vol. 12, no. 3, p. 1193, 1993.

[56] X. Niu, L. Renshaw-Gegg, L. Miller, and M. J. Guiltinan, "Bipartite determinants of DNA-binding specificity of plant basic leucine zipper proteins," Plant Molecular Biology, vol. 41, no. 1, pp. 1-13, 1999.

[57] M. Suckow, K. Schwamborn, B. Kisters-Woike, B. von Wilcken-Bergmann, and B. Muller-Hill, "Replacement of invariant bzip residuces within the basic region of the yeast transcriptional activator GCN4 can change its DNA binding specificity," Nucleic Acids Research, vol. 22, no. 21, pp. 4395-4404, 1994.

[58] E. K. O'Shea, R. Rutkowski, and P. S. Kim, "Mechanism of specificity in the Fos-Jun oncoprotein heterodimer," Cell, vol. 68, no. 4, pp. 699-708, 1992.

[59] A. Acharya, S. B. Ruvinov, J. Gal, J. R. Moll, and C. Vinson, "A heterodimerizing leucine zipper coiled coil system for examining the specificity of a position interactions: amino acids I, V, L, N, A, and K," Biochemistry, vol. 41, no. 48, pp. 14122-14131, 2002.

[60] J. Moitra, L. Szilak, D. Krylov, and C. Vinson, "Leucine is the most stabilizing aliphatic amino acid in the d position of a dimeric leucine zipper coiled coil," Biochemistry, vol. 36, no. 41, pp. 12567-12573, 1997.

[61] D. Krylov, I. Mikhailenko, and C. Vinson, "A thermodynamic scale for leucine zipper stability and dimerization specificity: e and g interhelical interactions," The EMBO Journal, vol. 13, no. 12, p. 2849, 1994.

[62] C. Vinson, T. Hai, and S. Boyd, "Dimerization specificity of the leucine zipper-containing bZIP motif on DNA binding: prediction and rational design," Genes & Development, vol. 7, no. 6, pp. 1047-1058, 1993.

[63] K. Schutze, K. Harter, and C. Chaban, "Post-translational regulation of plant bZIP factors," Trends in Plant Science, vol. 13, no. 5, pp. 247-255, 2008.

[64] D. Li, F. Fu, H. Zhang, and F. Song, "Genome-wide systematic characterization of the bZIP transcriptional factor family in tomato (Solanum lycopersicum L.)," BMC Genomics, vol. 16, no. 1, p. 771, 2015.

[65] J. Zhao, R. Guo, C. Guo, H. Hou, X. Wang, and H. Gao, "Evolutionary and expression analyses of the apple basic leucine zipper transcription factor family," Frontiers in Plant Science, vol. 7, p. 376, 2016.

[66] L. G. G. Correa, D. M. Riano-Pachon, C. G. Schrago, R. V. dos Santos, B. Mueller-Roeber, and M. Vincentz, "The role of bZIP transcription factors in green plant evolution: adaptive features emerging from four founder genes," PLoS One, vol. 3, no. 8, article e2944, 2008.

[67] T. Gu, Y. Han, R. Huang, R. J. McAvoy, and Y. Li, "Identification and characterization of histone lysine methylation modifiers in Fragaria vesca," Scientific Reports, vol. 6, p. 23581, 2016.

[68] M. Gao, H. Zhang, C. Guo et al., "Evolutionary and expression analyses of basic zipper transcription factors in the highly homozygous model grape PN40024 (Vitis vinifera L.)," Plant Molecular Biology Reporter, vol. 32, no. 5, pp. 1085-1102, 2014.

[69] J. Wang, J. Zhou, B. Zhang, J. Vanitha, S. Ramachandran, and S. Y. Jiang, "Genome-wide expansion and expression divergence of the basic leucine zipper transcription factors in higher plants with an emphasis on SorghumF," Journal of Integrative Plant Biology, vol. 53, no. 3, pp. 212-231, 2011.

[70] S.-J. Oh, S. I. Song, Y. S. Kim et al., "Arabidopsis CBF3/ DREB1A and ABF3 in transgenic rice increased tolerance to abiotic stress without stunting growth," Plant Physiology, vol. 138, no. 1, pp. 341-351, 2005.

[71] K. J-y, C. H-i, I. M-y, and S. Y. Kim, "Arabidopsis basic leucine zipper proteins that mediate stress-responsive abscisic acid signaling," The Plant Cell, vol. 14, no. 2, pp. 343-357, 2002.

Xiao-Long Wang, (1,2) Xinlu Chen, (2) Tian-Bao Yang, (3) Qunkang Cheng, (4) and Zong-Ming Cheng (1,2)

(1) College of Horticulture, Nanjing Agricultural University, Nanjing 210095, China

(2) Department of Plant Sciences, University of Tennessee, Knoxville, TN 37996-4560, USA

(3) Food Quality Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, United States Department of Agriculture, Beltsville, MD 20705, USA

(4) Department of Entomology and Plant Pathology, University of Tennessee, Knoxville, TN 37996-4560, USA

Correspondence should be addressed to Zong-Ming Cheng;

Received 11 October 2016; Revised 24 January 2017; Accepted 12 February 2017; Published 11 April 2017

Academic Editor: Ferenc Olasz

Caption: Figure 1: Exon/intron organization of FvbZIP genes was depicted for each group. The exons and introns are represented by green boxes and black lines, respectively. The red boxes denote the basic and hinge region of bZIP domain. The numbers above line denote the splicing phase. The numbers "0" and "2" denote different splicing phases, "0" means splicing occurred after the third nucleotide of the codon, and "2" means splicing occurred after the second nucleotide.

Caption: Figure 2: Intron patterns (a-e) within the basic and hinge regions of bZIP proteins. P0 and P2 stand for the intron splicing phase. P0 means splicing occurred after the third nucleotide of the codon, and P2 means splicing occurred after the second nucleotide. The black bars represent the sequence of the basic and hinge regions. The vertical lines denote the positions of intron splicing phases.

Caption: Figure 3: Summery for the distribution of conserved motifs identified from FvbZIP proteins by each clade given separately. Each motif is represented by a number in colored box. See Supplementary Table S1 for detailed motif information.

Caption: Figure 4: Prediction of the dimerization properties of FvbZIP proteins. (a) Pie charts depicting the frequency of different amino acids at the a, d, e, and g positions of the leucine zippers of all FvbZIP proteins (see Supplementary Figure S3 for the amino acid positions within the leucine zipper regions). (b) Histogram of the frequency of the Asn residues present at the a position of each heptad within leucine zipper for all bZIP proteins. (c) Histogram of the frequency of attractive or repulsive g [left and right arrow] e' pairs per heptad within leucine zipper for all FvbZIP proteins.

Caption: Figure 5: Heatmap of all FvbZIP genes expressed in different tissues (a) and red/yellow fruits (b). This heatmap was generated based on the RT-PCR data (normalized) using the MeV software. Green indicates low expression, black indicates intermediate expression, and red indicates high expression.

Caption: Figure 6: Heatmap of clade A FvbZIP genes expressed under drought (a) and heat (b) treatments. This heatmap was generated based on the RT-PCR data (normalized) using the MeV software. Green indicates low expression, black indicates intermediate expression, and red indicates high expression.
Table 1: bZIP genes identified in strawberry and their detailed

                                           Number of
               Accession                     amino      Molecular
ID              number        Chromosome     acids       weight

mrna00393   XP_011468749.1       LG7          320        35828.4
mrna00517   XP_004309443.1       LG6          536        60527.7
mrna01680   XP_004290316.1       LG2          585        64330.9
mrna02177   XP_004289625.1       LG2          304        34175.6
mrna02284   XP_004300452.2       LG5          148        17187.2
mrna02614   XP_004292766.1       LG2          158        18178.4
mrna03633   XP_004299018.1       LG4          711        76012.1
mrna03778   XP_004297132.1       LG4          1106      122711.1
mrna04187   XP_004304176.1       LG6          219        25024.5
mrna04504   XP_004298330.2       LG4          585        64648.6
mrna07554   XP_004289872.1    Unanchored      270         29591
mrna07844   XP_011460972.1       LG3          324        36127.8
mrna08154   XP_004290511.1       LG2          500        53629.2
mrna08186   XP_004289896.1       LG2          425        46306.8
mrna08484   XP_004291397.1       LG2          534        57868.1
mrna08566   XP_004291338.1       LG2          216        23379.1
mrna08757   XP_004291263.1       LG2          366        39534.3
mrna09110   XP_004291101.1       LG2          487        52096.4
mrna11666   XP_004291469.1    Unanchored      166        18034.9
mrna11837   XP_004301547.1       LG5          277        29886.4
mrna11979   XP_004303083.1       LG6          535        58896.3
mrna13716   XP_011465109.1       LG5          889        96685.2
mrna14220   XP_004307087.1       LG7          334        37188.7
mrna14556   XP_004289074.1       LG1          227        25468.2
mrna14942   XP_004287866.1       LG1          157         17845
mrna15193   XP_004291844.1       LG2          171        19140.1
mrna16561   XP_004302051.1       LG1          426         46739
mrna17796   XP_004303124.1       LG6          329        37026.8
mrna18282   XP_011467186.1       LG6          162        18593.9
mrna18928   XP_004306702.1       LG7          427        45946.1
mrna21344   XP_004307733.1       LG7          364        40093.2
mrna21797   XP_004291604.1    Unanchored      543         60038
mrna21832   XP_004289640.1       LG2          159        17974.9
mrna21882   XP_011467317.1       LG6          631         70226
mrna22776   XP_011462861.1       LG4          194        22165.9
mrna23487   XP_004309615.1       LG4          445        48838.1
mrna26148   XP_004301090.1       LG5          114        13507.3
mrna27194   XP_011460752.1       LG3          848        93710.6
mrna28103   XP_004296017.2       LG3          426        46283.3
mrna28250   XP_004294799.1       LG3          432        46987.6
mrna29159   XP_004300086.1       LG5          320        35003.2
mrna29546   XP_004300611.1    Unanchored      734        80035.7
mrna30252   XP_011466783.1       LG6          532        58072.2
mrna30280   XP_011466768.1       LG6          437        47924.5
mrna31321   XP_004301946.1       LG5          471        52672.4
mrna31322   XP_011465538.1       LG5          351        39638.4
mrna31621   XP_011466109.1       LG1          362        40736.2
mrna32022   XP_011463954.1       LG5          262        28695.4
mrna32024   XP_004298867.1       LG5          272        29939.3
mrna32629   XP_004292529.1       LG2          336        37040.8

            Theoretical   Instability   Aliphatic
ID              pI           index        index     GRAVY

mrna00393      8.79          55.97        67.66     -0.853
mrna00517      6.67          67.1         68.51     -0.76
mrna01680      5.96          62.76        54.74     -0.909
mrna02177      5.05          71.42        66.38     -0.994
mrna02284      5.49          60.81        58.65     -0.943
mrna02614       8.9          68.19        67.34     -0.852
mrna03633      6.28          52.99        71.31     -0.488
mrna03778      7.64          43.21        89.76     -0.056
mrna04187      5.63          64.63        65.89     -0.886
mrna04504      6.74          44.74        78.62     -0.458
mrna07554      6.01          33.29        57.48     -0.776
mrna07844      6.36          45.93        74.1      -0.667
mrna08154      6.25          54.88        69.92     -0.41
mrna08186      6.07          65.47        70.52     -0.723
mrna08484      5.95          60.87        61.22     -0.753
mrna08566      9.48          55.37        62.96     -0.634
mrna08757      5.49          46.89        64.56     -0.655
mrna09110      9.32          44.04        62.05     -0.64
mrna11666      9.56          68.79        70.6      -0.965
mrna11837      6.21          52.67        70.83     -0.739
mrna11979      8.66          51.21        88.77     -0.393
mrna13716      5.97          55.69        72.56     -0.57
mrna14220      7.83          55.28        81.92     -0.558
mrna14556      7.15          71.85        61.5      -0.899
mrna14942      4.93          44.46        79.49     -0.58
mrna15193      6.51          58.36        71.4      -0.794
mrna16561      5.18          63.18        57.98     -0.895
mrna17796      5.85          60.67        62.34     -1.004
mrna18282      7.08          60.44        84.81     -0.733
mrna18928      8.79          55.68        67.73     -0.725
mrna21344      5.91          59.84        59.84     -0.881
mrna21797      6.52          51.39        72.14     -0.545
mrna21832      6.29          66.22        66.98     -0.701
mrna21882      6.66          52.69        78.56     -0.428
mrna22776      5.34          53.68        54.74     -1.346
mrna23487      6.08          62.25        61.08     -0.783
mrna26148        9           46.47        94.04     -0.712
mrna27194       6.1          61.08        68.82     -0.549
mrna28103      6.35          59.02        56.92     -0.754
mrna28250      9.54          49.8         65.46     -0.663
mrna29159      4.73          59.96        80.22     -0.443
mrna29546      5.74          49.93        64.75     -0.662
mrna30252      7.49           61          67.74     -0.619
mrna30280      8.45          62.09        60.09     -0.778
mrna31321      8.48          47.98        67.94     -0.652
mrna31322      6.12          58.95        77.92     -0.572
mrna31621      6.35          43.26        79.59     -0.439
mrna32022       4.8          44.98        68.93     -0.842
mrna32024      4.97          47.36        82.1      -0.544
mrna32629      6.64          56.65        61.64     -0.838
COPYRIGHT 2017 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Wang, Xiao-Long; Chen, Xinlu; Yang, Tian-Bao; Cheng, Qunkang; Cheng, Zong-Ming
Publication:International Journal of Genomics
Article Type:Report
Date:Jan 1, 2017
Previous Article:Discovery of MicroRNAs and Their Target Genes Related to Drought in Paulownia "Yuza 1" by High-Throughput Sequencing.
Next Article:Evidence of Divergent Amino Acid Usage in Comparative Analyses of R5- and X4-Associated HIV-1 Vpr Sequences.

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters