Flavonoids are secondary metabolites present in all terrestrial plants. The flavonoid pathway has been extensively studied, and many of the involved genes and metabolites have been described in the literature. Despite this extensive knowledge, the functioning of the pathway in vivo is still poorly understood. Here, we study the flavonoid pathway using both experiments and mathematical models. We measured flavonoid metabolite dynamics in two tissues, hypocotyls and cotyledons, during tomato seedling development. Interestingly, the same backbone of interactions leads to very different accumulation patterns in the different tissues. Initially, we developed a mathematical model with constant enzyme concentrations that described the metabolic networks separately in both tissues. This model was unable to fit the measured flavonoid dynamics in the hypocotyls, even if we allowed unrealistic parameter values. This suggested us to investigate the effect of transcript abundance on flavonoid accumulation. We found that the expression of candidate flavonoid genes varies considerably with time. Variation in transcript abundance results in enzymatic variation, which could have a large effect on metabolite accumulation. Candidate transcript abundance was included in the mathematical model as representative for enzyme concentration. We fitted the resulting model to the flavonoid dynamics in the cotyledons, and tested it by applying it to the data from hypocotyls. When transcript abundance is included, we are indeed able to explain flavonoid dynamics in both tissues. Importantly, this is possible under the biologically relevant restriction that the enzymatic properties estimated by the model are conserved between the tissues.
References
[1]
Koes R, Quattrocchio F, Mol J (1994) The avonoid biosynthetic pathway in plants: function and evolution. BioEssays 16: 123–132.
[2]
Prasad S, Phromnoi K, Yadav VR, Chaturvedi MM, Aggarwal BB (2010) Targeting inammatory pathways by avonoids for prevention and treatment of cancer. Planta Med 76: 1044–1063.
[3]
Hodgson JM, Croft KD (2010) Tea avonoids and cardiovascular health. Mol Aspects Med 31: 495–502.
[4]
Grotewold E (2006) The genetics and biochemistry of oral pigments. Annu Rev Plant Biol 57: 761–780.
[5]
Taylor LP, Grotewold E (2005) Flavonoids as developmental regulators. Curr Opin Plant Biol 8: 317–323.
[6]
Bieza K, Lois R (2001) An arabidopsis mutant tolerant to lethal ultraviolet-b levels shows consti-tutively elevated accumulation of avonoids and other phenolics. Plant Physiol 126: 1105–1115.
[7]
Li J, Ou-Lee TM, Raba R, Amundson RG, Last RL (1993) Arabidopsis avonoid mutants are hypersensitive to uv-b irradiation. Plant Cell 5: 171–179.
[8]
Treutter D (2006) Significance of avonoids in plant resistance: a review. Environmental Chemistry Letters 4: 147–157.
[9]
Harborne JB, Williams CA (2000) Advances in avonoid research since 1992. Phytochemistry 55: 481–504.
[10]
Olsen KM, Hehn A, Jugdé H, Slimestad R, Larbat R, et al. (2010) Identification and characteri-sation of cyp75a31, a new avonoid 3'5'-hydroxylase, isolated from solanum lycopersicum. BMC Plant Biol 10: 21.
[11]
Zifkin M, Jin A, Ozga JA, Zaharia LI, Schernthaner JP, et al. (2012) Gene expression and metabo-lite profiling of developing highbush blueberry fruit indicates transcriptional regulation of avonoid metabolism and activation of abscisic acid metabolism. Plant Physiol 158: 200–224.
[12]
Castellarin SD, Matthews MA, Gaspero GD, Gambetta GA (2007) Water deficits accelerate ripen-ing and induce changes in gene expression regulating avonoid biosynthesis in grape berries. Planta 227: 101–112.
[13]
Bovy A, Gomez-Roldan V, Hall R (2010) Strategies to optimize the avonoid content of tomato fruit. Recent Advances in Polyphenol Research 2.
[14]
Bovy A, de Vos R, Kemper M, Schijlen E, Pertejo MA, et al. (2002) High-avonol tomatoes resulting from the heterologous expression of the maize transcription factor genes lc and c1. Plant Cell 14: 2509–2526.
[15]
Kubasek WL, Shirley BW, McKillop A, Goodman HM, Briggs W, et al. (1992) Regulation of avonoid biosynthetic genes in germinating arabidopsis seedlings. Plant Cell 4: 1229–1236.
[16]
Pelletier MK, Burbulis IE, Winkel-Shirley B (1999) Disruption of specific avonoid genes enhances the accumulation of avonoid enzymes and end-products in arabidopsis seedlings. Plant Mol Biol 40: 45–54.
[17]
Trabelsi N, d'Estaintot BL, Sigaud G, Gallois B, Chaudire J (2011) Kinetic and binding equilib-rium studies of dihydroavonol 4-reductase from Vitis vinifera and its unusually strong substrate inhibition. Journa of Biophysical Chemistry 2: 332–344.
[18]
Rodriguez A, Infante D (2011) Characterization in silico of avonoids biosynthesis in Theobroma cacao l. Network Biology 1(1): 34–45.
[19]
Clark ST, Verwoerd WS (2011) A systems approach to identifying correlated gene targets for the loss of colour pigmentation in plants. BMC Bioinformatics 12: 343.
[20]
Kitano H (2002) Systems biology: a brief overview. Science 295: 1662–1664.
[21]
Covert MW, Schilling CH, Palsson B (2001) Regulation of gene expression in ux balance models of metabolism. J Theor Biol 213: 73–88.
[22]
Nakabayashi R, Yamazaki M, Saito K (2010) A polyhedral approach for understanding avonoid biosynthesis in arabidopsis. N Biotechnol 27: 829–836.
[23]
Gutenkunst RN, Waterfall JJ, Casey FP, Brown KS, Myers CR, et al. (2007) Universally sloppy parameter sensitivities in systems biology models. PLoS Comput Biol 3: 1871–1878.
[24]
Apri M, de Gee M, Molenaar J (2012) Complexity reduction preserving dynamical behavior of biochemical networks. J Theor Biol 304C: 16–26.
[25]
Crosby KC, Pietraszewska-Bogiel A, Gadella TWJ, Winkel BSJ (2011) Frster resonance energy transfer demonstrates a avonoid metabolon in living plant cells that displays competitive inter-actions between enzymes. FEBS Lett 585: 2193–2198.
[26]
Preu？ A, Stracke R, Weisshaar B, Hillebrecht A (2009) Arabidopsis thaliana expresses a second functional avonol synthase. FEBS Lett 583: 1981–1986.
[27]
Groenenboom MAC, Astola L, Gomez-Roldan M, van Eeuwijk F, et al. Function prediction of flavonoid glycosyltransferase in tomato. In preparation.
[28]
Astola L, Groenenboom MAC, Gomez-Roldan M, Bovy A, Hall R, et al. Tree graphs and identifiable parameter estimation in metabolic networks. In preparation.
[29]
Gomez-Roldan M, Groenenboom MAC, Astola L, de Vos R, van Eeuwijk F, et al. Tomato seedlings as a valuable model to understand dynamic metabolic networks. In preparation.
[30]
Bolstad BM, Irizarry RA, Astrand M, Speed TP (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19: 185–193.
[31]
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, et al. (2003) Summaries of affymetrix genechip probe level data. Nucleic Acids Res 31: e15.
[32]
Sato S, Tabata S, Hirakawa H, Asamizu E, Shirasawa K, et al. (2012) The tomato genome sequence provides insights into eshy fruit evolution. Nature 485: 635–641.
[33]
Butelli E, Titta L, Giorgio M, Mock HP, Matros A, et al. (2008) Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors. Nat Biotechnol 26: 1301–1308.
