18 figures, 3 tables and 23 additional files

Figures

PIN1::PIN1:GFP expression in young WT and kan1kan2 leaf primordia.

(A) WT primordium of leaf 1, showing abaxial epidermis (a total of 10 leaves were imaged over 2 separate experiments). (B) As for A, but for a kan1kan2 primordium (a total of 15 leaves were imaged …

https://doi.org/10.7554/eLife.18165.003
PIN1::PIN1:GFP polarity patterns in WT and kan1kan2 leaf development.

Confocal images of the PIN1::PIN1:GFP marker in the abaxial epidermis of the same WT leaf primordium imaged over a period of 2 days. Approximate leaf widths (measured from projections of the z …

https://doi.org/10.7554/eLife.18165.004
DR5::GFP expression in a kan1kan2 leaf during outgrowth development.

Confocal images of DR5::GFP in the same kan1kan2 leaf imaged over a period of 3 days. Times relative to the first observation of an outgrowth, and leaf widths, are given above images. B ii, C ii and …

https://doi.org/10.7554/eLife.18165.005
Figure 4 with 1 supplement
Leaf phenotype of the kan1kan2cuc2 mutant.

(A) Whole kan1kan2 plant. (B) OPT images of a kan1kan2 leaf showing abaxial leaf surface. (C) Whole kan1kan2cuc2 plant. (D) OPT images of a kan1kan2cuc2 leaf. Scale bars in A and C = 1 cm, scale …

https://doi.org/10.7554/eLife.18165.006
Figure 4—figure supplement 1
Abaxial ridges and serrations produced in kan1kan2cuc2 mutants.

(A) Example of a kan1kan2cuc2 leaf with an abaxial ridge-like thickening (indicated by blue arrow). (B) Example of a kan1kan2cuc2 leaf with serrations. Scale bars = 1 mm.

https://doi.org/10.7554/eLife.18165.007
PIN1::PIN1:GFP in leaf one of the kan1kan2cuc2 mutant.

(A) Confocal image of PIN1::PIN1:GFP in the abaxial epidermis of the first leaf primordium of a kan1kan2cuc2 mutant. (B) Time-lapse confocal images of the abaxial side of the first leaf of a kan1kan2…

https://doi.org/10.7554/eLife.18165.008
Expression of CUC2::RFPer and PIN1::PIN1:GFP during kan1kan2 outgrowth development.

(A) Confocal images of a CUC2::RFPer reporter in the abaxial epidermis of a young kan1kan2 first leaf primordium. (B) Combined CUC2::RFPer (red) and PIN1::PIN1:GFP confocal channels for the same …

https://doi.org/10.7554/eLife.18165.009
Figure 7 with 1 supplement
Comparison of basic behaviours of indirect coupling, flux-based and up-the-gradient models.

(AC) Investigating the ability of cells to polarise in the absence of pre-established external asymmetries or polarisable neighbours for up-the-gradient (A), flux-based (B) and indirect coupling (C)…

https://doi.org/10.7554/eLife.18165.010
Figure 7—figure supplement 1
Comparison of models’ ability to generate polarity for a cell surrounded by a fixed environment.

In each case the central cell is initialized with noise in the concentrations of PIN at each cell edge and all seven cells have the same auxin concentration. The auxin concentration in each outside …

https://doi.org/10.7554/eLife.18165.011
Formation of a proximo-distal polarity field and centres of convergence in the convergent alignment model.

(A) Formation of a proximo-distal polarity field due to the presence of an elevated initial auxin concentration at the leaf tip (orange cells) and an elevated rate of auxin removal from the leaf …

https://doi.org/10.7554/eLife.18165.012
Formation of a proximo-distal polarity field in the flux-based and indirect coupling models.

(A) Formation of a proximo-distal polarity field in the flux-based model due to the presence of elevated auxin biosynthesis at the leaf base (orange cells) and an elevated rate of auxin removal from …

https://doi.org/10.7554/eLife.18165.013
Figure 10 with 1 supplement
Formation of centres of polarity convergence in the flux-based and indirect coupling models.

(A) A proximo-distal polarity field is initially established with the flux based model due to the presence of elevated auxin biosynthesis in the proximal half of the leaf (orange cells) and elevated …

https://doi.org/10.7554/eLife.18165.014
Figure 10—figure supplement 1
Incorporation of a D6 protein kinase-like activity into the indirect coupling model.

D6 kinase activity was added into simulations of proximodistal polarity establishment (A), (compare with Figure 9D) and convergence formation (B) (compare with Figure 10D). We assume that D6 protein …

https://doi.org/10.7554/eLife.18165.015
Figure 11 with 1 supplement
Expression of LAX1::GUS and AUX1::AUX1:YFP in WT and kan1kan2 leaves.

(A) Expression pattern of LAX1::GUS in WT leaves one and two. LAX1::GUS was expressed at the tips of developing primordia (arrows in (i), black dashed lines indicate leaf outlines, arrow heads …

https://doi.org/10.7554/eLife.18165.016
Figure 11—source data 1

Counts of outgrowths used to generate Figure 11G.

https://doi.org/10.7554/eLife.18165.017
Figure 11—figure supplement 1
Expression patterns of LAX2::GUS and LAX3::GUS.

(A) Expression of LAX2::GUS in a kan1kan2+/- leaf, showing expression in leaf vascular tissue. Scale bar = 100 μm. (B) Expression of LAX3::GUS in a WT seedling, showing expression is absent from …

https://doi.org/10.7554/eLife.18165.018
LAX1 and AUX1 expression in kan1kan2cuc2 mutants.

(A) Expression of LAX1::GUS in the abaxial surface of leaf 1 of a kan1kan2cuc2 mutant. Data representative of images from ten seedlings in two separate experiments. (B) AUX1::AUX1:YFP expression in …

https://doi.org/10.7554/eLife.18165.019
PIN1 immuno-localisation in transverse cross-sections of kan1kan2 leaves.

(A) PIN1 immuno-localisation in a kan1kan2 leaf before outgrowth emergence. i) and ii) are consecutive sections through the tissue (each section is 8 µm thick). A centre of PIN1 polarity convergence …

https://doi.org/10.7554/eLife.18165.020
YUCCA1::GUS expression in WT leaves.

(A) Expression of YUC1::GUS in a WT leaf 1 primordium. i) abaxial epidermis (10 plants were imaged across two separate experiments) ii) transverse cross-sectional view of leaf 1 at a similar …

https://doi.org/10.7554/eLife.18165.021
Expression of YUCCA1::GUS in kan1kan2 leaves.

(A) Expression of YUC1::GUS in kan1kan2 leaf one primordia. i) abaxial epidermis (a total of 15 leaves imaged, in three separate experiments), ii) transverse cross-sectional view of leaf 1 at a …

https://doi.org/10.7554/eLife.18165.022
Effect of a band of locally elevated auxin biosynthesis on flux-based, indirect coupling and up-the-gradient models.

(A) Flux-based model. A proximodistal polarity field is initially established due to elevated auxin production at the leaf base (orange cells) and elevated auxin import and removal at the leaf tip …

https://doi.org/10.7554/eLife.18165.023
Time-lapse imaging of YUC4::GFP in kan1kan2 leaves.

(A) Confocal images of YUC4::GFP in the same kan1kan2 leaf imaged over a period of 3 days. Times relative to outgrowth emergence, and leaf widths, are given above images. Yellow circles indicate …

https://doi.org/10.7554/eLife.18165.024
Classification of plant and animal polarity models.

Models of plant and animal epidermal polarity may be classified according to whether they require the presence of pre-established asymmetries or polarisable neighbours (referred to as a polarising …

https://doi.org/10.7554/eLife.18165.025

Tables

Table 1

Parameter values used in simulations of the indirect coupling model.

https://doi.org/10.7554/eLife.18165.026
Symboldescriptionunitvalue
Δtnumerical time steps (seconds)0.05
Rcarea of cytoplasmic compartmentµm2260
Rwarea of cell wall compartmentµm22.88 (for long cell edges); 3.0 (for short cell edges)
lnlength of membrane compartmentsµm2.88; 3.0
lwlength of wall compartmentsµm2.88; 3.0
cPINdefault initial concentrations of PIN1 in cytoplasmic compartmentsAu/ µm20.003
ρPINPIN1 default membrane binding rateµm/s0.03
τA*-dependent promotion of PIN1 bindingµm2/Au.s2
µPINPIN1 default membrane unbinding rate/s0.004
DPINPIN1 diffusion in cell membraneµm2/s0.1
ΨPINPIN-dependent active auxin efflux rateµm2/Au.s40
ɛlimit for noise addition during initialisation of A* and B* concentrationsdimensionless0.0166
γAuxAuxin-dependent promotion of A* to A conversionµm2/Au.s0.75*
ρAuxproduction rate of AuxinAu/µm2.s1.0 × 10−4*
µAuxdegradation rate of Auxin/s0.01*
vininflux auxin permeabilityµm/s0.75*
  1. All parameters relate to either equations specified here or in Abley et al., 2013.

  2. * In the simulations used to generate Figure 7C and Figure 7—figure supplement 1C, γAux = 0.9 µm2/Au.s, ρAux = 1.3 × 10−4 Au/µm2.s and µAux = 0.02/s. Also, in some simulations, at tissue boundaries or centres of polarity convergence, auxin production, degradation and influx auxin permeability rates may vary from the background rates. Modulation of the influx auxin permeability rate is used to simulate elevated rates of auxin import. In Figure 9C, the auxin production rate in the bottom-most row of cells = 2 × 10−4 Au/µm2.s, and the auxin degradation rate in the top most row of cells = 0.07/s. In Figure 9D, and at the beginning of the simulation used to generate Figure 10C and D the auxin production rate in the bottom-most row of cells = 2 × 10−4Au/µm2.s. Throughout both simulations, in the top-most row, the auxin degradation rate = 0.015/s and the auxin influx permeability = 4.5 µm/s. In the simulation used to generate Figure 10C and D, changes in auxin production occur during the simulation and are described below. In this simulation, in the cell that forms the centre of convergence, the inwards permeation of auxin = 37.5 µm/s and the auxin degradation rate = 0.15/s. In Figure 16C and D, the proximo-distal polarity field is established as described for Figure 9D. In the three cells with elevated auxin production, the auxin production rate is 1.4 × 10−3Au/µm2.s. In the cell with elevated auxin import and removal (Figure 16D), the auxin degradation rate = 0.1 /s and the auxin influx permeability = 37.5 µm/s. In the simulations used to generate Figure 10—figure supplement 1A and B, parameter values are as for Figures 9D and 10D, respectively except ΨPIN = 80 µm2/Au.s and γAux = 0.85 µm2/Au.s. Here, c, the concentration of A* at which half of the PINs are activated = 0.22 Au/µm and the hill coefficient, h = 2.

Table 2

Parameter values used in simulations of the flux-based model.

https://doi.org/10.7554/eLife.18165.027
Symboldescriptionunitvalue
Δtnumerical time steps (seconds)0.01
Rarea of cytoplasmic compartmentµm2260
llength of cell edge compartmentµm15 for short cells edges; 8.6 for the each of the two compartments of a long edge *
CAdefault initial concentration of auxin in cytoplasmic compartmentsAu/ µm20.01
CPINdefault initial concentration of PIN1 at cell edge compartmentsAu /µm0.01
ɛPINlimit for noise addition during initialisation of PIN1 concentrations at cell edgesdimensionless0.025**
ɛAuxlimit for noise addition during initialisation of Auxin concentrationsdimensionless0.025
TPIN-dependent active auxin efflux rateµm2/Au.s1
γUnbinding rate of PIN1 from the cell edge/s0.1
ρproduction rate of AuxinAu/µm2.s0.0001§
µdegradation rate of Auxin/s0.02
αFlux-dependent promotion of PIN1 allocation to a cell edgedimensionless1#
DPassive permeation rate of auxinµm/s5
PmaxMaximum concentration of PIN1 at a cell edgeAu /µm0.01
IAuxin import rateµm / s0
  1. * For cells with hexagonal geometries, each cell edge has a length of 10 µm.

  2.  Value given is that used in simulations used to generate Figure 7.B and Figure 7—figure supplement 1B. In all other simulations, CPIN = 0.

  3. ‡ Values apply to simulations used to generate Figure 7B (ɛPIN) and Figure 7—figure supplement 1B, and Figure 7E (ɛAux). All other simulations are initialised without noise addition.

  4. § The value given for the auxin production rate applies to simulations used to generate Figures 9A,B, 10A, B and 16A,B. In the simulation used to generate Figure 7B and Figure 7—figure supplement 1B, ρ = 0.0025 Au/µm2.s and in the simulation used to generate Figure 7E, ρ = 0.002 Au/µm2.s.

  5. # The value given for, α, the flux-dependent promotion of PIN allocation to cell edges, applies to simulations used to generate Figures 9A,B, 10A,B and 16A,B. In the simulation used to generate Figure 7B, α = 4 × 10−3 and in the simulation used to generate Figure 7E, α = 3.2 × 10−3.

  6. In the simulation used to generate Figures 7B, Figure 7—figure supplement 1B and 7E, Pmax = 0.04Au /µm.

Table 3

Parameter values used in up-the-gradient simulations.

https://doi.org/10.7554/eLife.18165.028
Symboldescriptionunitvalue
Δtnumerical time steps (seconds)0.01
Rarea of cytoplasmic compartmentµm2260
llength of cell edge compartmentµm15 for short cells edges; 8.6 for the each of the two compartments of a long edge *
CAdefault initial concentration of auxin in cytoplasmic compartmentsAu/ µm20.01 (Figure 7A,D, Figure 7—figure supplment 1A.)
0.005 (other Figs.)
CPINdefault initial concentration of PIN1 at cell edge compartmentsAu /µm0.1
ɛPINlimit for noise addition during initialisation of PIN1 concentrations at cell edgesdimensionless0.025
ɛAuxlimit for noise addition during initialisation of Auxin concentrationsdimensionless0.025
TPIN-dependent active auxin efflux rateµm2/Au.s80
ρproduction rate of AuxinAu /µm2.s0.0003
µdegradation rate of Auxin/s0.005
PINiTotal amount of PIN1 in a cell available for binding to edge compartmentsAu /µm0.1
DPassive permeation rate of auxinµm/s10
bExponentiation base for PIN1 allocation to the membranedimensionless6
IAuxin import rateµm / s0
  1. * For cells with hexagonal geometries, each cell edge has a length of 10 µm.

  2.  Values given for the initial concentration of PIN1 at cell edge compartments, and the limit for noise addition to this concentration, are those used in simulations used to generate Figure 7A and Figure 7—figure supplement 1A. In all other simulations, CPIN and ɛPIN = 0. The value given for ɛAux applies only to Figure 7D, in all other simulations, ɛAux = 0.

  3. Value given for the production rate of auxin applies to the simulations used to generate Figures 7D, 8A, B, 11H, 14D and 16E,F. In the simulation used to generate Figure 7A and Figure 7—figure supplement 1A, ρ = 0.0025 Au /µm2.s.

Additional files

Source code 1

Code used to generate Figure 7A.

https://doi.org/10.7554/eLife.18165.029
Source code 2

Code used to generate Figure 7B.

https://doi.org/10.7554/eLife.18165.030
Source code 3

Code used to generate Figure 7C.

https://doi.org/10.7554/eLife.18165.031
Source code 4

Code used to generate Figure 7D.

https://doi.org/10.7554/eLife.18165.032
Source code 5

Code used to generate Figure 7E.

https://doi.org/10.7554/eLife.18165.033
Source code 6

Code used to generate Figure 7F.

https://doi.org/10.7554/eLife.18165.034
Source code 7

Code used to generate Figure 8A.

https://doi.org/10.7554/eLife.18165.035
Source code 8

Code used to generate Figure 8B.

https://doi.org/10.7554/eLife.18165.036
Source code 9

Code used to generate Figure 9A.

https://doi.org/10.7554/eLife.18165.037
Source code 10

Code used to generate Figure 9B.

https://doi.org/10.7554/eLife.18165.038
Source code 11

Code used to generate Figure 9C.

https://doi.org/10.7554/eLife.18165.039
Source code 12

Code used to generate Figure 9D.

https://doi.org/10.7554/eLife.18165.040
Source code 13

Code used to generate Figure 10A and B.

https://doi.org/10.7554/eLife.18165.041
Source code 14

Code used to generate Figure 10C and D.

https://doi.org/10.7554/eLife.18165.042
Source code 15

Code used to generate Figure 11H.

https://doi.org/10.7554/eLife.18165.043
Source code 16

Code used to generate Figure 14D.

https://doi.org/10.7554/eLife.18165.044
Source code 17

Code used to generate Figure 16A.

https://doi.org/10.7554/eLife.18165.045
Source code 18

Code used to generate Figure 16B.

https://doi.org/10.7554/eLife.18165.046
Source code 19

Code used to generate Figure 16C.

https://doi.org/10.7554/eLife.18165.047
Source code 20

Code used to generate Figure 16D.

https://doi.org/10.7554/eLife.18165.048
Source code 21

Code used to generate Figure 16E.

https://doi.org/10.7554/eLife.18165.049
Source code 22

Code used to generate Figure 16F.

https://doi.org/10.7554/eLife.18165.050
Supplementary file 1

Supplementary model information.

Instructions on how to run models and explanation of the code for each model.

https://doi.org/10.7554/eLife.18165.051

Download links