Plants are highly plastic in their potential to adapt to changing environmental conditions. For example, they can selectively promote the relative growth of the root and the shoot in response to limiting supply of mineral nutrients and light, respectively, a phenomenon that is referred to as balanced growth or functional equilibrium. To gain insight into the regulatory network that controls this phenomenon, we took a systems biology approach that combines experimental work with mathematical modeling. We developed a mathematical model representing the activities of the root (nutrient and water uptake) and the shoot (photosynthesis), and their interactions through the exchange of the substrates sugar and phosphate (Pi). The model has been calibrated and validated with two independent experimental data sets obtained with Petunia hybrida. It involves a realistic environment with a day-and-night cycle, which necessitated the introduction of a transitory carbohydrate storage pool and an endogenous clock for coordination of metabolism with the environment. Our main goal was to grasp the dynamic adaptation of shoot:root ratio as a result of changes in light and Pi supply. The results of our study are in agreement with balanced growth hypothesis, suggesting that plants maintain a functional equilibrium between shoot and root activity based on differential growth of these two compartments. Furthermore, our results indicate that resource partitioning can be understood as the emergent property of many local physiological processes in the shoot and the root without explicit partitioning functions. Based on its encouraging predictive power, the model will be further developed as a tool to analyze resource partitioning in shoot and root crops.

The plant hormone auxin plays a central role in growth and morphogenesis. In shoot apical meristems, auxin flux is polarized through its interplay with PIN proteins. Concentration-based mathematical models of the flux can explain some aspects of phyllotaxis for the L1 surface layer, where auxin accumulation points act as sinks and develop into primordia. The picture differs in the interior of the meristem, where the primordia act as auxin sources, leading to the initiation of the vascular system. Self-organization of the auxin flux involves large numbers of molecules and is difficult to treat by intuitive reasoning alone; mathematical models are therefore vital to understand these phenomena. We consider a leading computational model based on the so-called flux hypothesis. This model has been criticized and extended in various ways. One of the basic counter-arguments is that simulations yield auxin concentrations inside canals that are lower than those seen experimentally. Contrary to what is claimed in the literature, we show that the model can lead to higher concentrations within canals for significant parameter regimes. We then study the model in the usual case where the response function Φ defining the model is quadratic and unbounded, and show that the steady state vascular patterns are formed of loopless directed trees. Moreover, we show that PIN concentrations can diverge in finite time, thus explaining why previous simulation studies introduced cut-offs which force the system to have bounded PIN concentrations. Hence, contrary to previous claims, extreme PIN concentrations are not due to numerical problems but are intrinsic to the model. On the other hand, we show that PIN concentrations remain bounded for bounded Φ, and simulations show that in this case, loops can emerge at steady state.

The plant hormone auxin is fundamental for plant growth, and its spatial distribution in plant tissues is critical for plant morphogenesis. We consider a leading model of the polar auxin flux, and study in full detail the stability of the possible equilibrium configurations. We show that the critical states of the auxin transport process are composed of basic building blocks, which are isolated in a background of auxin depleted cells, and are not geometrically regular in general. The same model was considered recently through a continuous limit and a coupling to the von Karman equations, to model the interplay of biochemistry and mechanics during plant growth. Our conclusions might be of interest in this setting, since, for example, we establish the existence of Lyapunov functions for the auxin flux, proving in this way the convergence of pure transport processes toward the set of critical configurations.