Impact Factor 3.677

The world's most-cited Plant Sciences journal

Original Research ARTICLE

Front. Plant Sci., 22 September 2016 | https://doi.org/10.3389/fpls.2016.01425

Transcriptome Analysis of Sucrose Metabolism during Bulb Swelling and Development in Onion (Allium cepa L.)

  • 1Beijing Vegetable Research Center, Beijing Academy of Agriculture and Forestry Sciences/Key Laboratory of Biology and Genetic Improvement of Horticultural Crops (North China), Beijing, China
  • 2National Key Laboratory of Crop Genetic Improvement, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan, China
  • 3Key Laboratory for Biology of Greenhouse Vegetable of Shandong Province, Institute of Vegetables and Flowers, Shandong Academy of Agricultural Sciences, Jinan, China
  • 4Vegetable Research Laboratory, Institute of Agricultural Sciences of Xuhuai Region in Jiangsu, Lianyungang, China

Allium cepa L. is a widely cultivated and economically significant vegetable crop worldwide, with beneficial dietary and health-related properties, but its sucrose metabolism is still poorly understood. To analyze sucrose metabolism during bulb swelling, and the development of sweet taste in onion, a global transcriptome profile of onion bulbs was undertaken at three different developmental stages, using RNA-seq. A total of 79,376 unigenes, with a mean length of 678 bp, was obtained. In total, 7% of annotated Clusters of Orthologous Groups (COG) were involved in carbohydrate transport and metabolism. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, “starch and sucrose metabolism” (147, 2.40%) constituted the primary metabolism pathway in the integrated library. The expression of sucrose transporter genes was greatest during the early-swelling stage, suggesting that sucrose transporters (SUTs) participated in sucrose metabolism mainly at an early stage of bulb development. A gene-expression analysis of the key enzymes of sucrose metabolism suggested that sucrose synthase, cell wall invertase, and invertase were all likely to participate in the hydrolysis of sucrose, generating glucose, and fructose. In addition, trehalose was hydrolyzed to two molecules of glucose by trehalase. From 15 to 40 days after swelling (DAS), both the glucose and fructose contents of bulbs increased, whereas the sucrose content decreased. The growth rate between 15 and 30 DAS was slower than that between 30 and 40 DAS, suggesting that the latter was a period of rapid expansion. The dataset generated by our transcriptome profiling will provide valuable information for further research.

Introduction

Onion (Allium cepa L.), a group of monocotyledonous biennial herbs belonging to the Alliaceae family, is the most economically important vegetable plant (Jakše and Bohanec, 2003). It may have been the earliest cultivated form of any vegetable crop. Dating back 5000 years, onions were already an important food source in ancient Egypt. With many health-related benefits, onions are frequently recognized as having an important dietary role, especially in preventing cardiovascular disease and cancer (Havey et al., 2004). Onions can be classified as sweet or non-sweet. Their significance in cooking is determined by their taste characteristics (pungent and sweet) and flavor profile. About 80% of onion bulb dry matter consists of non-structural carbohydrates (Darbyshire and Henry, 1981). The main carbohydrate components are glucose, fructose, sucrose, and fructo-oligosaccharides. Glucose, fructose, and sucrose account for 65% of the dry matter content, which varies from ~5% of fresh weight in sweet onions to ~30% in dehydrated varieties (Darbyshire and Henry, 1979; McCallum et al., 2006).

Onion pungency is caused by a range of sulfur compounds. When onions are first cut, some of these compounds affect the eyes and produce tears (Tewari and Bandyopadhyay, 1975). A high degree of pungency can mask a high level of sugar, resulting in the onion not being considered to be sweet. Also, onions with low pungency and low sugar content can be regarded as bland. Ideally, a sweet onion will have a high sugar level and low pungency. Thus, the balance between the pungency and sugar levels determines the perception of sweetness in an onion. In all cases, the sweetness and pungency that are produced are important aspects of the formation and development of onion bulbs. Despite the literature on sucrose metabolism in plants, the genetic mechanisms involved in the formation and development of onion bulbs have not been reported.

The formation and development of onion bulbs are closely related to sucrose metabolism (Sinclair et al., 1995; Mallor et al., 2011). In the non-photosynthetic cells of higher plants, sucrose is transported from the photosynthetic apparatus and cleaved to its constituent monosaccharides, hexoses (Hexs) or phosphorylated Hexs, which can then be used either in catabolic or biosynthetic reactions (Ruan, 2014). The only known enzymatic processes of sucrose (Suc) cleavage in plants are catalyzed by invertases [Suc combines with H2O to generate glucose (Glc) and fructose (Fru)] and sucrose synthases (SuSys) [sucrose combines with uridine diphosphate (UDP) to generate fructose and UDP-glucose (UDP-Glc)] (Koch, 2004). These processes typically degrade sucrose in vivo, but the products of their reactions differ in important ways (Tymowska-Lalanne and Kreis, 1998; Winter and Huber, 2000). Invertases produce glucose instead of UDP-Glc, and thus also form twice as many Hexs. As an important transportation form of photosynthetic products, in common with fumaric acid, sorbitol, mannitol, and raffinose, Suc is moved from leaves to sink tissues through the phloem by sucrose transporters (SUTs) (Riesmeier et al., 1994; Patrick et al., 2013). These are clusters of transmembrane proteins that guide sugars across the membrane to implement the transport between phloem and source/sink cells by an apoplast upload and download pathway. Suc is loaded into the phloem by the respective sugar transporters, the sieve element/companion cell (SE/CC) complex or through plasmodesmata connecting the SE/CC complex. Suc is degraded by four different enzymatic mechanisms (Koch, 2004; Vandeputte and Delcour, 2004; Jiang et al., 2014). First, Suc is hydrolyzed into Hexs (Glc and Fru) by cell wall invertase (CWIN; EC 3.2.1.26) before entering the cytoplasm. Hexs are then transported into the cytoplasm by symporters. Secondly, Suc that is imported through plasmodesmata may be degraded by SuSy, producing UDP-Glc and Fru. Thirdly, Suc taken up the cytoplasm is hydrolyzed into Glc and Fru by a cytoplasmic invertase (CIN). Fourthly, cytosolic Suc may be carried into vacuoles for hydrolysis by vacuolar invertase (VIN). CWIN and VIN activities are limited by their respective inhibitors (INHs). The intracellular Hex is used for glycolysis and in the synthesis of sugar polymers. Additionally, the Hex level may also be sensed by the nucleus-localized hexokinase (HXK) or other proteins to regulate gene expression. Suc synthesis mainly involves a single pathway in which Hexs (Glc and Fru) are converted into hexose-6-phosphates by HXK (EC 2.7.1.1). Fru-6-phosphate and UDP-Glc are converted into Suc-phosphate (Suc-P) by Suc-P synthase (SPS), and then Suc-P is converted into Suc and phosphate by the action of Suc-P phosphatase. At the same time, a small amount of trehalose (T) is synthesized, which is a disaccharide and present in plants at low micromolar concentrations, suggesting a regulatory rather than a metabolic role (Paul et al., 2008); using Glc-6-phosphate and UDP-Glc as substrates, T-6-phosphate (T6P) synthase (TPS) first synthesizes T6P, which is then converted to T by T6P phosphatase (TPP). Finally, T is hydrolyzed by trehalase (Tl) into two molecules of Glc (Figure 1; Cabib and Leloir, 1958; Koch, 2004; O'Hara et al., 2013; Li et al., 2014; Ruan, 2014).

FIGURE 1
www.frontiersin.org

Figure 1. Sucrose (Suc) loading, transport and unloading pathways in a non-photosynthetic cell. SE/CC, sieve element/companion cell; PD, plasmodesmata; Suc, sucrose; Glc, glucose; Fru, fructose; SUT, sucrose transporter; INH, inhibitor; CWIN, cell wall invertase; CIN, cytoplasmic invertase; SPS, sucrose-phosphate synthase; Ti, trehalase; TPP, trehalose-6-phosphate phosphatase; TPS, trehalose-6-phosphate synthase; SuSy, sucrose synthase; HXK, hexokinase; VIN, vacuolar invertase.

Reports on bulbous plants, such as Gladiolus hybridus (He et al., 2008), Tulipa gesneriana (Yu et al., 2009), and Lilium (Zheng et al., 2012; Li et al., 2014), have focused on starch or Suc metabolism. These studies have indicated that starch or Suc metabolism play an important role in the formation and development of bulbs. As with other bulbs, Suc metabolism is crucial in the development of onion (A. cepa L.) bulbs. In the study reported here, RNA-Seq, which is a powerful approach for detecting both differentially expressed genes (DEGs) and novel expressed genes over a broad dynamic range (Blencowe et al., 2009; Wang et al., 2009), has been used to elucidate Suc metabolism in onion, with the following objectives: (i) to identify DEGs involved in the formation and development of onion bulbs; and (ii) to screen the critical genes that are responsible for the changes in Suc, Glc, and Fru metabolism during the swelling of onion bulbs.

Materials and Methods

Plant Material and Sample Collection

The Utah Yellow Sweet Spain cultivar “Y1351” was used in this study. Fresh, undamaged onion (A. cepa L) bulbs were obtained from Beijing Academy of Agriculture and Forestry Vegetable Research Center. The bulbs were monitored until they formed a definite shape (bulbous) and size (7-cm diameter, with weights of 260–300 g), and samples were collected at three time points for cDNA library construction. The number of samples collected was ~10 on the 15th day after swelling (DAS) of bulb (2-cm diameter and 3–4 g weight), ~6 on the 30th DAS (5-cm diameter and 100–110 g weight), and ~3 on the 40th DAS (7-cm diameter and 260–300 gram), with three biological replicates per time point. Fresh samples were washed in running water to remove dirt and excess water was removed with blotting paper. Using a scalpel, 1–2 cm was excised from the top and the bottom of each bulb, and the remainder of the bulb was divided into three portions and used as the experimental material. For real-time quantitative PCR, material from bulbs was harvested at the three different time points, frozen immediately in liquid nitrogen and stored at −80°C prior to use.

RNA Isolation and cDNA Library Construction

Total RNA from each sample was isolated using TRIzol Reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's instructions, with the modification that, to remove the supernatant fluid, precipitation with isovolumetric iso-propanol and high-salt solution (1.2 mol/L sodium chloride and 0.8 mol/L sodium citrate) was used instead of iso-propanol precipitation. The RNA concentration prior to gel electrophoresis was determined using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA); and to confirm the quality of RNA for further processing, the integrity of total RNA samples was determined using a RNA Nanodrop 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). The RIN (RNA integrity number) values of these samples were 9.3–9.6. The construction of the libraries and RNA-seq were performed by the Biomarker Biotechnology Corporation (Beijing, China).

The cDNA libraries from 15 DAS, 30 DAS, and 40 DAS samples were generated using an Illumina TruSeq™ RNA Sample Preparation Kit (Illumina, San Diego, CA, USA). Briefly, mRNA was enriched and purified from 3 μg of total RNA using oligo (dT)-attached magnetic beads and broken into short fragments using Illumina proprietary fragmentation buffer at 94°C for 15 min. Taking these cleaved mRNA fragments as templates, first-strand cDNA was synthesized using random hexamers and SuperScript II. Second-strand cDNA was then synthesized using DNA polymerase I, dNTPs, and RNaseH. These cDNA fragments were then subjected to an end-repair, poly(A) tailing was added, together with a ligation sequencing adapter, and the fragments were then purified and size-selected using AMPure XP beads. Next, PCR amplification was performed to enrich the purified cDNA template. Finally, the nine libraries were sequenced using an Illumina HisSeq™ 2500.

De novo Transcriptome Assembly and Annotation

Using the HisSeq™ 2500, transcriptome sequencing generated 100-bp pair-end (PE) raw reads. For further analysis, the raw reads were filtered by removing the adaptor reads and low-quality reads or unknown nucleotides. The rest of the clean reads of sufficiently high quality were generated and then assembled using the Trinity platform (http://trinityrnaseq.sourceforge.net/) for transcriptome assembly without a reference genome, with the parameters of “K-mer = 25, group pairs distance = 300,” and with other parameters set to default values (Grabherr et al., 2011). For each library, short reads were first assembled into longer contigs based on their overlap regions. The contigs were subjected to clustering to form the components. Different contigs from other transcripts and their distance were then recognized by the approach of De Bruijn graph assembly, and by mapping clean reads back to the corresponding contigs based on pair-end information, and thus the transcripts were produced. Finally, the potential transcript sequences were clustered using the TGI Clustering tool in order to obtain uni-transcripts (Pertea et al., 2003).

The unigene sequences were searched against the following public databases: Non-redundant protein (NR) database (Deng et al., 2006), Swiss-Prot (Apweiler et al., 2004) Gene Ontology (GO) database (Ashburner et al., 2000), Clusters of Orthologous Groups (COG) database (Tatusov et al., 2000), and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG, Kanehisa et al., 2004). If the results from different database searches conflicted, then the order of priority was as follows: NR, KEGG, Swiss-Prot, GO, and COG. The threshold values for NR and Swiss-Prot were cut-off e-values of 10−5 and for COG it was a cut-off e-value of 10−3. The Blast2GO program (Conesa et al., 2005) was used to obtain GO annotation of the unigenes. WEGO and Top GO software (https://www.bioconductor.org/packages/2.12/bioc/html/topGO.html) was then used to perform GO functional classification and enrichment analysis of all the unigenes, in order to view the distribution of gene functions.

Differential Gene Expression Analysis

To compare the abundantly expressed genes in different samples, Bowite software (Langmead et al., 2009) was used to align between unigenes and reads obtained from each sample, and then associated with RSEM (Li and Colin, 2011) to estimate the expression of unigenes. Transcript count information for sequences corresponding to each unigene was calculated and normalized to the fragments per kilobase of transcript per million mapped reads (FPKM; Trapnell et al., 2010) values, in order to eliminate the influence of gene length and sequencing errors on the gene expression calculation. Significant DEGs were screened using DESeq software (Anders and Huber, 2010). The corrected P-values from this method accounted for multiple tests used the key factor, which was false discovery rate (FDR). FDR < 0.01 and |log2 (fold change)| > 1 or < −1 were set as the thresholds for differential gene expression. Fold changes in the expression levels between samples were used as the criteria in the screening process.

Suc, Glc, and Fru Assays

The levels of Suc, Glc, and Fru in bulbs at the three different developmental stages were determined by high-performance liquid chromatography (HPLC). Briefly, 1.0 g fresh onion bulb tissue was lyophilized and extracted in 40.0 mL distilled water by sonication for 30 min. The sonicated solution was then centrifuged at 13,000 × g for 5 min and filtered through a 0.45 μm Millipore filter prior to analysis. The sugar contents at the three different bulb developmental stages were determined in three independent experiments and carried out using three biological replicates.

Quantitative Real-Time PCR Analysis

To confirm RNA-seq results (DEGs selected) and to determine the roles of key enzymes involved in sucrose, glucose and fructose metabolism, quantitative real-time PCR (qRT-PCR) was performed using SYBR® green (TaKaRa, Dalian, China). The 1st-strand cDNA was synthesized using a PrimeScript II 1st Strand cDNA Synthesis Kit (Takara, Dalian, China). β-Actin was used as an internal reference gene. SYBR® green primers for quantitative real-time PCR were designed used Primer Premier 5.0 software. Quantitative real-time PCR was performed on a Roche Light Cycler 480 system (Bio-Rad, USA). The reaction mixture total volume of 20 μL contained 10 μL of 2 × SuperReal PreMix Plus, 2 μL of cDNA mix, 0.6 μL of each primer, and 6.8 μL of RNase-free double-distilled H2O. The program for RT-PCR was 95°C for 15 min, followed by 40 cycles of 95°C for 10 s, 60°C for 20 s, and 72°C for 20 s. A melting curve analysis was performed at the end of each PCR reaction at 95°C for 5 s, 65°C for 60 s, and 97°C continuously, prior to termination at 40°C for 30 s. Relative mRNA levels were calculated using the 2−ΔΔCT method to normalize and calibrate gene expression levels relative to the internal reference β-actin (Livak and Schmittgen, 2001). Genes and primer sequences can be found in Table S1.

Results

Collection of Bulbs

The transcriptome of onion (A. cepa L.) cultivar “Y1351” bulbs was evaluated during development, with reference to characteristic physiological and morphological changes. The first developmental stage, at ~15 DAS, has elliptical-shaped bulbs; the second stage, at ~30 DAS, has suborbicular-shaped bulbs; and the third stage, at 40 DAS, has bulbous-shaped bulbs and is approaching maturity (Figures S1A–C).

RNA-Seq Analysis

A total of 72.53 million clean reads from the samples were obtained following quality assessment and data filtering. The GC% of sequenced data from the three libraries was ~45%, and the percentage of reads with an average quality score >30 was >85%. This indicated that the accuracy and quality of the sequencing data were sufficient for further analysis. The general sequencing statistics are shown in Table S2. The distribution of the expression of the unigenes was similar in the three different samples, which indicated that there was no bias in the process of constructing the cDNA libraries (Figure S2).

Using Trinity, the reads were assembled into 79,376 unigenes with a mean length of 678 bp. The median contig (N50) length was 1093 bp and 28,724 unigenes had sequences of 200–300-bp in length (Figure S3).

Gene Annotation and Function Classification

The BLASTX algorithm was used to query the assembled sequences against Nr, Swiss-Prot, GO, COG, and KEGG databases, resulting in 27,955 (35.22%) perfectly annotated unigenes. Of these, 27,699 (99.08%) of the annotated unigenes were matched to the Nr database, while the smallest number of annotated unigenes was found in the KEGG database (Table S3). To further appraise the completeness of the RNA-seq data, COG classifications were performed with the annotated gene sequences. In total, 8288 unigenes were annotated according to 11,713 functions in 25 COG categories (Figure 2). Among them, “general function prediction only” (2217, 18.93%) ranked highest. Moreover, “carbohydrate transport and metabolism” was represented by only 585 unigenes (i.e., 7% of the annotated COG).

FIGURE 2
www.frontiersin.org

Figure 2. Histogram of clusters of orthologous groups (COG) classification. The 8288 unigenes were grouped into 25 COG categories.

Unigenes that were successfully annotated to the GO database were classified into three principal categories for the three different stages of bulb development, comprising “cellular component,” “molecular function,” and “biological process,” which were further subdivided into 56 categories (Table S4; Figure 3). Of these, the categories that were most represented in the “biological process” principal category were “metabolic process” [15 DAS: 11,737 (64.82%); 30 DAS: 11,522 (64.81%); and 40 DAS: 10,979 (64.5%)], which was the largest group, followed by “cellular process” [15 DAS: 11,084 (61.21%); 30 DAS: 10,894 (61.27%); and 40 DAS: 10,396 (61.08%)] and “response to stimulus” [15 DAS: 4912 (27.13%); 30 DAS: 4869 (27.39%); and 40 DAS: 4737 (27.83%)]. Within the “cellular component” principal category, unigenes [15 DAS: 12,866 (24.56%); 30 DAS: 12,678 (24.54%); 40 DAS: 12,293 (24.52%)] belonged to the “cell” category. In the principal category of “molecular function,” the two categories most represented were “binding” [15 DAS: 8461 (46.73%); 30 DAS: 8216 (46.21%); and 40 DAS: 7734 (45.44%)], and “catalytic activity” [15 DAS: 8339 (46.05%); 30 DAS: 8226 (46.27%); and 40 DAS: 7831 (46.01%)].

FIGURE 3
www.frontiersin.org

Figure 3. Histogram of gene ontology (GO) classification for unigenes derived from bulbs in three developmental stages. The unigenes corresponded to three main categories: “biological process,” “cellular component,” and “molecular function.” The left-hand y-axes indicate the percentage of genes.

Pathway Enrichment Analysis of DEGs

To identify the active pathways in the expansion stage of the onion bulb, the unigenes obtained were mapped to the canonical reference pathways in the KEGG database. KEGG pathway enrichment analysis was performed to categorize the biological functions of DEGs. Specific enrichment of genes was obtained for 95 pathways at 15 DAS vs. 30 DAS, for 79 pathways at 15 DAS vs. 40 DAS and for 86 pathways at 30 DAS vs. 40 DAS. Between library pairs, a total of 6113 unigenes were annotated in 120 KEGG pathways. The top 20 enriched pathways are listed in Table 1. “Ribosome” was the most common term and contained 401 (6.56%) DEGs, followed by “oxidative phosphorylation” (235, 3.84%), “protein processing in endoplasmic reticulum” (226, 3.70%), “plant hormone signal transduction” (202, 3.30%), and “RNA transport” (186, 3.04%). Notably, the other pathways were mainly involved in secondary metabolism, energy metabolism, protein processing, and RNA transport and splicing. But “starch and sucrose metabolism” (147, 2.40%) constituted the primary metabolism pathway in the integrated library and might be more active during bulb swelling and development.

TABLE 1
www.frontiersin.org

Table 1. The top 20 enriched KEGG pathways.

Gene Expression during Bulb Development

Setting a P < 0.01 and |log2 (fold change)| ≥ 1, DEGs involved in the development of bulbs were identified as up-regulated or down-regulated genes. As shown in the Venn diagram, a total of 146 DEGs were detected across all three pairs of developmental stages (Figure 4). The numbers of up- or down-regulated genes are shown in Figure 5. The numbers of up- or down-regulated genes changed markedly between 15 and 30 DAS, (2120 up-regulated vs. 1334 down-regulated genes), while the gene-expression pattern changed steadily between 15 and 40 DAS (1197 up-regulated vs. 1111 down-regulated genes). Between 30 and 40 DAS, 2308 DEGs were detected with 1197 up-regulated and 1111 down-regulated genes (Figure 5; Table S5).

FIGURE 4
www.frontiersin.org

Figure 4. Venn diagram showing genes differentially expressed in any two stages of onion bulb development.

FIGURE 5
www.frontiersin.org

Figure 5. Profile of the number of differentially expressed genes at different developmental stages. The numbers of up- and down-regulated genes in comparisons of 15 DAS vs. 30 DAS, 15 DAS vs. 40 DAS, and 30 DAS vs. 40 DAS bulb samples are shown.

Determination of Levels of Sugars

To investigate the accumulation and metabolism of sugars in swelling onion bulbs, the contents of Fru, Glc, and Suc at the three different developmental stages were measured by HPLC (Figure 6). At 15 DAS, the Glc and Fru contents were 172.93 mg/g fresh weight (fw) and 164 mg/g fw, respectively. The Suc content at this stage was at its maximum value, 120.22 mg/g fw. Interestingly, the Glc content (from 194.92 to 252.87 mg/g fw) in bulbs was higher than that (from 181.87 to 212.21 mg/g fw) of Fru and increased quickly from 30 to 40 DAS. At the same time, the Suc content decreased gradually from 15 DAS to reach 119.51 mg/g fw at 30 DAS and 95.63 mg/g fw at 40 DAS. This indicated that, after 30 DAS, bulbs started swelling rapidly, entering the rapid expansion period. After 40 DAS, the bulb had become essentially fully swollen and the Suc, Glc, and Fru contents barely changed (Grzelak et al., 2009).

FIGURE 6
www.frontiersin.org

Figure 6. Changes in the contents of Suc, Fru, and Glc at three different stages of onion bulb development.

Critical Genes Involved in Suc Metabolism in Onion Bulb Swelling

To gain a deeper insight into Suc metabolism during onion bulb swelling, a list of 147 differentially regulated starch and Suc metabolism genes was assembled by combining data from the literature and a keyword search in the onion RNA-seq annotation (Table S6). In the light of the changes in Fru, Glc, and Suc levels, we focused on six DEGs which were assumed to play important roles in Suc metabolism during the early stages of onion bulb swelling.

Verification of Gene Expression by qRT-PCR

To confirm the RNA-seq results, six DEGs involved in Suc metabolism were selected for qRT-PCR analysis, comprising CWIN, SUT, SuSy (SuSy1, SuSy2), and INV (INV1, INV2). The results of the agarose gel electrophoresis indicated that all six primer pairs amplified a single band (Figure S4). The expression levels of the selected DEGs revealed by qRT-PCR were generally consistent with those from the RNA-seq analysis at the three developmental stages. Linear regression was [(qRT-PCR value) = a (RNA-Seq value) + b] and the correlation coefficient (R2) was 0.7079** (Figure 7), which indicated that the results of the RNA-Seq analysis showed a high degree of correlation with those of qRT-PCR. The qRT-PCR values and FPKM values for the six genes are shown in Table S7.

FIGURE 7
www.frontiersin.org

Figure 7. Coefficient analysis of gene expression levels determined from RNA-seq and qRT-PCR data. The qRT-PCR log2 values (x-axis) were plotted against onion bulb developmental stages (y-axis). **Indicates a significant difference at p ≤ 0.01.

The expression levels of the SuSy2 and SUT genes were highest at 15 DAS, and then decreased gradually with time, which was consistent with the decline in Suc content after 15 DAS, during which period the Glc and Fru contents increased. This would be consistent with cleavage of Suc by INV and CWIN, in accordance with the relatively higher expression of the corresponding genes at 30 DAS (Figures 6, 8). After 30 DAS, the Suc content sharply decreased, which was consistent with an increase in the expression of INV1, whereas the expression of CWIN and SuSy1 declined (Figure 8).

FIGURE 8
www.frontiersin.org

Figure 8. Relative expression levels of six DEGs involved in Suc metabolism at three different bulb developmental stages. The expression levels on the y-axis are relative to β-actin protein as control. Data are mean values of three biological replicates and error bars represent standard deviations.

Discussion

The Transcriptome and Bulb Development

RNA-seq, a powerful and efficient approach for gene expression analysis at the genome level, is especially useful in the case of non-model plants for which genomic sequence information is lacking. In this study, a total of 5416 DEGs were determined for the three different developmental stages of onion bulbs. Among them, more genes were found to be differentially expressed between 15 and 30 DAS than between the other pairs of developmental stages, with 2120 up-regulated genes and 1334 down-regulated genes. Correspondingly, the genes (CWIN, SuSy1, and INV1) associated with bulbs rapidly expanding were highly expressed at 30 DAS. These temporal gene expression patterns showed that the 30 DAS stage was the most active and rapidly expanding stage of bulb development; however, it was not clear why more DEGs were detected at 30 DAS.

The transcript distributions in the GO categories present a clear compositional depiction of the categories of genes involved in bulb formation and development. Among them, the principal category “biological process” accounts for 52,392 unigenes (42.20%) and the principal category “cellular component” accounts for 40.74%. Within “biological process,” 22.69% belonged to the category “metabolic process,” consistent with the involvement of genes of Suc metabolism in bulb development. According to the COG classifications, the category “carbohydrate transport and metabolism” was prominent, indicating that the morphogenesis and growth of bulbs requires the participation of carbohydrates. Furthermore, in the KEGG pathways, “starch and sucrose metabolism” (147, 2.40%) constituted the primary metabolism pathway in the integrated library.

Suc Metabolism in Various Plants

Besides its involvement in bulb development, Suc metabolism is also vital to fruit development. In the flavedo of Citrus paradise, there are ~22 significant enzymes involved in the synthesis of the major sugars during fruit ripening (Patel et al., 2014). A physiological analysis revealed that Solanum lycopersicum inhibitor mediated Suc metabolism by regulating CWIN activity, leading to an increase in fruit sugar content (Zhang et al., 2015). During pineapple (Ananas comosus) fruit development, the large rise in Suc was accompanied by dramatic up-regulated changes in SPS, and the cycle of Suc breakdown in the cytosol of sink tissues could be mediated by SuSy and INV (Zhang et al., 2012). The development of other fruits, such as grapes (Xin et al., 2013; Jiang et al., 2014) and peach (Amygdalus persica L.) (Desnoues et al., 2013) involves Suc metabolism. Moreover, Suc metabolism is prominent during flower blooming in Rosa chinensis “Pallida” (Yan et al., 2014).

Suc metabolism is crucial in tuber vegetables (which, though not bulbs, are examples of subterranean organ formation). Potato tubers (Solanum tuberosum L.) regulate Suc and starch metabolism in response to short-term water deficit, and there is a cycle of Suc degradation and re-synthesis in tuber discs of potato (Geigenberger, 1993; Geigenberger et al., 1997; Ferreira et al., 2010). Carrots (Daucus carota L.) expressing antisense mRNA for SuSy show markedly reduced SuSy activity in taproots (Tang and Sturm, 1999), and similar observations have been made in sweet potato (Ipomoea batatas L) (Firon et al., 2013) and turnip (Brassica rapa L.) (Gupta et al., 2001). In the current study, RNA-seq revealed 147 differentially regulated sucrose metabolism genes (Table S6). In the KEGG pathway analysis, “starch and sucrose metabolism” constituted the primary metabolism pathway (Table 1). All of the studies in the literature support the view that Suc metabolism is essential for bulb development in bulbous ornamentals (Zheng et al., 2012).

Gene Expression during Bulb Development

SUTs play a central role, as they orchestrate sucrose allocation both intracellularly and at the whole plant level. Suc produced in mesophyll cells of leaves may be effluxed into the apoplasm of mesophyll or phloem parenchyma cells by SUTs or by sucrose/H+ antiport (Kühn and Grof, 2010). The expression of SUT genes has been found to be very high in young tissues, consistent with our present studies (Figure 8). CWINs are typically expressed in sinks in which phloem unloading or subsequent postphloem transport follows an apoplasmic pathway. For example, in young tomato fruit, CWIN is abundant in the placenta phloem parenchyma, where Suc is unloaded apoplasmically (Jin et al., 2009). Interestingly, there is little or only weak CWIN expression in the phloem of the tomato fruit pericarp, which exhibits symplasmic unloading early in development (Ruan and Patrick, 1995; Jin et al., 2009). These studies showed that apoplasmic and symplasmic pathways could operate simultaneously in different cellular sites within an organ, matched by high and low CWIN levels, respectively. In the present study, the expression level of CWIN was highest at 30 DAS (Figure 8), which was consistent with the subsequent decrease in the content of Suc and with the fact that the bulbs were undergoing rapid expansion.

SuSy can reversibly catalyze Suc breakdown, whereas INV catalyzes Suc breakdown irreversibly (Li et al., 2014). Nevertheless (despite its name), SuSy is generally deemed to be involved in Suc hydrolysis rather than in Suc synthesis (Angeles-Núñez and Tiessen, 2010). In our study, the expression levels of SuSy2 and SUT genes were relatively higher at 15 DAS compared to 30 DAS, or 40 DAS, whereas CWIN, INV1, and INV2 genes were expressed at relatively low levels (Figure 8). It is possible that the Suc content increased gradually before 15 DAS; in this case, Suc would be mainly involved in metabolism prior to 15 DAS, providing energy for bulb swelling and development, and also playing a crucial role given that Suc serves as an important signaling molecule in relation to cellular metabolic status (Smeekens and Hellmann, 2014). Suc metabolism controls the generation of the sugar signaling molecules Suc, Glc, and T6P, and then regulates plant developmental processes, including meristem activation, by modulating the expression levels of transcription factors (Ruan, 2014). Consistently, the expression of some SuSy isoforms has been found to be induced or enhanced under hypoxia or anoxia in Arabidopsis (Bieniawska et al., 2007), in maize roots (Chourey et al., 1998), and in potato tubers (Fu and Park, 1995). The presence of the simple disaccharide, T, which is the major disaccharide of various bacteria, algae, fungi, and insects and which generates two Glc molecules on hydrolysis (Patrick et al., 2013), might explain why the content of Glc was still higher than that of Fru (Figure 6). VINs have long been considered to play a major role in cell expansion through their osmotic effect, and their activity is high in many rapidly expanding tissues and in tissues that are accumulating sugar, such as elongating Arabidopsis roots, expanding tomato fruits, and carrot taproots (Klann et al., 1993; Tang et al., 1999; Sergeeva et al., 2006; Ruan et al., 2010). The pattern of gene expression in the current work indicated that Suc content was positively correlated with the expression of SuSy and INV (Schäfer et al., 2004; Uys et al., 2007), which would be consistent with bulb swelling and the development of a sweet taste on account of increasing Glc and Fru content.

Conclusions

In this study, we used three different developmental stages of onion (A. cepa L.) bulb, without access to a reference genome, and conducted a comparative gene expression analysis at the transcriptional level using RNA-seq. A total of 5416 DEGs were determined across the three different stages; of these, according to the KEGG database, 147 DEGs participated in starch, and Suc metabolism. Taking into account the content of Suc, Glc, and Fru at the different developmental stages, six critical DEGs involved in Suc metabolism, comprising CWIN, SUT, SuSy (SuSy1, SuSy2), and INV (INV1, INV2), were selected. The expression levels of SUT and SuSy2 were highest at the early stage (15 DAS), and were lowest at 30 DAS; in contrast, the expression levels of CWIN and INV (INV1, INV2) increased. The findings were consistent with a decrease in the content of Suc and increases in the levels of Fru and Glc prior to 30 DAS. The expression levels of CWIN, SuSy1, and INV2 were all maximal at 30 DAS, after which the content of Suc decreased, consistent with the onset of the rapid expansion period. After 30 DAS, the contents of Fru and Glc and the expression level of INV1 sharply increased and reached maximum values at 40 DAS. These results suggest that programmed gene expression, together with changes in the contents of Suc, Fru, and Glc, are vital to bulb swelling. The findings will be of value in improving our understanding of the molecular mechanisms involved in bulb swelling and development in Allium spp. and closely related plants.

Author Contributions

CZ analyzed the data, performed the experiments, and drafted the manuscript. HZ participated in the data analysis and helped to draft the manuscript. ZZ helped to analyze the data. BL and ZC contributed analysis tools and helped to draft the manuscript. YL conceived the study, participated in its coordination, and helped to draft the manuscript. All of the authors have read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was funded by grants from the National Key Technology Research and Development Program of China (2012BAD01B00, 2012BAD50G01, and 2014BAD01B08). Our research was also supported by the Technological Innovation Capacity Program of the Beijing Academy of Agricultural and Forestry Sciences (KJCX20150111).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01425

References

Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106

PubMed Abstract | CrossRef Full Text | Google Scholar

Angeles-Núñez, J. G., and Tiessen, A. (2010). Arabidopsis sucrose synthase 2 and 3 modulate metabolic homeostasis and direct carbon towards starch synthesis in developing seeds. Planta 232, 701–718. doi: 10.1007/s00425-010-1207-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Apweiler, R., Bairoch, A., Wu, C. H., Barker, W. C., Boeckman, B., Ferro, S., et al. (2004). UniProt: the universal protein knowledgebase. Nucleic Acids Res. 32, D115–D119. doi: 10.1093/nar/gkh131

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. Nat. Genet. 25, 25–29. doi: 10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Bieniawska, Z., Paul Barratt, D. H., Garlick, A. P., Thole, V., and Kruger, N. J. (2007). Analysis of the sucrose synthase gene family in Arabidopsis. Plant J. 49, 810–828. doi: 10.1111/j.1365-313X.2006.03011.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Blencowe, B. J., Ahmad, S., and Lee, L. J. (2009). Current-generation high-throughput sequencing: deepening insights into mammalian transcriptomes. Genes Dev. 23, 1379–1386. doi: 10.1101/gad.1788009

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabib, E., and Leloir, L. F. (1958). The biosynthesis of trehalose phosphate. J. Biol. Chem. 231, 259–275.

PubMed Abstract | Google Scholar

Chourey, P. S., Taliercio, E. W., Carlson, S. J., and Ruan, Y.-L. (1998). Genetic evidence that the two isozymes of sucrose synthase present in developing maize endosperm are critical, one for cell wall integrity and the other for starch biosynthesis. Mol. Gen. Genet. 259, 88–96. doi: 10.1007/s004380050792

PubMed Abstract | CrossRef Full Text | Google Scholar

Conesa, A., Götz, S., Garcia-Gomez, J. M., Terol, J., Talón, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610

PubMed Abstract | CrossRef Full Text | Google Scholar

Darbyshire, B., and Henry, R. J. (1979). The association of fructans with high percentage dry weight in onion cultivars suitable for dehydrating. J. Sci. Food Agric. 30, 1035–1038. doi: 10.1002/jsfa.2740301103

CrossRef Full Text | Google Scholar

Darbyshire, B., and Henry, R. J. (1981). Differences in fructan content and synthesis in some Allium species. New Phytol. 87, 249–256. doi: 10.1111/j.1469-8137.1981.tb03196.x

CrossRef Full Text | Google Scholar

Deng, Y., Li, J. Q., Wu, S. F., Zhu, Y. P., Chen, Y. W., and He, F. C. (2006). Integrated nr database in protein annotation system and its localization. Comput. Eng. 32, 71–72.

Google Scholar

Desnoues, E., Gibon, Y., Baldazzi, V., Signoret, V., Génard, M., and Quilot-Turion, B. (2013). Profiling sugar metabolism during fruit development in a peach progeny with different fructose-to-glucose ratios. BMC Plant Biol. 14:336. doi: 10.1186/s12870-014-0336-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira, S. J., Senning, M., Sonnewald, S., Keβling, P. M., Goldstein, R., and Sonnewald, U. (2010). Comparative transcriptome analysis coupled to X-ray CT reveals sucrose supply and growth velocity as major determinants of potato tuber starch biosynthesis. BMC Genomics 11:93. doi: 10.1186/1471-2164-11-93

PubMed Abstract | CrossRef Full Text | Google Scholar

Firon, N., LaBonte, D., Villordon, A., Kfir, Y., Solis, J., Lapis, E., et al. (2013). Transcriptional profiling of sweetpotato (Ipomoea batatas) roots indicates down-regulation of lignin biosynthesis and up-regulation of starch biosynthesis at an early stage of storage root formation. BMC Genomics 14:460. doi: 10.1186/1471-2164-14-460

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, H., and Park, W. D. (1995). Sink- and vascular-associated sucrose synthase functions are encoded by different gene classes in potato. Plant Cell 7, 1369–1385. doi: 10.1105/tpc.7.9.1369

PubMed Abstract | CrossRef Full Text | Google Scholar

Geigenberger, P., Reimholz, R., Geiger, M., Merlo, L., Canale, V., and Stitt, M. (1997). Regulation of sucrose and starch metabolism in potato tubers in response to short-term water deficit. Planta 201, 502–518. doi: 10.1007/s004250050095

CrossRef Full Text | Google Scholar

Geigenberger, P. S. M. (1993). Sucrose synthase catalyses a readily reversible reaction in vivo in developing potato tubers and other plant tissues. Planta 189, 329–339. doi: 10.1007/BF00194429

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from rna-seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Grzelak, K., Milala, J., Król, B., Adamicki, F., and Badełek, E. (2009). Content of quercetin glycosides and fructooligosaccharides in onion stored in a cold room. Eur. Food Res. Technol. 228, 1001–1007. doi: 10.1007/s00217-009-1018-z

CrossRef Full Text | Google Scholar

Gupta, A. K., Singh, J., and Kaur, N. (2001). Sink development, sucrose metabolising enzymes and carbohydrate status in turnip (Brassica rapa L.). Acta Physiol. Plant. 23, 31–36. doi: 10.1007/s11738-001-0019-8

CrossRef Full Text | Google Scholar

Havey, M. J., Galmarini, C. R., Gökçe, A. F., and Henson, C. (2004). QTL affecting soluble carbohydrate concentrations in stored onion bulbs and their association with flavor and health-enhancing attributes. Genome 47, 463–468. doi: 10.1139/g04-005

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Shi, L., Yuan, Z., Xu, Z., Zhang, Z., and Yi, M. (2008). Effects of lipoxygenase on the corm formation and enlargement in Gladiolus hybridus. Sci. Hortic. 118, 60–69. doi: 10.1016/j.scienta.2008.05.019

CrossRef Full Text | Google Scholar

Jakše, M., and Bohanec, B. (2003). “Doubled haploid production in crop plants,” in Haploid Induction in Onion via Gynogenesis, eds M. Maluszynski, K. J. Kasha, B. P. Forster, and I. Szarejko (New York, NY: Springer Netherlands), 281–285.

Google Scholar

Jiang, H. Y., Li, W., He, B. J., Gao, Y. H., and Lu, J. X. (2014). Sucrose metabolism in grape (Vitis vinifera L) branches under low temperature during overwintering covered with soil. Plant Growth Regul. 72, 229–238. doi: 10.1007/s10725-013-9854-z

CrossRef Full Text | Google Scholar

Jin, Y., Ni, D. A., and Ruan, Y.-L. (2009). Posttranslational elevation of cell wall invertase activity by silencing its inhibitor in tomato delays leaf senescence and increases seed weight and fruit hexose level. Plant Cell 21, 2072–2089. doi: 10.1105/tpc.108.063719

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Goto, S., Kawashima, S., Okuno, Y., and Hattori, M. (2004). The KEGG resource for deciphering the genome. Nucleic Acids Res. 32, D277–D280. doi: 10.1093/nar/gkh063

PubMed Abstract | CrossRef Full Text | Google Scholar

Klann, E. M., Chetelat, R. T., and Bennett, A. B. (1993). Expression of acid invertase gene controls sugar composition in tomato (Lycopersicon) fruit. Plant Physiol. 103, 863–870.

PubMed Abstract | Google Scholar

Koch, K. (2004). Sucrose metabolism: regulatory mechanisms and pivotal roles in sugar sensing and plant development. Curr. Opin. Plant Biol. 7, 235–246. doi: 10.1016/j.pbi.2004.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kühn, C., and Grof, C. P. (2010). Sucrose transporters of higher plants. Curr. Opin. Plant Biol. 13, 287–297. doi: 10.1016/j.pbi.2010.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., Trapnell, C., and Pop, M. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25. doi: 10.1186/gb-2009-10-3-r25

PubMed Abstract | CrossRef Full Text

Li, B., and Colin, N. D. (2011). RSEM: accurate transcript quantification from RNA Seq data with or without a reference genome. BMC Bioinform. 12:323. doi: 10.1186/1471-2105-12-323

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X. Y., Wang, C. X., Cheng, J. Y., Jing, Z., Silva, J. A., Liu, X. Y., et al. (2014). Transcriptome analysis of carbohydrate metabolism during bulblet formation and development in Lilium davidii var. unicolor. BMC Plant Biol. 14:358. doi: 10.1186/s12870-014-0358-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative pcr and the 2 -ΔΔCT, method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Mallor, C., Balcells, M., Mallor, F., and Sales, E. (2011). Genetic variation for bulb size, soluble solids content and pungency in the Spanish sweet onion variety Fuentes de Ebro. Response to selection for low pungency. Plant Breed. 130, 55–59. doi: 10.1111/j.1439-0523.2009.01737.x

CrossRef Full Text | Google Scholar

McCallum, J., Clarke, A., Pither-Joyce, M., Shaw, M., Butler, R., Brash, D., et al. (2006). Genetic mapping of a major gene affecting onion bulb fructan content. Theor. Appl. Genet. 112, 958–967. doi: 10.1007/s00122-005-0199-5

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Hara, L. E., Paul, M. J., and Wingler, A. (2013). How do sugars regulate plant growth and development? New insight into the role of trehalose-6-phosphate. Mol. Plant 6, 261–274. doi: 10.1093/mp/sss120

PubMed Abstract | CrossRef Full Text | Google Scholar

Patel, M., Manvar, T., Apurwa, S., Ghosh, A., Tiwari, T., and Chikara, S. K. (2014). Comparative de novo transcriptome analysis and metabolic pathway studies of Citrus paradisi flavedo from naive stage to ripened stage. Mol. Biol. Rep. 41, 3071–3080. doi: 10.1007/s11033-014-3166-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Patrick, J. W., Botha, F. C., and Birch, R. G. (2013). Metabolic engineering of sugars and simple sugar derivatives in plants. Plant Biotechnol. J. 11, 142–156. doi: 10.1111/pbi.12002

PubMed Abstract | CrossRef Full Text | Google Scholar

Paul, M. J., Primavesi, L. F., Jhurreea, D., and Zhang, Y. (2008). Trehalose metabolism and signaling. Annu. Rev. Plant Biol. 59, 417–441. doi: 10.1146/annurev.arplant.59.032607.092945

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertea, G., Huang, X. Q., and Liang, F. (2003). TIGR gene indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics 19, 651–652. doi: 10.1093/bioinformatics/btg034

PubMed Abstract | CrossRef Full Text | Google Scholar

Riesmeier, J. W., Willmitzer, L., and Frommer, W. B. (1994). Evidence for an essential role of the sucrose transporter in phloem loading and assimilate partitioning. EMBO J. 13, 1–7.

PubMed Abstract | Google Scholar

Ruan, Y.-L. (2014). Sucrose metabolism: gateway to diverse carbon use and sugar signaling. Annu. Rev. Plant Biol. 65, 33–67. doi: 10.1146/annurev-arplant-050213-040251

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruan, Y.-L., Jin, Y., Yang, Y. J., Li, G. J., and Boyer, J. S. (2010). Sugar input, metabolism, and signaling mediated by invertase: roles in development, yield potential, and response to drought and heat. Mol. Plant 3, 942–955. doi: 10.1093/mp/ssq044

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruan, Y.-L., and Patrick, J. W. (1995). The cellular pathway of postphloem sugar transport in developing tomato fruit. Planta 196, 434–444. doi: 10.1007/BF00203641

CrossRef Full Text | Google Scholar

Schäfer, W. E., Rohwer, J. M., and Botha, F. C. (2004). A kinetic study of sugarcane sucrose synthase. Eur. J. Biochem. 271, 3971–3977. doi: 10.1111/j.1432-1033.2004.04288.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sergeeva, L. I., Keurentjes, J. J., Bentsink, L., Vonk, J., van der Plas, L. H., Koornneef, M., et al. (2006). Vacuolar invertase regulates elongation of Arabidopsis thaliana roots as revealed by QTL and mutant analysis. Proc. Natl. Acad. Sci. U.S.A. 103, 2994–2999. doi: 10.1073/pnas.0511015103

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinclair, P. J., Blakeney, A. B., and Barlow, E. W. R. (1995). Relationships between bulb dry matter content, soluble solids concentration and non-structural carbohydrate composition in the onion (Allium cepa). J. Sci. Food Agric. 69, 203–209. doi: 10.1002/jsfa.2740690210

CrossRef Full Text | Google Scholar

Smeekens, S., and Hellmann, H. A. (2014). Sugar sensing and signaling in plants. Front. Plant Sci. 5:113. doi: 10.3389/fpls.2014.00113

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, G. Q., Lüscher, M., and Sturm, A. (1999). Antisense repression and vacuolar and cell wall invertase in transgenic carrot alters early plant development and sucrose partitioning. Plant Cell 11, 177–189. doi: 10.1105/tpc.11.2.177

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, G. Q., and Sturm, A. (1999). Antisense repression of sucrose synthase in carrot (Daucus carota L.) affects growth rather than sucrose partitioning. Plant Mol. Biol. 41, 465–479. doi: 10.1023/A:1006327606696

PubMed Abstract | CrossRef Full Text | Google Scholar

Tatusov, R. L., Galperin, M. Y., Natale, D. A., and Koonin, E. V. (2000). The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 28, 33–36. doi: 10.1093/nar/28.1.33

PubMed Abstract | CrossRef Full Text | Google Scholar

Tewari, G. M., and Bandyopadhyay, C. (1975). Quantitative evaluation of lachrymatory factor in onion by thin-layer chromatography. J. Agric. Food Chem. 23, 645–647. doi: 10.1021/jf60200a044

CrossRef Full Text | Google Scholar

Trapnell, C., Williams, B. A., Pertea, G., and Mortazavi, A. (2010). Transcript assembly and quantification by RNA Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

Tymowska-Lalanne, Z., and Kreis, M. (1998). The plant invertases: physiology, biochemistry, and molecular biology. Adv. Bot. Res. 28, 71–117. doi: 10.1016/S0065-2296(08)60294-3

CrossRef Full Text | Google Scholar

Uys, L., Botha, F. C., Hofmeyr, J. H. S., and Rohwer, J. M. (2007). Kinetic model of sucrose accumulation in maturing sugarcane culm tissue. Phytochemistry 68, 2375–2392. doi: 10.1016/j.phytochem.2007.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Vandeputte, G. E., and Delcour, J. A. (2004). From sucrose to starch granule to starch physical behaviour: a focus on rice starch. Carbohyd. Polym. 58, 245–266. doi: 10.1016/j.carbpol.2004.06.003

CrossRef Full Text | Google Scholar

Wang, Z., Gerstein, M., and Snyder, M. (2009). RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet. 10, 57–63. doi: 10.1038/nrg2484

PubMed Abstract | CrossRef Full Text | Google Scholar

Winter, H., and Huber, S. C. (2000). Regulation of sucrose metabolism in higher plants: localization and regulation of activity of key enzymes. Crit. Rev. Biochem. Mol. Biol. 35, 253–289. doi: 10.1080/10409230008984165

PubMed Abstract | CrossRef Full Text | Google Scholar

Xin, H., Zhang, J., Zhu, W., Wang, N., Fang, P., Han, Y., et al. (2013). The effects of artificial selection on sugar metabolism and transporter genes in grape. Tree Genet. Genomes 9, 1343–1349. doi: 10.1007/s11295-013-0643-7

CrossRef Full Text | Google Scholar

Yan, H. J., Zhang, H., Chen, M., Jian, H. Y., Baudino, S., Caissard, J. C., et al. (2014). Transcriptome and gene expression analysis during flower blooming in Rosa chinensis ‘Pallida’. Gene 540, 96–103. doi: 10.1016/j.gene.2014.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Z., Chen, L. C., Suzuki, H., Ariyada, O., Erra-Balsells, R., Nonami, H., et al. (2009). Direct profiling of phytochemicals in tulip tissues and in vivo monitoring of the change of carbohydrate content in tulip bulbs by probe electrospray ionization mass spectrometry. J. Am. Soc. Spectrom. 20, 2304–2311. doi: 10.1016/j.jasms.2009.08.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, N., Jiang, J., Yang, Y. L., and Wang, Z. H. (2015). Functional characterization of an invertase inhibitor gene involved in sucrose metabolism in tomato fruit. J. Zhejiang Univ. Sci. B 16, 845–856. doi: 10.1631/jzus.B1400319

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X. M., Wang, W., Du, L. Q., Xie, J. H., Yao, Y. L., and Sun, G. M. (2012). Expression patterns, activities and carbohydrate-metabolizing regulation of sucrose phosphate synthase, sucrose synthase and neutral invertase in pineapple fruit during development and ripening. Int. J. Mol. Sci. 13, 9460–9477. doi: 10.3390/ijms13089460

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, R. R., Wu, Y., and Xia, Y. P. (2012). Chlorocholine chloride and paclobutrazol treatments promote carbohydrate accumulation in bulbs of Lilium Oriental hybrids ‘Sorbonne’. J. Zhejiang Univ. Sci. B 13, 136–144. doi: 10.1631/jzus.B1000425

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Allium cepa L., bulb swelling, RNA-seq, sucrose metabolism, gene expression

Citation: Zhang C, Zhang H, Zhan Z, Liu B, Chen Z and Liang Y (2016) Transcriptome Analysis of Sucrose Metabolism during Bulb Swelling and Development in Onion (Allium cepa L.). Front. Plant Sci. 7:1425. doi: 10.3389/fpls.2016.01425

Received: 02 June 2016; Accepted: 07 September 2016;
Published: 22 September 2016.

Edited by:

Rudy Dolferus, Commonwealth Scientific and Industrial Research Organisation, Australia

Reviewed by:

Trine (Anne Kathrine) Hvoslef-Eide, Norwegian University of Life Sciences, Norway
Xiao-Li Tan, Jiangsu University, China

Copyright © 2016 Zhang, Zhang, Zhan, Liu, Chen and Liang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yi Liang, liangyi@nercv.org

These authors have contributed equally to this work