- Research article
- Open Access
Basic leucine zipper (bZIP) transcription factor genes and their responses to drought stress in ginseng, Panax ginseng C.A. Meyer
BMC Genomics volume 22, Article number: 316 (2021)
Ginseng is an important medicinal herb in Asia and Northern America. The basic leucine zipper (bZIP) transcription factor genes play important roles in many biological processes and plant responses to abiotic and biotic stresses, such as drought stress. Nevertheless, the genes remain unknown in ginseng.
Here, we report 91 bZIP genes identified from ginseng, designated PgbZIP genes. These PgbZIP genes were alternatively spliced into 273 transcripts. Phylogenetic analysis grouped the PgbZIP genes into ten groups, including A, B, C, D, E, F, G, H, I and S. Gene Ontology (GO) categorized the PgbZIP genes into five functional subcategories, suggesting that they have diversified in functionality, even though their putative proteins share a number of conserved motifs. These 273 PgbZIP transcripts expressed differentially across 14 tissues, the roots of different ages and the roots of different genotypes. However, the transcripts of the genes expressed coordinately and were more likely to form a co-expression network. Furthermore, we studied the responses of the PgbZIP genes to drought stress in ginseng using a random selection of five PgbZIP genes, including PgbZIP25, PgbZIP38, PgbZIP39, PgbZIP53 and PgbZIP54. The results showed that all five PgbZIP genes responded to drought stress in ginseng, indicating that the PgbZIP genes play important roles in ginseng responses to drought stress.
These results provide knowledge and gene resources for deeper functional analysis of the PgbZIP genes and molecular tools for enhanced drought tolerance breeding in ginseng.
Ginseng (Panax ginseng C.A. Meyer) is an important medicinal herb in Asia and Northern America. In China, ginseng has a long cultivation history and is mainly cultivated in Jilin Province where it is known as Jilin ginseng. Ginsenosides, present in most tissues of ginseng, are recognized as the most valuable active components of ginseng [1, 2]. Ginseng has a lot of benefits for human, such as relieving pain, improving brain function, and increasing anti-tumor activity [3,4,5]. However, ginseng is frequently suffering from different biotic and abiotic stresses, including, but not limited to, diseases, insect pests, drought, cold, heat, and daylight intensity, greatly threatening its production. Therefore, it is necessary to comprehensively investigate the genes involved in plant defense to the stresses in ginseng.
Transcription factors (TFs) have been shown to play a vital role in plant responses to various biotic or abiotic stresses. The basic leucine zipper (bZIP) transcription factor containing a conserved bZIP domain that is composed of 60–80 amino acids is known as one of the largest TF families [6, 7]. The conserved bZIP domain is composed of two important functional regions: the basic region and the leucine zipper region, linked by one hinge [8, 9]. The basic region usually contains an invariant N-x7-R/K motif (approximately 16 amino acids) and is responsible for both nuclear localization and DNA binding. The leucine zipper region mediates the homo- and/or hetero-dimerization as it contains a less conserved dimerization motif [10,11,12,13].
The bZIP genes have been documented to play a vital role in a number of biological processes, including plant tissue and organ differentiation and vascular development [14, 15], embryogenesis , and seed maturation . Studies have also shown that the bZIP genes code key components in plant regulation of biotic and abiotic stresses, e.g., pathogens [18, 19], osmosis [20, 21], salinity [22, 23], cold [13, 24], and drought [25, 26]. It has been reported that AtbZIP28 was activated by thermal stress, and then regulated the expressions of heat-responsive genes to protect plants from heat stress . In rice, OsbZIP23 and OsbZIP72 were reported to attenuate drought stress by activating ABA signaling [28, 29]. Knockouting SlbZIP1 and SlAREB1 that belong to Group A of the SlbZIP gene family increased salt stress tolerance, while the over-expressions of SlbZIP1 and SlAREB1 decreased salt stress tolerance in tomato [30, 31].
The bZIP gene family has been analyzed in several plant species and shown to vary in size. For example, 75 bZIP genes were identified in Arabidopsis , 69 in tomato , and 89 in rice . However, no research on the bZIP gene family has been reported yet in ginseng. The present study first identified the 91 bZIP genes from ginseng, which were designated PgbZIP genes. We then examined their conserved protein motifs, phylogeny, putative functionality, and expression characteristics and co-expression networks in different tissues, different year-old roots, and in the roots of different genotypes. Because drought stress restricts plant growth and development , influencing ginseng production, PgbZIP genes were further studied in response to drought stress in ginseng.
Identification of PgbZIP genes
A total of 1,957 transcript sequences containing the bZIP domain were identified from Database A consisting of 248,993 transcripts. As the conserved domains of 1,684 of the 1,957 transcripts were incomplete or out of ORFs (open-reading frames), the remaining 273 transcripts that contain complete bZIP domains in their ORFs were identified as the PgbZIP gene transcripts for Jilin ginseng. These 273 transcripts had a sequence length of 210 to 3,651 bp, with an average length of 1,449 bp (Additional file 1: Table S1). Analysis showed that these 273 transcripts were alternatively spliced from 91 PgbZIP genes . Of these 273 PgbZIP transcripts, 190 contained full-length ORFs that were derived from 62 PgbZIP genes. The full-length proteins encoded by the 190 PgbZIP transcripts contained amino acids varying from 46 (PgbZIP84) to 785 (PgbZIP63-1), with an average of 294 amino acids (Additional file 2: Table S2). In comparison, 45 (49 %) of the 91 Jilin ginseng PgbZIP genes identified in this study were orthologous to 111 (76 %) of the 146 Korean ginseng PgbZIP genes. The remaining 46 (51 %) of the Jilin ginseng PgbZIP genes were newly discovered or Jilin ginseng-specific (Additional file 3: Table S3; Additional file 4: Fig. S1).
Phylogeny and conserved motifs of the PgbZIP gene family
The longest transcript for each of the 62 genes containing a full-length ORF was used to construct the NJ (neighbor-joining) phylogenetic tree of the PgbZIP gene family. Fifty-six bZIP genes that were identified from Arabidopsis (20 AtbZIP genes), rice (19 OsbZIP genes) and tomato (17 SlbZIP genes) (Additional file 5: Table S4) were used as outgroups. The 62 PgbZIP genes were clustered into ten clades, defined ten groups in this study, with the AtbZIP, OsbZIP and SlbZIP genes from Arabidopsis, rice and tomato (Fig. 1a). This result suggested that the PgbZIP gene family is an ancient gene family that originated before splitting between the monocot (rice) and dicot (Arabidopsis and tomato) plants. The PgbZIP gene family has the same number of groups as the AtbZIP, OsbZIP or SlbZIP gene family , but consists of more groups than the bZIP gene family of castor bean , cucumber  or sorghum . Group A of the PgbZIP gene family has the largest number of PgbZIP genes, with 11 PgbZIP genes and Group H has only one PgbZIP gene (PgbZIP13). Similarly, we also constructed the MP (maximum parsimony) tree for the PgbZIP gene family (Additional file 6: Fig. S2). The MP tree was essentially the same as the NJ tree, with a difference from the NJ tree in grouping of only PgbZIP08, PgbZIP09 and PgbZIP84, which was likely due to their low bootstrap confidences for both the NJ and MP trees.
Twenty conserved motifs were identified from the putative proteins encoded by the 62 PgbZIP genes that were spliced into transcripts with full-length ORFs (Additional file 2: Table S2). The distribution of these 20 conserved motifs in the 62 PgbZIP genes is shown in Fig. 1b. Motif 1, annotated as the bZIP domain, is presented in all 62 PgbZIP genes of the gene family. The putative proteins of most PgbZIP genes in a group of the gene family usually have a similar set of motifs (Fig. 1b). For instance, the putative proteins of most PgbZIP genes of Group A contains Motifs 3, 12, 15, 16 and 17; those of Group D harbor Motifs 2, 4, 6, 7, 8, 13 and 20; and those of Group G share Motif 11. Moreover, some groups of the PgbZIP gene family also share the same motif. For instance, Groups G and S both possess Motif 10, and Groups F and G are common in Motif 14. These results suggested the similarities of the PgbZIP genes in functionality.
Functional differentiation of the PgbZIP gene family
We examined the functional differentiation of the PgbZIP gene family by categorizing the 273 PgbZIP gene transcripts using Gene Ontology (GO). Two hundred fifty-one (91.9 %) of the 273 PgbZIP gene transcripts were annotated and categorized into all three primary categories, Biological Process (BP), Molecular Function (MF) and Cellular Component (CC) (Fig. 2a; Additional file 7: Table S5). BP contained 235 PgbZIP transcripts, MP contained 249 PgbZIP transcripts and CC had two PgbZIP transcripts (Fig. 2a). At Level 2, these 251 PgbZIP gene transcripts were categorized into five subcategories, including two BP subcategories (transcription DNA-templated and regulation of gene expression), two MF subcategories (nucleic acid binding transcription factor activity and DNA binding), and one CC subcategory (cytosol) (Fig. 2b). Of these 5 subcategories, all except cytosol were enriched in number of PgbZIP transcripts (P ≤ 0.01).
The PgbZIP gene transcripts expressed in 14 tissues (Fig. 3a) of a four-year-old plant, the four-year-old roots of 42 genotypes (Fig. 3b), and the roots of four differently aged plants (Fig. 3c) were also categorized, respectively. These PgbZIP transcripts were also categorized into the above mentioned five subcategories, but the numbers of the PgbZIP gene transcripts categorized into the five subcategories varied substantially across tissues, genotypes and developmental stages. These results together demonstrated the functional differentiation of the PgbZIP gene family and also confirmed their functional consistency across tissues, genotypes and developmental stages.
Expression characteristics of the PgbZIP gene transcripts
To characterize the expressions of the PgbZIP genes, the 273 PgbZIP transcripts were investigated in expression in different tissues, different year-old plant roots, and the four-year-old roots of different genotypes (Additional file 8: Table S6). The PgbZIP transcripts were used for this experiment because different transcripts spliced from one gene may have different biological functions . The expressions of the PgbZIP transcripts varied dramatically, ranging from 0.0 TPM to 307.7 TPM, to 178.5 TPM, and to 169.9 TPM among different tissues, different year-old plant roots, and different genotypes, respectively.
For different tissues of a four-year-old plant, 248 (91 %) of the 273 PgbZIP transcripts expressed in at least one tissue (relative to the reference transcriptome), but most (77 %) expressed in two or more tissues. Sixty-eight (25 %) PgbZIP transcripts expressed in all 14 tissues and 37 PgbZIP transcripts (14 %) showed tissue-specific expressions (Fig. 4a). For different developmental stages, 167 (61 %) of the 273 PgbZIP transcripts expressed at least in one of 5-, 12-, 18- and 25-year-old plant roots, 93 (34 %) expressed at all four developmental stages, 35 (13 %) were developmental stage expression-specific, and 39 (14 %) expressed at two or three of the developmental stages (Fig. 4b). For the four-year-old roots of different genotypes, 208 (76 %) of the 273 PgbZIP transcripts expressed in at least one genotype, 55 (20 %) expressed in all 42 genotypes, 12 (4 %) were genotype expression-specific, and 141 (52 %) expressed in 2–41 genotypes (Fig. 4c). For individual tissues, root developmental stages or genotypes, from 40 − 60 % of the 273 PgbZIP transcripts expressed (Additional file 9: Fig. S3).
Furthermore, gene expression heatmaps were constructed for the 273 PgbZIP transcripts that expressed in different tissues (Fig. 5a), different year-old plant roots (Fig. 5b) and the four-year-old roots of different genotypes (Fig. 5c) to estimate their co-regulations and expression patterns. The results showed that a number of the PgbZIP transcripts had identical expression patterns across the 14 tissues, four different developmental stages or 42 genotypes, suggesting that they were co-regulated.
The co-expression network of the PgbZIP genes
To assess the potential functional relationships among different members of the PgbZIP gene family, the co-expression network of its 91 PgbZIP genes was constructed based on the expressions of their 273 transcripts in the four-year-old plant roots of 42 genotypes at a P ≤ 0.05 (Fig. 6). Two hundred seventy-three functionally unknown Jilin ginseng gene transcripts were randomly selected from Database A as the negative controls. Consequently, 208 of the 273 PgbZIP transcripts formed a co-expression network that was composed of 208 nodes, 1,994 edges (Fig. 6a), and 18 clusters (Fig. 6b). In comparison, the co-expression network of the PgbZIP transcripts was much more robust than that constructed from the 273 randomly-selected unknown ginseng transcripts at all significance levels from P ≤ 5.0E-02 through P ≤ 1.0E-08 (Fig. 6c,d). Statistical analysis confirmed the tendency that PgbZIP transcripts was more likely to form a co-expression network than the randomly-selected unknown ginseng transcripts (Fig. 6e,f). These results concluded that the members of the PgbZIP gene family were more likely to form a co-expression network, suggesting that they likely function correlatively.
Response of the PgbZIP gene family to drought stress
To test whether the PgbZIP gene family functions in plant response to drought stresses, five PgbZIP genes, PgbZIP25, PgbZIP38, PgbZIP39, PgbZIP53 and PgbZIP54, were randomly selected from the PgbZIP gene family and examined in plant response to drought stress. Ginseng seedlings were treated with 20 % (w/v) PEG-6000 that is widely used to drought-stress plants for 3, 6, 12, 24 and 48 h. The RWCs (relative water contents) of the seedlings treated with and without PEG-6000 were determined and compared. The RWCs of the seedlings treated with PEG-6000 for 24 h were significantly reduced, relative to the control seedlings not treated with PEG-6000. No significant difference in RWC was observed for the seedlings treated with PEG-6000 for other time points (Additional file 10: Fig. S4). Nevertheless, the expressions of all five PgbZIP genes studied were up-regulated in the seedlings treated with PEG-6000. Specifically, the expressions of PgbZIP25, PgbZIP38, PgbZIP39, PgbZIP53 and PgbZIP54 in the seedlings treated with PEG-6000 for 3 h were up-regulated by 4.1-, 12.6-, 8.6-, 4.7- and 18.9-fold, respectively, over those in the untreated seedlings (P ≤ 0.01). However, as the treatment time increased the expressions of these five genes varied differently. For instance, the expression of PgbZIP25 reached the peak after the seedlings stressed with PEG-6000 for 6 h, while no significant different expression of the gene was observed in the treated seedlings from that in the control seedlings after stressed for 12 h, 24 and 48 h. PgbZIP38 was continuously up-regulated in the seedlings treated by PEG-6000 for 3 h through 24 h, but returned to the expression level as in the control seedlings at 48 h. The expression of PgbZIP53 in the seedlings stressed with PEG-6000 for 6 h through 48 h showed no significant difference from that in the control seedlings. PgbZIP39 and PgbZIP54 showed irregular expression variation in the seedlings stressed at different time points (Fig. 7).
We have identified 91 PgbZIP genes that expressed in a four-year-old Jilin ginseng plant. The size of the PgbZIP gene family is comparable to those identified in rice (89 bZIP genes) , Brachypodium distachyon (96) , barley (89) , and sorghum (92) , but larger than those identified in Arabidopsis (75) , castor bean (49) , cucumber (64)  and grapevine (55) , and smaller than those identified in maize (125)  and soybean (160) . This result suggests that the PgbZIP gene family is an intermediate-sized transcription factor gene family. It seems independent of its genome size and perennial growth nature.
Phylogenetic analysis in the present study shows that the PgbZIP gene family is made up of 10 groups, which is the same as the bZIP gene families identified in Arabidopsis, rice, tomato and maize. This result suggests that the PgbZIP gene family has a similar evolutionary trajectory to those of Arabidopsis, rice, tomato and maize. Notably, all ten groups of the PgbZIP gene family are grouped with the bZIP genes of Arabidopsis, tomato and rice, implying that the PgbZIP gene family has existed before splitting between monocots (rice) and dicots (ginseng, Arabidopsis, tomato). The clustering of the PgbZIP gene family and its significant sharing of conserved motifs with those of Arabidopsis, rice and tomato indicate that they have a common ancestor.
It has been reported that the bZIP genes are involved in a variety of biological processes, including drought/osmotic stress response [9, 36, 37], growth and development and cell elongation [40, 41], organ and tissue differentiation , and seed storage protein gene regulation . This study categorizes the PgbZIP genes into five subcategories that belong to all three primary categories. This result suggests that the functionality of the PgbZIP gene family has been substantially differentiated. However, the functional differentiation of the gene family was far smaller than those of the PgRLK gene family (23 subcategories)  and the PgNBS gene family (36 subcategories)  identified in Jilin ginseng. Our result indicates that the PgbZIP gene family mainly functions in DNA-templated transcription, regulation of gene expression, and nucleic acid binding transcription factor activity.
The PgbZIP gene family actively expresses in all tissues, at all developmental stages and in all genotypes examined in this study, but only approximately 50 % of the genes in the family expressed in a tissue, at a developmental stage and in a genotype. Although the expression activities of the genes in the PgbZIP gene family vary dramatically across tissues, at different developmental stages and in different genotypes, most of them expressed in multiple tissues and only a small portion are tissue-, developmental stage- or genotype-specific. Expression heatmap analysis reveals that co-regulation of the PgbZIP gene expressions exists across tissues, developmental stages and genotypes, but the co-regulation is observed only for a limited number of the PgbZIP genes. The tendency of co-expression network formation for the majority of the PgbZIP genes indicates the functional correlation of the PgbZIP genes and also their functional differentiation.
Previous studies showed that the bZIP genes isolated from mung bean, adzuki bean, Arabidopsis, wheat, rice, and Tamarix hispida were involved in plant response to drought and salt stresses [45,46,47,48,49]. This study shows that all five PgbZIP genes, PgbZIP39, PgbZIP25, PgbZIP38, PgbZIP53 and PgbZIP54, randomly selected from the A, S, G, S and I groups of the PgbZIP gene family, respectively, responded to the drought stress stimulated by PEG-6000. This result confirms that the PgbZIP gene family also plays a role in plant response to drought stress. Nevertheless, additional research will be needed to further investigate the tolerance of the PgbZIP genes to drought stress and their underlying molecular mechanisms.
Ninety-one PgbZIP genes were identified from Jilin ginseng and systematically analyzed in phylogeny, conservation, functional differentiation, expression, and network interaction. The PgbZIP gene family is an ancient gene family and has substantially differentiated in functionality. The expressions of the PgbZIP genes varied spatially, temporally and across genotypes, but they were more likely to form a co-expression network, suggesting their functional correlation. It is observed that the PgbZIP gene family is involved in plant response to drought stress in ginseng. Together, the results of this study conclude that the PgbZIP gene family consists of at least 91 gene members, its functionality has diverged but its members remain functionally correlated at a substantial degree, and it likely plays a significant role in plant response to drought stress.
The Jilin Ginseng Transcriptome Databases developed from 14 tissues of a four-year-old plant (Database A), the roots of four different year-old plants (Database B) , and the four-year-old roots of 42 genotypes collected from Jilin province, China, defined from S1 to S42 (Database C) (Additional file 11: Table S7) [44, 51] were used for this study. All the samples for the databases were collected at the fruiting stage of the plants.
Identification of the bZIP gene family in ginseng
Three steps were conducted to identify the bZIP genes from Jilin ginseng. First, the local Hidden Markov Model (HMM) search was carried out with the bZIP domain HMM profiles, including PF07716, PF00170, PF03131, and PF012498 (http://pfam.sanger.ac.uk/), using the HMMER3.0 software (HMMER: http://hmmer.wustl.edu/). Then, BLAST search was performed at a threshold of 1.0E-06 to identify putative bZIP genes from Database A. Third, all putative bZIP genes were subjected to filtration by the Online Conserved Domain Search. The resulting bZIP genes were identified as PgbZIP genes that were defined from PgbZIP01 through PgbZIP91 and whose transcripts were indicated with suffix, such as “_1, _2, and so on” (Additional file 1: Table S1).
Comparison of the PgbZIP genes between Jilin ginseng and Korean ginseng
To have an overview of the Jilin ginseng PgbZIP genes with those identified from the genome database of Korean ginseng, we compared these two sets of PgbZIP genes. The sequences of the Korean ginseng PgbZIP genes were downloaded from the Korean Ginseng Genome Database (http://ginsengdb.snu.ac.kr/index.php)  and aligned with the Jilin ginseng PgbZIP genes identified from the Jilin Ginseng Transcriptome Database (see above). The sequences with an identity of ≥ 95 % and an alignment length of ≥ 240 bp (the bZIP domain maximum length) were identified as the orthologous genes of the Jilin and Korean ginsengs in this study .
Phylogenetic analysis of the PgbZIP gene family in Jilin ginseng
The Jilin ginseng PgbZIP genes were translated into putative proteins and those having complete ORFs were selected as representatives for phylogenetic analysis of the PgbZIP gene family, with one longest bZIP protein sequence per gene. The full-length protein sequences of the bZIP genes of Arabidopsis thaliana (Arabidopsis), Solanum lycopersicum (tomato) and Oryza sativa (rice) were downloaded from the Plant TF Database (PlantTFDB, v3.0) , the Sol Genomics Network (SGN; http://solgenomics.net/) and the Rice Genome Annotation Project RGAP (http://rice.plantbiology.msu.edu/) , respectively, and used as outgroups. Then, the PgbZIP-coding proteins of ginseng, along with those of Arabidopsis, rice and tomato, were subjected to multiple sequence alignment. Finally, the phylogenetic tree of the PgbZIP gene family was constructed using MEGA version 7.0 (http://www.megasoftware.net), with the Neighbor-Joining (NJ) and Maximum Parsimony (MP) algorithm, respectively and 1,000 replications.
Conserved motifs of the PgbZIP genes
The ORFs of the full-length PgbZIP transcripts were first identified using the online ORF Finder at NCBI (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). Then, the conserved motifs of the PgbZIP proteins were searched using the Motif Elicitation Tool (version 4.9.1, http://meme.sdsc.edu/meme/cgi-bin/meme.cgi)  with a maximum number of 20. The minimum and maximum lengths of the conserved motifs were 10 and 50 amino acids, respectively, and other parameters were used as default .
Functional categorization of the PgbZIP transcripts
The 273 PgbZIP transcripts were submitted to Blast2GO (Version 4.1.5) to perform GO analysis and categorization according to their GO terms. The number of transcripts categorized into each functional subcategory was subjected to enrichment analysis  using the GO functional categorization of 1,000 unknown ginseng gene transcripts randomly selected from Database A  as controls. The difference between the observed number of PgbZIP transcripts and the expected number of randomly-selected unknown ginseng transcripts categorized into each functional subcategory (Level 2) were examined by Chi-square test.
Expressions and network analysis of the PgbZIP transcripts
The expressions of the 273 PgbZIP transcripts were extracted from Databases A , B  and C , respectively. The R programming language was used to construct the expression heatmaps of the 273 PgbZIP transcripts and the BioLayout Express3D software (Version 3.2) was used to construct their co-expression networks .
Plant growth and drought stress
The seeds of Jilin ginseng, Damaya, were grown in the pots containing nutritional soil for 4 weeks. Then, the four-week-old seedlings were treated with 20 % PEG-6000 (polyethylene glycol 6000) to simulate drought stress. PEG-6000 has been widely used to mimic drought stress for study of plant response to drought stress [59,60,61]. After treated with PEG-6000 for 3 h, 6 h, 12 h 24 and 48 h, the seedlings were harvested and weighed immediately (fresh weight). The seedlings were then placed in a vial containing distilled water at 4℃ for 24 h, and the saturated weight of the seedlings was measured. Third, these seedlings were dried at 70℃ for 72 h, and their dry weights were measured. Finally, the relative water contents (RWCs) of the control and drought-stressed seedlings were measured using the following formula: [(fresh weight - dry weight) / (saturated weight - dry weight)] x 100 .
RNA isolation and analysis by quantitative real‐time PCR
After treated with PEG-6000 for 0 h, 3 h, 6 h, 12 h 24 and 48 h, the total RNAs of seedlings were isolated by TRIzol reagent (Biotake, Beijing, China) and the first-strand cDNAs were synthesized using the PrimeScript™ RT reagent Kit with gDNA Eraser (TaKaRa, Tokyo, Japan). Ubiquinol-cytochrome C reductase (QCR) gene of ginseng was selected as the reference gene . Five PgbZIP genes, PgbZIP25, PgbZIP38, PgbZIP39, PgbZIP53 and PgbZIP54, were randomly selected from the 91 PgbZIP genes and the primers specific for these genes were designed and synthesized (Additional file 12: Table S8). qRT-PCR was conducted using an Applied Biosystems 7500 Real Time PCR System (Thermo Fisher Scientific,Waltham, USA) and SYBR Premix Ex Taq™ II (TaKaRa, Tokyo, Japan). The PCR conditions were 30 s at 95℃, and followed by 40 cycles of 5 s at 95℃, 34 s at 60℃, and one cycle of 15 s at 95℃, 60 s at 60℃, and finally, 15 s at 95℃. Three biological replicates were conducted and the gene relative expression was calculated using the 2−ΔΔCT formula.
Availability of data and materials
The data used for this study are available at Sequence Read Archive (SRA) of National Center for Biotechnology Information (NCBI), BioProject PRJNA302556; and at Gene Expression Omnibus (GEO) of NCBI, SRP066368 and SRR13131364 - SRR13131405. The plant materials are available through corresponding authors, upon request.
basic leucine zipper.
- PgbZIP :
Panax ginseng bZIP
open reading frame
multiple EM for motif elicitation
Hidden Markov Model
polyethylene glycol 6000
relative water content
Choi KT. Botanical characteristics, pharmacological effects and medicinal components of Korean Panax ginseng C.A. Meyer. Acta Pharmacolo Sin. 2008;29:1109–18.
Zhang YC, Li G, Jiang C, Yang B, Yang HJ, Xu HY, et al. Tissue-specific distribution of ginsenosides in different aged ginseng and antioxidant activity of ginseng leaf. Molecules. 2014;19:17381–99.
Popovich DG, Yeo SY, Zhang W. Ginseng (Panax quinquefolius) and licorice (Glycyrrhiza uralensis) root extract combinations increase hepatocarcinoma cell (hep-g2) viability. Evid Based Complement Alternat Med. 2011;2011:408273.
Lee ST, Chu K, Sim JY, Heo JH, Kim M. Panax ginseng enhances cognitive performance in Alzheimer disease. Alzheimer Dis Assoc Disord. 2008;22:222–6.
Xie X, Wang HT, Li CL, Gao XH, Ding JL, Zhao HH, et al. Ginsenoside Rb1 protects PC12 cells against β-amyloid-induced cell injury. Mol Med Rep. 2010;3:635–9.
Alves MS, Dadalto SP, Gonçalves AB, De Souza GB, Barros VA, Fietto LG. Plant bZIP transcription factors responsive to pathogens: a review. Int J Mol Sci. 2013;14:7815–28.
Nijhawan A, Jain M, Tyagi AK, Khurana JP. Genomic survey and gene expression analysis of the basic leucine zipper transcription factor family in rice. Plant Physiol. 2008;146:333–50.
Wei K, Chen J, Wang Y, Chen Y, Chen S, Lin Y, et al. Genome-wide analysis of bZIP-encoding genes in maize. DNA Res. 2012;19:463–76.
Liu X, Chu Z. Genome-wide evolutionary characterization and analysis of bZIP transcription factors and their expression profiles in response to multiple abiotic stresses in Brachypodium distachyon. BMC Genom. 2015;16:227.
Jakoby M, Weisshaar B, Dröge-Laser W, Vicente-Carbajosa J, Tiedemann J, Kroj T, et al. bZIP transcription factors in Arabidopsis. Trends Plant Sci. 2002;7:106–11.
Foster R, Izawa T, Chua NH. Plant bZIP proteins gather at ACGT elements. Faseb J. 1994;8:192.
Pourabed E, Golmohamadi FG, Monfared PS, Razavi SM, Shobbar ZS. Basic leucine zipper family in barley: genome-wide characterization of members and expression analysis. Mol Biotechnol. 2015;57:12–26.
Liao Y, Zou HF, Wei W, Hao YJ, Tian AG, Huang J, 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. 2008;228:225–40.
Abe M, Kobayashi Y, Yamamoto S, Daimon Y, Yamaguchi A, Ikeda Y, et al. FD, a bZIP protein mediating signals from the floral pathway integrator FT at the shoot apex. Science. 2005;309:1052–6.
Silveira AB, Gauer L, Tomaz JP, Cardoso PR, Carmello-Guerreiro S, Vincentz M. The Arabidopsis AtbZIP9 protein fused to the VP16 transcriptional activation domain alters leaf and vascular development. Plant Sci. 2007;172:1148–56.
Guan Y, Ren H, Xie H, Ma Z, Chen F. Identification and characterization of bZIP-type transcription factors involved in carrot (Daucus carota L.) somatic embryogenesis. Plant J. 2010;60:207–17.
Lara P, Onate-Sanchez L, Abraham Z, Ferrandiz C, Diaz I, Carbonero P, et al. Synergistic activation of seed storage protein gene expression in Arabidopsis by ABI3 and two bZIPs related to OPAQUE2. J Biol Chem. 2003;278:21003–11.
Thurow C, Schiermeyer A, Krawczyk S, Butterbrodt T, Nickolov K, Gatz C. Tobacco bZIP transcription factor TGA2.2 and related factor TGA2.1 have distinct roles in plant defense responses and plant development. Plant J. 2010;44:100–13.
Kaminaka H, NaKe C, Epple P, Dittgen J, Schütze K, Chaban C, et al. bZIP10-LSD1 antagonism modulates basal defense and cell death in Arabidopsis following infection. Embo J. 2006;25:4400–11.
Weltmeier F, Ehlert A, Mayer CS, Dietrich K, Wang X, Schütze K, et al. Combinatorial control of Arabidopsis proline dehydrogenase transcription by specific heterodimerisation of bZIP transcription factors. EMBO J. 2014;25:3133–43.
Balázs A, Pócsi I, Hamari Z, Leiter E, Emri T, Miskei M, et al. AtfA bZIP-type transcription factor regulates oxidative and osmotic stress responses in Aspergillus nidulans. Mol Genet Genomics. 2010;283:289–303.
Wang Z, Yan L, Wan L, Huai D, Kang Y, Shi L. Genome-wide systematic characterization of bZIP transcription factors and their expression profiles during seed development and in response to salt stress in peanut. BMC Genom. 2019;20:51.
Huang C, Zhou J, Jie Y, Xing H, Zhong Y, She W, et al. A ramie (Boehmeria nivea) bZIP transcription factor BnbZIP3 positively regulates drought, salinity and heavy metal tolerance. Mol Breed. 2016;36:120.
Yang Y, Li J, Li H, Yang Y, Guang Y, Zhou Y. The bZIP gene family in watermelon: genome-wide identification and expression analysis under cold stress and root-knot nematode infection. Peer J. 2019;7:e7878.
Ying S, Zhang DF, Fu J, Shi YS, Song YC, Wang TY. Cloning and characterization of a maize bZIP transcription factor, ZmbZIP72, confers drought and salt tolerance in transgenic Arabidopsis. Planta. 2012;235:253–66.
Liu C, Wu Y, Wang X. bZIP transcription factor RISBZ5: a potential negative regulator of cold and drought stress response in rice. Planta. 2012;235:1157–69.
Gao H, Brandizzi F, Benning C, Larkin RM. A membrane-tethered transcription factor defines a branch of the heat stress response in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2008;105:16398–403.
Xiang Y, Tang N, Du H, Ye H, Xiong L. 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 Physiol. 2008;148:1938–52.
Lu G, Gao C, Zheng X, Han B. Identification of OsbZIP72 as a positive regulator of ABA response and drought tolerance in rice. Planta. 2009;229:605–15.
Orellana S, Yañez M, Espinoza A, Verdugo I, González E, Ruiz-Lara S, et al. The transcription factor SlAREB1 confers drought, salt stress tolerance and regulates biotic and abiotic stress-related genes in tomato. Plant Cell Environ. 2010;33:2191–208.
Zhu M, Meng X, Cai J, Li G, Dong T, Li Z. Basic leucine zipper transcription factor SlbZIP1 mediates salt and drought stress tolerance in tomato. BMC Plant Biol. 2018;18:83.
Li D, Fu F, Zhang H, Song F. Genome-wide systematic characterization of the bZIP transcriptional factor family in tomato ( Solanum lycopersicum L.). BMC Genom. 2015;16:771.
Godwin J, Farrona S. Plant epigenetic stress memory induced by drought: a physiological and molecular perspective. Methods Mol Biol. 2020;2093:243–59.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Jin Z, Xu W, Liu A. Genomic surveys and expression analysis of bZIP gene family in castor bean (Ricinus communis L.). Planta. 2014;239:299–312.
Baloglu MC, Eldem V, Hajyzadeh M, Unver T. Genome-wide analysis of the bZIP transcription factors in cucumber. PloS One. 2014;9:e96014.
Wang J, Zhou J, Zhang B, Vanitha J, Ramachandran S, Jiang SY. Genome-wide expansion and expression divergence of the basic leucine zipper transcription factors in higher plants with an emphasis on sorghum. J Integr Plant Biol. 2011;53:212–31.
Syed NH, Kalyna M, Marquez Y, Barta A, Brown JW. Alternative splicing in plants – coming of age. Trends Plant Sci. 2012;17:616–23.
Liu J, Chen N, Chen F, Cai B, Santo SD, Tornielli GB, et al. Genome-wide analysis and expression profile of the bZIP transcription factor gene family in grapevine (Vitis vinifera). BMC Genom. 2014;15:281–1.
Chen X, Yao Q, Gao X, Jiang C, Nicholas P, Fu X. Shoot-to-root mobile transcription factor HY5 coordinates plant carbon and nitrogen acquisition. Curr Biol. 2016;26:640–6.
Fukazawa J, Sakai T, Ishida S, Yamaguchi I, Takahashi Y. REPRESSION OF SHOOT GROWTH, a bZIP transcriptional activator, regulates cell elongation by controlling the level of gibberellins. Plant Cell. 2000;12:901–15.
Yan Q, Wu F, Ma T, Zong X, Ma Q, Li J. Comprehensive analysis of bZIP transcription factors uncovers their roles during dimorphic floret differentiation and stress response in Cleistogenes songorica. BMC Genom. 2019;20:760.
Lin Y, Wang K, Li X, Sun C, Yin R, Wang Y, et al. Evolution, functional differentiation, and co-expression of the RLK gene family revealed in Jilin ginseng, Panax ginseng C.A. Meyer. Mol Genet Genomics. 2018;293:845.
Yin R, Zhao M, Wang K, Lin Y, Wang Y, Sun C, et al. Functional differentiation and spatial-temporal co-expression networks of the NBS-encoding gene family in Jilin ginseng, Panax ginseng C.A. Meyer. PloS One. 2017;12:e0181596.
Wang L, Zhu J, Li X, Wang S, Wu J. Salt and drought stress and ABA responses related to bZIP genes from V. radiata and V. angularis. Gene. 2018;651:152–60.
Yoshida T, Fujita Y, Maruyama K, Mogami J, Todaka D, Shinozaki K, et al. Four ArabidopsisAREB/ABF transcription factors function predominantly in gene expression downstream of SnRK2 kinases in abscisic acid signalling in response to osmotic stress. Plant Cell Environ. 2015;38:35–49.
Bi C, Yu Y, Dong C, Yang Y, Zhai Y, Du F. The bZIP transcription factor TabZIP15 improves salt stress tolerance in wheat. Plant Biotechnol J. 2021;19:209–11.
Tang N, Zhang H, Li X, Xiao J, Xiong L. Constitutive activation of transcription factor OsbZIP46 improves drought tolerance in rice. Plant Physiol. 2012;158:1755–68.
Ji X, Liu G, Liu Y, Zheng L, Nie X, Wang Y. The bZIP protein from Tamarix hispida, ThbZIP1, is ACGT elements binding factor that enhances abiotic stress signaling in transgenic Arabidopsis. BMC Plant Biol. 2013;13:1–13.
Wang K, Jiang S, Sun C, Lin Y, Yin R, Wang Y, et al. The spatial and temporal transcriptomic landscapes of ginseng, Panax ginseng C.A. Meyer. Sci Rep. 2015;5:18283.
Wang Y, Li X, Lin Y, Wang Y, Wang K, Sun C, et al. Structural variation, functional differentiation, and activity correlation of the cytochrome P450 gene superfamily revealed in ginseng. Plant Genome. 2018;11:0106.
Kim NH, Jayakodi M, Lee SC, Choi BS, Jang W, Lee J, et al. Genome and evolution of the shade-requiring medicinal herb Panax ginseng. Plant Biotechnol J. 2018;16:1904–17.
Zhu L, Zhao M, Chen M, Li L, Jiang Y, Liu S, et al. The bHLH gene family and its response to saline stress in Jilin ginseng, Panax ginseng C.A. Meyer. Mol Genet Genomics. 2020;295:877–90.
Jin J, Zhang H, Kong L, Gao G, Luo J. PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 2014;42:D1182–7.
Bailey TL, Johnson J, Grant CE, Noble WS. The MEME Suite. Nucleic Acids Res. 2015;43:W39–49.
Brown P, Baxter L, Hickman R, Beynon J, Moore JD, Ott S. MEME-LaB: motif analysis in clusters. Bioinformatics. 2013;29:1696–7.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
Theocharidis A, van Dongen S, Enright AJ, Freeman TC. Network visualization and analysis of gene expression data using BioLayout Express3D. Nat Protoc. 2009;4(10):1535–50.
Nasrollahi V, Mirzaie-Asl A, Piri K, Nazeri S, Mehrabi R. The effect of drought stress on the expression of key genes involved in the biosynthesis of triterpenoid saponins in liquorice (Glycyrrhiza glabra). Phytochemistry. 2014;103:32–7.
Muscolo A, Sidari M, Anastasi U, Santonoceto C, Maggio A. Effect of PEG-induced drought stress on seed germination of four lentil genotypes. J Plant Interactions. 2014;9:354–63.
Ahmad MA, Javed R, Adeel M, Rizwan M, Yang Y. PEG 6000-stimulated drought stress improves the attributes of in vitro growth, steviol glycosides production, and antioxidant activities in Stevia rebaudiana Bertoni. Plants. 2020;9:1552.
Liu J, Wang Q, Sun M, Zhu L, Yang M, Zhao Y. Selection of reference genes for quantitative real-time PCR normalization in Panax ginseng at different stages of growth and in different organs. PloS One. 2014;9:e112177.
The authors thank the editor and reviewers for their valuable comments and suggestions for this manuscript.
This research was supported by an award from China 863 Project (2013AA102604-3), the Bureau of Science and Technology of Jilin Province (20200801063GH, 20190201264JC, 20180414077GH, 20180101027JC), and the Development and Reform Commission of Jilin Province (2018C047-3).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The PgbZIP genes identified in this study and their transcript sequences.
Table S2. The 62 PgbZIP genes that were spliced into transacripts having full-length ORFs and their putative protein sequences.
Table S3. Sequence alignment of the Jilin ginseng PgbZIP genes to the Korean Ginseng Genome Database.
Fig. S1. Comparison of the number of PgbZIP genes between the Korean Ginseng Genome Database and the Jilin Ginseng Transcriptome Database.
Table S4. The published bZIP genes of Arabidopsis, rice and tomato used as outgroups for phylogenetic analysis of the PgbZIP gene family.
Fig. S2. Phylogenetic tree of PgbZIP proteins constructed with the MP algorithm using the bZIP proteins of Arabidopsis, tomato and rice as the outgroups.
Table S5. Annotation and GO categorization of the 273 PgbZIP transcripts.
Table S6. Expressions of the 273 PgbZIP transcripts in the four-year-old roots of 42 genotypes, 14 tissues of a four-year-old plant and 4 different year-old roots.
Fig. S3. Percentage of the 273 PgbZIP transcripts expressed in different tissues (a), in the roots of differently aged plants (b), in the four-year-old roots of different genotypes (c).
Fig. S4. Relative water content (RWC) of ginseng seedlings subjected to PEG stress. The control plants were maintained under normal water irrigation conditions. The values are presented as the means of three replicates. “*”, P ≤ 0.05.
Table S7. Plant materials used for this study.
Table S8. Primers used for qRT-PCR analysis.
About this article
Cite this article
Li, H., Chen, J., Zhao, Q. et al. Basic leucine zipper (bZIP) transcription factor genes and their responses to drought stress in ginseng, Panax ginseng C.A. Meyer. BMC Genomics 22, 316 (2021). https://doi.org/10.1186/s12864-021-07624-z
- Panax ginseng
- bZIP transcription factor
- Functional differentiation
- Drought stress tolerance