- Research article
- Open Access
Transcriptome changes during fruit development and ripening of sweet orange (Citrus sinensis)
BMC Genomics volume 13, Article number: 10 (2012)
The transcriptome of the fruit pulp of the sweet orange variety Anliu (WT) and that of its red fleshed mutant Hong Anliu (MT) were compared to understand the dynamics and differential expression of genes expressed during fruit development and ripening.
The transcriptomes of WT and MT were sampled at four developmental stages using an Illumina sequencing platform. A total of 19,440 and 18,829 genes were detected in MT and WT, respectively. Hierarchical clustering analysis revealed 24 expression patterns for the set of all genes detected, of which 20 were in common between MT and WT. Over 89% of the genes showed differential expression during fruit development and ripening in the WT. Functional categorization of the differentially expressed genes revealed that cell wall biosynthesis, carbohydrate and citric acid metabolism, carotenoid metabolism, and the response to stress were the most differentially regulated processes occurring during fruit development and ripening.
A description of the transcriptomic changes occurring during fruit development and ripening was obtained in sweet orange, along with a dynamic view of the gene expression differences between the wild type and a red fleshed mutant.
The typical course of fruit development involves expansion, sweetening and increasing pigmentation . From the consumers' point of view, the appearance, texture and taste of the fruit are all of high importance. These properties involve attaining a suitable composition of sugars, organic acids, amino acids and carotenoids. The underlying mechanisms of fruit development and ripening have been extensively studied in tomato , but are not well explored in non-climacteric fruits. Citrus is a widely grown fruit crops, which exhibits non-climacteric ripening behaviour. Its fruit contains a juicy pulp made of vesicles within segments . The growth and development of the citrus fruit can be divided into three stages: cell division, an expansion phase involving cell enlargement and water accumulation, and the ripening stage . In the latter stage, carotenoids and other soluble solids are accumulated, chlorophyll is lost, the cell wall is extensively modified, the organic acid content is reduced, and the concentration of a number of volatiles increases.
Citrus fruit provides a convenient vehicle to study gene regulation during non-climacteric fruit development and ripening. In grape (Vitis vinifera), another non-climacteric fruit, several comprehensive mRNA expression profiling studies have been presented to describe fruit development and ripening [5–8], while in citrus, an EST sequencing project showed that 20% of the sequences were metallothionein . Application of a citrus cDNA microarray has suggested that unique genetic regulatory networks arise during fruit development , while a more global transcriptome analysis was able to identify ethylene-responsive genes in the mandarin fruit . More recently, a comprehensive study of the clementine fruit transcriptome has proposed a mechanism for citrate utilization . However, these transcriptomic studies of fruit development in citrus have mainly been based on microarray analysis, the next generation sequencing technology provides new opportunities for more accurate and powerful deep transcriptome analysis of fruit development.
The ripening of citrus fruit is accompanied by carbohydrate build-up, acid reduction, and carotenoid accumulation. Citrus fruits accumulate most of their sucrose in the juice cells . Glucose is decomposed via glycolysis and the pentose phosphate pathway (OPP). As the major acid present in citrus fruit, citric acid contributes > 90% of the total organic acid content. Citric acid synthesis is thought to take place in the mitochondria via the tricarboxylic acid (TCA) cycle. Although the genes encoding some of the key enzymes (specially, citrate synthase (CS), aconitase , and NADP-isocitrate dehydrogenase (NADP-IDH) ) have been isolated, their activity cannot completely account for variation in the level of citric acid in citrus fruit . Carotenoids are also important components of the citrus fruit, and their composition and content in sweet orange fruits have been extensively studied [17, 18]. Cross-talk between carotenoid, sugar and organic acid metabolism has been documented [19, 20]. The presence of sucrose may promote colour break , while its deficiency delays lycopene accumulation in tomato . The down-regulation of CS and NAD-dependent IDH results in a decrease in the levels of both organic acids and carotenoids .
The fruit of the bud mutant Hong Anliu is characterized by a high sucrose and low citric acid level and its tendency to accumulate lycopene is responsible for the red pigmentation of its pulp . The mutant appears to be highly isogenic with its progenitor wild type variety Anliu, as revealed by a genotypic analysis based on microsatellites and AFLPs. The mutant has therefore provided a platform to study cross-talk between primary (sugar and organic acid synthesis) and secondary (carotenoid synthesis) metabolism enzymes, which are fundamental in the determination of citrus fruit quality. Suppressive subtraction hybridization combined with cDNA microarray analysis has been applied to determine that the differentially expressed genes were mainly enriched in the stage of 170 DAF (days after flowering) in the mutant fruit . At this stage, a total of 582 genes were found to be differentially expressed between the wild type Anliu (hereafter WT) and Hong Anliu (hereafter MT) as revealed by RNA-seq analysis . However, how genes are dynamically and differentially expressed during fruit development and ripening has not yet been determined. Here, the developmental changes of fruit transcriptome of sweet orange were investigated.
Plant material and RNA preparation
WT and MT plants were both cultivated in the same orchard at the Institute of Citrus Research (Guilin, Guangxi Province, China), with the same climatic conditions. Fruit samples were harvested at 120, 150, 190, and 220 DAF (days after flowering) from three different trees in 2009. At each developmental stage, ten representative fruits were sampled from each tree. The pulp was separated from the peel, and the pulp was sliced. The sliced WT pulps samples were combined with one another (similarly for the MT samples), snap-frozen in liquid nitrogen and kept at -80°C until required [24, 25]. One aliquot was used to extract RNA isolation, as described previously . The remainder of the powder was used for the determination of sugar and organic acid composition and concentration, and the content of H2O2.
RNA-seq and functional assignment
The WT and MT fruit pulp harvested at 120, 150, 190, and 220 DAF was subjected to RNA-seq using an Illumina Genome Analyzer at Beijing Genomics Institute (Shenzhen) in 2009. The abundance of each tag was normalized to one transcript per million (TPM) for between sample comparison purposes. The raw data was filtered to remove low quality sequences including ambiguous nucleotides, adaptor sequences, and below 3 TPM, as described previously . The sequencing data can be accessed at the website: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=dxqjxoygumyauzm&acc=GSE22505. To link the expressed signatures to known genes from orange, the TIGR unigene dataset (http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=orange) was used as a reference database. The Z-score method using the p-value as a statistical significance index  was applied to identify differentially expressed genes. A cluster analysis was performed according to Eisen et al. ; the log2 of TPM for each gene was used for the hierarchical clustering analysis. Gene Ontology (GO) categorization was carried out as described previously . The ultra-geometric test was applied to perform GO enrichment analysis. In the significance analysis of the enrichment of a GO item, the p-value represents the probability of satisfying the hypothesis that the designated genes involved in the GO item has not been enriched (statistical significance at P = 0.05).
Real-time quantitative RT-PCR
The differential expression of a selection of the genes identified as being differentially expressed was validated by applying real-time quantitative RT-PCR (qRT-PCR). The sequences of the primer pairs (designed using Primer Express 3.0 (Applied Biosystems, Foster City, CA, USA)) are listed in additional file 1. All qRT-PCRs were performed using an ABI 7500 Real Time System (Applied Biosystems) using the actin gene as the reference . Primers for both the target gene and the reference were diluted in SYBR GREEN PCR Master Mix (Applied Biosystems) and 20 μL of the reaction mix were added to each well. Reactions were performed via an initial incubation at 50°C for 2 min and at 95°C for 10 min, and then cycled at 95°C for 15 s, and 60°C for 60 s for 40 cycles. The resulting data were handled by the instrument on-board software Sequence Detector Version 1.3.1 (Applied Biosystems).
Analysis of sugar, organic acid and H2O2
Soluble sugar and organic acid composition and concentrations were determined by gas chromatography (GC) using 3 g of the powdered pulp as described previously  with minor modifications. The powder was suspended in chilled 80% methanol and then held in a 75°C water bath for 30 min. After a 2 h ultrasonic extraction and centrifugation at 4000g for 10 min, the supernatant was collected and 1 mL internal standard (2.5% w/v phenyl-β-D-glucopyranoside, 2.5% w/v methyl-α-D-glucopyranoside) was added. The solution was made up to 50 mL with 80% methanol, and a 2 mL aliquot was centrifuged at 12000g for 15 min. A 0.5 mL aliquot of this final supernatant was vacuum-dried and then re-dissolved in 800 μL 2% w/v hydroxylamine hydrochloride in pyridine at 75°C for 1 h. Then 400 μL hexamethyldisilazane and 200 μL trimethylchlorosilane were added and the sample was held at 75°C for 2 h. A 0.5 μL aliquot was used for GC analysis in an Agilent 6890N device (Santa Clara, CA, USA) equipped with a flame ionization detector. A capillary column (HP-5, 5% phenyl-methylpolysiloxane, 30 m × 25 μm i.d. × 0.1 μm) was employed, with nitrogen as the carrier gas at a flow rate of 45 mL/min, and flow-rates of hydrogen and air set to 40 mL/min and 450 mL/min, respectively. Sugars and organic acids were identified through a comparison of retention times using standard compounds from Sigma (St. Louis, MO, USA). The concentration of H2O2 was measured using a hydrogen peroxide detection kit supplied by Nanjing Jiancheng Institute of Biological Technology (Nanjing, China). A 0.8 g sample of powdered pulp was suspended in 7.2 ml saline (0.90% w/v of NaCl) and centrifuged for 10 min at 10, 000g. The intensity of yellow complex formed by the reaction of molybdate and H2O2, as measured spectrophotometrically at 405 nm, was used to assess the concentration of H2O2. Three replicates were conducted for each sample.
The fruit transcriptome sampled at four developmental stages
In total, eight cDNA preparations were sequenced from fruit pulp sampled at 120, 150, 190, and 220 DAF from WT and MT. The average number of tags produced for each library was 4.01 million (Table 1). The raw data were submitted and available from the NCBI/GEO repository (accession number GEO: GSE22505; website: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=dxqjxoygumyauzm&acc=GSE22505). After filtering, the number of robust tags per library ranged from 2.5 to 4.3 million, and the number of distinct tags from 61,000 to 113,000. A saturation analysis (additional file 2) demonstrated that as sequencing depth was increased, the number of new distinct tags decreased, but only until the number of sequences had reached 2.5 million. We concluded therefore that the libraries were all fully saturated and hence large enough for gene expression analysis. The distribution of distinct tag abundance and total tag number exhibited very similar tendencies for all eight libraries (additional file 3). Transcripts which accounted for nearly 60% of the total number were in less than 7% of the categories, and transcripts that accounted for 40% of the categories were less than 5% of the total number, indicating that only a few genes were expressed at a high level.
Transcriptome changes during fruit development and ripening
To map tags to known genes, a reference citrus unigene dataset containing 26,826 contigs and 73,607 singletons was used. The procedure identified between 68.1% and 76.2% of the tags (additional file 4), of which 20,155 to 36,173 (31.8% to 33.1%) produced unambiguous identifications (one tag mapping to one gene). The libraries were relatively uniform with respect to mapping efficiency. A total of 18,829 genes were detected in at least one of the four stages in the wild type sweet orange, of which 8,825 genes were expressed in all the four stages. In this study, we solely used the wild type sweet orange as a model to demonstrate the transcriptome changes during fruit development and ripening. Three genes were most highly expressed in wild type, two of which were encoded a stress-response protein (one 22 kDa polypeptide response to low temperature stress) and a heat shock protein, while the function of the third one (F28C11.8) is unknown. Changes in the transcriptome during fruit development and ripening were examined by cluster analysis of gene expression patterns, which arranged the 18,829 genes into 22 groups (Figure 1); the 10,005 (53.1%) genes expressed in three or less of the three stages fell into groups 1 to 11. The largest group (20) comprised the 3,075 (16.3%) genes whose expression increased continuously during fruit development and ripening; this group included the genes encoding sucrose phosphate synthase, cysteine proteinase and a sucrose transporter. The second largest group (10) contained the 2,970 (15.8%) genes which were not expressed at 120 DAF but maintained a stable expression level at other three developmental stages. The 2,618 (13.9%) genes in group 2 were not expressed at 120 and 190 DAF. The cluster analysis also revealed that the abundance of 89.7% of the transcripts detected in the WT pulp varied over the course of fruit development and ripening (Figure 1). Many of the transcripts were single stage-specific (additional file 5). A comparison of expression patterns between WT (Figure 1) and MT (additional file 6) revealed that 20 of the groups were common to both, while 97.5% of the genes expressed in MT showed a similar expression pattern to that in WT (the exceptions belonged to groups 13 and 19).
Differentially expressed genes during fruit development and ripening
Of the 18,829 genes detected in the WT, 9,377, 7,886, and 7,757 were differentially expressed between 120 DAF and 150 DAF, 150 DAF and 190 DAF, and 190 DAF and 220 DAF, respectively. Of these, 36.7% were assigned to one of 18 GO categories (Figure 2). The categories "metabolic process", "cellular process", "establishment of localization", "localization", "biological regulation", "pigmentation", and "response to stimulus" based on biological process captured most of these genes. Several of the differentially expressed genes were associated with cell wall metabolism and softening (Table 2). A gene encoding pectinesterase (TC14614) showed a decrease in expression during fruit development and ripening, consistent with the modification of pectin which occurs during ripening. A xyloglucan endotransglycosylase gene (TC9277) was up-regulated during fruit development and ripening, while another (TC25537) was only detectable at the early stage of fruit development. The abundance of three expansin (EY741042, TC20787, and EY738078) and three alpha-expansin (EY727139, EY743186, and EY699153) transcripts increased during fruit development and ripening. Some genes involved in sucrose metabolism, the TCA cycle and carotenoid biosynthesis were found to be differentially expressed (Table 2). The expression of a sucrose phosphate synthase gene family member (EY677217) increased during ripening, and its expression profile correlated well with the rise in total soluble solids. Related genes encoded CS (EY703799), aconitase-iron regulated protein 1 (CB293814), phytoene synthase (PSY) (EY722043) and phytoene desaturase (PDS) (TC3). Some of the differentially expressed genes were associated with stress response (Table 2): these included genes encoding a heat shock protein and various metallothioneins. Most of the metallothionein genes were up-regulated in the later stages of ripening. Three cysteine protease genes (EY718215, TC22740, and EY707198), were also highly expressed at the late stage of fruit development and ripening. Also, the expression of a valencene synthase gene (TC363) increased during fruit development and ripening.
Developmental difference between fruit transcriptome of MT and WT
The comparison between the transcriptomes of WT and MT identified 6,540 genes showing a two-fold or greater difference in expression level at 120 DAF, 3,529 at 150 DAF, 4,601 at 190 DAF, and 3,289 at 220 DAF. Application of the Z-score method  suggested that the expression of 634, 568, 540, and 616 of these was significantly different at p < 0.05 and |log2Ratio| ≥ 1 in the four developmental stages, respectively (additional file 7). The top ten differentially expressed genes in each of the four developmental stages have been listed in additional file 8. Many encode stress-related products, such as cysteine protease Cp5 and metallothionein-like protein. At all the four developmental stages the number of up-regulated genes was less than that of down-regulated genes (Figure 3).
When the expression profiles of the genes differentially expressed between MT and WT were subjected to a cluster analysis (additional file 9), over one half (492/883) turned out to be up-regulated in MT at all the developmental stages except 150 DAF. Some examples of this large cluster included genes encoding Cu/Zn superoxide dismutase (TC12069), ascorbate peroxidase (EY685405, TC22775 and TC14669), MYB1 (EY649968) and plastidic glucose 6-phosphate/phosphate translocator (EY722703). Only five genes were detected as differentially expressed at all four stages; one encoded a cysteine protease (Cp5) (TC5370), which had been reported exhibiting the acid-activatable cysteine protease forms , three had no assigned functions, and one shared no homology with any entry in. The five genes could be classified into two groups (two genes in one group, three genes in the other), with their expression patterns of being opposite to each other but with both groups having a turning point at 150 DAF (additional file 10).
The GO categories of the set of differentially expressed genes (additional file 11) revealed that most encoded products associated with "protein binding", "hydrolase activity", "transferase activity" and "transporter activity". At the 150 DAF stage, the most common categories were "signal transducer activity" and "electron carrier activity". GO enrichment analysis was also carried out to identify which key processes were altered in MT (Table 3). The genes encoding phytoene synthase (EY722043) and ζ-carotene desaturase (ZDS) (TC15628) involved in carotenoid metabolic process were enriched at 150 DAF in MT compared to WT, and those encoding capsanthin/capsorubin synthase (CCS) (TC5834 and TC9367) were also enriched at 150 DAF in MT.
Verification of differentially expressed genes during fruit development and ripening
RNA sampled from fruit of MT and WT at the four selected stages of fruit development and ripening provided the template for qRT-PCR validation of the sequence-based transcription profiles for 22 of the differentially expressed genes (Figure 4). Linear regression [(RNA-seq value) = a(RT-PCR value) + b] analysis showed an overall correlation coefficient of 0.8379, indicating a good correlation between transcript abundance assayed by real-time PCR and the transcription profile revealed by RNA-seq data (Figure 4E).
Changes in fruit soluble sugars, organic acids, carotenoids and H2O2 content
Since the expression of a number of genes implicated in carbohydrate metabolism and mitochondria-related citric acid metabolism varied greatly during fruit development and ripening, an attempt was made to monitor the dynamics of pulp soluble sugar and organic acid content. The content of soluble sugars increased markedly during the late stages of fruit development and ripening in both WT and MT (Figure 5A). The concentrations of glucose and fructose remained rather constant up to 150 DAF, and thereafter tended to be higher in MT than in WT. The concentration of sucrose was higher in MT than in WT throughout fruit development and ripening. The citric acid content fluctuated, but was lower in MT than in WT at all developmental stages (Figure 5B). The concentration of quinic acid decreased substantially over the course of fruit development and ripening, but did not differ between WT and MT. Carotenoids and lycopene both accumulated over time in MT, but remained at a low level in WT (Figure 5C). Finally, H2O2 content fell as the fruit developed and ripened, but was higher in MT than in WT at 120 DAF (Figure 5D).
Transcriptome dynamics during fruit development and ripening
Fruit ripening is a highly coordinated, genetically programmed and irreversible process which involves a series of physiological, biochemical, and organoleptic changes allowing for the development of an edible ripe fruit . Fruit ripening in citrus is accompanied by carbohydrate build-up, acid reduction, carotenoid accumulation and chlorophyll degradation . The bud mutant MT produces fruit with high sucrose and lycopene, but low citric acid content . Here, transcriptome changes over the course of fruit development and ripening in MT and WT were monitored and annotated. We must point out that we collected samples from three different trees and pooled for RNA-seq analysis for each developmental stage. We did not include biological replicates considering that we used a pair of genotypes (WT and MT) and the data from WT and MT could corroborate to each other to some extent. We sequenced each pool once technically since the next generation sequencing data are highly replicable with relatively little technical variation . We also used real-time qRT-PCR to verify the transcription profile revealed by RNA-seq data (Figure 4), with an overall correlation coefficient of 0.8379 which indicated the RNA-seq data was reliable. The RNA-seq approach detected a similar number of genes (18,829 genes in WT and 19,440 genes in MT) in both genotypes, with no apparent difference between the global number of genes expressed in the two genotypes at any of the four stages of fruit development and ripening sampled (additional file 4). The WT and MT patterns of gene expression were also largely alike (20 of the 24 expression patterns were similar). Thus, we believed that overall picture of the transcriptome captured by RNA-seq is robust.
An overwhelming proportion of the genes identified (84.8% in MT and 89.7% in WT) varied in their level of expression over the course of fruit development and ripening, reflecting the occurrence of a massive genetic re-programming. A large number of these genes were expressed in a stage-specific manner, which implicates their involvement in physiological processes which take place only at a specific developmental stage(s). A major group of the differentially expressed genes was involved in cell wall modification, which is not surprising since the major textural changes associated with the softening of fruit are due to enzyme-mediated alterations in the structure and composition of the cell wall , especially the cell walls of juice sacs in citrus. Changes in the activity of several cell wall-related genes were known to result in the abnormal development of juice sac granulation [35, 36], while modifications in cell wall structure or in the components of the membranes of the segments and juice sacs during fruit development and ripening clearly influenced the formation of the fruit pulp melting characteristic . One of the cell wall-related genes revealed here encoded a pectinesterase, an enzyme which modifies the assembly and disassembly of pectin, a common component of the primary cell wall. In tomato, the gene for pectinesterase was highly expressed prior to ripening, and was down-regulated by ethylene as ripening begins . Here, the expression of the gene encoding the pectinesterase PPE8B precursor decreased as the WT fruit matured (Figure 4A), while in the blood orange, the expression of a pectinesterase gene has been shown to increase during fruit development and ripening . It might be due to the different members identified in the two studies. The hemicellulose xyloglucan is a common component of the cell wall, and is hydrolysed and transglycosylated by xyloglucan endotransglycosylase in growing tissues and ripening fruits . Here, a gene encoding this enzyme was up-regulated during fruit development and ripening (Figure 4A), indicating its probable role in cell wall degradation during ripening. Both the early enlargement of the citrus fruit driven by cell expansion and the later ripening process require the presence of expansins to loosen the cell walls , and several genes encoding expansins were detected in the present study and their expression was higher at the later stages of fruit development and ripening.
The accumulation of carbohydrates represents one of the most obvious changes which occur during citrus fruit development and ripening. The perceived changes in expression of genes involved in carbohydrate metabolism here were consistent with the findings of other transcriptomic analyses in citrus . The type of sugar deposited to a high level in the cell vacuole in citrus is predominantly sucrose, unlike in grape, where it is glucose and fructose . Reflecting this difference, the expression profile of the gene encoding sucrose synthase in citrus  was rather different from that in the grape . Nevertheless, in both species, sugar is important for the regulation of colour development, and perhaps also for other ripening processes . It was notable that the gene encoding a sucrose phosphate synthase was induced during fruit development and ripening, perhaps because the activity of this enzyme may be required for the re-synthesis of sucrose to allow its further transport to the vacuole . Of some interest also was the behaviour of a gene encoding valencene synthase, an enzyme the activity of which was known to be induced by ethylene, and which was part of the natural ripening process in citrus . Valencene is an important component of the aroma of the ripe sweet orange fruit. Both in the present experiments (Figure 4A) and elsewhere , valencene synthase transcript accumulated in the ripening fruit.
Molecular processes involved in the formation of the fruit traits of red flesh sweet orange
Suppressive subtraction hybridization, in combination with cDNA microarray analysis, has identified a set of 267 genes which were differentially expressed between MT and WT . RNA-seq technology is more capable of identifying a nearly complete inventory of transcripts and by this method a total of 582 genes were found to be differentially expressed between WT and MT at the stage of 170 DAF . In the current study, we tested on four developmental stages of WT and MT and extended our understanding of the global and dynamic changes during fruit development and ripening in MT and WT. Almost all of the members of the 267 gene set revealed by SSH strategy were also identified in the present study, although some discrepant expression patterns were apparent between this new data set and previously assembled microarray-based set. For example, the abundance of the transcript encoding cysteine protease, which appeared to differ between MT and WT in both studies, was documented by the microarray analysis as being lower at all developmental stages, whereas it appeared to be higher at 150 DAF in MT in the present study (additional file 10).
The major biological processes occurring in the mitochondria (TCA cycle, coupling electron transfer, and oxidative phosphorylation) were remarkably altered in MT. The intermediates of the TCA cycle can be channelled into the syntheses of fats, terpenoids, porphyrins, nucleotides, and amino acids. In MT, the level of citric acid, the major organic acid present in citrus fruit, was consistently around 25% that present in WT fruit (Figure 5). However, no major difference was detected in the expression of the genes encoding CS and aconitase-iron regulated protein (Figure 4A), two predominant enzymes involved in the TCA cycle. Five differentially expressed genes, all associated with mitochondria-related processes (coupling electron transfer and oxidative phosphorylation) were down-regulated in MT compared with WT (Figure 4B). These included the genes encoding a NADH-ubiquinone oxidoreductase 75kDa subunit (TC7962), cytochrome C oxidase (TC18256), and cytochrome oxidase subunit 3 (TC156), suggesting that MT mitochrondria were capable of less efficient electron transport than that WT ones. If, as a result, flux through the TCA cycle is decreased, the accumulation of citric acid is likely to be compromised. In addition, PDS and ZDS catalyzing desaturation of phytoene to lycopene involve net electron transfer . In tomato, a NAD(P)H dehydrogenase complex which participating in electron transfer was involved in carotenoid biosynthetic pathway , suggesting the possible exist of cross-talk between electron transfer and carotenoid accumulation in sweet orange.
In the plant cell, the mitochondrial electron transport chain is a major site of reactive oxygen species (ROS) production . Here, the concentration of one of the primary ROS molecules (H2O2) was higher in MT than that in WT at 120 DAF (Figure 5D). The delicate balance between antioxidant defence and ROS production can be disrupted by either compromised antioxidant defence or the inhibition of electron flow . Here, the primary anti-oxidant enzymes (SOD, APX and GR) were more active in MT than in WT pulp, suggesting that the level of oxidative stress may be greater in MT than in WT . The expression of a large number of stress-related genes was also substantially different in MT and WT (additional file 8). Functioning as an important ROS scavenging pathway, the ascorbate-glutathione cycle requires the cofactor NADPH, which is provided by the OPP pathway . This pathway is also a major source of NADPH for many biosynthetic processes, including carotenoid biosynthesis , in non-photosynthetic organs such as the fruit. The gene encoding glucose 6-phosphate dehydrogenase, which is considered as the first and rate-limiting enzyme of the OPP pathway in all cells , was up-regulated in MT. A statistical analysis of the qRT-PCR result confirmed that the level of transcription of this gene was significantly higher in MT than in WT, especially at 120 DAF (Figure 4C). Carotenoid biosynthesis, another effective anti-oxidative process , was also higher in MT than in WT. Lycopene is the most potent antioxidant among the carotenoids . The expression level of several carotenoid biosynthesis genes, encoding namely PSY, ZDS, lycopene β-cyclase (LCYb) and CCS, was greatly changed in MT. PSY is generally accepted to be a rate-limiting enzyme in carotenoid biosynthesis pathway. The CCS product is an enzyme which is mechanistically similar to LCYb, and the low transcript level of CCS may well be responsible for the accumulation of lycopene in red grapefruits (Citrus paradisi) . The qRT-PCR analysis confirmed that both upstream genes (PSY and ZDS) were up-regulated and both downstream ones (LCYb and CCS) down-regulated in MT (Figure 4D), consistent with the mechanism regulating lycopene accumulation in tomato .
The present study has provided a dynamic view of the transcriptome during fruit development and ripening of a sweet orange red-fleshed mutant and its progenitor wild type. Cell wall biosynthesis, carbohydrate metabolism, the TCA cycle, and carotenoid biosynthesis were all differentially regulated during fruit development and ripening. These differentially regulated processes may well be important for the formation of the pleiotropic fruit trait of Hong Anliu sweet orange.
Coombe BG: The Development of Fleshy Fruits. Annu Rev Plant Physiol. 1976, 27: 207-228. 10.1146/annurev.pp.27.060176.001231.
Fei ZJ, Tang XM, Alba RM, White JA, Ronning CM, Martin GB, Tanksley SD, Giovannoni JJ: Comprehensive EST analysis of tomato and comparative genomics of fruit ripening. Plant J. 2004, 40 (1): 47-59. 10.1111/j.1365-313X.2004.02188.x.
Iglesias DJ, Cercós M, Colmenero-Flores JM, Naranjo MA, Ríos G, Carrera E, Ruiz-Rivero O, Lliso I, Morillon R, Tadeo FR, Talon M: Physiology of citrus fruiting. Braz J Plant Physiol. 2007, 19: 333-362.
Bain JM: Morphological, anatomical, and physiological changes in the developing fruit of the Valencia orange, Citrus sinensis (L) Osbeck. Aust J Bot. 1958, 6 (1): 1-23. 10.1071/BT9580001.
Mohammad A, Shangwu C, Yechun W, Oliver Y, Laszlo K, Wenping Q: Berry skin development in Norton grape: Distinct patterns of transcriptional regulation and flavonoid biosynthesis. BMC Plant Biol. 2011, 11: 7-10.1186/1471-2229-11-7.
Deluc LG, Grimplet J, Wheatley MD, Tillett RL, Quilici DR, Osborne C, Schooley DA, Schlauch KA, Cushman JC, Cramer GR: Transcriptomic and metabolite analyses of Cabernet Sauvignon grape berry development. BMC genomics. 2007, 8 (1): 429-10.1186/1471-2164-8-429.
Pilati S, Perazzolli M, Malossini A, Cestaro A, Demattè L, Fontana P, Dal Ri A, Viola R, Velasco R, Moser C: Genome-wide transcriptional analysis of grapevine berry ripening reveals a set of genes similarly modulated during three seasons and the occurrence of an oxidative burst at veraison. BMC genomics. 2007, 8 (1): 428-10.1186/1471-2164-8-428.
Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiol. 2010, 152 (4): 1787-10.1104/pp.109.149716.
Moriguchi T, Kita M, Hisada S, Endo-Inagaki T, Omura M: Characterization of gene repertoires at mature stage of citrus fruits through random sequencing and analysis of redundant metallothionein-like genes expressed during fruit development1. Gene. 1998, 211 (2): 221-227. 10.1016/S0378-1119(98)00138-3.
Shimada T, Fuiii H, Endo T, Yazaki J, Kishimoto N, Shimbo K, Kikuchi S, Omura M: Toward comprehensive expression profiling by microarray analysis in citrus: monitoring the expression profiles of 2213 genes during fruit development. Plant Sci. 2005, 168 (5): 1383-1385. 10.1016/j.plantsci.2005.01.013.
Fujii H, Shimada T, Sugiyama A, Nishikawa F, Endo T, Nakano M, Ikoma Y, Shimizu T, Omura M: Profiling ethylene-responsive genes in mature mandarin fruit using a citrus 22K oligoarray. Plant Sci. 2007, 173 (3): 340-348. 10.1016/j.plantsci.2007.06.006.
Cercós M, Soler G, Iglesias DJ, Gadea J, Forment J, Talón M: Global analysis of gene expression during development and ripening of citrus fruit flesh. A proposed mechanism for citric acid utilization. Plant Mol Biol. 2006, 62 (4): 513-527. 10.1007/s11103-006-9037-7.
Soule J, Grierson W: Anatomy and physiology. Fresh citrus fruits. Edited by: W. Wardowski, S. Nagy, and W. Grierson. 1986, AVI, Westport, Conn, 1-22.
Sadka A, Dahan E, Cohen L, Marsh KB: Aconitase activity and expression during the development of lemon fruit. Physiol Plantarum. 2000, 108 (3): 255-262. 10.1034/j.1399-3054.2000.108003255.x.
Sadka A, Dahan E, Or E, Cohen L: NADP(+)-isocitrate dehydrogenase gene expression and isozyme activity during citrus fruit development. Plant Sci. 2000, 158 (1-2): 173-181. 10.1016/S0168-9452(00)00328-9.
Katz E, Fon M, Lee YJ, Phinney BS, Sadka A, Blumwald E: The citrus fruit proteome: insights into citrus fruit metabolism. Planta. 2007, 226 (4): 989-1005. 10.1007/s00425-007-0545-8.
Wang YC, Chuang YC, Hsu HW: The flavonoid, carotenoid and pectin content in peels of citrus cultivated in Taiwan. Food Chem. 2008, 106 (1): 277-284. 10.1016/j.foodchem.2007.05.086.
Fanciullino AL, Dhuique-Mayer C, Froelicher Y, Talón M, Ollitrault P, Morillon R: Changes in carotenoid content and biosynthetic gene expression in juice sacs of four orange varieties (Citrus sinensis) differing in flesh fruit color. J Agr Food Chem. 2008, 56 (10): 3628-3638. 10.1021/jf0732051.
Alquézar B, Zacarías L, Rodrigo MJ: Molecular and functional characterization of a novel chromoplast-specific lycopene β-cyclase from Citrus and its relation to lycopene accumulation. J Exp Bot. 2009, 60 (6): 1783-1797. 10.1093/jxb/erp048.
Liu Q, Xu J, Liu YZ, Zhao XL, Deng XX, Guo LL, Gu JQ: A novel bud mutation that confers abnormal patterns of lycopene accumulation in sweet orange fruit (Citrus sinensis L. Osbeck). J Exp Bot. 2007, 58 (15-16): 4161-4171. 10.1093/jxb/erm273.
Huff A: Sugar Regulation of Plastid Interconversions in Epicarp of Citrus Fruit. Plant Physiol. 1984, 76: 307-312. 10.1104/pp.76.2.307.
Télef N, Stammitti-Bert L, Mortain-Bertrand A, Maucourt M, Carde JP, Rolin D, Gallusci P: Sucrose deficiency delays lycopene accumulation in tomato fruit pericarp discs. Plant Mol Biol. 2006, 62: 453-469. 10.1007/s11103-006-9033-y.
Sienkiewicz-Porzucek A, Nunes-Nesi A, Sulpice R, Lisec J, Centeno DC, Carillo P, Leisse A, Urbanczyk-Wochniak E, Fernie AR: Mild reductions in mitochondrial citrate synthase activity result in a compromised nitrate assimilation and reduced leaf pigmentation but have no effect on photosynthetic performance or growth. Plant Physiol. 2008, 147 (1): 115-127. 10.1104/pp.108.117978.
Liu Q, Zhu AD, Chai LJ, Zhou WJ, Yu KQ, Ding J, Xu J, Deng X: Transcriptome analysis of a spontaneous mutant in sweet orange [Citrus sinensis (L.) Osbeck] during fruit development. J Exp Bot. 2009, 60 (3): 801-813. 10.1093/jxb/ern329.
Xu Q, Yu KQ, Zhu AD, Ye JL, Liu Q, Zhang JC, Deng XX: Comparative transcripts profiling reveals new insight into molecular processes regulating lycopene accumulation in a sweet orange (Citrus sinensis) red-flesh mutant. BMC Genomics. 2009, 10: 540-10.1186/1471-2164-10-540.
Liu YZ, Liu Q, Tao NG, Deng XX: Efficient isolation of RNA from fruit peel and pulp of ripening navel orange (Citrus sinensis Osbeck). J HAU. 2006, 25: 300-304.
Kal AJ, van Zonneveld AJ, Benes V, van den Berg M, Koerkamp MG, Albermann K, Strack N, Ruijter JM, Richter A, Dujon B: Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources. Mol Biol Cell. 1999, 10 (6): 1859-1872.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.
Bartolozzi F, Bertazza G, Bassi D, Cristoferi G: Simultaneous determination of soluble sugars and organic acids as their trimethylsilyl derivatives in apricot fruits by gas-liquid chromatography. J Chromatogr A. 1997, 758 (1): 99-107. 10.1016/S0021-9673(96)00709-1.
North MJ, Nicol K, Sands TW, Cotter DA: Acid-activatable cysteine proteinases in the cellular slime mold Dictyostelium discoideum. J Biol Chem. 1996, 271 (24): 14462-10.1074/jbc.271.24.14462.
Prasanna V, Prabha TN, Tharanathan RN: Fruit Ripening Phenomena-An Overview. Crit Rev Food Sci Nutr. 2007, 47 (1): 1-19. 10.1080/10408390600976841.
Talon M, Gmitter FG: Citrus genomics. Int J Plant Genomics. 2008, 2008: 528361-
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome research. 2008, 18 (9): 1509-1517. 10.1101/gr.079558.108.
Goulao LF, Oliveira CM: Cell wall modifications during fruit ripening: when a fruit is not the fruit. Trends Food Sci Tech. 2008, 19 (1): 4-25. 10.1016/j.tifs.2007.07.002.
Singh R: 65-year research on citrus granulation. Ind J Hort. 2001, 58 (1-2): 112-144.
Sharma RR, Saxena SK: Rootstocks influence granulation in Kinnow mandarin (Citrus nobilis C. deliciosa). Sci Hort. 2004, 101 (3): 235-242. 10.1016/j.scienta.2003.10.010.
Waldron KW, Parker ML, Smith AC: Plant cell walls and food quality. Compr Rev Food Sci Food saf. 2003, 2 (4): 128-146. 10.1111/j.1541-4337.2003.tb00019.x.
Giovannoni JJ: Genetic regulation of fruit development and ripening. Plant Cell. 2004, 16: S170-10.1105/tpc.019158.
Licciardello C, Russo MP, Vale' G, Recupero RG: Identification of differentially expressed genes in the flesh of blood and common oranges. Tree Genet Genomes. 2008, 4 (2): 315-331. 10.1007/s11295-007-0111-3.
Catala C, Rose JKC, York WS, Albersheim P, Darvill AG, Bennett AB: Characterization of a tomato xyloglucan endotransglycosylase gene that is down-regulated by auxin in etiolated hypocotyls. Plant Physiol. 2001, 127 (3): 1180-10.1104/pp.010481.
Brummell DA, Harpster MH, Dunsmuir P: Differential expression of expansin gene family members during growth and ripening of tomato fruit. Plant Mol Biol. 1999, 39 (1): 161-169. 10.1023/A:1006130018931.
Castellarin SD, Gambetta GA, Wada H, Shackel KA, Matthews MA: Fruit ripening in Vitis vinifera: spatiotemporal relationships among turgor, sugar accumulation, and anthocyanin biosynthesis. J Exp Bot. 2011
Komatsu A, Moriguchi T, Koyama K, Omura M, Akihama T: Analysis of sucrose synthase genes in citrus suggests different roles and phylogenetic relationships. J Exp Bot. 2002, 53 (366): 61-10.1093/jexbot/53.366.61.
Sharon-Asa L, Shalit M, Frydman A, Bar E, Holland D, Or E, Lavi U, Lewinsohn E, Eyal Y: Citrus fruit flavor and aroma biosynthesis: isolation, functional characterization, and developmental regulation of Cstps1, a key gene in the production of the sesquiterpene aroma compound valencene. Plant J. 2003, 36 (5): 664-674. 10.1046/j.1365-313X.2003.01910.x.
Mayer MP, Beyer P, Kleinig H: Quinone compounds are able to replace molecular oxygen as terminal electron acceptor in phytoene desaturation in chromoplasts of Narcissus pseudonarcissus L. Eur J Biochem. 1990, 191 (2): 359-363. 10.1111/j.1432-1033.1990.tb19130.x.
Nashilevitz S, Melamed-Bessudo C, Izkovich Y, Rogachev I, Osorio S, Itkin M, Adato A, Pankratov I, Hirschberg J, Fernie AR: An Orange Ripening Mutant Links Plastid NAD (P) H Dehydrogenase Complex Activity to Central and Specialized Metabolism during Tomato Fruit Maturation. Plant Cell. 2010, 22 (6): 1977-1997. 10.1105/tpc.110.074716.
Moller IM: Plant mitochondria and oxidative stress: Electron transport, NADPH turnover, and metabolism of reactive oxygen species. Annu Rev Plant Physiol Plant Mol Biol. 2001, 52: 561-591. 10.1146/annurev.arplant.52.1.561.
Turrens JF: Mitochondrial formation of reactive oxygen species. J Physiol London. 2003, 552 (2): 335-344. 10.1113/jphysiol.2003.049478.
Pan ZY, Liu Q, Yun Z, Guan R, Zeng WF, Xu Q, Deng XX: Comparative proteomics of a lycopene-accumulating mutant reveals the important role of oxidative stress on carotenogenesis in sweet orange (Citrus sinensis [L.] osbeck). Proteomics. 2009, 9 (24): 5455-5470. 10.1002/pmic.200900092.
May MJ, Vernous T, Leaver C, Van Montagu M, Inzé D: Glutathione homeostasis in plants: implications for environmental sensing and plant development. J Exp Bot. 1998, 49: 649-667. 10.1093/jexbot/49.321.649.
Neuhaus HE, Emes MJ: Nonphotosynthetic metabolism in plastids. Annu Rev Plant Physiol Plant Mol Biol. 2000, 51: 111-140. 10.1146/annurev.arplant.51.1.111.
Stryer L: Biochemistry. 1989, W.H.Freeman and Co., San Francisco, CA, 3
Couee I, Sulmon C, Gouesbet G, El Amrani A: Involvement of soluble sugars in reactive oxygen species balance and responses to oxidative stress in plants. J Exp Bot. 2006, 57 (3): 449-459. 10.1093/jxb/erj027.
Clinton SK: Lycopene: Chemistry, Biology, and Implications for Human Health and Disease. Nutr Rev. 1998, 56 (2): 35-51.
Giuliano G, Bartley GE, Scolnik PA: Regulation of carotenoid biosynthesis during tomato development. Plant Cell. 1993, 5 (4): 379-
This work was supported by the National Basic research program of China (973 program; No. 2011CB100601), and the National Natural Science Foundation of China (No. 30830078 and 31071779).
KQY, QX and XLD were responsible for generating the RNA-seq data and for the interpretation of the data. KQY carried out qRT-PCR experiments and measured sugar, organic acid and H2O2 content, and drafted the manuscript. QX conceived the study, participated in its design and helped to draft the manuscript. FG and YDD participated in the statistical analyses. XXD proposed and supervised the research. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: The primer sequence information. This file listed the primers sequences used for real-time quantitative RT-PCR validation of RNA-seq data. (DOC 53 KB)
Additional file 2: The saturation evaluations of the eight libraries in this study. This file contained the information of the saturation evaluations of the RNA-seq tags in the eight libraries (MT and WT at four selected fruit developmental stages) against sequencing depth. The results revealed that with the increase of total sequence number (sequencing depth), the number of genes identified increased, but the number stabilized once the number of sequences reached 2.5 million, indicating enough information has been included in the RNA-seq data. (DOC 115 KB)
Additional file 3: Distribution of total tags number (A) and distinct tags number (B) in MT and WT at 120, 150, 190, and 220 DAF. This file showed the distribution of the number of total tags and distinct tags obtained in MT and WT at different developmental stages. (DOC 66 KB)
Additional file 4: Summary of tags mapped against a reference set of sweet orange unigenes. This file contained the summary result of tags mapping against a reference set of sweet orange unigenes. (DOC 42 KB)
Additional file 5: Number of stage-specific genes expressed in MT and WT. This file contained the summary result of stage-specific genes number in MT and WT. (DOC 38 KB)
Additional file 6: . This file contained the result of the hierarchical cluster analysis of genes expression profiles in MT. The log2 of transcripts per million (TPM) for each gene was used for the hierarchical clustering analysis at four developmental stages (120, 150, 190 and 220 DAF). In all, 19,440 genes were classified into 22 regulatory patterns, designated groups 1-22. (DOC 1 MB)
Additional file 7: List of differentially expressed genes between MT and WT. The table contained information of the differentially expressed genes with expression difference > 2, and genes differentially expressed at 0.05 significance level at each of the four fruit developmental stages. (XLS 700 KB)
Additional file 8: The ten most differentially expressed genes between MT and WT at each of the four selected fruit developmental stages. This file listed the ten most differentially expressed genes between MT and WT at different developmental stages, with their expression ratios between MT and WT, also containing simple annotation information. (DOC 106 KB)
Dynamics patterns of gene expression of a set of genes differentially expressed between MT and WT at each of the four selected fruit developmental stages
Additional file 9: . This file contained the result of the hierarchical cluster analysis of expression profiles of differentially expressed genes between MT and WT at different developmental stages. The log2 of the ratio between the MT and the WT TPM for each gene was used to perform the cluster analysis. (DOC 56 KB)
Additional file 10: . This file contained the pattern of genes which were differentially expressed at all selected stages. At each stage (120, 150, 190 and 220 DAF), the log2 of the ratio between the MT and the WT TPM for each gene is represented. (DOC 49 KB)
Additional file 11: Functional categorization of genes differentially expressed between WT and MT. This file showed the distribution of GO categories of differentially expressed genes between WT and MT at the four selected stages of fruit development and ripening. The categorization was based on molecular activity of Gene Ontology items. Percentages are based on the proportion of the number of genes in each set. (DOC 211 KB)
About this article
Cite this article
Yu, K., Xu, Q., Da, X. et al. Transcriptome changes during fruit development and ripening of sweet orange (Citrus sinensis). BMC Genomics 13, 10 (2012). https://doi.org/10.1186/1471-2164-13-10
- Fruit Development
- Citrus Fruit
- Sweet Orange