Tiago Moreira Machado

Understanding stress-induced metabolic reprogramming in crop plants can inform breeding strategies and support the development of stress-resilient varieties. Genome-scale metabolic modelling has shown promise in elucidating network-level responses to changing environments, yet as an optimality-based approach it relies on the definition of an objective function, which is far from trivial for non-optimal conditions. To address this uncertainty, we used a time-resolved, data-informed metabolic model of rice ( Oryza sativa L.) cold stress response as a test case, and explored two complementary approaches. We used sampling of the solution space combined with machine learning to identify reactions and pathways best characterizing the stress-induced metabolic shift, and used this information to perform Pareto analysis, placing growth and a stress-related objective in competition. This trade-off analysis identified key branch points in carbohydrate, amino acid, phenylpropanoid, nucleotide, and fatty acid biosynthesis, where resource reallocation towards stress-protection comes at the expense of growth. It further revealed differential flux modes across subcellular compartments and shifts in reducing equivalent provision as distinguishing features of the stress response. Together, these results provide a mechanistic understanding of the metabolic trade-offs and branch points governing cold stress response, and identify potential targets to optimize the cold response-growth trade-off in rice.

C4 photosynthesis evolved from the ancestral C3 pathway through coordinated leaf anatomical and metabolic reorganization that concentrates CO ₂ to reduce photorespiration. Quantitative understanding of these structure-function relationships remains limited. Here we used anatomy-aware metabolic modeling of a mesophyll-bundle sheath cell system to analyze the interdependence between leaf anatomy and photosynthetic metabolism on the C3-C4 spectrum. Our model faithfully recapitulates the transitory steps from C3 to C4 photosynthesis, reveals a crucial role for plasmodesmata in enabling the C3 to C4 transition, and points at potential pre-C2 metabolic states that provide benefits under conditions that promote elevated photorespiration. Incorporating bundle cell suberisation with our model predicts reduction of PSII activity and dominance of the NADP-ME C4 subtype in leaves with suberized bundle sheath cells and proposes a role for oxygen evolution at PSII as a potential driver for this mechanism. Varying bundle sheath leakage and photorespiratory conditions along the C3-C4 spectrum identify conditions under which C3-C4 intermediacy provides energetic benefits and underlines the notion of intermediate photosynthesis as a stable evolutionary state. Overall, our study sheds new light on the quantitative relationship between leaf anatomy and metabolism and its interaction with the environment and suggests targets for climate-adaptation in C3 plants.