# A Mechanical Model to Interpret Cell-Scale Indentation Experiments on Plant Tissues in Terms of Cell Wall Elasticity and Turgor Pressure

Richard Malgat, François Faure, Arezki Boudaoud

► To cite this version:

Richard Malgat, François Faure, Arezki Boudaoud. A Mechanical Model to Interpret Cell-Scale Indentation Experiments on Plant Tissues in Terms of Cell Wall Elasticity and Turgor Pressure. Frontiers in Plant Science, 2016, 7, pp.1351. 10.3389/fpls.2016.01351. hal-01367360

HAL Id: hal-01367360

https://hal.science/hal-01367360

Submitted on 16 Sep 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Mechanical Model to Interpret Cell-Scale Indentation Experiments on Plant Tissues in Terms of Cell Wall Elasticity and Turgor Pressure

Richard Malgat<fnref n="1" />,²,³, François Faure<fnref n="1" />,² and Arezki Boudaoud³*

Morphogenesis in plants is directly linked to the mechanical elements of growing tissues, namely cell wall and inner cell pressure. Studies of these structural elements are now often performed using indentation methods such as atomic force microscopy. In these methods, a probe applies a force to the tissue surface at a subcellular scale and its displacement is monitored, yielding force-displacement curves that reflect tissue mechanics. However, the interpretation of these curves is challenging as they may depend not only on the cell probed, but also on neighboring cells, or even on the whole tissue. Here, we build a realistic three-dimensional model of the indentation of a flower bud using SOFA (Simulation Open Framework Architecture), in order to provide a framework for the analysis of force-displacement curves obtained experimentally. We find that the shape of indentation curves mostly depends on the ratio between cell pressure and wall modulus. Hysteresis in force-displacement curves can be accounted for by a viscoelastic behavior of the cell wall. We consider differences in elastic modulus between cell layers and we show that, according to the location of indentation and to the size of the probe, force-displacement curves are sensitive with different weights to the mechanical components of the two most external cell layers. Our results confirm most of the interpretations of previous experiments and provide a guide to future experimental work.

**Keywords:** mechanical model, cell wall, turgor, physically-based simulation, atomic force microscope, indentation, shoot apical meristem, floral meristem

## 1. INTRODUCTION

Morphogenesis relies on well-defined patterns of growth determined by gene expression (Coen et al., 2004). Mechanistically, it is thought that gene activity influences the mechanical properties of tissues at cellular and sub-cellular level (Traas and Monéger, 2010; Mirabet et al., 2011; Robinson et al., 2013; Ali et al., 2014), thus controlling patterns of growth. Plants are well-suited to investigate the mechanistic basis of morphogenesis: Growth is limited by the polysaccharide-made walls that surround cells and is driven by the osmotically-generated turgor pressure, which puts walls in tension (Schopfer, 2006).

**OPEN ACCESS**

**Edited by:**
Joshua L. Heazlewood,
University of Melbourne, Australia

**Reviewed by:**
Rosemary Dyson,
University of Birmingham, UK
Abraham D. Stroock,
Cornell University, USA

*Correspondence:
Arezki Boudaoud
arezki.boudaoud@ens-lyon.fr

**Specialty section:**
This article was submitted to Plant Biophysics and Modeling, a section of the journal Frontiers in Plant Science

**Received:** 27 January 2016
**Accepted:** 23 August 2016
**Published:** 07 September 2016

**Citation:** Malgat R, Faure F and Boudaoud A (2016) A Mechanical Model to Interpret Cell-Scale Indentation Experiments on Plant Tissues in Terms of Cell Wall Elasticity and Turgor Pressure. Front. Plant Sci. 7:1351. doi: 10.3389/fpls.2016.01351

<footnote n="1">¹ Institut National de Recherche en Informatique et en Automatique, Grenoble, France, ² Laboratoire Jean Kuntzmann, Centre National de la Recherche Scientifique, Grenoble, France, ³ Reproduction et Développement des Plantes, Université de Lyon, Ecole Normale Supérieure de Lyon, Université Claude Bernard Lyon 1, Institut National de la Recherche Agronomique, Centre National de la Recherche Scientifique, Lyon, France</footnote>
Several lines of evidence indicate links between gene expression, cell wall properties and growth, notably at the shoot apex. The activity of cell wall remodeling proteins, expansins (Fleming et al., 1997; Cho and Cosgrove, 2000; Pien et al., 2001) and pectinmethyl esterases (Peaucelle et al., 2008), enables outgrowth or enhances growth. The first steps of flower organogenesis are associated with a reduction in cell wall stiffness (Peaucelle et al., 2011; Braybrook and Peaucelle, 2013). More generally, there is an inverse correspondence between wall stiffness and growth rate (Milani et al., 2013). Finally, cell stiffness correlates with the expression of CLAVATA3, the glycopeptide associated with the central zone (Milani et al., 2014).

Here, we aim at providing a theoretical framework to interpret measurements of plant mechanics at cellular resolution. Frameworks at organ scale have a long history (Niklas, 1992) and enable the deduction of average properties of cell layers from experiments where the whole organ is stretched or bent (Niklas, 1992), based for instance on models that link cell wall elastic modulus (the higher the modulus, the stiffer the wall), turgor pressure, and average elastic modulus of the tissue (Nilsson et al., 1958).

Experimental approaches to measure mechanical properties have been recently scaled down thanks to nano-indentation systems, such as atomic force microscopes, whereby a nanometric to micrometric probe is used to apply a force in the nN-µN range on the sample of interest while its displacement is monitored, yielding force-displacement curves at well-defined locations. Measurements were performed on the shoot apical meristem of Arabidopsis (Milani et al., 2011; Peaucelle et al., 2011; Braybrook and Peaucelle, 2013; Milani et al., 2014), on roots (Fernandes et al., 2012), on cotyledons and leaves (Hayot et al., 2012; Forouzesh et al., 2013; Sampathkumar et al., 2014), on onion scales (Lintilhac et al., 2000; Routier-Kierzkowska et al., 2012; Beauzamy et al., 2015a), on pollen tubes (Geitmann and Parre, 2004; Vogler et al., 2013), or on culture cells (Radotić et al., 2012; Weber et al., 2015). These approaches (reviewed in Geitmann, 2006; Milani et al., 2013; Routier-Kierzkowska and Smith, 2013; Beauzamy et al., 2014; Vogler et al., 2015) have contributed to a renewal of plant biomechanics (Moulia, 2013).

The interpretation of nano-indentation experiments raises a number of questions. Indeed, when a piece of wall is stretched, the elastic modulus that quantifies wall stiffness can be obtained directly from the slope of force vs. displacement curve because the surface on which the force is applied remains constant (Boudaoud, 2010), whereas in indentation experiments on live tissues, the surface of contact between the probe and the sample increases with depth and, in addition, turgor pressure may contribute to the mechanical response. What information is revealed by force-displacement curves? Are measurements sensitive to the cell wall elasticity, cell wall viscosity, and/or to turgor pressure? Are measurements sensitive to the properties of the cell indented or of a larger group of cells? Do probes of different sizes reveal properties at different scales (Peaucelle et al., 2011)? These are typical questions that we tackle in the present study. To do so, we need to build realistic mechanical models of plant tissues.

Such models have been worked out at the subcellular scale. When indentation depth is small with respect to wall thickness and the wall is assumed to be locally homogeneous, standard models from contact mechanics (Johnson, 1987) enable the deduction of the transverse elastic modulus (Milani et al., 2011), which quantifies the wall stiffness in the direction of its thickness. The assumption of homogeneity can be relieved by assuming that the elastic modulus varies smoothly with the distance normal to the surface (Lee et al., 2009) or that the wall is made of two types of homogeneous materials (Roduit et al., 2009). The assumption of linear elastic behavior of the material can be relieved in models that include nonlinear elasticity (Valero et al., 2016) or viscoelastoplastic behavior (Tvergaard and Needleman, 2011). Cell-scale models can account for the cell wall and for turgor (Routier-Kierzkowska et al., 2012), leading to the deduction of wall elastic modulus (Hayot et al., 2012) or turgor pressure (Forouzesh et al., 2013) from indentation experiments. Physical models can also be combined to extract many physical parameters from a single force-displacement curve (Beauzamy et al., 2015a; Bonilla et al., 2015). Incidentally, inferring turgor provides a useful alternative to the more standard pressure probe (Tomos and Leigh, 1999), in which cells are impaled for measurements. In tissue-scale models (Hamant et al., 2008; Bassel et al., 2014; Sampathkumar et al., 2014), too much spatial detail—a resolution smaller than wall thickness—would lead to a huge computational time; therefore cell walls are modeled as thin plates or shells and wall thickness is accounted for only through stretching and bending moduli, which quantify the stiffness of the plate when stretched or bent, respectively. In such shell models, indentation has been explored only in geometries corresponding to single cells in the shape of a sphere (Vella et al., 2012b), an ellipsoid (Vella et al., 2012a), or a capped cylinder (Vogler et al., 2013; Weber et al., 2015), which helped the deduction of turgor pressure in various systems (Vogler et al., 2013; Beauzamy et al., 2015a,b; Weber et al., 2015).

Here we aim at generalizing shell models to tissues; accordingly, we do not make predictions at the sub-wall scale, but rather aim at providing a framework to guide indentation experiments at the cell scale. Our approach (Malgat et al., 2014) is based on the Simulation Open Framework Architecture (SOFA, Faure et al., 2012). We simulate the indentation of a floral meristem (Fernandez et al., 2010; Boudon et al., 2015) and vary turgor pressure and the elastic moduli of the different types of walls. We finally discuss our results in relation with recent experiments at cell scale (Peaucelle et al., 2011; Fernandes et al., 2012; Braybrook and Peaucelle, 2013; Milani et al., 2014; Beauzamy et al., 2015b; Peaucelle et al., 2015).

## 2. MATERIALS AND METHODS

In Malgat et al. (2014), we developed a mechanical model for the indentation of a plant tissue and we validated it for a square-shaped cell wall; we also illustrated the model with a template originating from a flower bud, assuming all cell walls to have the same elastic modulus. Here, we briefly recall the main modeling ingredients and explain the new ones.
## 2.1. Flower Bud Template

The template was obtained previously (Fernandez et al., 2010; Boudon et al., 2015). A flower bud was imaged by confocal microscopy and the 3D image was processed so as to extract the center of mass of each cell and the surface of the bud. The bud was partitioned into virtual cells using a Voronoi tessellation closed by the surface of the bud. The facets of cells were meshed with triangles and the mesh was refined near the probe in order to improve precision, see Figure 1A. For certain simulations, we defined categories of walls, based on the organization of the bud into cell layers: the L1 epidermal layer, the L2 subepidermal layer, and other cells (L3).

We also deformed this template by a small random displacement of all vertices and found force-displacement curves to be rather insensitive to such changes in the mesh (Supplementary Figure 1).

## 2.2. Constitutive Elements of the Model

Each cell wall (facet) was assumed to behave as a thin plate of thickness $h$, elastic modulus $E$, and Poisson ratio $\nu$. Accordingly

each wall has a stretching modulus $Eh$ and a bending modulus $Eh^3/12(1 - \nu^2)$. In addition to our previous work (Malgat et al., 2014), we also implemented anisotropy in the stretching modulus, which is modeled by the generalization of Hooke's law to an orthotropic material. Each thin plate is in plane stress, so that the stress tensor $\sigma$ reduces to its 3 components in the $(x, y)$ coordinates of the plate. Taking the notation $\epsilon$ for the strain tensor, Hooke's law can then be written as

$$ \begin{bmatrix} \sigma_{xx} \\ \sigma_{yy} \\ \sigma_{xy} \end{bmatrix} = \frac{1}{1 - \nu_{xy}\nu_{yx}} \begin{bmatrix} E_x & \nu_{yx}E_x & 0 \\ \nu_{xy}E_y & E_y & 0 \\ 0 & 0 & G_{xy}(1 - \nu_{xy}\nu_{yx}) \end{bmatrix} \begin{bmatrix} \epsilon_{xx} \\ \epsilon_{yy} \\ 2\epsilon_{xy} \end{bmatrix} $$

with $E_i$ the elastic modulus along axis $i$, $G_{xy}$ the shear modulus in direction and $\nu_{ij}$ is the Poisson's ratio that corresponds to a contraction in direction $j$ when an extension is applied in direction $i$. Symmetry gives the constraint $\nu_{xy}/E_x = \nu_{yx}/E_y$. In the simulations with mechanical anisotropy, we varied the ratio $\alpha = E_x/E_y$, keeping constant the following quantities: average modulus $E = (E_x + E_y)/2$, first Poisson's ratio $\nu_{xy} = \nu$, and the shear modulus $G_{xy} = E/2(1 + \nu)$. Like other

<figure><img src="image_2.png" /></figure>
mechanical properties, mechanical anisotropy was assumed to be stationary. The direction of the stiffest direction was fixed based on qualitative observations of the cytoskeleton (see below).

The mechanical model was solved using the finite element method as implemented in the Simulation Open Framework Architecture (SOFA, Faure et al., 2012) and detailed in Malgat et al. (2014).

Also in addition to our previous work (Malgat et al., 2014), we considered a linear viscoelastic behavior, restricting ourselves to isotropic materials, whereby $E$ is replaced by $E + \eta \frac{d}{dt}$. The corresponding viscoelastic time scale is $\tau = \eta/E$. This constitutive law describes cell wall behavior at time scales that are small with respect to growth.

## 2.3. Physical Parameters Explored

The range of physical parameters is given in Table 1. Indentations were performed in the middle of a cell (around which the mesh was refined) and at its periphery (above an anticlinal wall). Except for the anisotropic cases, we show hereafter results obtained on a single cell from the top of the flower bud, though these results are typical of such locations, as our preliminary exploration revealed that the force-displacement curves did not vary qualitatively according to the exact location in the cell or to the cell probed.

For each static force-displacement curve, equilibrium is computed for 30–50 successive values of the force (enough values to obtain a smooth curve), with force increments in the 10 nN range, using at each step a full implicit Euler solver coupled with a conjugate gradient solver (Malgat et al., 2014) until the convergence criterion is reached. As a consequence, each force-displacement curve requires a computation time ranging from 6 to 9 h on a standard desktop computer. The results presented here involve about a hundred such force-displacement curves.

In the case of the viscoelastic model, the force was increased according to the equation $F = F_m t/T$, where $T$ is the time needed to reach the maximal value $F_m$, and then decreased according to $F = F_m(2 - t/T)$. As we do not wait for equilibrium, simulations are significantly faster, yielding a force-displacement curve in less than an hour.

# 3. RESULTS

## 3.1. Homogeneous Stiffness of Cell Walls

We simulated the indentation of a flower bud by a spherical probe (Figures 1B–D) and we first assumed that all cell walls had the same elastic modulus (Figure 2A). We initially took a probe diameter of 5 µm (Peaucelle et al., 2011; Braybrook and Peaucelle, 2013). In all cases, we found strain to be smaller than about 5%, consistent with our use of a linearly elastic material (Supplementary Figure 2).

We started with no pressure, as would be the case when tissues are plasmolyzed (Peaucelle et al., 2011; Braybrook and Peaucelle, 2013; Peaucelle et al., 2015). Force-displacement curves differ qualitatively according to whether indentation is performed near the cell middle or near an anticlinal cell wall (wall roughly perpendicular to the surface), as seen in Figure 2B. The vicinity of anticlinal cell walls appear effectively stiffer as about three times higher forces are needed to reach the same depth;

<table><caption>TABLE 1 | Range of physical parameters.</caption><tr><td>Typical cell size <i>R</i></td><td>5 µm</td></tr><tr><td>Modulus (<i>E</i>)</td><td>range 5–500 MPa</td></tr><tr><td>Poisson's ratio (<i>ν</i>)</td><td>0.49</td></tr><tr><td>Thickness (<i>h</i>)</td><td>0.4 µm</td></tr><tr><td>Mechanical anisotropy (<i>α</i> = <i>E</i><sub>x</sub>/<i>E</i><sub>y</sub>)</td><td>1:1 to 5:1</td></tr><tr><td>Pressure (<i>P</i>)</td><td>range 0.1–1 MPa</td></tr><tr><td>Probe diameter</td><td>0.6–5 µm</td></tr><tr><td>Maximum displacement (<i>δ</i>)</td><td>10% of cell size or 0.5 µm</td></tr><tr><td>Viscoelastic timescale (<i>τ</i> = <i>η</i>/<i>E</i>)</td><td>0.1–10 s</td></tr><tr><td>Time for force increase (<i>T</i>)</td><td>1 s</td></tr></table>

the corresponding force-displacement, $F(\delta)$, curves exhibit two regimes, roughly linear at small displacement and superlinear (concave) at larger displacement. We fitted all curves to power-law equations with prefactor $a$ and exponent $b$, $F = a\delta^b$, taking the origin of displacement at $\delta = 3.5\%$ for the second regime at anticlinal walls. The exponents were $b = 2$ at cell middle and $b = 1$ (small depth) or $b = 1.3$ (large depth) at anticlinal walls. The superlinear curves are consistent with experimental observations (Peaucelle et al., 2011; Braybrook and Peaucelle, 2013). The linear behavior has not been observed, but it occurs at displacements smaller than wall thickness, which falls out of the scope of our model. All prefactors are proportional to the elastic modulus (Figure 2C), which yields the only force scale in the absence of turgor pressure.

We then switched to a pressurized state corresponding to experiments on turgid tissues (Milani et al., 2014; Beauzamy et al., 2015b) and took a reasonable value of pressure of $P = 0.5$ MPa. Curves are rather linear as can be seen in Figure 2D, whether indenting near an anticlinal wall or near the middle of a cell. We fitted these curves with linear equations, $F = a\delta$, and found that the prefactor $a$ depends linearly on the elastic modulus (Figure 2E), though this dependence is weak when indenting at the cell middle. Consequently, turgor pressure is the main parameter influencing indentation on periclinal walls, whereas cell wall and turgor contribute with comparable magnitude to these curves near an anticlinal wall. Altogether, it is the ratio of pressure to modulus, $P/E$, that matters for qualitative changes in behavior upon indentation.

Finally, we studied the impact of probe size (Figure 3A) for different values of pressure, when indenting in the middle of a cell (Figure 3B) or near an anticlinal wall (Figure 3C). In all cases, probe size has little influence at small displacements and matters more at higher displacements. At the cell middle, a small probe yields higher forces than a large probe because the cell wall is locally more bent, as can be predicted in simple geometries (Vella et al., 2012b). In contrast, at the cell periphery, a small probe yields a smaller force: This might be explained by a smaller contact with the anticlinal wall, which has a larger mechanical resistance to indentation. At lower or zero pressure, probe size has a minimal influence at the cell's middle, whereas it has a relatively strong influence at cell periphery. Overall, probe size does not affect the qualitative behavior of force-displacement curves when pressure is varied.
A

<figure><img src="image_3.png" /></figure>

B

<figure><img src="image_4.png" /></figure>

C

<figure><img src="image_5.png" /></figure>

D

<figure><img src="image_6.png" /></figure>

E

<figure><img src="image_7.png" /><figcaption>FIGURE 2 | Uniform elasticity: Indentation by a sphere (radius 5 µm) of a flower bud with walls having uniform elastic modulus; sensitivity to modulus.</figcaption></figure>

(A) Section showing fixed cell walls in blue. (B–E) Indentation with no turgor pressure (B,C) and with turgor $P = 0.5$ MPa (D,E). (B,D) Force-displacement curves and power-law fits to the curves, near the cell middle and near an anticlinal wall, shown for three values of elastic modulus $E$. (C,E) Prefactors $a$ of the fits to force-displacement curves as a function of the elastic modulus—note that in (C) the values for cell middle are multiplied by 20. With no pressure the exponent $b$ takes the values 2 near the cell middle and 1 (small displacement) or 1.3 (large) near an anticlinal wall (B,C), whereas with pressure the exponent $b$ is always close to 1 consistently with the linear curves (D,E).

## 3.2. Viscoelasticity

Experiments (Peaucelle et al., 2011; Fernandes et al., 2012; Braybrook and Peaucelle, 2013; Milani et al., 2014; Beauzamy

et al., 2015b; Peaucelle et al., 2015) show hysteresis in force-displacement curves: approach (indentation) and retract (detindentation) do not superimpose. We therefore explored the
<figure><img src="image_8.png" /><figcaption>**FIGURE 3** | Probe size: Indentation by a sphere of a flower bud with walls having uniform elastic modulus; sensitivity to probe size. (A) Two probe sizes are considered: 0.6 and 5 µm; the modulus is fixed at 50 MPa. (B,C) Force-displacement curves obtained near cell middle (B) or near anticlinal wall (C). The values of pressure range from 0.1 MPa to 1 MPa; the lighter colors correspond to smaller values. The red and blue curves correspond to the small and large probe, respectively.</figcaption></figure>

possibility of a linear viscoelastic behavior of the cell walls. We considered a Kelvin-Voigt constitutive law, where stress is the sum of an elastic stress that depends linearly on strain and of a viscous stress that depends linearly on strain rate. This constitutive law describes wall relaxation at time scales much smaller than the time scales of irreversible wall deformation associated with growth. The Kelvin-Voigt law yields a viscoelastic timescale $\tau = \eta/E$, which is the ratio between viscosity $\eta$ and elastic modulus $E$. The rate of force application becomes important when viscosity is accounted for. We considered that the force increases at constant rate up to its maximal value and then decreases to 0, the two phases having the same duration $T$. Dimensional analysis indicates that the results depend on the ratio, $\rho = \tau/T$, of the two times.

We first considered $\rho = 1$, and retrieved hysteresis in the force-displacement curves without (Figure 4A) and with (Figure 4B) turgor pressure, as in experiments; however unlike in experiments, indentation depth slightly increased after the force started to decrease, which might be ascribed to the sudden, unrealistic shift in force rate that we imposed. Hysteresis is smaller in a turgid bud, because pressure contributes significantly to apparent stiffness with no dependence on the velocity of loading. We then varied $\rho$ from small to large values (Figures 4C,D). A higher viscosity always leads to a higher apparent stiffness upon approach. In experiments, the hysteresis is relatively small, suggesting that the viscous time scale is smaller than the force ramp duration and that wall viscosity has a

smaller contribution than wall elasticity to tissue mechanical behavior.

### 3.3. Anisotropy

The cell wall can be mechanically anisotropic when cellulose fibrils are aligned (Cosgrove, 2015). For instance, the elastic modulus of epidermal onion peels is 5 times larger in the main direction of cellulose fibrils than in the perpendicular direction (Kerstens et al., 2001). Based on the orientation of cortical microtubules at the shoot apex (Hamant et al., 2008), surface walls are expected to be mechanically isotropic at the tip of the floral bud and anisotropic at the flanks with the stiffest direction circumferential to the tip. We simulated this configuration with the stiffest direction parallel to the red band shown in Figure 5A, and indented an anisotropic cell from this region, either at its middle or at its periphery. The cell indented had a refined mesh as in the isotropic case. The mean elastic modulus was kept constant, whereas the degree of anisotropy was varied from 1 (isotropic) to 5 (the elastic modulus in one direction is 5 times higher than in the perpendicular direction).

We first considered the case of no pressure. When indenting in the middle (Figure 5B), anisotropy has negligible influence on force-displacement curves at small depth (up to 6% displacement), and anisotropy makes the wall appear slightly softer (by about 5%) at larger depths. When indenting at the cell periphery, anisotropy has no influence at all because anticlinal walls are dominant there (Figure 5C).
<figure><img src="image_9.png" /></figure>

FIGURE 4 | Uniform viscoelasticity: Indentation by a sphere (radius 5 µm) of a flower bud with viscoelastic walls having uniform elastic modulus *E* = 60 MPa and uniform viscosity *η*; sensitivity to the ratio *ρ* = *τ*/*T* of the viscoelastic time scale, *τ* = *η*/*E*, to the duration of the force ramp, *T*; same geometry as in Figure 2A. (A,B) Approach and retract force-displacement curves for *ρ* = 1. (C,D) Approach curves for values of *ρ* = 0, 1/10, 1/3, 1, 3, and 10. (A,C) no pressure and (B,D) turgor pressure *P* = 0.5 MPa.

With pressure, the trends are similar to the plasmolyzed case, though they are somewhat amplified. The changes between the unpressurized case and the pressurized case might be ascribed to the small anisotropic inflation of cells. Mechanical anisotropy makes the wall appear softer (Figures 5D,E). The strongest effect of anisotropy occurs in the cell middle at larger depths: a 5:1 anisotropy makes cells appear softer by 15%.

Overall, mechanical anisotropy of the cell wall slightly reduces
the apparent stiffness of the tissue.

3.4. Variations in Stiffness of Cell Walls
across Layers

Next, we released the assumption that cell wall stiffness is
homogeneous, because there are many indications that, in the
aerial part of the plant, wall stiffness decreases from the surface
toward the inside (see e.g., Beauzamy et al., 2015b). We therefore
subdivided cell walls in 4 categories (Figure 6A): surface walls
(A, elastic modulus E<sub>A</sub>), anticlinal walls in the epidermal L1 cell
layer (B, modulus E<sub>B</sub>), cell walls of the subepidermal L2 layer (C,
modulus E<sub>C</sub>), and all other walls (D, modulus E<sub>D</sub>). We started
with the common value 100 MPa for the elastic modulus and we
gradually reduced the values in the different layers to finally get
to a modulus that is half the initial value in all layers.

We first considered the case of no pressure. In the middle
of the cell (Figure 6B), we observe that we need to reduce
the modulus of anticlinal cell walls by a factor of 2 in order
to reduce the resisting force by about 20%. Surface walls have
a strong influence on force-displacement curves, whereas the
L2 layer and deeper walls have little influence. Near the cell
periphery (Figure 6C), the elastic modulus of anticlinal walls has
the strongest influence on indentation curves, though surface
walls have a comparable influence. D-walls do not contribute
mechanically. Interestingly, reducing L2 cell walls by a factor of
2 can lead to a reduction of the resisting force by about 20%.
This could be expected because the L2 is mechanically closer to
anticlinal cell walls of the L1 than to the middle of surface cell
walls.

In the pressurized case, indentation in the middle of the cell (Figure 6D) is only influenced by the elastic modulus of the cell wall, whereas indentation at the cell periphery (Figure 6E) is influenced by the moduli of both surface and anticlinal cell walls.

4. DISCUSSION

Using a rather realistic mechanical model of a plant tissue
indented by a spherical probe, we explored the behavior of force-
displacement curves when varying probe size, contact location
<figure><img src="image_10.png" /><figcaption>**FIGURE 5** | Mechanical anisotropy: Indentation by a sphere of a flower bud with anisotropic surface walls; sensitivity to mechanical anisotropy. (A) The flanks (red) are mechanically anisotropic; surface walls are stiffer in the circumferential direction (along the red band) and indentation is performed on a cell from the flanks. The radius of the probe is 5 µm. (B–E) Force-displacement curves obtained with zero pressure (B,C) and a pressure P = 0.5 MPa (D,E), near the cell middle (B,D) or above an anticlinal wall (C,E). The anisotropy ratio takes the values 1 : 1 (isotropic), 3 : 1, and 5 : 1. The average elastic modulus is E = 100 MPa.</figcaption></figure>

(near cell center or near anticlinal wall), elastic modulus of cell walls, viscosity of walls, mechanical anisotropy of walls, or turgor pressure. This mechanical model is directly relevant to many recent experimental studies (Peaucelle et al., 2011; Fernandes et al., 2012; Braybrook and Peaucelle, 2013; Milani et al., 2014; Beauzamy et al., 2015b; Peaucelle et al., 2015), though it does not apply to indentation depths significantly smaller than wall thickness (Milani et al., 2011; Hayot et al., 2012; Routier-Kierzkowska et al., 2012; Radotic et al., 2012; Forouzesh et al., 2013).

As a first approximation, indentation curves are mostly sensitive to the elastic modulus of surface walls and to turgor, while mechanical anisotropy has little influence on

these curves. This is consistent with previous theoretical work
in configurations corresponding to single cells (Vella et al.,
2012a,b; Weber et al., 2015). Therefore the main parameter
controlling the qualitative behavior of curves is the ratio of
turgor to average elastic modulus, *P*/*E*, as well as the spatial
variations in this ratio. If cell size, *R*, or wall thickness,
*h*, were varied, work on thin shells (Vella et al., 2012a,b;
Weber et al., 2015) suggests that *PR*/**Eh* should be the
relevant parameter. The vicinity of anticlinal walls appears
stiffer (Peaucelle et al., 2011; Braybrook and Peaucelle, 2013;
Milani et al., 2014; Peaucelle et al., 2015), in the absence of
other features such as a softer material above anticlinal cell
walls in onion epidermis (Routier-Kierzkowska et al., 2012).
<figure><img src="image_11.png" /></figure>

**FIGURE 6** | Spatially varying elasticity: Indentation by a sphere of a flower bud where elastic modulus varies according to layers. (A) Definition of the layers: A-surface walls, B-walls anticlinal to the surface, C-all walls of the second cell layer, and D-all other walls. (B–E) Force-displacement curves obtained with zero pressure (B,C) and a pressure *P* = 0.5 MPa (D,E), near the cell middle (B,D) or above an anticlinal wall (C,E). Elastic moduli (*E*<sub>A</sub>, *E*<sub>B</sub>, *E*<sub>C</sub>, *E*<sub>D</sub>) take the values (100, 100, 100, 100), (100, 100, 100, 50), (100, 100, 50, 50), (100, 50, 50, 50), or (50, 50, 50, 50) MPa, respectively.

We also find that a viscoelastic behavior of the cell walls may account for the hysteresis in force-displacement curves observed in experiments; however, with turgor pressure, such hysteresis might also originate in water movement in the tissue (Beauzamy et al., 2015b).

In turgid tissues, probe size does not affect the results qualitatively. Apparent stiffness in the cell middle is mostly sensitive to turgor, making it possible to infer turgor as performed by Beauzamy et al. (2015b). At the cell periphery, indentation is primarily sensitive to turgor and to the modulus of anticlinal walls and secondarily to surface walls. Milani et al. (2014) ascribed the difference in cell peripheral stiffness between the central zone and peripheral zone of aerial meristems to cell walls,

though they did not exclude variations in turgor. This restriction is relevant based on the present results and future work should address whether there are spatial differences in turgor at the shoot apex.

In plasmolyzed tissues, probe size matters. With a small probe, apparent stiffness is mostly sensitive to walls of the epidermis - surface walls or anticlinal walls according to the location. With a probe comparable to cell size, apparent stiffness of anticlinal walls is also sensitive to the elastic modulus of the sub-epidermal L2 cell layer. This gives theoretical grounds to how Peaucelle et al. (2011) interpreted their results on the difference in stiffness between flower initia and meristem proper: finding no difference with a small probe and different values with a large probe, they
concluded that cell walls of the L2 and L3 soften at initia. Our
results indicate that this conclusion mostly holds for the L2.

An apparent limitation of our study is that we assumed
that cell walls are linearly elastic, knowing that cell walls can
exhibit nonlinear elasticity, with e.g., walls that are stiffer when
under tension (Kierzkowski et al., 2012; Lipchinsky et al., 2013).
Nevertheless, cell walls undergo only small deformations under
indentation at the depths explored here. Given that any nonlinear
material has an elastic modulus for each value of tension,
the moduli considered in the present study for turgid tissues
could be interpreted as the elastic moduli of the cell walls
in the turgid state. Therefore, our work is broadly applicable
to the interpretation of indentation experiments on plant
tissues.

A last question is the relevance of our work to morphogenesis
and more specifically to growth control. Strictly speaking, the
expansion of cell walls is determined by their extensibility
and not by their elasticity (Cosgrove, 2015). However, studies
of mechanics of the shoot apex (as reviewed in Milani
et al., 2013) suggest a correlation between elasticity and
extensibility. In this context, our results can help constraining
three-dimensional models of growing plant tissues (Boudon
et al., 2015) that address the mechanical regulation of
morphogenesis.

REFERENCES

Ali, O., Mirabet, V., Godin, C., and Traas, J. (2014). Physical models of plant
development. *Annu. Rev. Cell Dev. Biol.* 30, 59–78. doi: 10.1146/annurev-
cellbio-101512-122410

Bassel, G. W., Stamm, P., Mosca, G., de Reuille, P. B., Gibbs, D. J., Winter, R.,
et al. (2014). Mechanical constraints imposed by 3d cellular geometry and
arrangement modulate growth patterns in the arabidopsis embryo. *Proc. Natl.
Acad. Sci. U.S.A.* 111, 8685–8690. doi: 10.1073/pnas.1404616111

Beauzamy, L., Derr, J., and Boudaoud, A. (2015a). Quantifying hydrostatic pressure
in plant cells by using indentation with an atomic force microscope. *Biophys. J.*
108, 2448–2456. doi: 10.1016/j.bpj.2015.03.035

Beauzamy, L., Louveaux, M., Hamant, O., and Boudaoud, A. (2015b). Mechanically, the shoot apical meristem of arabidopsis behaves like a shell inflated by a pressure of about 1 mpa. *Front. Plant Sci.* 6:1038. doi: 10.3389/fpls.2015.01038

Beauzamy, L., Nakayama, N., and Boudaoud, A. (2014). Flowers under pressure:
ins and outs of turgor regulation in development. *Ann. Bot.* 114, 1517–1533.
doi: 10.1093/aob/mcu187

Bonilla, M. R., Stokes, J. R., Gidley, M. J., and Yakubov, G. E. (2015). Interpreting atomic force microscopy nanoindentation of hierarchical biological materials using multi-regime analysis. *Soft Matter* 11, 1281–1292. doi: 10.1039/C4SM02440K

Boudaoud, A. (2010). An introduction to the mechanics of morphogenesis for plant biologists. *Trends Plant Sci.* 15, 353–360. doi: 10.1016/j.tplants.2010.04.002

Boudon, F., Chopard, J., Ali, O., Gilles, B., Hamant, O., Boudaoud, A., et al. (2015). A computational framework for 3D mechanical modeling of plant morphogenesis with cellular resolution. *PLoS Comput. Biol.* 11:e1003950. doi: 10.1371/journal.pcbi.1003950

Braybrook, S. A., and Peaucelle, A. (2013). Mechano-chemical aspects of organ
formation in *Arabidopsis thaliana*: the relationship between auxin and pectin.
*PLoS ONE* 8:e57813. doi: 10.1371/journal.pone.0057813

Cho, H. T., and Cosgrove, D. J. (2000). Altered expression of expansin modulates
leaf growth and pedicel abscission in arabidopsis thaliana. *Proc. Natl. Acad. Sci.
U.S.A.* 97, 9783–9788. doi: 10.1073/pnas.160276997

AUTHOR CONTRIBUTIONS

RM, FF, and AB designed research; RM performed research; AB
supervised the project.

FUNDING

The work was funded by the INRIA Morphogenetics project, by
the ANR (#12-BSV2-023-02), and by the ERC (Starting Grant
Phymorph #307387).

ACKNOWLEDGMENTS

We thank the whole Sofa team for their everyday support,
especially Benjamin Gilles for outstanding plugins, P. Das, J.
Chopard, and F. Boudon for the realistic structure of meristem,
O. Ali for suggestions about mechanical anisotropy, and C. Godin
and J. Traas for initiating the morphogenetics project.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found
online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.
01351

Coen, E., Rolland-Lagan, A.-G., Matthews, M., Bangham, J. A., and Prusinkiewicz,
P. (2004). The genetics of geometry. *Proc. Natl. Acad. Sci. U.S.A.* 101, 4728–
4735. doi: 10.1073/pnas.0306308101

Cosgrove, D. J. (2015). Plant cell wall extensibility: connecting plant cell growth with cell wall structure, mechanics, and the action of wall-modifying enzymes. *J. Exp. Bot.* 67, 463–476. doi:10.1093/jxb/erv511

Faure, F., Duriez, C., Delingette, H., Allard, J., Gilles, B., Marchesseau, S., et al. (2012). "SOFA: a multi-model framework for interactive physical simulation," in *Soft Tissue Biomechanical Modeling for Computer Assisted Surgery, volume 11 of Studies in Mechanobiology, Tissue Engineering and Biomaterials*, ed Y. Payan (Berlin: Springer), 283–321.

Fernandes, A. N., Chen, X., Scotchford, C. A., Walker, J., Wells, D. M., Roberts, C. J., et al. (2012). Mechanical properties of epidermal cells of whole living roots of *Arabidopsis thaliana*: an atomic force microscopy study. *Phys. Rev. E.* 85:021916. doi: 10.1103/PhysRevE.85.021916

Fernandez, R., Das, P., Mirabet, V., Moscardi, E., Traas, J., Verdeil, J.-L., et al. (2010). Imaging plant growth in 4d: robust tissue reconstruction and lineaging at cell resolution. *Nat. Methods* 7, 547–553. doi: 10.1073/pnas.0306308101

Fleming, A. J., McQueen-Mason, S., Mandel, T., and Kuhlemeier, C. (1997). Induction of leaf primordia by the cell wall protein expansin. *Science* 276, 1415–1418. doi: 10.1126/science.276.5317.1415

Forouzesh, E., Goel, A., Mackenzie, S. A., and Turner, J. A. (2013). *In vivo* extraction of arabidopsis cell turgor pressure using nanoindentation in conjunction with finite element modeling. *Plant J.* 73, 509–520. doi: 10.1111/tpj.12042

Geitmann, A. (2006). Experimental approaches used to quantify physical
parameters at cellular and subcellular levels. *Am. J. Bot.* 93, 1380–1390. doi:
10.3732/ajb.93.10.1380

Geitmann, A., and Parre, E. (2004). The local cytomechanical properties of growing pollen tubes correspond to the axial distribution of structural cellular elements. *Sexual Plant Reproduct.* 17, 9–16. doi: 10.1007/s00497-004-0210-3

Hamant, O., Heisler, M. G., Jönsson, H., Krupinski, P., Uyttewaal, M., Bokov, P., et al. (2008). Developmental patterning by mechanical signals in arabidopsis. *Science* 322, 1650–1655. doi: 10.1126/science.1165594

Hayot, C. M., Forouzesh, E., Goel, A., Avramova, Z., and Turner, J. A. (2012). Viscoelastic properties of cell walls of single living plant cells determined
by dynamic nanoindentation. *J. Exp. Bot.* 63, 2525–2540. doi: 10.1093/jxb/err428

Johnson, K. (1987). *Contact Mechanics*. Cambridge: Cambridge University Press.

Kerstens, S., Decraemer, W. F., and Verbelen, J.-P. (2001). Cell walls at the plant surface behave mechanically like fiber-reinforced composite materials. *Plant Physiol.* 127, 381–385. doi: 10.1104/pp.010423

Kierzkowski, D., Nakayama, N., Routier-Kierzkowska, A.-L., Weber, A., Bayer, E., Schorderet, M., et al. (2012). Elastic domains regulate growth and organogenesis in the plant shoot apical meristem. *Science* 335, 1096–1099. doi: 10.1126/science.1213100

Lee, D., Barber, J. R., and Thouless, M. D. (2009). Indentation of an elastic half space with material properties varying with depth. *Int. J. Eng. Sci.* 47, 1274–1283. doi: 10.1016/j.ijengsci.2008.08.005

Lintilhac, P. M., Wei, C., Tanguay, J. J., and Outwater, J. O. (2000). Ball tonometry: a rapid, nondestructive method for measuring cell turgor pressure in thin-walled plant cells. *J. Plant Growth Regul.* 19, 90–97. doi: 10.1007/s00344000009

Lipchinsky, A., Sharova, E., and Medvedev, S. (2013). Elastic properties of the growth-controlling outer cell walls of maize coleoptile epidermis. *Acta Physiol. Plant.* 35, 2183–2191. doi: 10.1007/s11738-013-1255-4

Malgat, R., Boudaoud, A., and Faure, F. (2014). “Mechanical modeling of three-dimensional plant tissue indented by a probe,” in *Workshop on Virtual Reality Interaction and Physical Simulation*, eds J. Bender, C. Duriez, F. Jaillet, and G. Zachmann (Geneva: The Eurographics Association), 59–68. doi: 10.2312/vriphys.20141224

Milani, P., Braybrook, S. A., and Boudaoud, A. (2013). Shrinking the hammer: micromechanical approaches to morphogenesis. *J. Exp. Bot.* 64, 4651–4662. doi: 10.1093/jxb/ert619

Milani, P., Gholamirad, M., Traas, J., Arnéodo, A., Boudaoud, A., Argoul, F., et al. (2011). *In vivo* analysis of local wall stiffness at the shoot apical meristem in arabidopsis using atomic force microscopy. *Plant J.* 67, 1116–1123. doi: 10.1111/j.1365-313X.2011.04649.x

Milani, P., Mirabet, V., Cellier, C., Rozier, F., Hamant, O., Das, P., et al. (2014). Matching patterns of gene expression to mechanical stiffness at cell resolution through quantitative tandem epifluorescence and nanoindentation. *Plant Physiol.* 165, 1399–1408. doi: 10.1104/pp.114.237115

Mirabet, V., Das, P., Boudaoud, A., and Hamant, O. (2011). The role of mechanical forces in plant morphogenesis. *Ann. Rev. Plant Biol.* 62, 365–385. doi: 10.1146/annurev-arplant-042110-103852

Moulia, B. (2013). Plant biomechanics and mechanobiology are convergent paths to flourishing interdisciplinary research. *J. Exp. Bot.* 64, 4617–4633. doi: 10.1093/jxb/ert320

Niklas, K. (1992). *Plant Biomechanics, An Engineering Approach to Plant Form and Function*. Chicago, IL: The University of Chicago Press.

Nilsson, S. B., Hertz, C. H., and Falk, S. (1958). On the relation between turgor pressure and tissue rigidity. *2. Physiol. Plant.* 11, 818–837. doi: 10.1111/j.1399-3054.1958.tb08275.x

Peaucelle, A., Braybrook, S. A., Guillou, L. L., Bron, E., Kuhlemeier, C., and Höfte, H. (2011). Pectin-induces changes in cell wall mechanics underlie organ initiation in *Arabidopsis*. *Curr. Biol.* 21, 1720–1726. doi: 10.1016/j.cub.2011.08.057

Peaucelle, A., Louvet, R., Johansen, J. N., Höfte, H., Laufs, P., Pelloux, J., et al. (2008). Arabidopsis phyllotaxis is controlled by the methyl-esterification status of cell-wall pectins. *Curr. Biol.* 18, 1943–1948. doi: 10.1016/j.cub.2008.10.065

Peaucelle, A., Wightman, R., and Höfte, H. (2015). The control of growth symmetry breaking in the *Arabidopsis* hypocotyl. *Curr. Biol.* 25, 1746–1752. doi: 10.1016/j.cub.2015.05.022

Pien, S., Wyrzykowska, J., McQueen-Mason, S., Smart, C., and Fleming, A. (2001). Local expression of expansin induces the entire process of leaf development and modifies leaf shape. *Proc. Natl. Acad. Sci. U.S.A.* 98, 11812–11817. doi: 10.1073/pnas.191380498

Radotic, K., Roduit, C., Simonovic, J., Hornitschek, P., Frankhauser, C., Mutavžić, D., Steinbach, G., et al. (2012). Atomic force microscopy stiffness tomography

on living arabidopsis thaliana cells reveals the mechanical properties of surface and deep cell-wall layers during growth. *Biophys. J.* 103, 386–394. doi: 10.1016/j.bpj.2012.06.046

Robinson, S., Burian, A., Couturier, E., Landrein, B., Louveaux, M., Neumann, E. D., et al. (2013). Mechanical control of morphogenesis at the shoot apex. *J. Exp. Bot.* 64, 4729–4744. doi: 10.1093/jxb/ert199

Roduit, C., Sekatski, S., Dietler, G., Catsicas, S., Lafont, F., and Kasas, S. (2009). Stiffness tomography by atomic force microscopy. *Biophys. J.* 97, 674–677. doi: 10.1016/j.bpj.2009.05.010

Routier-Kierzkowska, A.-L., Weber, A., Kochova, P., Felekis, D., Nelson B. J., Kuhlemeier, C., et al. (2012). Cellular force microscopy for *in vivo* measurements of plant tissue mechanics. *Plant Physiol.* 158, 1514–1522. doi: 10.1104/pp.111.191460

Routier-Kierzkowska, A.-L., and Smith, R. S. (2013). Measuring the mechanics of morphogenesis. *Curr. Opin. Plant Biol.* 16, 25–32. doi: 10.1016/j.pbi.2012.11.002

Sampathkumar, A., Krupinski, P., Wightman, R., Milani, P., Berquand, A., Boudaoud, A., et al. (2014). Subcellular and supracellular mechanical stress prescribes cytoskeleton behavior in *Arabidopsis* cotyledon pavement cells. *eLife* 3:e01967. doi: 10.7554/eLife.01967

Schopfer, P. (2006). Biomechanics of plant growth. *Am. J. Bot.* 93, 1415–1425. doi: 10.3732/ajb.93.10.1415

Tomos, A. D., and Leigh, R. A. (1999). The pressure probe: a versatile tool in plant cell physiology. *Ann. Rev. Plant Physiol. Plant Mol. Biol.* 50:447.

Traas, J., and Monéger, F. (2010). Systems biology of organ initiation at the shoot apex. *Plant Physiol.* 152, 420–427. doi: 10.1104/pp.109.150409

Tvergaard, V., and Needleman, A. (2011). Polymer indentation: numerical analysis and comparison with a spherical cavity model. *J. Mech. Phys. Solids* 59, 1669–1684. doi: 10.1016/j.jmps.2011.06.006

Valero, C., Navarro, B., Navajas, D., and García-Aznar, J. M. (2016). Finite element simulation for the mechanical characterization of soft biological materials by atomic force microscopy. *J. Mech. Behav. Biomed. Mater.* 62, 222–235. doi: 10.1016/j.jmbbm.2016.05.006

Vella, D., Ajdari, A., Vaziri, A., and Boudaoud, A. (2012a). Indentation of ellipsoidal and cylindrical elastic shells. *Phys. Rev. Lett.* 109:144302. doi: 10.1103/PhysRevLett.109.144302

Vella, D., Ajdari, A., Vaziri, A., and Boudaoud, A. (2012b). The indentation of pressurized elastic shells: from polymeric capsules to yeast cells. *J. R. Soc. Interf.* 9, 448–455. doi: 10.1098/rsif.2011.0352

Vogler, H., Draeger, C., Weber, A., Felekis, D., Eichenberger, C., Routier-Kierzkowska, A.-L., et al. (2013). The pollen tube: a soft shell with a hard core. *Plant J.* 73, 617–627. doi: 10.1111/tpj.12061

Vogler, H., Felekis, D., Nelson, B. J., and Grossniklaus, U. (2015). Measuring the mechanical properties of plant cell walls. *Plants* 4, 167–182. doi: 10.3390/plants4020167

Weber, A., Braybrook, S., Huflejt, M., Mosca, G., Routier-Kierzkowska, A.-L., and Smith, R. S. (2015). Measuring the mechanical properties of plant cells by combining micro-indentation with osmotic treatments. *J. Exp. Bot.* 66, 3229–3241. doi: 10.1093/jxb/erv135

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2016 Malgat, Faure and Boudaoud. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.