Pattern formation in auxin flux

Fluctuations in auxin levels depend upon synchronicity of cell divisions in a one-dimensional model of auxin transport

Auxin is a well-studied plant hormone, the spatial distribution of which remains incompletely understood. Here, we investigate the effects of cell growth and divisions on the dynamics of auxin patterning, using a combination of mathematical modelling and experimental observations. In contrast to most prior work, models are not designed or tuned with the aim to produce a specific auxin pattern. Instead, we use well-established techniques from dynamical systems theory to uncover and classify ranges of auxin patterns as exhaustively as possible as parameters are varied. Previous work using these techniques has shown how a multitude of stable auxin patterns may coexist, each attainable from a specific ensemble of initial conditions. When a key parameter spans a range of values, these steady patterns form a geometric curve with successive folds, often nicknamed a snaking diagram. As we introduce growth and cell division into a one-dimensional model of auxin distribution, we observe new behaviour which can be explained in terms of this diagram. Cell growth changes the shape of the snaking diagram, and this corresponds in turn to deformations in the patterns of auxin distribution. As divisions occur this can lead to abrupt creation or annihilation of auxin peaks. We term this phenomenon 'snake-jumping'. Under rhythmic cell divisions, we show how this can lead to stable oscillations of auxin. We also show that this requires a high level of synchronisation between cell divisions. Using 18 hour time-lapse imaging of the auxin reporter DII:Venus in roots of Arabidopsis thaliana, we show auxin fluctuates greatly, both in terms of amplitude and periodicity, consistent with the snake-jumping events observed with non-synchronised cell divisions. Periodic signals downstream of the auxin signalling pathway have previously been recorded in plant roots. The present work shows that auxin alone is unlikely to play the role of a pacemaker in this context.

Auxin-driven patterning with unidirectional fluxes

The plant hormone auxin plays an essential role in the patterning of plant structures. Biological hypotheses supported by computational models suggest that auxin may fulfil this role by regulating its own transport, but the plausibility of previously proposed models has been questioned. We applied the notion of unidirectional fluxes and the formalism of Petri nets to show that the key modes of auxin-driven patterning-the formation of convergence points and the formation of canals-can be implemented by biochemically plausible networks, with the fluxes measured by dedicated tally molecules or by efflux and influx carriers themselves. Common elements of these networks include a positive feedback of auxin efflux on the allocation of membrane-bound auxin efflux carriers (PIN proteins), and a modulation of this allocation by auxin in the extracellular space. Auxin concentration in the extracellular space is the only information exchanged by the cells. Canalization patterns are produced when auxin efflux and influx act antagonistically: an increase in auxin influx or concentration in the extracellular space decreases the abundance of efflux carriers in the adjacent segment of the membrane. In contrast, convergence points emerge in networks in which auxin efflux and influx act synergistically. A change in a single reaction rate may result in a dynamic switch between these modes, suggesting plausible molecular implementations of coordinated patterning of organ initials and vascular strands predicted by the dual polarization theory.

Self-organization of plant vascular systems: claims and counter-claims about the flux-based auxin transport model

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.