Open Access

De novo production of the flavonoid naringenin in engineered Saccharomyces cerevisiae

  • Frank Koopman1, 2, 3,
  • Jules Beekwilder2, 4, 5,
  • Barbara Crimi1, 2, 3,
  • Adele van Houwelingen4,
  • Robert D Hall2, 4, 5,
  • Dirk Bosch2, 4, 5,
  • Antonius JA van Maris1, 3,
  • Jack T Pronk1, 2, 3 and
  • Jean-Marc Daran1, 2, 3Email author
Microbial Cell Factories201211:155

DOI: 10.1186/1475-2859-11-155

Received: 4 September 2012

Accepted: 15 November 2012

Published: 8 December 2012

Abstract

Background

Flavonoids comprise a large family of secondary plant metabolic intermediates that exhibit a wide variety of antioxidant and human health-related properties. Plant production of flavonoids is limited by the low productivity and the complexity of the recovered flavonoids. Thus to overcome these limitations, metabolic engineering of specific pathway in microbial systems have been envisaged to produce high quantity of a single molecules.

Result

Saccharomyces cerevisiae was engineered to produce the key intermediate flavonoid, naringenin, solely from glucose. For this, specific naringenin biosynthesis genes from Arabidopsis thaliana were selected by comparative expression profiling and introduced in S. cerevisiae. The sole expression of these A. thaliana genes yielded low extracellular naringenin concentrations (<5.5 μM). To optimize naringenin titers, a yeast chassis strain was developed. Synthesis of aromatic amino acids was deregulated by alleviating feedback inhibition of 3-deoxy-d-arabinose-heptulosonate-7-phosphate synthase (Aro3, Aro4) and byproduct formation was reduced by eliminating phenylpyruvate decarboxylase (Aro10, Pdc5, Pdc6). Together with an increased copy number of the chalcone synthase gene and expression of a heterologous tyrosine ammonia lyase, these modifications resulted in a 40-fold increase of extracellular naringenin titers (to approximately 200 μM) in glucose-grown shake-flask cultures. In aerated, pH controlled batch reactors, extracellular naringenin concentrations of over 400 μM were reached.

Conclusion

The results reported in this study demonstrate that S. cerevisiae is capable of de novo production of naringenin by coexpressing the naringenin production genes from A. thaliana and optimization of the flux towards the naringenin pathway. The engineered yeast naringenin production host provides a metabolic chassis for production of a wide range of flavonoids and exploration of their biological functions.

Keywords

Saccharomyces cerevisiae Naringenin de novo Flavonoids Metabolic engineering

Background

In recent years, plant flavonoids, which comprise a family of over 9000 compounds, have attracted a tremendous increase in research interest[13]. This interest is mainly attributed to highly promising human health applications of specific flavonoids[48]. The biological activities of flavonoid compounds have been investigated in relation to a multitude of human pathological conditions, including cancer, diabetes, obesity and Parkinson’s disease[6, 914]. The identified mechanisms of action include scavenging of oxygen radicals, anti-inflammatory, antiviral and antitumor activities[15, 16].

Both for health related research and commercial nutritional applications, availability of sufficient amounts of defined flavonoid preparations is important. To date, flavonoid production mostly relies on isolation from plants. However, investigation and subsequential industrialization from plants is hampered by their low production efficiency. In addition to the low growth rate of some of the producing plants, extraction and separation of flavonoids with highly related structures complicate plant-based production, thereby impeding progress in the exploration of the biological activities of flavonoids[13, 14]. Although flavonoids can be produced chemically, efficient production of flavonoids by organic synthesis is severely hindered by the complexity of the molecules, as well as by the necessity of utilizing toxic chemicals and extreme reaction conditions[13, 14].

In response to the poor production efficiency from plants and chemical synthesis, research groups have directed their attention to the heterologous production of flavonoids in microorganisms such as Escherichia coli and Saccharomyces cerevisiae, using metabolic engineering and synthetic biology[14, 1720]. In plants, formation of the central metabolite naringenin involves six steps catalyzed by phenylalanine ammonia lyase (PAL), cinnamate 4-hydroxylase (C4H) and its associated cytochrome P450 reductase (CPR), 4-coumaric acid-CoA ligase (4CL), chalcone synthase (CHS) and chalcone isomerase (CHI) (Figure1)[21, 22]. Alternatively, a tyrosine ammonia lyase (TAL) can convert tyrosine directly to coumaric acid and circumvent the use of membrane bound P450 related enzymes, which may pose challenges in E. coli[14, 18]. Additionally, in plants, these enzymes and their corresponding genes are often represented by several isoforms, which may differ in substrate preference or kinetic properties[23]. Moreover, the different isoforms are proposed to be organized in one or more enzyme complexes that may promote substrate channeling. In Arabidopsis thaliana, interaction between CHS, CHI and an enzyme downstream of naringenin, flavonol 3-hydroxylase were demonstrated based on affinity chromatography and immunoprecipitation assays[2426].
https://static-content.springer.com/image/art%3A10.1186%2F1475-2859-11-155/MediaObjects/12934_2012_Article_781_Fig1_HTML.jpg
Figure 1

Schematic representation of the engineered naringenin production pathway in S. cerevisiae. Six A. thaliana genes were overexpressed: PAL1 (phenylalanine ammonia lyase), C4H (Cinnamate 4-hydroxylase), CPR1 (cytochrome P450 reductase), 4CL3 (4-coumaric acid-CoA ligase), CHS3 (chalcone synthase) and CHI1 (chalcone isomerase), and one gene from Rhodobacter capsulatus; Tal1 (tyrosine ammonia lyase). Dashed lines indicate feedback inhibition. Grey arrows indicate the S. cerevisiae pathway for phenylethanol production. Bold dark grey arrows indicate the naringenin production pathway as described for A. thaliana[22]. Aro3/Aro4: 3-deoxy-D-arabino-heptulosonate-7-phosphate (DAHP) synthase, Pdc1, 5, 6; pyruvate decarboxylases, Aro10; phenylpyruvate decarboxylase.

Expression of enzyme combinations, originating from a variety of host organisms has yielded microbial strains capable of producing the key flavonoid precursor naringenin (Figure1). Several reports describe successful biotransformation processes in which a phenylpropanoid precursor, such as coumaric acid, is converted into naringenin by metabolically engineered E. coli or S. cerevisiae[1, 14, 17, 19, 20, 27, 28]. The highest naringenin titers obtained through biotransformation were achieved in E. coli, with naringenin titers reaching 1.74 mM (474 mg·l-1) from 2.6 mM coumaric acid supplied to the medium[1]. The reliance on ‘expensive’ phenylpropanoid precursors might represent a major hurdle for economically feasible flavonoid production[18]. Hitherto, only one study reports de novo naringenin production from glucose. Using an engineered E. coli strain, naringenin titers of up to 106.5 μM (29 mg·l-1) were obtained[18].

S. cerevisiae has several attractive characteristics as a metabolic engineering platform for flavonoid production. In addition to its excellent accessibility to molecular and synthetic biology techniques[29, 30], its eukaryotic nature may facilitate functional expression of plant-derived flavonoid-biosynthetic genes. For example, S. cerevisiae can functionally express cytochrome P450-containing enzymes and its subcellular compartmentation is comparable to that of plant cells[31]. Finally, its GRAS (generally recognized as safe) status facilitates subsequent application for the production of pharma- and nutraceuticals.

The goal of the present study was to define a metabolic engineering strategy for de novo production of naringenin by S. cerevisiae, using glucose as sole carbon source. For optimal synergistic activity, the flavonoid biosynthetic genes PAL1, C4H, CPR1, 4CL3, CHS3 and CHI1 used in this study were derived from a single plant species, A. thaliana and selected for in planta co-expression profiles. After expression of the plant pathway genes, optimization of naringenin production was explored by engineering of precursor supply to the naringenin pathway and by reducing the formation of byproducts derived from yeast metabolism.

Results

Selection of naringenin biosynthetic genes from A. thaliana

As a first step towards heterologous expression of naringenin in S. cerevisiae, flavonoid biosynthetic genes were selected from A. thaliana. It has been proposed that, in plants, flavonoid biosynthetic enzymes can be organized in protein complexes, where the proteins involved in the pathway are spatially co-localized. The combination of isoenzymes in such complexes is hypothesized to facilitate the synthesis of different flavonoids[24]. In addition, complex formation might enhance pathway activity and carbon flux via metabolic channeling of intermediates[24]. Moreover, formation of toxic intermediates is kept locally and are sequentially converted[32]. Protein complex formation is likely to require specific protein-protein interactions and therefore, co-evolution of the protein structures of the enzymes that form these complexes. In A. thaliana, several isoenzymes exist for most of the reactions leading to naringenin formation (Table1). In order to identify the best set of isoenzymes for naringenin production and maintain the potential benefit of plant-specific protein-protein interactions upon expression in S. cerevisiae, a subset of isoenzymes from one single species, A. thaliana, was identified and subsequently subjected to expression profile correlation analysis. The subset of isoenzymes was selected based on previous reports indicating either the involvement in the lignin biosynthetic pathway[33] and/or the flavonoid biosynthetic pathway (Figure1), since these pathways share the same reactions up to the metabolic branch-point coumaroyl-CoA. The specific genes from this isoenzyme subset were then further assessed for the correlation of their expression levels in a dataset comprising 392 transcriptome studies of A. thaliana[34].
Table 1

A. thaliana genes included in the expression correlation studies

Gene

Gene product

A. thaliana locus

PAL1

phenylalanine ammonia lyase

at2g37040

PAL2

phenylalanine ammonia lyase

at3g53260

PAL3

phenylalanine ammonia lyase

at5g04230

PAL4

phenylalanine ammonia lyase

at3g10340

C4H

trans-cinnamate 4-monooxygenase

at2g30490

4CL1

4-coumaric acid-CoA ligase

at1g51680

4CL2

4-coumaric acid-CoA ligase

at3g21240

4CL3

4-coumaric acid-CoA ligase

at1g65060

4CL4

4-coumaric acid-CoA ligase

at3g21230

CHS1

Chalcone synthase

at4g00040

CHS2

Chalcone synthase

at4g34850

CHS3

Chalcone synthase

at5g13930

CHI1

Chalcone isomerase

at5g05270

CHI2

Chalcone isomerase

at3g63170

CHI3

Chalcone isomerase

at3g55120

F3H

Flavanone 3-hydroxylase

at3g51240

FLS

Flavonol synthase

at5g08640

HCT

Quinate O-hydroxycinnamoyltransferase

at5g48930

CCoOMT

Caffeoyl-CoA 3-O-methyltransferase

at4g34050

C3H

p-coumaric acid 3-hydroxylase

at2g40890

Bold letters indicate genes used in the construction of engineered S. cerevisiae strains.

From the expression correlation analysis, using the BAR Expression Angler[35], two clear modules could be distinguished (Figure2). The first module comprises flavonoid biosynthetic genes, including the 4-coumaric acid-CoA ligase gene (4CL3), the chalcone synthase gene (CHS3), two chalcone isomerase genes (CHI1, CHI3) and a number of genes that mediate the further modification of naringenin. The second, a lignin module, comprised two phenylalanine ammonia lyase (PAL1, PAL2) genes, the trans-cinnamate 4-monooxygenase (C4H) gene and three 4-coumaric acid-CoA ligase genes (4CL1, 4CL2, 4CL4), in addition to a number of genes known to be involved in formation of lignin and phenolic esters.
https://static-content.springer.com/image/art%3A10.1186%2F1475-2859-11-155/MediaObjects/12934_2012_Article_781_Fig2_HTML.jpg
Figure 2

Co-expression analysis of candidate genes for flavonoid biosynthesis in A. thaliana . Correlations between the expression of a set of target genes (see Table1) were established by bioinformatic analysis of 392 A. thaliana transcriptome datasets[34]. Circles represent genes as described in Table1. Light grey circles represent genes where the expression correlates with that of genes involved in lignin metabolism (represented by C3H, HCT, CCoOMT), dark grey circles represent genes with expression correlating with genes of the flavonoid pathway (represented by FLS, F3H and several UGTs). White circles represent genes with expression profiles which do not correlate with genes involved in either the lignin or the flavonoid pathway. Solid lines represent Pearson correlations of expression above 0.7, dotted lines represent Pearson correlations between 0.6 and 0.7.

Since these results are in good agreement with a previous study on expression correlation in plants[36], the following genes for the pathway towards naringenin were selected for heterologous expression in S. cerevisiae: the flavonoid biosynthetic genes 4CL3 (at1g65060), CHS1 (at4g00040) and CHI1 (at5g05270), as well as, C4H (at2g30490) and PAL1 (at2g37040) from the lignin biosynthetic genes[37]. PAL1 and CHI1 were preferred over PAL2 and CHI3, as their expression profiles showed stronger correlations with C4H and 4CL3 respectively (Figure2).

Construction and evaluation of a naringenin producing strain

To enable naringenin production in S. cerevisiae, one episomal and one integrative expression vector were constructed which together, harbor the five flavonoid biosynthetic genes. Additionally, activation of the cytochrome P450 C4H requires a cytochrome P450 reductase (CPR). To choose the best candidate gene, the two A. thaliana CPR variants (CPR1 or CPR2) were separately included in the pathway engineering strategy. First, the centromeric episomal plasmid pUDE172, carrying PAL1 and yeast codon-optimized versions of the A. thaliana C4H (coC4H) and CPR1 (coCPR1) genes was constructed. In a same manner, separate construct, PAL1 and coC4H were combined with the coCPR2 gene. Subsequently, the integration plasmid pUDI065 was constructed, which carried the non yeast optimized A. thaliana genes at4CL3, atCHS1 and atCHI1. Only the, C4H and CPR1 in the naringenin production pathway were codon optimized since these genes hold the majority of rare codons for yeast (8), compared to the other genes (PAL1, 4CL3, CHS3 and CHI1) (3). When coCPR2 was expressed instead of coCPR1, significantly lower naringenin titers were observed (data not shown); coCPR2 was therefore not used in further experiments. Introduction of the coCPR1 version of the centromeric expression vector and the integration vector yielded the S. cerevisiae strain IMU011. In shake-flask cultures containing synthetic medium and with glucose as the sole carbon source, this strain produced naringenin to a concentration of 5.4 μM (Figure3A). Interestingly, naringenin was measured extracellularly, although its export mechanism has not yet been elucidated.
https://static-content.springer.com/image/art%3A10.1186%2F1475-2859-11-155/MediaObjects/12934_2012_Article_781_Fig3_HTML.jpg
Figure 3

Naringenin production in S. cerevisiae. A) Increase in heterologous production of naringenin. B) Phenylethanol production in engineered strains. () IMU011 (atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,at4CL3↑), () IMX183 (aro3 Δ ARO4 G226S atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,at4CL3↑), (□) IMX185 (aro3 Δ ARO4 G226S pdc5 Δ, pdc6 Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,at4CL3↑), (■) IMX106 (aro3 Δ ARO4 G226S aro10 Δ, pdc5 Δ, pdc6 Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,at4CL3↑), (Δ) non-naringenin producing reference strain CEN.PK2-1C. Cultures were grown in shake flasks on synthetic medium containing 20 g·l-1 glucose and appropriate growth factors to supplement the auxotrophic requirements of the strains. All cultures were performed in triplicate. Error bars denote standard deviation.

Alleviation of tyrosine feedback inhibition of yeast 3-deoxy-D-arabinose-heptulosonate-7-phosphate (DAHP) synthase

Since phenylalanine is a key intermediate in naringenin biosynthesis (Figure1), further metabolic engineering was directed at improving its intracellular availability. Aromatic amino acid biosynthesis in S. cerevisiae is subject to strong feedback inhibition by phenylalanine and tyrosine. Introduction of a tyrosine insensitive ARO4 allele (ARO4 G226S )[38, 39] in conjunction with the deletion of the other allele of the 3-deoxy-d-arabinose-heptulosonate-7-phosphate (DAHP) synthase (ARO3) has previously been shown to cause a 4-fold increase in flux through the aromatic amino acid pathway[40]. To test the potential impact of this genetic intervention on naringenin synthesis, S. cerevisiae strain IMX183 (ARO4 G226S aro3Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,at4CL3↑) was constructed. The deregulation of DAHP synthase in this strain lead to a 2-fold increase of the naringenin titer in shake flask cultures (10.4 μM) (Figure3). However, this increase was accompanied by a strongly enhanced accumulation of phenylethanol, whose extracellular concentration was increased 100-fold to a level ca. 20 fold higher than the improved naringenin concentration (Figure3B). Therefore, to improve further naringenin production, we adopted a strategy to reduce the diversion of aromatic amino acid biosynthesis into this byproduct.

Elimination of competing phenylpyruvate decarboxylase activity

In S. cerevisiae, phenylethanol is produced via the Ehrlich pathway. Decarboxylation of phenylpyruvate, the 2-oxo acid associated to phenylalanine, yields phenylacetaldehyde, which is subsequently reduced into phenylethanol and/or oxidized to phenylacetate[4143]. Decarboxylation of phenylpyruvate can be catalyzed by four different thiamine pyrophosphate-dependent 2-oxo acid decarboxylases encoded by ARO10, PDC1, PDC5, and PDC6[43, 44]. Recent work in our group has demonstrated that, among the four decarboxylases capable of phenylpyruvate decarboxylase, Pdc1 and Pdc6 showed a much lower affinity and decarboxylation rate of phenylpyruvate than Pdc5 and, in particular, Aro10[44]. Since absence of all three pyruvate decarboxylase genes (PDC1, PDC 5 and PDC6) abolishes growth on glucose in synthetic media, strains only retaining PDC1 (aro10Δ, pcd5Δ, pdc6Δ) were constructed[45]. The intermediate strain IMX185 (ARO4 G226S , aro3Δ, pdc5Δ, pdc6Δ) did not show a reduced phenylethanol titer compared to its ancestor IMX183 (ARO4 G226S , aro3Δ) (Figure3B). In contrast, the naringenin producing strain IMX106 (ARO4 G226S , aro3Δ, pdc5Δ, pdc6Δ, aro10Δ) that also carries the ARO10 deletion exhibited a 22-fold lower concentration of extracellular phenylethanol (44 μM) compared to both ancestor strains IMX183 and IMX185 (955 μM). This strong reduction in phenylethanol production coincided with a 3-fold increase of extracellular naringenin (up to 46.5 μM) (Figure3A), indicating that reduction of the flux through the Ehrlich pathway had indeed led to substantial improvement in naringenin production. Analysis of culture supernatants revealed that strain IMX106 also produced coumaric acid up to a concentration of 200 μM (Figure4). Moreover, this strain showed increased production of another aromatic metabolite that, based on HPLC spectra and LC-MS analysis, was identified as phloretic acid (Figure4). Production of phloretic acid was only observed in cultures that also produced coumaric acid and which expressed 4CL3 (data not shown). This suggests that phloretic acid is most likely produced from coumaroyl-CoA (Figure4).
https://static-content.springer.com/image/art%3A10.1186%2F1475-2859-11-155/MediaObjects/12934_2012_Article_781_Fig4_HTML.jpg
Figure 4

Product formation by S. cerevisiae IMX106 ( aro3 Δ ARO4 G226S aro10 Δ, pdc5 Δ, pdc6 Δ, atPAL1, coC4H, coCPR1, atCHI1, atCHS3 , at4CL3). Formation of naringenin (),coumaric acid () and phloretic acid (□). Cultures were grown in shake flasks on synthetic medium containing 20 g·l-1 glucose. All cultures were performed in triplicate. Error bars denote standard deviation.

Alleviate the bottleneck downstream of coumaric acid

Naringenin chalcone synthase, which catalyzes the formation of chalcone by condensing coumaroyl-CoA with three molecules of malonyl-CoA, is known to be an enzyme with low catalytic activity[46]. The transient accumulation and later reconsumption of coumaric acid in shake flask cultures suggested that reactions downstream of coumaric acid were limiting naringenin production. To test whether the capacity of chalcone synthase was indeed controlling flux through the pathway, two additional copies of the coCHS3 gene were introduced into strain IMX106 on an episomal plasmid (pUDE188), yielding strain IMX197. The two additional copies resulted in a 2.5 fold increase in naringenin accumulation (134.5 μM) in shake flask cultures (Figure5A), indicating that coCHS3 was indeed a limiting step in the naringenin production. Additional coCHS3 copies also caused a decreased production of phloretic acid (Figure5B), consistent with the hypothesis that phloretic acid production occurs when coumaroyl-CoA cannot be efficiently converted to naringenin.
https://static-content.springer.com/image/art%3A10.1186%2F1475-2859-11-155/MediaObjects/12934_2012_Article_781_Fig5_HTML.jpg
Figure 5

Stepwise increase in naringenin formation by S cerevisiae. Formation of A) naringenin and B) phloretic acid in the engineered strains IMX106 ()(aro3 Δ ARO4 G226S aro10 Δ, pdc5 Δ, pdc6 Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,at4CL3↑), IMX197 ()(aro3 Δ ARO4 G226S aro10 Δ, pdc5 Δ, pdc6 Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3↑, at4CL3↑, cotal1↑) and IMX198 (□) (aro3 Δ ARO4 G226S aro10 Δ, pdc5 Δ, pdc6 Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,coCHS3↑, at4CL3↑, cotal1↑). Cultures were grown in shake flasks on synthetic medium containing 20 g·l-1 glucose. All cultures were performed in duplicate. Error bar denotes deviation of the mean.

Expression of a tyrosine ammonia lyase leads to increased naringenin synthesis

The engineered naringenin production pathway directs carbon via the aromatic amino acid phenylalanine to coumaric acid through phenylalanine ammonia lyase, trans-cinnamate 4-monooxygenase and subsequently to naringenin through the 4-coumaric acid-CoA ligase and chalcone synthase reactions (Figure1). Deregulation of aromatic amino acid synthesis in S. cerevisiae by eliminating feed-back inhibition of DAHP synthase not only increases intracellular phenylalanine levels, but also the intracellular concentration of tyrosine[40]. Since S. cerevisiae cannot interconvert tyrosine and phenylalanine, this increased intracellular availability of tyrosine could not be exploited for naringenin production. In other organisms, deamination of tyrosine by tyrosine ammonia lyase does provide an alternative route to coumaric acid. This has previously been shown to increase naringenin production by E coli when fed with either glucose or tyrosine[18, 20, 47]. To investigate the added benefit of the direct conversion of tyrosine into coumaric acid in S. cerevisiae IMX197, a cotal1 gene from Rhodobacter capsulatus was expressed. This R. capsulatus gene was chosen based on its known high catalytic efficiency for tyrosine[48]. Integration of this yeast codon-optimized cotal1 gene resulted in increased naringenin production in the resulting strain IMX198 (aro3 Δ ARO4 G226S aro10 Δ, pdc5 Δ, pdc6 Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,coCHS3↑, at4CL3↑, cotal1↑) (Figure5B), to a titer of 200 μM in shake-flask cultures. Additionally, the phloretic acid concentrations in this strain again increased compared to the parental strain IMX197 (623,9 ± 11,4 μM and 396,1 ± 6,9 respectively), while the coumaric acid concentrations in cultures of the two strains remained the same.

Naringenin production in controlled aerobic batch cultures

To further characterize strain IMX198 (aro3 Δ ARO4 G226S aro10 Δ, pdc5 Δ, pdc6 Δ, atPAL1↑, coC4H↑, coCPR1↑, atCHI1↑, atCHS3,coCHS3↑, at4CL3↑, cotal1↑) under controlled conditions, this strain was cultured in a 2L batch bioreactor with 20 g·l-1 glucose at pH 5.0 (Figure6A, B). When S. cerevisiae is grown aerobically in batch cultures on glucose, alcoholic fermentation is the predominant mode of glucose metabolism[49, 50] and is characterized by a diauxic growth profile. During the glucose consumption phase the specific growth rate of IMX198 was 0.2 h-1, which is approximately 50% of the specific growth rate of the reference strain S. cerevisiae CEN.PK113-7D under the same conditions[51]. Besides the expected formation of ethanol, a specific naringenin production rate of 12.545 ± 0.333 μmol.g-1CDWh-1 was obtained. After complete consumption of glucose, a naringenin titer 148.06 ± 5.67 μM at a naringenin yield on glucose of 0.002 ± 0.000 (g·g-1) was obtained. When all the glucose was consumed, ethanol, acetate and glycerol that were produced during the first phase were subsequently consumed (Figure6A). During this reconsumption phase, naringenin titers increased to 414.63 ± 1.60 μM, indicating that most naringenin is produced during this second phase (Figure6B). However, it must be taken into account that during this reconsumption phase, naringenin production is also facilitated by the presence of available coumarate that was previously formed during the glucose consumption phase and also by the higher amount of biomass, compared to the glucose phase. When only calculating the product yield over the total glucose and ethanol consumption phase, a yield of 0.006 ± 0.000 (g·g-1) was obtained, which is approximately triple of what is obtained during solely the glucose consumption phase.
https://static-content.springer.com/image/art%3A10.1186%2F1475-2859-11-155/MediaObjects/12934_2012_Article_781_Fig6_HTML.jpg
Figure 6

Product formation of the naringenin-producing strain S. cerevisiae IMX198 in bioreactors. Growth and extracellular metabolite formation were studied in aerobic and pH controlled (pH 5.0) batch cultures of IMX198 on glucose. The results shown are from a single representative experiment. A) Concentrations of glucose (), ethanol (), acetate (□), glycerol (■) and optical density (OD660) (Δ). B) Concentrations of naringenin (▲), coumaric acid (♦) and phloretic acid ().

In summary, the naringenin production of IMX198 in a well aerated and pH controlled batch bioreactor cultivation required only 38 hours and yielded more than twice the naringenin titer that was achieved in equivalent shake flasks cultivations (200 vs. 415 μM).

Discussion

The present study, for the first time, demonstrates de novo production of the key flavonoid intermediate naringenin, from glucose, by an engineered S. cerevisiae strain. Combined expression of the product pathway, codon optimization, improvement of precursor supply and reduction of byproduct formation led to concentrations of over 400 μM naringenin in aerobic, glucose-grown batch cultures. These concentrations are over 4-fold higher than reported in a previous study on de novo biosynthesis of naringenin by an engineered E. coli strain[18].

In previous studies on microbial biotransformation of aromatic precursors to naringenin, the core biosynthetic pathway for naringenin synthesis (phenylalanine and/or tyrosine ammonia lyase, trans-cinnamate 4-monooxygenase, P450 cytochrome reductase, 4-coumaric acid-CoA ligase, chalcone synthase) was, in most cases, reconstituted by co-expressing genes from different plants and/or bacteria[1719, 27, 28]. Only one of those studies used naringenin biosynthetic genes originating from a single donor organism (A. thaliana), which resulted in the first demonstration of de novo naringenin biosynthesis in E. coli[20]. In line with this observation, co-expression of a set of biosynthetic genes from A. thaliana, that were selected based on correlation of in planta expression profiles, was sufficient for low-level de novo production of naringenin in S. cerevisiae (Figure3). The observations that heterologous co-expression of genes originating from single plant species is beneficial for naringenin production, support the hypothesis that maintaining donor specific protein-protein interactions can improve pathway kinetics and specificity[52]. Besides using genes from a single plant species, also the recently developed synthetic protein scaffolds may be used to explore the impact of the physical association of pathway enzymes[5355]. The potential of this approach is illustrated by the recent demonstration that protein scaffolding in yeast led to a five-fold increase of the titer of the stilbenoid resveratrol[55]. The optimized supply of precursors and decreased production of aromatic byproducts make the yeast strains developed in this study an interesting models to explore the impact of pathway topology on flavonoid production in a heterologous host.

The improvement of naringenin production in this study followed a typical metabolic engineering approach consisting of repeated cycles of design, construction, and analysis[56]. After introduction of the pathway-encoding genes and observation of low levels of naringenin, alleviation of feedback inhibition of aromatic amino acid biosynthesis by phenylalanine and tyrosine further improved naringenin production. As seen previously[57]. removal of feedback inhibition also resulted in an overflow of phenylalanine into the Ehrlich pathway[40] and, consequently, production of phenylethanol as a major byproduct. Restricting the activity of the Ehrlich pathway by deletion of three decarboxylase-encoding genes strongly reduced formation of phenylethanol and had a strong positive impact on naringenin production (Figure3). Precursor supply was further improved by expression of a tyrosine ammonia lyase, which introduced a second route towards coumarate besides phenylalanine ammonia lyase. The subsequent introduction of two additional copies of the A. thaliana CHS3 gene further increased the flux from coumaric acid to naringenin (Figure5B). This observation may be attributed to the catalytic properties of CHS3. This enzyme combines 4 substrates (coumaroyl CoA and three copies of malonyl CoA) to synthesise naringenin, and it is well known to exhibit a very low catalytic activity. Although originating from a different plant, the catalytic constant kcat of the chalcone synthase CHS3 from alfalfa is 2 min-1,[58], which is roughly 25-fold lower than the kcat of the protein preceding the chalcone synthase in the metabolic pathway, the coumaroyl-CoA ligase (4CL3) from A.thaliana which has a kcat of 50 min-1[23].

In the naringenin-producing S. cerevisiae strains, the by-product phloretic acid was produced in concentrations similar to those achieved for naringenin. Phloretic acid production was only produced in cultures that both accumulated high concentrations of coumaric acid and expressed A. thaliana 4CL3. Phloretic acid has previously been detected in engineered S. cerevisiae, but the mechanism of its formation remains unknown[59]. Even though it was observed that increasing the flux to naringenin by increasing the copy number of the A. thaliana CHS3 gene, further analysis of the mechanism of phloretic acid production by these engineered S. cerevisiae, e.g. via systematic gene deletion studies, provides a good opportunity to further improve naringenin production. Also improvement of the provision of malonyl-CoA, a key precursor for both naringenin production as well as for lipid production, might further improve naringenin production. It has previously been shown that overexpression of acetyl-CoA carboxylase, which carboxylates acetyl-CoA to malonyl-CoA, leads to increased production of fatty acids in the yeast Hansenula polymorpha[60] and of biodiesel or 6-methylsalicylic acid production in S. cerevisiae[61, 62]. Moreover, malonyl-CoA synthase in E. coli combined with malonate supplementation, improved biotransformation of coumaric acid into flavonoids[63]. Increasing the availability of malonyl-CoA therefore provides an interesting target for further research.

Demonstration of de novo production of naringenin in S. cerevisiae represents an important step towards commercial production of plant-derived flavonoids from glucose. Even though commercial production of naringenin for direct consumer applications will require substantial further improvement of product yields and titers. Furthermore, this engineered naringenin-producing yeast strain provides an attractive platform for expression of plant-derived pathways that convert naringenin into other flavonoids, either by reconstruction of individual plant pathways or by combinatorial approaches based on the powerful yeast-based methods for in vivo recombination[29, 64]. Availability of a robust yeast platform for such studies paves the way for identification, characterization and eventual industrialization of many plant-derived flavonoids with pharma- and nutraceutical properties.

Conclusions

Here we provide the first report on the production of the flavonoid naringenin from glucose as sole carbon source using the yeast S. cerevisiae, and the optimization of the strain to reach a naringenin titer of over 400 μM. Several obstacles were overcome when expressing a foreign pathway in S. cerevisiae, which include the choice of enzymes, the expression systems the optimization of flux and the inhibition of competing pathways. This study successfully describes the reduction of the phenyl ethanol side product formation. All were done to increase the final titer of flavonoid production.

This strain may be used as an ideal platform for production of a wide variety of flavonoids originating from naringenin, thereby facilitating the production, characterization and possible implementation of such compounds as pharmaceutical and neutraceutical.

Methods

Strains and maintenance

All strains used in this study (Table2) are derived from the S. cerevisiae CEN.PK strain family background[65, 66]. Stock cultures were grown at 30°C in 500 ml shake flask cultures containing 100 ml synthetic medium (according to[67]) with 20 g·l-1 glucose and appropriate growth factors to supplement the specific auxotrophic requirements of the strains[68]. After overnight growth, a final concentration of 20% glycerol was added and 1 ml aliquots were stored at -80°C.
Table 2

Strains and plasmids used in this study

Strain

Genotype

Source

CEN.PK2-1c

MATalpha ura3 his3 leu2 trp1 MAL2-8cSUC2

[69]

CEN.PK717.5A

MATa aro3 Δ ARO4 G226S

This study

IMU011

MATalpha ura3 his3 leu2 trp1 MAL2-8cSUC2, pUDE172, PYK2::pUDI065

This study

IMK328

MATalpha ura3 his3 leu2 trp1 MAL2-8cSUC2 Δaro3(46,1065)::loxP ARO4 G226S

This study

IMK389

IMK328 pdc6Δ(-6,-2)::loxP pdc5Δ(-6,-2)::loxP

This study

IMK393

IMK389 aro10Δ(-6,-2)::loxP

This study

IMX183

IMK328 pUDE172, PYK2::pUDI065

This study

IMX185

IMK389 pUDE172, PYK2::pUDI065

This study

IMX106

IMK393 pUDE172, PYK2::pUDI065

This study

IMX193

IMK393 pUDE172, PYK2::pUDI065, TRP1::pUDI069

This study

IMX197

IMK393 pUDE172, PYK2::pUDI065, pUDE188

This study

IMX198

IMK393 pUDE172, PYK2::pUDI065, TRP1::pUDI069, pUDE188

This study

Plasmids

  

pUG6

PCR template for loxp-KanMX4-loxp cassette

[70]

pSH47

Centromeric plasmid, URA3, PGAL1-Cre- TCYC1

[70]

pAG416GPD-ccdB

Centromeric plasmid, URA3, PTDH3-ccdb-TCYC1, Addgene plasmid 14148

[71]

pAG306GPD-ccdB

Integration plasmid, URA3, PTDH3-ccdb-TCYC1, Addgene plasmid 14140

[71]

pAG305GPD-ccdB

Integration plasmid, LEU2, PTDH3-ccdb-TCYC1, Addgene plasmid 14138

[71]

pAG304GPD-ccdB

Integration plasmid, TRP1, PTDH3-ccdb-TCYC1, Addgene plasmid 14136

[71]

pAG423GPD-ccdB

2 μm ori, HIS3, PTDH3-ccdb-TCYC1, Addgene plasmid 14150

[71]

p426TEF

2 μm ori, URA3, PTEF –TTEF

[72]

p426TPI

2 μm ori, URA3 PTPI –TADH

[73]

pUNI.U21762

Cloning vector, Ka, atCHS3 cDNA#

[74]

pUNI.U11924

Cloning vector, Ka, atCHI1 cDNA

[74]

pDZ.4CL3

Cloning vector, Bla, at4CL3 cDNA

[75]

pUD168

pMA cloning vector, Bla, A-PTDH3-atPAL1-TCYC1-B*

This study

pUD169

pMA cloning vector, Bla, B-PPGI-coCPR1-TPGI-C

This study

pUD170

pMA cloning vector, Bla, C-PTPI-coC4H-TADH-D

This study

pUD175

pMA cloning vector, Bla, E-PPGK-coCHS3-TPGK-F

This study

pUDE185

2 μm ori, HIS3, PTDH3-coCHS3-TCYC1

This study

pUDE186

2 μm ori, URA3, PTEF-coCHS3-TTEF

This study

pUDE188

2 μm ori, HIS3, PTDH3-coCHS3-TCYC1, PTEF-coCHS3-TTEF

This study

pUDE103

2 μm ori, URA3, PTEF-at4CL3-TTEF

This study

pUDE172

Centromeric plasmid, URA3, PTDH3-atPAL1-TCYC1, PTPI-coC4H-TADH, PPGI-coCPR1-TPGI

This study

pUDI053

Integration plasmid, LEU2, PTDH3-atCHI1-TCYC1

This study

pUDI060

Integration plasmid, LEU2, PTDH3-atCHI1-TCYC1, PTEF-at4CL3-TTEF

This study

pUDI061

Integration plasmid, LEU2, PTDH3-atCHI1-TCYC1, PTPI-atCHS3-TADH, PTEF-at4CL3-TTEF

This study

pUDI065

Integration plasmid, LEU2, PTDH3-atCHI1-TCYC1, PTPI-atCHS3-TADH, PTEF-at4CL3-TTEF, PYK2(1-710)

This study

pUDI069

Integration plasmid, TRP1, PTDH3-coTAL1-TCYC1

This study

* co: codon optimized, #; at: A. thaliana.

Introduction of chromosomal gene alteration

Plasmids and oligonucleotide primers used in this study are listed in Tables2 and3, respectively. The geneticin (G418) resistance cassette (KanMX4) was amplified by using the pUG6 vector as the template[70]. The KanMX4 PCR product was transformed to the appropriate host (Table2) using the lithium acetate method[76]. Synthetic medium agar plates containing 2% glucose and 200 mg·l-1 G418 (Sigma-Aldrich, Zwijdrecht, The Netherlands) and appropriate growth factors to supplement the auxotrophic requirements of the strains, were used to select for the presence of the KanMX4 gene. The KanMX4 marker was removed by the Cre loxP recombination system using the plasmid pSH47 as described previously[70].
Table 3

Oligonucleotide primers used in this study

Name

Sequence (5’-3’)

Primers for knockouts

 

ARO3 KO for

ATGTTCATTAAAAACGATCACGCCGGTGACAGGAAACGCTTGGAAGACTGCAGCTGAAGCTTCGTACG

ARO3 KO rev

CTATTTTTTCAAGGCCTTTCTTCTGTTTCTAACACCTTCTGCCAATAGCTGCATAGGCCACTAGTGGATCTG

PDC5 KO for

CATAATCAATCTCAAAGAGAACAACACAATACAATAACAAGAAGAACAAACAGCTGAAGCTTCGTACGC

PDC5 KO rev

AAAGTAAAAAAATACACAAACGTTGAATCATGAGTTTTATGTTAATTAGCGCATAGGCCACTAGTGGATCTG

PDC6 KO for

GGCGGCTGTTTGAAGCCATTCTATCTTAATCTTGTGCTATTGCAGTCCTCCAGCTGAAGCTTCGTACGC

PDC6 KO rev

GTAAGTTTTATTTGCAACAATAATTCGTTTGAGTACACTACTAATGGCGCATAGGCCACTAGTGGATCTG

ARO10 KO for

GATACTCAAAACAAGTTGACGCGACTTCTGTAAAGTTTATTTACAAGATAACAAAGAAACTCCCTTAAGCCAGCTGAAGCTTCGTACGC

ARO10 KO rev

GGGTTTTTTATGTGTTAATGAACAGAAAACGAACAATTGGTAGCAGTGTTTTATAATTGCGCCCACAAGTCATAGGCCACTAGTGGATCTG

ARO4 for

GACGCATTGTTAGCATTG

ARO4 rev

CGAATTGGCAGTGGTAGAG

Primers for cloning

 

FK7

GCGACTAGTATGTCTTCATCCAACGCC

FK8

GCGGTCGACTCAGTTCTCTTTGGCTAGTT

FK9

GCGACTAGTATGATCACTGCAGCTCTAC

FK10

GCGGTCGACTCAACAAAGCTTAGCTTTGAG

FK11

GAGGCCGGCGCAATTAACCCTCACTAAAG

FK12

TGGACTCCAACGTCAAAG

FK13

GCGACTAGTATGGTGATGGCTGGTGC

FK14

GCGCTCGAGTTAGAGAGGAACGCTGTGC

FK15

CGCGCAATTAACCCTCAC

FK16

GCGGAGCTCCACTATAGGGCGAATTGGG

FK29

GCGTCCCGGAATGCCAGAGTCCAGATTG

FK30

GCGTCCCGGAGATATCTTGCCCTTCAGAACCC

FK52

GCGACTAGTATGACCTTACAATCCCAAAC

FK53

GCGCTCGAGTTAGGCTGGAG

FK86

ATCATGAACTTGCGCTCAATTCCGCGCAGAAGGCAATGCTATAAGACCTACGTCCACGGATTGCGCCTAAGACCGGATAAAGCACCGCATAGGGTAATAACTG

FK87

ATCGGAAATTCGACCGTGTGCTAGTGCCTATTGATGATCTGGCGGAATGTCTGCCGTGCCATAGCCATGCCTTCACATATAGTGTAATACGGTTATCCACAGAATC

FK169

GCGACTAGTATGGTTATGGCTGGTGC

FK170

GCGCTCGAGAATTACAATGGAACAGAGTGC

FK171

GCGGAGCTCGCTGGAGCTCATAGCTTC

FK172

GCGGAGCTCGTACCCAGTATAGCGACC

Primers for verification

 

FK105

TCTTTCCTGCGTTATCCC

FK106

GGCATGTACGGGTTACAG

FK107

CGCGTGTACGCATGTAAC

FK108

TCCCGTTAGGAACATTGG

FK109

GCAAATGCCTGCAAATCG

FK110

AACGTGCAGATGGTGATG

FK111

CATTATTGAACAGCGTCCAAG

FK112

AGAACCGTGGATGATGTG

The introduction of a feed-back-insensitive DAHP synthase activity in S.cerevisiae was accomplished in two steps. Firstly, the ARO3 allele was deleted in strain CEN.PK 2.1C using primers ARO3 KO for and ARO3 KO rev as described above. In this resulting strain, the chromosomal ARO4 allele was replaced by the feed-back-insensitive ARO4 G226S allele[38, 39]. Genomic DNA of strain CEN.PK717.5A was prepared using the genomic DNA isolation kit (Zymo Research, Orange, CA), according to the manufacturer’s recommendations. The ARO4 G226S allele was then amplified using primers Aro4 for/Aro4 rev, and transformed. Selection was performed by screening for tyrosine tolerance by plating on selective agar plates containing synthetic medium, 2% glucose and 1 g·l-1 tyrosine. The resulting strain was named IMK328 (Table2).

Plasmid construction

Plasmid pUDE172 was obtained by transformation-associated recombination methods as described previously[29, 30, 77]. The backbone fragment harboring the yeast and bacterial origins of replication and markers was amplified from pAG416GPD-ccdB (Table2) using primers FK86/87. Genes were optimized for S. cerevisiae using the JCat algorithm[78]), based on the A. thaliana gene sequence. The DNA fragments were synthesized (Life Technologies, Bleiswijk, The Netherlands) harboring the yeast codon-optimized genes (codon, promoter, terminator and 2 unique 80 bp flanks that included an Eco RV sites at each end). The codon-optimized plant genes were flanked by different promoter and terminator combinations to prevent plasmid instability due to homologous recombination. Prior to transformation-associated recombination, the synthesized fragments were digested from the cloning vector using Eco RV and gel purified. The pAG416GPD-ccdB backbone fragment and the 3 gene cassettes atPAL1, coC4H and coCPR were mixed at a ratio of 1:2:2:2 (50ng :100ng; 100ng: 100ng) and transformed to CEN.PK2-1C. Colonies were selected on glucose synthetic medium plates in the absence of uracil and checked by multiplex PCR using the primers FK105 – FK112. After plasmid isolation, retransformation to E. coli and subsequent plasmid purification, the structure of the resulting plasmid pUDE172 (Genbank accession number: JX268037) was verified using restriction analysis.

The A. thaliana CHI1 gene was amplified using primers FK7/8 from plasmid pUNI.U21762 (Table2), and then agarose gel purified and digested using Spe I and Sal I. The restricted fragment was ligated into an Spe I and Sal I linearized pAG305-GPD-ccdB plasmid. The A. thaliana 4CL3 gene was amplified using primers FK9/10 from pDZ.4CL3, purified and then ligated into p426TEF (Table2) using Spe I and Sal I. Subsequently, the PTEF-at4CL3-TTEF cassette was amplified using primer FK11/12, purified and ligated into a Nae I digested and FastAP (Thermo Scientific, Etten-Leur, The Netherlands) treated pUDI053, yielding pUDI060. The A. thaliana CHS3 gene was amplified using primer FK13/14, purified and ligated into p426TPI (Table2) using Spe I and Xho I yielding pUDE101. The PTPI-atCHS3-TADH cassette was amplified using primer FK15/16, purified and digested by Sac I. The resulting fragment was ligated into Sac I digested and FastAP treated pUDI060 plasmid, yielding pUDI061. A fragment of the PYK2 gene (1-710bp) was amplified from the CEN.PK2-1C genome using primer FK29/30, and then purified and digested using Pfo I. This fragment was ligated into pUDI061 plasmid that had been digested by Pfo I and treated with alkaline phosphatase (FastAP, Thermo Scientific) yielding pUDI065 (Genbank accession number: JX268039). The resulting plasmid was linearized using Nco I and integrated into the chromosomal PYK2 locus.

The tal1 gene from R. capsulatus[48, 79], was codon optimized for yeast using JCat[78], and synthesized (GeneScript, Piscataway, NJ). The synthetic cotal1 was amplified using primers FK52/53 and ligated into pAG304-ccdB using Spe I and Xho I resulting in pUDI069 (Genbank accession number: JX268036). Integration of this plasmid into the trp1 locus was preceded by linearization using Eco RV.

A DNA fragment containing the custom-synthesized, yeast codon-optimized CHS3 sequence (Life Technologies, Bleiswijk, The Netherlands) was amplified using primers FK169/170 and ligated into both a pAG325GPD-ccdB and p426TEF (Table2) using Spe I and Xho I yielding respectively, pUDE185 and 186. Subsequently, the PTEF-coCHS3-TTEF cassette was amplified using primers FK 171/172 and restricted using Sac I. This cassette was ligated into a Sac I digested and FastAP (Thermo Scientific) treated pUDE185. Plasmids were screened for opposing gene orientation by restriction analysis, yielding pUDE188 (Genbank accession number: JX268038).

Molecular biology procedures

PCR amplification was performed using Phusion® Hot Start II High Fidelity Polymerase (Thermo Scientific) according to the manufacturer's instructions in a Biometra TGradient Thermocycler (Biometra, Gottingen, Germany). Agarose gel separation was performed using 1% (w/v) agarose (Sigma-Aldrich, Zwijdrecht, The Netherlands) gel in 1×TAE (40 mM Tris-acetate pH 8.0 and 1 mM EDTA). Isolation of agarose trapped fragments was performed using Zymoclean Gel DNA Recovery kit (Zymo Research, Orange, CA). Restriction endonucleases, DNA ligases and FastAP (Thermo Scientific) were used according to the manufacturer's instructions. Transformation and amplification of plasmids were performed in E. coli DH5α electrocompetent cells (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. Plasmids were isolated from E. coli with the Sigma GenElute Plasmid Miniprep Kit (Sigma-Aldrich). Alternatively, plasmid purification was performed using a Zymoprep™™ Yeast Plasmid Miniprep (Zymo Research). Sequencing of constructs was performed by Baseclear BV (Baseclear, Leiden, The Netherlands).

Cultivation and media

E. coli was grown at 37°C in Luria Broth medium containing the appropriate antibiotic, 100 μg·l-1 ampicillin or 50 μg·l-1 kanamycin (Sigma-Aldrich, Zwijndrecht, The Netherlands). S. cerevisiae was grown at 30°C in medium containing demineralized water, 20 g·l-1 glucose, 5 g·l-1 (NH4)2SO4, 3 g·l-1 KH2PO4, 0.5 g·l-1 MgSO4.7H2O, vitamins and trace elements[67]. The pH of the medium was set to 6.0 using KOH. Appropriate growth factors to supplement the auxotrophic requirements of the strains were added at 150 mg·l-1 for uracil, 500 mg·l-1 for leucine, 75 mg·l-1 for tryptophan and 125 mg·l-1 for histidine[68]. Unless otherwise stated, cultures were grown in 500 ml shake flasks with 100 ml medium by adding 1 ml frozen stock culture and incubating at 30°C in an Innova incubator shaker (New Brunswick Scientific, Edison, NJ) set at 200 rpm. Controlled aerobic batch cultures were grown at 30°C in 2l bioreactors (Applikon, Schiedam, The Netherlands), using a working volume of 1l. In the bioreactor experiments, the (NH4)2SO4 concentration in the synthetic medium was increased to 10 g·l-1 to avoid nitrogen depletion towards the end of the culture. Emulsion C antifoam at 0.05% (w/v) (Sigma-Aldrich, Zwijndrecht, The Netherlands) was added separately after autoclaving. The pH was maintained at pH 5.0 by automatic addition of either 2 M KOH or 2 M H2SO4, the stirrer speed was fixed at 800 rpm and the aeration rate was set at 500 ml·min-1.

Analytical methods

Optical density was measured at 660 nm using a Libra S11 spectrophotometer (Biochrom, Cambridge, UK). Biomass dry weights were determined by filtration of 10 ml culture over dry, preweighed 0.45 μm nitrocellulose filters (Gelman Laboratory, Ann Arbor, USA). After removal of the medium, the filter was washed twice with demineralized water and dried for 20 min using a microwave set at 350 W. Glucose, ethanol, glycerol and acetate were analyzed using an Aminex HPX-87H ion exchange column (BioRad, Veenendaal, The Netherlands) operated at 60°C with 5 mM H2SO4 as mobile phase at a flow rate of 0.6 ml·min-1. For measurement of phenylpropanoic compounds, culture samples were diluted with an equal volume of 100% ethanol. After vigorous mixing, cells were spun down at 13000 rpm in a Sorval Heraeus #3328 rotor for 2 min. The supernatant was analysed using HPLC. For phenylpropanoic compounds, measurements were made using a Waters 2695 separation module and a Waters 996 photodiode array detector. Phloretic acid, coumaric acid, phenyl ethanol, cinnamic acid and naringenin were measured respectively at 275, 309, 214, 277, 289 nm using an Agilent Zorbax SB-C18 Column (4.6 × 5.0, 3.5 micron) operating at 30oC. A gradient of acetonitrile and 20 mM KH2PO4 (pH 2) with 1% acetonitrile was used as eluent, at a flow rate of 1 ml·min-1, increasing from 0 to 10% acetonitrile in 6 min followed by an increase to 40% acetonitrile until 23 min. From 23 min to 27 min, the flow was set to 100% KH2PO4. Naringenin, coumaric acid, cinnamic acid, phloretic acid and phenylethanol standards were obtained from Sigma Aldrich (Sigma-Aldrich, Zwijndrecht, The Netherlands).

For analysis of carbon dioxide production in bioreactor cultures, the off-gas was first cooled in a condenser (2°C) and dried with a Perma Pure Dryer (Permapure, Toms River, NJ). CO2 concentrations in the off-gas were then measured with a NGA 2000 Rosemount gas analyzer (Rosemount Analytical, Orrville,OH).

Identification of phenylpropanoid intermediates and naringenin was performed using liquid chromatography coupled to both photodiode-array detection and accurate mass quadrupole time-of-flight mass spectrometry (LC–PDA-QTOF MS) was performed using a Waters Alliance 2795 HPLC connected to a Waters 2996 PDA detector and subsequently a QTOF Ultima V4.00.00 mass spectrometer (Waters, MS technologies, UK) operating in negative ionization mode, an analytical column (Luna 3 μ C18/2 100A; 2.0 × 150 mm) attached to a C18 pre-column (2.0 × 4 mm; AJO-4286; both Phenomenex, USA). Eluents A (ultrapure water:formic acid (1000:1, v/v) ) and B (acetonitrile : formic acid (1000 : 1, v/v) ) were used at 0.19 ml.min. The gradient started at 5% B and increased linearly to 35% B in 45 min, after which the column was washed and equilibrated for 15 min before the next injection. The injection volume was 5 μl. Leucine enkephalin ([M-H]- = 554.2620) was used as a lock mass for on-line accurate mass correction[80]. Data were recorded using MassLynx 4.0 software (Waters).

Co-expression correlation analysis of naringenin biosynthetic genes in A. thaliana

Expression correlation analysis was performed using the BAR Expression Angler[35]. Each candidate gene listed in Table1 was used as “bait”, and expression was correlated over a set of 392 micro-array experiments generated using the ATH1 Affymetrix Whole Genome GeneChip, from the Nottingham Arabidopsis Stock Centre's microarray database[34]. Correlations between expression levels of genes from Table1 were recorded when the Pearson correlation coefficient was above 0.6.

Declarations

Acknowledgements

We thank Marinka Almering, Erik de Hulster, Bart Oud and Victor Guadalupe Medina for their technical assistance. The research group of JTP is part of the Kluyver Centre for Genomics of Industrial Fermentation (http://www.kluyvercentre.nl/pro1/general/home.asp/); the research group of RDH is part of the Centre for BioSYstems Genomics, both of which are supported by The Netherlands Genomics Initiative (http://www.genomics.nl/). FK, JTP and J-MD. were also supported by the “Platform Green Synthetic Biology” programme (http://www.pgsb.nl/) funded by NGI. FK, J-MD, AJAM were also financially supported by NGI through the Horizon breakthrough programme (project 93519011). JB, AH and RDH acknowledge support by the EU 7th Framework ATHENA project (FP7-KBBE-2009-3-245121-ATHENA). This research was also supported by BIOFLAVOUR, COST Action FA0907 (http://www.bioflavour.insa-toulouse.fr).

Authors’ Affiliations

(1)
Department of Biotechnology, Delft University of Technology
(2)
Platform for Green Synthetic Biology
(3)
Kluyver Centre for Genomics of Industrial Fermentation
(4)
Plant Research International (PRI)
(5)
Centre for Biosystems Genomics

References

  1. Xu P, Ranganathan S, Fowler ZL, Maranas CD, Koffas MA: Genome-scale metabolic network modeling results in minimal interventions that cooperatively force carbon flux towards malonyl-CoA. Metab Eng. 2011, 13: 578-587. 10.1016/j.ymben.2011.06.008View Article
  2. Facchini PJ, Bohlmann J, Covello PS, De Luca V, Mahadevan R, Page JE, Ro DK, Sensen CW, Storms R, Martin VJ: Synthetic biosystems for the production of high-value plant metabolites. Trends Biotechnol. 2012, 30: 127-131. 10.1016/j.tibtech.2011.10.001View Article
  3. Xiao ZP, Peng ZY, Peng MJ, Yan WB, Ouyang YZ, Zhu HL: Flavonoids Health Benefits and Their Molecular Mechanism. Mini-Rev Med Chem. 2011, 11: 169-177. 10.2174/138955711794519546View Article
  4. Wedick NM, Pan A, Cassidy A, Rimm EB, Sampson L, Rosner B, Willett W, Hu FB, Sun Q, van Dam RM: Dietary flavonoid intakes and risk of type 2 diabetes in US men and women. Am J Clin Nutr. 2012, 95: 925-933. 10.3945/ajcn.111.028894View Article
  5. Knekt P, Kumpulainen J, Jarvinen R, Rissanen H, Heliovaara M, Reunanen A, Hakulinen T, Aromaa A: Flavonoid intake and risk of chronic diseases. Am J Clin Nutr. 2002, 76: 560-568.
  6. Wang B, Zhang X: Inhibitory effects of Broccolini leaf flavonoids on human cancer cells. Scanning. 2012, 34: 1-5. 10.1002/sca.20278View Article
  7. McCullough ML, Peterson JJ, Patel R, Jacques PF, Shah R, Dwyer JT: Flavonoid intake and cardiovascular disease mortality in a prospective cohort of US adults. Am J Clin Nutr. 2012, 95: 454-464. 10.3945/ajcn.111.016634View Article
  8. Martin C, Butelli E, Petroni K, Tonelli C: How can research on plants contribute to promoting human health?. Plant Cell. 2011, 23: 1685-1699. 10.1105/tpc.111.083279View Article
  9. Zava DT, Duwe G: Estrogenic and antiproliferative properties of genistein and other flavonoids in human breast cancer cells in vitro. Nutr Cancer. 1997, 27: 31-40. 10.1080/01635589709514498View Article
  10. Greenwald P: Clinical trials in cancer prevention: current results and perspectives for the future. J Nutr. 2004, 134: 3507S-3512S.
  11. Hou DX, Fujii M, Terahara N, Yoshimoto M: Molecular mechanisms behind the chemopreventive effects of anthocyanidins. J Biomed Biotechnol. 2004, 2004: 321-325. 10.1155/S1110724304403040View Article
  12. Allister EM, Borradaile NM, Edwards JY, Huff MW: Inhibition of microsomal triglyceride transfer protein expression and apolipoprotein B100 secretion by the citrus flavonoid naringenin and by insulin involves activation of the mitogen-activated protein kinase pathway in hepatocytes. Diabetes. 2005, 54: 1676-1683. 10.2337/diabetes.54.6.1676View Article
  13. Fowler ZL, Koffas MA: Biosynthesis and biotechnological production of flavanones: current state and perspectives. Appl Microbiol Biotechnol. 2009, 83: 799-808. 10.1007/s00253-009-2039-zView Article
  14. Wang Y, Chen S, Yu O: Metabolic engineering of flavonoids in plants and microorganisms. Appl Microbiol Biotechnol. 2011, 91: 949-956. 10.1007/s00253-011-3449-2View Article
  15. Nijveldt RJ, van Nood E, van Hoorn DE, Boelens PG, van Norren K, van Leeuwen PA: Flavonoids: a review of probable mechanisms of action and potential applications. Am J Clin Nutr. 2001, 74: 418-425.
  16. Limem I, Guedonc E, Hehn A, Bourgaud F, Ghedira LC, Engasser JM, Ghoul M: Production of phenylpropanoid compounds by recombinant microorganisms expressing plant-specific biosynthesis genes. Process Biochem. 2008, 43: 463-479. 10.1016/j.procbio.2008.02.001. 10.1016/j.procbio.2008.02.001View Article
  17. Leonard E, Lim KH, Saw PN, Koffas MA: Engineering central metabolic pathways for high-level flavonoid production in Escherichia coli. Appl Environ Microbiol. 2007, 73: 3877-3886. 10.1128/AEM.00200-07View Article
  18. Santos CNS, Koffas M, Stephanopoulos G: Optimization of a heterologous pathway for the production of flavonoids from glucose. Metab Eng. 2011, 13: 392-400. 10.1016/j.ymben.2011.02.002View Article
  19. Trantas E, Panopoulos N, Ververidis F: Metabolic engineering of the complete pathway leading to heterologous biosynthesis of various flavonoids and stilbenoids in Saccharomyces cerevisiae. Metab Eng. 2009, 11: 355-366. 10.1016/j.ymben.2009.07.004View Article
  20. Watts KT, Lee PC, Schmidt-Dannert C: Exploring recombinant flavonoid biosynthesis in metabolically engineered Escherichia coli. Chem Biochem. 2004, 5: 500-507.
  21. Koopmann E, Hahlbrock K: Differentially regulated NADPH: cytochrome P450 oxidoreductases in parsley. Proc Natl Acad Sci USA. 1997, 94: 14954-14959. 10.1073/pnas.94.26.14954View Article
  22. Winkel-Shirley B: Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001, 126: 485-493. 10.1104/pp.126.2.485View Article
  23. Costa MA, Bedgar DL, Moinuddin SG, Kim KW, Cardenas CL, Cochrane FC, Shockey JM, Helms GL, Amakura Y, Takahashi H, et al: Characterization in vitro and in vivo of the putative multigene 4-coumarate:CoA ligase network in Arabidopsis: syringyl lignin and sinapate/sinapyl alcohol derivative formation. Phytochem. 2005, 66: 2072-2091. 10.1016/j.phytochem.2005.06.022. 10.1016/j.phytochem.2005.06.022View Article
  24. Burbulis IE, Winkel-Shirley B: Interactions among enzymes of the Arabidopsis flavonoid biosynthetic pathway. Proc Natl Acad Sci USA. 1999, 96: 12929-12934. 10.1073/pnas.96.22.12929View Article
  25. Winkel BSJ: Metabolic channeling in plants. Annu Rev Plant Biol. 2004, 55: 85-107. 10.1146/annurev.arplant.55.031903.141714View Article
  26. Yu O, Shi J, Hession AO, Maxwell CA, McGonigle B, Odell JT: Metabolic engineering to increase isoflavone biosynthesis in soybean seed. Phytochem. 2003, 63: 753-763. 10.1016/S0031-9422(03)00345-5. 10.1016/S0031-9422(03)00345-5View Article
  27. Jiang H, Wood KV, Morgan JA: Metabolic engineering of the phenylpropanoid pathway in Saccharomyces cerevisiae. Appl Environ Microbiol. 2005, 71: 2962-2969. 10.1128/AEM.71.6.2962-2969.2005View Article
  28. Yan Y, Chemler J, Huang L, Martens S, Koffas MA: Metabolic engineering of anthocyanin biosynthesis in Escherichia coli. Appl Environ Microbiol. 2005, 71: 3617-3623. 10.1128/AEM.71.7.3617-3623.2005View Article
  29. Gibson DG, Benders GA, Andrews-Pfannkoch C, Denisova EA, Baden-Tillson H, Zaveri J, Stockwell TB, Brownley A, Thomas DW, Algire MA, et al: Complete chemical synthesis, assembly, and cloning of a Mycoplasma genitalium genome. Science. 2008, 319: 1215-1220. 10.1126/science.1151721View Article
  30. Shao Z, Zhao H: DNA assembler, an in vivo genetic method for rapid construction of biochemical pathways. Nucleic Acids Res. 2009, 37: e16-, 10.1093/nar/gkn724View Article
  31. Jiang H, Morgan JA: Optimization of an in vivo plant P450 monooxygenase system in Saccharomyces cerevisiae. Biotechnol Bioeng. 2004, 85: 130-137. 10.1002/bit.10867View Article
  32. Zhang YH: Substrate channeling and enzyme complexes for biotechnological applications. Biotechnol Adv. 2011, 29: 715-725. 10.1016/j.biotechadv.2011.05.020View Article
  33. Raes J, Rohde A, Christensen JH, Van de Peer Y, Boerjan W: Genome-wide characterization of the lignification toolbox in Arabidopsis. Plant Physiol. 2003, 133: 1051-1071. 10.1104/pp.103.026484View Article
  34. Craigon DJ, James N, Okyere J, Higgins J, Jotham J, May S: NASCArrays: a repository for microarray data generated by NASC's transcriptomics service. Nucleic Acids Res. 2004, 32: D575-D577. 10.1093/nar/gkh133View Article
  35. Toufighi K, Brady SM, Austin R, Ly E, Provart NJ: The botany array resource: e-northerns, expression angling, and promoter analyses. Plant J. 2005, 43: 153-163. 10.1111/j.1365-313X.2005.02437.xView Article
  36. Saito K, Hirai MY, Yonekura-Sakakibara K: Decoding genes with coexpression networks and metabolomics - 'majority report by precogs'. Trends Plant Sci. 2008, 13: 36-43. 10.1016/j.tplants.2007.10.006View Article
  37. Ehlting J, Buttner D, Wang Q, Douglas CJ, Somssich IE, Kombrink E: Three 4-coumarate:coenzyme A ligases in Arabidopsis thaliana represent two evolutionarily divergent classes in angiosperms. Plant J. 1999, 19: 9-20. 10.1046/j.1365-313X.1999.00491.xView Article
  38. Schnappauf G, Hartmann M, Kunzler M, Braus GH: The two 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase isoenzymes from Saccharomyces cerevisiae show different kinetic modes of inhibition. Arch Microbiol. 1998, 169: 517-524. 10.1007/s002030050605View Article
  39. Hartmann M, Schneider TR, Pfeil A, Heinrich G, Lipscomb WN, Braus GH: Evolution of feedback-inhibited beta/alpha barrel isoenzymes by gene duplication and a single mutation. Proc Natl Acad Sci USA. 2003, 100: 862-867. 10.1073/pnas.0337566100View Article
  40. Luttik MA, Vuralhan Z, Suir E, Braus GH, Pronk JT, Daran JM: Alleviation of feedback inhibition in Saccharomyces cerevisiae aromatic amino acid biosynthesis: quantification of metabolic impact. Metab Eng. 2008, 10: 141-153. 10.1016/j.ymben.2008.02.002View Article
  41. Vuralhan Z, Luttik MA, Tai SL, Boer VM, Morais MA, Schipper D, Almering MJ, Kotter P, Dickinson JR, Daran JM, Pronk JT: Physiological characterization of the ARO10-dependent, broad-substrate-specificity 2-oxo acid decarboxylase activity of Saccharomyces cerevisiae. Appl Environ Microbiol. 2005, 71: 3276-3284. 10.1128/AEM.71.6.3276-3284.2005View Article
  42. Vuralhan Z, Morais MA, Tai SL, Piper MD, Pronk JT: Identification and characterization of phenylpyruvate decarboxylase genes in Saccharomyces cerevisiae. Appl Environ Microbiol. 2003, 69: 4534-4541. 10.1128/AEM.69.8.4534-4541.2003View Article
  43. Hazelwood LA, Daran JM, van Maris AJ, Pronk JT, Dickinson JR: The Ehrlich pathway for fusel alcohol production: a century of research on Saccharomyces cerevisiae metabolism. Appl Environ Microbiol. 2008, 74: 2259-2266. 10.1128/AEM.02625-07View Article
  44. Romagnoli G, Luttik MA, Kotter P, Pronk JT, Daran JM: Substrate specificity of thiamine-pyrophosphate-dependent 2-oxo-acid decarboxylases in Saccharomyces cerevisiae. Appl Environ Microbiol. 78: 7538-7548.
  45. Flikweert MT, de Swaaf M, van Dijken JP, Pronk JT: Growth requirements of pyruvate-decarboxylase-negative Saccharomyces cerevisiae. FEMS Microbiol Lett. 1999, 174: 73-79. 10.1111/j.1574-6968.1999.tb13551.xView Article
  46. Zuurbier KWM, Fung SY, Scheffer JJC, Verpoorte R: Assay of chalcon synthase activity by high-performance liquid-chromatography. Phytochem. 1993, 34: 1225-1229. 10.1016/0031-9422(91)80005-L. 10.1016/0031-9422(91)80005-LView Article
  47. Hwang EI, Kaneko M, Ohnishi Y, Horinouchi S: Production of plant-specific flavanones by Escherichia coli containing an artificial gene cluster. Appl Environ Microbiol. 2003, 69: 2699-2706. 10.1128/AEM.69.5.2699-2706.2003View Article
  48. Kyndt JA, Meyer TE, Cusanovich MA, Van Beeumen JJ: Characterization of a bacterial tyrosine ammonia lyase, a biosynthetic enzyme for the photoactive yellow protein. FEBS Lett. 2002, 512: 240-244. 10.1016/S0014-5793(02)02272-XView Article
  49. De Deken RH: The Crabtree effect: a regulatory system in yeast. J Gen Microbiol. 1966, 44: 149-156.View Article
  50. Van Hoek P, Van Dijken JP, Pronk JT: Effect of specific growth rate on fermentative capacity of baker's yeast. Appl Environ Microbiol. 1998, 64: 4226-4233.
  51. van Maris AJA, Bakker BM, Brandt M, Boorsma A, de Mattos MJT, Grivell LA, Pronk JT, Blom J: Modulating the distribution of fluxes among respiration and fermentation by overexpression of HAP4 in Saccharomyces cerevisiae. FEMS Yeast Research. 2001, 1: 139-149.View Article
  52. Winkel-Shirley B: Evidence for enzyme complexes in the phenylpropanoid and flavonoid pathways. Physiol Plant. 1999, 107: 142-149. 10.1034/j.1399-3054.1999.100119.x. 10.1034/j.1399-3054.1999.100119.xView Article
  53. Dueber JE, Wu GC, Malmirchegini GR, Moon TS, Petzold CJ, Ullal AV, Prather KL, Keasling JD: Synthetic protein scaffolds provide modular control over metabolic flux. Nat Biotechnol. 2009, 27: 753-759. 10.1038/nbt.1557View Article
  54. Moon TS, Dueber JE, Shiue E, Prather KLJ: Use of modular, synthetic scaffolds for improved production of glucaric acid in engineered E. coli. Metab Eng. 2010, 12: 298-305. 10.1016/j.ymben.2010.01.003View Article
  55. Wang Y, Yu O: Synthetic scaffolds increased resveratrol biosynthesis in engineered yeast cells. J Biotechnol. 2011, 157: 258-260.View Article
  56. Kim IK, Roldao A, Siewers V, Nielsen J: A systems-level approach for metabolic engineering of yeast cell factories. FEMS Yeast Res. 2012, 12: 228-248. 10.1111/j.1567-1364.2011.00779.xView Article
  57. Santos CNS: Combinatorial search strategies for the metabolic engineering of microorganisms. 2010, Cambridge, MA: Massachusetts Institute of Technology, PhD thesis,2010.
  58. Jez JM, Ferrer JL, Bowman ME, Dixon RA, Noel JP: Dissection of malonyl-coenzyme A decarboxylation from polyketide formation in the reaction mechanism of a plant polyketide synthase. Biochemistry. 2000, 39: 890-902. 10.1021/bi991489fView Article
  59. Beekwilder J, Wolswinkel R, Jonker H, Hall R, de Vos CH, Bovy A: Production of resveratrol in recombinant microorganisms. Appl Environ Microbiol. 2006, 72: 5670-5672. 10.1128/AEM.00609-06View Article
  60. Ruenwai R, Cheevadhanarak S, Laoteng K: Overexpression of acetyl-CoA carboxylase gene of Mucor rouxii enhanced fatty acid content in Hansenula polymorpha. Mol Biotechnol. 2009, 42: 327-332. 10.1007/s12033-009-9155-yView Article
  61. Shi S, Valle-Rodriguez JO, Khoomrung S, Siewers V, Nielsen J: Functional expression and characterization of five wax ester synthases in Saccharomyces cerevisiae and their utility for biodiesel production. Biotechnol Biofuels. 2012, 5: 7-, 10.1186/PREACCEPT-1932279820621895View Article
  62. Wattanachaisaereekul S, Lantz AE, Nielsen ML, Nielsen J: Production of the polyketide 6-MSA in yeast engineered for increased malonyl-CoA supply. Metab Eng. 2008, 10: 246-254. 10.1016/j.ymben.2008.04.005View Article
  63. Leonard E, Yan Y, Fowler ZL, Li Z, Lim CG, Lim KH, Koffas MA: Strain improvement of recombinant Escherichia coli for efficient production of plant flavonoids. Mol Pharm. 2008, 5: 257-265. 10.1021/mp7001472View Article
  64. Naesby M, Nielsen SV, Nielsen CA, Green T, Tange TO, Simon E, Knechtle P, Hansson A, Schwab MS, Titiz O, et al: Yeast artificial chromosomes employed for random assembly of biosynthetic pathways and production of diverse compounds in Saccharomyces cerevisiae. Microb Cell Fact. 2009, 8: 45-, 10.1186/1475-2859-8-45View Article
  65. van Dijken JP, Bauer J, Brambilla L, Duboc P, Francois JM, Gancedo C, Giuseppin ML, Heijnen JJ, Hoare M, Lange HC, et al: An interlaboratory comparison of physiological and genetic properties of four Saccharomyces cerevisiae strains. Enzyme Microb Technol. 2000, 26: 706-714. 10.1016/S0141-0229(00)00162-9View Article
  66. Nijkamp JF, van den Broek M, Datema E, de Kok S, Bosman L, Luttik MA, Daran-Lapujade P, Vongsangnak W, Nielsen J, Heijne WH, et al: De novo sequencing, assembly and analysis of the genome of the laboratory strain Saccharomyces cerevisiae CEN.PK113–7D, a model for modern industrial biotechnology. Microb Cell Fact. 2012, 11: 36-, 10.1186/1475-2859-11-36View Article
  67. Verduyn C, Postma E, Scheffers WA, Vandijken JP: Physiology of Saccharomyces cerevisiae in Anaerobic Glucose-Limited Chemostat Cultures. J Gen Microbiol. 1990, 136: 395-403. 10.1099/00221287-136-3-395View Article
  68. Pronk JT: Auxotrophic yeast strains in fundamental and applied research. Appl Environ Microbiol. 2002, 68: 2095-2100. 10.1128/AEM.68.5.2095-2100.2002View Article
  69. Entian KD, Kotter P: 25 Yeast genetic strain and plasmid collections. Method Microbiol. 2007, 36: 629-666.View Article
  70. Gueldener U, Heinisch J, Koehler GJ, Voss D, Hegemann JH: A second set of loxP marker cassettes for Cre-mediated multiple gene knockouts in budding yeast. Nucleic Acids Res. 2002, 30: e23-, 10.1093/nar/30.6.e23View Article
  71. Alberti S, Gitler AD, Lindquist S: A suite of Gateway (R) cloning vectors for high-throughput genetic analysis in Saccharomyces cerevisiae. Yeast. 2007, 24: 913-919. 10.1002/yea.1502View Article
  72. Mumberg D, Muller R, Funk M: Yeast Vectors for the Controlled Expression of Heterologous Proteins in Different Genetic Backgrounds. Gene. 1995, 156: 119-122. 10.1016/0378-1119(95)00037-7View Article
  73. de Kok S, Yilmaz D, Suir E, Pronk JT, Daran JM, van Maris AJ: Increasing free-energy (ATP) conservation in maltose-grown Saccharomyces cerevisiae by expression of a heterologous maltose phosphorylase. Metab Eng. 2011, 13: 518-526. 10.1016/j.ymben.2011.06.001View Article
  74. Yamada K, Lim J, Dale JM, Chen H, Shinn P, Palm CJ, Southwick AM, Wu HC, Kim C, Nguyen M, et al: Empirical analysis of transcriptional activity in the Arabidopsis genome. Science. 2003, 302: 842-846. 10.1126/science.1088305View Article
  75. Seki M, Narusaka M, Kamiya A, Ishida J, Satou M, Sakurai T, Nakajima M, Enju A, Akiyama K, Oono Y, et al: Functional annotation of a full-length Arabidopsis cDNA collection. Science. 2002, 296: 141-145. 10.1126/science.1071006View Article
  76. Gietz RD, Schiestl RH: High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method. Nat Protoc. 2007, 2: 31-34. 10.1038/nprot.2007.13View Article
  77. Noskov VN, Koriabine M, Solomon G, Randolph M, Barrett JC, Leem SH, Stubbs L, Kouprina N, Larionov V: Defining the minimal length of sequence homology required for selective gene isolation by TAR cloning. Nucleic Acids Res. 2001, 29: E32-, 10.1093/nar/29.6.e32View Article
  78. Grote A, Hiller K, Scheer M, Munch R, Nortemann B, Hempel DC, Jahn D: JCat: a novel tool to adapt codon usage of a target gene to its potential expression host. Nucleic Acids Res. 2005, 33: W526-W531. 10.1093/nar/gki376View Article
  79. Watts KT, Mijts BN, Lee PC, Manning AJ, Schmidt-Dannert C: Discovery of a substrate selectivity switch in tyrosine ammonia-lyase, a member of the aromatic amino acid lyase family. Chem Biol. 2006, 13: 1317-1326. 10.1016/j.chembiol.2006.10.008View Article
  80. De Vos RC, Moco S, Lommen A, Keurentjes JJ, Bino RJ, Hall RD: Untargeted large-scale plant metabolomics using liquid chromatography coupled to mass spectrometry. Nat Protoc. 2007, 2: 778-791. 10.1038/nprot.2007.95View Article

Copyright

© Koopman et al.; licensee BioMed Central Ltd. 2012

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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.