Distinct expression and function of carotenoid metabolic genes and homoeologs in developing wheat grains
BMC Plant Biology volume 16, Article number: 155 (2016)
β-carotene, the most active provitamin A molecule produced by plants, plays important roles in human nutrition and health. β-carotene does not usually accumulate in the endosperm (i.e. flour) of mature wheat grains, which is a major food source of calories for humans. Therefore, enriching β-carotene accumulation in wheat grain endosperm will enable a sustainable dietary supplementation of provitamin A. Several metabolic genes affecting β-carotene accumulation have already been isolated from wheat, including phytoene synthase 1 (PSY1), lycopene ε-cyclase (LCYe) and carotenoid β-ring hydroxylase1/2 (HYD1/2).
In this work, we cloned and biochemically characterized two carotenoid cleavage dioxygenases (CCDs), CCD1 and CCD4, from wheat. While CCD1 homoeologs cleaved β-apo-8′-carotenal, β-carotene, lutein and zeaxanthin into apocarotenoid products, CCD4 homoeologs were inactive towards these substrates in in vitro assays. When analyzed by real-time qPCR, PSY1, LCYe, HYD1/2 and CCD1/4 homoeologs showed distinct expression patterns in vegetative tissues and sections of developing tetraploid and hexaploid wheat grains, suggesting that carotenoid metabolic genes and homoeologs are differentially regulated at the transcriptional level in wheat.
The CCD1/4 enzyme activity and the spatial-temporal gene expression data provide critical insights into the specific carotenoid metabolic gene homoeologs that control β-carotene accumulation in wheat grain endosperm, thus establishing the knowledge base for generation of wheat varieties with enhanced β-carotene in the endosperm through breeding and genome editing approaches.
Wheat (Triticum spp.), rice and maize constitute the three most widely consumed cereal grains worldwide. The color of wheat grain endosperm (i.e. flour) is determined largely by carotenoid pigments and has been selected according to consumer’s preference during the history of wheat breeding. Tetraploid durum wheat (T. turgidum) is used for making pasta and couscous and has been selected for increased concentrations of yellow pigments. Lutein, a non-provitamin A carotenoid exhibiting yellow color, is the major carotenoid molecule present in the endosperm of the mature tetraploid wheat grains. By contrast, hexaploid bread wheat (T. aestivum) has been selected for white flour, resulting in low concentrations of total carotenoids in the endosperm of the mature hexaploid wheat grains [1–3]. Therefore, wheat grain endosperm generally lacks provitamin A carotenoids for conversion into vitamin A in mammals.
Since humans cannot synthesize vitamin A from basic hydrocarbon building blocks, this essential micronutrient must be obtained from dietary sources, in the form of preformed vitamin A or provitamin A [4–7]. Considering the indispensable role of vitamin A in human nutrition and the importance of wheat flour in supplying dietary calories, enhancing provitamin A accumulation in wheat grain endosperm holds great promise for staple food-based provision/supplementation of vitamin A nutrition. Among the carotenoid molecules that possess provitamin A activities, β-carotene can be converted into vitamin A most efficiently due to its possession of two unmodified β-ionone rings (a structural feature that is essential for vitamin A activity)  (Fig. 1). As such, understanding how carotenoid, particularly β-carotene, accumulation is controlled in wheat grain endosperm is crucial for vitamin A biofortification in this tissue.
Several metabolic genes controlling β-carotene accumulation have been isolated from wheat, including phytoene synthase 1 (PSY1) that channels carbon fluxes to carotenoid biosynthesis, lycopene ε-cyclase (LCYe) that diverts the common biosynthetic precursor lycopene to a competing pathway of lutein biosynthesis, and β-carotene hydroxylase1/2 (HYD1/2) that enable turnover of β-carotene to downstream hydroxylated products [8–10] (Fig. 1).
Besides hydroxylation, β-carotene and its derivatives can also be modified via conversion into smaller molecules (aka. apocarotenoids) by carotenoid cleavage dioxygenases (CCDs) (Fig. 1). Of the four plant CCDs (including CCD1, CCD4, CCD7 and CCD8), coupled actions of CCD7 and CCD8 are required for biosynthesis of strigolactones from 9-cis-β-carotene in root tissues . CCD1s identified from different plant species localize to the cytosol and have demonstrated cleavage activities towards all-trans-β-carotene (β-carotene hereafter) and additional linear or cyclic (apo)carotenoid substrates at different double bond positions (e.g. C9-C10/C9′-C10′, C7-C8/C7′-C8′ and C5-C6/C5′-C6′) in in vitro enzyme assays [12–32]. Among the CCD1s identified from monocots, rice and maize CCD1s (OsCCD1 and ZmCCD1) have been functionally characterized [15, 16, 30, 33]. Both enzymes cleave linear (e.g. lycopene) and cyclic (e.g. β-carotene and zeaxanthin) carotenoids symmetrically at C9-C10/C9′-C10′ or C7-C8/C7′-C8′ positions [15, 16, 30]. However, overexpression of OsCCD1 in Golden Rice endosperm (engineered for high β-carotene accumulation) did not change its β-carotene content significantly, arguing against a role of OsCCD1 in carotenoid cleavage in rice endosperm . It remains to be determined whether CCD1 enzymes can transform carotenoids into apocarotenoid products in other cereal grains, such as wheat.
In contrast to the cytosolic localization of CCD1s, CCD4 proteins contain transit peptides for targeting to the plastids [24, 35, 36]. Correlations between reduced CCD4 gene expression/activity and increased carotenoid accumulation or reduced apocarotenoid volatile emission in floral, fruit or storage tissues have been observed in Arabidopsis, chrysanthemum, grape, peach, potato and saffron, suggesting a role of CCD4s in carotenoid cleavage in these plants [19, 21, 24, 37–40]. However, increased violaxanthin accumulation was detected in RNAi knockdown lines of StCCD4 (potato CCD4), even though StCCD4 was unable to cleave violaxanthin in enzyme assays, suggesting that a reduction in CCD4 activity may not lead to the accumulation of its immediate substrates in planta . When assayed in vitro, the apple and chrysanthemum CCD4s showed substantial cleavage of β-carotene at C9-C10/C9′-C10′, while the Arabidopsis, rose and osmanthus CCD4s had little conversion of β-carotene into the apocarotenoid products [18, 42]. Although CsCCD4a/b/c of saffron (Crocus sativus) all cut β-carotene at C9-C10/C9′-C10′, minor cleavage activities of β-carotene at C7-C8/C7′-C8′ were only observed for CsCCD4c [24, 43]. Lutein, neoxanthin and violaxanthin may also serve as substrates for CsCCD4c as suggested by the carotenoid profiles of Nicotiana benthamiana plants transiently overexpressing CsCCD4c .
Since polyploid wheat has several subgenomes (AA, BB and DD genomes), multiple homoeologs are present for the carotenoid metabolic genes in wheat. While tetraploid wheat (genomes AABB) was formed from the hybridization of T. urartu (AA) and an unknown species from the Sitopsis group (BB) about 500,000 years ago, hexaploid wheat (genomes AABBDD) originated from the hybridization of tetraploid wheat (AABB) and diploid wheat (Aegilops tauschii) (DD) approximately 8000 to 10,000 years ago [42, 44]. It was proposed that genome duplication and evolution in tetraploid and hexaploid wheat may lead to subfunctionalization or neofunctionalization of gene homoeologs as has been reported for other polyploid plant species [45–47]. However, neofunctionalization occurs less often in polyploid wheat species due to their relatively short evolutionary time spans . Moreover, differences among wheat subgenomes are often manifested as differential expression of homoeologous genes (i.e. different homoeologs of the same gene). Consistent with this notion, we recently observed distinct expression patterns of HYD1/2 and LCYe genes and homoeologs in developing tetraploid and hexaploid wheat grains (whole grains were analyzed in this study) . Particularly, HYD-B1 homoeolog expression strongly resembles those of embryo-specific genes . These data raised the possibility that carotenoid metabolic gene homoeologs could be differentially regulated and may have specific functions in different sections of wheat grains.
The overall objective of this work is to obtain a comprehensive understanding of metabolic genes and homoeologs controlling carotenoid, particularly β-carotene, accumulation in wheat. Since CCD1/4 had not been isolated and functionally characterized in wheat, we first cloned CCD1 and CCD4 homoeologs from wheat and determined the in vitro enzyme activities of their encoded proteins. We then analyzed and compared carotenoid content as well as expression of carotenoid metabolic gene homoeologs (including PSY1, LCYe, HYD1/2 and CCD1/4) in vegetative tissues and three sections of developing tetraploid and hexaploid wheat grains.
Plant growth and tissue collection
Seeds of tetraploid wheat var. Kronos and hexaploid wheat breeding line UC1041 were germinated on pre-wetted filter paper and transferred, at 7 d post-germination, to either vermiculite followed by growing in a temperature-controlled growth chamber (for collection of vegetative tissues), or soil followed by growing in a temperature-controlled greenhouse (for collection of grains). Leaf, stem and root tissues of wheat seedlings were harvested after growing in the chamber for three weeks; the collected tissues were immediately frozen in liquid nitrogen. The greenhouse-grown wheat plants were tagged when first anthers became visible from the middle florets of an ear. Developing grains were harvested according to the six defined developmental stages, including watery/1, early milk/2, late milk/3, soft dough/4, hard dough/5 and ripening/6, as previously described . Grains of developmental stages 3–5 were dissected, using forceps and scalpels, into pericarp, endosperm and embryo sections. Grains at stages 1, 2 and 6 are not suitable for dissection as stages 1 and 2 grains are quite small and watery while stage 6 grains are very dry. The dissected grain tissues were immediately frozen in liquid nitrogen. Three biological replicates, each containing pooled sections from multiple grains harvested from several plants, were used for gene expression and metabolite analyses. The frozen vegetative tissues and grain sections were ground into fine powder in liquid nitrogen using mortar and pestle and stored at −80 °C until further analysis.
RNA extraction and cloning of wheat CCD1 and CCD4 homoeologs
Total RNA was extracted from wheat tissues using TRI reagent (Invitrogen, Carlsbad, CA). Quantity and quality of the RNA samples were determined using a Nanodrop® spectrophotometer, according to absorption at 260 nm (RNA quantity) as well as the ratios of A260/A280 and A260/A230 (RNA quality). RNA integrity was also evaluated by agarose gel electrophoresis. First strand cDNA synthesis was performed using the BioRad iScript cDNA synthesis kit with mixed random hexamers and oligo(dT)20 primers (Hercules, CA).
The coding sequences of wheat CCD1 or CCD4 homoeologs are highly similar at the 5′ and 3′ ends. Therefore, the same set of primers, except for CCD-B1, was used for amplification of wheat CCD1 or CCD4 homoeologs and cloning into the pENTR/D-TOPO vector (Invitrogen). Plasmids extracted from multiple colonies were sequenced in each cloning experiment to identify different homoeologs of CCD1 or CCD4. Upon sequence confirmation, the CCD1/4 gene homoeologs in pENTR/D-TOPO were then recombined into pDEST17 (Invitrogen) for expression as His-tagged proteins in E. coli. Arabidopsis thaliana CCD1 (AtCCD1) and β-glucuronidase (GUS) were also cloned into pDEST17 and used as positive and negative controls for carotenoid cleavage activities, respectively. To improve the solubility of the recombinant CCD4 proteins, His-tagged wheat CCD4 homoeologs were also subcloned into the pMAL-c2x vector (New England BioLabs, Ipswich, MA) for tagging of the maltose binding protein (MBP). The His-tagged wheat CCD-A1 was cloned into pMAL-c2x in parallel and used as a control for comparing the cleavage activities between His-tagged and MBP-His-tagged CCD proteins.
When analyzed by TargetP , wheat CCD4 homoeologs were predicted to contain N-terminal plastid transit peptides of various lengths (50 aa for CCD-A4, 65 aa for CCD-B4, and 49 aa for CCD-D4). Since it was suggested that removal of subcellular targeting sequences may improve recombinant protein expression in E. coli [49, 50], truncated CCD4 homoeologs (i.e. without the transit peptide-encoding DNA sequences) were also cloned into pENTR/D-TOPO and then pDEST17 for protein expression in E. coli. Primers used for cloning wheat CCD1 and CCD4 homoeologs are listed in Additional file 1: Table S1.
Purification of recombinant proteins and CCD enzyme assays
Recombinant plasmids carrying AtCCD1, GUS, wheat CCD1 or CCD4 (full length or truncated) homoeologs were transformed into chemically competent cells of E. coli strain Rosetta (DE3)pLysS. The E. coli cells were grown at 37 °C until OD600 reached 0.6-0.8. Isopropyl β-D-1-thiogalactopyranoside (IPTG) was added to the bacterial cell culture to a final concentration of 0.4 mM. The cells continued to grow at 16 °C for 20 h and were collected by centrifugation at 3000 x g for 20 min. Purification of His-tagged proteins was carried out using Ni-NTA agarose beads (Qiagen, Valencia, CA). MBP-His-tagged proteins were purified using amylose resins (New England BioLabs). Induction and purification of the recombinant proteins were examined by SDS-PAGE.
β-apo-8′-carotenal (Sigma Aldrich, St. Louis, MO), lutein (Chromadex, Irvine, CA) and zeaxanthin (purified from E. coli cells expressing pAC-ZEAX) were dissolved in ethanol containing 1.5 % n-octyl-β-D-glucopyranoside (Sigma Aldrich). β-carotene (Sigma Aldrich) was solubilized as previously described with slight modifications . Briefly, β-carotene was dissolved in 20 μl of chloroform with 2 % Tween-80. The mixture was dried under vacuum, resuspended in 20 μl of ethanol, and then added to the reaction mixture. For in vitro CCD enzyme assays, the 200 μl reaction mixture included 20 μg purified recombinant proteins, 100 mM Tris-HCl (pH 7.0), 0.5 mM FeSO4, 5 mM ascorbic acid, and the (apo)carotenoid substrate (120 μM β-apo-8′-carotenal, 40 μM β-carotene, 20 μM zeaxanthin or 100 μM lutein). The reaction was allowed to proceed at 30 °C for 2 h and stopped by extraction with ethyl acetate; 20 μl of the ethyl acetate extract was injected on a reverse phase HPLC. The enzyme assays were carried out in triplicate.
Recombinant plasmids carrying wheat CCD1 or CCD4 (full length) homoeologs as well as the GUS control were also transformed into E. coli JM109 (DE3) cells that express pAC-BETA (containing genes for β-carotene production). An overnight culture initiated from a single colony was used to inoculate 25-ml Luria Bertani (LB) with 1 % glucose and appropriate antibiotics. The bacterial culture was grown at room temperature with shaking until OD600 reached 0.6. IPTG was added to a final concentration of 0.02 mM and the cells continued to grow at room temperature for an additional 12 h. The growth media were then extracted with ether, dried under N2, and resuspended in ethyl acetate; 20 μl of the ethyl acetate resuspension was separated on HPLC. The medium extractions were repeated at least three times from independently transformed and grown bacterial cell cultures.
HPLC and MS analyses
The HPLC program for analysis of CCD enzyme assay products consisted of three solvents: (A) 20 % acetonitrile, (B) acetonitrile:H2O:triethylamine (900:99:1, v/v/v), and (C) ethyl acetate. The gradient elution included: 0–2 min, 100 % A; 2–5 min, 100-50 % A and 0–50 % B; 5–10 min, 50-0 % A and 50–100 % B; 10–13 min, 100-75 % B and 0–25 % C; 13–16 min, 75-30 % B and 25–70 % C; 16–19 min, 30-0 % B and 70–100 % C; 19–21 min, 100 % C; 21–22 min, 0–100 % A and 100-0 % C. Carotenoid extraction from wheat grain tissues as well as HPLC separation and quantification were carried out as previously described . The HPLC program for analysis of bacterial growth media included (A) 50 % acetonitrile, (B) acetonitrile:H2O:triethylamine (900:99:1, v/v/v), and (C) ethyl acetate, with the following gradient: 0–10 min, 100-0 % A and 0–100 % B; 10–13 min, 100-75 % B and 0–25 % C; 13–16 min, 75-30 % B and 25–70 % C; 16–19 min, 30-0 % B and 70–100 % C; 19–21 min, 100 % C; 21–22 min, 0–100 % B and 100-0 % C. The flow rate was maintained at 1 ml min−1 for all of the HPLC analysis.
Mass spectrometry (MS) analysis was performed on a Thermo Electron LTQ-Orbitrap Hybrid mass spectrometer (Thermo Scientific, Waltham, MA). Product peaks from the CCD enzyme assays were collected from the HPLC runs and used for injection onto the mass spectrometer. An isocratic flow of 50 % (A) H2O and 50 % (B) acetonitrile was maintained at 0.2 ml min−1 for MS analysis. The mass spectra were acquired by electrospray ionization (ESI) in the positive mode with a mass range of m/z 50–1000 Da.
Representative CCD and nine-cis-epoxycarotenoid dioxygenase (NCED) proteins from different plant species were selected for multiple sequence alignment using Multiple Sequence Comparison by Log-Expectation (MUSCLE) . The aligned sequences were then used for construction of a neighbor-joining (NJ) tree by MEGA5 with the pairwise deletion method . Robustness of the NJ tree was tested with 1000 rounds of bootstrapping. The rice (Os), sorghum (Sb) and maize (Zm) CCD and NCED sequences were obtained from . The GenBank (non-A. thaliana sequences) and AGI (A. thaliana sequences) accession numbers are: AtCCD1, At3g63520; AtCCD4, At4g19170; AtCCD7, At2g44990; AtCCD8, At4g32810; AtNCED2, At4g18350; AtNCED3, At3g14440; AtNCED5, At1g30100; AtNCED6, At3g24220; AtNCED9, At1g78390; CaCCD1, DQ157170; CcCCD1, DQ157166; CisCCD1, AB219165; ClCCD1, AB219168; CmCCD1, DQ269467; CsCCD1a, AJ132927; CsCCD1b, EU523661; CuCCD1, AB219164; DcCCD1, DQ192203; FaCCD1, ACA13522; OfCCD1, BAJ05401; PhCCD1, AY576003; RdCCD1, ABY47994; SlCCD1A, AY576001; SlCCD1B, AY576002; TaCCD-A1, KU975445; TaCCD-B1, KU975446; TaCCD-D1, KU975447; VvCCD1, AY856353; CsCCD4a, EU523662; CsCCD4b, EU523663; CsCCD4c, JN131499; CuCCD4, BAO18774; CxmCCD4a, ABY60885; CxmCCD4b, BAF36656; MdCCD4, ABY47995; OfCCD4, ABY60887; PpCCD4, JX309999; RdCCD4, ABY60886; StCCD4, XP_006359966; TaCCD-A4, KU975448; TaCCD-B4, KU975449; TaCCD-D4, KU975450; BoLCD, AJ489277; CsZCD, AJ489276. At, Arabidopsis thaliana; Bo, Bixa orellana; Ca, Coffea arabica; Cc, Coffea canephora; Cis, Citrus sinensis; Cl, Citrus limon; Cm, Cucumis melo; Cs, Crocus sativus; Cu, Citrus unshiu; Cxm, Chrysanthemum x morifolium; Dc, Daucas carota; Fa, Fragaria ananassa; Md, Malus x domestica; Of, Osmanthus fragans; Os, Oryza sativa; Ph, Petunia hybrida; Pp, Prunus persica; Rd, Rosa x damascene; Sb, Sorghum bicolor; Sl, Solanum lycopersicum; St, Solanum tuberosum; Ta, Triticum aestivum; Vv, Vitis vinifera; Zm, Zea mays. LCD, Lycopene Cleavage Dioxygenase; ZCD, Zeaxanthin Cleavage Dioxygenase.
Real-time qPCR analysis
Total RNA was treated with RNase-free DNase I (Fermentas, Glen Burnie, MD). cDNA synthesis was carried out using 5 μg total RNA, random hexamers, and the iScript cDNA synthesis kit (BioRad). Primers specific for wheat HYD1 and HYD2 homoeologs as well as the reference genes were reported previously . Primers specific for PSY1, LCYe, CCD1 and CCD4 homoeologs were designed and verified using nullisomic-tetrasomic and ditelosomic lines of hexaploid wheat var. Chinese Spring (Additional file 2: Table S2; Additional file 3: Figure S1). The qPCR products were cloned and sequenced to verify amplification of the target gene homoeologs.
Real-time qPCR analysis was carried out as previously described using the iTaq SYBR® Green Supermix . qPCR amplification efficiency for the target gene homoeologs was in the range of 90 % to 121 %. The relative standard curve method was used for quantification of the transcripts . cDNAs synthesized from whole grains were used for construction of the standard curves in order to allow comparison of homoeologs expressed in different grain sections. Normalization of gene expression was carried out using the geometric mean of two reference genes, Ta2291 and Ta54227, which are stably expressed in different wheat tissues .
One way Analysis of Variation (ANOVA) followed by unpaired, two-tailed t-test (carotenoid content as well as gene expression in tetraploid wheat) or Tukey’s test (gene expression in hexaploid wheat) were performed. All statistical analysis was carried out using JMP (SAS Institute, Cary, NC).
Cloning and biochemical characterization of wheat CCD1 and CCD4 homoeologs
Wheat CCD1 and CCD4 homoeologs were identified by searching TIGR gene indices as well as the genomic sequences of Ae. tauschii and hexaploid wheat var. Chinese Spring using the Arabidopsis CCD1 and CCD4 sequences as queries. CCD1 and CCD4 homoeologs were assigned to different wheat subgenomes based on the wheat rice synteny and the subgenome assignment was verified using nullisomic-tetrasomic and ditelosomic lines of Chinese Spring (Additional file 3: Figure S1). When compared with CCDs and NCEDs that were previously identified from other plant species, wheat CCD1 and CCD4 homoeologs fell within the respective CCD1 and CCD4 clades with strong support (Fig. 2). Moreover, wheat CCD1 and CCD4 homoeologs clustered closely with their monocot relatives (Fig. 2).
To determine the catalytic activities of wheat CCD1 and CCD4 homoeologs towards (apo)carotenoid substrates, purified recombinant proteins of CCD1/4 were used for in vitro enzyme assays (Additional file 4: Figure S2). The three wheat CCD1 homoeologs cleaved β-apo-8′-carotenal (C30 apocarotenoid) at the C9-C10 position to form β-ionone (m/z 193.17, [M + H]+) and β-apo-10,8′-carotendial (C17 dialdehyde; m/z 257.17, [M + H]+) (Fig. 3; Additional files 5: Figure S3 and Additional file 6: Figure S4). In addition to β-apo-8′-carotenal, wheat CCD1 homoeologs also exhibited significant cleavage activities towards lutein, while relatively small amounts of products were generated when β-carotene and zeaxanthin were provided as substrates (Fig. 3; Additional files 5: Figure S3 and Additional file 6: Figure S4). It should be noted that β-carotene was only partially soluble in the reaction buffer, which may have limited the accessibility of this substrate to the enzymes.
The C14 dialdehyde β-apo-10,10′-carotendial (m/z 217, [M + H]+) was produced in reactions with β-carotene, zeaxanthin or lutein, suggesting symmetrical cleavage of these carotenoids at C9-C10/C9′-C10′ double bonds by wheat CCD1 homoeologs (Fig. 3; Additional files 5: Figure S3 and Additional file 6: Figure S4). Symmetrical cleavage at C9-C10/C9′-C10′ of the carotenoid substrates was also supported by the release of β-ionone (from β-carotene; m/z 193.17, [M + H]+), 3-OH-β-ionone (from zeaxanthin; major fragment m/z 191.08, [M + H]+), and 3′-OH-α-ionone + 3-OH-β-ionone (from lutein; m/z 209.17, [M + H]+) in these reactions (Fig. 3; Additional files 5: Figure S3 and Additional file 6: Figure S4). 3-OH-β-ionone and 3′-OH-α-ionone co-eluted at 6.97 min in the reaction with lutein as substrate, and showed different fragmentation patterns in MS analysis than 3-OH-β-ionone obtained from the zeaxanthin cleavage reaction (Fig. 3) [31, 55]. CCD1 activities towards β-carotene were further verified by analyzing the growth media of E. coli cells co-transformed with pAC-BETA (a plasmid that contains β-carotene biosynthetic genes) and CCD1 homoeologs. β-ionone was also found in the growth media of β-carotene producing E. coli cells that express wheat CCD1 homoeologs (Fig. 4; Additional file 7: Figure S5).
Unlike the His-tagged CCD1 proteins that partially partitioned to the soluble protein fraction, His-tagged CCD4 proteins were largely insoluble (Additional file 4: Figure S2). CCD4 homoeologs were then expressed as MBP-His-tagged proteins and showed improved protein solubility (Additional file 4: Figure S2). Overall, the purified MBP-His-CCD4 proteins did not exhibit any cleavage activity towards β-apo-8′-carotenal, β-carotene, zeaxanthin and lutein in enzyme assays (Fig. 3; Additional files 5: Figure S3 and Additional file 6: Figure S4). In addition, β-ionone was not detected in growth media of β-carotene producing E. coli cells co-expressing wheat CCD4 homoeologs (Fig. 4; Additional file 7: Figure S5). To demine whether the MBP tag may affect protein activities, MBP-His-CCD-A1 was also expressed in E. coli and showed improved protein solubility as compared to His-CCD-A1 (Additional file 4: Figure S2). Comparable cleavage activities of the (apo)carotenoid substrates were demonstrated by MBP-His-CCD-A1 and His-CCD-A1, indicating that the MBP tag does not affect the activity of CCD proteins (Figs. 3 and 4). Upon removal of N-terminal plastid transit peptides, solubility of His-truncated CCD-4A/B/D proteins was moderately improved (Additional file 4: Figure S2). However, His-truncated CCD-4A/B/D proteins were not active towards β-apo-8′-carotenal, β-carotene, zeaxanthin and lutein (Additional file 8: Figure S6).
Carotenoid metabolite profiles in different sections of developing wheat grains
Our previous analysis of whole grains revealed progressively decreased carotenoid accumulation during the 6 defined stages of tetraploid and hexaploid wheat grain development . To further investigate spatial carotenoid accumulation in developing grains, total carotenoids from endosperm, embryo and pericarp sections of tetraploid (var. Kronos) and hexaploid (breeding line UC1041) wheat grains, at stages 3–5 (representing late milk, soft and hard dough stages), were analyzed and compared (Tables 1 and 2).
The pericarp of wheat grains contain β-carotene, neoxanthin, violaxanthin and lutein, resembling the carotenoid composition of photosynthetic tissues (Tables 1 and 2) . In the pericarp tissue, similar individual and total carotenoid levels were observed between stages 3 and 4 of tetraploid and hexaploid wheat; however, there was a consistent drop in carotenoid content at stage 5, which was significant for β-carotene, neoxanthin and total carotenoids (Tables 1 and 2). On the other hand, an increasing trend in the ratio of β,ε-/β,β-carotene branch (β,ε/β,β) carotenoids was evident for stages 3–5 of grain pericarp (Tables 1 and 2), indicating a shift in carotenoid accumulation favoring the α-carotene/lutein branch relative to the β-carotene/zeaxanthin branch (Fig. 1).
In contrast to pericarp, lutein and violaxanthin were the only detectable carotenoids in the endosperm of tetraploid and hexaploid wheat grains (Tables 1 and 2). The endosperm of tetraploid wheat consistently accumulated higher concentrations of lutein and violaxanthin than that of hexaploid wheat at each grain developmental stage. As grain endosperm matures, violaxanthin decreased more rapidly than lutein and led to a one-fold increase in the β,ε/β,β ratio at stage 5 of tetraploid and hexaploid wheat grains (Tables 1 and 2).
In addition to the carotenoids present in pericarp and endosperm tissues, antheraxanthin and zeaxanthin were also found in embryo (Tables 1 and 2). Particularly, antheraxanthin accounted for about 20 % of total carotenoids in the embryo tissue. As with pericarp and endosperm, total carotenoids in embryo decreased towards the late stage of grain development, contributed mainly by reduced lutein. While the β,ε/β,β ratios in maturing embryo decreased in tetraploid wheat, they remained mostly constant in hexaploid wheat (Tables 1 and 2).
Expression of carotenoid metabolic gene homoeologs in different sections of developing wheat grains
To assess the spatial contribution of carotenoid metabolic gene homoeologs to grain carotenoid accumulation, expression of PSY1, LCYe, HYD1/2 and CCD1/4 homoeologs in pericarp, endosperm and embryo tissues of tetraploid and hexaploid wheat grains was analyzed using real-time qPCR (Figs. 5 and 6).
In the pericarp tissue, PSY-A1 and PSY-B1 showed comparable and similarly decreased expression during tetraploid wheat grain development, except in stage 5 where the PSY-B1 transcript level was significantly higher than that of PSY-A1 (Fig. 5). The preponderance of the B-genome PSY1 homoeolog was more evident in hexaploid wheat; PSY-B1 was the major PSY1 homoeolog expressed in stage 3 and the only PSY1 transcript detected in stages 4 and 5 (Figs. 5 and 6). While LCYe, HYD1 and HYD2 homoeologs exhibited similar transcript accumulation through grain development in tetraploid wheat, HYD1 homoeologs showed increased expression at stage 5 of grain development in hexaploid wheat. In addition, LCYe-A was the only LCYe homoeolog expressed in the pericarp of hexaploid wheat. Of the CCD1/4 homoeologs, CCD-B4 expression in pericarp increased gradually towards stage 5 of tetraploid and hexaploid wheat grains. The other CCD1/4 homoeologs showed similar or slightly decreased expression (CCD-A1) during tetraploid and hexaploid wheat grain development (Figs. 5 and 6).
In the endosperm tissue, transcripts of all three PSY1 homoeologs were under the detection limit of real-time qPCR in hexaploid wheat grains, contrasting to the high PSY1 homoeolog expression found in tetraploid wheat where decreased PSY-A1 and PSY-B1 transcript levels were observed in stages 3–5 (Figs. 5 and 6). Interestingly, LCYe-A was the only LCYe homoeolog expressed in the endosperm of both tetraploid and hexaploid wheat grains. Of the HYD1/2 genes, HYD1 homoeolog expression was not detectable, while the HYD2 homoeologs showed decreased expression in tetraploid and hexaploid wheat (Figs. 5 and 6). Only CCD-A1 and CCD-B4 expression was observed in tetraploid wheat grain endosperm. In contrast, all three CCD1 homoeologs were expressed in the hexaploid wheat grain endosperm, with CCD-A1 and CCD-D1 showing increased expression during grain development; on the other hand, none of the CCD4 homoeologs had detectable transcript accumulation in hexaploid wheat grain endosperm (Figs. 5 and 6).
In the embryo tissue, PSY-B1 and PSY-D1 were the most abundant PSY1 homoeolog in tetraploid and hexaploid wheat, respectively (Figs. 5 and 6). LCYe-B in tetraploid wheat as well as LCYe-A and LCYe-B in hexaploid wheat showed increased expression during grain maturation, particularly at stage 5. While HYD2 homoeologs had mostly consistent expression in the three grain developmental stages analyzed, enhanced HYD1 homoeolog expression was observed at the late stage of tetraploid and hexaploid wheat grain development (Figs. 5 and 6). Interestingly, transcripts of CCD1, but not CCD4, homoeologs were detectable in the embryo tissue of tetraploid and hexaploid wheat grains (Figs. 5 and 6).
Expression of carotenoid metabolic gene homoeologs in vegetative tissues
To determine the genetic factors controlling carotenoid accumulation in vegetative tissues of tetraploid and hexaploid wheat, expression of carotenoid metabolic gene homoeologs in three vegetative tissues, including leaf, stem and root, was analyzed by real-time qPCR (Figs. 7 and 8).
In leaves, HYD1/2 and CCD1/4 homoeologs showed similar expression in both tetraploid and hexaploid wheat. On the other hand, PSY-A1 and LCYe-A transcripts were relatively more abundant than their corresponding homoeologs in tetraploid wheat, whereas PSY-D1 and LCYe-B transcripts were relatively more abundant than their corresponding homoeologs in hexaploid wheat (Figs. 7 and 8).
In stems of tetraploid and hexaploid wheat, the relative expression patterns among PSY1, LCYe, HYD1 and CCD1/4 homoeologs generally resembled those in leaves. However, expression levels of PSY1 and CCD1/4 homoeologs were much reduced in stems as compared to leaves (Figs. 7 and 8). Unlike the comparable transcript accumulation of HYD2 homoeologs in leaves, HYD-A2 appeared to be the major HYD2 homoeolog expressed in stems of tetraploid and hexaploid wheat (Figs. 7 and 8).
In roots of tetraploid and hexaploid wheat, PSY1 and CCD4 homoeologs were under the limit of detection, and LCYe homoeologs were absent or very lowly expressed (Figs. 7 and 8). In tetraploid wheat roots, HYD-A1 and HYD-A2 were clearly the major HYD1 and HYD2 homoeologs, respectively. In contrast, comparable expression was observed for the three homoeologs of HYD2, while HYD-D1 had only baseline expression in hexaploid wheat roots (Figs. 7 and 8). CCD-A1 and CCD-B1 homoeolog expression did not differ significantly for both tetraploid and hexaploid wheat roots, whereas CCD-D1 showed two-fold higher expression than the other two CCD1 homoeologs in roots of hexaploid wheat seedlings (Figs. 7 and 8).
Lutein, but not β-carotene, accumulates in the endosperm (flour) of tetraploid and hexaploid wheat grains at varied levels (Tables 1 and 2). While re-directing the carbon flux from lutein to β-carotene formation may be accomplished in wheat grains by blocking the LCYe catalyzed reaction as demonstrated in maize grains and potato tubers [56, 57], it is crucial to understand whether the β-carotene produced could accumulate or will be turned over via hydroxylation by HYDs or cleavage by CCDs (Fig. 1). Our in vitro enzyme assay results showed that wheat CCD1 homoeologs could cleave β-carotene symmetrically, suggesting that CCD1 homoeologs can potentially contribute to the degradation of this provitamin A molecule in grains modified with increased β-carotene production (Figs. 3 and 4; Additional files 5: Figure S3 and Additional file 6: Figure S4). Since lutein can also serve as substrate to CCD1 homoeologs in vitro and it naturally accumulates in wheat grain endosperm, it is possible that the rate of lutein biosynthesis is higher than that of potential degradation by CCD1 homoeologs in this tissue. Alternatively, there could be unknown mechanisms that protect lutein from degradation by CCD1 homoeologs in wheat grain endosperm.
In contrast to CCD1 homoeologs, wheat CCD4 homoeologs did not act on β-carotene, lutein and zeaxanthin (Figs. 3 and 4; Additional files 5: Figure S3, Additional file 6: Figure S4, Additional file 7: Figure S5 and Additional file 8: Figures S6). It remains to be determined whether they could function towards other carotenoid molecules, such as violaxanthin, that were not used for in vitro testing. AtCCD4 (Arabidopsis CCD4) was previously found in the plastid-localized plastoglobules; in particular, it interacts with a zinc finger protein VAR3 in a protein complex required for normal chloroplast and palisade cell development [35, 58]. It will be interesting to determine whether wheat CCD4 homoeologs are also bound to plastoglobules and cooperates with VAR3-like proteins. Further reverse genetic analysis of wheat CCD4 homoeologs will help elucidate their enzymatic functions in planta and the role of their cleavage products in plastid, cell and tissue development.
It was previously shown that maize genotypes with low levels of CHY (i.e. HYD) and LCYe transcripts correlate with high levels of β-carotene accumulation in grains [57, 59], which suggests that carotenoid metabolic gene expression may be indicative of their functions in plants. We therefore examined the expression profiles of carotenoid metabolic gene homoeologs to determine the specific gene homoeologs that control β-carotene accumulation in different sections of developing wheat grains. Overall, the gene expression data indicated that homoeologous carotenoid metabolic genes are regulated at the transcriptional level in sections of wheat grains (Figs. 5, 6, 7 and 8). Additionally, tissue specific expression was more prominent than genome specific expression for the carotenoid metabolic gene homoeologs in wheat (Figs. 5, 6, 7 and 8). Recently, a genome-wide study examined homoeolog expression and genome interplay among three cell types (including endosperm, aleurone layer and transfer cells) of hexaploid wheat grains at 10-, 20- and 30-day after anthesis . Similar to that observed in our gene expression analysis, this study also reported differential expression of homoeologous genes without an overall dominance of a particular subgenome . However, embryo and pericarp tissues were not examined and expression of carotenoid metabolic gene homoeologs were not analyzed in this report .
In developing embryo of tetraploid and hexaploid wheat grains, LCYe-A and LCYe-B homoeologs showed increased expression contrasting to decreased lutein accumulation , suggesting that mechanisms other than transcriptional control could also be involved in the regulation of lutein metabolism in wheat grains. Accompanying largely reduced individual and total carotenoids in pericarp at the late stage of wheat grain development, CCD-B4 expression increased sharply, raising the possibility that it may be responsible for carotenoid degradation in pericarp during grain dehydration. This observation is consistent with a recent study in Arabidopsis where AtCCD4 expression rose during seed drying (comparable to wheat grain stages 4–5) and was shown genetically as a major contributor of β-carotene degradation during Arabidopsis seed desiccation . However, since whole Arabidopsis seeds were used in this study, the spatial location of AtCCD4 expression and β-carotene accumulation within Arabidopsis seeds were not determined.
Our gene expression results also provided important insights for future attempts to increase β-carotene content in wheat grain endosperm. In the endosperm of tetraploid and hexaploid wheat grains, HYD1 homoeologs were not expressed; CCD1 and CCD4 homoeologs were either non-detectable or had low level expression, thus may not impact carotenoid accumulation in this tissue. On the other hand, there were significant transcript accumulation of LCYe-A (but not LCYe-B or LCYe-D) and HYD2 homoeologs, suggesting that downregulation or loss-of-function of these homoeologs may be sufficient to result in more β-carotene in the endosperm of tetraploid and hexaploid wheat grains. Since the total carotenoid content in hexaploid wheat grain endosperm is low, more carbons will also need to the directed to the carotenoid biosynthetic pathway, such as introducing PSY1 genes with higher expression or activity, to increase β-carotene levels in hexaploid wheat grain endosperm.
Considering the critical roles that carotenoids play in light harvesting and photoprotection, it is important to ensure that downregulation of specific carotenoid metabolic gene homoeologs in grain endosperm will not compromise carotenoid biosynthesis and function in the photosynthetic tissue. LCYe, HYD1/2 and CCD1/4 homoeologs showed comparable expression in leaves (Figs. 7 and 8), suggesting that loss/reduction of activities of specific carotenoid metabolic gene homoeologs in the photosynthetic tissue could be compensated for by the overlapping activities of homologous (e.g. HYD1 and HYD2) and/or homoeologous (e.g. different homoeologs of HYD2) genes.
Taken together, the CCD1/4 enzyme activity and the spatial gene expression analyses suggested that reduced expression of LCYe-A and one or more of the loci, including HYD-A2, HYD-B2, and CCD1 homoeologs (only for hexaploid wheat), could lead to β-carotene enrichment in the endosperm of wheat grains without compromising carotenoid metabolism in leaves. Identification of the specific carotenoid metabolic gene homoeologs controlling β-carotene accumulation will facilitate efficient and effective provitamin A biofortification of wheat grains through plant breeding and genome editing technologies.
ANOVA, analysis of variation; CCD, carotenoid cleavage dioxygenase; ESI, electrospray ionization; GUS, β-glucuronidase; HYD, carotenoid β-ring hydroxylase; IPTG, isopropyl β-D-1-thiogalactopyranoside; LCD, lycopene cleavage dioxygenase; LCYe, lycopene ε-cyclase; MBP, maltose binding protein; MS, mass spectrometry; MUSCLE, multiple sequence comparison by log-expectation; NCED, nine-cis-epoxycarotenoid dioxygenase; NJ, neighbor-joining; PSY, phytoene synthase; ZCD, zeaxanthin cleavage dioxygenase
Panfili G, Fratianni A, Irano M. Improved normal-phase high-performance liquid chromatography procedure for the determination of carotenoids in cereals. J Agric Food Chem. 2004;52:6373–7.
Hentschel V, Kranl K, Hollmann J, Lindhauer M, Bohm V, Bitsch R. Spectrophotometric determination of yellow pigment content and evaluation of carotenoids by high-performance liquid chromatography in durum wheat grain. J Agric Food Chem. 2002;50:6663–8.
Calucci L, Capocchi A, Galleschi L, Ghiringhelli S, Pinzino C, Saviozzi F, Zandomeneghi M. Antioxidants, free radicals, storage proteins, puroindolines, and proteolytic activities in bread wheat (Triticum aestivum) seeds during accelerated aging. J Agric Food Chem. 2004;52:4274–81.
Britton G. Functions of intact carotenoids. In: Britton G, Liaaen-Jensen S, Pfander H, editors. Carotenoids Volume 4: Natural Functions. Basel: Birkhäuser Verlag; 2008.
Britton G. Vitamin A and vitamin A deficiency. In: Britton G, Liaaen-Jensen S, Pfander H, editors. Carotenoids Volume 5: Nutrition and Health. Basel: Birkhäuser Verlag; 2009.
Bendich A. From 1989 to 2001: what have we learned about the “biological actions of beta-carotene”? J Nutr. 2004;134:225S–30.
Mora JR, Iwata M, von Andrian UH. Vitamin effects on the immune system: vitamins A and D take centre stage. Nat Rev Immunol. 2008;8:685–98.
Howitt C, Cavanagh C, Bowerman A, Cazzonelli C, Rampling L, Mimica J, Pogson B. Alternative splicing, activation of cryptic exons and amino acid substitutions in carotenoid biosynthetic genes are associated with lutein accumulation in wheat endosperm. Funct Integr Genomics. 2009;9:363–76.
Zhang W, Dubcovsky J. Association between allelic variation at the Phytoene synthase 1 gene and yellow pigment content in the wheat grain. Theor Appl Genet. 2008;116:635–45.
Qin X, Zhang W, Dubcovsky J, Tian L. Cloning and comparative analysis of carotenoid β-hydroxylase genes provides new insights into carotenoid metabolism in tetraploid (Triticum turgidum ssp. durum) and hexaploid (Triticum aestivum) wheat grains. Plant Mol Biol. 2012;80:631–46.
Alder A, Jamil M, Marzorati M, Bruno M, Vermathen M, Bigler P, Ghisla S, Bouwmeester H, Beyer P, Al-Babili S. The path from β-carotene to carlactone, a strigolactone-like plant hormone. Science. 2012;335:1348–51.
Auldridge ME, Block A, Vogel JT, Dabney-Smith C, Mila I, Bouzayen M, Magallanes-Lundback M, DellaPenna D, McCarty DR, Klee HJ. Characterization of three members of the Arabidopsis carotenoid cleavage dioxygenase family demonstrates the divergent roles of this multifunctional enzyme family. Plant J. 2006;45:982–93.
Auldridge ME, McCarty DR, Klee HJ. Plant carotenoid cleavage oxygenases and their apocarotenoid products. Curr Opin Plant Biol. 2006;9:315–21.
Mathieu S, Terrier N, Procureur J, Bigey F, Günata Z. A carotenoid cleavage doxygenase from Vitis vinifera L.: functional characterization and expression during grape berry development in relation to C13-norisoprenoid accumulation. J Exp Bot. 2005;56:2721–31.
Vogel J, Tan B, McCarty D, Klee H. The carotenoid cleavage dioxygenase 1 enzyme has broad substrate specificity, cleaving multiple carotenoids at two different bond positions. J Biol Chem. 2008;283:11364–73.
Ilg A, Beyer P, Al-Babili S. Characterization of the rice carotenoid cleavage dioxygenase 1 reveals a novel route for geranial biosynthesis. FEBS J. 2009;276:736–47.
Gang D. Evolution of flavors and scents. Annu Rev Plant Biol. 2005;56:301–25.
Huang F, Molnár P, Schwab W. Cloning and functional characterization of carotenoid cleavage dioxygenase 4 genes. J Exp Bot. 2009;60:3011–22.
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–201.
Ohmiya A. Carotenoid cleavage dioxygenases and their apocarotenoid products in plants. Plant Biotechnol. 2009;26:351–8.
Gonzalez-Jorge S, Ha S, Magallanes-Lundback M, Gilliland L, Zhou A, Lipka A, Nguyen Y, Angelovici R, Lin H, Cepela J, et al. Carotenoid cleavage dioxygenase4 is a negative regulator of β-carotene content in Arabidopsis seeds. Plant Cell. 2013;25:4812–26.
Rodrigo M, Alquézar B, Alós E, Medina V, Carmona L, Bruno M, Al-Babili S, Zacarías L. A novel carotenoid cleavage activity involved in the biosynthesis of citrus fruit-specific apocarotenoid pigments. J Exp Bot. 2013;64:4461–78.
Bouvier F, Suire C, Mutterer J, Camara B. Oxidative remodeling of chromoplast carotenoids: identification of the carotenoid dioxygenase CsCCD and CsZCD genes involved in crocus secondary metabolite biogenesis. Plant Cell. 2003;15:47–62.
Rubio A, Rambla JL, Santaella M, Gómez MD, Orzaez D, Granell A, Gómez-Gómez L. Cytosolic and plastoglobule-targeted carotenoid dioxygenases from Crocus sativus are both involved in β-ionone release. J Biol Chem. 2008;283:24816–25.
Simkin AJ, Underwood BA, Auldridge M, Loucas HM, Shibuya K, Schmelz E, Clark DG, Klee HJ. Circadian regulation of the PhCCD1 carotenoid cleavage dioxygenase controls emission of β-ionone, a fragrance volatile of petunia flowers. Plant Physiol. 2004;136:3504–14.
Simkin AJ, Schwartz SH, Auldridge M, Taylor M, Klee HJ. The tomato carotenoid cleavage dioxygenase 1 genes contribute to the formation of the flavor volatiles β-ionone, pseudoionone, and geranylacetone. Plant J. 2004;40:882–92.
Kato M, Matsumoto H, Ikoma Y, Okuda H, Yano M. The role of carotenoid cleavage dioxygenases in the regulation of carotenoid profiles during maturation in citrus fruit. J Exp Bot. 2006;57:2153–64.
Ibdah M, Azulay Y, Portnoy V, Wasserman B, Bar E, Meir A, Burger Y, Hirschberg J, Schaffer A, Katzir N, et al. Functional characterization of CmCCD1, a carotenoid cleavage dioxygenase from melon. Phytochemistry. 2006;67:1579–89.
Simkin AJ, Moreau H, Kuntz M, Pagny G, Lin C, Tanksley S, McCarthy J. An investigation of carotenoid biosynthesis in Coffea canephora and Coffea arabica. J Plant Physiol. 2008;165:1087–106.
Sun Z, Hans J, Walter M, Matusova R, Beekwilder J, Verstappen F, Ming Z, van Echtelt E, Strack D, Bisseling T, et al. Cloning and characterisation of a maize carotenoid cleavage dioxygenase (ZmCCD1) and its involvement in the biosynthesis of apocarotenoids with various roles in mutualistic and parasitic interactions. Planta. 2008;228:789–801.
García-Limones C, Schnäbele K, Blanco-Portales R, Luz Bellido M, Caballero JL, Schwab W, Muñoz-Blanco J. Functional characterization of FaCCD1: a carotenoid cleavage dioxygenase from strawberry involved in lutein degradation during fruit ripening. J Agric Food Chem. 2008;56:9277–85.
Baldermann S, Kato M, Kurosawa M, Kurobayashi Y, Fujita A, Fleischmann P, Watanabe N. Functional characterization of a carotenoid cleavage dioxygenase 1 and its relation to the carotenoid accumulation and volatile emission during the floral development of Osmanthus fragrans Lour. J Exp Bot. 2010;61:2967–77.
Vallabhaneni R, Bradbury LMT, Wurtzel ET. The carotenoid dioxygenase gene family in maize, sorghum, and rice. Arch Biochem Biophys. 2010;504:104–11.
Ilg A, Yu Q, Schaub P, Beyer P, Al-Babili S. Overexpression of the rice carotenoid cleavage dioxygenase 1 gene in Golden Rice endosperm suggests apocarotenoids as substrates in planta. Planta. 2010;232:691–9.
Ytterberg AJ, Peltier J-B, van Wijk KJ. Protein profiling of plastoglobules in chloroplasts and chromoplasts. A surprising site for differential accumulation of metabolic enzymes. Plant Physiol. 2006;140:984–97.
Ma G, Zhang L, Matsuta A, Matsutani K, Yamawaki K, Yahata M, Wahyudi A, Motohashi R, Kato M. Enzymatic formation of β-citraurin from β-cryptoxanthin and zeaxanthin by Carotenoid Cleavage Dioxygenase4 in the flavedo of citrus fruit. Plant Physiol. 2013;163:682–95.
Campbell R, Ducreux LJM, 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–64.
Brandi F, Bar E, Mourgues F, Horvath 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.
Falchi R, Vendramin E, Zanon L, Scalabrin S, Cipriani G, Verde I, Vizzotto G, Morgante M. Three distinct mutational mechanisms acting on a single gene underpin the origin of yellow flesh in peach. Plant J. 2013;76:175–87.
Lashbrooke JG, Young PR, Dockrall SJ, Vasanth K, Vivier MA. Functional characterisation of three members of the Vitis vinifera L. carotenoid cleavage dioxygenase gene family. BMC Plant Biol. 2013;13:156.
Bruno M, Beyer P, Al-Babili S. The potato carotenoid cleavage dioxygenase 4 catalyzes a single cleavage of β-ionone ring-containing carotenes and non-epoxidated xanthophylls. Arch Biochem Biophys. 2015;572:126–33.
Dubcovsky J, Dvorak J. Genome plasticity a key factor in the success of polyploid wheat under domestication. Science. 2007;316:1862–6.
Rubio-Moraga A, Rambla J, Fernández-de-Carmen A, Trapero-Mozos A, Ahrazem O, Orzáez D, Granell A, Gómez-Gómez L. New target carotenoids for CCD4 enzymes are revealed with the characterization of a novel stress-induced carotenoid cleavage dioxygenase gene from Crocus sativus. Plant Mol Biol. 2014;86:555–69.
Huang S, Sirikhachornkit A, Su X, Faris J, Gill B, Haselkorn R, Gornicki P. Genes encoding plastid acetyl-CoA carboxylase and 3-phosphoglycerate kinase of the Triticum/Aegilops complex and the evolutionary history of polyploid wheat. Proc Natl Acad Sci U S A. 2002;99:8133–8.
Akhunova A, Matniyazov R, Liang H, Akhunov E. Homoeolog-specific transcriptional bias in allopolyploid wheat. BMC Genomics. 2010;11:505.
Nomura T, Ishihara A, Yanagita R, Endo T, Iwamura H. Three genomes differentially contribute to the biosynthesis of benzoxazinones in hexaploid wheat. Proc Natl Acad Sci U S A. 2005;102:16490–5.
Bottley A, Koebner RMD. Variation for homoeologous gene silencing in hexaploid wheat. Plant J. 2008;56:297–302.
Emanuelsson O, Brunak S, von Heijne G, Nielsen H. Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. 2007;2:953–71.
Koch B, Sibbesen O, Halkier BA, Svendsen I, Lindberg Møller B. The primary sequence of cytochrome P450tyr, the multifunctional N-hydroxylase catalyzing the conversion of L-tyrosine to p-hydroxyphenylacetaldehyde oxime in the biosynthesis of the cyanogenic glucoside dhurrin in Sorghum bicolor (L.) Moench. Arch Biochem Biophys. 1995;323:177–86.
Hull AK, Celenza JL. Bacterial epression and purification of the Arabidopsis NADPH–cytochrome P450 reductase ATR2. Protein Express Purif. 2000;18:310–5.
Edgar R. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res. 2004;32:1792–7.
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–9.
Applied Biosystems. Guide to performing relative quantification of gene expression using real-time quantitative PCR. 2008.
Paolacci A, Tanzarella O, Porceddu E, Ciaffi M. Identification and validation of reference genes for quantitative RT-PCR normalization in wheat. BMC Mol Biol. 2009;10:11.
Mein JR, Dolnikowski GG, Ernst H, Russell RM, Wang X-D. Enzymatic formation of apo-carotenoids from the xanthophyll carotenoids lutein, zeaxanthin and β-cryptoxanthin by ferret carotene-9′,10′-monooxygenase. Arch Biochem Biophys. 2011;506:109–21.
Diretto G, Tavazza R, Welsch R, Pizzichini D, Mourgues F, Papacchioli V, Beyer P, Giuliano G. Metabolic engineering of potato tuber carotenoids through tuber-specific silencing of lycopene epsilon cyclase. BMC Plant Biol. 2006;6:13.
Harjes C, Rocheford T, Bai L, Brutnell T, Kandianis C, Sowinski S, Stapleton A, Vallabhaneni R, Williams M, Wurtzel E, et al. Natural genetic variation in lycopene epsilon cyclase tapped for maize biofortification. Science. 2008;319:330–3.
Næsted H, Holm A, Jenkins T, Nielsen HB, Harris CA, Beale MH, Andersen M, Mant A, Scheller H, Camara B, et al. Arabidopsis VARIEGATED 3 encodes a chloroplast-targeted, zinc-finger protein required for chloroplast and palisade cell development. J Cell Sci. 2004;117:4807–18.
Yan J, Kandianis C, Harjes C, Bai L, Kim E, Yang X, Skinner D, Fu Z, Mitchell S, Li Q, et al. Rare genetic variation at Zea mays crtRB1 increases beta-carotene in maize grain. Nat Genet. 2010;42:322–7.
Pfeifer M, Kugler KG, Sandve SR, Zhan B, Rudi H, Hvidsten TR, Consortium IWGS, Mayer KFX, Olsen O-A. Genome interplay in the grain transcriptome of hexaploid bread wheat. Science. 2014;345:1250091.
We thank Saarah Kuzay, Quoc Thiem and Hana Leist for assistance with wheat grain collection and dissection, and Dr. Francis Cunningham (University of Maryland) for providing the pAC-BETA and pAC-ZEAX plasmids. This work was supported by the startup fund from University of California, Davis to LT. SY was supported by Chinese Scholarship Council Fellowship and UC Davis, Department of Plant Sciences Graduate Research Fellowship. JD acknowledges support from the Howard Hughes Medical Institute and the Gordon and Betty Moore Foundation.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
LT and XQ designed the experiment. XQ and KF collected and dissected wheat grains; XQ performed real-time qPCR analysis; XQ and SY conducted CCD enzyme assays; JD provided expertise in wheat genetics and sequence analysis and contributed to the revision of the manuscript; XQ, SY and LT analyzed the data; LT wrote the paper. All authors have read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Primers used for cloning of wheat CCD1 and CCD4 homoeologs. (DOCX 15 kb)
Homoeolog-specific primers used for real-time qPCR analysis. (PDF 17 kb)
Verification of homoeolog-specific primers using tetrasomic- nullisomic and ditelosomic lines of hexaploid wheat var. Chinese Spring. (PDF 377 kb)
E. coli expression and purification of wheat CCD1 and CCD4 proteins. (PDF 187 kb)
In vitro activities of wheat CCD-B1 and CCD-B4 homoeologs. (PDF 148 kb)
In vitro activities of wheat CCD-D1 and CCD-D4 homoeologs. (PDF 148 kb)
Characterization of wheat CCD-B1, CCD-B4, CCD-D1 and CCD-D4 activities in β-carotene producing E. coli cells. (PDF 45 kb)
In vitro activities of truncated wheat CCD4 homoeologs. (PDF 259 kb)
About this article
Cite this article
Qin, X., Fischer, K., Yu, S. et al. Distinct expression and function of carotenoid metabolic genes and homoeologs in developing wheat grains. BMC Plant Biol 16, 155 (2016). https://doi.org/10.1186/s12870-016-0848-7