- Research article
- Open Access
Transcriptome and metabolite analysis reveal the drought tolerance of foxtail millet significantly correlated with phenylpropanoids-related pathways during germination process under PEG stress
BMC Plant Biology volume 20, Article number: 274 (2020)
Foxtail millet [Setaria italica (L.) P. Beauv.] is an excellent crop known for its superior level of drought tolerance across the world. Especially, less water is needed during its germination period than the other cereal crops. However, the knowledge of the mechanisms underlying the abiotic stress effects on seed germination of foxtail millet is largely unknown.
The water uptake pattern of foxtail millet seeds was ploted during germination period, according to which the germination time course of millet was separated into three phases. We sequenced the transcriptome of foxtail millet seeds, which were treated by PEG during different germination phases after sowing. The transcriptional studies revealed that more DEGs were identified during the further increase in water uptake period (phase III) than during the rapid initial uptake period (phase I) and the plateau period (phase II) under PEG stress. The pathway analysis of DEGs showed that the highly enriched categories were related to phenylpropanoid biosynthesis, plant hormone signal transduction and phenylalanine metabolism during phase III. The 20 phenylpropanoids-related genes of germinating foxtail millet were found to be down-regulated during the further increase in water uptake period under PEG stress. Further expression analysis identified 4 genes of phenylalanine ammonia-lyase, 4-coumarate-CoA ligase 3, cinnamoyl-CoA reductase 1, cationic peroxidase SPC4 in phenylpropanoids-related pathway, which played important roles in foxtail millet in response to PEG stress during different germination periods. The studies of metabolites in phenylpropanoid biosynthesis pathway revealed that higher amount of cinnamic acid was accumulated in germinating seeds under PEG stress, while the contents of p-coumaric acid, caffeic acid, ferulic acid and sinapic acid were decreased. And the effects of five phenolic compounds on germination and growth of foxtail millet showed that 1 mM concentration of cinnamic acid inhibited shoot and root growth, especially root development. Ferulic acid, caffeic acid, sinapic acid and p-coumaric acid could increase the root length and root/sprout in lower concentration.
These findings suggest that key genes and metabolites of foxtail millet related with phenylpropanoids pathway may play prominent roles in the regulation of resistance to drought during germination. Foxtail millet can probably avoid drought by regulating the levels of endogenous allelochemicals.
Drought is one of the most frequent and severe abiotic stress factors, which adversely affects plant growth and crop productivity in many arid and semiarid regions [1,2,3]. Additionally, seasonal droughts often occur unevenly in the non-arid regions . In particular, droughts in spring severely impact the germination of seeds.
Seed germination commences with the uptake of water by the dry seed (imbibition), and is completed when usually the radicle extends to penetrate the structures that surround it . Uptake of water by a mature dry seed is triphasic, including a rapid initial uptake (phase I), a followed plateau phase (phase II), and a further increase in water uptake (phase III). Phase III occurs only after germination is completed, while the embryonic axes elongate and the radicles protrude, often called visible germination . During germination process, the influx of water into the cells of dry seeds results in rapidly resuming metabolic activity of the quiescent dry seed, and a series of complex structures, as well as physiological and molecular changes occur, such as temporary membrane structural perturbations, mobilization of the major storage reserves, recruitment of polysomes, the translation of preformed mRNAs, mRNA de novo synthesis [6, 7]. All the cellular and metabolic events occur in the nondormant seeds (before the completion of germination) and the imbibed dormant seeds, but the metabolic activities of the latter are subtly different from those of the former . The imbibed mature seed is sensitive to different environmental factors during germination process. However, the knowledge about the physiological and molecular mechanisms underlying the environmental effects on germination was largely lacking [8, 9].
Foxtail millet [Setaria italica (L.) P. Beauv.] is one of the oldest cereals and is thought to have played an important role in ancient civilization as a staple crop . Foxtail millet is known as a relatively drought-tolerant crop across the world, and grows in arid and semi-arid regions. Its morphological and anatomical characteristics give it strong drought resistance, such as thick cell wall, dense reticulate root system and small leaf area . Foxtail millet not only has strong drought tolerance, but also possesses small diploid genome sequence (about 515 Mb), strong inbreeding, short growth period and abundant germplasm resources. These features make millet an ideal model system for studying abiotic challenges [12, 13]. Moreover, it was reported that foxtail millet had significantly high water use efficiency compared with wheat, maize and sorghum . Especially, it needs less water in germination period than the other cereal crops, which can germinate when water absorption accounts for about 26% of seed weight . However, the metabolite changes and transcriptome reprogramming of foxtail millet in response to drought stress during germination process are rarely studied.
Phenylpropanoids, which are plant-specific natural products, play important functions during growth, development, and environmental interactions [15, 16]. Phenylpropanoids are precursors of a wide range of phenolic compounds, such as flavonoids, isoflavonoids, anthocyanins, plant hormones, phytoalexins, and lignins [17, 18]. Phenylpropanoids are synthesized from phenylalanine via the central phenylpropanoid pathway. The emergence of the phenylpropanoid pathway in plants is an important evolutionary adaptation that enables plant defense against abiotic and biotic stresses . Phenylpropanoids are in response to environmental cues and serve important functions in several different pathways including plant defense against pathogens and predators, protection from UV irradiation, signal transduction and communication with other organisms, and regulatory molecules [15, 17, 19,20,21,22]. And many phenylpropanoids and related metabolites have been reported to be allelochemicals, such as cinnamic acid and its hydroxylated derivatives, which can prohibit the germination of many plant seeds . However, the function of phenylpropanoids in seed germination of foxtail millet is still unclear.
RNA-Seq is an important technology that has been used to elucidate the complexity of regulation of gene expression during various stress conditions, as well as to obtain genome-wide estimations of relative gene expression. Rahman et al. sequenced the salinity responsive leaf transcriptome of the susceptible and tolerant finger millet. They found in the tolerant Trichy 1, the genes of several functional groups, such as transporters and transcription factors, were highly up-regulated, and genes involved in flavonoid biosynthesis were down-regulated specifically. Salinity stress inhibited photosynthesis related genes in the susceptible genotype . But the transcriptional profile analysis of the germinating seeds remains limited in foxtail millet under drought stress.
In this study, we aimed to understand the transcriptional and metabolic basis of drought responsiveness in foxtail millet during germination period. We plotted the figure associated with water uptake pattern of germinating seeds. The transcriptome analysis using RNA-Seq was performed for germination seeds in response to osmotic stress induced by polyethylene glycol (PEG). Many differentially expressed genes (DEGs) were identified separately in different phases of seed germination. Phenylpropanoid biosynthesis was found to be the highly enriched category. The related DEGs and metabolites of phenylpropanoid pathway were further analyzed in foxtail millet under PEG stress during different growth stages. These results provide an opportunity to elucidate the molecular mechanism underlying drought resistance of foxtail millet during germination process.
Phenotypic symptoms of foxtail millet during germination period under PEG stress
Total 8 foxtail millet cultivars were screened for their drought tolerances by germination percentage, relative sprout length and relative root length. Among them, cultivar Jingu 20 showed the highest germination percentage and relative root length, indicating its strongest dehydration tolerance (Shown in Additional file 1: Figure S1 and Additional file 2: Figure S2). Thus Jingu 20 was selected to understand the phenotypic and molecular basis of drought tolerance mechanisms.
The water uptake pattern and germination phenotype of Jingu 20 seeds were analyzed under 24 °C ~ 26 °C conditions in our culture room (Fig. 1). The pattern included three stages of a rapid initial uptake, a followed plateau phase and a further increase in water uptake (Fig. 1b), which were consistent with the previous results reported by Bewley . The rapid initial uptake period (phase I) was from 0 h to 6 h. The plateau phase (phase II) was from 6 h to 12 h. A further increase in water uptake (phase III) commenced at 12 h after sowing, at which the radicles of few seeds protruded. After sowing for 14 h, the radicles of all seeds penetrated the structures that surrounded them (Fig. 1a).
The phenotypic symptoms of Jingu 20 treated by PEG at 2 h (2H), 8 h (8H) and 14 h (14H) after sowing under water condition were shown in Fig. 1c. Using seeds grown at water condition (CK-1) as controls, PEG leaded to remarkable reduction in the growth of roots and shoots of 2H, 8H, 14H, CK-2 (seeds sowed under PEG stress) and CK-3 (seeds treated by PEG after absorbing water for 40 min). The length of roots and shoots were found to be affected more severely in CK-3 than in 2H, 8H, 14H. No significant differences were found in the length of the growth zone between of 2H, 8H, 14H and CK-2. And PEG stress slightly affected the growth of roots (Fig. 1d), but significantly reduced the length of shoot (Fig. 1e). These results indicated that the shoots suffered more inhibition than the roots in foxtail millet under PEG stress during germination.
Transcriptome analysis of foxtail millet in response to drought stress during seed germination
Drought responsive transcriptome profiling
To gain a global view on the drought induced changes at transcriptome level and drought responsive metabolic pathways in foxtail millet during germination, twelve cDNA libraries of Jingu 20 were sequenced from CK2H (collected during phase I, without PEG stress treatment), CK8H (collected during phase II, without PEG stress treatment)), CK14H (collected during phase III, without PEG stress treatment), P2H (collected during phase I, treated by PEG stress for 1 h and 3 h), P8H (collected during phase II, treated by PEG stress for 1 h and 3 h), P14H (collected during phase III, treated by PEG stress for 1 h and 3 h). After removing low quality reads and adapter sequences, a total of 38,845,760 and 39,679,812 (CK2H), 35,384,758 and 36,744,034 (P2H), 56,164,700 and 43,525,742 (CK8H), 52,411,398 and 47,479,886 (P8H), 46,838,670 and 49,151,796 (CK14H), 41,676,048 and 48,986,294 (P14H) clean reads were generated, which were listed in Additional file 8: Table S1. These libraries with Q30 > 85.02% were perfectly matched to the foxtail millet reference sequences from 78.50 to 82.84%.
Identification and functional annotation of differentially expressed genes
Differentially expressed genes (DEGs) during the progression of stresses show the stresses responsiveness and their putative roles in drought tolerance. DEGs were identified at a absolute threshold of fold change ≥2 and FDR ≤0.01 (FDR, false discovery rate). In phase I of foxtail millet germination, total 6 genes were detected in PEG stress samples (P2H) with CK2H (without PEG stress) as controls, of which 3 were up-regulated and 3 down-regulated (Additional file 3: Figure S3 and Table S3). At phase II, 42 genes were differentially regulated in the PEG stress samples (P8H) in comparison with CK8H (without PEG stress), including 24 up-regulated and 18 down-regulated genes. During phase III, a total of 302 up-regulated and 355 down-regulated genes were identified in PEG treated samples (P14H) when compared with the controls (CK14H, without PEG stress). The result showed that more DEGs were found during the further increase in water uptake period under drought stress.
COG refers to clusters of orthologous groups for eukaryotes. In this database, every protein is assumed to be evolved from a common ancestor protein. According to COG, 228 DEGs were classified into 21 different COG categories in the comparison of CK14H vs. P14H during phase III. Among these categories, the “general function prediction only”, “signal transduction mechanisms”and “transcription” were found to be the largest group in the compared groups (Additional file 4: Figure S4).
To further understand the biological functions of genes under drought stress, all DGEs were analyzed against the kyoto encyclopedia of genes and genomes (KEGG). In CK14H vs. P14H of phase III, 121 DEGs were mapped to the reference canonical pathways (Additional file: Table S3), which were classified into 70 functional categories (Additional file: Figure S5). Eight DEGs were found in CK8H vs. P8H, including 7 KEGG functional categories, and one DEG in CK2H vs. P2H (Additional file: Table S3). Using KEGG, the significantly enriched categories were identified, such as “phenylpropanoid biosynthesis, plant hormone signal transduction and phenylalanine metabolism” in CK14H vs. P14H (Fig. 2). Phenylpropanoid biosynthesis was dramatically enriched in CK8H vs. P8H (Additional file: Figure S6). The enrichment of other secondary metabolic pathways during different time point was not significant (Fig. 2 and Additional file: Figure S6 and Table S3). These findings suggested that many phenylpropanoids-related genes might play important roles in drought response during the further increase in water uptake period. Therefore, the following analysis mainly focused on phase III of foxtail millet germination under PEG stress.
Drought response transpoters in foxtail millet
Multiple transporters play important roles in plants in response to various environmental stimuli . In this study, 37 genes encoding transporters were identified in CK14H vs. P14H of Jingu 20 under PEG stress, and 1 gene in CK8H vs. P8H (Additional file: Table S4). Two outer envelope pore proteins and 2 S-adenosylmethionine carriers were up-regulated in CK14H vs. P14H. The expression levels of aquaporin PIP2–5, cation/calcium exchanger 1, glycerol-3-phosphate transporter 4, bidirectional sugar transporter SWEET1b and mannan endo-1,4-beta-mannosidase 1 were separately increased in the same sample. Amino-acid permease BAT1 homolog was found to be up-regulated in CK14H vs. P14H and CK8H vs. P8H. Cation/calcium exchanger NIP2–2 was down-regulated in CK14H vs. P14H under PEG stress.
Drought response signal transduction elements in foxtail millet
Protein kinases and phosphatases turn on or off stress responses by phosphorylation/ dephosphorylation to activate transcription factors and other genes related to stress tolerance [26, 27]. Total 38 protein kinases were found to be differentially regulated in CK14H vs. P14H (Additional file: Table S4). Among them, 6 genes were up-regulated, which encoded CBL-interacting protein kinase 16, SNF1-related protein kinase regulatory subunit gamma-like PV42a, SNF1-related protein kinase regulatory subunit gamma-1, two proline-rich receptor-like protein kinase PERK2, proline-rich receptor-like protein kinase PERK9. Out of 12 differentially expressed phosphatases, nine phosphatases 2C were found to be up-regulated in CK14H vs. P14H. The transcript of one phosphatase 2C was increased by PEG stress in CK8H vs. P8H of foxtail millet. Furthermore, 2 abscisic acid receptor PYL4 and 2 JAZ-like repressors of jasmonate signaling (Protein TIFY 10B and CASP-like protein 4A2) were found to be down-regulated under PEG stress in CK14H vs. P14H.
Drought response transcription factors in foxtail millet
Transcription factors regulate the expression of many downstream genes at the transcriptional level and control many biological processes such as cell division, growth, and response to environmental stress . In this study, 56 genes encoding transcription factors were identified in CK14H vs. P14H of Jingu 20 under PEG stress (Additional file: Table S4). These genes mainly included 14 ethylene-responsive transcription factor, 7 MYB, 6 heat stress transcription factor, 4 bZIP transcription factors, 3 zinc finger CCCH domain-containing protein, 3 WRKY, 3 transcription factor bHLH, 2 dehydration-responsive element-binding protein, 1 AP2/ERF and B3 domain-containing protein. The transcripts of 6 ethylene-responsive transcription factor genes were up-regulated, and 8 ones were down-regulated. Four bZIP transcription factors and three zinc finger CCCH domain-containing protein were found to be up-regulated. The levels of 4 heat stress transcription factor, 2 MYB, 1 WRKY and 1 dehydration-responsive element-binding protein also increased in CK14H vs. P14H. Three ethylene-responsive transcription factor and 1 homeobox-leucine zipper protein were identified in CK8H vs. P8H. Activation of ethylene-responsive transcription factor, heat stress transcription factor, MYB, WRKY, etc. showed that foxtail millet modulated stress tolerance during germination period under PEG conditions.
Drought responsiveness of genes related to phytohormone biosynthesis in foxtail millet
Phytohormones play important roles in response and adaptation to stress by reducing or mitigating the negative effects of stress . In CK14H vs. P14H, genes involved in auxin, cytokinine, ethylene, salicylic acid biosynthesis were found to be down-regulated (Additional file: Table S4). The expression levels of three 9-cis-epoxycarotenoid dioxygenase 1 and beta-carotene 3-hydroxylase related to abscisic acid biosynthesis were increased under PEG stress. Three genes encoding for gibberellin 2-beta-dioxygenase 1, ent-copalyl diphosphate synthase 1 and cytochrome P450 88A1 related to gibberellin biosynthesis were promoted by PEG stress. One gene of phenylalanine ammonia-lyase involved in salicylic acid biosynthesis was depressed by PEG treatment.
Regulation of genes involved in osmotic homeostasis under PEG stress
Plants often overcome the adverse effects of osmotic stress by accumulating metabolites or compatible solutes. Late embryogenesis abundant protein and dehydrin associated with tolerance against water stress . In CK14H vs. P14H, 9 genes encoding late embryogenesis abundant protein and 1 dehydrin were found to be up-regulated (Additional file: Table S4). Phosphoethanolamine N-methyltransferase involved in the biosynthesis of choline was up-regulated, which is the upstream gene of glycine betaine production. The transcript of phosphoethanolamine N-methyltransferase was increased in CK8H vs. P8H. Accumulation of proline is the adaptive response of plants against environmental stresses . One up-regulated gene of gamma-glutamyl phosphate reductase and one down-regulated gene of proline dehydrogenase 2 were found in CK14H vs. P14H, which might lead to the accumulation of proline. Sucrose and raffinose acts as compatible osmolytes in response to stress [30, 31]. In CK14H vs. P14H, one gene encoding alpha-galactosidase involved in production of galactose, glucose, sucrose and raffinose were found to be up-regulated in foxtail millet under PEG stress.
DEGs involved in phenylpropanoids-related pathway and their down-regulated expression pattern in response to drought
The metabolism of phenylalanine is an upstream pathway of phenylpropanoid. Phenylpropanoid biosynthesis provides the precursors for a wide range of phenolic compounds, such as ferulic acid, p-coumaric acid [15, 19]. During phase III, seventeen genes involved in phenylalanine metabolism were down-regulated in CK14H vs. P14H. The transcripts of 20 genes were depressed by PEG in the phenylpropanoid biosynthesis pathway (Fig. 3). After the duplication between two pathways was removed, 20 phenylpropanoids-related genes were found to be down-regulated in Jingu 20 under drought stress in phase III, which included 4-coumarate-CoA ligase 3, shikimate O-hydroxycinnamoyltransferase, cinnamoyl-CoA reductase 1, beta-glucosidase 6, phenylalanine ammonia-lyase and 15 peroxidases (Table 1). According to Fig. 3, the lower expression levels of these phenylpropanoids-related genes might affect the accumulation of phenolic compounds in the seeds under drought stress. And many phenylpropanoids and related metabolites have been reported to be allelochemicals and influence germination and growth of many plant, such as coumaric acid, ferulic acid . Therefore, the expression patterns of phenylpropanoids-related DEGs, the contents and roles of different phenolic compounds deserve further analysis in foxtail milllet.
Expression analysis of DEGs related with phenylpropanoids pathway using qRT-PCR
Since the most enriched pathway was “biosynthesis of phenylpropanoid”, eight DEGs encoding key enzymes in this pathway were selected for qRT-PCR validation (Table 1 and Fig. 4). In the sample of 14P1 (PEG treated for 1 h after germinating for 14 h under normal condition), the expression levels of Seita.1G240500, Seita.1G065800, Seita.1G361000, Seita.9G492600 and Seita.3G004800, decreased significantly compared with those in 14CK1 (germinating for 15 h under normal condition). While Seita.1G023100, Seita.4G047200 and Seita.5G145500 only showed similar expression trends. In 14P3 (PEG treated for 3 h after germinating for 14 h under normal condition), the genes, including Seita.1G065800, Seita.9G492600, Seita.3G004800, Seita.1G023100, Seita.4G047200 and Seita.5G145500, were dramatically down-regulated with 14CK3 (germinating for 17 h under normal condition) as control. While Seita.1G240500 and Seita.1G361000 just had decreased expression tendency. The comparison of qRT-PCR and RNA-seq assay data indicated a similar expression patterns in most of the selected DEGs except for the difference of fold change between them. The correlation between RNA-seq and qRT-PCR was analyzed in terms of fold changes. The Pearson coefficient was 0.54 (p value 0.03). These data supported the reliability of the sequencing results, and confirmed the down-regulated transcriptional expression of phenylpropanoids-related genes during seed germination under drought stress (Fig. 4).
To further verifying the down-regulated trends of DEGs in phenylpropanoids-related pathway during the germination period, eight key genes were analyzed by qRT-PCR at phase I (germinating for 2 h at normal condition), phase II (germinating for 8 h) and phase III (germinating for 14 h) under drought treatment (− 0.50 MPa PEG for 0 h, 1 h, 3 h, 12 h and 48 h) (Fig. 5). The results of qRT-PCR showed that the expression of phenylalanine ammonia-lyase (Seita.1G240500), 4-coumarate-CoA ligase 3 (Seita.1G065800), cinnamoyl-CoA reductase 1 (Seita.1G361000), beta-glucosidase 6 (Seita.9G492600), cationic peroxidase SPC4 (Seita.3G004800) and peroxidase 52 (Seita.1G023100) were up-regulated during phase I, while the expression of genes encoding shikimate O-hydroxycinnamoyl-transferase (Seita.4G047200) and peroxidase 4 (Seita.5G145500) only increased significantly at 48 h under PEG stress at phase I.
In phase II, the expression levels of four DEGs encoding 4-coumarate--CoA ligase 3 (Seita.1G065800), cinnamoyl-CoA reductase 1 (Seita.1G361000), cationic peroxidase SPC4 (Seita.3G004800) and peroxidase 52 (Seita.1G023100) decreased, while those of beta-glucosidase 6 (Seita.9G492600) were down-regulated only under PEG stress at 48 h. The expression of peroxidase 4 (Seita.5G145500) increased, whereas the expression levels of shikimate O-hydroxy- cinnamoyl-transferase (Seita.4G047200) were up-regulated only under PEG stress at 48 h. Additionally, the expression of gene encoding phenylalanine ammonia-lyase (Seita.1G240500) did not change obviously.
During phase III, the expression of five DEGs encoding phenylalanine ammonia-lyase (Seita.1G240500), 4-coumarate-CoA ligase 3 (Seita.1G065800), cinnamoyl-CoA reductase 1 (Seita.1G361000), beta-glucosidase 6 (Seita.9G492600) and cationic peroxidase SPC4 (Seita.3G004800) decreased under PEG stress, while the expression levels of shikimate O-hydroxycinnamoyl-transferase (Seita.4G047200) and peroxidase 4 (Seita.5G145500) were up-regulated, and the expression of peroxidase 52 (Seita.1G023100) remained unchanged. These results suggested that, the expression trends of phenylpropanoids related genes were different under drought stress during different germination period. And these genes probably play important roles in the regulation of germinating foxtail millet in response to drought stress.
Metabolite accumulation pattern in phenylpropanoid pathway of foxtail millet
Five phenylpropanoids-related metabolites were analyzed in germination seeds of millet treated by PEG at 14 h after growing under water condition. The results shown that the higher amount of cinnamic acid was accumulated in germinating seeds under PEG than that in the control. The levels of sinapic acid increased at 3 h, but decreased at 12 h, 48 h and 7 days under PEG conditions. However, the contents of p-coumaric acid and ferulic acid decreased during stress. The contents of caffeic acid declined significantly under PEG stress (Fig. 6).
The effects of phenolic compounds on foxtail millet germination
The five phenolic compounds were bioassayed for their effects on germination and growth of foxtail millet (Fig. 7). The analysis results indicated that the germination of foxtail millet was completely inhibited at 10 mM concentration of ferulic acid, p-coumaric acid, cinnamic acid and caffeic acid. The 10 mM concentration of sinapic acid dramatically reduced shoot and root growth. Ferulic acid, p-coumaric acid, cinnamic acid and caffeic acid remarkably inhibited shoots growth at 1 mM concentration. In particular, cinnamic acid more obviously restrained the growth of roots than the growth of shoots at 1 mM concentration (Fig. 7b). By contrast, cinnamic acid slightly promoted shoot growth at at 10− 6 M, 10− 5 M and 10− 4 M concentration. The 10− 6 M concentration of ferulic acid apparently promoted root growth and the root/sprout in this experiment. Sinapic acid, caffeic acid and p-coumaric acid separately promoted root growth at 0.01 mM concentration. Sinapic acid and caffeic acid helped noticeably the growth of root/sprout at concentrations of 0.01 M, 0.1 M and 1 mM. p-Coumaric acid promoted the root/sprout at 10− 6 M, 10− 5 M and 10− 4 M concentration.
Foxtail millet is known as a relatively drought-tolerant crop, which needs less water in germination period than the other cereal crops . Germination and growth of seeds depend on the germination conditions. Under 24 °C ~ 26 °C conditions in our culture room, the water uptake pattern of foxtail millet seeds was explored during germination period, which included phase I (from 0 to 6 h after sowing), phase II (from 6 to 12 h after sowing) and phase III (commencing at 12 h after sowing)(Fig. 1), which was consistent with the previous report .
According to three phases, twelve cDNA libraries have been sequenced in germinating seeds to gain a global view of the drought induced changes and germination regulation pattern at transcriptome level, and responsive metabolic pathways. Moreover, the stress time was very short, such as 1 h and 3 h. In response to PEG stress, total 302 up-regulated and 355 down-regulated genes were identified between CK14H and P14H in phase III, 24 up-regulated and 18 down-regulated genes between CK8H and P8H in phase II, and 3 up-regulated and 3 down-regulated genes between CK2H and P2H in phase I (Additional file 3: Figure S3). These results indicated that foxtail millet was more sensitive to drought during the further increase in water uptake period than in the rapid initial uptake and plateau period. Thus the subsequent analysis emphases were laid on the further increase in water uptake period.
Using KEGG enrichment, the pathways of phenylpropanoid biosynthesis, plant hormone signal transduction and phenylalanine metabolism were highly enriched in CK14H vs. P14H (Fig. 2). These findings suggested that many phenylpropanoids-related genes played important roles in drought response during phase III of foxtail millet. Similarly, the DEGs of germinating millet under drought stress for 10 h and 18 h were found to be associated with phenylpropanoids metabolism and plant hormone signal transduction . The changes of transcript levels of phenylpropanoids-related genes were found in Ricinus communis under low temperature . Soybean genes in the phenylpropanoid synthesis pathway were up-regulated in both alkaline and low Fe conditions .
Phenylalanine metabolism provides the precursors for the synthesis of phenylpropanoid . In our study, 17 genes involved in phenylalanine metabolism were down-regulated under PEG stress, and the transcripts of 20 genes were decreased in the phenylpropanoid biosynthesis pathway in CK14H vs. P14H (Fig. 3). A total of 20 phenylpropanoids-related genes were found to be down-regulated in foxtail millet under stress (Table 1). The identified phenylpropanoids-related genes were phenylalanine ammonia-lyase, 4-coumarate--CoA ligase 3, shikimate O-hydroxycinnamoyltransferase, cinnamoyl-CoA reductase 1, beta-glucosidase 6 and fifteen peroxidase. Phenylalanine ammonia lyase (PAL, EC126.96.36.199), the first enzyme in the phenylpropanoid biosynthesis pathway, catalyzes the conversion of L-phenylalanine to ammonia and trans-cinnamic acid . PAL is an inducible enzyme that responds to biotic and abiotic stresses, including pathogen infection, wounding, nutrient depletion, UV irradiation, and extreme temperatures [17, 20, 22, 34,35,36,37,38,39]. PAL plays an important role in plant defense, which is involved in the biosynthesis of salicylic acid (SA) [40,41,42]. 4-Coumarate--CoA ligase and cinnamoyl-CoA reductase (CCR) are important enzymes in catalyzing the conversion of coumaric acid, caffeic acid, ferulic acid, 5-hydroxyferulic acid and sinapic acid to lignin [43,44,45,46,47,48,49,50,51,52]. And many phenylpropanoids and their metabolites are known allelochemicals, for example cinnamic acid, p-coumaric acid, caffeic acid, ferulic acid and sinapic acid, which have been reported to prohibit the germination of many plant seeds . Therefore, these results indicated that the lower expression levels of genes in phenylpropanoid biosynthesis pathway might lead to the accumulation changes of allelochemicals in foxtail millet seeds, thus regulating the seed germination and seedling growth under drought stress.
To verify the expression pattern of these phenylpropanoids-related genes, eight key gene in phenylpropanoids-related pathway were analyzed by qRT-PCR at phase I (germinating for 2 h), phase II (germinating for 8 h), phase III (germinating for 14 h) under PEG treatment. These results showed that during different germination period, the expression trends of phenylpropanoids related genes were different under drought stress, respectively. For instance, when the seeds suffered from drought during phase I, the expression levels of these genes encoding phenylalanine ammonia-lyase (Seita.1G240500), 4-coumarate--CoA ligase 3 (Seita.1G065800), cinnamoyl-CoA reductase 1 (Seita.1G361000), beta-glucosidase 6 (Seita.9G492600) and cationic peroxidase SPC4 (Seita.3G004800) would be up-regulated, while the expression levels of these genes were down-regulated during phase III. In phase II, only 4-coumarate--CoA ligase 3 (Seita.1G065800), Cinnamoyl-CoA reductase 1 (Seita.1G361000) and Cationic peroxidase SPC4 (Seita.3G004800) were decreased. These results implied that the phenylpropanoids-related pathway played different roles in the regulation of foxtail millet in response to drought stress during different germination periods. In other words, the expression changes of phenylpropanoids related genes might increase or reduce the concentrations of phenylpropanoids related metabolites in the germination seeds, and thus result in promoting and inhibiting the germination and growth of foxtail millet under drought stress in different germination stages.
Water stress caused by osmotic stress or drought, can stimulate accumulation of allelochemicals . To confirm the accumulation changes of these phenylpropanoids-related metabolites, five metabolites related to phenylpropanoids were analyzed in germination seeds of foxtail millet. The analysis shown that the higher amount of cinnamic acid was accumulated in germinating seeds under PEG stress than that in the control. However, the contents of p-coumaric acid, caffeic acid, ferulic acid and sinapic acid were decreased. This is consistent with the conclusion that drought caused the increase of ferulic acid content in wheat .
The effects of the five phenolic compounds on germination and growth of foxtail millet were tested by external application. The highest concentrations of five phenolic compounds remarkably inhibited the germination or the growth of shoots and roots. In particular, cinnamic acid apparently restrained the growth of roots more than the growth of shoots at 1 mM concentration. The lower concentrations of ferulic acid, sinapic acid, caffeic acid and p-coumaric acid promoted root growth and root/sprout. They would jointly play roles at the same times and keep a coordination among themselves. Combined with the concentration changes of five metabolites in phenylpropanoids-related pathway under PEG stress, the content of cinnamic acid was dramatically increased and the amounts of ferulic acid, sinapic acid, caffeic acid and p-coumaric acid decreased in germinating seeds. These results implied that phenolic compounds would regulate the growth of foxtail millet under drought conditions. Accordingly, from phenotypic symptoms, drought stress slightly affected the growth of roots, while significantly reduced the length of shoots. Similar observations have been reported in Arabidopsis thaliana . Ferulic acid, p-coumaric acid and sinapic acid showed inhibitory effect on the germination and radicle length of A. thaliana. Ferulic acids inhibited total germination of A. thaliana above 500 and 750 μM. Wu et al. found that ferulic acids had no effect on annual bluegrass germination . Reigosa et al. reported that lower concentrations of phenolic compounds stimulated the growth of six weeds or had no actions . Ferulic and p-coumaric acids had inhibitory effects on root growth of A. thaliana without nutrients . And it was found that in tobacco p-coumaric inhibited root growth while ferulic acid had stimulatory effects . These results all indicate that various plant show different response to allelochemicals .
Due to lack of motility and immune system, plants have developed their defense strategies in the process of evolution, for example, the production of secondary metabolites as a tool in order to adapt to the changing environment and overcome stress constraints . These allelochemicals act as first line of defense against abiotic stresses [60, 61]. In particular, many desert plants have adapted to the harsh environment through secondary metabolites. Allelopathy of these chemicals can also affect neighboring plants . Allelopathy activity from plants grown in dry soils was greater than that provoked from plants grown in well-watered soils. Under water stress conditions, the donor plants contained a greater amount of allelopathic chemicals per dry weight than in absence of water stress, and the growth of target plants was reduced. This fact could cause autotoxicity problems in natural conditions [53, 63]. From the above results, foxtail millet has probably developed defense strategies in response to drought stress based on phenylpropanoids-related metabolites as desert plant, and avoided drought by their allelopathy and autotoxicity roles. Of course, these need further detailed experimental verification in the future.
In addition, in this study, many DEGs related to transport, signal transduction, phytohormone and osmotic homeostasis were identified in Jingu 20 during different germination periods under PEG stress. In CK14H vs. P14H, 37 transporters, 38 protein kinases, 2 PYL, 12 phosphatases, 2 JAZ, 56 transcription factors, 13 genes involved in phytohormone biosynthesis and 18 genes related to osmo-protention were differentially expressed under drought stress. These results showed that foxtail millet was in response to drought stress and modulating stress tolerance during germination period. With regard to the relation between phenylpropanoids-related pathway and the above pathways, it was only known that phenylalanine metabolism was upstream of salicylic acid and phenylpropanoid biosynthesis (Additional file: Figure S7). The knowledge of its mechanisms was largely unknown and needed to be further studied in foxtail millet.
The transcriptomic results demonstrate differentially expressed genes related to phenylpropanoids pathway involved in foxtail millet drought resisting. There were remarkable decreases in the expressions of phenylpropanoids-related genes during the further increase in water uptake period (phase III) of germination under PEG stress, accompanied by an increase of cinnamic acid, and reductions in p-coumaric acid, caffeic acid, ferulic acid and sinapic acid. The higher concentrations of external five phenolic acid inhibited seeds germination or the growth of shoots and roots, but the lower concentrations promoted the growth of seedlings or had no effects. The 1 mM concentration of cinnamic acid more obviously restrained the growth of roots than the growth of shoots. Ferulic acid, caffeic acid, sinapic acid and p-coumaric acid could increase the root length and root/sprout in lower concentration. Thus these findings indicate that foxtail millet could avoid drought by regulation of allelochemical concentrations. Allelochemicals would be accumulated and reduced by regulating the gene expressions of phenylpropanoids-related pathway during germination. These integrative results would provide valuable information for understanding the process of drought resistance during in foxtail millet germination, which would be benefit to foxtail millet cultivar breeding and innovation, promoting germination under drought stress, and seed priming.
Materials and methods
Plant materials and cultivar screening by PEG treatment
All foxtail millet cultivars were obtained from and preserved in Millet Research Institute of Shanxi Academy of Agricultural Sciences, which were planted every three years. We undertook the formal identification of the plant materials by phenotype and stress resistance. The voucher specimens of foxtail millet cultivars in this research were deposited in Physiological and Biochemical Laboratory. The deposition numbers of Yugu1, An04–4783, Changnong35, Jingu20, Jingu34, Jingu9, Tieganhan, Lugu6 were si-2015-b13, si-2015-358, si-2015-g54, si-2015-g1, si-2015-g3, si-2015-b15, si-2015-b19, si-2015-b21, respectively. We evaluated the cultivar variations of the 8 foxtail millet genotypes in the context of early drought tolerance using polyethylene glycol (PEG) induced osmotic stress at the germination and early seedling growth stage. The stress was imposed by exposing the germinating grain to polyethylene glycol. The osmotic potential of polyethylene glycol 6000 was − 0.50 MPa . The germination percentage, the relative sprout length and the relative root length were observed against the controls. From these screening, the cultivar of highest drought resistance was selected to understand the phenotypic and molecular basis of tolerance mechanisms.
Water uptake pattern of foxtail mellit seeds during germination stage
Fifty seeds of the selected cultivar were allowed for germination under water in pertri-dishes, each was replicated three times. Seeds germinated during June–August when the temperature remained between 24 °C ~ 26 °C in culture room. Their fresh weights were measured every one hour by drying them with absorbent paper during germination period in foxtail millet. Then the growth rate of seed fresh weight was calculated against dry seed weight. The map was draw by the growth rate of seed fresh weight.
Effect of PEG stress during germination period and sample preparation
Seeds of the selected cultivar were surface-sterilized in 3% sodium hypochloride for 20 min and rinsed 5 times (4 min/each) with distilled water. These seeds were allowed for germination in 9 cm petri-dishes at room temperature ranged from 24 °C to 26 °C and relative humidty ranged from 40 to 50% in the culture room. Each was replicated three times. After germinating for 2 h, 8 h and 14 h (including the time of sterilizing and growing) under normal conditions, the seeds were transplanted in plates containing polyethylene glycol 6000 solution. The osmotic potential of polyethylene glycol 6000 was − 0.50 MPa. The control plants were cultured and transferred in the same ways as the dehydration treatments, but without the addition of PEG. The shoot length and root length of foxtail millet stressed by PEG were measured on the 7th day of germination. Under the same conditions, the germinating seeds used for sequencing were collected, which were CK2H, CK8H, CK14H (germinating respectively for 2 h, 8 h, 14 h under normal condition respectively, then growing for 1 h and 3 h without PEG stress treatment), and P2H, P8H, P14H (germinating separately for 2 h, 8 h, 14 h under normal condition respectively, then growing for 1 h and 3 h under PEG treatment). Three biological replicates were set at the same time. Each biological duplication was carried out every 48 h. Three biological repeats were mixed into one sample. The experiments were repeated twice for collecting the sequencing samples. The first repeated samples were − 1 of sequencing samples, and the second repeated ones were − 2 of sequencing samples. The seeds treated by PEG were sampled at the times (1 h, 3 h, 12 h and 48 h) respectively. At the same time, the seeds growing at normal condition were collected during germination phases for expression analysis by qRT-PCR. The collected samples were frozen in liquid nitrogen and stored at − 80 °C for RNA Extraction.
RNA extraction, library preparation, and transcriptome sequencing
Frozen seed samples were ground in liquid nitrogen and total RNA was extracted using CTAB methods. RNA degradation and contamination were monitored on 1% agarose gels. RNA purity was checked using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA). The transcriptome cDNA libraries were prepared using NEBNext®Ultra™ RNA Library Prep Kit (NEB, USA) and purified with AMPure XP system (Beckman Coulter, Beverly, USA) according to manufacturer’s protocols. Sequencing was performed on an Illumina Hiseq 2500 following manufacturer’s recommendations (Beijing, China).
After RNA sequencing, the raw data of fasta format were firstly processed through in-house perl scripts. Clean reads were obtained by removing reads containing adapter, reads containing ploy-N and low quality reads from raw data. Q30 and GC-content of the clean data were calculated. The clean reads were mapped and annotated based on the reference genome. Tophat2 was used as the tool soft for maping . Read count for each gene was obtained from the mapping results. Gene expression levels were estimated by RSEM for each sample . DESeq R package was used to analyze differential expression of two groups . P value corrections were performed using the Benjamini and Hochberg’s approach. The genes with corrected P-value less than 0.05 were considered differentially expressed.
Gene functional annotation
Gene function was annotated based on the following databases: Nr (NCBI non-redundant protein sequences), Nt (NCBI non-redundant nucleotide sequences), Pfam (Protein family), COG (Clusters of orthologous groups of proteins), Swiss-Prot (A manually annotated and reviewed protein sequence database), KO (KEGG ortholog database) and GO (Gene ontology). The enrichment of differential expression genes in KEGG pathways was analyzed using KOBAS software .
Validation of differential expression using qRT-PCR
Expression of eight genes involved in phenylpropanoids-related pathway under PEG stress were examined by qRT-PCR analysis using SYBR Premix Ex Taq II (TaKaRa, Dalian, China) on Thermal Cycler Dice Real Time System (TaKaRa Code.TP800, Japanese). The PCR reactions were performed according to the manufacturer’s protocol. Each 25 μL amplification reaction contained 12.5 μL of SYBR Premix Ex Taq II (TaKaRa, Dalian, China), 1 μL of each primer (10 μM), 8.5 μL of sterile distilled water, and 2 μL of cDNA template. The PCR reaction contained an initial denaturation (95 °C for 30 s) followed by 40 cycles of 95 °C denaturation for 5 s and 60 °C annealing for 30 s. Gene specific primers were designed according to non-conserved region sequences of each gene using Primer Premier 5.0 software. Gene-specific primers were listed in Table S2. β-Actin gene of foxtail millet (Seita.7G294000) was used as internal controls for normalization as described by Xu . The relative expression levels of the genes under PEG stress during different stages were calculated via the double standard curve method .
Determination of metabolite content in phenylpropanoid biosynthesis of foxtail millet
The contents of cinnamic acid, p-coumaric acid, caffeic acid, ferulic acid and sinapic acid in germinating seeds were analyzed by HPLC . The seeds of foxtail millet germinated for 14 h under water conditions, then grew for 3 h, 12 h, 48 h and 7 days under PEG treatment (− 0.50 MPa).
The effects of phenolic compounds on the germination of foxtail millet
Cinnamic acid, p-coumaric acid, caffeic acid, ferulic acid and sinapic acid were tested for their effects on the germination and growth of foxtail millet. Five phenolic compounds were prepared in concentrations of 10 mM, 1 mM, 0.1 mM, 0.01 mM and 0.001 mM . Solutions were bioassayed on seeds of foxtail millet. Twenty five seeds were placed on Whatman 3 MM paper in 9 cm diameter Petri dishes, to which 7 ml of solutions were add at the beginning and two additional ml of each solution were added on the third day. Three replicates of each treatment were incubated in the germination chamber under the same conditions as mentioned above.
Availability of data and materials
The RNA-Seq data has been deposited in the Sequence Read Archive (SRA) at the National Center for Biotechnology Information (NCBI). The accession number is PRJNA573803, which includes 12 accession items (SRX6967585, SRX6967584, SRX6967583, SRX6967582, SRX6967581, SRX6967580, SRX6967579, SRX6967578, SRX6967577, SRX6967576, SRX6967575, SRX69675794).
Differentially expressed genes,
False discovery rate
Clusters of orthologous groups
Kyoto encyclopedia of genes and genomes
Quantitative real time polymerase chain reaction
Phenylalanine ammonia lyase
Cetyltrimethylammonium ammonium bromide
Fragments per kilobase of exon model per million mapped reads
Polymerase chain reaction
High performance liquid chromatography
Qiao G, Wen XP, Yu LF, Ji XB. Identification of differentially expressed genes preferably related to drought response in pigeon pea (Cajanus cajan) inoculated by arbuscular mycorrhizae fungi (AMF). Acta Physiol Plant. 2012;34:711–1721. https://doi.org/10.1007/s11738-012-0966-2..
Zhu XC, Song FB, Liu SQ, Liu TD. Arbuscular mycorrhizae improves photosynthesis and water status of Zea mays L. under drought stress. Plant Soil Environ. 2012;58:186–91. https://doi.org/10.1007/s11032-011-9671-x.
Omirou M, Ioannis MI, Constantinos E. Mycorrhizal inoculation affects arbuscular mycorrhizal diversity in watermelon roots, but leads to improved colonization and plant response under water stress only. Appl Soil Ecol. 2013;63:112–9. https://doi.org/10.1016/j.apsoil.2012.09.013.
Bewley JD, Black M. Seeds: physiology of development and germination. New York: Plenum Press; 1994. https://doi.org/10.1007/978-1-4899-1002-8.
Bewley JD. Seed Germination and Dormancy. Plant Cell. 1997;9:1055–66. https://doi.org/10.1105/tpc.9.7.1055.
Catusse J, Job C, Job D. Transcriptome- and proteome-wide analyses of seed germination. CR Biol. 2008;331:815–22. https://doi.org/10.1016/j.crvi.2008.07.023.
Rajjou L, Duval M, Gallardo K, Catusse J, Bally J, Job C, Job D. Seed germination and vigor. Annu Rev Plant Biol. 2012;63:507–33. https://doi.org/10.1146/annurev-arplant-042811-105550.
Wang X, Wang L, Yan X, Wang L, Tan M, Geng X, Wei W. Transcriptome analysis of the germinated seeds identifies low-temperature responsive genes involved in germination process in Ricinus communis. Acta Physiol Plant. 2016;38(1):6. https://doi.org/10.1007/s11738-015-1994-5.
Tan L, Chen S, Wang T, Dai S. Proteomic insights into seed germination in response to environmental factors. Proteomics. 2013;13:1850–70. https://doi.org/10.1002/pmic.201200394.
Fukunaga K, Wang ZM, Kato K, Kawase M. Geographical variation of nuclear genome RFLPs and genetic differentiation in foxtail millet, Setaria italica (L.) P. Beauv. Genet Res Crop Evol. 2002;49:95–101. https://doi.org/10.1023/a:1013852007770.
Diao X. Foxtail millet production and future development direction in China. In: Chai Y, Wan SH, editors. Reports on minor grain development in China. Beijing: Chinese Agri Sci Tech Press; 2007. p. 32–43.
Lata C, Gupta S, Prasad M. Foxtail millet: a model crop for genetic and genomic studies in bioenergy grasses. Crit Rev Biotechnol. 2013;33:328–43. https://doi.org/10.3109/07388551.2012.716809.
Muthamilarasan M, Prasad M. Advances in Setaria genomics for genetic improvement of cereals and bioenergy grasses. Theor Appl Genet. 2015;128:1–14. https://doi.org/10.1007/s00122-014-2399-3.
Li R, Zhao Z, Mao X. Principle of foxtail millet genetics and breeding. Beijing: China Meteorological Press; 1997. p. 26.
Vogt T. Phenylpropanoid biosynthesis. Mol Plant. 2010;3:2–20. https://doi.org/10.1093/mp/ssp106.
Liu XY, Wang PP, Wu YF, Cheng AX, Lou HX. Cloning and functional characterization of two 4-Coumarate: CoA ligase genes from Selaginella moellendorffii. Molecules. 2018;23:595. https://doi.org/10.3390/molecules23030595.
Dixon RA, Paiva N. Stress-induced phenylpropanoid metabolism. Plant Cell. 1995;7:1085–97. https://doi.org/10.1105/tpc.7.7.1085.
La CS, Gouzerh G, Dhondt S, Hoffmann L, Fritig B, Legrand M, Heitz T. Metabolic reprogramming in plant innate immunity: the contributions of phenylpropanoid and oxylipin pathways. Immunol Rev. 2004;198:267–84. https://doi.org/10.1111/j.0105-2896.2004.0129.x.
Ferrer JL, Austin MB, Stewart C Jr, Noel JP. Structure and function of enzymes involved in the biosynthesis of phenylpropanoids. Plant Physiol Biochem. 2008;46:356–70. https://doi.org/10.1016/j.plaphy.2007.12.009.
Payyavula RS, Navarre DA, Kuhl JC, Pantoja A, Pillai SS. Differential effects of environment on potato phenylpropanoid and carotenoid expression. BMC Plant Biol. 2012;12:39. https://doi.org/10.1186/1471-2229-12-39.
Hahlbrock K, Scheel D. Physiology and molecular biology of phenylpropanoid metabolism. Annu Rev Plant Physiol Plant Mol Biol. 1989;40:347–69. https://doi.org/10.1146/annurev.pp.40.060189.002023.
Kim DS, Hwang BK. An important role of the pepper phenylalanine ammonia-lyase gene (PAL1) in salicylic acid-dependent signalling of the defence response to microbial pathogens. J Exp Bot. 2014;65:2295–306. https://doi.org/10.1093/jxb/eru109.
Reigosa MJ, Malvido-Pazos E. Phytotoxic effects of 21 plant secondary metabolites on Arabidopsis thaliana germination and root growth. J Chem Ecol. 2007;33:1456–66. https://doi.org/10.1007/s10886-007-9318-x.
Rahman H, Jagadeeshselvam N, Valarmathi R, Sachin B, Sasikala R, Senthil N, Sudhakar D, Robin S, Muthurajan R. Transcriptome analysis of salinity responsiveness in contrasting genotypes of finger millet (Eleusine coracana L.) through RNA-sequencing. Plant Mol Biol. 2014;85:485–503. https://doi.org/10.1007/s11103-014-0199-4.
Sun X, Xu L, Wang Y, Luo X, Zhu X, Kinuthia KB, Nie S, Feng H, Li C, Liu L. Transcriptome-based gene expression profiling identifies differentially expressed genes critical for salt stress response in radish (Raphanus sativus L.). Plant Cell Rep. 2015;10. https://doi.org/10.1007/s00299-015-1887-5.
Ozturk ZN, Talamé V, Deyholos M, Michalowski CB, Galbraith DW, Gozukirmizi N, Tuberosa R, Bohnert HJ. Monitoring large-scale changes in transcript abundance in drought-and salt-stressed barley. Plant Mol Biol. 2002;48:551–73. https://doi.org/10.1023/a:1014875215580.
Zhang X, Zhen J, Li Z, Kang D, Yang Y, Kong J, Hua J. Expression profile of early responsive genes under salt stress in upland cotton (Gossypium hirsutum L.). Plant Mol Biol Rep. 2011;29:626–37. https://doi.org/10.1007/s11105-010-0269-y.
Bartels D, Sunkar R. Drought and salt tolerance in plants. Crit Rev Plant Sci. 2005;24:23–58. https://doi.org/10.1080/07352680590910410.
Delauney AJ, Verma DPS. Proline biosynthesis and osmoregulation in plants. Plant J. 1993;4:215–23. https://doi.org/10.1046/j.1365-313X.1993.04020215.x.
Ashraf M, Harris P. Potential biochemical indicators of salinity tolerance in plants. Plant Sci. 2004;166:3–16. https://doi.org/10.1016/j.plantsci.2003.10.024.
Nishizawa A, Yabuta Y, Shigeoka S. Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 2008;147:1251–63. https://doi.org/10.1104/pp.108.122465.
Xu BX, Yin MQ, Wen YY, Pei SS, Ke ZJ, Zhang B, Yuan XY. Gene expression profiling of foxtail millet (Setaria italica L.) under drought stress during germination. Sci Agric Sin. 2018;51:1431–47.
Waters BM, Amundsen K, Graef G. Gene expression profiling of iron deficiency chlorosis sensitive and tolerant soybean indicates key roles for phenylpropanoids under alkalinity stress. Front Plant Sci. 2018;9:10. https://doi.org/10.3389/fpls.2018.00010.
Edwards K, Cramer CL, Bolwell GP, Dixon RA, Schuch W, Lamb CJ. Rapid transient induction of phenylalanine ammonia-lyase mRNA in elicitor-treated bean cells. Proc Natl Acad Sci USA. 1985;82:6731–5. https://doi.org/10.1073/pnas.82.20.6731.
Liang XW, Dron M, Cramer CL, Dixon RA, Lamb CJ. Differential regulation of phenylalanine ammonia-lyase genes during plant development and by environmental cues. J Biol Chem. 1989;264:14486–92. https://doi.org/10.1016/0092-8674(89)90112-8.
Liang XW, Dron M, Schmid J, Dixon RA, Lamb CJ. Developmental and environmental regulation of a phenyl- alanine ammonia-lyase-β-glucuronidase gene fusion in transgenic tobacco plants. Proc Natl Acad Sci USA. 1989;86:9284–8. https://doi.org/10.1073/pnas.86.23.9284.
Huang J, Gu M, Lai Z, Fan B, Shi K, Zhou YH, Yu JQ, Chen Z. Functional analysis of the Arabidopsis PAL gene family in plant growth, development, and response to environmental stress. Plant Physiol. 2010;153:1526–38. https://doi.org/10.1104/pp.110.157370.
Jin Q, Yao Y, Cai Y, Lin Y. Molecular cloning and sequence analysis of a phenylalanine ammonia-lyase gene from Dendrobium. PLoS One. 2013;8:e62352. https://doi.org/10.1371/journal.pone.0062352.
MacDonald MJ, D’Cunha GB. A modern view of phenylalanine ammonia lyase. Biochem Cell Biol. 2007;85:273–82. https://doi.org/10.1139/o07-018.
MauchMani B, Slusarenko AJ. Production of salicylic acid precursors is a major function of phenylalanine mmonia-lyase in the resistance of Arabidopsis to Peronospora parasitica. Plant Cell. 1996;8:203–12. https://doi.org/10.2307/3870265.
Nugroho LH, Verberne MC, Verpoorte R. Activities of enzymes involved in the phenylpropanoid pathway in constitutively salicylic acid-producing tobacco plants. Plant Physiol Biochem. 2002;40:775–60. https://doi.org/10.1016/s0981-9428(02)01437-7.
Chaman ME, Copaja SV, Argandona VH. Relationships between salicylic acid content, phenylalanine ammonia-lyase (PAL) activity, and resistance of barley to aphid infestation. J Agric Food Chem. 2003;51:2227–31. https://doi.org/10.1021/jf020953b.
Dixon RA, Achnine L, Kota P, Liu CJ, Reddy MS, Wang L. The phenylpropanoid pathway and plant defence-a genomics perspective. Mol Plant Pathol. 2002;3:371–90. https://doi.org/10.1046/j.1364-3703.2002.00131.x.
Boerjan W, Ralph J, Baucher M. Lignin biosynthesis. Annu Rev Plant Biol. 2003;54:519–46. https://doi.org/10.1146/annurev.arplant.54.031902.134938.
Piquemal J, Lapierre C, Myton K, O’Connell A, Schuch W, Grima-Pettenati J, Boudet A. Downregulation of cinnamoyl CoA reductase induces significant changes of lignin profiles in transgenic tobacco plants. Plant J. 1998;13:71–83. https://doi.org/10.1046/j.1365-313X.1998.00014.x.
Jones L, Ennos AR, Turner SR. Cloning and characterization of irregular xylem4 (irx4): a severely lignin-deficient mutant of Arabidopsis. Plant J. 2001;26:205–16. https://doi.org/10.1046/j.1365-313x.2001.01021.x.
Goujon T, Ferret V, Mila I, Pollet B, Ruel K, Burlat V, Joseleau JP, Barriere Y, Lapierre C, Jouanin L. Down- regulation of the AtCCR1 gene in Arabidopsis thaliana: effects on phenotype, lignins and cell wall degradability. Planta. 2003;217:218–28. https://doi.org/10.2307/23387842.
Kawasaki T, Koita H, Nakatsubo T, Hasegawa K, Wakabayashi K, Takahashi H, Umemura K, Umezawa T, Shimamoto K. Cinnamoyl-CoA reductase, a key enzyme in lignin biosynthesis, is an effector of small GTPase Rac in defense signaling in rice. Proc Natl Acad Sci USA. 2006;103:230–5 www.pnas.org/cgi/doi/10.1073/pnas.0509875103.
Leple JCI. Downregulation of cinnamoyl-coenzyme a reductase in poplar: multiple-level phenotyping reveals effects on cell wall polymer metabolism and structure. Plant Cell. 2007;19:3669–91. https://doi.org/10.1105/tpc.107.054148.
Wadenback J, von Arnold S, Egertsdotter U, Walter MH, Grima-Pettenati J, Goffner D, Gellerstedt G, Gullion T, Clapham D. Lignin biosynthesis in transgenic Norway spruce plants harboring an antisense construct for cinnamoyl CoA reductase (CCR). Transgenic Res. 2008;17:379–92. https://doi.org/10.1007/s11248-007-9113-z.
Zhou R, Jackson L, Shadle G, Nakashima J, Temple S, Chen F, Dixon RA. Distinct cinnamoyl CoA reductases involved in parallel routes to lignin in Medicago truncatula. Proc Natl Acad Sci USA. 2010;107:17803–8. https://doi.org/10.2307/20780539.
Tamasloukht B, Wong Quai Lam MS, Martinez Y, Tozo K, Barbier O, Jourda C, Jauneau A, Borderies G, Balzergue S, Renou JP, Huguet S, Martinant JP, Tatout C, Lapierre C, Barrière Y, Goffner D, Pichon M. Characterization of a cinnamoyl-CoA reductase 1(CCR1) mutant in maize:effects on lignification,fibre development,and global gene expression. J Exp Bot. 2011;62:3837–48. https://doi.org/10.1093/jxb/err077.
Reigosa M J, Pedrol N, González, Luís. Allelopathy || Allelopathy and abiotic stress. 2006. https://doi.org/10.1007/1-4020-4280-9_9.
Hura T, Grzesiak S, Hura K, Thimet E, Tokarz K, Wedzony M. Physiological and biochemical tools useful in drought-tolerance detection in genotypes of winter triticale: accumulation of ferulic acid correlates with drought tolerance. Ann Bot. 2007;100:767–75. https://doi.org/10.1093/aob/mcm162.
Wu L, Guo X, Harivandi MA. Allelopathic effects of phenolic acids detected in buffalograss (Buchloe dactyloides) clippings on growth of annual bluegrass (Poa annua) and buffalograss seedlings. Environ Exp Bot. 1998;39:159–67. https://doi.org/10.1016/s0098-8472(97)00040-3.
Reigosa MJ, Souto XC, González L. Effect of phenolic compounds on the germination of six weeds species. Plant Growth Reg. 1999;28:83–9. https://doi.org/10.1023/a:1006269716762.
Reddy S, Rao SSR. Effect of hydroxy cinnamic acids on growth, nucleic acid, carbohydrate and nicotine content of tobacco (Nicotiana tabacum L.). Indian J Plant Physiol. 1999;4:167–70.
Baleroni CSS, Ferrarese MLL, Braccini AL, Scapim CA, Ferrarese-Filho O. Efects of ferulic and p-coumaric acids on canola (Brassica napus L. cv. Hyola 401) seed germination. Seed Sci Tech. 2000;28:333–40.
Maqbool N, Wahid A, Farooq M, Cheema ZA, Siddique KHM. Allelopathy and abiotic stress interaction in crop plants. Berlin Heidelberg: Springer; 2013. https://doi.org/10.1007/978-3-642-30595-5_19.
Bacilio-Jiménez M, Aguilar-Flores S, Ventura-Zapata E, Pérez-Campos E, Bouquelet S, Zenteno E. Chemmical characterization of root exudates from rice (Oryza sativa) and their effect on the chemotatic response of endophytic bacteria. Plant Soil. 2003;249:271–7. https://doi.org/10.1023/a:1022888900465.
Asao T, Kitawaza H, Tomita K, Suyama K, Yamamoto H, Hosoki T, Pramanik MHR. Mitigation of cucumber autotoxicity in hydroponic cultureusing microbial strain. Sci Hort. 2004;99:207–14. https://doi.org/10.1016/s0304-4238(03)00098-0.
Friedman J. Allelopathy in desert ecosystems; 1987. https://doi.org/10.1021/bk-1987-0330.ch006.
Tongma S, Kobayashi K, Usui K. Allelopathic activity of Mexican sunflower [Tithonia diversifolia (Hemsl.) a. gray] in soil under natural field conditions and different moisture conditions. Weed Biol Manage. 2001;1:115–9. https://doi.org/10.1046/j.1445-6664.2001.00020.x.
Michel BE, Kaufmann MR. The osmotic potential of polyethylene glycol 6000. Plant Physiol. 1973;51:914–6. https://doi.org/10.1104/pp.51.5.914.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36. https://doi.org/10.1186/gb-2013-14-4-r36.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323. https://doi.org/10.1186/1471-2105-12-323.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106. https://doi.org/10.1186/gb-2010-11-10-r106.
Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21:3787–93. https://doi.org/10.2307/1592215.
Xu Y, Hui L, Li X, Jing L, Wang Z, Yang Q, Chang Y. Systematic selection and validation of appropriate reference genes for gene expression studies by quantitative real-time PCR in pear. Acta Physiol Plant. 2015;37:1–16. https://doi.org/10.1007/s11738-015-1784-0.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001;29:900. https://doi.org/10.1093/nar/29.9.e45.
Molnár-Perl I, Füzfai Z. Chromatographic, capillary electrophoretic and capillary electrochromatographic techniques in the analysis of flavonoids. J Chromatogr A. 2005;1073:201–27. https://doi.org/10.1016/j.chroma.2004.10.068.
We would like to thank Xianmin Diao for good advice in the experiments.
This work was supported by China Agriculture Research System [CARS-06-13.5-A23], the Project Plan of Shanxi Academy of Agricultural Sciences [YGG17021 and YCX2018206]. The funder had no role in this study design, data collection and analysis, decision to publish, or preparation of the manuscript.
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.
Additional file 1: Figure S1. Differential responses of foxtail millet cultivars to PEG stress during germination stage. Phenotypes of eight foxtail millet cultivars were observed on the seventh day under PEG conditions during germination period. The cultivars in first line from left to right were Yugu1, An04–4783, Changnong35 and Jingu20, respectively. The genotypes in second line from left to right were Jingu34, Jingu9, Tieganhan and Lugu6, respectively.
Additional file 2: Figure S2. Germination rates of foxtail millet cultivars under PEG stress during germination period. Germination rates of Yugu1, An04–4783, Changnong35, Jingu20, Jingu34, Jingu9, Tieganhan and Lugu6 were counted on the 7th day under PEG conditions. The mean values and SD were calculated using one-way ANOVA followed by Tukey HSD multiple comparison.
Additional file 3: Figure S3. Volcano plot of genes with differential expression between control and treated sample libraries. The X-axis indicated the log2(FC) of DEGs (FC, fold change). The Y-axis represented the -log10(FDR) of differential expression genes (FDR, False Discovery Rate). Red dots indicated up-regulated genes. Green dots represented down-regulated genes, and black dots indicated non-different expression genes. DEGs between control and PEG stress were identified in different germination periods according to a threshold of fold change ≥2 and FDR ≤0.01 (FDR, false discovery rate). A, DEGs of CK2H vs. P2H. (n = 6; 3 up and 3 down) B, DEGs of CK8H vs. P8H. (n = 42; 24 up and 18 down) C, DEGs of CK14H vs. P14H (n = 657; 302 up and 355 down).
Additional file 4: Figure S4. COG classifications of DEGs of foxtail millet in CK14H vs. P14H. A total of 228 DEGs in CK14H vs. P14H were assigned to 21 COG categories. The capital letters on the x-axis indicated the COG categories as listed on the right of the histogram. The y-axis represented the percentages of the corresponding category among all categories. The number of DEGs belong to the category, and the percentages of the corresponding category among all categories were separately indicated in square brackets.
Additional file 5: Figure S5. KEGG classifications of DEGs of foxtail millet in CK14H vs. P14H. The x-axis indicated the percentages of DEGs among the total annotated genes. The left y-axis represented the pathways. The number of DEGs belonged to the annotated pathway were represented after the histograms.
Additional file 6: Figure S6. The pathways with the most significant Q value in CK8H vs. P8H. The x-axis indicated rich factor of DEGs belong to the corresponding pathway. The left y-axis represented the pathways. The sizes of bubble represented the number of DEGs in the corresponding pathway, and the colors of the bubble represented the enrichment Q value of the corresponding pathway.
Additional file 8: Table S1 - Summary of the sequencing and the reads mapping from the control groups (CK) and the PEG stress groups (P). The sequencing samples were CK2H (germinating for 2 h under normal conditions, then growing for 1 h and 3 h without PEG stress treatment), P2H(germinating for 2 h, then growing for 1 h and 3 h under PEG treatment), CK8H (germinating for 8 h under normal conditions, then growing for 1 h and 3 h without PEG stress treatment), P8H(germinating for 8 h under normal conditions, then growing for 1 h and 3 h under PEG treatment), CK14H(germinating for 14 h under normal conditions, then growing for 1 h and 3 h without PEG stress treatment), P14H (germinating for 14 h under normal conditions, then growing for 1 h and 3 h under PEG treatment). -1 represented the first repeat. -2 represented the second repeat.
Additional file 11: Table S4. Significant DEGs involved in PEG stress during different germination periods. DEGs related with transporters, signal transduction components, transcription factors, phytohormones biosynthesis and osmotic homeostasis were selected according to gene functions and metabolic pathways.
About this article
Cite this article
Yu, A., Zhao, J., Wang, Z. et al. Transcriptome and metabolite analysis reveal the drought tolerance of foxtail millet significantly correlated with phenylpropanoids-related pathways during germination process under PEG stress. BMC Plant Biol 20, 274 (2020). https://doi.org/10.1186/s12870-020-02483-4
- Foxtail millet
- Transcriptome analysis
- Phenylpropanoids-related pathways and metabolites
- PEG stress