Comparative analysis of carotenoid accumulation in two goji (Lycium barbarum L. and L. ruthenicumMurr.) fruits
© Liu et al.; licensee BioMed Central Ltd. 2014
Received: 14 February 2014
Accepted: 29 September 2014
Published: 16 December 2014
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 ,. When accumulated in the chromoplasts of flowers and fruits, carotenoids act as visual attractants for pollinating insects and seed-dispersing animals ,. Furthermore, carotenoids are the precursors of important apocarotenoids, such as volatile flavor/aroma terpenes, and the growth regulators abscisic acid (ABA) and strigolactone -. Recently, oxidized products from plant carotenoids have been implicated as signals induced by environmental stressors . In addition to these biological functions, carotenoids serve as major micronutrients in the human diet ,. In particular, β-carotene, β-carotene and β-cryptoxanthin are precursors for vitamin A biosynthesis , while lutein and zeaxanthin slow aging-related damage to the retina .
During the past two decades, the carotenoid biosynthetic pathway in plants has been well elucidated -. Previous studies have shown that there are at least three mechanisms that regulate carotenoid accumulation in the chromoplasts . 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 . 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 -.
Second, carotenoid degradation by carotenoid cleavage dioxygenases (CCDs) may be central to determining the final carotenoid concentrations in chromoplasts . 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 . 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 -.
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 . When the Or gene was transformed into Arabidopsis, the Arabidopsis calli exhibited an orange color with chromoplast formation . In tomatoes, the perturbed activity of several light-signal-related genes, including UV-Damaged DNA Binding Protein1 (DDB1) , De-Etiolated1 (DET1) ,, Cullin4 (CUL4) , HY5, COP1LIKE , Cryptochrome 2 (CRY2) , Golden 2-Like (GLK2)  and Arabidopsis Pseudo Response Regulator2-Like (APRR2-Like) , caused changes in the plastid number and size, which indirectly affected the concentrations of the carotenoids and other phytonutrients in the ripening fruits .
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 . 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 . 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 ,. L. ruthenicum is another TCM used for the treatment of heart disease, abnormal menstruation and menopause . The functional compounds in the black fruit (BF) of L. ruthenicum are primarily comprised of anthocyanins, essential oils and polysaccharides -. 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 ,. 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
Light microscopy of the green fruits (S1) and ripe fruits (S4)
Isolation of the carotenogenesis-related genes from L. barbarum and L. ruthenicum
Sequence information of the carotenogenesis-related genes from L. barbarum and L. ruthenicum
Protein identity Lb vs Sl c
(bp a, Lb/Lr)
Lb vs Lr (%)
(accession number of Sl) d
Lr vs Sl
The comparative RNA-seq profile of the carotenogenesis-related genes in the ripening fruits of the two Lyciumspecies
Comparative analysis of the carotenogenesis-related gene expression in L. barbarum and L. ruthenicumfruits via qRT-PCR
As the first enzyme in the MEP pathway, DXS was shown to be a regulatory enzyme in tomato fruit carotenogenesis . 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
Construct design for enzymatic assays in E. coli
Genes to be analysed
Plasmids contained in E. coli (BL21) and the carotenoid abeing produced (in parentheses)
pACCRT-E + pEASY-E1
pACCRT-E + pEASY-LbPSY1 or pACCRT-E + pEASY-LrPSY1
pACCRT-EIB + pEASY-E1
pACCRT-EIB + pEAST-LbCYC-B or pACCRT-EIB + pEASY-LrCYC-B
pACCAR16ΔcrtX + pEASY-E1
pACCAR16ΔcrtX + pEASY-LbCRTR-B2 or pACCAR16ΔcrtX + pEASY-LrCRTR-B2
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 ,. Unlike RF, the ripe fruit of L. ruthenicum is deep purple in color with a high petunidin content produced by the anthocyanin pathway . 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) . 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 . In the Or cauliflower mutant, the failure in chromoplast formation blocked the biosynthesis and accumulation of carotenoids . In transgenic Or-overexpressing Arabidopsis and rice, the chromoplast differentiation occurred in the calli of both species and induced the biosynthesis of the carotenoids ,. In addition, CHRC and HSP21 were also shown to play significant roles in chromoplast development and carotenoid storage ,.
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 . Distinct from the plastoglobules in the chloroplasts, these fibrils are characterized by a high homogeneity of apolar compounds, most of which are esterified xanthophylls . 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 . Recently, this viewpoint was also verified in apples (Malus x domestica Borkh) . The majority of the zeaxanthin that accumulates in the ripe, red fruits of L. barbarum is esterified to zeaxanthin-dipalmitate . 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 . 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 . 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 . 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 . 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 . 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 ,,. 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 . Of these, CCD4 has been proven to play a decisive role in the regulation of the carotenoid content in some plant organs , including chrysanthemum petals , peach fruits , and potato tubers . 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.
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.
The carotenoids were extracted from the fruits as previously described . 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 . 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.
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.
The tissues were treated for microscopy as previously described . 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 (http://www.softberry.com/berry.phtml) (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 ; 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) . All of the data were normalized to Lycium Actin1.
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.
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.
- Cazzonelli CI: Goldacre review: carotenoids in nature: insights from plants and beyond. Funct Plant Biol. 2011, 38: 833-847. 10.1071/FP11192.View ArticleGoogle Scholar
- Frank HA, Cogdell RJ: Carotenoids in photosynthesis. Photochem Photobiol. 2008, 63: 257-264. 10.1111/j.1751-1097.1996.tb03022.x.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Grotewold E: The genetics and biochemistry of floral pigments. Annu Rev Plant Biol. 2006, 57: 761-780. 10.1146/annurev.arplant.57.032905.105248.View ArticlePubMedGoogle Scholar
- Xie X, Yoneyama K: The strigolactone story. Annu Rev Phytopathol. 2010, 48: 93-117. 10.1146/annurev-phyto-073009-114453.View ArticlePubMedGoogle Scholar
- Walter MH, Floss DS, Strack D: Apocarotenoids: hormones, mycorrhizal metabolites and aroma volatiles. Planta. 2010, 232: 1-17. 10.1007/s00425-010-1156-3.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Havaux M: Carotenoid oxidation products as stress signals in plants. Plant J. 2014, 79 (4): 597-606. 10.1111/tpj.12386.View ArticlePubMedGoogle Scholar
- Rao A, Rao L: Carotenoids and human health. Pharmacol Res. 2007, 55: 207-216. 10.1016/j.phrs.2007.01.012.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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
- R Mer S, Fraser PD: Recent advances in carotenoid biosynthesis, regulation and manipulation. Planta. 2005, 221: 305-308. 10.1007/s00425-005-1533-5.View ArticleGoogle Scholar
- Hirschberg J: Carotenoid biosynthesis in flowering plants. Curr Opin Plant Biol. 2001, 4: 210-218. 10.1016/S1369-5266(00)00163-1.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Hannoufa A, Hossain Z: Regulation of carotenoid accumulation in plants. Biocatal Agric Biotechnol. 2012, 1: 198-202.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohmiya A: Carotenoid cleavage dioxygenases and their apocarotenoid products in plants. Plant Biotechnol. 2009, 26: 351-358. 10.5511/plantbiotechnology.26.351.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Li L, Yuan H: Chromoplast biogenesis and carotenoid accumulation. Arch Biochem Biophys. 2013, 539: 102-109. 10.1016/j.abb.2013.07.002.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.,Google Scholar
- Bramley PM: Regulation of carotenoid formation during tomato fruit ripening and development. J Exp Bot. 2002, 53: 2107-2113. 10.1093/jxb/erf059.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.