Skip to main content

Comparative analysis of carotenoid accumulation in two goji (Lycium barbarum L. and L. ruthenicumMurr.) fruits



The traditional Chinese medicinal plants Lycium barbarum L. and L. ruthenicum Murr. are valued for the abundance of bioactive carotenoids and anthocyanins in their fruits, respectively. However, the cellular and molecular mechanisms contributing to their species-specific bioactive profiles remain poorly understood.


In this study, the red fruit (RF) of L. barbarum was found to accumulate high levels of carotenoids (primarily zeaxanthin), while they were undetectable in the black fruit (BF) of L. ruthenicum. Cytological and gene transcriptional analyses revealed that the chromoplast differentiation that occurs in the chloroplast during fruit ripening only occurs in RF, indicating that the lack of chromoplast biogenesis in BF leads to no sink for carotenoid storage and the failure to synthesize carotenoids. Similar enzyme activities of phytoene synthase 1 (PSY1), chromoplast-specific lycopene β-cyclase (CYC-B) and β-carotene hydroxylase 2 (CRTR-B2) were observed in both L. ruthenicum and L. barbarum, suggesting that the undetectable carotenoid levels in BF were not due to the inactivation of carotenoid biosynthetic enzymes. The transcript levels of the carotenoid biosynthetic genes, particularly PSY1, phytoene desaturase (PDS), β-carotene desaturase (ZDS), CYC-B and CRTR-B2, were greatly increased during RF ripening, indicating increased zeaxanthin biosynthesis. Additionally, carotenoid cleavage dioxygenase 4 (CCD4) was expressed at much higher levels in BF than in RF, suggesting continuous carotenoid degradation in BF.


The failure of the chromoplast development in BF causes low carotenoid biosynthesis levels and continuous carotenoid degradation, which ultimately leads to undetectable carotenoid levels in ripe BF. In contrast, the successful chromoplast biogenesis in RF furnishes the sink necessary for carotenoid storage. Based on this observation, the abundant zeaxanthin accumulation in RF is primarily determined via both the large carotenoid biosynthesis levels and the lack of carotenoid degradation, which are regulated at the transcriptional level.


Carotenoids are isoprenoids that are synthesized by all photosynthetic organisms as well as some non-photosynthetic bacteria and fungi. In plants, chloroplastic carotenoids are constituents of light-harvesting complexes and the photosynthetic reaction center, where they also play important roles in protecting tissues against photo-oxidative damage [1],[2]. When accumulated in the chromoplasts of flowers and fruits, carotenoids act as visual attractants for pollinating insects and seed-dispersing animals [3],[4]. Furthermore, carotenoids are the precursors of important apocarotenoids, such as volatile flavor/aroma terpenes, and the growth regulators abscisic acid (ABA) and strigolactone [5]-[7]. Recently, oxidized products from plant carotenoids have been implicated as signals induced by environmental stressors [8]. In addition to these biological functions, carotenoids serve as major micronutrients in the human diet [9],[10]. In particular, β-carotene, β-carotene and β-cryptoxanthin are precursors for vitamin A biosynthesis [11], while lutein and zeaxanthin slow aging-related damage to the retina [12].

During the past two decades, the carotenoid biosynthetic pathway in plants has been well elucidated [13]-[15]. Previous studies have shown that there are at least three mechanisms that regulate carotenoid accumulation in the chromoplasts [16]. First, the transcript abundance of rate-limiting structural genes is predicted to be the primary mechanism controlling the carotenoid content and composition in the chromoplasts [17]. During tomato (Solanum lycopersicum) fruit development, increasing expression of phytoene synthase 1 (PSY1) and phytoene desaturase (PDS), diminishing expression of the chloroplast-related lycopene β-cyclase (LCY-B) and lycopene β-cyclase (LCY-E), and low transcript levels of the chromoplast-specific LCY-B and β-ring hydroxylase (CRTR-B) (corresponding to CYC-B and CRTR-B2, respectively), lead to the accumulation of lycopene as the major carotenoid [18]-[21].

Second, carotenoid degradation by carotenoid cleavage dioxygenases (CCDs) may be central to determining the final carotenoid concentrations in chromoplasts [22]. For example, despite active carotenoid biosynthesis in both the yellow and white petals of chrysanthemums (Chrysanthemum morifolium), the carotenoids are degraded by CmCCD4a into colorless compounds in the white petals [23]. In potatoes and peaches, the different carotenoid content among the cultivars can also be attributed to the distinct enzymatic activity of the CCD4-degrading carotenoids [24]-[26].

Finally, the sink capacity of carotenoid-accumulating tissues has recently been implicated in the control of carotenoid levels. The characterization of the Orange (Or) protein, which is involved in chromoplast biogenesis, revealed the importance of the carotenoid storage sink for carotenoid accumulation. Due to a failure in chromoplast formation, the cauliflower (Brassica oleracea) Or mutant lacks carotenoid accumulation [27]. When the Or gene was transformed into Arabidopsis, the Arabidopsis calli exhibited an orange color with chromoplast formation [28]. In tomatoes, the perturbed activity of several light-signal-related genes, including UV-Damaged DNA Binding Protein1 (DDB1) [29], De-Etiolated1 (DET1) [30],[31], Cullin4 (CUL4) [32], HY5, COP1LIKE [33], Cryptochrome 2 (CRY2) [34], Golden 2-Like (GLK2) [35] and Arabidopsis Pseudo Response Regulator2-Like (APRR2-Like) [36], caused changes in the plastid number and size, which indirectly affected the concentrations of the carotenoids and other phytonutrients in the ripening fruits [37].

Lycium barbarum L. (Chinese: gouqi or ningxiagouqi) and L. ruthenicum Murr. (Chinese: heigouqi or heiguogouqi), which are two shrub plants belonging to the Solanaceae family, have been used as traditional medicinal plants in China and other Asian countries for centuries [38]. L. barbarum, in particular, has high economic significance in Northwest China, with its red fruit (RF; gouqizi in Chinese, also known as goji berry or wolfberry) being used for both traditional Chinese medicine (TCM) and nutritional purposes [39]. Modern pharmacological studies have begun to investigate the biochemical mechanisms of the medicinal effects of wolfberry, including the antioxidant, immunomodulatory and neuroprotective properties, which are primarily attributed to the polysaccharides (LBP), flavonoids and carotenoids [40],[41]. L. ruthenicum is another TCM used for the treatment of heart disease, abnormal menstruation and menopause [42]. The functional compounds in the black fruit (BF) of L. ruthenicum are primarily comprised of anthocyanins, essential oils and polysaccharides [42]-[45]. As the primary pigment in RF, carotenoids have been extensively studied, and zeaxanthin and esterified zeaxanthin were reported to be the major bioactive compounds that accumulate in RF, especially for its traditional use in eyesight improvement [41],[46]. However, the content and composition of the carotenoids in BF have not been comprehensively reported, and the mechanisms controlling the species differences in the carotenoid biosynthesis between RF and BF remain unknown. Analyses of these differences may provide novel insights into the regulation of carotenoid accumulation in goji fruits, with important implications for their medicinal and nutritional value.


The carotenoid accumulation differs between the RF and BF from different Lyciumspecies

The analysis of the total carotenoid content in the red fruits of L. barbarum and the black fruits of L. barbarum at four developmental stages (S1-S4, Figure 1) revealed an increase in the carotenoid content of RF from S2 to S4, reaching a maximum of 508.90 μg g-1 fresh weight (FW) (Additional file 1). On the converse, the amount of carotenoids in BF declined from 34.46 μg g-1 FW in the S1 fruit to undetectable levels in the S4 fruit (Additional file 1).

Figure 1
figure 1

Photographs of L. barbarum and L. ruthenicum fruits (RF and BF, respectively) at different developmental stages (S1-S4). S1, green fruit stage; S2, color break stage; S3, light color stage; S4, ripe fruit stage.

The carotenoid composition and content in ripe RF were previously reported [46], with zeaxanthin accounting for the highest carotenoid proportion, followed by β-cryptoxanthin and β-carotene, and with most of the xanthophylls esterified. In this study, to detect the carotenoid accumulation regardless of esterification, the carotenoid content in the four developmental stages (S1-S4, Figure 1) of RF and BF was analyzed after saponification. Xanthophyll esters were undetectable in BF (data not shown). In S1 of RF and BF, the chloroplastic carotenoids, violaxanthin, lutein and β-carotene comprised the majority of the carotenoids (Additional files 2 and 3). Two additional compounds (both unidentified) were detected in BF. As both the red and black fruits developed, the amount of chloroplastic carotenoids (lutein, violaxanthin and β-carotene) declined (Additional files 2 and 3). During BF ripening, no additional carotenoids showed rising levels, and all of the existing carotenoids gradually decreased to undetectable levels (Figure 2B; Additional files 1 and 3). Meanwhile, in RF, several other carotenoids, especially zeaxanthin, increased dramatically from S2 to S4 (Figure 2A; Additional files 1 and 2). Specifically, zeaxanthin reached 381.6 μg g-1 FW, and β-cryptoxanthin and β-carotene reached 17.59 μg g-1 FW and 28.99 μg g-1 FW, respectively (Additional file 1).

Figure 2
figure 2

Content of zeaxanthin, β-cryptoxanthin and β-carotene in RF (A) and BF (B).

Light microscopy of the green fruits (S1) and ripe fruits (S4)

The chromoplast differentiation in the fruits of the two Lycium species was comparatively determined by examining the plastids in the mesocarp of the green fruits (S1) and the ripe fruits (S4) under a light microscope. The chloroplasts were observed in the green fruits of both species (Figures 3A and C). However, orange, globular chromoplasts were only observed in the ripe RF of L. barbarum (Figure 3B). Consistent with the absence of carotenoid accumulation, a failure in chromoplast formation was observed in the ripe BF of L. ruthenicum (Figure 3D).

Figure 3
figure 3

Light micrographs of plastids in RF and BF. (A) Green RF (S1) cell with chloroplasts. (B) Ripe RF (S4) cell with orange globular chromoplasts. (C) Green BF (S1) cell with chloroplasts. (D) Ripe BF (S4) cell without colour chromoplasts. Fruits are not stained to show the natural colour of plastids.

Isolation of the carotenogenesis-related genes from L. barbarum and L. ruthenicum

To compare the gene sequences encoding the enzymes responsible for the biosynthesis, degradation, and storage of carotenoids in both species, the full-length open reading frames (ORF) of twenty-five putative carotenogenesis-related genes were isolated (Table 1). The phylogenetic relationship of each putative protein (with known functional proteins in other organisms) was used to confirm the orthologous relationships of these proteins with the clearly defined proteins in tomatoes (Additional file 4). Two isoforms of 1-deoxy-D-xylulose 5-phosphate synthase (DXS1 and DXS2), which are involved in the first step of the 2-C-methyl-Derythritol4-phosphate (MEP) pathway, were isolated in each species. In particular, three pairs of carotenoid biosynthetic genes (PSY1/PSY2, LCY-B/CYC-B, CRTR-B1/CRTR-B2) were isolated; PSY1, CYC-B and CRTR-B2 are putatively specific for carotenoid biosynthesis in chromoplasts [18]. The ORF length of the majority of the genes [except for 15-cis-β-carotene isomerase (Z-ISO), P450-type ε-ring hydroxylase (CYP97C11), CRTR-B1 and zeaxanthin epoxidase (ZEP)] was identical between the two Lycium species. The average identity of the protein sequences between L. barbarum and L. ruthenicum was 98.42%, while L. barbarum and S. lycopersicum shared 89.03% identity and L. ruthenicum and S. lycopersicum shared 88.83% identity. Consistent with other species, most of the proteins were predicted to localize to the chloroplast by ProtComp (Table 1).

Table 1 Sequence information of the carotenogenesis-related genes from L. barbarum and L. ruthenicum

The comparative RNA-seq profile of the carotenogenesis-related genes in the ripening fruits of the two Lyciumspecies

To comparatively overview the expression of the carotenogenesis-related genes in both species, RNA-seq data derived from the fruit s S1 to S3 stages were profiled in this study. Generally, the parameters, transcriptional read amounts and reads per kilobase of coding sequence per million reads (RPKM) are used for assessing the gene expression levels when analyzing RNA-seq data. As shown in Additional file 5, the RPKMs of the 25 carotenogenesis-related genes were calculated. The RPKM of the chromoplast-related genes (CHRC, Or1 and HSP21) in RF were much higher than those in BF, suggesting that these genes are more active in RF than in BF. Particularly, the RPKM of CHRC reached nearly 30,000 in RF (Figure 4). In RF, the RPKMs of some of the carotenoid biosynthetic genes (PDS, ZDS, CYC-B and CRTR-B2) showed increasing trends during fruit ripening and approached the hundreds in S2 and S3 (Figure 4, Additional file 5). In contrast, in all three BF stages, the RPKMs of all of the carotenoid biosynthetic genes were less than fifty (Additional file 5). During BF development, only the LrCCD4 transcripts increased and sharply reached 2,000 RPKM in S3 (Figure 4, Additional file 5). However, the LbCCD4 expression obviously declined from S1 to S3 (Figure 4). These results suggest that more carotenoids are degraded in BF than in RF.

Figure 4
figure 4

The RPKM values, calculating from the RNA-seq data of L. barbarum and L. ruthenicum , of eight carotenoid-ralated genes which are obviously changing during fruit development. LrCCD4, carotenoid cleavage dioxygenase 4 from L. ruthenicum. Other seven genes are from L. barbarum: PDS, phytoene desaturase; ZDS, β-carotene desaturase; CYC-B, chromoplast-specific lycopene β-cyclase; CRTR-B2, non-heme di-iron carotenoid ε-ring hydroxylase 2; CHRC, Chromoplast-specific carotenoid-associated protein; HSP21, heat shock protein 21; CCD4, carotenoid cleavage dioxygenase 4.

Comparative analysis of the carotenogenesis-related gene expression in L. barbarum and L. ruthenicumfruits via qRT-PCR

To confirm the expression patterns of all 25 of the carotenogenesis-related genes during RF and BF ripening (S1-S4), qRT-PCR was used (Figure 5). Consistent with the findings from the RNA-seq data, the LrCHRC transcripts were very low in BF, while the LbCHRC transcripts were abundant in RF, particularly in the S2 and S3 stages. Both Or genes displayed constant expression during BF ripening, while the LbOr1 transcript level was much higher (5.9- to 9.2-fold) in RF S1-S3 than in S4. The transcript abundance of HSP21 was increased by 7.6-fold from S1 to S4 of RF, while it decreased to undetectable levels in S3 and S4 of BF.

Figure 5
figure 5

Expression patterns of carotenoid-related genes in ripening RF ( L. barbarum ; gray bars) and BF ( L. ruthenicum ; black bars). The expression of Actin1 was used to normalize the mRNA levels for each sample. Three replicates were performed for each sample. LF, leaf samples; the fruit developmental stages (S1-S4) shown are identical to those depicted in Figure 1. DXS, 1-deoxy-D-xylulose 5-phosphate-synthase; PSY, phytoene synthase; PDS, phytoene desaturase; Z-ISO, 15-cis-β-carotene isomerase; ZDS, β-carotene desaturase; CRTISO, carotene isomerase; LCY-B, lycopene β-cyclase; CYC-B, chromoplast-specific lycopene β-cyclase; CRTR-B, non-heme di-iron carotenoid ε-ring hydroxylase; CYP97A29, P450 carotenoid ε-ring hydroxylase; CYP97C11, P450 carotenoid ε-ring hydroxylase; NCED, 9-cis-epoxycarotenoid dioxygenase; CCD, carotenoid cleavage dioxygenase; CHRC, Chromoplast-specific carotenoid-associated protein; Or, orange; HSP21, heat shock protein 21; ZEP, zeaxanthin epoxidase; VDE, violaxanthin de-epoxidase.

As the first enzyme in the MEP pathway, DXS was shown to be a regulatory enzyme in tomato fruit carotenogenesis [47]. In Lycium, DXS1 showed similar expression profiles in the two species examined here, with much higher mRNA levels (10-fold) in the leaves than in the fruit. Specifically, the LrDXS2 transcripts were equally abundant in the leaves during S1 and S2, before decreasing in S3 and S4, while the LbDXS2 transcripts increased during the color-break stage (S2) by 10-fold and gradually decreased thereafter (Figure 5). The transcripts of the putative chromoplast-specific genes (PSY1, CYC-B and CRTR-B2) increased dramatically (by 64-, 17- and 53-fold, respectively) during RF ripening, whereas they were consistently expressed at relatively low levels throughout BF ripening (Figure 5). Similarly, the PDS, ZDS, and CRTISO transcript abundance was low during BF ripening but increased (by 30-, 6.7- and 6.3-fold, respectively) at the color-break stage (S2) in RF and remained high thereafter (Figure 5). The Z-ISO transcript was not detected in the leaves or fruits of either species via qRT-PCR. Similarly, the lutein synthesis genes LCY-E, CYP97A29 and CYP97C11 were expressed at low levels or were undetectable in the leaves and fruits of both species (Figure 5). The ZEP and VDE genes, which act downstream of zeaxanthin, showed similar expression profiles, with moderate transcript abundance in the leaves of both L. barbarum and L. ruthenicum, where they likely participate in the xanthophyll cycle. In the fruits however, both genes were consistently expressed during the ripening of BF but were barely detected in RF (Figure 5).

The transcript abundance of NCED1 increased gradually throughout BF ripening, while in RF, it peaked at S2 and decreased thereafter (Figure 5). No NCED6 transcripts were detected in the leaves or fruits of either species, except for extremely low levels in S1 and S2. In RF, the expression profile of LbCCD1A was similar to LbNCED1, while in BF, LrCCD1A was anti-correlated with LbNCED1. LrCCD4 was highly expressed in BF, particularly in the late developmental stages, while LbCCD4 was expressed at relatively low levels in RF and gradually decreased throughout the ripening process (Figure 5).

Functional analysis of the key carotenoid biosynthesis enzymes from L. barbarum and L. ruthenicum

To verify the functionality of the key carotenoid biosynthesis enzymes identified in the two Lycium species, the bioactivities of PSY1, CYC-B and CRTR-B2 were tested in E. coli. Previous studies have demonstrated that these are rate-controlling enzymes in chromoplast-specific carotenoid biosynthesis [18],[19],[21]. The carotenoids in transformed E. coli cells were detected using high-performance liquid chromatography (HPLC) (Figure 6). The positive controls (Table 2) showed the expected absorbance spectra corresponding to phytoene, β-carotene and zeaxanthin. For the functional assays, the peaks of the carotenoids isolated from the bacteria containing the substrate synthesizing plasmids (pACCRT-E for GGPP, pACCRT-EIB for lycopene, pACCAR16ΔcrtX for β-carotene), coupled with the vectors containing the Lycium enzymes, represented the peaks for phytoene, β-carotene and β-cryptoxanthin. Beta-cryptoxanthin is an intermediate for zeaxanthin, and therefore is an indicator of insufficient CRTR-B2 hydroxylase activity E. coli (Figure 6) [48],[49]. Overall, these results suggest similar PSY1, CYC-B and CRTR-B2 bioactivities between L. barbarum and L. ruthenicum.

Figure 6
figure 6

Pigments produced in E. coli in the functional analysis for PSY1, CYC-B and CRTR-B2 of both species. In the functional analysis of PSY1 (A), pigments extracted from E. coli cells harboring pACCRT-E, the engineered plasmid producing GGPP, and pEASY-E1, the empty vector; plasmids pACCRT-E and pEASY-LbPSY1, which encodes LbPSY1; plasmids pACCRT-E and pEASY-LrPSY1, which encodes LrPSY1; and pACCRT-EB, the engineered plasmid producing phytoene (peak 1) as a positive control. The absorption spectra of phytoene are presented in the boxes with retention times of 10.0 min. In the functional analysis of CYC-B (B), the plasmids are pACCRT-EIB, the engineered plasmid producing lycopene (peak 2), and pEASY-E1; pACCRT-EIB and pEASY-LbCYC-B, which encodes LbCYC-B; pACCRT-EIB and pEASY-LrCYC-B, which encodes LrCYC-B; and pACCAR16ΔcrtX, the engineered plasmid producing β-carotene (peak 3) as a positive control. The absorption spectra of β-carotene are presented in the boxes with retention times of 9.4 min. In the functional analysis of CRTR-B2 (C), the plasmids are pACCAR16ΔcrtX (producing β-carotene, peak 3) and pEASY-E1; pACCAR16ΔcrtX and pEASY-LbCRTR-B2, which encodes LbCRTR-B2; pACCAR16ΔcrtX and pEASY-LrCRTR-B2, which encodes LrCRTR-B2; and pACCAR25ΔcrtX, the engineered plasmid producing zeaxanthin (peak 5). Beta-cryptoxanthin standard was used as the indicator of peak 4. The absorption spectra of β-cryptoxanthin are presented in the boxes with retention times of 10.5 min.

Table 2 Construct design for enzymatic assays in E. coli


The fruits of two valuable Lyciumspecies show opposite carotenoid accumulation patterns

The red pigmentation of the ripe L. barbarum fruit is due to the high accumulation of specific carotenoids [41],[46]. Unlike RF, the ripe fruit of L. ruthenicum is deep purple in color with a high petunidin content produced by the anthocyanin pathway [42]. To unravel the molecular regulatory basis for the differences in carotenoid accumulation between RF and BF, we first characterized the carotenoid compositional changes during fruit ripening in both species.

The phytochemical analysis revealed that the carotenoid accumulation increased during RF ripening (Additional file 1). As shown in Figure 2, the zeaxanthin precursor β-cryptoxanthin accumulated at consistently low levels during the RF ripening process, which was accompanied by a high level of zeaxanthin accumulation, consistent with previous studies (Additional files 1 and 2) [46]. These results suggest that the flux through the carotenoid pathway in RF is primarily directed into the β, β-carotene branch to produce zeaxanthin. At the same time, the chloroplastic carotenoids lutein and violaxanthin, present in pre-ripe S1 green RF, gradually decreased during fruit development (Additional file 2). Beta-carotene is also a chloroplastic carotenoid, and it declined from S1 to S2 in RF; however, as another intermediate of zeaxanthin biosynthesis, it increased from S2 to S4, consistent with the zeaxanthin accumulation (Additional file 1).

In BF, the products of both the β, β-carotene and β, β-carotene branches of the carotenoid pathway, namely lutein and violaxanthin, gradually decreased to undetectable levels during ripening (Additional files 1 and 3). In contrast to the results for RF, the carotenoid compositions did not change during BF ripening, and the content of all of the existing compositions gradually declined to undetectable levels (Additional file 3). It is interesting to reveal the mechanisms underlying the different carotenoid accumulation patterns between RF and BF.

The failure in chromoplast development results in no carotenoid accumulation in ripe BF

In fruits and flowers, a large abundance of carotenoids can be stored in the chromoplasts. Following the research on the Or gene, the formation of the chromoplast was recognized as a vital factor for carotenoid accumulation [50]. In the Or cauliflower mutant, the failure in chromoplast formation blocked the biosynthesis and accumulation of carotenoids [27]. In transgenic Or-overexpressing Arabidopsis and rice, the chromoplast differentiation occurred in the calli of both species and induced the biosynthesis of the carotenoids [28],[51]. In addition, CHRC and HSP21 were also shown to play significant roles in chromoplast development and carotenoid storage [52],[53].

In this study, we observed that orange, globular chromoplasts existed in the cells of the ripe RF (Figure 3B), consistent with the abundant carotenoid accumulation. Likewise, consistent with the poor accumulation of carotenoids in BF, these organelles were not observed in the ripe BF (Figure 3D). Therefore, the development of the chromoplasts may be the primary cause of the differences in carotenoid accumulation between RF and BF. The expression profiles of the chromoplast-related genes (Or, CHRC and HSP21) during RF and BF development also supported this speculation. Or1, CHRC and HSP21 all showed much higher expression levels in RF compared to BF (Figures 4 and 5).

Within the chromoplast, the carotenoids and the CHRC protein are predominantly stored in lipoprotein fibrils [54]. Distinct from the plastoglobules in the chloroplasts, these fibrils are characterized by a high homogeneity of apolar compounds, most of which are esterified xanthophylls [55]. In potato tubers, a positive correlation between the total carotenoid content and the esterified xanthophyll fraction was observed, suggesting that esterification facilitates the accumulation of these lipophilic compounds within the plastids [56]. Recently, this viewpoint was also verified in apples (Malus x domestica Borkh) [57]. The majority of the zeaxanthin that accumulates in the ripe, red fruits of L. barbarum is esterified to zeaxanthin-dipalmitate [46]. Given the very high content of zeaxanthin in RF determined here, it is possible that, similar to potatoes and apples, the esterification of zeaxanthin may be a key regulatory step in carotenoid accumulation in L. barbarum fruit.

The species differences in the carotenoid accumulation are not due to differences in the functions of key enzymes

In addition to the chromoplast development, another explanation for the large species-specific differences in the total carotenoid content of Lycium fruits could be altered by the functionality of the carotenoid biosynthetic enzymes. Given that chloroplastic carotenoids are indispensable for plant survival and to investigate this possibility, we focused on key enzymes that may be not necessary for chloroplastic carotenoid biosynthesis. These are the chromoplast-specific PSY1, CYC-B and CRTR-B2 enzymes. Protein sequence comparisons revealed up to 97% identity between species, with no insertion/deletion or frame shift mutations indicative of non-functional proteins in L. ruthenicum (Table 1). Furthermore, the protein expression in E. coli revealed that all of the enzymes from both species were equally functional in catalyzing their respective carotenoid substrates (Figure 6). Therefore, the low carotenoid content in the BF of L. ruthenicum is unlikely to be due to reduced activities of the carotenoid biosynthetic enzymes.

The biosynthesis of zeaxanthin in RF is regulated at the transcriptional level

Increasing evidence suggests that the carotenoid content in the chromoplasts is predominantly regulated at the transcriptional level [16]. The ripening of the tomato fruit is one of the best studied systems for the regulation of carotenoid biosynthesis and accumulation in the chromoplasts. Changes in the production of carotenoids associated with tomato fruit ripening are mainly controlled via the transcriptional regulation of biosynthetic genes [58]. During RF, the zeaxanthin accumulation was significantly correlated with the upregulated expression of the upstream biosynthetic genes DXS2, PSY1, PDS, ZDS, CRTISO, CYC-B, and CRTR-B2 and the very low expression of the downstream gene ZEP (Figures 4 and 5; Additional file 5). In tomatoes, DXS1 is ubiquitously expressed and shows the highest expression levels during fruit ripening, while the DXS2 transcripts are not detected in the fruit [59]. Here, the DXS1 transcript was detected at very low levels in RF, whereas DXS2 was expressed more highly, especially in S2 (Figure 5). This difference between tomatoes and L. barbarum may indicate a functional divergence of paralogous genes in the different species [60]. During RF development, the upregulation of the chromoplast-specific genes LbPSY1, LbCYC-B and LbCRTR-B2 indicated the presence of chromoplast-specific carotenoid biosynthesis in RF (Figures 4 and 5). The PDS, ZDS and CRTISO genes showed similar expression patterns during RF ripening, with a sharp increase from S1 to S2 that was maintained until S4 (Figure 5). Therefore, it is possible that the transcription of these three genes is controlled by the same mechanism or regulatory factor(s). Due to the failure in chromoplast formation during BF development, all of the abovementioned genes were expressed with generally lower levels in BF than those in RF (Figure 5; Additional file 5).

LCY-E and two P450 family hydroxylases (CYP97A29 and CYP97C11) are primarily involved in the biosynthesis of lutein [61]. The qRT-PCR results showed that these three genes exhibit a similar expression pattern in both Lycium species, being highly expressed in the leaves and expressed at low levels in the fruit (Figure 5). This was consistent with the low or no lutein accumulation in RF and BF.

Carotenoid degradation may continuously occur in the BF of L. ruthenicum

The relatively low transcript levels (but not no transcript) of the carotenoid biosynthetic genes in all four BF stages (Figure 5, Additional file 5) suggested that although the chromoplasts were not well formed in BF, the biosynthesis of the carotenoids still occurred to a small extent. However, the content of the carotenoids decreased to undetectable levels during BF ripening (Additional file 2). Therefore, we speculated that the carotenoid degradation occurred during the development of BF. In plants, the degradation of carotenoids is catalyzed by a family of CCDs, which contribute to the overall control of the cellular carotenoid content [6],[22],[62]. Arabidopsis has nine CCD family members, five of which have been classified as ABA-related AtNCEDs, and the remaining are AtCCD1, AtCCD4, AtCCD7 and AtCCD8 [63]. Of these, CCD4 has been proven to play a decisive role in the regulation of the carotenoid content in some plant organs [64], including chrysanthemum petals [23], peach fruits [24],[25] and potato tubers [26]. Interestingly, our results confirmed that LrCCD4 was highly expressed during BF development (Figures 4 and 5; Additional file 5). These results suggest that the low activity biosynthetic carotenoids were gradually degraded by the highly active LrCCD4 during BF ripening, further resulting in almost undetectable carotenoid levels in the ripe BF. In contrast, the transcripts of LbCCD4 showed a decreasing trend during RF ripening, and its expression level was much lower in RF than in BF (Figure 5, Additional file 5). Therefore, the lower rate of carotenoid degradation may be another factor for the increased carotenoid accumulation in RF compared to BF.


In conclusion, a regulatory model for the species-specific differences in carotenoid accumulation in L. barbarum and L. ruthenicum fruits has been proposed. The development of carotenoid sink organelles (chromoplasts) is likely the primary cause of the differences in carotenoid accumulation between RF and BF. In RF, based on the formation of chromoplasts, a high flux towards zeaxanthin, which is regulated at the transcriptional level, combined with a low rate of carotenoid degradation concurrently determine the observed accumulation of high levels of zeaxanthin. In BF, where the chromoplasts are not formed, small amounts of carotenoids are biosynthesized, but they are mostly degraded by LrCCD4; therefore, no carotenoids can be detected in ripe BF. The esterification of zeaxanthin in RF may be a possible regulatory step for carotenoid biosynthesis, which still requires further investigation. This study has improved our understanding of the regulatory mechanisms controlling the levels of important medicinal and nutritional compounds in Lycium.


Plant material

The L. barbarum and L. ruthenicum samples (mature leaves and fruits at four developmental stages) used in this study were collected from Zhongning County, the Ningxia Hui Autonomous Region and the Turpan Desert Botanical Garden of the Chinese Academy of Sciences, China. The samples of the four fruit developmental stages (for both RF and BF) were harvested based on the phenotype of the fruit epidermis (Figure 1): the green fruit stage (S1, 3 days before color break), the color-break stage (S2), the light-color stage (S3, 3 days after break) and the ripe fruit stage (S4, 6 days after break). For each developmental stage, more than twenty fruits were collected randomly and were then separated into three replicate groups. After harvest, each group of fruits was weighed, frozen in liquid nitrogen, and stored at −80C until further use.

Carotenoid extraction

The carotenoids were extracted from the fruits as previously described [56]. Briefly, the fruits were ground into a fine powder with liquid nitrogen and extracted three times using 5 ml of hexane/acetone/ethanol (2:1:1, v/v/v; with 0.1% butylated hydroxytoluene) via an ultrasonic treatment for 30 min until the sample was colorless. After centrifugation (4000 g for 10 min at 4C), the extracts were combined into a 50-ml tube, followed by shaking with 5 ml of NaCl-saturated solution for 1 min, and the supernatant was collected. The residue was partitioned with 5 ml of hexane and repeated three times, and all of the supernatants were combined and dried in a Vacufuge Plus vacuum concentrator (Eppendorf, Germany). Dichloromethane (2 ml) was added for the HPLC analysis of the samples that did not require saponification. For the samples that required saponification, the residue was dissolved in 2 ml of methyltert-butylether (MTBE), after which, 2 ml of a 15% (w/v) KOH/methanol solution was added for the saponification for 6 h in the dark under nitrogen [46]. After the saponification, the solutions were partitioned with 2 ml of MTBE and 4 ml of NaCl-saturated solution, and the supernatant was collected. The lower aqueous layer was repeatedly partitioned three times with 2 ml of MTBE. The supernatants were pooled, vacuum dried, and dissolved in 2 ml of dichloromethane for the HPLC analysis and the total carotenoid quantification.

HPLC analysis

The HPLC analysis was carried out on an LC-20A liquid chromatograph (Shimadzu, Japan) with two LC-20AT pumps and an SPD-M20A UV/VIS detector. All of the separations were performed using a reverse-phase C30 carotenoid column (250 4.6-mm i.d., 5-μm particle size, YMC, Japan) coupled to a 23 4.0-mm guard column. The data were acquired and processed using Shimadzu LC solution software. A binary mobile phase of methanol/acetonitrile (3:8, v/v) (A) and dichloromethane/hexane (1:1, v/v) (B) was used with the following gradient elution: 95% A and 5% B initially, decreased to 50% A in 15 min and returned to 95% A in 20 min, then maintained until 25 min. The column temperature was maintained at 30C, with a flow rate of 1 ml/min and a detection wavelength of 450nm. The quantification was performed using a calibration curve generated with commercially available β-carotene, β-cryptoxanthin and zeaxanthin standards (Sigma-Aldrich) (Additional file 6). For the preparation of the standard curves, a mixture of zeaxanthin, β-cryptoxanthin and β-carotene was dissolved into five concentrations (zeaxanthin: 2.5, 5, 25, 50 and 100 μg ml-1; β-cryptoxanthin: 0.5, 1, 5, 10 and 20 μg ml-1; β-carotene: 0.5, 1, 5, 10 and 20 μg ml-1). Three standard curves were each prepared by plotting the concentration of the carotenoid standard to its area. The regression equations and correlation coefficients (R 2) of the standard curves are shown in Additional file 7. Lutein and violaxanthin were identified via their absorption spectra and based on previous reports. All of the spectra of the identified carotenoids are listed in Additional file 8. The total carotenoid content was estimated using a spectrophotometer using zeaxanthin for the standard curve drawing.

Light microscopy

The tissues were treated for microscopy as previously described [36]. Briefly, immediately after excision with a sterile razor blade, the young fruits (S1) and ripe fruits (S4) of L. barbarum and L. ruthenicum (mesocarp only, 1-mm2 sections) were fixed in 3.5% glutaraldehyde solution for 1 h in darkness. The young fruit tissue was disrupted at 65C in a solution of disodium EDTA (EDTA-Na2; 0.1 M, pH 9.0) for 20 min, followed by maceration with clean forceps on glass microscope slides. The ripe-fruit mesocarp was disrupted in a solution of EDTA-Na2 (0.1 M, pH 9.0) at room temperature. The samples were imaged using a ZEISS Axioplan2 imaging microscope (Carl Zeiss, Germany). The images were captured with a ZEISS Axiocam MRC digital camera and Axiovision 4.6 software.

Gene cloning and sequence analysis

Total RNA was isolated from the L. barbarum and L. ruthenicum tissues using TRIzol Reagent (Invitrogen, USA). For the gene cloning for each species, the five RNA samples (from the leaf and four developmental fruits) were combined and 20 ng of the RNA was reverse transcribed into cDNA using SuperScript II (Invitrogen, USA) according to the manufacturer's instructions. A unigene library of L. barbarum (unpublished data) was searched using the carotenogenesis-related genes of tomatoes, and 25 putative carotenogenesis-related unigenes from L. barbarum were identified, namely DXS1, DXS2, PSY1, PSY2, PDS, Z-ISO, ZDS, CRTISO, LCY-B, CYC-B, LCY-E, CYP97A29, CYP97C11, CRTR-B1, CRTR-B2, ZEP, VDE, NCED1, NCED6, CCD1A, CCD4, CHRC, Or1, Or2 and HSP21 (Table 1). The primers (Additional file 9) were designed to amplify all of the genes containing a complete ORF in L. barbarum and L. ruthenicum. The PCR was performed in 25-μl mixtures containing 1 U of PrimeSTAR HS DNA Polymerase (Takara, China), 5 μl of PrimeSTAR buffer (5), 0.2 MM of each dNTP, 0.3 μM of each primer and 50 ng of cDNA. The following conditions were used: denaturation at 98C for 5 min; 30 Cycles of 98C for 10 s, 60C for 15 s, 72C for 2 min, and a final 10 min extension at 72C. The products were gel purified, ligated into the pMD19-T vector (Takara, China) and verified via sequencing (Sangon, China). The predicted coding and amino acid sequences for each gene were reciprocally compared between L. barbarum, L. ruthenicum and S. lycopersicum using DNAMAN (Table 1). The sub-cellular localization of the protein sequences was predicted using ProtComp ( (Table 1).

RNA-seq expression profiling of the carotenogenesis-related genes during the fruit ripening process (S1-S3) of RF and BF

The RNA-seq data of two replicates for the three developmental stages (S1-S3) of RF and BF were used to characterize the expression profiles of the carotenogenesis-related genes. The gene expression levels were quantified in RPKM. Twenty-five carotenogenesis-related genes were analyzed, including two genes from the MEP pathway: DXS1 and DXS2; fifteen carotenoid biosynthetic genes: PSY1, PSY2, PDS, Z-ISO, ZDS, CRTISO, LCY-B, CYC-B, LCY-E, CRTR-B1, CRTR-B2, CYP97A29, CYP97C11, ZEP and VDE; four carotenoid cleavage genes: NCED1, NCED6, CCD1A and CCD4; and four chromoplast-related genes: CHRC, Or1, Or2 and HSP21 (Additional file 5).

Functional analyses of PSY1, CYC-B and CRTR-B2 in E. coli

To test the differences between the functions of the key genes in L. barbarum and L. ruthenicum, three rate-controlling carotenoid biosynthesis genes from each species (LbPSY1/LrPSY1, LbCYC-B/LrCYC-B and LbCRTR-B2/LrCRTR-B2) were cloned into pEASY-E1 (Trans, China). The carotenoid producing plasmids were kindly provided by Dr. Norihiko Misawa. The positive controls, negative controls and gene functional assays are described in Table 2. The plasmids were transformed into E. coli BL21 (DE3) and grown overnight at 37C on Luria-Bertani (LB) solid medium with appropriate antibiotics (30 μg/ml chloramphenicol for the positive controls containing only one plasmid or 100 μg/ml ampicillin and 30 μg/ml chloramphenicol for selecting two plasmids). The selected positive colonies were shaken overnight at 37C in LB liquid medium with the appropriate antibiotics and then isopropyl β-D-1-thiogalactopyranoside (IPTG) was added to a concentration of 1.0 MM. The bacterial cultures were incubated at room temperature in darkness for 24 h for the gene expression. The cells were harvested by centrifugation (3000 g for 15 min at 4C), washed in sterile distilled water, lyophilized and used for the carotenoid extraction and HPLC analysis as described above.

Quantitative reverse transcription-PCR (qRT-PCR)

qRT-PCR was performed in five samples (the leaf and four developmental fruits) for each species. The cDNA was synthesized from 20 ng of RNA using a PrimeScript RT reagent Kit with gDNA Eraser (Perfect Real Time) (Takara, China). The qRT-PCR primers were designed using PRIMER3 software (Rozen and Skaletsky, 2000 [65]; listed in Additional file 9). The primers were designed to amplify the conserved sequences between L. barbarum and L. ruthenicum. The PCR was performed on a LightCycler@ 480 (Roche, USA) using SYBR Premix Ex Taq (Perfect Real Time) (Takara, China). The reaction mixture contained 10.0 μl of 2 SYBR Premix Ex Taq II, 2 μl of the cDNA solution (40 ng/μl), 0.8 μl of each primer (10 μM) and 6.4 μl of ddH2O, in a final volume of 20 μl. The amplification conditions were as follows: 95C for 30 s, 40Cycles of 95C for 5 s, 50C for 30 s, and 72C for 30 s. The relative expression levels were calculated using the ΔΔCT method as described by Lyi et al. (2007) [66]. All of the data were normalized to Lycium Actin1.

Accession numbers

The sequence data reported in this article can be found in GenBank under the following accession numbers: KC190187 and KC190188 for LbPDS and LrPDS, respectively; KF957678-KF957701 for DXS1, DXS2, PSY1, PSY2, Z-ISO, ZDS, CRTISO, LCY-B, CYC-B, LCY-E, CRTR-B1, CRTR-B2, CYP97A29, CYP97C11, ZEP, VDE, NCED1, NCED6, CCD1A, CCD4, CHRC, Or1, Or2 and HSP21 isolated from L. barbarum; KF957702-KF957725 for these 24 genes isolated from L. ruthenicum.

Additional files


  1. Cazzonelli CI: Goldacre review: carotenoids in nature: insights from plants and beyond. Funct Plant Biol. 2011, 38: 833-847. 10.1071/FP11192.

    Article  CAS  Google Scholar 

  2. Frank HA, Cogdell RJ: Carotenoids in photosynthesis. Photochem Photobiol. 2008, 63: 257-264. 10.1111/j.1751-1097.1996.tb03022.x.

    Article  Google Scholar 

  3. Tanaka Y, Sasaki N, Ohmiya A: Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008, 54: 733-749. 10.1111/j.1365-313X.2008.03447.x.

    Article  CAS  PubMed  Google Scholar 

  4. Grotewold E: The genetics and biochemistry of floral pigments. Annu Rev Plant Biol. 2006, 57: 761-780. 10.1146/annurev.arplant.57.032905.105248.

    Article  CAS  PubMed  Google Scholar 

  5. Xie X, Yoneyama K: The strigolactone story. Annu Rev Phytopathol. 2010, 48: 93-117. 10.1146/annurev-phyto-073009-114453.

    Article  CAS  PubMed  Google Scholar 

  6. Walter MH, Floss DS, Strack D: Apocarotenoids: hormones, mycorrhizal metabolites and aroma volatiles. Planta. 2010, 232: 1-17. 10.1007/s00425-010-1156-3.

    Article  CAS  PubMed  Google Scholar 

  7. Parry AD, Horgan R: Carotenoids and abscisic acid (ABA) biosynthesis in higher plants. Physiol Plant. 2006, 82: 320-326. 10.1111/j.1399-3054.1991.tb00100.x.

    Article  Google Scholar 

  8. Havaux M: Carotenoid oxidation products as stress signals in plants. Plant J. 2014, 79 (4): 597-606. 10.1111/tpj.12386.

    Article  CAS  PubMed  Google Scholar 

  9. Rao A, Rao L: Carotenoids and human health. Pharmacol Res. 2007, 55: 207-216. 10.1016/j.phrs.2007.01.012.

    Article  CAS  PubMed  Google Scholar 

  10. Fraser PD, Bramley PM: The biosynthesis and nutritional uses of carotenoids. Prog Lipid Res. 2004, 43: 228-265. 10.1016/j.plipres.2003.10.002.

    Article  CAS  PubMed  Google Scholar 

  11. Weber D, Grune T: The contribution of β-carotene to vitamin A supply of humans. Mol Nutr Food Res. 2012, 56: 251-258. 10.1002/mnfr.201100230.

    Article  CAS  PubMed  Google Scholar 

  12. Sabour-Pickett S, Nolan JM, Loughman J, Beatty S: A review of the evidence germane to the putative protective role of the macular carotenoids for age-related macular degeneration. Mol Nutr Food Res. 2012, 56: 270-286. 10.1002/mnfr.201100219.

    Article  CAS  PubMed  Google Scholar 

  13. Ruiz-Sola M, Rodrguez-Concepcin M: Carotenoid biosynthesis in Arabidopsis: A colorful pathway. The Arabidopsis Book. Edited by: Torii K. The American Society of Plant Biologists, Rockville. 2012, e0158-

    Google Scholar 

  14. R Mer S, Fraser PD: Recent advances in carotenoid biosynthesis, regulation and manipulation. Planta. 2005, 221: 305-308. 10.1007/s00425-005-1533-5.

    Article  Google Scholar 

  15. Hirschberg J: Carotenoid biosynthesis in flowering plants. Curr Opin Plant Biol. 2001, 4: 210-218. 10.1016/S1369-5266(00)00163-1.

    Article  CAS  PubMed  Google Scholar 

  16. Ohmiya A: Qualitative and quantitative control of carotenoid accumulation in flower petals. Sci Hortic. 2013, 163: 10-19. 10.1016/j.scienta.2013.06.018.

    Article  CAS  Google Scholar 

  17. Hannoufa A, Hossain Z: Regulation of carotenoid accumulation in plants. Biocatal Agric Biotechnol. 2012, 1: 198-202.

    CAS  Google Scholar 

  18. Galpaz N, Ronen G, Khalfa Z, Zamir D, Hirschberg J: A chromoplast-specific carotenoid biosynthesis pathway is revealed by cloning of the tomato white-flower locus. Plant Cell. 2006, 18: 1947-1960. 10.1105/tpc.105.039966.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Ronen G, Carmel-Goren L, Zamir D, Hirschberg J: An alternative pathway to β-carotene formation in plant chromoplasts discovered by map-based cloning of Beta and old-gold color mutations in tomato. Proc Natl Acad Sci U S A. 2000, 97: 11102-11107. 10.1073/pnas.190177497.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Ronen G, Cohen M, Zamir D, Hirschberg J: Regulation of carotenoid biosynthesis during tomato fruit development: expression of the gene for lycopene epsilon-cyclase is down-regulated during ripening and is elevated in the mutant Delta . Plant J. 1999, 17: 341-351. 10.1046/j.1365-313X.1999.00381.x.

    Article  CAS  PubMed  Google Scholar 

  21. Giuliano G, Bartley GE, Scolnik PA: Regulation of carotenoid biosynthesis during tomato development. Plant Cell. 1993, 5: 379-387. 10.1105/tpc.5.4.379.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Ohmiya A: Carotenoid cleavage dioxygenases and their apocarotenoid products in plants. Plant Biotechnol. 2009, 26: 351-358. 10.5511/plantbiotechnology.26.351.

    Article  CAS  Google Scholar 

  23. Ohmiya A, Kishimoto S, Aida R, Yoshioka S, Sumitomo K: Carotenoid cleavage dioxygenase (CmCCD4a) contributes to white color formation in chrysanthemum petals. Plant Physiol. 2006, 142: 1193-1201. 10.1104/pp.106.087130.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Adami M, De Franceschi P, Brandi F, Liverani A, Giovannini D, Rosati C, Dondini L, Tartarini S: Identifying a carotenoid cleavage dioxygenase (ccd4) gene controlling Yellow/White fruit flesh color of peach. Plant Mol Biol Rep. 2013, 31: 1166-1175. 10.1007/s11105-013-0628-6.

    Article  CAS  Google Scholar 

  25. Brandi F, Bar E, Mourgues F, Horvth G, Turcsi E, Giuliano G, Liverani A, Tartarini S, Lewinsohn E, Rosati C: Study of Redhaven peach and its white-fleshed mutant suggests a key role of CCD4 carotenoid dioxygenase in carotenoid and norisoprenoid volatile metabolism. BMC Plant Biol. 2011, 11: 24-10.1186/1471-2229-11-24.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Campbell R, Ducreux LJ, Morris WL, Morris JA, Suttle JC, Ramsay G, Bryan GJ, Hedley PE, Taylor MA: The metabolic and developmental roles of carotenoid cleavage dioxygenase4 from potato. Plant Physiol. 2010, 154: 656-664. 10.1104/pp.110.158733.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Lu S, Van Eck J, Zhou X, Lopez AB, O halloran DM, Cosman KM, Conlin BJ, Paolillo DJ, Garvin DF, Vrebalov J: The cauliflower Or gene encodes a DnaJ cysteine-rich domain-containing protein that mediates high levels of β-carotene accumulation. Plant Cell. 2006, 18: 3594-3605. 10.1105/tpc.106.046417.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Wang Y, Yang Y, Fei Z, Yuan H, Fish T, Thannhauser TW, Mazourek M, Kochian LV, Wang X, Li L: Proteomic analysis of chromoplasts from six crop species reveals insights into chromoplast function and development. J Exp Bot. 2013, 64 (4): 949-961. 10.1093/jxb/ers375.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Cookson P, Kiano J, Shipton C, Fraser P, Romer S, Schuch W, Bramley P, Pyke K: Increases in cell elongation, plastid compartment size and phytoene synthase activity underlie the phenotype of the high pigment-1 mutant of tomato. Planta. 2003, 217: 896-903. 10.1007/s00425-003-1065-9.

    Article  CAS  PubMed  Google Scholar 

  30. Davuluri GR, Van Tuinen A, Fraser PD, Manfredonia A, Newman R, Burgess D, Brummell DA, King SR, Palys J, Uhlig J: Fruit-specific RNAi-mediated suppression of DET1 enhances carotenoid and flavonoid content in tomatoes. Nat Biotechnol. 2005, 23: 890-895. 10.1038/nbt1108.

    Article  CAS  PubMed  Google Scholar 

  31. Mustilli AC, Fenzi F, Ciliento R, Alfano F, Bowler C: Phenotype of the tomato high pigment-2 mutant is caused by a mutation in the tomato homolog of DEETIOLATED1 . Plant Cell. 1999, 11: 145-157. 10.1105/tpc.11.2.145.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Wang S, Liu J, Feng Y, Niu X, Giovannoni J, Liu Y: Altered plastid levels and potential for improved fruit nutrient content by downregulation of the tomato DDB1-interacting protein CUL4. Plant J. 2008, 55: 89-103. 10.1111/j.1365-313X.2008.03489.x.

    Article  CAS  PubMed  Google Scholar 

  33. Liu Y, Roof S, Ye Z, Barry C, Van Tuinen A, Vrebalov J, Bowler C, Giovannoni J: Manipulation of light signal transduction as a means of modifying fruit nutritional quality in tomato. Proc Natl Acad Sci U S A. 2004, 101: 9897-9902. 10.1073/pnas.0400935101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Giliberto L, Perrotta G, Pallara P, Weller JL, Fraser PD, Bramley PM, Fiore A, Tavazza M, Giuliano G: Manipulation of the blue light photoreceptor cryptochrome 2 in tomato affects vegetative development, flowering time, and fruit antioxidant content. Plant Physiol. 2005, 137: 199-208. 10.1104/pp.104.051987.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Powell AL, Nguyen CV, Hill T, Cheng KL, Figueroa-Balderas R, Aktas H, Ashrafi H, Pons C, Fernndez-Muoz R, Vicente A: Uniform ripening encodes a Golden 2-like transcription factor regulating tomato fruit chloroplast development. Science. 2012, 336: 1711-1715. 10.1126/science.1222218.

    Article  CAS  PubMed  Google Scholar 

  36. Pan Y, Bradley G, Pyke K, Ball G, Lu C, Fray R, Marshall A, Jayasuta S, Baxter C, Van Wijk R: Network inference analysis identifies an APRR2-like gene linked to pigment accumulation in tomato and pepper fruits. Plant Physiol. 2013, 161: 1476-1485. 10.1104/pp.112.212654.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Azari R, Tadmor Y, Meir A, Reuveni M, Evenor D, Nahon S, Shlomo H, Chen L, Levin I: Light signaling genes and their manipulation towards modulation of phytonutrient content in tomato fruits. Biotechnol Adv. 2010, 28: 108-118. 10.1016/j.biotechadv.2009.10.003.

    Article  CAS  PubMed  Google Scholar 

  38. Potterat O: Goji (Lycium barbarum and L. chinense): phytochemistry, pharmacology and safety in the perspective of traditional uses and recent popularity. Planta Med. 2010, 76: 7-19. 10.1055/s-0029-1186218.

    Article  CAS  PubMed  Google Scholar 

  39. Chang RC-C, So K-F: Use of anti-aging herbal medicine, Lycium barbarum, against aging-associated diseases. What do we know so far?. Cell Mol Neurobiol. 2008, 28: 643-652. 10.1007/s10571-007-9181-x.

    Article  CAS  PubMed  Google Scholar 

  40. Yao X, Peng Y, Xu LJ, Li L, Wu QL, Xiao PG: Phytochemical and biological studies of Lycium medicinal plants. Chem Biodivers. 2011, 8: 976-1010. 10.1002/cbdv.201000018.

    Article  CAS  PubMed  Google Scholar 

  41. Wang C, Chang S, Inbaraj BS, Chen B: Isolation of carotenoids, flavonoids and polysaccharides from Lycium barbarum L. and evaluation of antioxidant activity. Food Chem. 2010, 120: 184-192. 10.1016/j.foodchem.2009.10.005.

    Article  CAS  Google Scholar 

  42. Zheng J, Ding C, Wang L, Li G, Shi J, Li H, Wang H, Suo Y: Anthocyanins composition and antioxidant activity of wild Lycium ruthenicum Murr. from Qinghai-Tibet Plateau. Food Chem. 2011, 126: 859-865. 10.1016/j.foodchem.2010.11.052.

    Article  CAS  Google Scholar 

  43. Peng Q, Song J, Lv X, Wang Z, Huang L, Du Y: Structural Characterization of an Arabinogalactan-Protein from the Fruits of Lycium ruthenicum . J Agric Food Chem. 2012, 60: 9424-9429. 10.1021/jf302619c.

    Article  CAS  PubMed  Google Scholar 

  44. Altintas A, Kosar M, Kirimer N, Baser K, Demirci B: Composition of the essential oils of Lycium barbarum and L. ruthenicum fruits. Chem Nat Compd. 2006, 42: 24-25. 10.1007/s10600-006-0028-3.

    Article  CAS  Google Scholar 

  45. Zeng S, Wu M, Zou C, Liu X, Shen X, Hayward A, Liu C, Wang Y: Comparative analysis of anthocyanin biosynthesis during fruit development in two Lycium species. Physiol Plant. 2014, 150 (4): 505-516. 10.1111/ppl.12131.

    Article  CAS  PubMed  Google Scholar 

  46. Inbaraj BS, Lu H, Hung CF, Wu WB, Lin CL, Chen BH: Determination of carotenoids and their esters in fruits of Lycium barbarum Linnaeus by HPLC-DAD-APCI-MS. J Pharm Biomed Anal. 2008, 47: 812-818. 10.1016/j.jpba.2008.04.001.

    Article  CAS  PubMed  Google Scholar 

  47. Lois LM, Rodrguez-Concepcin M, Gallego F, Campos N, Boronat A: Carotenoid biosynthesis during tomato fruit development: regulatory role of 1-deoxy-d-xylulose 5-phosphate synthase. Plant J. 2000, 22: 503-513. 10.1046/j.1365-313x.2000.00764.x.

    Article  CAS  PubMed  Google Scholar 

  48. Ikoma Y, Matsumoto H, Kato M: The characteristics of carotenoid biosynthesis in Citrus fruit. Jpn Agric Res Q. 2014, 48: 9-16. 10.6090/jarq.48.9.

    Article  CAS  Google Scholar 

  49. Li Q, Farre G, Naqvi S, Breitenbach J, Sanahuja G, Bai C, Sandmann G, Capell T, Christou P, Zhu C: Cloning and functional characterization of the maize carotenoid isomerase and β-carotene hydroxylase genes and their regulation during endosperm maturation. Transgenic Res. 2010, 19: 1053-1068. 10.1007/s11248-010-9381-x.

    Article  CAS  PubMed  Google Scholar 

  50. Li L, Yuan H: Chromoplast biogenesis and carotenoid accumulation. Arch Biochem Biophys. 2013, 539: 102-109. 10.1016/

    Article  CAS  PubMed  Google Scholar 

  51. Bai C, Rivera SM, Medina V, Alves R, Vilaprinyo E, Sorribas A, Canela R, Capell T, Sandmann G, Christou P: An in vitro system for the rapid functional characterization of genes involved in carotenoid biosynthesis and accumulation. Plant J. 2014, 77: 464-475. 10.1111/tpj.12384.

    Article  CAS  PubMed  Google Scholar 

  52. Kilambi HV, Kumar R, Sharma R, Sreelakshmi Y: Chromoplast specific carotenoid associated protein appears to be important for enhanced accumulation of carotenoids in hp1 tomato fruits. Plant Physiol. 2013, 161: 2085-2101. 10.1104/pp.112.212191.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Neta-Sharir I, Isaacson T, Lurie S, Weiss D: Dual role for tomato heat shock protein 21: protecting photosystem II from oxidative stress and promoting color changes during fruit maturation. Plant Cell. 2005, 17: 1829-1838. 10.1105/tpc.105.031914.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Derure J, Romer S, D harlingue A, Backhaus RA, Kuntz M, Camara B: Fibril assembly and carotenoid overaccumulation in chromoplasts: a model for supramolecular lipoprotein structures. Plant Cell. 1994, 6: 119-133. 10.1105/tpc.6.1.119.

    Article  Google Scholar 

  55. Vishnevetsky M, Ovadis M, Vainstein A: Carotenoid sequestration in plants: the role of carotenoid-associated proteins. Trends Plant Sci. 1999, 4: 232-235. 10.1016/S1360-1385(99)01414-4.

    Article  PubMed  Google Scholar 

  56. Fernndez-Orozco R, Gallardo-Guerrero L, Hornero-Mndez D: Carotenoid profiling in tubers of different potato (Solanum sp) cultivars: Accumulation of carotenoids mediated by xanthophyll esterification. Food Chem. 2013, 141: 2864-2872. 10.1016/j.foodchem.2013.05.016.

    Article  Google Scholar 

  57. Delgado-Pelayo R, Gallardo-Guerrero L, Hornero-Mndez D: Chlorophyll and carotenoid pigments in the peel and flesh of commercial apple fruits varieties. Food Res Int 2014. doi:10.1016/j.foodres.2014.03.025.,

  58. Bramley PM: Regulation of carotenoid formation during tomato fruit ripening and development. J Exp Bot. 2002, 53: 2107-2113. 10.1093/jxb/erf059.

    Article  CAS  PubMed  Google Scholar 

  59. Paetzold H, Garms S, Bartram S, Wieczorek J, Urs-Gracia E-M, Rodrguez-Concepcin M, Boland W, Strack D, Hause B, Walter MH: The isogene 1-deoxy-D-xylulose 5-phosphate synthase 2 controls isoprenoid profiles, precursor pathway allocation, and density of tomato trichomes. Mol Plant. 2010, 3: 904-916. 10.1093/mp/ssq032.

    Article  CAS  PubMed  Google Scholar 

  60. Walter MH, Hans J, Strack D: Two distantly related genes encoding 1-deoxy-d-xylulose 5-phosphate synthases: differential regulation in shoots and apocarotenoid-accumulating mycorrhizal roots. Plant J. 2002, 31: 243-254. 10.1046/j.1365-313X.2002.01352.x.

    Article  CAS  PubMed  Google Scholar 

  61. Stigliani AL, Giorio G, DAmbrosio C: Characterization of P450 carotenoid β-and ε-hydroxylases of tomato and transcriptional regulation of xanthophyll biosynthesis in root, leaf, petal and fruit. Plant Cell Physiol. 2011, 52: 851-865. 10.1093/pcp/pcr037.

    Article  CAS  PubMed  Google Scholar 

  62. Auldridge ME, McCarty DR, Klee HJ: Plant carotenoid cleavage oxygenases and their apocarotenoid products. Curr Opin Plant Biol. 2006, 9: 315-321. 10.1016/j.pbi.2006.03.005.

    Article  CAS  PubMed  Google Scholar 

  63. Tan BC, Joseph LM, Deng WT, Liu L, Li QB, Cline K, McCarty DR: Molecular characterization of the Arabidopsis 9-cis epoxycarotenoid dioxygenase gene family. Plant J. 2003, 35: 44-56. 10.1046/j.1365-313X.2003.01786.x.

    Article  CAS  PubMed  Google Scholar 

  64. Huang F-C, Molnr P, Schwab W: Cloning and functional characterization of carotenoid cleavage dioxygenase 4 genes. J Exp Bot. 2009, 60: 3011-3022. 10.1093/jxb/erp137.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    CAS  PubMed  Google Scholar 

  66. Lyi SM, Zhou X, Kochian LV, Li L: Biochemical and molecular characterization of the homocysteine S-methyltransferase from broccoli (Brassica oleracea var. italica). Phytochemistry. 2007, 68: 1112-1119. 10.1016/j.phytochem.2007.02.007.

    Article  CAS  PubMed  Google Scholar 

  67. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  68. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.

    CAS  PubMed  Google Scholar 

Download references


This work was supported by grants from the National Natural Science Foundation of China (31100223, 31200242, 31470391) and the Instrument Developing Project of the Chinese Academy of Sciences (Grant No. YZ201227). We would like to thank Dr. Norihiko Misawa for providing plasmids, Drs. Xiaoyi Wei and David Ow for helping with facilities, and Dr. Alice Hayward from The University of Queensland for polishing this manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ying Wang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

YL participated in the sampling, carried out the experiments and statistical analysis, and drafted the manuscript. SZ participated in the sampling and data analysis and drafted the manuscript. WS participated in the sequence alignment. MW participated in the HPLC analysis. WH participated in the sampling and helped in field experiment. XS helped in carotenoid extraction. YW conceived of the study, designed and coordinated the research. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Carotenoid concentrations (μg g -1fresh weight) in all four developmental stages (S1-S4) of L. barbarum and L. ruthenicum fruits. (PDF 108 KB)


Additional file 2: Accumulation profiles of carotenoids in ripening RF revealed by HPLC. Peak identification: (1) violaxanthin, (2) lutein, (3) β-carotene, (4) β-cryptoxanthin, (5) zeaxanthin. (PNG 82 KB)


Additional file 3: Accumulation profiles of carotenoids in ripening BF revealed by HPLC. Peak identification: (1) violaxanthin, (2) lutein, (3) β-carotene, (6) unidentified, (7) unidentified. (PNG 91 KB)


Additional file 4: Phylogenetic trees for the carotenogenesis-ralated proteins from L. barbarum (Lb), L. ruthenicum (Lr) and other organisms. The trees were constructed by MEGA 5.1 program [67] using the neighbor-joining method [68].The proteins in Solanum lycopersicum (Sl), Capsicum annuum (Ca), Arabidopsis thaliana (At), Nicotiana tabacum (Nt) and Brassica oleracea (Bo) include SlDXS1 (FN424051), SlDXS2 (FN424052), SlPSY1 (ABU40772), SlPSY2 (ABU40771), SlPDS (AGO05926), SlZ-ISO (XP_004252966), SlZDS (AGO05927), SlCRTISO (AAL91366), SlLCY-B (NP_001234226), SlCYC-B (AAG21133),SlLCY-E (Y14387), SlCRTR-B1 (CAB55625), SlCRTR-B2 (CAB55626), SlCYP97A29 (ACJ25969), SlCYP97C11 (ACJ25968), SlZEP (P93236), SlVDE (ACM92036), SlNCED1 (CAB10168), SlNCED6 (XP_004240215), SlCCD1A (AAT68187), SlCCD1B (AAT68188), SlCCD4 (XP_004246004), SlCHRC (ABC42191), SlHSP21 (AAB07023); CaTKT2 (CAA75778), CaPSY (P37272), CaPDS (P80093), CaZDS (Q9SMJ3), CaLCY-B (ADH04271), CaCCS (CAA54495), CaCRTR-B1 (CAA70888), CaCRTR-B2 (CAA70427), CaZEP (Q96375), CaFBN (CAA50750); AtCLA1 (AAC49368), AtPSY (AAM62787), AtPDS (AAL15300), AtZ-ISO (NP_001117264), AtZDS (AAM63349), AtCRTISO (NP_172167), AtLCY-B (AAB53337), AtLCY-E (AAB53336), AtCRTR-B2 (NP_200070), AtCRTR-B1 (AAC49443), AtCYP97A3 (AAL08302), AtCYP97C1 (AAM13903), AtZEP (BAB08942), AtVDE (AAC50032), AtNCED2 (AT4G18350), AtNCED3 (AT3G14440), AtNCED5 (AT1G30100), AtNCED6 (AT3G24220), AtNCED9 (AT1G78390), AtCCD1 (AT3G63520), AtCCD4 (AT4G19170), AtCCD7 (AT2G44990), AtCCD8 (AT4G32810), AtFBN (NP_192311), AtHSP21 (NP_194497); NtOr (AEV23056) and BoOr (ABH07405). (PDF 159 KB)


Additional file 5: RNA-seq expression profiles of twenty-five carotenogenesis-related genes during fruit ripening (S1-S3) of RF and BF, and the expression profiles of orthologous genes in tomato. (XLS 38 KB)

Additional file 6: Chromatogram of standard mixture. (PNG 124 KB)

Additional file 7: Standard curves of zeaxanthin, β-cryptoxanthin and β-carotene. (PNG 1015 KB)


Additional file 8: Spectra of the carotenoids discussed in the sutdy. (1) violaxanthin, (2) lutein, (3) β-carotene, (4) β-cryptoxanthin, (5) zeaxanthin, (6) unidentified, (7) unidentified. (PNG 226 KB)

Additional file 9: Specific primers used for gene cloning and qRT-PCR. (PDF 62 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Zeng, S., Sun, W. et al. Comparative analysis of carotenoid accumulation in two goji (Lycium barbarum L. and L. ruthenicumMurr.) fruits. BMC Plant Biol 14, 269 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Carotenoids
  • Chromoplast
  • Fruit development
  • Gene expression
  • Lycium barbarum
  • L. ruthenicum