The Coordinated Upregulated Expression of Genes Involved in MEP, Chlorophyll, Carotenoid and Tocopherol Pathways, Mirrored the Corresponding Metabolite Contents in Rice Leaves during De-Etiolation

Light is an essential regulator of many developmental processes in higher plants. We investigated the effect of 4-hydroxy-3-methylbut-2-enyl diphosphate reductase 1/2 genes (OsHDR1/2) and isopentenyl diphosphate isomerase 1/2 genes (OsIPPI1/2) on the biosynthesis of chlorophylls, carotenoids, and phytosterols in 14-day-old etiolated rice (Oyza sativa L.) leaves during de-etiolation. However, little is known about the effect of isoprenoid biosynthesis genes on the corresponding metabolites during the de-etiolation of etiolated rice leaves. The results showed that the levels of α-tocopherol were significantly increased in de-etiolated rice leaves. Similar to 1-deoxy-D-xylulose-5-phosphate synthase 3 gene (OsDXS3), both OsDXS1 and OsDXS2 genes encode functional 1-deoxy-D-xylulose-5-phosphate synthase (DXS) activities. Their expression patterns and the synthesis of chlorophyll, carotenoid, and tocopherol metabolites suggested that OsDXS1 is responsible for the biosynthesis of plastidial isoprenoids in de-etiolated rice leaves. The expression analysis of isoprenoid biosynthesis genes revealed that the coordinated expression of the MEP (2-C-methyl-D-erythritol 4-phosphate) pathway, chlorophyll, carotenoid, and tocopherol pathway genes mirrored the changes in the levels of the corresponding metabolites during de-etiolation. The underpinning mechanistic basis of coordinated light-upregulated gene expression was elucidated during the de-etiolation process, specifically the role of light-responsive cis-regulatory motifs in the promoter region of these genes. In silico promoter analysis showed that the light-responsive cis-regulatory elements presented in all the promoter regions of each light-upregulated gene, providing an important link between observed phenotype during de-etiolation and the molecular machinery controlling expression of these genes.

. Tocopherol composition and content in the etiolated leaves of 14-day-old rice plants during de-etiolation. Values represent the mean of five biological replicates and bars represent standard deviations. Asterisks indicate statistically significant differences compared with the control (0 h illumination) according to the Student's t-test: * p < 0.05.

Functional Analysis of OsDXS Gene Family
In the MEP pathway, pyruvate and glyceraldehyde 3-phosphate are condensed through the action of DXS. DXS is encoded in the rice nuclear genome by three OsDXS gene members [11]. OsDXS3 encoded protein has been demonstrated to have an enzymatic activity that complements the growth of Escherichia coli mutant disrupted in the corresponding DXS gene [12]. A functional E. coli color complementation assay was carried out in a β-carotene accumulating E. coli strain to verify enzymatic activity of OsDXS1 and OsDXS2 encoded proteins. β-Carotene can be synthesized from foreign carotenogenic genes of pACCAR16ΔcrtX which contains all the carotenogenic genes for the generation of β-carotene from farnesyl pyrophosphate (FPP) in non-carotenogenic E. coli as a heterologous host [13]. E. coli cells expressing pACCAR16∆crtX and pUC8 were able to accumulate β-carotene when induced by isopropyl β-D-1-thiogalactopyranoside (IPTG) and showed a light-yellow color (Figure 2A). The color of E. coli strains co-transformed with pACCAR16∆crtX and pUC8-OsDXS1/2/3 were deeper, reflecting the additional flux provided by OsDXS1/2/3 ( Figure 2A). Pigments were extracted from liquid cultures and the β-carotene levels were quantified, revealing that co-transformation with pACCAR16∆crtX and the pUC8 empty vector resulted in a β-carotene yield of 258.99 ± 66.03 μg g −1 DW. The addition of pUC8-OsDXS1 resulted in a 3.2-fold increase to 820.47 ± 81.67 μg g −1 DW, whereas the addition of pUC8-OsDXS2 and pUC8-OsDXS3 resulted in a 2.9-fold increase

Functional Analysis of OsDXS Gene Family
In the MEP pathway, pyruvate and glyceraldehyde 3-phosphate are condensed through the action of DXS. DXS is encoded in the rice nuclear genome by three OsDXS gene members [11]. OsDXS3 encoded protein has been demonstrated to have an enzymatic activity that complements the growth of Escherichia coli mutant disrupted in the corresponding DXS gene [12]. A functional E. coli color complementation assay was carried out in a β-carotene accumulating E. coli strain to verify enzymatic activity of OsDXS1 and OsDXS2 encoded proteins. β-Carotene can be synthesized from foreign carotenogenic genes of pACCAR16∆crtX which contains all the carotenogenic genes for the generation of β-carotene from farnesyl pyrophosphate (FPP) in non-carotenogenic E. coli as a heterologous host [13]. E. coli cells expressing pACCAR16∆crtX and pUC8 were able to accumulate β-carotene when induced by isopropyl β-D-1-thiogalactopyranoside (IPTG) and showed a light-yellow color (Figure 2A). The color of E. coli strains co-transformed with pACCAR16∆crtX and pUC8-OsDXS1/2/3 were deeper, reflecting the additional flux provided by OsDXS1/2/3 ( Figure 2A). Pigments were extracted from liquid cultures and the β-carotene levels were quantified, revealing that co-transformation with pACCAR16∆crtX and the pUC8 empty vector resulted in a β-carotene yield of 258.99 ± 66.03 µg g −1 DW. The addition of pUC8-OsDXS1 resulted in a 3.2-fold increase to 820.47 ± 81.67 µg g −1 DW, whereas the addition of pUC8-OsDXS2 and pUC8-OsDXS3 resulted in a 2.9-fold increase to 750.37 ± 75.36 µg g −1 DW and 3.3-fold increase to 847.37 ± 63.98 µg g −1 DW ( Figure 2B). The results indicated that OsDXS1 and OsDXS2 were similar to OsDXS3, encoding functional DXS activities.

Expression Analysis of Genes Involved in MEP and MVA Pathways
The function of OsDXS1/2/3, OsDXR (DXP reductoisomerase gene), 4-(cytidine 5diphospho)-2-C-methyl-D-erythritol kinase gene (OsCMK), 2-C-methyl-D-erythritol 2,4cyclodiphosphate synthase gene (OsMDS), and OsIPPI1/2 encoded proteins has been experimentally confirmed ( Figure 2) [4,10,12,14]. The expression profiles of MEP and MVA pathway genes were compared to determine whether the light-triggered the isoprenoid biosynthesis gene expression in rice leaves during the de-etiolation process. The changes in expression levels during illumination were investigated by qRT-PCR separately. Figure 3 illustrates the heat map of comparisons of these gene expressions in rice leaves. The genes in MEP and MVA pathways had different responses to light. The majority of the MEP pathway genes were decreased slightly in the first 0.5 h of illumination, then started to increase until peaks, and after that decreased again. There were exceptions such as OsDXS2 and OsDXS3, and the HMBPP synthase gene (OsHDS) among MEP pathway genes. The expression level of OsDXS1 reached the peak after 6 h illumination, which was four-fold greater compared with that in the dark control. Unlike OsDXS1, OsDXS2 and OsDXS3 were downregulated after exposure to light, the highest expression level of OsDXS2 was detected after 1 h (the mRNA level of OsDXS2 was doubled compared with that at 0.5 h), then kept decreasing until 12 h of illumination and slightly increased from 12 to 24 h of illumination. The transcript abundance of OsDXS3 was lower than the other two OsDXS genes in rice leaves and the mRNA level of OsDXS3 after 1 h of illumination dropped to one-tenth of the expression level in dark control. From 1 to 8 h, the expression level of OsDXS3 increased and then decreased sharply ( Figure 3). error bars represent standard deviations. Double asterisks indicate highly significant difference compared with the control (p < 0.01, Student's t-test).

Expression Analysis of Genes Involved in MEP and MVA Pathways
The function of OsDXS1/2/3, OsDXR (DXP reductoisomerase gene), 4-(cytidine 5 diphospho)-2-C-methyl-D-erythritol kinase gene (OsCMK), 2-C-methyl-D-erythritol 2,4 cyclodiphosphate synthase gene (OsMDS), and OsIPPI1/2 encoded proteins has been ex perimentally confirmed ( Figure 2) [4,10,12,14]. The expression profiles of MEP and MVA pathway genes were compared to determine whether the light-triggered the isoprenoi biosynthesis gene expression in rice leaves during the de-etiolation process. The change in expression levels during illumination were investigated by qRT-PCR separately. Figur 3 illustrates the heat map of comparisons of these gene expressions in rice leaves. Th genes in MEP and MVA pathways had different responses to light. The majority of th MEP pathway genes were decreased slightly in the first 0.5 h of illumination, then starte to increase until peaks, and after that decreased again. There were exceptions such a OsDXS2 and OsDXS3, and the HMBPP synthase gene (OsHDS) among MEP pathwa genes. The expression level of OsDXS1 reached the peak after 6 h illumination, which wa four-fold greater compared with that in the dark control. Unlike OsDXS1, OsDXS2 an OsDXS3 were downregulated after exposure to light, the highest expression level o OsDXS2 was detected after 1 h (the mRNA level of OsDXS2 was doubled compared wit that at 0.5 h), then kept decreasing until 12 h of illumination and slightly increased from 12 to 24 h of illumination. The transcript abundance of OsDXS3 was lower than the othe two OsDXS genes in rice leaves and the mRNA level of OsDXS3 after 1 h of illuminatio dropped to one-tenth of the expression level in dark control. From 1 to 8 h, the expressio level of OsDXS3 increased and then decreased sharply ( Figure 3).  The expression level of OsDXR (DXP reductoisomerase gene) rose rapidly after 2 h irradiation and reached a peak after 4 h of illumination, which was 7.6-fold higher compared with the dark control. From 4 to 10 h of illumination, the expression level of OsDXR stayed at a high level and then decreased sharply after 12 h irradiation. The highest expression level of OsMCT (2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase gene) was detected after 6 h of irradiation, which was 12-fold higher than that in dark control. After that, the OsMCT expression level decreased slightly and stayed at a high level from 8 to 12 h. After 24 h of irradiation, the expression level of OsMCT declined rapidly, which was 5.8-fold higher than that in the dark control and a half compared with the peak. The expression pattern of OsCMK response to continuous WL (white light) was similar to that of OsMCT. The maximum expression level of OsCMK was detected at 6 h of illumination (10-fold compared with that in dark control). The expression level of OsMDS kept increasing after exposure to WL, and reached the peak after 10 h of illumination, which was 9.9-fold compared to that of the dark control, then decreased to around 6.0-fold after 12 h irradiation. Unlike the rest of the MEP pathway genes, the highest expression level of OsHDS (HMBPP synthase gene) was detected after 24 h of illumination (9.8-fold compared to the dark control). The dynamic changes in mRNA levels of the OsHDR (HMBPP reductase gene) family have been described by Jin et al. (2020) [10].
The mRNA levels of genes in the MVA pathway during the de-etiolation process were opposite to those in the MEP pathway. The expression levels of MVA pathway genes had higher mRNA levels under dark conditions than those after exposure to light. MVA pathway genes had a similar response after exposure to light. The expression levels of the genes at 10 h of illumination were all lower than those in the dark control. The responses of OsHMGS (3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) synthase gene) family genes to WL illumination were not exactly the same. The expression level of OsHMGS1 was increased slightly in the first 2 h of illumination, then decreased sharply from 2 to 6 h, and the expression level of OsHMGS1 was stable at a low level until 10 h (less than half compared with that in dark control). Similar to OsHMGS1, the mRNA level of OsHMGS2 was stable and kept at a low level after 4 h of illumination. The expression level of OsHMGS3 started to decrease after 1 h of illumination, and the lowest expression level was detected after 12 h, which was only 0.08-fold compared with that in the dark control. The OsHMGR (HMG-CoA reductase gene) family also showed distinct reactions to continuous WL illumination. The maximum expression levels of OsHMGR1-3 were detected at 0.5, 2, and 0.5 h of illumination, respectively. The expression level of OsHMGR1 dropped by 50% compared with the dark control after 24 h of illumination. The mRNA level of OsHMGR2 and OsHMGR3 in the dark control was 7.9-fold and 18-fold greater than that after 24 h of illumination, respectively ( Figure 3). In the first 1 h, the expression level of OsMK (mevalonate kinase gene) increased slightly, then decreased until 24 h. The expression level of OsMK in the dark control was 3.4-fold higher than that after 24 h illumination. The mRNA level of OsPMK (phospho-MVA kinase gene) and OsMVD (MVA diphosphate decarboxylase gene) both increased slightly in the first 0.5 h of irradiation, and after that kept decreasing. The expression changes in the OsIPPI family, involved in both MVA and MEP pathways, were described by Jin et al. (2020) [10].

Expression Analysis of Genes Involved in Carotenoid, Tocopherol, and Chlorophyll Biosynthetic Pathways
The function of rice phytoene synthase 2 (OsPSY2), phytoene desaturase (OsPDS), ζcarotene desaturase (OsZDS), and lycopene β-cyclase (OsLYCB) has been demonstrated [15,16]. The expression levels of the selected genes related to the carotenoid biosynthetic pathway are shown in Figure 4. The expression level of both OsPSY1 and OsPSY2 increased immediately after exposure to WL illumination. The expression level of both OsPSY1 and OsPSY2 reached high levels after 8 and 10 h of illumination (approximately 5-and 25-fold, respectively, greater than the expression level in the dark control), respectively, but both decreased sharply after 10 h of illumination. The expression of OsPSY1 remained nearly constant from 12 to 24 h of illumination, while the expression of OsPSY2 increased rapidly again after 12 h of illumination and reached its maximum after 24 h. The expression level of OsPDS decreased slightly during the first 1 h of illumination, then increased to the peak after 8 h of illumination (4.1-fold compared with that in the dark control). After the peak, the expression level of OsPDS declined until 24 h, when its expression level was 2.4-fold higher than that of the dark control. The expression level of OsZDS kept increasing in the first 12 h of illumination and then declined at 24 h (8.5-and 4.5-fold higher than the dark control). The mRNA levels of both OsLYCB and OsLYCE (lycopene ε-cyclase gene) reached their peaks after 8 h of illumination (4.7-and 2.3-fold compared to the dark control, respectively). Unlike the other genes, OsBCH2 (β-carotene hydroxylase 2 gene) showed a decrease in expression level after exposure to WL, which decreased to about 50% of that of the dark control after 24 h of illumination ( Figure 4). increased rapidly again after 12 h of illumination and reached its maximum after 24 h. The expression level of OsPDS decreased slightly during the first 1 h of illumination, then increased to the peak after 8 h of illumination (4.1-fold compared with that in the dark control). After the peak, the expression level of OsPDS declined until 24 h, when its expression level was 2.4-fold higher than that of the dark control. The expression level of OsZDS kept increasing in the first 12 h of illumination and then declined at 24 h (8.5-and 4.5-fold higher than the dark control). The mRNA levels of both OsLYCB and OsLYCE (lycopene ε-cyclase gene) reached their peaks after 8 h of illumination (4.7-and 2.3-fold compared to the dark control, respectively). Unlike the other genes, OsBCH2 (β-carotene hydroxylase 2 gene) showed a decrease in expression level after exposure to WL, which decreased to about 50% of that of the dark control after 24 h of illumination ( Figure 4). The key tocochromanol biosynthesis genes in rice have been identified [17], while the function of these genes remained undetermined. The expression levels of the selected tocopherol biosynthetic pathway genes in de-etiolated rice leaves are shown in Figure 5. The expression levels of OsGGR (geranylgeranyl reductase gene), OsVTE5 (phytol kinase The key tocochromanol biosynthesis genes in rice have been identified [17], while the function of these genes remained undetermined. The expression levels of the selected tocopherol biosynthetic pathway genes in de-etiolated rice leaves are shown in Figure 5. The expression levels of OsGGR (geranylgeranyl reductase gene), OsVTE5 (phytol kinase gene), OsVTE2 (homogentisate phytyltransferase gene), OsVTE1 (tocopherol cyclase gene), and OsVTE4 (γ-tocopherol methyltransferase gene) were all upregulated by WL. The expression levels of OsGGR increased rapidly during the first 8 h of illumination (10-fold greater than that of dark control) and then plummeted after 12 h of irradiation. The expression level of OsGGR mRNA remained on the rise from 12 to 24 h of illumination. The mRNA of OsVTE1 remained at low levels during the first 2 h of illumination, then started to increase to its maximum after 12 h of illumination (4.9-fold higher than that of dark control), and remained at this high level until 24 h irradiation. There were no significant differences in expression levels among 10-, 12-, and 24-h illuminated samples. Both OsVTE5 and OsVTE2 mRNA levels exhibited a downward trend at the beginning, then started to increase and reached a peak respectively at 8 and 12 h (2.5-and 2.0-fold higher than that in the dark control, respectively). Afterward, the expression levels started to decline, OsVTE5 mRNA level decreased sharply after 24 h of illumination, and the OsVTE2 mRNA level dropped to a level similar to that in the dark control after 24 h irradiation. The expression level of OsVTE4 increased rapidly during the first 4 h of illumination and decreased until 24 h of illumination. However, OsHGGT (homogentisate geranylgeranyltransferase gene) was downregulated by WL. The expression level decreased to three-fold greater than the dark control after 24 h of illumination. gene), OsVTE2 (homogentisate phytyltransferase gene), OsVTE1 (tocopherol cyclase gene), and OsVTE4 (γ-tocopherol methyltransferase gene) were all upregulated by WL. The expression levels of OsGGR increased rapidly during the first 8 h of illumination (10fold greater than that of dark control) and then plummeted after 12 h of irradiation. The expression level of OsGGR mRNA remained on the rise from 12 to 24 h of illumination. The mRNA of OsVTE1 remained at low levels during the first 2 h of illumination, then started to increase to its maximum after 12 h of illumination (4.9-fold higher than that of dark control), and remained at this high level until 24 h irradiation. There were no significant differences in expression levels among 10-, 12-, and 24-h illuminated samples. Both OsVTE5 and OsVTE2 mRNA levels exhibited a downward trend at the beginning, then started to increase and reached a peak respectively at 8 and 12 h (2.5-and 2.0-fold higher than that in the dark control, respectively). Afterward, the expression levels started to decline, OsVTE5 mRNA level decreased sharply after 24 h of illumination, and the OsVTE2 mRNA level dropped to a level similar to that in the dark control after 24 h irradiation. The expression level of OsVTE4 increased rapidly during the first 4 h of illumination and decreased until 24 h of illumination. However, OsHGGT (homogentisate geranylgeranyltransferase gene) was downregulated by WL. The expression level decreased to three-fold greater than the dark control after 24 h of illumination. Figure 5. Relative expression levels of tocopherol pathway genes in etiolated leaves during de-etiolation at different times after the onset of irradiation with white light. Abbreviations: DMGGBQ, 2,3-dimethyl-6-geranylgeranylbenzoquinol; DMPBQ, 2,3-dimethyl-6-phytylbenzoquinol; MGGBQ, 2-methyl-6-geranylgeranylbenzoquinol; MPBQ, 2-methyl-6-phytylbenzoquinol; OsGGR, geranylgeranyl reductase gene; OsHGGT, rice homogentisate geranylgeranyltransferase gene; OsVTE1, tocopherol cyclase gene; OsVTE2, homogentisate phytyltransferase gene; OsVTE4, γ-tocopherol methyltransferase gene; OsVTE5, phytol kinase gene; PPP, phytyl pyrophosphate; α-Toc, α-Tocopherol; β-Toc, β-Tocopherol; γ-Toc, γ-Tocopherol; δ-Toc, δ-Tocopherol; α-T3, α-Tocotrienol; β-T3, β-Tocotrienol; γ-T3, γ-Tocotrienol; δ-T3, δ-Tocotrienol.
So far, the function of a few chlorophyll biosynthesis genes (rice chlorophyll synthase gene (OsCHLG) and rice chlorophyll a oxygenase 1 gene (OsCAO1)) has been demonstrated [4,18,19]. Changes in the expression levels of the selected chlorophyll biosynthetic pathway genes during rice de-etiolation are shown in Figure 6. The expression levels of OsPORA (protochlorophyllide oxidoreductase A gene) and OsPORB reacted differently to Plants 2021, 10, 1456 9 of 16 WL illumination. The OsPORB mRNA level increased significantly in the first 8 h (15-fold higher than that in dark control) and decreased rapidly after 24 h of illumination. The expression level at 24 h was approximately three-fold greater than in the dark control. However, the expression of OsPORA started to decline after exposure to WL, indicating that OsPORA was downregulated by light. The expression of OsCHLG increased in the first 0.5 h of illumination, decreased slightly after 2 h of illumination, increased again after 4 h of irradiation, and remained stable from 8 to 12 h. The maximum expression level was 4.5-fold higher than that in the dark control, which was detected after 12 h irradiation. After that, OsCHLG expression decreased immediately to the same level as that in the dark control. The expression level of OsCAO1 reached a peak after 6 h of illumination (8.4-fold), decreased until 12 h, and slightly increased again from 12 to 24 h. pathway genes during rice de-etiolation are shown in Figure 6. The expression lev OsPORA (protochlorophyllide oxidoreductase A gene) and OsPORB reacted diffe to WL illumination. The OsPORB mRNA level increased significantly in the first 8 fold higher than that in dark control) and decreased rapidly after 24 h of illuminatio expression level at 24 h was approximately three-fold greater than in the dark co However, the expression of OsPORA started to decline after exposure to WL, indi that OsPORA was downregulated by light. The expression of OsCHLG increased first 0.5 h of illumination, decreased slightly after 2 h of illumination, increased agai 4 h of irradiation, and remained stable from 8 to 12 h. The maximum expression lev 4.5-fold higher than that in the dark control, which was detected after 12 h irrad After that, OsCHLG expression decreased immediately to the same level as that in th control. The expression level of OsCAO1 reached a peak after 6 h of illumination (8.4 decreased until 12 h, and slightly increased again from 12 to 24 h.

In Silico Promoter Analysis of Light-Responsive Cis-Regulatory Elements in the Prom Regions of Light-Upregulated Genes
The full-length CDS (coding sequence) of each rice (O. sativa cultivar: Nippo gene was used to identify its corresponding promoter region from NCBI. The mRNA the longest 5′-untranslated region and its corresponding promoter sequence (up to upstream of transcription start site) of each upregulated gene during de-etiolation olated rice leaves were retrieved from NCBI ( Table 1). The online tools PlantCare [2 PLACE [21] were used for the identification of light-responsive cis-acting regulato ments in the promoter of each gene. The promoters of the upregulated genes (Os

In Silico Promoter Analysis of Light-Responsive Cis-Regulatory Elements in the Promoter Regions of Light-Upregulated Genes
The full-length CDS (coding sequence) of each rice (O. sativa cultivar: Nipponbare) gene was used to identify its corresponding promoter region from NCBI. The mRNA with the longest 5 -untranslated region and its corresponding promoter sequence (up to 1.0 kb upstream of transcription start site) of each upregulated gene during de-etiolation of etiolated rice leaves were retrieved from NCBI ( Table 1). The online tools PlantCare [20] and PLACE [21] were used for the identification of light-responsive cis-acting regulatory elements in the promoter of each gene. The promoters of the upregulated genes (OsDXS1,  Table 1). Table 1. In-silico analysis of light-responsive cis-acting regulatory elements of upregulated genes during de-etiolation of etiolated rice leaves. The positions of the cis-acting regulatory elements are relative to those of cDNA (+1 is the first nucleotide of cDNA).

GenBank Accession Numbers
OsZDS The analysis of Arabidopsis revealed that the genes involved in the MVA and MEP pathways behaved differently in response to light stimulation. The expression levels of MVA pathway genes were inhibited by light, whereas those of MEP pathway genes were upregulated [8]. In rice, Kim et al. (2005) [11] reported that, in the case of light illumination, only the relative expression of OsDXS1 increased relative to the unilluminated control. OsDXS3 has been proposed to be specifically involved in the synthesis of diterpenoid phytoalexins (such as momilactones and phytocassanes) in suspension-cultured rice cells after treatment with a chitin elicitor [12,22].
Carotenoids, tocopherols, and chlorophylls are synthesized primarily from precursors derived from the MEP pathway [8]. The number of carotenoids and chlorophylls increased significantly in response to the illumination of etiolated rice seedlings compared with that of the dark control [10]. The quantity of α-tocopherol significantly increased (Figure 1). The first step in the MEP pathway is catalyzed by DXS. There are three OsDXS genes in the rice genome. The expression of OsDXS1 increased four-fold after 4 h of illumination, whereas the expression levels of OsDXS2 and OsDXS3 decreased after exposure to light (Figure 3). Similar to OsDXS3 [12], both the encoded OsDXS1 and OsDXS2 are shown to be functional enzymes (Figure 2) in engineered E. coli. The appearance time of MEP-derived isoprenoids in the de-etiolated rice leaves coincided with the time of the positive response of OsDXS1 to light (Figure 3; Supplementary Figure S1), suggesting that the biosynthesis of chlorophylls, carotenoids, and tocopherols in the de-etiolated rice leaves is mainly generated by OsDXS1. The rest of the MEP pathway gene (OsDXR, OsMCT, OsCMK, OsMDS, OsHDS, and OsHDR1) expression (Figure 3; Jin et al., 2020 [10] for OsHDR1) had positive correlations with light illumination correlated with the MEP pathway-derived carotenoid, chlorophyll, and tocopherol accumulation. A study by Huang et al. (2018) reported photosynthetic pigments in an OsMDS mutant (IspF), which displayed a yellowgreen leaf phenotype, and the chlorophyll and carotenoid levels were significantly lower than those in WT rice. Therefore, OsMDS is essential for the biosynthesis of chlorophylls and carotenoids [14]. The genes encoding enzymes involved in chlorophyll biosynthesis, such as OsCHLG, OsPORB, and OsCAO1, were upregulated by light, except OsPORA ( Figure 6), which was associated with the accumulation of chlorophylls.
Phytoene synthase (PSY) is the rate-limiting enzyme of carotenogenesis. In Arabidopsis seedlings, the expression level of PSY increased under the continuous illumination of both FR and red light [23]. In white mustard (Sinapis alba L.), the upregulation of PSY expression was observed after light exposure, while PDS and GGPS (geranylgeranyl pyrophosphate synthase gene) expression levels remained constant [24]. Both rice OsPSY1 and OsPSY2 are subject to phytochrome and involved in carotenoid biosynthesis in photosynthetically active tissues with different expression patterns during chloroplast development [25]. Carotenoid biosynthesis genes exhibit high expression in light-grown rice tissues [26]. In agreement with reports by Welsch et al. (2008) and Chaudhary et al. (2010) [25,26], the expression levels of OsPSY1/2, OsPDS, OsZDS, OsLYCB, and OsLYCE increased after exposure to WL in this study, whereas the OsBCH2 expression level decreased relative to that of the dark control (Figure 4), reflecting the approximately two-fold increases in the abundance of β-carotene after 24 h of illumination compared to that of the dark control [10].
Studies on Arabidopsis found that AtVTE2 showed limited activities in tocopherol synthesis [27]. OsVTE2 was upregulated by light. Both OsVTE1 and OsVTE5 were also upregulated by WL in de-etiolated rice leaves ( Figure 5). Meanwhile, OsVTE4 increased rapidly in the early stage of rice de-etiolation. The upregulation of tocopherol biosynthetic genes in de-etiolated rice leaves was reflected by the increase in tocopherol content in rice leaves during de-etiolation.
Phytosterols are derived from the MVA pathway. In a previous study [10], the level of β-sitosterol increased significantly in the first 6 h of light illumination and then decreased after 12 h of illumination to the same level as that of the dark control. After 24 h of light illumination, the levels of campesterol and stigmasterol showed no difference compared with that of the dark control [10]. The expression levels of the selected MVA pathway genes were inhibited with no obvious changes during the de-etiolated rice leaves (Figure 3), consistent with the MVA pathway-derived phytosterol changes.

The Underpinning Mechanistic Basis of Coordinated Light-Responsive Gene Expression during the De-Etiolation Process
The coordinated expression of light-responsive genes during the de-etiolation can reflect similarities between their promoters. A range of light-responsive elements (LREs) has been documented in different promoters, many of which positively or negatively mediate gene expression in response to light [3]. Although many LREs and their transcription factors have been identified, no single element is found in all light-upregulated gene promoters, suggesting a complex light-regulation network and a lack of a universal switch. It has been suggested that combinations of LREs, rather than individual elements, could confer proper light-responsiveness to a light-insensitive basal promoter [3,28]. To investigate the basis of coordinated light-responsive gene expression during the de-etiolation of rice etiolated leaves, in silico promoter analysis of upregulated genes during de-etiolation was performed. Most promoters contained at least two motifs involved in the response to light, but both OsCMK and OsHDR1 promoters (Table 1) had only one light-responsive motif (IBOXCORE, GATAA), supporting that these genes were upregulated by light. However, surprisingly, except for OsCMT, OsLYCE, and OsVTE2 promoters, other light-upregulated gene promoters all presented a cis-regulatory element with a consensus sequence GATAA (IBOXCORE) ( Table 1). IBOXCORE (I-box) is a conserved sequence upstream of light-upregulated genes for both monocots and dicots [21,28]. (http://www.dna.affrc.go.jp/htdocs/PLACE/; 18 March 2020 (See IBOX, S000199)).

Plant Material
Wild type rice (Oryza sativa L. cv. EYI105) seeds were sterilized in 3.7% NaClO for 20 min and washed with sterilized water thoroughly, then were aseptically grown on Murashige and Skoog (MS) medium containing 0.3% phytagel in pots (diameter 10 cm, height 20 cm, 9 seeds in each pot) in a controlled growth chamber at 25 • C in complete darkness. Fourteen-day-old etiolated rice seedlings were illuminated with continuous WL (100 µmol m −2 s −1 ) for 0.5, 1, 2, 4, 6, 8, 10, 12, or 24 h (hour). Leaves of de-etiolated rice seedlings were rapidly frozen in liquid nitrogen after collection and stored at −80 • C. The mixture of all the rice leaves from the same pot was collected and stored as one sample.

Tocochromanols Analysis
Rice leaf samples were freeze-dried before analysis, and 10 mg of powdered tissue was extracted using a chloroform-methanol mixture, as described [10]. The gas chromatographymass spectrometry (GC-MS) method was applied to determine tocochromanols as described in detail previously [29]. Compositions of tocochromanols were identified using a mass spectral library built from in-house standards [30]. The samples used for tocochromanols analysis were from the same batch [10] as those used for carotenoids, chlorophylls, phytosterols, and pathway gene expression analysis.

Rice DXS cDNA Cloning and Functional Complementation Analysis
The OsDXS1, OsDXS2, and OsDXS3 cDNAs were cloned directly from rice leaf mRNA by RT-PCR using two pairs of oligonucleotide primers with terminal EcoRI and BamHI restriction sites for both OsDXS1 and OsDXS3, or EcoRI and HindIII restriction sites for OsDXS2 (Supplementary Table S2). The primers were designed to flank the open reading frame (ORF) of each gene according to sequences in GenBank (accession number XM_015785019 for OsDXS1, XM_015787004 for OsDXS2, and XM_015792458 for OsDXS3). The RT-PCR products were transferred to vector pCRII TOPO (Invitrogen, Carlsbad, CA, USA) to generate pCR-OsDXS1, pCR-OsDXS2, and pCROsDXS3 for sequencing. The Os-DXS1, OsDXS2, and OsDXS3 cDNAs were then introduced into pUC8 at the EcoRI/BamHI or EcoRI/HindIII sites to generate vectors pUC8-OsDXS1, pUC8-OsDXS2, and pUC8-OsDXS3 for functional complementation in engineered E. coli cells. The integrity of all intermediate and final constructs was confirmed by sequencing.

Gene Expression Analysis with Quantitative Real-Time PCR (qRT-PCR)
Total RNA was isolated using the RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) and DNA was removed with DNase I (RNase-Free DNase Set, Qiagen) [10]. First-strand cDNA was synthesized from 1 µg total RNA using Ominiscript Reverse Transcriptase (Qiagen) in a 20 µL total reaction, and quantitative real-time RT-PCR (qRT-PCR) was applied [10,32] using the primers listed in Supplementary Table S1. The amplified DNA fragments for each gene were confirmed by sequencing. The expression levels of the genes in different samples were normalized against rice actin mRNA (OsActin). Three biological replicates each comprising three technical replicates were tested for each sample.

Data and Statistical Analysis
Significant differences in quantitative dependent variables for any two alternative treatments were calculated using Student's t-test, with a p-value of 0.05 considered significant.

Conclusions
Both OsDXS1 and OsDXS2 genes were shown to encode functional DXS enzymes. The appearance of MEP-derived plastidial isoprenoids in the de-etiolated rice leaves correlates well with the observed upregulation of OsDXS1 expression, suggesting that OsDXS1 is responsible for the biosynthesis of chlorophylls, carotenoids, and tocopherols in the deetiolated rice leaves. The coordinated upregulated expression of genes involved in the MEP pathway, chlorophyll, carotenoid, and tocopherol biosynthetic pathways mirrored the increased levels of chlorophyll, carotenoid, and tocopherol contents in rice leaves during de-etiolation, demonstrating that the synthesis of these plastidial metabolites is under transcriptional control. Moreover, the work was extended to elucidate the underpinning mechanistic basis of coordinated light-responsive gene expression during the de-etiolation process. The light-responsive cis-regulatory elements were presented in all promoter regions (up to 1.0 kb upstream of transcription start site) of each light-upregulated gene, serving as a base for the further characterization of these light-induced genes during the de-etiolation process.