Transcriptome and metabolite analysis reveal the drought tolerance of foxtail millet significantly correlated with phenylpropanoids-related pathways during germination process under PEG stress

Background

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.

Results

水吸收的谷子种子模式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.

Conclusions

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 [2］.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 [4］.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 [5］.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 [5］.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 [10］.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 [11］.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 [11］.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 [14］.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 [19］.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 [23］.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 [24］.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 file1: Figure S1 and Additional file2: 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 [5］.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 file8: 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 file3:图S3和表3）.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 file4: 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）.使用KEGG,大大丰富了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.2and Additional file: Figure S6and 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 [25］.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）.两个外信封孔隙S-adenos蛋白质和2ylmethionine 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 [28］.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 [24］.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 [24］.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 [29］.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 (Table1）.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 [23］.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 (Table1and 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 (pvalue 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）.存在表明expressi的结果on 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 [14］.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 [5］.

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 file3: 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 [32］.The changes of transcript levels of phenylpropanoids-related genes were found inRicinus communisunder low temperature [8］.Soybean genes in the phenylpropanoid synthesis pathway were up-regulated in both alkaline and low Fe conditions [33］.

Phenylalanine metabolism provides the precursors for the synthesis of phenylpropanoid [15］.在我们的研究中,参与phenylalanin 17基因e 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 (Table1）.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, EC4.3.1.24), the first enzyme in the phenylpropanoid biosynthesis pathway, catalyzes the conversion of L-phenylalanine to ammonia and trans-cinnamic acid [19］.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 [23］.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 [53］.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 [54］.

五个酚类化合物对微生物的影响ination 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 inArabidopsis thaliana[23］.Ferulic acid, p-coumaric acid and sinapic acid showed inhibitory effect on the germination and radicle length ofA. thaliana. Ferulic acids inhibited total germination ofA. thalianaabove 500 and 750 μM. Wu et al. found that ferulic acids had no effect on annual bluegrass germination [55］.Reigosa et al. reported that lower concentrations of phenolic compounds stimulated the growth of six weeds or had no actions [56］.Ferulic and p-coumaric acids had inhibitory effects on root growth ofA. thalianawithout nutrients [23].And it was found that in tobacco p-coumaric inhibited root growth while ferulic acid had stimulatory effects [57］.These results all indicate that various plant show different response to allelochemicals [58］.

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 [59］.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 [62］.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.

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 [64］.相对发芽率、发芽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

选择的品种是surface-steril的种子ized 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).

Data analysis

核糖核酸测序后,fasta格式的原始数据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 [65］.Read count for each gene was obtained from the mapping results. Gene expression levels were estimated by RSEM for each sample [66］.DESeq R package was used to analyze differential expression of two groups [67］.Pvalue corrections were performed using the Benjamini and Hochberg’s approach. The genes with correctedP-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 [68］.

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 SYBRPremix Ex TaqII (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 SYBRPremix Ex Taq二世(豆类,大连,中国),每个引物(1 1μL0 μ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.β-Actingene of foxtail millet (Seita.7G294000) was used as internal controls for normalization as described by Xu [69］.The relative expression levels of the genes under PEG stress during different stages were calculated via the double standard curve method [70］.

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 [71］.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 [56］.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.

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).

PEG:

Polyethylene glycol

DEGs:

Differentially expressed genes,

FDR:

False discovery rate

COG:

Clusters of orthologous groups

KEGG:

Kyoto encyclopedia of genes and genomes

qRT-PCR:

Quantitative real time polymerase chain reaction

M:

mol/L

cDNA:

Complementary DNA

PAL:

Phenylalanine ammonia lyase

CCR:

Cinnamoyl-CoA reductase

RNA:

Ribonucleic acid

CTAB:

Cetyltrimethylammonium ammonium bromide

FC:

Fold-change

FPKM:

Fragments per kilobase of exon model per million mapped reads

GO:

Gene ontology

PCR:

Polymerase chain reaction

RNA-Seq:

RNA sequencing

HPLC:

High performance liquid chromatography

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.

Author notes

1. Aili Yu and Jinfeng Zhao contributed equally to this work.

Affiliations

1. Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, 046011, China

   Aili Yu, Jinfeng Zhao, Zhenhua Wang, Kai Cheng, Peng Zhang, Gang Tian, Xin Liu, Erhu Guo, Yanwei Du & Yuwen Wang

Contributions

JZ and AY designed the experiments. AY, JZ, ZW, KC, PZ, GT, XL, EG, YD, YW conducted the experiments, analyzed the results and prepared the manuscript. All authors have read and approved the manuscript for publication.

Corresponding author

Correspondence toJinfeng Zhao.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

关于文书期刊上施普林格自然保持中立isdictional claims in published maps and institutional affiliations.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions