Tissue-wide integration of mechanical cues promotes effective auxin patterning

New plant organs form by local accumulation of auxin, which is transported by PIN proteins that localize following mechanical stresses. As auxin itself modifies tissue mechanics, a feedback loop between tissue mechanics and auxin patterning unfolds—yet the impact of tissue-wide mechanical coupling on auxin pattern emergence remains unclear. Here, we use a model composed of a vertex model for plant tissue mechanics and a compartment model for auxin transport to explore the collective mechanical response of the tissue to auxin patterns and how it feeds back onto auxin transport. We compare a model accounting for a tissue-wide mechanical integration to a model that regards cells as mechanically isolated. We show that tissue-wide mechanical coupling not only leads to more focused auxin spots via stress redistribution, but that it also mitigates the disruption to patterning when considering noise in the mechanical properties of each cell of the tissue. We find that this mechanism predicts that a local turgor increase correlates with auxin concentration, and yet auxin spots can exist regardless of the exact local turgor distribution.


Introduction
Formation of organs entails an effective coordination of local cell growth typically initiated by patterns of one or more morphogenic factors. Understanding how these patterns of morphogenic agents robustly emerge is fundamental for predicting organ morphogenesis. Plants organ formation is interesting from a physical perspective due to the strong mechanical coupling between plant cells, and the fact that growth is driven by changes in the mechanical properties of the cell wall and internal pressure [1][2][3][4][5]. Evidence indicates that the morphogenic factors such as the plant hormone auxin change the mechanics of the tissue [6,7], with implications for the shaping of organs [8,9]. Interestingly, the transporters of auxin respond to mechanical cues [10,11], leading to an intertwining of chemical and mechanical cues.
The phytohormone auxin, Indole-3-Acetic Acid, is the key morphogenic agent in plants. Auxin accumulation drives a wide range of plant developmental processes including, but A difference in auxin, therefore, induces a stress difference between the two compartments separating both cells. This stress difference is such that PIN binds preferentially in the cell with lower auxin concentration, increasing the flow of auxin into the cell with higher auxin concentration Fig. 2 Schematic difference between the tissue-wide mechanical model (left) and the uncoupled tissue approximation (right). In the tissue-wide mechanical model, turgor pressure, T , and stiffness determine the vertex positions that minimize mechanical energy. Wall strain and stress are then inferred from the mechanical configuration. In the approximation, we prescribe average wall stress,σ , with a static geometry. This approximation disregards the effect of stiffness variations on strain. The prescription of stress in the approximation renders the mechanical interaction to be only between nearest neighbours and uncoupled from all other cells. In the tissue-wide mechanical model, the mechanical state is a function of all cells in the tissue

Methods
In order to investigate the interaction between auxin cell wall softening and collective tissue mechanics, we use a vertex model to describe the mechanical behaviour of the tissue and a compartment model to express auxin concentration and transport between adjacent cells.
α is that same quantity when the cell is at rest. σ i j refers to the longitudinal stress acting on the compartment of the wall. (bottom right) Cell α has an auxin concentration a α which is expressed, degraded and transported, both passively and actively. The active component of auxin transport relies on the density of membrane-bound efflux auxin carriers facing a particular wall compartment, p i j 2.1 Geometrical set-up of the tissue The tissue is described by a tiling of two-dimensional space into M cells surrounded by their cell walls. Walls are represented as edges connecting two vertices each, positioned at Here, we reserve Latin indices for vertex numbering and Greek ones for cells. Each cell wall segment has two compartments, one facing each cell. Therefore, we represent each cell wall with two edges of opposite direction, one for each compartment. The position of tissue vertices fully define geometrical quantities such as cell areas, A α , cell perimeters, L α , wall lengths, l i j = l ji , and cell centroids, X α (Fig. 3 top left). To simplify notation significantly, we also define for each cell the cyclically ordered set of all vertices around that cell, V α , arranged counterclockwise (ccw). Hence, we use i∈V α to signify the sum over all vertices surrounding cell α with an arbitrary start, where i + 1 and i − 1 mean, respectively, the next and previous ccw vertex. Similarly, we introduce N α as the cyclically ordered (counterclockwise) set of all neighbouring regions around cell α, one for each edge of α (Fig. 3 top right).

Tissue mechanics-tissue-wide coupling
Vertex models are a widely employed theoretical approach to describe mechanics of epithelial tissues and morphogenesis [9,[48][49][50][51][52][53]. The essence of vertex models is that cell geometry within a tissue is given as the mechanical equilibrium of the tissue. In the case of plant cells, the shape of a cell is a competition between the turgor pressure, T α , all cells exert on each other and the cell's resistance to deformation with stiffness, E α . Strain acting on each cell will be described using the second moment of area of the corresponding cell in reference to its centroid, M α , whose components are where the primed coordinates represent the translation transformation, α , we define cell strain as the normalized difference between both matrices, and stress with σ α = E α ε α . Having described the tissue mechanically ( Fig. 3 bottom left), we define the energy for a single cell as the sum of work done by turgor pressure and elastic deformation energy, resulting in the tissue mechanical energy, Using this model, we obtain the shape of the tissue by minimizing H with respect to vertex positions. After minimizing (Eq. 5), we quantify the stress acting on each wall through the average strain acting on each cell given by (Eq. 4). Assuming that cell wall rest length is the same between two adjacent wall compartments then it follows that they are under the same longitudinal strain, which is, to first approximation, the average between the two cells surrounding them. Therefore, longitudinal average strain acting on a specific wall used here is where t αβ is a unit vector along the wall separating cell α and cell β. Note that this interpolation assumes a continuous strain field. Then the stresses acting on each compartment are by the constitutive equation of a linear elastic isotropic material with Poisson ratio ν = 0, Note that we are only considering the longitudinal components with regards to the cell wall, which means thatε αβ and σ αβ are scalar quantities. More details on the mechanical model used can be found in the supporting text. As argued in the supporting material, our choice of ν = 0 does not impact the qualitative behaviour studied here.

Tissue mechanics-uncoupled tissue approximation
To assess the impact of collective mechanical behaviour within a tissue on auxin pattern selforganization, we approximate the tissue-wide mechanical model to a static tissue geometry where we approximate the effects of turgor pressure of each individual cells in the static tissue by a constant average stressσ acting on it [10]. Again assuming that both wall compartments have the same rest length, we infer that the stress acting on a particular wall depends only on σ and the stiffness of the adjacent cells. Effectively, the average longitudinal strain acting on a wall surrounded by cells α and β would simply bē This way, instead of minimizing the full mechanical model (Eq. 5) given a set of turgor pressures T α and rest shape matrices M α we can, in the static tissue, immediately compute stress with Eq. 7 yielding, Interestingly, Eq. 9 is valid for ν = 0 as demonstrated in the supporting material. In order to compare the two models, we choose the value ofσ to be the same as the stress obtained through minimisation of (Eq. 5), for a given set of T α and M (0) α , with the constraint of the same end geometry.
Note that not only can this approximation be interpreted as the tissue being mechanically coupled only to the nearest neighbours, disregarding the rest of the tissue, (Fig. 2), but also as an analogous non-mechanical auxin concentration feedback model.

Auxin transport-compartment model
Compartment models for auxin transport are well adapted to the context of plant development, since the prerequisite of a boundary of a plant cell is particularly well defined by courtesy of the cell wall.
Although passive diffusion occurs across cell walls, the dominant players in auxin transport are membrane-bound carriers [22,24]. Namely, efflux transporters of the PIN family are important due to their anisotropic positioning around a cell [16], which leads to a net auxin flow from one cell to the next. Let a α denote an non-dimensional and normalized average auxin concentration inside cell α. Following the model by [10], which is similar to previous mathematical models [25,26,29], auxin evolves according to auxin metabolism in the cell, passive diffusion between cells and active transport across cell walls via PIN, where γ * is the auxin production rate, δ * is the auxin decay rate, W αβ = l αβ /A α , with K , P, and D as adjustable parameters. D is the passive permeability of plant cells, whereas P is permeability of the cell wall due to PIN-mediated transport of auxin, and K is the Michaelis-Menten constant for the efflux of auxin. More information on how this expression is derived can be found in the supporting text. Although this description ignores the auxin present within the extracellular domain and inside the cell wall, it has been shown that under physiological assumptions, this is a valid approximation [29]. The active transport term depends on the amount of bound PIN in each cell wall, where f αβ , β ∈ N α expresses the ratio between binding and unbinding rates of a particular wall ( Fig. 3 bottom right). Note that p αβ is different from wall to wall and from cell to cell. This means that in general, p αβ = p βα , or equivalently, p i j = p ji . This is consistent with the fact that there are two compartments to a cell wall shared by two adjacent cells. Expression (Eq. 11) is based on the assumption that cell walls around a particular cell compete for the same pool of PIN molecules and that the amount of PIN scales with cell perimeter. This competition has been shown to be important in the polarization of PIN [29]. Alternatively, one could also scale the amount of PIN with cell size or not scale it at all. In the former case, smaller cells would be slightly preferred for auxin accumulation, whereas in the latter, larger cells would be preferred instead. Since we want to study the impact of stress patterns on the tissue, we want to decouple it from this effect as much as possible, choosing instead to scale the amount of PIN with perimeter. The trivial fixed point of these dynamical equations is given by a α = μ * /δ * , ∀α, which also results in equal PIN density across all walls, provided turgor pressure T α and stiffness E α are the same across the tissue.
The feedback between tissue mechanics and auxin pattern unfolds as auxin transport affects tissue mechanics due to auxin, a α , controlling cell wall stiffness, E α , and in reverse tissue stress, σ α , affects auxin transport by regulating PIN binding rates, f αβ , as hypothesized by [10,11].

Mechanical regulation of PIN binding
According to the hypothesis presented by [10,11], mechanical cues up-regulate PIN binding. The distinction between whether these mechanical cues are strain or stress has been studied recently by [54], yet the exact nature remains unclear. Following the model presented by [10], we model the binding-unbinding ratio, f αβ , as being a power law on positive stress, where the stresses, σ αβ , follow from tissue mechanics after minimization of the full mechanical model (Eq. 5), or, in the averaged stress approximation, it is the stress load on that particular compartment given by (Eq. 9). Furthermore, n is the exponent of this power law, and η captures the coupling between stress and PIN. Effectively, this mechanical coupling to PIN parameter corresponds to the sensing and subsequent response to stress, loosely translating into how much resources the cell needs to spend for processing stress cues.

Auxin-mediated cell wall softening
Auxin affects the mechanical properties of a cell wall via methyl esterification of pectin [6,7], resulting in a decrease of the stiffness of the cell wall. We assume that all cell wall  Then we integrate auxin dynamics for a time step and update the stiffness of each cell. This process knocks the system out of the previous mechanical energy minimum, and it has to be relaxed again. Alternatively, we can shortcut energy minimization using the averaged stress approximation for a static tissue. This procedure is repeated until t = t max . The parameters r , wall loosening effect, and η, stress coupling, interface both models and are, therefore, of critical importance to the mechanism studied compartments surrounding cell α share the same stiffness, E α . To capture this effect, we model stiffness with a Hill function [10], where r ∈ [0, 1[ which we define as the cell wall loosening effect, m is the Hill exponent of this interaction, and E 0 is the stiffness of the cell walls when its auxin concentration is a α = 1. At low values of auxin, E α approaches the value (1 + r ) E 0 , whereas at high auxin concentration, E α approaches (1 − r ) E 0 . Given a distribution of auxin, we can compute the wall stiffness in (Eq. 5) from (Eq. 13), or the stress acting on a specific compartment in (Eq. 9) for the approximated model.

Integrating auxin transport and tissue mechanics
At each time step, Δt, starting from an auxin distribution, we compute the stiffness of each cell according to (Eq. 13). Then, with the input of all turgor pressures, we minimize (Eq. 5) to obtain tissue geometry and stresses acting on each wall. Auxin concentration in each cell will evolve according to (Eq. 10), where the active transport term will be regulated by stress according to (Eq. 12) via (Eq. 11). A new auxin distribution will result at the end of this iteration, and we will be ready to take another time step (Fig. 4). We repeat this process until t = t max .

Implementation
We implemented this model with C++ programming language, where we have used the Quad-Edge data structure for geometry and topology of the tissue [55], implemented in the library Quad-Edge [56]. In order to minimize the mechanical energy of the tissue, we have used a limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm (L-BFGS) [57,58], implemented in the library NLopt [59]. For solving the set of ODEs presented in the compartment model, we used the explicit embedded Runge-Kutta-Fehlberg method (often referred to as RKF45) implemented in the GNU Scientific Library (GSL) [60]. We wrapped the resulting classes into a python module with SWIG. For additional details regarding the parameters used for the simulations of the following section, consult Table S1 in the supporting material.

Observables
In order to quantify the existence of auxin patterns, we compute the difference between an emerging auxin concentration pattern and the trivial steady state of uniform auxin concentration pattern defined as a α = γ * /δ * , ∀α. To account for a large range of orders of magnitude of auxin concentration, we consider as an order parameter, where · M denotes an average over all cells within the tissue. This way, ϕ ≈ 0 means that there are no discernible patterns, whereas ϕ ≈ 1 implies prominent auxin patterning. The term δ 2 defines the sensitivity of this measure, such that an average deviation of δ yields ϕ ≈ 1/2 (for small δ). We will choose δ = 0.1, i. e. , a 10% deviation from the trivial steady state.
We also keep track of the average of auxin above basal levels in order to gauge the potential degree of modulation of auxin-mediated cell behaviour.
Furthermore, to characterize cells with regards to PIN localization we introduce the magnitude of the average PIN efflux direction, wheren ii+1 is the unit vector normal to the wall pointing outwards from α. Aside from a global measure of auxin patterning, it is also important to locally relate auxin to tissue mechanics. Namely, for auxin we are interested in auxin concentration, a α , and auxin local gradient, obtained by interpolation, where In fact, the quantity |∇a α | can be used as an indicator of whether there is an interface between auxin spots and the rest of the tissue. With regards to tissue mechanics, the local quantities we quantify are the isotropic component of stress, and the stress deviator tensor projected along the direction of the auxin gradient, where σ α = σ α − I P α , and I is the identity matrix. Therefore, P α is a measure if a cell is being compressed (P α < 0), or pulled apart (P α > 0), and D α translates into if a cell is more compressed along the auxin gradient than perpendicular to it (D α < 0), or vice-versa (D α > 0). Finally, to measure the disruption of an auxin pattern we approximate entropy by means of a Riemann sum, where Π (a) is the probability density function of auxin and Δa the partition size. Note that it is only meaningful to compare entropy measures obtained with the same partition size Δa. Here, the probability density function of auxin concentration is obtained by applying a kernel density estimation on the resulting tissue auxin values. Note that Π (a) is a continuous function. In order to infer it from simulation data, for each auxin value in the tissue, a α , we add Kernel functions K w (a), obeying where w is a smoothing parameter defining the width of the Kernel, this parameter is sometimes called bandwidth. This statistical tool is called kernel density estimation (KDE) [61]. We use the Epanechnikov kernel because it is bounded and we can force Π (a) = 0, a ≤ 0.

The tissue-wide mechanical model captures stress patterns after ablation
First we verify that the tissue-wide mechanical model captures the expected mechanical behaviour and auxin patterning when a cell is ablated. To model ablation, we set the stiffness of the ablated cell walls to E 0 = 0, block all auxin transport to and from it, block PIN transporters of adjacent cells from binding to the shared wall with the ablated cell, and, finally, we lower the turgor pressure to only 10% of the original value. This remnant of pressure represents the surface tension emerging from pressure of the inner layers of the shoot apical meristem acting on a curved surface, as required by the Young-Laplace equation. This is necessary since the model only simulates the epidermal layer in a plane. We observe that the region neighbouring the ablation site gets depleted of auxin due to PIN binding preferentially to the walls circumferentially aligned around the ablated cell in accordance with the stress principal directions (Fig. 5a). This stress pattern is in agreement with calculations performed by [41] in this setting and PIN aligns according to the ablation experiments in [10].
We also simulated different wound shapes. The resulting stress patterns are shown the supporting material. Stress directions align along the shape of the ablation wound.

Conditions for auxin patterns emergence
The uncoupled tissue approximation allows to analytically compute the conditions for spontaneous auxin pattern emergence in a general regular lattice (Fig. 6). Effectively, for a regular grid, the condition for pattern formation is, where M = nmr, W = 4/ √ 3 is a geometrical factor specific to the used grid, and p 0 = f (σ ) /(1 + f (σ )) (see supporting material for the linear stability analysis details). Equation 22 is the closed form of more general expressions presented by [10,29] tailored to our system and parameters.
To quantify the existence of auxin patterns in the model with tissue-wide stress patterning, we computed the order parameter ϕ defined in Eq. 14 for simulations with different values of wall loosening effect r and stress coupling η (Fig. 6). These two parameters are conceptually important since the former is the cause for stiffness inhomogeneity of the tissue, and the latter represents a plant cell's sensitivity to mechanical cues.
We observe a very good agreement between the conditions for pattern emergence (Eq. 22) analytically predicted in the case of the uncoupled tissue approximation and the transition of ϕ in the case of tissue-wide stress patterning (Fig. 6). This means that at the onset of patterns emergence the auxin concentrations are similar enough to make the assumption that the effect of turgor pressure is simply an isotropic stress across the entire tissue, validating the approximation near the transition. This observation is in agreement with the auxin pattern emergence mechanism hypothesis by [10] (Fig. 1). The agreement between the two models does not necessarily apply after patterns emerge. This poses the question of the role of mechanics in potentially enhancing or hindering auxin flows.

Global mechanical response reinforces PIN polarity
To understand the role of tissue-wide stress patterning on the emergence of PIN-driven auxin patterns, we quantify how PIN rearranges in the model with tissue-wide stress patterning ver sus the uncoupled tissue approximation.
We compute the average PIN efflux direction, i.e., average PIN polarity for each combination of the parameters r (auxin-induced cell wall loosening) and η (coupling of PIN to stress) under the approximated (Fig. 7 top left) and tissue-wide ( Fig. 7 top right) stress coupling regimes.
We observe an overall increase in PIN polarity in the tissue-wide stress coupling regime compared with the uncoupled tissue approximation. PIN polarity also becomes more sensitive to r . For very low values of r , tissue stress patterns are slightly detrimental to auxin patterning. These data show that saturation of PIN polarity happens earlier with respect to η for intermediate values of r . For high values of r , we observe a non-monotonic dependence of polarity on η, effectively translating into an optimal value of η.
Visual inspection of the simulations results reveals higher PIN density in proximity of auxin spots and an increase in magnitude of these auxin peaks upon tissue-wide stress patterning ( Fig. 7a-d). Moreover, we observe a severe alteration of pattern size and wavelength between both models (Fig. 7 bottom left).
These results show that tissue-wide stress patterning reinforces PIN polarity and that auxin spots are sharper. Next we will quantify how much sharper these auxin spots become.

Tissue-wide coupling induces efficient emergence of auxin spots
Auxin levels in the shoot apical meristem have been shown to affect cell fate reliably [20], even if the flexibility of the auxin signalling mechanism allows for many potential outcomes [62]. We explore auxin spot concentration achieved by both models in order to gauge the impact of tissue-wide stress patterns on the distinguishability of primordium cells.
For this, we first characterize quantitatively the auxin spot average concentration measured for each simulation of the uncoupled tissue approximation (Fig. 8 left) and tissue-wide for the tissue-wide stress patterning. PIN polarity magnitude increases when considering the mechanics of the whole tissue, with a particularly strong dependence on the wall loosening affect r of auxin. The labels represent the parameters plotted for (a, b, c, d) comparison between example results of auxin concentration and PIN density of simulations using the uncoupled tissue approximation (a, c) and the tissue-wide stress patterning (b, d), for the same value of η = 5.5, and r = 0.65 (bottom left) or r = 0.90 (bottom right). In both instances, we observe that PIN polarity and auxin concentration are higher upon tissue-wide stress patterning (b, d) Fig. 8 Characterization of auxin spot concentration reveals more focused auxin spots due to tissue-wide integration via mechanical coupling. Average auxin concentration for cells above basal auxin concentration (a α > 1), for the uncoupled tissue approximation (left), and upon tissue-wide stress patterning (right), as a function of stress-PIN coupling, η, and wall loosening effect, r . Spot auxin concentration increases with both η and r in (left); however, in (right), it increases predominantly with r . For medium to high values of r , auxin concentration jumps to several times immediately after emergence stress coupling (Fig. 8 right) regimes. We use, as a proxy, the average of cells with auxin concentration a α > 1 to identify auxin spots.
We observe that the dependence on the parameter r recognized for PIN polarity translates into auxin spot concentration. For medium to high values of r , auxin concentration is several Fig. 9 Map of auxin distribution and PIN density aligns with stress direction. Green lines represent principal direction of stress, measured as γ = , where λ ± are the largest and lowest eigenvalues of the stress tensor. We observe that stress directions in part congruent with auxin spots. r = 0.90 and η = 5.5 times higher when accounting for tissue-wide behaviour than when considering the uncoupled tissue approximation.
Additionally, at the onset of pattern formation for medium to high values of r , we observe a considerable jump in average auxin spot concentration for a small change in η. This increase in sensitivity to a change in η of the system, under the aforementioned conditions, implies a boost in mechanosensing capabilities when considering tissue-wide stress patterning.
Our results point to stress patterns being responsible for the enhancement of auxin spot concentration and flows. In order to make sure we understand why, we decided to observe and quantify stress patterns and their connection to auxin distribution.

Part of wall stress within spots is borne by walls at the interface
In order to analyse tissue-wide stress patterns, we choose an example that has simple auxin patterns that allow for a straightforward interpretation. Under this condition, we choose the parameters r = 0.90 and η = 5.5 already presented in Fig. 7d, for which we plot on it a measure of anisotropy along the largest principal stress direction (Fig. 9). Here it becomes apparent that stress patterns are related, even if not absolutely, to auxin spot patterns.
We can section the results according to high or low auxin concentration (Fig. 10 bottom  left), and high or low auxin gradient (Fig. 10 bottom right). Here, high auxin cells are a proxy for auxin spot cells, and high auxin gradient cells are a proxy for cells neighbouring auxin spots. Taking into account that in the uncoupled tissue approximation P α =σ and D α = 0 by construction, we can get a better picture of tissue-wide stress patterns.
We observe from data ( Fig. 10 bottom left) that P α in cells of auxin spots is lower than in the uncoupled tissue approximation and accompanied by a slight shift in the opposite direction of the P α of the remaining cells. Additionally, we register a noticeable shift towards negative D α for high auxin gradient cells (Fig. 10 bottom right). , since P α =σ signifies the stress that would be expected in the uncoupled tissue approximation, high auxin concentration cell expansion is constrained by the remaining cells which are, in turn, under a larger amount of stress. On the other hand (bottom right) we observe that the auxin spot neighbours have, on average, negative values of D α , indicating that the largest principal stress direction is perpendicular to auxin gradients, i.e., circumferentially aligned around auxin spots, as suggested by Fig. 9. r = 0.90 and η = 5.5 Taken together, these data suggests that cell walls at the interface of a spot are under a larger amount of stress whereas the cells within auxin spots have decreased stress. This leads to reinforced polar auxin transport towards the spot and hence higher auxin concentration. The lower isotropic stress component inside the auxin spot suggests that the diffusive term inside auxin spots increases in importance relative to the active transport term.

Tissue-wide stress coupling mitigates disruption by noise
Up until now, our simulations were performed on hexagonal tissues in the absence of noise. This also raises the question of how tissue-wide stress patterns impact pattern emergence robustness against noise.
In plant tissue as any biological entity, noise prevails. As such cells within a tissue differ in their mechanical parameters. In order to inspect how parameter noise disrupts pattern emergence, we choose to sample reference stiffness, E 0 , from a normal distribution for each cell. As outlined in the supplementary material, we expect this parameter to be the most disruptive to the active term and it is reasonable to assume it changes from cell to cell. We then simulate the resulting tissue with the uncoupled tissue approximation and tissue-wide stress coupling.
We simulate tissues with r = 0.65 and η = 5.5 for both models by promoting E 0 to a random variable sampled from Gaussian distribution with meanĒ 0 = 300 MPa and standard deviation of αĒ 0 , α ∈ {0.03, 0.06, 0.09, 0.12, 0.15}, where α is the noise strength. For each value of α, five simulations were performed per model. We fit the resulting auxin distributions for a noise strength of 9%. Even though patterns are heavily disrupted, we can still discern more clearly high auxin concentration spikes upon tissue-wide coupling. (bottom right) Entropy difference between the resulting distributions for a given noise strength and in the absence of noise. In the presence of tissue-wide stress patterns disruption of tissue patterning is consistently lower than in the uncoupled tissue approximation to a probability density function (Fig. 11 top left and top right). We observe that noise in reference stiffness impacts the patterning behaviour in a severe manner ( Fig. 11 bottom left). Yet, with tissue-wide stress coupling spots of noticeable auxin accumulation are preserved. In order to quantify the disruption, we compute the entropy (Eq. 20) of a fitted auxin probability density function by means of a kernel density estimation on the resulting auxin distributions. The kernel used for all fits was the Epanechnikov kernel with a bandwidth of about 0.202. This number arises in the rule-of-thumb estimate for the Gaussian kernel for the sample size and dimension of this system. The partition size used for the numerical approximation of the entropy is the same for all instances. Afterwards, we measure the entropy difference between the expected auxin distribution for each value of α, of each model and for all simulations (Fig. 11 bottom right). The reference entropy is taken to be the average of the uncoupled tissue approximation at α = 0. We can infer from these results that tissue-wide coupling helps to rescue auxin accumulation spots despite its heavy disruption in comparison to the uncoupled tissue approximation.

High turgor preferred but not required for sustaining auxin maxima
It is of interest to the experimental community at this point in time how auxin spots and turgor pressure correlates [63,64]. To explore how the the tissue-wide mechanical model Regardless of initial auxin concentration, for high turgor, we observe an auxin maximum predictably emerges correlated with a turgor maximum (Fig. 12 right column). Nevertheless, if an initial auxin concentration exists, we also predict that the spot remains regardless of whether this position is a turgor minimum or not (Fig. 12 bottom row). Low turgor regions can still exhibit patterns adding to the complexity of this simple measure ( Fig. 12 middle row, left).
From these data, we can conclude is that developmental history is as important as turgor pressure for predicting auxin maxima positioning. We can predict high turgor leads to auxin accumulation, yet low turgor gives us little insight on auxin distribution. We can also predict that a high auxin concentration region.

Discussion
Here, we used a model composed of a vertex model for plant tissue mechanics, and a compartment model for auxin transport to uncover the role of tissue-wide mechanical coupling on auxin redistribution. We first verified that the tissue-wide mechanical model successfully captures the behaviour of plant tissue upon ablation experiments and the conditions for emergence of auxin patterns. We then compared the behaviour of our model featuring tissue-wide mechanical coupling to an approximation which regards cells as mechanically isolated. We observe the emergence of focused auxin spots with high auxin concentration when tissuewide mechanical coupling is implemented. Notably, depending on the parameters of the tissue-wide stress model, auxin spot concentration is more sensitive to stress than what could be predicted from the approximation. We observe that tissue-wide mechanical effects unaccounted for by the approximation have a positive impact on PIN polarity. Furthermore, we show that stress patterning of the tissue mitigates the disruption caused by noise, increasing robustness of the system. Finally, we observe auxin concentration correlating with high local turgor pressure. This behaviour coexists with the possibility of having auxin maxima anti-correlating with turgor.
The auxin-induced cell wall loosening effect (r parameter in this work) is an important determinant of the feedback of auxin on tissue mechanics. The range of values of r for which substantial pattern focusing occurs is around r ∼ 0.60 and above in our model. This translates into a variation of stiffness from a minimum value E min up to E max = 4E min (see supplementary material). Although high, this range is within biological expectation and supported by AFM measurements on auxin treated tissues [7] and comparable to previous simulations of this mechanism [10] where E max /E min = 5 which translates into r = 2/3.
Comparison of the tissue-wide stress patterning case to the uncoupled tissue approximation reveals that auxin spot concentration has a very steep transition in the former case (Fig. 8). This results in a several-fold increase in auxin concentration at values of η close to the threshold for pattern formation. What was once a relatively subtle graded response of auxin spot concentration on stress behaves now as an on-off switch by virtue of tissue mechanical relaxation. Since the mechanical perturbations being highlighted through the comparison are purely passive, this improvement in sensing comes at no additional cost for the plant and therefore has the potential to increase efficiency.
In the present work, we explored the parameter space (η, r ) exclusively. We observe consistently that pattern wavelengths shorten from the uncoupled tissue approximation and the tissue-wide coupling model. It would be interesting to systematically probe the diversity of patterns and how they change upon tissue-wide mechanical coupling. For our simulations, we used the parameters n, m, K from [10], parameters on which we have little empirical information. Yet, the sensitivity analysis from [29] suggests that n and K especially should affect patterning the most. We speculate the parameter m, specific to wall loosening, to be of similar importance. We expect that a study focusing on these three parameters would yield more interesting patterns shapes.
This work focused exclusively in the hypothesis that PIN is mechanically regulated. However, competing chemical feedback mechanisms have been proposed. Recently, mechanics and ARF-mediated PIN expression have been modelled together by [40] and show promising pattern formation capabilities. Other factors we have not taken into account is the family of auxin importers of the AUX family, which have been shown to be present in the epidermal layer of the shoot apical meristem [16]. Auxin binding proteins have also been hypothesized to promote auxin flow polarization [65]. Another observed interaction is cytokinin action controlling PIN polarity during lateral root formation [66].
The PIN regulation used in the auxin transport compartment model was specifically stressbased. In the supplementary text, we show results using strain-based PIN binding instead. We observe the same overall auxin spot focusing behaviour. It is still unclear whether the PIN density change due to mechanics is a result of strain or stress [11]. In fact, this question has been tackled recently by [54] concluding that in most simulated experiments both strain and stress-based models behave similarly. A notable exception is the experimentally observed correlation of PIN polarity and auxin concentration [37,67]. On one hand, this observation is not captured by the stress-based PIN binding model. On the other hand, available experimental data and simulations suggest stress sensing being easier to explain [54]. Furthermore, the polarity difference could be rescued by the observation that ARF-mediated PIN expression is higher at the tip of the primordium [68].
The specific distribution of emergent stress patterns is remarkable in the sense that it coincides with the shape-induced stress patterns, as indicated by microtubule orientation, around the tip of the primordium as it emerges from the meristem [41]. Therefore, tissue-wide stress patterning sets the stage for primordium outgrowth by focusing efficiently auxin, forming local circumferential stress that in turn may re-orient microtubules and prefigure the shape of the primordium. This process could, in turn, be capable of reinforcing auxin transport to the tip of the newly forming organ. Yet, quantifying this requires further modelling. Therefore, it would be interesting to include auxin transport in already existing models for primordium outgrowth [8,9].
Even though the analysed numerical simulations were limited to noise in the parameter E 0 , it showcases the power of the aforementioned auxin peak focusing that happens upon tissue-wide mechanical coupling. In this instance, we show here the power of tissue-wide stress patterns to mitigate the information loss due to noise by inspecting entropy of auxin concentration distribution. This result, especially when paired with the increase in sensitivity mentioned above, is indeed remarkable. This is due to the fact that in a wide range of optimized systems, biological or otherwise, robustness and efficiency are thought to be in opposition to each other, as illustrated, for example, by [69]. This is because robustness is usually brought upon by additional systems which would be considered clutter by a system geared towards efficiency. This opens a novel line of argumentation in the discourse of the evolution of mechanical signalling in multicellular organisms.
Lastly, we probed the behaviour of the used tissue-wide mechanical model when faced with local turgor variations. Our results indicate that once established auxin spots can endure low turgor scenarios, even if they would prefer high turgor regions all else being equal. Maintaining a turgor pressure difference for so long, however, might not be feasible for the plant. To answer how this setting could be achieved would require modelling water transport between plant cells along the lines of [63]. Nevertheless, our model can explain, at least in part, why these two quantities do not correlate in a straightforward manner.

Conclusion
Even though the mechanisms by which PIN preferentially associate with stressed cell walls is unclear, here we show that there are substantial advantages by intertwining tissue-wide mechanics and auxin patterning. Even if auxin patterning is possible by chemical processes and local mechanical coupling, tissue-wide mechanics may provide a way for patterning to still occur at a lower energy cost for the tissue. Moreover, this process can also provide robustness to the patterning, factoring in tissue-wide stress pattern, a sort of proprioceptive mechanism.