Genetic Regulation of GA Metabolism during Vernalization, Floral Bud Initiation and Development in Pak Choi (Brassica rapa ssp. chinensis Makino)

Pak choi (Brassica rapa ssp. chinensis Makino) is a representative seed vernalization vegetable and premature bolting in spring can cause significant economic loss. Thus, it is critical to elucidate the mechanism of molecular regulation of vernalization and floral bud initiation to prevent premature bolting. Gibberellin (GA) is the key plant hormone involved in regulating plant development. To gain a better understanding of GA metabolism in pak choi, the content of GA in pak choi was measured at different stages of plant development using enzyme-linked immunosorbent assay. The results showed that the GA content increased significantly after low-temperature treatment (4°C) and then decreased rapidly with vegetative growth. During floral bud initiation, the GA content increased rapidly until it peaked upon floral bud differentiation. To elucidate these changes in GA content, the expression of homologous genes encoding enzymes directly involved in GA metabolism were analyzed. The results showed that the changes in the expression of four genes involved in GA synthesis (Bra035120 encoding ent-kaurene synthase, Bra009868 encoding ent-kaurene oxidase, Bra015394 encoding ent-kaurenoic acid oxidase, and Bra013890 encoding GA20-oxidase) were correlated with the changes in GA content. In addition, by comparing the expression of genes involved in GA metabolism at different growth stages, seven differentially expressed genes (Bra005596, Bra009285, Bra022565, Bra008362, Bra033324, Bra010802, and Bra030500) were identified. The differential expression of these genes were directly correlated with changes in GA content, suggesting that these genes were directly related to vernalization, floral bud initiation and development. These results contribute to the understanding of the molecular mechanism of changes in GA content during different developmental phases in pak choi.

Pak choi (Brassica rapa ssp. chinensis Makino) is a representative seed vernalization vegetable and premature bolting in spring can cause significant economic loss. Thus, it is critical to elucidate the mechanism of molecular regulation of vernalization and floral bud initiation to prevent premature bolting. Gibberellin (GA) is the key plant hormone involved in regulating plant development. To gain a better understanding of GA metabolism in pak choi, the content of GA in pak choi was measured at different stages of plant development using enzyme-linked immunosorbent assay. The results showed that the GA content increased significantly after low-temperature treatment (4 • C) and then decreased rapidly with vegetative growth. During floral bud initiation, the GA content increased rapidly until it peaked upon floral bud differentiation. To elucidate these changes in GA content, the expression of homologous genes encoding enzymes directly involved in GA metabolism were analyzed. The results showed that the changes in the expression of four genes involved in GA synthesis (Bra035120 encoding ent-kaurene synthase, Bra009868 encoding ent-kaurene oxidase, Bra015394 encoding ent-kaurenoic acid oxidase, and Bra013890 encoding GA20-oxidase) were correlated with the changes in GA content. In addition, by comparing the expression of genes involved in GA metabolism at different growth stages, seven differentially expressed genes (Bra005596, Bra009285, Bra022565, Bra008362, Bra033324, Bra010802, and Bra030500) were identified. The differential expression of these genes were directly correlated with changes in GA content, suggesting that these genes were directly related to vernalization, floral bud initiation and development. These results contribute to the understanding of the molecular mechanism of changes in GA content during different developmental phases in pak choi.
Keywords: pak choi, vernalization, gibberellins metabolism, expression profiles, gene INTRODUCTION Pak choi (Brassica rapa ssp. chinensis Makino) is a cruciferous vegetable and is representative plant that require vernalization. After exposure to low temperatures for a period of time, the plant blossoms under high temperatures and long days. A cultivar with low chilling requirements bolts easily; thus, premature bolting occurs frequently in spring and causes significant reduction in yield and quality (Liu et al., 2009). Therefore, it is crucial to study the molecular regulation of vernalization to prevent bolting in pak choi and possibly other cruciferous vegetables.
Plant hormones play an important role in vernalization and flowering (Aryal and Ming, 2014). Among the six known types plant hormones, gibberellin (GA) is particularly significant (Galvão et al., 2012). GA can partially replace long-day or low temperature to promote flowering in plants (Pearce and Dubcovsky, 2013;Hu et al., 2016). However, the phytohormone GA was shown to play duel, opposing roles in Arabidopsis. GA promoted the termination of vegetative development, but it inhibited flower formation (Yamaguchi et al., 2014). Treatment with GA has been shown to promote the flowering of pak choi Yu et al., 2009). Therefore, studying the regulation of GA on flowering in pak choi is important.
The genes involved in GA metabolism have been isolated in a few plant species, including grape (Wang et al., 2012), Arabidopsis thaliana (Phillips et al., 1995), tomato (Rebers, 1999), Camellia lipoensis (Xiao et al., 2016), tomato (Shen et al., 2016), and pumpkin (Yamaguchi et al., 1996). Although these genes were found to affect the synthesis of GA, their expression patterns were different. In tomato, the expression of the GA1homologous gene LeCPS increased during the early stage of floral bud development (Rebers, 1999). In pumpkin, GA2 was detected high transcription in the vegetative tissues (Yamaguchi et al., 1996). GA3 from Arabidopsis was expressed in inflorescence tissues (Helliwell et al., 1998). Nty, the GA4 ortholog in tobacco, was expressed in specific parts of the floral buds and floral organs (e.g., tapetum) (Yamaguchi et al., 1996;Itoh et al., 1999). GA5 was expressed in the stem and inflorescence of Arabidopsis during pod development (Phillips et al., 1995). These findings indicated that GA plays an important role in different stages of plant flowering. However, the process of GA metabolism during vernalization, floral bud initiation and development in pak choi remains poorly understood. In this study, the GA content in pak choi was measured at different growth stages, and the expression of gene encoding enzymes directly involved in GA metabolism was analyzed, which elucidated the molecular regulation mechanism of GA metabolism during vernalization, floral bud initiation and development in pak choi.

Plant Materials and Growth Conditions
Early-bolting pak choi inbred line "75 # " was used in this study and was provided by the Institute of Vegetable Research, Shanxi Academy of Agricultural Sciences. The experiment included two treatments: low-temperature treatment and control. For the former (designated as V), the germinating seeds were kept at 4 • C for 20 days. In the control (designated as CK), the seeds were germinated at 25 • C. Seedlings with the same size were transplanted to trays with 50 holes full of substrate at the same time. The plants were then cultivated using traditional methods.

Sample Collection
After low-temperature treatment, the shoot apices were collected and designated as V0 and CK0 for the low-temperature treatment and control groups, respectively. On 10 days (vegetative growth), 15 days (immediately differentiation) and 16 days (floral bud differentiation stage 1) after transplantation, shoot apices were also collected from both groups, the low-temperature treatment and control samples collected at these times were named as V10, CK10, V15, CK15, V16, and CK16, respectively. Each 0.2 g sample was used to determine the GA content with three biological replicates. Each 0.1 g sample was used for RNA extraction. All samples were frozen in liquid nitrogen and stored at −80 • C.

Measurement of GA Content
The GA content at different developmental stages of pak choi were measured by enzyme-linked immunosorbent assay using the method described by Deng et al. (2008). Two types of GA (bioactive GA 1 and GA 3 ) were measured.
Total RNA Extraction and cDNA Library Construction, Sequencing, Gene Expression Analysis, and Functional Annotation Total RNA was extracted using an RNeasy Plant Mini Kit according to the manufacturer's instructions (QIAGEN, 74903). The samples used for sequencing were V0, CK0, V10, V15, CK16, and V16. All reads for each sequencing sample were deposited into GenBank, and the SRA accession numbers are SRP111574 (for V10), SRP064332 (for V0, CK0, V15, and CK16) (Sun et al., 2015) and SRP075755 (for V16) (Song et al., 2017; http://www.ncbi.nlm.nih.gov/sra). cDNA library construction and sequencing were performed by Biomarker Technologies Co., Ltd, Beijing, China. The methods of sequencing, expression analysis and functional annotation were previously described by Sun et al. (2015).

Quantitative Real-Time PCR
To validate the RNA sequencing results, quantitative real-time PCR (qRT-PCR) was performed using gene-specific primers for four randomly selected genes. Primer3 software was used to design specific primers. The ACTIN gene of pak choi was used as a reference gene. The RNA was extracted using an RNeasy Plant Mini Kit (QIAGEN, 74903). First-strand cDNA was obtained using a PrimeScript R RT reagent kit (Perfect Real Time; TaKaRa, RR037A). qRT-PCR was carried out using SYBR R Premix Ex TaqTM II (Tli RNaseH Plus; TaKaRa, RR820A) with an ABI 7500 instrument. The 25-µl reaction mixture included 20 ng cDNA. The program was as follows: 94 • C for 1 min; and 40 cycles of 94 • C for 30 s, 55 • C for 30 s, and 72 • C for 30 s. The relative expression levels were calculated using the 2 − Ct method (Stanko et al., 2014).

GA Contents of Pak Choi Shoot Apices at Different Growth Stages
The GA contents of the pak choi shoot apices at different growth stages were measured (Figure 1). After low-temperature treatment (V0), the GA content was 7.36 ng/g·FW, then decreased rapidly with vegetative growth and was approximately 6.80 ng/g·FW 10 d after transplanting (V10). Before floral bud differentiation (V15), the GA content increased rapidly to 9.34 ng/g·FW. Subsequently, the GA content increased slowly and reached 9.86 ng/g·FW at floral bud differentiation stage 1 (V16). Thus, the GA content first decreased and then increased after transplantation, peaking at floral bud differentiation stage 1. However, the GA content of the shoot apices in the control group did not show the same trend. In addition, the GA content in the low-temperature treatment group was significantly higher than in the control group, with the exception of 10 d after transplanting (V10 and CK10). This suggested that GA accumulation in the shoot apices can initiate floral bud differentiation, and low-temperature treatment can increase GA content.

Quality Assessment of the Sequencing Results
The expression profiles of genes in the pak choi shoot apices at different growth stages were analyzed by RNA sequencing. As shown in Table 1, a total of 11,111,314, 11,268,910, 10,849,774, 11,072,438, 10,709,754, and 10,169,936 reads were produced from the V0, CK0, V10, V15, V16, and CK16 samples, respectively; clean reads (those remaining after filtering out the adaptor sequences, contaminating sequences and low-quality reads) accounted for 87. 09, 99.86, 99.79, 99.90, 99.86, and 99.86% of the total reads, respectively, indicating that the sequencing quality was excellent.
To obtain gene expression information, high-quality clean reads from six libraries were aligned with the reference genome database. A total of 82. 94, 84.80, 84.57, 85.08, 84.25, and 84.77% of the reads, including unique mapped reads and multiple mapped reads, were mapped to the reference genome for samples V0, CK0, V10, V15, V16, and CK16, respectively. The proportion

Expression Analysis of Genes Encoding GA-Metabolizing Enzymes
Unlike the majority of GA, which is inactive, GA 1 and GA 3 have biological activity (Pan, 2008). GA is mainly synthesized in the stems and roots of plants, and a small amount is synthesized in mature leaves. It is not transported to other tissues. In cells, GA is synthesized in the plastids, endoplasmic reticulum and cytoplasmic matrix (Reinecke and Pharis, 2013). The KEGG pathway map of diterpenoid biosynthesis (ko00904) from RNA sequencing was shown in Supplementary Figure S1, in which the GA metabolism is marked in red. GA is synthesized in three steps from geranylgeranyl pyrophospha (GGPP) as a precursor. The first step occurs in the protoplast. GGPP forms copal pyrophosphate under the catalysis of CPS (5.5.1.13) and then forms entkaurene with catalysis by KS (4.2.3.19). The second step takes place in the endoplasmic reticulum with O 2 as the substrate, NADPH as the cofactor (Rademacher, 2003), and KO (1.14.13.78) as the catalyst. Ent-kaurene undergoes constant oxidation, generating ent-Kaur-16-en-19-ol, ent-Kaur-16-en-19-al, and ent-Kaur-16-en-19-oate. Finally, GA12-aldehyde is produced under the catalysis of KAO (1.14.13.79). The third step is carried out in the cytoplasm. GA12-aldehyde forms different structures of GA under the catalysis of different enzymes. During this process, under the action of GA20ox (1.14.11.12), GA3ox (1.14.11.15) and cytochrome P450 monooxygenase (P450-1, 2, 3), GA12-aldehyde generates bioactive GA 1 and GA 3 . Under the action of GA2ox (1.14.11.13), GA 1 undergoes 2β -hydroxylation to form inactive GA 8 and GA 34 . GA20ox and GA2ox play key roles in the entire process of GA metabolism. Some genes encoding enzymes related to GA metabolism in pak choi were listed in Table 2. Seven enzymes involving 18 genes were found to participate in GA synthesis. The enzyme involved in the catabolism of GA is mainly GA2-oxidase, 11 genes were found encoding the enzyme in pak choi.
Based on the results of the gene expression profiles, the expression of genes encoding enzymes involved in GA metabolism were analyzed. The results showed that the expression of Bra000864 encoding CPS, Bra035120 encoding KS, Bra009868 encoding KO, Bra015394 encoding KAO, Bra019165 and Bra013890 encoding GA20ox were all upregulated in CK0 vs. V0, V10 vs. V15, and CK16 vs. V16; that is, the gene expression in cold treatment sample were higher than the control after cold treatment and at floral bud differentiation stage 1. Furthermore, the gene expression in shoot apices was higher immediately prior to floral bud differentiation than in the vegetative growth stage. The expression of Bra009285 encoding GA20ox was up-regulated after low-temperature treatment and immediately prior to floral bud differentiation; however, its expression was not detected in shoot apices at floral bud differentiation stage 1, indicating that low temperature promoted its expression. Bra026757 encoding GA3ox was expressed only in immediately prior to floral bud differentiation; it was not expressed in CK0, V0, V16 or CK16. Thus, the expression of Bra026757 was lower in shoot apices at the vegetative growth stage than in the shoot apices immediately prior to floral bud differentiation.
The expression of some genes that participate in GA biosynthesis were irregular. For example, the expression of Bra036239 encoding CPS, Bra005596 encoding KAO and Bra008480 encoding GA3ox were up-regulated after lowtemperature treatment, indicating that the low temperature promoted their expression. Bra028706 and Bra022565 encoding GA20ox were up-regulated after low-temperature treatment and during floral bud differentiation; that is, their expression were higher in the low-temperature treatment group compared to in the control after treatment and at floral bud differentiation stage 1. Bra028277 encoding GA20ox and Bra026122 encoding GA3ox were up-regulated in floral bud differentiation and down-regulated in CK0 vs. V0 and V10 vs. V15, indicating that their expression at floral bud differentiation stage 1 were greater than in the control. These results indicated that the high expression of Bra028706, Bra022565, Bra028277, and Bra026122 are beneficial for the flowering of pak choi.
The above results showed that the expression of genes involved in GA biosynthesis were up-regulated, indicating that GA synthesis was accelerated. The measured GA contents were higher in V0, V15, and V16 than in CK0, V10, and CK16, respectively, and the changes in gene expression were consistent with the changes in GA content.
Among the genes encoding GA2ox involved in GA catabolism, the changes in gene expression were irregular. For example, the expression of Bra030500, Bra010802, and Bra030187 were down-regulated after low-temperature treatment, indicating that low temperature inhibited their expression. The expression of Bra008362 was down-regulated after low-temperature treatment and in floral bud differentiation. The expression of Bra033324, Bra020890, and Bra038780 were down-regulated in floral bud differentiation but up-regulated in CK0 vs. V0 and V10 vs. V15. These results indicated that the expression of these genes were down-regulated, causing GA to decompose more slowly, thereby increasing GA content. This further explains the change in GA content observed in pak choi.
In addition, the genes Bra027106 encoding GA20ox, Bra008481, and Bra020909 encoding GA3ox involved in GA synthesis, and Bra005415 encoding GA2ox involved in GA decomposition were not expressed during vernalization, floral bud initiation and development.

Identification of Differentially Expressed Genes Related to GA Metabolism
According to the sequencing results, the expression of genes related to GA metabolism were analyzed in different periods ( Table 2). After low-temperature treatment, 12 genes related to GA biosynthesis were up-regulated compared to the control, while two were down-regulated; three genes were up-regulated and four were down-regulated in the catabolic process. When the shoot apices immediately enter the floral bud differentiation, 10 genes related to GA biosynthesis were up-regulated compared to in the vegetative stage, whereas five were down-regulated; four were up-regulated, and six were down-regulated in the catabolic process. At floral bud differentiation stage 1, 10 genes related to GA biosynthesis were up-regulated compared to in the vegetative stage, while three genes were down-regulated; six gene were up-regulated and four were down-regulated in the catabolic process. In summary, after low-temperature treatment, more genes involved in GA synthesis were up-regulated than down-regulated during vegetative growth and at floral bud differentiation stage 1; for the genes taking part in GA catabolism, slightly more genes were down-regulated than up-regulated, except for at floral bud differentiation stage 1. This may have led to the higher GA contents in samples V0, V15, and V16 compared to in samples CK0, V10, and CK16. These results were coincide with the measured GA contents. According to FDR < 0.01 and Fold Change ≥ 2, the differentially expressed genes (DEGs) were identified ( Table 3). The DEGs related to GA metabolism were not found in floral bud differentiation stage 1 (V16) vs. the vegetative shoot apices (CK16). Four DEGs after low-temperature treatment (V0 vs. CK0) and three DEGs between V15 and V10 were identified by gene ontology (GO); all of them were involved in GA biosynthesis (GO:0009686) or catabolism (GO:0045487).
Bra005596, the KAO2 homologous gene, was up-regulated after low-temperature treatment. In the second step of GA biosynthesis, in which kaurene generates ent-Kaur-16-en-19ol, ent-Kaur-16-en-19-al, and ent-Kaur-16-en-19-oate followed  (Helliwell et al., 2001a,b). The subcellular localization of AtKAO1 and AtKAO2 revealed that they located in the endoplasmic reticulum (Helliwell et al., 2001b), where the second stage of GA synthesis is carried out. In addition, AtKAO2 was highly expressed in germinating seeds, flowers and fruit pods (Regnault et al., 2014). The expression of KAO2 was up-regulated in pak choi, which was beneficial for early flowering. Both Bra009285 and Bra022565, which encode GA20ox, were up-regulated after low-temperature treatment. Five homologous genes encode GA20ox in Arabidopsis: AtGA20ox1, AtGA20ox2, AtGA20ox3, AtGA20ox4, and AtGA20ox5. Bra009285 and Bra022565 corresponded to AtGA20ox3 (YAP169) and AtGA20ox2, respectively. AtGA20ox3 was expressed in dry seeds, imbibed seeds and silique (Yan et al., 2014) and was found to aid GA biosynthesis in multiple biological processes (Plackett et al., 2012). GA20ox2 promotes floral bud formation and fruit pod growth. The over-expression of GA20ox can lead to early flowering (Sun, 2008). The up-regulation of Bra009285 and Bra022565 was speculated to contribute to early flowering in pak choi.
Among the genes encoding GA2ox involved in GA catabolism, the GA2ox6 homologous genes Bra030500 and Bra033324 were down-regulated after low-temperature treatment and immediate differentiation, respectively. The ATGA2ox2 and ATGA2ox1 homologous genes were Bra010802 and Bra008362, respectively, which were down-regulated and up-regulated after transplantation, respectively. Active GA2oxidase can be inactivated by 2β -hydroxylation, which is one of the main methods of regulating GA activity in cells (Rieu et al., 2008;Wiesen et al., 2016). In Arabidopsis, five genes encode GA2oxidase: GA2ox1, GA2ox2, GA2ox3, GA2ox4, and GA2ox6. The over-expression of these genes in Arabidopsis causes dwarfing, reduces the GA content, and delays flowering (Wang et al., 2004). Similar results were also obtained in transgenic poplar (Populus L), rice and bean (Sakamoto et al., 2001;Busov et al., 2003;Appleford et al., 2007). In this study, the genes encoding GA2ox expressed were down-regulated after low-temperature treatment, which may have increased the content of active GA. The expression of Bra030500, Bra033324 and Bra010802 were down-regulated, which was consistent with the GA content. However, the expression of Bra008362 was up-regulated, which did not coincide with the GA measurements. Further analysis is required to explain this discrepancy.

qRT-PCR Verification
To verify the reliability of the gene expression profiling, Bra000393, Bra004928, Bra033324, and Bra010802 were selected for qRT-qPCR verification in six samples (CK0, V0, V10, V15, CK16, and V16). The qRT-PCR results were consistent with the sequencing results, confirming the reliability of the sequencing results (Figure 2).

Relationship between GA and Vernalization, Floral Bud Initiation and Development
Hormones and GA in particular play an important role in plant vernalization, floral bud initiation and development. GA has been shown to induce the flowering of long-day plants under shortday conditions (Kinet, 1993;Buchanan et al., 2002). The GA content in germinating pak choi seeds treated at low temperature was approximately three times higher than that of the control (Li et al., 2002a). This suggested that vernalization treatment can promote the biosynthesis of endogenous GA, in agreement with the report of Burn et al. (1993). The results also suggest that vernalization can increase the activities of GA-biosynthetic enzymes and thus increase the content of GA. In broccoli, the GA content after low-temperature treatment was three times that in the control; the GA content then gradually decreased with green vernalization, reached a minimum at the critical period of floral bud differentiation, and then increased sharply (Jiang and Yu, 2004). The GA content was low in the vegetative growth stage of rape, increased rapidly during floral bud differentiation, and decreased sharply after floral bud differentiation (Li, 2010). For Chinese flowering cabbage, the GA content of the stem tip was highest in the vegetative period, decreased when floral bud differentiation was about to start, and then increased slightly when the floral bud began to differentiate (Li et al., 2002b). The variation of GA content in Chinese cabbage during floral bud differentiation was similar to that observed in Chinese flowering cabbage (Li, 2002). In summary, in brassicaceous vegetables, the GA content increases after cold vernalization, decreases when floral bud differentiation is imminent, and then gradually increases after floral bud differentiation. In this study, the GA content of pak choi was significantly higher after low-temperature treatment than in the control. The GA content decreased with vegetative growth and increased after floral bud differentiation, consistent with the results of previous studies. However, right before floral bud differentiation, the GA content was inconsistent. One possible explanation for this is that pak choi needs more GA at the beginning of floral bud differentiation.

Expression of Some Important Genes Involved in GA Metabolism
To further clarify the cause of the variation in GA content in pak choi during different growth stages, the GA metabolic pathway was analyzed. In the synthesis of GA, the conversion of GGPP to ent-kaurene is catalyzed CPS and KS. CPS is the key enzyme in GA biosynthesis. The oxidation of entkaurene to GA12-aldehyde is catalyzed by cytochrome P450 monooxygenase. The first three steps from ent-kaurene to entkaurene acid are catalyzed by KO (Helliwell et al., 1999). The following three steps from ent-kaurene acid to GA12aldehyde are catalyzed by KAO (Helliwell et al., 2001b). Both KO and KAO belong to the cytochrome P450 monooxygenase family. The final step in GA biosynthesis is the formation of bioactive GA catalyzed by GA20-oxidase and GA3-oxidase. GA2oxidase catalyzes the formation of inactive GA. These enzymes play important roles in GA metabolism and affect the level of active GA in plants, thereby affecting plant growth and development.
By analyzing the expression profiles, the changes in the expression of four genes encoding enzymes related to GA metabolism were found to be consistent with the changes in GA content. These genes were Bra035120, Bra009868, Bra015394 and Bra015394, and their Arabidopsis homolog genes encode KS1, KO, KAO1, and GA20ox1, respectively. In the first stage of GA synthesis, KS catalyzes the synthesis of ent-kaurene from copal pyrophosphate. In Arabidopsis, AtKS1 was located in the chloroplast stroma by subcellular localization (Helliwell et al., 2001b). The loss of function of KS will lead to serious plant dwarfing and loss of reproductive capacity (Regnault et al., 2014). In this study, Bra035120 was up-regulated with the growth and development of pak choi, and its expression was consistent with the GA content of the shoot tips, in agreement with previous results. KO and KAO, which are members of the cytochrome P450 family of enzymes, participated in the second phase of GA biosynthesis. The process from ent-kaurene to ent-kaurene acid is catalyzed by AtKO1, which belongs to the cytochrome P450 CYP701A subfamily (Regnault et al., 2014). The up-regulation of AtKO1 in transgenic Arabidopsis resulted in an increase in GA 4 content . KAO1 is highly expressed in the seeds and flowers of Arabidopsis, while its expression is low in seedlings and vegetative tissues . Bra009868 and Bra015394 have been suggested to be beneficial for the synthesis of active GA, which agrees with the results of this study. In Arabidopsis, the homologous gene of Bra013890 is GA20ox1, which belongs to the GA20ox family and is highly expressed in stems and flowers (Yan et al., 2014). GA20ox1 plays an important role in the internode elongation of stamen filaments (Sun, 2008). The silencing of SlGA20ox1 does not contribute to pollen production in tomato (Olimpieri et al., 2011). Thus, Bra035120, Bra009868, Bra015394, and Bra013890 are thought to play important roles in the vernalization, floral bud initiation and development of pak choi.
Seven DEGs related to GA metabolism (Bra005596, Bra009285, Bra030500, Bra022565, Bra033324, Bra010802 and Bra008362) were identified by analyzing expression profiles. In Arabidopsis, their respective homologous genes are KAO2, YAP169 (GA20ox3), GA2ox6, GA20ox2, GA2ox6, ATGA2OX2, and ATGA2OX1, respectively. KAO mutations have been shown to cause severe stunting in a variety of plants including barley, rice, corn, peas and sunflowers (Helentjaris, 1995;Davidson et al., 2003;Sakamoto et al., 2004;Fambrini et al., 2011). KAO is thought to play an important role in the synthesis of GAs. YAP169 (GA20ox3) and GA20ox2 are members of the GA20ox family. The over-expression of GA20ox leads to early flowering and stem elongation (Sun, 2008). GA2ox6, ATGA2OX2, and ATGA2OX1 encode GA2ox, and the over-expression of GA2ox can delay plant flowering and cause dwarfing (Sun, 2008). The over-expression of SlGA2ox1 leads to reduction in endogenous GA in tomato (Shen et al., 2016). GA3ox is also a key enzyme in GA biosynthesis, and its mutation leads to dwarfing and flower hypoplasia (Sun, 2008). DEGs encoding GA3-oxidase were not found in this study. It may be that GA20-oxidase and GA2-oxidase play important roles in regulating the level of active GA in plants, while GA3-oxidase has little effect on the late stage of GA biosynthesis in pak choi (Supplementary Figure S1). These results improve our understanding of the molecular flowering mechanism in pak choi.

CONCLUSION
In this study, the GA content of the shoot tips at different stages in pak choi were measured. The results showed that low-temperature treatment increased the GA content, and enhanced GA accumulation initiated floral bud differentiation. To elucidate the molecular mechanism responsible for the changes in GA content, the expression of genes encoding enzymes involved in GA metabolism were analyzed based on the RNA sequencing. The results showed that the expression patterns of most genes involved in GA metabolism, particularly those of four genes (Bra035120, Bra009868, Bra015394, and Bra013890) coincided with the observed changes in GA content. In addition, by analyzing differentially expressed genes at different growth stages, seven genes (Bra005596, Bra009285, Bra022565, Bra008362, Bra033324, Bra010802, and Bra030500) were identified, and their expression were also consistent with the measured GA contents. These results contribute to the understanding of the molecular mechanism of changes in GA content during different developmental phases in pak choi.

AUTHOR CONTRIBUTIONS
XW and MS performed the experiment, data analysis and prepared the manuscript. XQ and AP participated in performed the experiments. LH, GX, and GL participated in the experimental design. ML conceived the idea and participated in the interpretation of results and preparation of manuscript. JZ revised the manuscript. All authors read and approved the final manuscript.