BackgroundCassava is a well-known starchy root crop utilized for food, feed and biofuel production. However, the comprehension underlying the process of starch production in cassava is not yet available.ResultsIn this work, we exploited the recently released genome information and utilized the post-genomic approaches to reconstruct the metabolic pathway of starch biosynthesis in cassava using multiple plant templates. The quality of pathway reconstruction was assured by the employed parsimonious reconstruction framework and the collective validation steps. Our reconstructed pathway is presented in the form of an informative map, which describes all important information of the pathway, and an interactive map, which facilitates the integration of omics data into the metabolic pathway. Additionally, to demonstrate the advantage of the reconstructed pathways beyond just the schematic presentation, the pathway could be used for incorporating the gene expression data obtained from various developmental stages of cassava roots. Our results exhibited the distinct activities of the starch biosynthesis pathway in different stages of root development at the transcriptional level whereby the activity of the pathway is higher toward the development of mature storage roots.ConclusionsTo expand its applications, the interactive map of the reconstructed starch biosynthesis pathway is available for download at the SBI group’s website (http://sbi.pdti.kmutt.ac.th/?page_id=33). This work is considered a big step in the quantitative modeling pipeline aiming to investigate the dynamic regulation of starch biosynthesis in cassava roots.

Analysis of metabolic flux was used for system level assessment of carbon partitioning in Kasetsart 50 (KU50) and Hanatee (HN) cassava cultivars to understand the metabolic routes for their distinct phenotypes. First, the c onstraint- b ased metabolic m odel of cassava storage r oots, rMeCBM, was developed based on the carbon assimilation pathway of cassava. Following the subcellular compartmentalization and curation to ensure full network connectivity and reflect the complexity of eukaryotic cells, cultivar specific data on sucrose uptake and biomass synthesis were input, and rMeCBM model was used to simulate storage root growth in KU50 and HN. Results showed that rMeCBM-KU50 and rMeCBM-HN models well imitated the storage root growth. The flux-sum analysis revealed that both cultivars utilized different metabolic precursors to produce energy in plastid. More carbon flux was invested in the syntheses of carbohydrates and amino acids in KU50 than in HN. Also, KU50 utilized less flux for respiration and less energy to synthesize one gram of dry storage root. These results may disclose metabolic potential of KU50 underlying its higher storage root and starch yield over HN. Moreover, sensitivity analysis indicated the robustness of rMeCBM model. The knowledge gained might be useful for identifying engineering targets for cassava yield improvement.

The existing genome-scale metabolic model of carbon metabolism in cassava storage roots, rMeCBM, has proven particularly resourceful in exploring the metabolic basis for the phenotypic differences between high and low-yield cassava cultivars. However, experimental validation of predicted metabolic fluxes by carbon labeling is quite challenging. Here, we incorporated gene expression data of developing storage roots into the basic flux-balance model to minimize infeasible metabolic fluxes, denoted as rMeCBMx, thereby improving the plausibility of the simulation and predictive power. Three different conceptual algorithms, GIMME, E-Flux, and HPCOF were evaluated. The rMeCBMx-HPCOF model outperformed others in predicting carbon fluxes in the metabolism of storage roots and, in particular, was highly consistent with transcriptome of high-yield cultivars. The flux prediction was improved through the oxidative pentose phosphate pathway in cytosol, as has been reported in various studies on root metabolism, but hardly captured by simple FBA models. Moreover, the presence of fluxes through cytosolic glycolysis and alanine biosynthesis pathways were predicted with high consistency with gene expression levels. This study sheds light on the importance of prediction power in the modeling of complex plant metabolism. Integration of multi-omics data would further help mitigate the ill-posed problem of constraint-based modeling, allowing more realistic simulation.

Increasing cassava production could mitigate one of the global food insecurity challenges by providing a sustainable food source. To improve the yield potential, physiological strategies (i.e., the photosynthetic efficiency, source-to-sink carbon partitioning, and intracellular carbon metabolism) can be applied in breeding to screen for superior genotypes. However, the influences of source-to-sink carbon partitioning and carbon metabolism on the storage root development of cassava are relatively little understood. We hypothesized that carbon partitioning and utilization vary modulating the distinctive storage root yields of high and low-yielding cassava varieties, represented in this study by varieties Kasetsart 50 (KU50) and Hanatee (HN), respectively. Plant growth, photosynthesis measurements, soluble sugars, and starch contents of individual tissues were analyzed at different developmental stages. Also, the diurnal patterns of starch accumulation and degradation in leaves were investigated through iodine staining. Despite a comparable photosynthetic rate, KU50 grew better and yielded greater storage roots than HN. Interestingly, both varieties differed in their carbon partitioning strategies. KU50 had a high photosynthetic capacity and was better efficient in converting photoassimilates to carbon substrates and allocating them to sink organs for their growth. In contrast, HN utilized the photoassimilates at a high metabolic cost, in terms of respiration, and inefficiently allocated carbon to stems rather than storage roots. These results highlighted that carbon assimilation and allocation are genetic potential characteristics of individual varieties, which in effect determine plant growth and storage root yield of cassava. The knowledge gained from this study sheds light on potential strategies for developing new high-yielding genotypes in cassava breeding programs.

The constraint-based rMeCBM-KU50 model of cassava storage root growth was analyzed to evaluate its sensitivity, with respect to reaction flux distribution and storage root growth rate, to changes in model inputted data and constraints, including sucrose uptake rate-related data—photosynthetic rate, total leaf area, total photosynthetic rate, storage root dry weight, and biomass function-related data. These mainly varied within ±90% of the model default values, although exceptions were made for the carbohydrate (−90% to 8%) and starch (−90% to 9%) contents. The results indicated that the predicted storage root growth rate was highly affected by specific sucrose uptake rates through the total photosynthetic rate and storage root dry weight variations; whereas the carbon flux distribution, direction and partitioning inclusive, was more sensitive to the variation in biomass content, particularly the carbohydrate content. This study showed that the specific sucrose uptake rate based on the total photosynthetic rate, storage root dry weight, and carbohydrate content were critical to the constraint-based metabolic modeling and deepened our understanding of the input–output relationship—specifically regarding the rMeCBM-KU50 model—providing a valuable platform for the modeling of plant metabolic systems, especially long-growing crops.

No abstract