Comparative constraint-based modelling of fruit development across species highlights nitrogen metabolism in the growth-defence trade-off.

Although primary metabolism is well conserved across species, it is useful to explore the specificity of its network to assess the extent to which some pathways may contribute to particular outcomes. Constraint-based metabolic modelling is an established framework for predicting metabolic fluxes and phenotypes and helps to explore how the plant metabolic network delivers specific outcomes from temporal series. After describing the main physiological traits during fruit development, we confirmed the correlations between fruit relative growth rate (RGR), protein content and time to maturity. Then a constraint-based method is applied to a panel of eight fruit species with a knowledge-based metabolic model of heterotrophic cells describing a generic metabolic network of primary metabolism. The metabolic fluxes are estimated by constraining the model using a large set of metabolites and compounds quantified throughout fruit development. Multivariate analyses showed a clear common pattern of flux distribution during fruit development with differences between fast- and slow-growing fruits. Only the latter fruits mobilise the tricarboxylic acid cycle in addition to glycolysis, leading to a higher rate of respiration. More surprisingly, to balance nitrogen, the model suggests, on the one hand, nitrogen uptake by nitrate reductase to support a high RGR at early stages of cucumber and, on the other hand, the accumulation of alkaloids during ripening of pepper and eggplant. Finally, building virtual fruits by combining 12 biomass compounds shows that the growth-defence trade-off is supported mainly by cell wall synthesis for fast-growing fruits and by total polyphenols accumulation for slow-growing fruits.

Model-assisted analysis for tuning anthocyanin composition in grape berries.

Anthocyanin composition is responsible for the red colour of grape berries and wines, and contributes to their organoleptic quality. However, anthocyanin biosynthesis is under genetic, developmental and environmental regulation, making its targeted fine-tuning challenging. We constructed a mechanistic model to simulate the dynamics of anthocyanin composition throughout grape ripening in Vitis vinifera, employing a consensus anthocyanin biosynthesis pathway. The model was calibrated and validated using six datasets from eight cultivars and 37 growth conditions. Tuning the transformation and degradation parameters allowed us to accurately simulate the accumulation process of each individual anthocyanin under different environmental conditions. The model parameters were robust across environments for each genotype. The coefficients of determination (R2) for the simulated versus observed values for the six datasets ranged from 0.92 to 0.99, while the relative root mean square errors (RRMSEs) were between 16.8 and 42.1 %. The leave-one-out cross-validation for three datasets showed R2 values of 0.99, 0.96 and 0.91, and RRMSE values of 28.8, 32.9 and 26.4 %, respectively, suggesting a high prediction quality of the model. Model analysis showed that the anthocyanin profiles of diverse genotypes are relatively stable in response to parameter perturbations. Virtual experiments further suggested that targeted anthocyanin profiles may be reached by manipulating a minimum of three parameters, in a genotype-dependent manner. This model presents a promising methodology for characterizing the temporal progression of anthocyanin composition, while also offering a logical foundation for bioengineering endeavours focused on precisely adjusting the anthocyanin composition of grapes.