Journal of Biomedicine and BiotechnologyVolume 2011 (2011), Article ID 731542, 7 pagesdoi:10.1155/2011/731542

Research Article

Cloning and Characterization of the -Carotene Desaturase Gene from Chlorella protothecoides CS-41

Meiya Li, Zhibing Gan, Yan Cui, Chunlei Shi, and Xianming Shi

MOST-USDA Joint Research Center for Food Safety and Bor Luh Food Safety Center, Department of Food Science & Technology, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai 200240, China

Received 31 May 2011; Revised 4 August 2011; Accepted 8 August 2011

Academic Editor: J. R. Botella

Copyright © 2011 Meiya Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

To elucidate the lutein biosynthesis pathway in the lutein-producing alga, Chlorella protothecoides CS-41, the -carotene desaturase gene (zds) was isolated from Chlorella protothecoides using the approach of rapid amplification of cDNA ends. The full-length cDNA sequence was 2031 bp and contained 1755 bp putative open reading frame which encodes a 584 amino acid deduced polypeptide whose computed molecular weight was 63.7 kDa. Sequence homology research indicated that the nucleotide and putative protein had sequence identities of 72.5% and 69.5% with those of the green alga Chlamydomonas reinhardtii, respectively. Phylogenetic analysis demonstrated that the ZDS from C. protothecoides CS-41 had a closer relationship with those of chlorophyta and higher plants than with those of other species. In addition, we also found that the zds gene expression was upregulated in response to light.

1. Introduction

Green algae in the genus Chlorella are widely distributed in freshwater environments. They are typically small, nonmotile, unicellular organisms with a single chloroplast and a rigid cell wall [1]. Chlorella protothecoides CS-41 is a strain that produces a high yield of lutein, which triggered our interest to investigate it.

Carotenoids are widely distributed in the world. They are a structurally diverse class of isoprenoids which are synthesized by all photosynthetic organisms and many nonphotosynthetic bacteria and fungi. They play a critical role in human nutrition and health [2]. In addition, they protect the photosynthetic apparatus in plants [3] from photo-oxidation. So far, more than 700 types of carotenoids have been found from natural sources [4]. For their extensive commercial and industrial uses, carotenoids especially lutein and astaxanthin are in high demand around the world.

Currently, for lutein production, the most widely used source is French marigold (Tagetes patula). However, mass plantation of marigold occupies a large land area and is easily influenced by season and climate. Previous research in our laboratory showed that Chlorella protothecoides CS-41 cultivated heterotrophically can produce a considerable amount of lutein [5]. By optimizing their cultivation conditions and the composition of the medium, one can improve the production of lutein for a higher yield [6, 7]. However, there are no reports about enhancing lutein production by genetic modification, although the genetic engineering technology has become more and more popular in the area of molecular biology. Therefore, the purpose of this research was to characterize the genes involved in lutein biosynthesis, with the expectation of increasing the lutein content in Chlorella protothecoides CS-41 by genetic engineering in the near future.

The carotenoid biosynthesis pathway in Chlorella protothecoides CS-41 is similar to that of higher plants, which is a complicated secondary metabolic system [8, 9] and is highly conserved in all carotenogenic organisms [10] (Figure 1). The process starts with the condensation of two geranylgeranyl diphosphate molecules to form phytoene. Phytoene is converted into colored -carotene in a two-step desaturation reaction by phytoene desaturase. Further desaturation and cyclization reactions lead to the biosynthesis of - and -carotene, which are subsequently converted into xanthophylls. ZDS is a key enzyme in carotenogenesis serving downstream of -carotene in the carotenoid biosynthetic pathway. In order to have a fuller understanding of the carotenogenic pathway of Chlorella protothecoides, we decided to isolate the cDNA of the zds gene and characterize this essential enzyme as well as other key enzymes such as PSY, PDS, and LCYE using a series of bioinformatics tools and functional analysis methods.

Figure 1: Putative biosynthesis pathway of carotenoids. The names of enzymes are abbreviated as follows: CrtB, bacterial phytoene synthase; CrtI, bacterial phytoene desaturase; PSY, phytoene synthase; PDS/CrtP, phytoene desaturase; ZDS/CrtQ, -carotene desaturase; CRTISO/CrtH, carotenoid isomerase.

2. Materials and Methods

2.1. Microbial Strains and Culture Conditions

The microalgal strain used in this study was Chlorella protothecoides CS-41, obtained from the CSIRO Marine Laboratories (Hobart, Australia). They were grown in modified Basal medium [11], which contained 10 g L⁻¹ glucose, at 28°C and shaken at 180 rpm. The cultures were later collected at the log or late log phase of growth.

Escherichia coli strain DH5α was used as the host for the multiplication of plasmids.

2.2. Genomic DNA and RNA Isolation

Genomic DNA was extracted using a modified cetyltrimethylammonium bromide (CTAB) method [12]. The total RNA was isolated from 15 mL of C. protothecoides CS-41 cells grown to the late log phase using RNAprep Plant Kit (Tiangen Biotech, China). DNA contamination was removed with DNaseI (Tiangen Biotech, China) according to the manufacturer’s instruction. The quantity of the total DNA and RNA was determined using a Nanodrop1000 spectrophotometer (Thermo scientific, USA).

2.3. Cloning of Full-Length zds cDNA and Genomic DNA

Degenerate primers were designed for the amplification of the partial cDNA from C. protothecoides CS-41. The primers were derived from the highly conserved nucleotide and amino acid sequences reported for the zds gene from four kinds of algae: Volvox carteri (XM_002954182), Chlorella NC_64A (gm_36_00005), Chlamydomonas reinhardtii (XM_001700734), and Dunaliella salina (HM754265); and four kinds of higher plants: Arabidopsis thaliana (NM_202493), Solanum lycopersicum (EF650012), Carica papaya (FJ812088), and Zea mays (EU974617).

One-step reverse transcriptase-polymerase chain reaction (RT-PCR) with the primers ZDSF and ZDSR (Table 1) was performed to amplify a portion of the zds cDNA. RT-PCR was performed using an RNA PCR Kit (AMV) ver.3.0 (TaKaRa, Dalian, China) according to the manufacturer’s protocol for cDNA synthesis, followed by PCR amplification (1 cycle of 94°C, 5 min; 20 cycles of 94°C, 30 s, 65–50°C, 30 s (decrease 0.5°C per cycle), and 72°C, 2 min; 5 cycles of 94°C, 30 s, 44°C, 30 s, and 72°C, 2 min; 1 cycle of 72°C, 10 min). The PCR product was purified with an AxyPrep PCR column (Axygen Biosciences, USA) to remove excess primers. Then the product was ligated into pMD-18T vector (TaKaRa, Dalian, China) according to the manufacturer’s instructions and chemically transformed into E. coli strain DH5α competent cells. Transformants were selected on LB-agar plates containing 100 mg/L ampicillin. Positive clones were selected and sequenced in both directions using M13 sequencing primers.

Table 1: PCR primers and targeted fragments of the zds gene in C. protothecoides CS-41.

Primers
Primer sequence (5′-3′)
Product size (bp)

Core fragment

413

    ZDSF
TCTTTGGCTGCTATTTCAACCTG

    ZDSR
ATCGTCACAGTTCAGGAAGC

5′-RACE

691

    ZRO1
TTCAGGAAGCCCAGAGCGTAG

    ZRI1
TTGCCGCCGAGGAAGAAGC

3′-RACE

1390

    ZFO1
ACCTGCTGCTCAAGGACCAC

    ZFI1
AGACGTGCGTGAGCTGGACTTC

    ZFI2
CCGATTGCTGACTACATCACTGCCA

ZDS-ORF

1755 for cDNA; 4587 for DNA

    ZDSF1
ATGCAGGCATCAGCGGTGGCC

    ZDSR1
TCACGCAGCCACCGCCTGCT

qRT-PCR

171

    ZDSRTF
TGGGCACCTCTCCCGTTG

    ZDSRTR
AGTTCAGGAAGCCCAGAGCG

Chlorella β-actin

166

    c-actinF
TGATGAAGATCCTCACAGAGC

    c-actinR
TCTGGCAGTTCAAAGGTCTT

Sets of specific primers were synthesized based on the sequence of putative insert for 5′ and 3′ rapid amplification of cDNA ends (RACE) [13]. ZRO1 and ZRI1 were used for 5′ RACE of zds gene, ZFO1 and ZFI1, ZFO1 and ZFI2 (Table 1) for 3′RACE of zds gene. RACE was performed using the 5′-Full RACE Kit and 3′-Full RACE Core Set Ver.2.0 (TaKaRa, Dalian, China) according to the manufacturer’s protocol. The RACE products were gel purified and sequenced as previously described. One pair of specific primer, ZDSF1 and ZDSR1 (Table 1), was designed from the sequences of the 5′ and 3′RACE fragments to amplify a full-length zds cDNA and its corresponding gene.

2.4. Sequence Alignment and Phylogenetic Analysis

Comparative and bioinformatics analysis of the nucleotide sequences and deduced amino acid sequences were carried out online at the websites (http://www.ncbi.nlm.nih.gov/) and (http://cn.expasy.org/). The nucleotide sequence, deduced amino acid sequence, and open reading frame [14] were analyzed, and the sequence comparison was conducted through a database search using BLAST programs. The phylogenetic analysis of zds gene from other plant species was aligned with Clustal X program version 1.83 using default parameters [15] and manual adjustments when necessary. Phylogenetic trees were constructed using MEGA (Molecular Evolutionary Genetics Analysis) version 4.0 [16]. The neighbor-joining (NJ) [17] method was used to construct the tree. In the NJ method, the distance was used to analyze the amino acid sequences. A total of 1000 repetitions were performed using the bootstrap method to determine the reliability of each node of the tree.

2.5. zds Expression Analysis in Response to Light

To analyze the light regulation pattern of the zds gene in C. protothecoides, the algal cells from late log phase were cultivated in the dark for more than 2 days, then collected by centrifugation at 5 500 g for 15 min under dark conditions. The pellet was resuspended in fresh medium without glucose, and then subjected to light treatment under a light intensity of 120 μmol m⁻² s⁻¹ for different induction times (0 h, 0.5 h, 1 h, 2 h, and 4 h). Each treatment was carried out in three replicates.

QPCR experiments were performed on an iQ Cycler (Bio-Rad, Watford, UK) with the specific primers (Table 1) and the SYBR ExScript RT-PCR kit (TaKaRa, Dalian, China). For each reaction, 1 μL of cDNA was added as template, 12.5 μL of 2 × SYBR Premix Ex Taq, 0.25 μM of forward primer, and 0.25 μM of reverse primer, and nuclease-free water were added to a final volume of 25 μL. The thermal cycling conditions were 3 min at 95°C for activation of the hot start Taq polymerase and 40 cycles of 95°C, 15 s, 60°C, 15 s, and 72°C, 30 s. Melting curve analysis following each qPCR was performed to assess product specificity. Each qPCR measurement was carried out independently at least three times, and the mean value was used for quantification. The method was used to analyze the relative changes in gene expression, and the expression of -actin gene was used as a normalized control, and the expression of the untreated samples was used as a negative control.

3. Results and Discussion

3.1. Cloning of C. protothecoides CS-41 zds Gene

Carotenoids comprise a diverse group of multifunctional, taxonomically widespread, and biotechnologically important pigments. Their biosynthetic pathway is a model system for understanding the evolution of secondary metabolism. ZDS is an essential carotenogenic enzyme that carries out the conversion from -carotene to lycopene in eukaryotes. So far, nuclear genome sequences of several eukaryotic algae, including C. reinhardtii [18, 19] and Chlorella NC_64A [1, 20] have been completed, but the characteristics of ZDS in these eukaryotic algae have not yet been released by the authors. As a result of the incomplete molecular investigation of carotenogenesis of eukaryotic algae, the carotenogenic enzymes of C. protothecoides CS-41 have not been isolated and characterized thoroughly, and the carotenogenic pathway of C. protothecoides CS-41 is still unclear. Nonetheless, several DNA sequences of carotenogenic enzymes of C. protothecoides CS-41 have been determined, such as the psy, pds, and lcye genes. In this study, the cDNA of zds was isolated using the RACE-PCR approach.

TD-PCR with primers ZDSF and ZDSR (Table 1) generated a predicted 413-bp fragment (Figure 2, lane 1). BLAST analysis showed that the nucleotide sequence of this fragment shared about 76% identity with that of C. reinhardtii, which suggested that the sequence of the fragment was derived from a putative ζ-carotene desaturase.

Figure 2: The PCR products of the zds gene of C. protothecoides CS-41. Lane 1: The core fragment of zds; lane 2: 5′-RACE fragment of zds; lanes 3 and 4: 3′-RACE fragments of zds; lane 5: the ORF fragment of zds; M: DL 2000.

With this sequence information, specific primers were designed for 5′ and 3′RACE of the related gene. A 691-bp fragment (Figure 2, lane 2) was generated with the primers ZRO1/ZRI1; a 486-bp fragment (Figure 2, lane 3) generated at the first cycle of 3′RACE with the primers ZFO1/ZFI1; a 975-bp fragment (Figure 2, lane 4) at the second cycle of 3′RACE with ZFO1/ZFI2. They were confirmed by sequencing to be the 5′ and 3′ regions of a phytoene desaturase gene, respectively. RT-PCR with the primer pair, ZDSF1/ZDSR1, generated a 1755-bp fragment (Figure 2, lane 5), which was identified as a full-length zds cDNA (submitted to GenBank, accession no. GU269622).

To characterize the corresponding gene of the zds cDNA, genomic PCR was performed. A 4587-bp fragment (submitted to GenBank, accession No. GU269623) was generated and sequenced. Analysis of the obtained nucleotide sequence revealed that the product was the zds gene. After analysis using online software (http://www.ncbi.nlm.nih.gov/IEB/Research/Ostell/Spidey/), this gene was found to consist of fourteen exons and thirteen introns (Figure 3).

Figure 3: Diagram of exons and introns of zds gene in C. protothecoides CS-41. Ten exons were (1) 1st bp to 112nd bp; (2) 217th bp to 332nd bp; (3) 509th bp to 586th bp; (4) 893rd bp to 1001st bp; (5) 1194th bp to 1285th bp; (6) 1500th bp to 1598th bp; (7) 1812nd bp to 1995th bp; (8) 2243rd bp to 2438th bp; (9) 2707th bp to 2851st bp; (10) 3070th bp to 3171st bp; (11) 3366th bp to 3427th bp; (12) 3589th bp to 3748th bp; (13) 4060th bp to 4238th bp; and (14) 4467th bp to 4587th bp.

3.2. Characterization of the Deduced ZDS Protein

By using the Compute pI/Mw tool at the website http://web.expasy.org/compute_pi/, the open reading frame of zds cDNA encodes a protein of 584 amino acids, which has a calculated molecular mass of 63.67 kDa and an isoelectric point of 8.68. The protein contained 64 strongly basic amino acids (Lys and Arg), 58 strongly acidic amino acids (Asp and Glu), 219 hydrophobic amino acids (Ala, Ile, Leu, Phe, Trp and Val), and 145 polar amino acids (Asn, Cys, Gln, Ser, Thr and Tyr).

A sequence homology search revealed that the putative protein had high homology with ZDS sequences of other Chlorella species, such as 64.8% identities and 69.1% positives with that from Chlorella NC_64A. ZDS was also highly similar to that from C. reinhardtii (69.5% identities, 75.4% positives) and D. salina (57.4% identities, 67.0% positives), further confirming that ZDS belongs to the algal ZDS family. The deduced ZDS similar to the PDS protein contained a conserved dinucleotide-binding site (GXGX₂GX₃AX₂LX₃GX₆EX₅GG) at the N-terminus and a carotenoid-binding domain (AGDYTKQRYLASMEGATFSG) at the C-terminus [21] (Figure 4).

(a)

(b)

Figure 4: Sequence alignment of deduced PDS and ZDS with those from other organisms. Species designations are A. thaliana (Ath), Z. mays (Zma), C. reinhardtii (Cre), C. protothecoides (Cpr), and D. salina (Dsa). Amino acids in black boxes indicate identical residues, while gray boxes indicate conserved residues. The blank frame shows the consensus motif of dinucleotide-binding domain. The red frame indicates the postulated carotenoid-binding domain.

BLASTP analysis suggested that this protein contains a FAD binding domain belonging to the Pyridine nucleotide-disulphide oxidoreductase superfamily that is involved in oxidative metabolism of a variety of hydrocarbons [21]. Prediction of a transmembrane region was performed using the TMHMM server version 2.0, which suggested that the C. protothecoides ZDS was not a potential transmembrane protein. The SignalP-HMM analysis results suggested that there is a signal peptide in this polypeptide with the maximum cleavage site probability of 0.307 between positions 16 and 17 (Figure 5), which implies that it belongs to the outside protein. The function of this protein in C. protothecoides CS-41 needs further investigation.

Figure 5: Analysis of the signal peptide of the ZDS from C. protothecoides CS-41.

3.3. Molecular Evolution Analysis

The zds gene was first cloned from C. protothecoides CS-41, prompting an investigation of their evolutionary position in the phylogenetic tree of various genes. Using MEGA 4.0, phylogenetic analysis of 20 ZDS protein sequences demonstrated that four clusters were displayed on the phylogenetic tree as follows: higher plants, chlorophyta, cyanobacteria, and bacteria. As displayed in the phylogenetic tree, the ZDS from C. protothecoides had a closer relationship with those of chlorophyta and higher plants than with those of cyanobacteria (Figure 6). This suggests that ZDS was derived from a common ancestor gene from bacteria via cyanobacteria to plants. This is consistent with the previous conclusion from the phylogenetic analysis that PSY and PDS from C. protothecoides are more closely related to those from other algae and higher plants.

Figure 6: Phylogenic tree of the zds sequences from various species. The phylogeny was derived using Neighbor-joining analysis. Bootstrap values >50% are indicated near the base of each branch (1000 replicates). Horizontal branch lengths represent relative evolutionary distances with the scale bar corresponding to 0.1 amino acid substitutions per site.

3.4. Gene Expression Response to Light

Light, as a particularly important environmental factor for all photosynthetic organisms, is essential for photosynthesis, and also involved in the regulation of many aspects of plant growth and development [22]. Carotenoids are ubiquitous and essential components of photosynthetic tissues in unicellular microalgae because they participate in the light harvesting process and prevent photooxidative damage of the photosynthetic apparatus [23]. Hence, investigating gene expression in the biosynthetic pathway of carotenoids is very important. As a key enzyme in the pathway of carotenoids, the transcription level of the zds gene was investigated after treatment under light by using real-time PCR. The result showed that the expression of zds gene increased to 5.3 times than that in the dark, which indicated that the expression of zds gene was upregulated in response to light (Figure 7). In addition to zds gene, we also investigated the expression of the psy, pds, and lcye genes in the unicellular green alga C. protothecoides (data not shown). Our results revealed that all of these genes in the biosynthesis pathway showed upregulation in response to light.

Figure 7: Gene expression levels determined from cells exposed to light for different lengths of time in C. protothecoides CS-41. Data (mean + SEM) are combined from three independent experiments. *indicated that gene expression levels were highly significantly different from that at 0 h (); #showed that gene expression levels were significantly different from that at 0 h ().

Since metabolically engineered microorganisms [24] and transgenic plants [25–28] are widely used in carotenoid production and carotenoid breeding to meet the high global demand, large carotenogenic gene pool should be provided first from the targeted engineered microorganisms for carotenoid production. The carotenoid biosynthesis pathway in higher plants and C. reinhardtii has been described in great detail using genetic, biochemical, and molecular data; however, this is the first research on the carotenoid pathway of a lutein-producing microalga, Chlorella protothecoides CS-41. Once these genes of key enzymes have been successfully isolated and characterized including the zds gene, it will lead to the elucidation of the carotenogenic pathway of C. protothecoides as well as the evolution of secondary metabolism.

Acknowledgments

This work was supported by National Natural Science Foundation of China (Grant nos. 31000779 and U1031003) and National 863 program of China (Grant no. 2006AA02Z226), and Shanghai Jiao Tong University Innovation Fund for Postgraduates. The authors also thank Professor James T. Staley at University of Washington for the critical review of this paper.

References

N. A. Eckardt, “The Chlorella genome: big surprises from a small package,” Plant Cell, vol. 22, no. 9, p. 2924, 2010.

N. I. Krinsky and E. J. Johnson, “Carotenoid actions and their relation to health and disease,” Molecular Aspects of Medicine, vol. 26, no. 6, pp. 459–516, 2005.

J. Hirschberg, “Carotenoid biosynthesis in flowering plants,” Current Opinion in Plant Biology, vol. 4, no. 3, pp. 210–218, 2001.

F. Ladislav, P. Vere, S. Karel, and V. Karel, “Reliability of carotenoid analyses: a review,” Current Analytical Chemistry, vol. 1, no. 1, pp. 93–102, 2005.

X. M. Shi, Y. Jiang, and F. Chen, “High-yield production of lutein by the green microalga Chlorella protothecoides in heterotrophic fed-batch culture,” Biotechnology Progress, vol. 18, no. 4, pp. 723–727, 2002.

X. M. Shi, H. J. Liu, X. W. Zhang, and F. Chen, “Production of biomass and lutein by Chlorella protothecoides at various glucose concentrations in heterotrophic cultures,” Process Biochemistry, vol. 34, no. 4, pp. 341–347, 1999.

X. Shi, Z. Wu, and F. Chen, “Kinetic modeling of lutein production by heterotrophic Chlorella at various pH and temperatures,” Molecular Nutrition and Food Research, vol. 50, no. 8, pp. 763–768, 2006.

C. Liang, F. Zhao, W. Wei, Z. Wen, and S. Qin, “Carotenoid biosynthesis in cyanobacteria: structural and evolutionary scenarios based on comparative genomics,” International Journal of Biological Sciences, vol. 1, no. 4, pp. 197–207, 2006.

F. Wang, J. G. Jiang, and Q. Chen, “Progress on molecular breeding and metabolic engineering of biosynthesis pathways of C(30), C(35), C(40), C(45), C(50) carotenoids,” Biotechnology Advances, vol. 25, no. 3, pp. 211–222, 2007.

C. Sieiro, M. Poza, T. de Miguel, and T. G. Villa, “Genetic basis of microbial carotenogenesis,” International Microbiology, vol. 6, no. 1, pp. 11–16, 2003.

X. M. Shi, F. Chen, J. P. Yuan, and H. Chen, “Heterotrophic production of lutein by selected Chlorella strains,” Journal of Applied Phycology, vol. 9, no. 5, pp. 445–450, 1997.

C. N. Stewart Jr. and L. E. Via, “A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications,” BioTechniques, vol. 14, no. 5, pp. 748–750, 1993.

M. A. Frohman, M. K. Dush, and G. R. Martin, “Rapid production of full-length cDNAs from rare transcripts: amplification using a single gene-specific oligonucleotide primer,” Proceedings of the National Academy of Sciences of the United States of America, vol. 85, no. 23, pp. 8998–9002, 1988.

G. Peers, T. B. Truong, E. Ostendorf et al., “An ancient light-harvesting protein is critical for the regulation of algal photosynthesis,” Nature, vol. 462, no. 7272, pp. 518–521, 2009.

J. D. Thompson, T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins, “The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools,” Nucleic Acids Research, vol. 25, no. 24, pp. 4876–4882, 1997.

K. Tamura, J. Dudley, M. Nei, and S. Kumar, “MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0,” Molecular Biology and Evolution, vol. 24, no. 8, pp. 1596–1599, 2007.

W. Zhang and Z. Sun, “Random local neighbor joining: a new method for reconstructing phylogenetic trees,” Molecular Phylogenetics and Evolution, vol. 47, no. 1, pp. 117–128, 2008.

G. L. Wheeler, D. Miranda-Saavedra, and G. J. Barton, “Genome analysis of the unicellular green alga Chlamydomonas reinhardtii indicates an ancient evolutionary origin for key pattern recognition and cell-signaling protein families,” Genetics, vol. 179, no. 1, pp. 193–197, 2008.

J. Godman and J. Balk, “Genome analysis of Chlamydomonas reinhardtii reveals the existence of multiple, compartmentalized iron-sulfur protein assembly machineries of different evolutionary origins,” Genetics, vol. 179, no. 1, pp. 59–68, 2008.

G. Blanc, G. Duncan, I. Agarkova et al., “The Chlorella variabilis NC64A genome reveals adaptation to photosymbiosis, coevolution with viruses, and cryptic sex,” Plant Cell, vol. 22, no. 9, pp. 2943–2955, 2010.

F. X. Cunningham Jr., Z. Sun, D. Chamovitz, J. Hirschberg, and E. Gantt, “Molecular structure and enzymatic function of lycopene cyclase from the cyanobacterium Synechococcus sp strain PCC7942,” Plant Cell, vol. 6, no. 8, pp. 1107–1121, 1994.

M. M. Neff, C. Fankhauser, and J. Chory, “Eight: an indicator of time and place,” Genes and Development, vol. 14, no. 3, pp. 257–271, 2000.

H. A. Frank and R. J. Cogdell, “Carotenoids in photosynthesis,” Photochemistry and Photobiology, vol. 63, no. 3, pp. 257–264, 1996.

C. Schmidt-Dannert, D. Umeno, and F. H. Arnold, “Molecular breeding of carotenoid biosynthetic pathways,” Nature Biotechnology, vol. 18, no. 7, pp. 750–753, 2000.

J. A. Paine, C. A. Shipton, S. Chaggar et al., “Improving the nutritional value of golden rice through increased pro-vitamin A content,” Nature Biotechnology, vol. 23, no. 4, pp. 482–487, 2005.

M. Aluru, Y. Xu, R. Guo et al., “Generation of transgenic maize with enhanced provitamin A content,” Journal of Experimental Botany, vol. 59, no. 13, pp. 3551–3562, 2008.

C. Bai, R. M. Twyman, G. Farré et al., “A golden era-pro-vitamin A enhancement in diverse crops,” In Vitro Cellular and Developmental Biology—Plant, vol. 47, no. 2, pp. 205–221, 2011.

L. J. M. Ducreux, W. L. Morris, P. E. Hedley et al., “Metabolic engineering of high carotenoid potato tubers containing enhanced levels of -carotene and lutein,” Journal of Experimental Botany, vol. 56, no. 409, pp. 81–89, 2005.