# Mechanosensitive control of plant growth: Bearing the load, sensing, transducing and responding

Bruno B. Moulia, Catherine Coutand, Jean-Louis J.-L. Julien

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

## To cite this version:

Bruno B. Moulia, Catherine Coutand, Jean-Louis J.-L. Julien. Mechanosensitive control of plant growth: Bearing the load, sensing, transducing and responding. Frontiers in Plant Science, 2015, 6, 20 p. 10.3389/fpls.2015.00052. hal-01123388

HAL Id: hal-01123388

https://hal.science/hal-01123388v1

Submitted on 28 May 2020

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_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

# Mechanosensitive control of plant growth: bearing the load, sensing, transducing, and responding

Bruno Moulia<fnref n="1" />,<fnref n="2" />*†, Catherine Coutand<fnref n="1" />,<fnref n="2" />† and Jean-Louis Julien<fnref n="1" />,<fnref n="2" />

**Edited by:**
Burkhard Schulz, Purdue University, USA

**Reviewed by:**
Wei-Hua Tang, Chinese Academy of Sciences, China
Nabil I. Elsheery, Tanta University, Egypt

**\*Correspondence:**
Bruno Moulia, UMR, PIAF
Integrative Physics and Physiology of Trees, Institut National de la Recherche Agronomique, 5 chemin de Beaulieu, F63039
Clermont-Ferrand, France
e-mail: bruno.moulia@clermont.inra.fr

As land plants grow and develop, they encounter complex mechanical challenges, especially from winds and turgor pressure. Mechanosensitive control over growth and morphogenesis is an adaptive trait, reducing the risks of breakage or explosion. This control has been mostly studied through experiments with artificial mechanical loads, often focusing on cellular or molecular mechanotransduction pathway. However, some important aspects of mechanosensing are often neglected. (i) What are the mechanical characteristics of different loads and how are loads distributed within different organs? (ii) What is the relevant mechanical stimulus in the cell? Is it stress, strain, or energy? (iii) How do mechanosensing cells signal to meristematic cells? Without answers to these questions we cannot make progress analyzing the mechanobiological effects of plant size, plant shape, tissue distribution and stiffness, or the magnitude of stimuli. This situation is rapidly changing however, as systems mechanobiology is being developed, using specific biomechanical and/or mechanobiological models. These models are instrumental in comparing loads and responses between experiments and make it possible to quantitatively test biological hypotheses describing the mechanotransduction networks. This review is designed for a general plant science audience and aims to help biologists master the models they need for mechanobiological studies. Analysis and modeling is broken down into four steps looking at how the structure bears the load, how the distributed load is sensed, how the mechanical signal is transduced, and then how the plant responds through growth. Throughout, two examples of adaptive responses are used to illustrate this approach: the thigmorphogenetic syndrome of plant shoots bending and the mechanosensitive control of shoot apical meristem (SAM) morphogenesis. Overall this should provide a generic understanding of systems mechanobiology at work.

**Keywords:** mechanobiology, biomechanics, thigmomorphogenesis, wind, turgor pressure, curvature, mechanotransduction, stress

## INTRODUCTION

Land plants continuously encounter mechanical challenges from without and within. External mechanical loads are imposed by the wind, rain, neighboring plants or solid substrates. The external bending loads imposed by winds induce a syndrome of mechanosensitive growth responses in the aerial stems of plants known as thigmomorphogenesis. The activity of their meristems is modulated to stunt vertical growth and stimulate an increase in girth, thereby making the plant more wind-resistant (see Telewski, 2006; Coutand, 2010; Monshausen and Haswell, 2013 for reviews). Internal loads may be imposed by the plant's own weight, inertial forces and the large hydro-static turgor pressure in cells. Even meristems, although protected from many external mechanical loads by young leaves in the shoot apical bud or by the bark in the lateral camium, are under considerable direct mechanical stress due to the inner turgor pressure and the mechanical barriers imposed by neighboring organs or tissues (e.g., Couturier et al., 2012; Baskin and Jensen, 2013). Therefore, precise mechanosensitive

control of the magnitude and direction of growth is required so that the size, shape, and edges of the growing organs and tissues are produced in a regular and stable way (Hamant, 2013). It follows that acclimation responses of growth and morphogenesis have been naturally selected to reduce the risk of breakage or explosion of plant parts during growth and development.

These two adaptive responses, stem thigmomorphogenesis and meristem growth, ultimately rely on mechanosensing of the internal mechanical state of the living cells of the plant as a cue for the regulation of growth and morphogenesis (Coutand, 2010; Moulia et al., 2011; Hamant, 2013; Monshausen and Haswell, 2013). Mechanosensing occurs at the cell level, yet mechanical stimulation involves loads that act on the whole organ, either at its boundaries (e.g., for wind-drag) or across its full volume (e.g., weight, inertial forces or turgor pressure). Therefore, changes in the mechanical state of tissues and cells that trigger cell mechanosensing depend on the load, on the mechanical structure of the organ, and on the mechanosensitive structure.

<footnote n="1">¹ NRA, UMR 547 PIAF, Clermont-Ferrand, France</footnote>

<footnote n="2">² Clermont Université, Université Blaise Pascal, UMR 547 PIAF, Clermont-Ferrand, France</footnote>

<footnote>† These authors have contributed equally to this work.</footnote>
The mechanosensitive structure is defined as the location and amount of mechanosensitive tissues involved in a response. Some mechanosensored modulations of growth and morphogenesis are triggered through long-distance internal signaling so the connection between the mechanosensitive structure and the responding structure also needs to be borne in mind (Coutand, 2010; Moulia et al., 2011).

Analyzing and modeling the biology of mechanosensing and of mechanosensitive growth responses thus involves three phases (Moulia et al., 2011). (i) Biomechanical analysis reveals how mechanical loads are distributed over the constitutive plant tissues and cells. (ii) The local mechanosensitive pathways are analyzed in the sensing cells. (iii) Mechanobiological integration combines the models describing how the plant's local mechanosensing relates to global growth responses. Other comprehensive reviews have focused on the local analysis of the mechanosensitive pathway or on the global responses of growth and morphogenesis (e.g., Braam, 2005; Telewski, 2006; Coutand, 2010; Monshausen and Haswell, 2013 to cite just a few). Our purpose instead is to review the integrative aspects, tracing them down the scale from the effect of the load on the plant to the effects on tissue elements and cells, and then up the scale from mechanosensitive gene expression to the growth and morphogenetic responses of the organ. Two mechanosensitive growth responses have been particularly extensively studied in the last two decades: thigmomorphogenesis of stems responding to external bending loads, and growth and morphogenesis of the shoot apical meristem (SAM). In particular, we aim to illustrate how integrative models combining structural mechanics with mechanosensory biology have been instrumental in understanding how mechanical loads are distributed within the plant, defining the heterogeneous stress and strain fields. We explain how the models become key experimental tools to qualitatively and quantitatively assess hypotheses about sensory mechanisms (e.g., does sensing occur through stretch-activated channels or wall-associated transmembrane proteins? Is strain sensed or is stress sensed or both?) or the influence of organ geometry and tissue distribution on the magnitude of mechanosensitive responses. This review is designed for a general biologist audience and aims to help biologists master the mechanical models they need for mechanobiological studies. There is no need for an advanced background in mechanics, mathematics or modeling as the crucial equations are introduced both verbally and graphically. The list of the models and of their acronyms can be found in Table 1.

## MECHANICAL CHARACTERISTICS OF LOADS AND THEIR HETEROGENEOUS DISTRIBUTION WITHIN THE PLANT

Some central concepts are introduced briefly here that are essential when taking a mechanical view of plant structure. More complete primers in plant mechanics, including lists of definitions, are available (Boudaoud, 2010; Moulia et al., 2011; Moulia, 2013).

## CRUCIAL MECHANICAL CONCEPTS

Mechanical loads may involve the action of localized forces (e.g., intermittent contact with a neighboring stone or organ) or

<table><caption>Table 1 | List of models with their acronyms and references.</caption><thead><tr><td>Acronym of the model</td><td>Full name</td><td>References</td></tr></thead><tbody><tr><td>CBmS</td><td>Composite beam model of the Stem (in flexion)</td><td>Moulia and Fournier, 1997; Gibson et al., 1988; Coutand and Moulia, 2000</td></tr><tr><td>PVm</td><td>Pressurized vessel model of the SAM</td><td>Hamant et al., 2008; Traas and Hamant, 2009</td></tr><tr><td>FEm</td><td>Finite elements model of two patches of the L1 + L2 tissues</td><td>Hamant et al., 2008</td></tr><tr><td>2D SFm</td><td>Two-dimensional cellular stress feedback model</td><td>Hamant et al., 2008</td></tr><tr><td>S³m</td><td>Sum of strain-sensing model</td><td>Coutand and Moulia, 2000; Coutand et al., 2009; Moulia et al., 2011</td></tr><tr><td>SAM SFm</td><td>Integrative SAM stress feedback model ()</td><td>Hamant et al., 2008</td></tr><tr><td>GSFm</td><td>Growth-strain feedback model</td><td>see replies to Hamant et al., 2008 by Schopfer and Meyerowitz in Science e-letters</td></tr></tbody></table>

distributed loads (e.g., external drag by wind flow, self-weight or the internal turgor pressure of the living cells). Some loads are static or quasi-static (their rate of change is slow), whereas others are dynamic so inertial forces due to the acceleration of mass need to be considered (e.g., wind-induced oscillations) (Rodriguez et al., 2008; Pivato et al., 2014).

Under the action of internal and/or external load(s) a solid body such as a plant organ can be globally displaced. This displacement involves a translation of the center of mass of the body and the body might rotate around the center of mass, described in terms of velocity. In addition, parts of the body may be displaced relative to one another, resulting in a change of shape, called deformation. These deformations are measured locally by strains, written as $\epsilon$. Strains may be linear or shear (angular) and are measured in relative units (i.e., strains are dimensionless). Straining stretches bonds and causes slide/shear of internal elements, thereby allowing internal reaction forces to build up. In this way the mechanical load is distributed through the material with the storage of elastic strain energy across the deformed body until an internal and external mechanical equilibrium is achieved, i.e., all the forces and moments acting on the body are balanced. The density of the resulting internal forces, i.e., the internal forces per unit of area, is called stress, written as $\sigma$ and measured in Pascal (Pa) which is equivalent to $N.m^{-2}$. Strains and stresses can be very heterogeneous across the body. Deformation is characterized by the strain and stress fields, i.e., the amount of strain and stress at every location of the body at a given time. The
strain and stress fields in a given load situation therefore, measure the mechanical state of a body like a cell, an organ or a plant.

The amount of stress produced by straining is linked to the rheological properties of the material. Rheology can be modeled in a so-called constitutive equation. If the stress increases linearly with increasing strain and linearly reverts during unloading, the material behaves as a linear elastic material. The slope of the stress-strain curve is called the Young's modulus, such that a stiff, rigid material has a high Young's modulus. Over a certain threshold, some materials may yield plastically (irreversibly). From this point, in pure plastic materials, stress does not increase further with strain, but in visco-plastic materials it varies depending on the strain rate. The growing cell wall has a visco-plastic rheology (called the Bingham flow model, see Dumais, 2013 for more details). Finally, when internal forces overcome the strength of the material, fracture occurs.

## POSING THE PROBLEM: EXPERIMENTAL SETUPS

Thigmomorphogenetic experiments are generally conducted in the lab by subjecting single plant stems to static bending. Typically a stem is bent by moving the top of the stem laterally while the base is anchored immobile in the soil (Telewski and Pruyn, 1998). Alternatively, the stem is first clamped in a vertical position with the roots bathed in a liquid medium, then the basal part of the stem is displaced (Figure 1). The latter setup decouples the effects of stem bending from tilting the apical growth zone (Coutand et al., 2000). Dynamic loading can also be imposed by vibrating the plant, for example (Der Loughian et al., 2014).

For most SAM experiments, the meristem is cut from the stem and the surrounding young leaves are removed (Figure 2). The isolated SAM is then cultured on a growth medium. Three types of mechanical perturbations have been used on excised SAM. The osmotic potential of the bathing solution can be changed to transiently manipulate the inner turgor pressure of the cells (Peaucelle et al., 2011). External loads, such as lateral compression of the whole meristem, can be applied mechanically (Hamant et al., 2008). Alternatively, outer cells can be ablated to create holes, thereby modifying the mechanical structure and state of the SAM (Hamant et al., 2008). To our knowledge, the mechanical states of intact SAM have not been studied so far.

In both experimental systems, the next step is to estimate the amount and spatial distribution of the changes in mechanical state (stress and strains). This can be done as (i) the external mechanical load and the mechanical structure of the organ are known, and (ii) the changes imposed by the experimenter or by the environmental conditions are measured. However, to estimate changes that occur at very different scales, we need to consider how a change in a unit of the cell wall or tissue affects the whole organ (and vice versa). For this, a mechanical model needs to be developed, using a scientific method originating from mechanical engineering called integrative structural mechanics (ISM) modeling. More complete coverage of ISM modeling can be found in Moulia et al. (2011).

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

**FIGURE 1** | Morphological and anatomical structure of a stem submitted to an external bending load from Coutand and Moulia (2000), Journal of Experimental Botany, by permission of the Society for Experimental Biology. (A) Side view of the basal part of the stem base before the application of bending. The stem is grown in hydroponics, and clamped below the primary growth zone, so that bending does not affect its position (Cl, metal clamp; IN 1-3, internodes; Hyp, hypocotyl; Cot, cotyledons. (B) Idealized geometrical scheme (A) as a cantilever beam. (C) Negative photograph of a cross-section (note the quasi-circular shape and the concentric layers of tissues (E, epidermis; Co, collenchyma; Pa, parenchyma; Ph.2, metaphloem; Ca, cambium; Xyl.2, metaxylem; Xyl.1, protoxylem; Ph.i, internal phloem). (D) Changes in the external diameter (●) and of the diameter of the pith (■) along the basal part of the stem.

## INTEGRATIVE MODELS IN MECHANICS

The general structure of an ISM model is shown in Figure 3, using schematic graphical conventions that will be used throughout this review. The major aspects of the models can be understood independently of the detailed model equations. First, the constitutive materials of the structure, the elementary “bricks” or units, are defined and the rheological properties of these elements specified. Do they behave elastic or do they undergo viscoplastic deformations? Are they isotropic, displaying the same properties in all directions? If they are anisotropic, in which direction is the anisotropy? (Coutand and Moulia, 2000; Hamant et al., 2008; Baskin and Jensen, 2013). It is very important to know the shape of these elements at rest (without a load) and whether the shape is dependent on other physical variables such as temperature, water status or time (e.g., Moulia, 2000; Hamant et al., 2008). The size of these elements is not prescribed as the mechanical structure of plants is multiscale and the scale at which to work is mostly a matter of informed choice (Boudaoud, 2010; Niklas and Spatz, 2012; Gibson, 2013). Next, the structure is defined by specifying how the elements are assembled (topology) and displayed geometrically. Note that to model heterogeneous organs, several materials may need to be considered (e.g., Moulia and Fournier, 1997; Coutand and Moulia, 2000; Routier-Kierzkowska et al.,
<figure><img src="image_5.png" /></figure>

**FIGURE 2** | Structure of the shoot apical meristem (SAM). (A) View of an *Arabidopsis thaliana* SAM from above Hamant et al. (2008), reprinted with permission from AAAS. (B) Side view of a tomato SAM (Robinson et al., 2013). (C-E) Generalized schematic representations of a typical dome-shaped shoot apex bearing a cylindrical young primordium. (C) Major structures. M, shoot apical meristem; P, organ primordium; B, boundary between the meristem and the primordium. (D) Morphological domains of the SAM. CZ, central zone, PZ, peripheral zone where new organs are generated. (E) Internal organization of the SAM. L1, presumed epidermis, L1 and L2, tunica layers, L3, corpus from Robinson et al. (2013), Journal of Experimental Botany, by permission of the Society for Experimental Biology.

2012). Finally, the mechanical loading applied to the structure is defined, as well as including any boundary conditions, which are displacement or force constraints at the boundaries of the structure.

With these three steps the model is now fully defined. When input values are known like the load applied or the structural change (e.g., making a hole in the structure), the state of the structure in the loaded state can be computed because mechanical laws specify (i) the conditions for equilibrium (static or dynamic), and (ii) the compatibility of strains between adjacent material elements or boundaries. Depending on the structure's geometry, some simplifying hypotheses can be used for calculations, e.g., beam or shell theories. In some cases the problem can even be solved analytically (e.g., Hamant et al., 2008). Mostly, however, numerical methods are required for computations. The outputs of such models can be multiple: knowledge of strain and stress fields, the velocity of the top of a plant, or bending rigidity, etc.

Plants are open systems. If cells grow or differentiate the amount and/or rheology of constitutive materials may change (e.g., cell wall maturation) and will need to be accounted for in a model. This has important implications in formulating the mechanical problems that are specific to biomechanical models (Moulia and Fournier, 2009). For example, for the tree gravitropic reaction, the problem can be solved by using beam theory hypotheses but requires an incremental formulation of the problem (e.g., Fournier et al., 1994; Fourcaud et al., 2008; Coutand et al., 2011) to take into account the growth and shrinkage of the cell wall rest-length during wood maturation (Coutand et al., 2011; Pot et al., 2014).

## LOAD DISTRIBUTION FROM THE PLANT TO MECHANOSENSITIVE CELLS

Two examples of analyzing load distribution from the scale of the whole plant down to the scale of mechanosensitive cells will be presented.

## COMPOSITE BEAM MODEL OF THE STEM SUBJECTED TO BENDING

The dicot stem is composed of several tissues of very variable stiffness, e.g., epidermis, parenchyma, sclerenchyma, and wood. Growth activity is concentrated in (i) the primary growth zone just below the SAM and (ii) the cambial zone, a thin shell of 1–20 layers of meristematic cells near the lateral surface of the stem, just beneath the bark. The primary growth zone is responsible for longitudinal growth, and the cambial zone for most radial expansion.

As the stem is generally a slender structure (the diameter:length ratio is less than 1/20) (see Figure 1) and its constitutive tissues are in transverse layers, the mechanical analysis can be simplified using the theory of heterogeneous composite beams or rods (Gibson et al., 1988; Moulia and Fournier, 1997), reframed in a mechanobiological context (Coutand and Moulia, 2000), and called the Composite Beam model of the Stem (CBmS) in the following.

This mechanical modeling is detailed step by step in **Figure 4**. Only longitudinal strains and stresses will be considered, noted by the subscript LL, as transverse shearing can be neglected when analyzing slender structures, which exhibit pure bending. The material element in the CBmS is a small volume of tissue. Two types of tissues, broadly organized into three concentric rings, were defined. Tissues such as parenchyma or phloem were treated as compliant materials, and tissues like wood as stiff materials. These elements are assumed to behave in the linear elastic range as has been confirmed experimentally (Coutand et al., 2000). These tissue elements are then assembled into infinitesimal slices of dS thickness according to the known anatomy of the stems. Finally, the stem can be viewed as a pile of infinitesimal slices, glued one next to another along an imaginary line inside the stem, called the neutral line.

During bending experiments one end of the stem is fixed and one end is free to move. The stem therefore, operates mechanically as a cantilever subjected to a local force F̄ transverse to the stem. The action of F̄ depends on the amount of the force F and on the lever arm L, i.e., the distance from the application point of the force to a given slice. This mechanical amplification effect can be modeled using a quantity called the bending moment M̂, the magnitude of which is M = F.L (Equation 1).

A central property of beam bending is that each cross section remains flat and orthogonal to the neutral line all along the deformed beam. A change in M̂ will thus induce a relative rotation (through an angle dθ) of two successive stem slice cross-sections. The effect of the rotation is an increment of length dl on the tensed side and a decrement of dl on the compressed side. The ratio C = $\frac{d\theta}{ds}$ (Equation 2) is called the curvature. It measures the spatial density (rate) of bending rotation and $\epsilon_{LL} = \frac{dl}{ds}$ (Equation 3) measures the longitudinal strain. The strain at any point located at a distance y from the cross-section center can
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | **Structure of an integrative structural mechanics (ISM) model from Moulia et al. (2011), by permission of Springer-Verlag Berlin Heidelberg.** The structure of an ISM model for use in plant biomechanics. ISM models consider (at least) two scales in the system: a scale of phenomenological empiricism called the material scale, and a scale of mechanistic spatial integration,</figcaption></figure>

the mechanical structure. The internal and boundary loads (inputs) result in a change in mechanical state that can be calculated using mechanical principles and robust simplifying theories. ISM models can produce various outputs characterizing mechanical state or dynamics, such as strain (ε) and stress (σ) fields, vibration modes, or rupture risk factors.

be computed as the product of the change in curvature and the distance y to the central line of the stem, ε<sub>LL,y</sub> = y·(C − C<sub>0</sub>) (Equation 4), where C<sub>0</sub> is the initial stem curvature before the load (if the initial configuration of the stem is straight). The value of the longitudinal strain thus varies with the position along the beam, being maximal at the beam periphery (y = R) along a radius aligned with ⃗F.

Straining allows internal reaction forces to build up to balance the effect of the external load. For elastic constitutive materials, the stress is calculated as the strain multiplied by the appropriate Young's modulus: σ<sub>LL,x,y,i</sub> = ε<sub>LL,x,y</sub>E<sub>LL,i</sub> (Equation 5) where E<sub>LL,i</sub> is the longitudinal elastic modulus of the i<sup>th</sup> material and x, y are the spatial coordinates within the cross section.

The amount and distribution of stresses and strains can then be calculated so that they balance out M̄ (as detailed in supplemental data). This yields ε<sub>LL</sub> = y(C − C<sub>0</sub>) = My/(E<sub>soft</sub>I<sub>1</sub> + E<sub>stiff</sub>I<sub>2</sub> + E<sub>soft</sub>I<sub>3</sub>) and σ<sub>LL</sub> = E<sub>i</sub>ε<sub>LL</sub>, (Equation 6). E<sub>i</sub> is the Young's modulus of the i<sup>th</sup> tissue slice. I<sub>j</sub> measures the effect of all the internal lever arms of the resisting stresses in a given tissue and works out as I<sub>i</sub> = π R<sup>4</sup><sub>Out,i</sub>/4 − R<sup>4</sup><sub>in,i</sub> where R<sub>out,i</sub> and R<sub>in,i</sub> are the outer and inner radii of the j<sup>th</sup> annulus of tissue. These formula specify the consequences of loading (F), stem geometry and anatomy (L, I<sub>j</sub>), and material stiffness (E<sub>i</sub>) on

the stress and strain fields. Note that if strains increase linearly from the center to the periphery, stress may be distributed non-continuously. Another striking property of beam bending apparent in the previous equation is that changes in cross-sectional geometry have much more effect on stresses and strains than changes in material stiffness. For example, for a given load F, doubling the elastic stiffness of all the tissues halves the longitudinal strain without changing longitudinal stress. Doubling the stem radius however (keeping the same proportion of concentric tissues) reduces both strains and stresses at the stem periphery 8-fold. This is another example of mechanical amplification by lever arms. We will see that this has profound consequences on the mechanical stability and the mechanosensitivity of a given stem. Finally, stresses and strains are highly anisotropic, with their principal component lying longitudinally along the length of the stem.

**Shell model of the SAM under internal pressure load**

The SAM, a group of continuously growing and dividing cells,
is a dome-shaped structure (Figure 2) composed of two outer
layers, named L1 (the outer epidermis) and L2, and an inner
bulk of cells named L3. Future definite lateral organs, like leaves,
sepals, or petals, are initiated as primordia, bulges at the periphery
(see reviews by Kwiatkowska, 2008; Burian et al., 2013; Robinson
<figure><img src="image_7.png" /></figure>

**FIGURE 4** | The ISM beam model of pure bending of a stem. ISM model used to analyze stem bending experiments, using the theory of composite heterogeneous beams in a cantilever setting. (A) Unloaded beam. The beam is composed of a pile of (virtual) slices of infinitesimal thickness delimited by (virtual) successive cross-sections, along a central line. (B) Loaded beam. Under bending moment $\vec{M}(\zeta)$, the beam curves. Each cross-section rotates by a small angle $d\theta(\zeta)$, with $\zeta$ being the position along the stem and x, y the coordinates within the current cross-section of the stem. (C-E) Detailed side (C,D) and top (E) views of a bent slice in a homogeneous stem. (F-H) Detailed side (F,G) and top (H) views of a bent slice in a heterogeneous stem made of one stiff (dark gray) and two compliant (light gray) concentric annuli of tissues. (C,F) Strain distribution across the cross-section. Note that the cross-section remains flat during the bending and only rotates respective to the previous cross-section at the bottom of the slice by an angle $d\theta$, irrespective of the anatomy of the stem. The spatial rate of change in angle of the successive cross-sections is the stem curvature $C = \frac{d\theta}{ds}$. Accordingly, the stem is elongated on the convex side by $dl(x, y) > 0$ and shortened on the

concave side by $dl(x,y) < 0$, with no change on the central (neutral) line. The longitudinal strain $\epsilon_{LL} = \frac{dl}{ds}$ is thus maximal at the periphery on the sides of the slice that face downwards and away from the orientation of the bending force. The heterogeneous anatomy of the stem has no effect on the relative distribution of strain across the cross-section, which remains linear and is given by $\epsilon_{LL,y} = y \cdot (C - C_0)$. Straining allows for internal reaction forces, which density is measured by stresses, to build up balancing the effect of the external load. Therefore, the amount of change in stem curvature (and hence the global amount of straining) only depends on the amount of external bending moment and on the overall bending stiffness of the stem. (E,F) Stress distribution in the cross-section. For elastic constituents, the stress is equal to the strain multiplied by the Young's modulus: $\sigma_{LL,x,y,i} = \epsilon_{LL,x,y} E_{LL,i}$ where $E_{LL,i}$ is the longitudinal elastic modulus of material i. In homogeneous stems stresses parallel strains. However, on a stem with a heterogeneous anatomy (F) the stresses also depend on the local stiffness of the tissue and they de-correlated with strains across the cross-section (with maximal stresses possibly occurring inside the stem).

et al., 2013). Between the primordium and the apical dome, a saddle-shaped boundary forms which later becomes a sharp crease that separates the growing primordium from the SAM. The thin-walled fully turgid cells in the SAM are under considerable mechanical load from turgor pressure and cell-to-cell mechanical interactions known as “tissue tensions.”

The mechanical analysis of the meristematic dome has been performed by Hamant et al. (2008) and Traas and Hamant (2009), first giving rise to the Pressurized Vessel model (PVm). The meristematic zone is modeled as a vessel according to thin-shell theory. The vessel “wall” corresponds to the outer wall of the L1 layer, which is thicker than the other walls and likely to bear much of

the load due to the turgor pressure of inner cells. The modeled vessel wall is built of thin shell elements of infinitesimal dimensions $ds$ and $dr$ and of thickness $t$. Their material properties are considered to be homogeneous across the SAM (for discussion see Baluska et al., 2003). Hamant et al. (2008) proposed that the elements should be linear elastic, but this is not necessary as the material could equally well be viscoelastic. The only restriction is that the material element should not show pure plastic properties as this would induce loss of rheological homogeneity during loading. These shell elements are smoothly assembled (i.e., essentially they are virtually “glued” together along their sides) into a typical SAM structure, modeled geometrically by
<figure><img src="image_8.png" /></figure>

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

**FIGURE 5** | The thin-walled pressurized vessel model of the shoot apical meristem. Integrative structural mechanics (ISM) model used to analyze the loading of the SAM by internal turgor pressure, using very thin shell theory (Hamant et al., 2008), reprinted with permission from AAAS. (A) The SAM modeled as a pressurized vessel. Each point has a coordinate in the orthoradial (r) and meridional (s) direction and P is pressure. (B) At the top of the apical dome, represented as a spherical dome, the stress is isotropic. If the flanks of the meristem are represented as a cylinder, the stress is greater in the circumferential (orthoradial) direction than along the meridian and strongly anisotropic stresses occurs on the flanks of the meristem. Maximal stress anisotropy occurs at the saddle-shaped boundary between the primordium and the central dome.

combining three simple adjoining structures (Figure 5). (i) The apical dome is represented as a spherical dome of radius R. (ii) The flanks of the meristem are represented by a cylinder of radius R. (iii) Where relevant an incipient primordium is represented by a smaller lateral dome. Each point of the vessel “wall” is characterized by its coordinates in orthoradial (r) and meridional (s) directions (Figure 5A). The difference between a beam slice and a shell element is that each shell element can display curvatures in two directions (i.e., $C_{rr}$ along an orthoradial line, $C_{ss}$ along a meridian) and may also display a twist ($C_{rs}$). Note that in the central or primordial domes, $C_{ss}$ and $C_{rr}$ have the same sign, the concave surface faces into the meristem and $C_{rs} = 0$), whereas in the boundary, $C_{ss}$ and $C_{rr}$ have opposite signs. Importantly, this geometry is assumed to be under static equilibrium, so that the model only aims at calculating the stresses required to achieve this static equilibrium in the specified geometrical configuration and strains are unknown. The load is considered to be a uniform and constant inner pressure P. The effect of this internal load in the model can be described verbally as follows. Each shell element builds up stresses in three directions, and its stress state is thus represented by a stress tensor ($\sigma_{ss}$, $\sigma_{rr}$, $\sigma_{rs}$), with $\sigma_{ss}$ and $\sigma_{rr}$ being tensions in the meridional and orthoradial directions, and $\sigma_{rs}$ is a shear stress within the vessel wall (a kind of internal friction stress). The values of each stress component can be fully estimated using the conditions of static equilibrium, and depends on the curvatures in each zone of the SAM. This was solved analytically at specific points of local symmetry.

In the central dome, equilibrium yields $\sigma_{ss} = \sigma_{rr} = \frac{PR}{2}$, $\sigma_{rs} = 0$ (Equation 7). The tensile stresses are isotropic as they have same value in the s and r directions. In the flanks of the meristem cylinder, $\sigma_{rr} = PR$, $\sigma_{ss} = \frac{PR}{2}$ (Equation 8), so the stress is highly

anisotropic, being twice as high in the circumferential direction as in the longitudinal direction.

In the saddle-shaped boundary between the apex and a primordium, one may assume that the orthoradial curvature is approximately the curvature of the dome $C_{rr} \approx \frac{1}{R}$ whereas the meridional curvature $C_{ss}$ is much higher (in absolute terms). The stresses matching static equilibrium are therefore: $\sigma_{ss} \approx \frac{P}{2C_{ss}}$, $\sigma_{rr} \approx \frac{P}{2C_{rr}} \approx \frac{PR}{2}$, $|\sigma_{rr}| > |\sigma_{ss}|$ (Equation 9). The outer wall of the SAM is under tension across the crease, but in compression along the crease. The amount of the stress depends on P and on one of the two curvatures $C_{rr}$ and $C_{ss}$, with higher curvature inducing lower stresses for a given P. The absolute amount of stress is much higher across the crease, and the stress distribution is highly anisotropic.

This model of the SAM as a thin-walled pressurized vessel was very instructive. However, it did not provide detailed insights into how stress is distributed in the cell walls of a given SAM zone. It was therefore, complemented by a second “zoom-in” model at the scale of a small patch of tissue in the L1 and L2 layers (Figure 6). This model is a Finite Elements model (FEm) of two patches of the L1 + L2 tissues, one in the primordium-boundary zone and one at the top of the dome. Detailed specification of the geometry of the cell walls of this patch was achieved by experimental measurements. Elements were meshed piecewise to form plates of finite size. Just as for the PVm, the material was assumed to be homogeneous, but also linear elastic (with no plastic deformation or growth) and the load resulted from a uniform internal pressure putting the L1 layer into a curved configuration. No internal pressure within the cells of the L1 and L2 layers was considered as the effect of uniform pressure among the neighboring cells cancels out within a cell layer. The boundaries of the patch were given the proper saddle shape of the primordial boundary or the hemispherical shape of the dome. The mechanical equilibrium state was then computed numerically. This was done both for the intact patch, and for a patch in which one or two cells were ablated (i.e., their outer walls were deleted from the model). Making one hole in the patch redirects the orientation of the main stresses to surround the hole, and increases the magnitude of these stresses. When holes are made in two adjacent elements, highly anisotropic stresses are induced between the two holes, even if they are positioned at the tip of the dome where stresses are normally isotropic. Note that as the elastic rheology of the cell wall was specified, the FEm could be used to estimate not only the stress distribution within the cell walls, but also the elastic strains of the walls (this was not possible using the PVm).

## MECHANOSENSING AND MECHANOTRANSDUCTION

Now that we have tracked the distribution of stresses and strains within the two types of organs and the two types of loads, we can study how the cells sense their local mechanical state. Models can be helpful tools at this stage too to define quantitative behaviors and to deduce which variable is sensed.

## STRAIN-SENSING OR STRESS-SENSING? DOES IT MATTER?

Mechanobiologists have paid relatively little attention to the issue of whether plant cells sense stress or strain, implicitly assuming that mechanical “stress” is the variable of interest perhaps due
<figure><img src="image_10.png" /><figcaption>**FIGURE 6** | Finite elements model (FEm) of a patch of the L1 + L2 layer of the SAM. ISM model for the numerical mechanical analysis of a small patch of the SAM with full cellular resolution. The example here displays the model (and the numerical simulation of its stress-field output) of a patch at the boundary between the primordium and the central dome, with the ablation of one L1 cell from Hamant et al. (2008), reprinted with permission from AAAS. (A) General view of the FEm of the patch in the primordium boundary zone from above indicating the simulated pattern of principal stress directions (red lines) on the outer surface of meristem tissue. Colors indicate relative values of stress (blue, low; green, medium; red, high). (B) Side view of the outermost cell layers L1 and L2. (C) Detail of the stress pattern around one hole due to cell ablation.</figcaption></figure>

to semantic confusion with physiological “stress” (Moulia et al., 2011). However, it is important to remember that stresses and strains do not parallel in heterogeneous constitutive materials (e.g., stems) or in materials behaving in the plastic range (e.g., during growth). Thus, a strain-sensing mechanism would not give the same output as a stress-sensing mechanism. Recently a “stress-sensing vs. strain-sensing controversy” has been stirred up [see replies to Hamant et al. (2008) by Schopfer and Meyerowitz in science e-letters, and (Moulia et al., 2011) and (Hamant, 2013)]. Addressing this issue requires a further step in the modeling.

## FROM CELLULAR MECHANISMS TO QUANTITATIVE LOCAL MECHANOSENSING

### Local mechanosensing of external loads: the “strain-sensing model”

Among the mechanisms involved in mechanosensing, mechanosensitive ionic channels, often known as stretch-activated channels (SAC) have been the subject of detailed quantitative studies using the patch-voltage-pressure-clamp technique on protoplasts, cells enzymatically stripped of their walls (e.g., Ding and Pickard, 1993; Haswell et al., 2008). Altering the turgor pressure induces strains and tensional stresses in the plasma membrane and in the channel. The ionic current passing through the channels can be monitored after clamping the voltage, thus quantifying their mechanosensitive responses. The general shape of these response curves is sigmoidal, and can easily be linearized in the range of small strains (Figure 7A). Based on these results, we assumed that the local mechanosensitive function of a tissue element can be approximated through a linear function over a threshold (Coutand and Moulia, 2000; Moulia et al., 2011):

$$dS_i = k_s \cdot (\varepsilon - \varepsilon_0) \cdot dV \quad \text{if } \varepsilon > \varepsilon_0, \text{ else } dS_i = 0 \quad (10)$$

where $dS_i$ is the local signal in the cell (in Figure 7A, $dS_i = dI$, where I is the ionic current), $k_s$ is a mechanosensitivity factor ($k_s = 0$ defines an insensitive tissue, while higher $k_s$ values equate to greater sensitivity), $\varepsilon$ is the local mechanical strain of the tissue element, $\varepsilon_0$ is a possible strain threshold or minimal effective strain ($\varepsilon_0 \ge 0$) (see Moulia et al., 2006 for a review), and $dV$ is the volume of the tissue element.

Equation (10) assumes that only tensile strains are sensed ($\varepsilon > \varepsilon_0 \ge 0$), but it can be extended straightforwardly to the case where both tensile and compressive strains are sensed in proportion to their absolute value, as is observed in animal bone tissues (Schriefer et al., 2005).

Equation (10) was assessed experimentally in Populus tremula × alba (Pta) (Coutand et al., 2009) by measuring the expression of the primary mechanosensitive gene ZFP2. ZFP2 codes for a zinc finger protein that is transiently over-expressed as early as 5 min after straining in the strained tissues, probably in a cell-autonomous manner (Leblanc-Fournier et al., 2008; Martin et al., 2009, 2010) The response of the cell mechanotransduction pathway—from the initial reaction in the cytoplasm to primary gene expression in the nucleus—could thus be assessed by measuring $Q_r$, the relative abundance of ZFP2 transcripts in small slices of the stem using quantitative real-time PCR (Coutand et al., 2009). The bending stresses and strains are highly heterogeneous across a stem element. An integrative model was thus necessary to express the prediction of Equation (10) at the level of a stem segment and to assess it experimentally. Combining Equation (10) with the strain field equation in bending (Equation 4), it was possible to derive

$$Q_{\text{organ}} = \frac{k_s \cdot k_{ds}}{C_0} \cdot \bar{\varepsilon} - \left( \frac{k_s \cdot k_{ds}}{C_0} \cdot \bar{\varepsilon}_0 - 1 \right) = k_r \cdot \bar{\varepsilon} - (k_r \cdot \bar{\varepsilon}_0 - 1)(11)$$

where $Q_r$ is the ratio between the abundances of Pta ZFP2 transcripts in the strained tissue elements and those in the unstrained control), $k_{ds}$ is the sensitivity of the pathway downstream of the primary sensory reaction, $C_0$ is the baseline transcript concentration in the unstrained control, $k_r = k_s k_{ds} / C_0$ is the apparent sensitivity of relative gene expression, and $\bar{\varepsilon}$ is the volume-averaged tensile strain (see Moulia et al., 2011 for details).

Our local mechanosensing model (Equation 11) thus predicts a linear increase in the relative expression of ZFP2 with an increase in the mean strain $\bar{\varepsilon}$, a prediction that can be tested experimentally. Indeed, the experimental relationship between measured $Q_r$ and volume-averaged strain $\bar{\varepsilon}$ was found to be linear (Figure 7B), with Equation (11) explaining 77% of the 1:500 variation in $Q_r$. This validates the hypothetical strain-sensing model stated in Equation (10) and gives the first in planta measurement of the mechanosensitivity of the mechanotransduction pathway. Under the conditions of this experiment, a 1% strain induces a transient 200-fold increment in transcription of Pta-ZFP2. It was surprising that the strain range in which this linear mechanosensing model holds true goes up to at least 5%, i.e., well beyond the range of elastic strains in cell walls. Cell internal components are likely to undergo a much larger range of elastic deformation than the cell wall alone, explaining the proportional sensing of strain even when wall stresses eventually plateau (see Sato et al., 2005).
<figure><img src="image_11.png" /><figcaption>**FIGURE 7** | Local mechanosensing of external loads. (A) Probability of mechanosensitive channel (MsC) opening and mean patch conductance as a function of patch depression (and hence membrane tension and MsC strain). Open and filled circles, two replicates. Dashed dotted line, linear fit. Modified from Ding and Pickard (1993), Copyright# 1993, The Plant Journal, John Wiley and Sons. (B) Relationship between the relative</figcaption></figure>

transcript abundance Qr of the primary mechanosensitive gene *Pta ZFP2* (measured by Q-RT-PCR) and predictions from the Strain-Sensing model through the volume-averaged strain in the bent stem segment *ε̄*, (i.e., Sum of the Strain-Sensing normalized to the volume of the bent tissue; Coutand et al., 2009, Journal of Experimental Botany, by permission of the Society for Experimental Biology).

**Local mechanosensing of internal loads: cellular stress-feedback model**

The mechanisms underlying the responses to internal loads have been investigated much less than those involved in sensing external loads. Mechanical signals control (i) the amount and distribution of the active PIN1 auxin transporters, possibly though Ca²⁺ influx acting on PINOID proteins via TOUCH3 proteins (Heisler and Lam, 2010; Nakayama et al., 2012), and (ii) the alignments of cortical microtubules (CMT) and the orientation of cell division planes. The calculated stress pattern in the SAM outer L1 layer, and the CMT distribution determined experimentally (Hamant et al., 2008) were very similar (Figure 8A).

However, this observation is only correlative. A step forward was made when it was confirmed that the distribution of microtubules changed to match the redistribution of the wall stresses as predicted by the local FEm when the meristematic dome was compressed or when two holes were made in the L1 layer, (Figures 8B,C). However, this still did not provide a mechanistic link. A putative sensing mechanism may involve wall-associated protein complexes linking the cell wall to CMT that would then be directly subject to cell wall stresses, but there is no direct experimental evidence for this at the subcellular level (see Baluska et al., 2003; Landrein and Hamant, 2013 for discussion).

In order to assess this mechanism of microtubule alignment by wall stresses with respect to the data on SAM dynamics, a new model was needed. Microtubule reorientation takes 4–12 h, long enough for growth to occur. Thus, a model was required that included cell geometry, growth, mechanosensing of load distribution, and microtubule orientation. This model was called the Two-Dimensional Cellular Stress Feedback model (2D SFm, Figure 9; Supplemental Material SM2).

**Mechanical structure.** The material element of this model is a piece of cell wall, assumed to display linear elastic properties, i.e., its stresses are proportional to its strains. The coefficient of proportionality is the stiffness of the wall material $E_w$, its Young's

<figure><img src="image_12.png" /><figcaption>**FIGURE 8** | Mechanosensing of internal loads in the SAM and microtubule re-orientation from Hamant et al. (2008), reprinted with permission from AAAS. (A) Schematic representation of stress directions and microtubule orientations in the different parts of an SAM bearing a cylindrical primordium. (B) Principal stress pattern at the outer surface of the meristem simulated in an FEm of a patch of SAM at the top of the dome with a two-cell ablation. The stress pattern is circumferential to each of the ablated regions and stress alignment is enhanced in the cell between the two ablated cells. (C) Cortical microtubule distribution in the L1 layer in the central zone after a two-cell ablation as visualized by the expression of a construct fusing the Green Fluorescent Protein and the Microtubule Binding Domain (GFP-MBD), Scalebar, 5 µm.</figcaption></figure>

modulus. Only one-dimensional (1D) stretching of the wall is considered. $E_w$ is under biological control, modeled as depending on the auxin concentrations in adjoining cells on both sides of the wall (see Supplemental data), and on microtubule orientation in the same cells ($\theta_{c1}$ and $\theta_{c2}$), according to:

$$E_w = E_{min} + E_{max} \left( \frac{\cos^2(\theta_1) + \cos^2(\theta_2)}{2} \right) \quad (12)$$
$E_{min}$ is the elastic stiffness of the isotropic cell wall matrix, and $E_{max}$ ($\frac{\cos^2(\theta_1)+\cos^2(\theta_2)}{2}$) is a stiffening term related to the directional angle $\Theta$ of CMTs (and hence microfibrils) relative to the wall direction on both sides of the cell wall. This $\cos^2(\Theta)$ angular dependency simply specifies mathematically the idea that parallel and antiparallel orientations both lead to the maximal longitudinal stiffening, whereas perpendicular orientation leads to no stiffening.

Interphasic expansion growth of the walls of meristematic cells is modeled by increasing the resting length of walls $l_{w0}$ at an absolute rate that is proportional to the amount of elastic strain of the wall above a yield threshold. Cell wall synthesis is assumed to follow wall extension (and cell division). The wall elements are assembled into a two-dimensional tissue model representing the cells of the L1 layer as hexagonal boxes, only considering the lateral walls of L1 cells. Each wall is considered to act as a 1D spring carrying a force $F_w$. The cell corners are assumed to behave like ball-joints in that there is no stiffness when the angle is changed. Finally, cell division is assumed to occur when cells reach a threshold size, and the new wall runs through the barycenter of the original cell and parallel to the direction of microtubules.

**Load modeling mechanical equilibrium.** Just as in the PVm, the load comes from the turgor pressure of the inner cells and is assumed to be homogeneous. This puts the L1 layer under 2D stress and the curved configuration allows local curvature to balance the inner force and the tensile reaction within the L1 wall.

**Mechanosensing and feed-back mechanisms.** The mechanosensitive reorientation of CMTs is modeled assuming that $\Theta_c$, the CMT direction for a cell, is sensitive to wall stresses. The model does not address the stress distribution within the wall-associated proteins and the cytoskeleton though. It is only assumed that the mean CMT orientation somehow follows the orientation of the net force resulting from the tensions in all the side walls:

$$ \theta_c = \theta \underset{\sum F_w}{\rightarrow} \qquad (13) $$

This alignment response is not instantaneous but occurs at a constant rate. This is the only mechanosensitive step in the model.

The feedback mechanism then comes from the following assumptions. The current mean CMT orientation is supposed to alter the elastic stiffness of each cell wall by modifying the direction of the cellulose microfibrils in the wall and hence the anisotropy of cell wall stiffness. This modifies the amount of growth in the different cell walls and the direction in which the new cell wall is laid down when a cell divides.

The inputs of the model are the turgor pressure, the shape of the L1 cells, and possibly the distribution of auxin. It predicts the elastic stresses and strains in the walls, the expansion growth, the mean orientation of CMT, and the orientation of phragmoplasts. These outputs depend on seven parameters ($E_{min}$, $E_{max}$, plastic growth extensibility, yield strain threshold for growth, the two parameters of auxin sensitivity, and rate of mechanosensitive reorientation of CMT, see Supplemental material for more details).

An important feature of this model is that different rates of expansion originate from the applied forces (pressure and

wall-wall interactions), and the dynamics of the cell wall elastic stiffness. Hence, there is no explicit relationship between maximal stress and maximal growth direction in the model. Depending on the stress patterns, the model can predict maximal growth along the maximal stress direction and perpendicular to it.

The validity of this model could not be assessed experimentally at the local level, but it was included in a model of the mechanosensitive behavior of the whole SAM (Section Whole-organ integration and experimental assessment).

# WHOLE-ORGAN INTEGRATION AND EXPERIMENTAL ASSESSMENT

## THE SUM OF STRAIN-SENSING MODEL, $S^3M$

### Quantifying global thigmomorphogenetic responses

To properly lay out the problem of integrated mechanosensing at this point, we need to consider the global growth responses of the plant to external bending loads in more detail. This has been made possible by using a quantitatively-controlled bending device while continuously monitoring primary elongation or secondary thickening using linear voltage displacement transducers (Coutand et al., 2000). It was found that elastic bending at the base of the stem induced a thigmomorphogenetic response in the distal primary growth zone, implying that a long-range internal secondary signal $S_{i,1}$ traveled from the bent tissues to the responding primary tissues (Coutand et al., 2000; also Brenner et al., 2006). The propagation of this signal to the apex is much faster than the typical reaction time of growth responses, and there is no obvious damping over longer distances (Moulia et al., 2011). The nature of the carrier of this long-distance signal is currently being investigated. Given its velocity it could be either an electric signal in the phloem, or more likely a pressure pulse in the xylem (Lopez-Rodriguez et al., 2014; Tixier et al., 2014). In contrast, the secondary growth response seems local to the bent zone (Mattheck and Bethge, 1998; Coutand et al., 2009). For both primary and secondary growth responses, the initially growth stops for one to a few hours, then growth restarts and eventually the growth rate returns to that of unstimulated controls. For primary growth, the recovery time is highly dependent on the amount of bending strain applied, typically ranging from 100 to 1000 min. No compensatory growth is observed so at the end of the experiment bent plants are shorter than control plants (e.g., 2 mm shorter per bending stimulus in the experiment by Coutand et al., 2000). Secondary growth though shows clear and long-lasting growth stimulation after the initial inhibition, with growth rate increasing over 3 days then decreasing to the control rate over the next 3–4 days. The effect of this stimulation of secondary growth (+0.35 mm per bending stimulus) was approximately 30 times higher than the effect of the initial inhibition, resulting in an overall stimulation of radial growth. Unlike primary growth, the timing of the response was much less dependent on the amount of bending strain than on the peak (and total) increment in growth rate (Coutand et al., 2009, 2010).

### Integrating local mechanotransduction into plant mechanosensing

Why do stems of different shape and structure respond differently to the same external load? We aimed to assess whether the strain-sensing hypothesis, combined with structural and geometrical effects on load distribution across the stem structure, can explain
<figure><img src="image_13.png" /></figure>

**FIGURE 9 |** Schematic representation of the SAM Stress Feedback model (SAM SFm). The SAM SFm (Hamant et al., 2008) incorporates an ISM biomechanical model of the mechanical load-bearing structure of the SAM, and a mechanobiological model of the responses to the mechanical state of a cell in terms of (i) cell-wall stress sensing by CMTs and (ii) the consequences on the longitudinal elastic stiffness of the cell wall due to the direction of the laying down of cellulose microfibrils with respect to the longitudinal direction of the cell wall. The elemental brick of the biomechanical model is a piece of cell wall Ⓐ (called the cell-wall element) which displays two rheological behaviors: (i) elastic straining and (ii) expansion growth changing the rest length $l_{w,0}$ of the element at a rate that is proportional to its elastic strain over a certain threshold. Cell wall growth is therefore analogous to visco-plastic creep. The transverse height (d) and thickness (t) of the wall element are assumed to be constant. Two levels of structure can then be assembled. At the first level, Ⓑ the side walls of a single hexagonal cell are assembled. The model can be run at this level, giving rise to the cell-level formulation of the SFm. Otherwise the cells can be assembled to form a surface mesh with typical SAM geometry，则 in the SAM-level formulation of the SFm. The load is the turgor pressure of inner tissues considered to be fully borne by the L1 cell(s) Ⓒ. This ISM biomechanical module outputs the (elastic) wall-stress field $\sigma_w^e$ at every position X on the different cell-wall elements composing the mechanical structure, at a given time, as well as the changes in rest-lengths of all the cell-wall elements due to expansion growth. This updates the geometry of the cellular structure for the next step and the outputs are transmitted to the

mechanobiological module. In the mechanobiological module of the SAM SFm, the mechanosensitive step occurs at the level of the cell Ⓒ, as it is an intrinsic cellular process. The central hypothesis of the module is that CMTs are re-aligned to the current direction of the direction of the resulting stress $\theta_c$, but this occurs at a constant pace with only some of the overall CMT population ($n_{cnew}$) being reoriented during each step (the rate is assumed to be independent of the mechanical state). The current mean orientations of CMTs in two cells sharing a given cell-wall ($\theta_1(t), \theta_2(t)$) determines the longitudinal elastic stiffness of the side cell wall $E_w$ presumably through the orientation of the laying down of the new cellulose microfibrils with respect to the existing wall, changing the anisotropy of the cell wall elastic rigidity and hence its longitudinal stiffness (here assumed to be instantaneous) Ⓓ (note that $(\theta_1(t), \theta_2(t))$ may differ from the targeted orientations $(\theta_{c1}, \theta_{c2})$ specified by the stress-feedback as CMT reorientation takes time). The latter process occurs at the level of each cell-wall element so different walls of the same cell differ in the amount of elastic stiffening they undergo. The new value of wall elastic stiffness $E_w$ is the output of the mechanosensitive module and is transferred to the mechanical module, in which the elastic stiffness is updated, immediately changing the constitutive law of the cell walls, and thus giving rise to a new mechanical equilibrium at the next step. This change in wall elastic rigidity is the way the current wall-stress state feeds back on the growth of the meristematic cells. Note that cell division may occur (not shown). A phragmoplast (new cell wall) is laid down parallel to the current CMT direction $\theta_1$ whenever the size of the stem passes a certain size threshold thereby changing the structure of the L1 cell wall “mesh.”

the variability in plant responses (Coutand and Moulia, 2000; Coutand et al., 2000, 2009). For this, we set out to build a *minimal model* of mechanosensing integration, from the level of the strained tissue element up to the thigmomorphogenetic growth responses in the entire stem. This model has been called the Sum of Strain-Sensing model ($S^3$m). Its development involved several steps (Coutand and Moulia, 2000; Coutand et al., 2009) and it was only completely assembled a few years ago (Moulia

et al., 2011). This model is designed to chart the effects on the global thigmomorphogenetic responses of both the mechanical and the mechanoperceptive structures of the organ. The model was made with the purpose of competitively assessing two possible candidate mechanisms for the mechanosensing of external loads, strain-sensing vs. stress-sensing. Building the model was analogous with the process of integrative modeling in structural mechanics as illustrated in **Figure 10** (presented in detail in the
<figure><img src="image_14.png" /><figcaption>FIGURE 10 | Schematic representation of thigmomorphogenetic model including the ISM beam model and the S³m. The ISM model of the mechanical load-bearing structure (left) is the CBmS designed by Coutand and Moulia (2000) to analyze stem bending experiments (see Figure 3). It is based on a validated composite-beam model of plant organ flexion (Moulia and Fournier, 1997). In its most simple configuration its inputs are the curvature field C(ζ) and the bending moment M(ζ) along the stem (measured as in Moulia et al., 1994). Its parameters are: (i) length, L, and diameters along the stem, D(ζ); (ii) estimates of tissue stiffness (longitudinal Young's moduli Coutand and Moulia, 2000); and (iii) the anatomical cross-sectional images processed using the model by Moulia and Fournier (1997). The elementary unit is a piece of tissue assumed to behave in the linear elastic range. Like all models based on beam-theory, this model defines two integration levels: the cross-section (which can be heterogeneous) and the stem. From the curvature field, it computes the strain field, ε(x, y, ζ, t) (and the stress field</figcaption></figure>

σ(x, y, ζ, t) if required) with ζ being the position along the stem and x, y the coordinates within the current cross-section of the stem, and t the time. The mechanosensitive model is the S³m model (Moulia et al., 2011). Its inputs are the strain fields ε(x, y, ζ, t) in each stem, and the stem geometry factors L and D(ζ), (all these data are received from the ISM beam model). S³m then generates the local sensing in a tissue element, and, if needed, the predicted amount of transcription Qr of a primary mechanosensitive gene. Then the integrated secondary signals Sᵢ,₁ and Sᵢ,₂ reaching the primary and secondary meristems, respectively, are computed. These signals are inputs of a module of thigmomorphogenetic growth responses (Coutand and Moulia, 2000; Coutand et al., 2009) outputting logarithmic dose-response modulations of primary and secondary growth. In a fully-coupled dynamic model of thigmomorphogenesis, the outputs of the thigmomorphogenetic growth response module can be used to update the size and geometry of the stem at the next step, so time integration can be simulated.

supplemental data), but extended it to include purely biological sensory responses.

The starting point was the local strain-sensing model (Equation 10) which states that the secondary signal output of each cell, dS₀, is proportional to the mechanotransduced signal in the mechanostimulated cell, and hence to dSᵢ (hypothesis H1).

The long-distance signal propagation was very fast compared with the growth response and was not damped down, so it could be neglected. The simplest model for the integration of the mechanical sensing is that the output signals, dS₀, of all the mechanosensitive cells simply sum up into a global secondary internal signal Sᵢ (hypothesis H2). In short, for a given strain amplitude, the more cells that are strained, the higher the Sᵢ is.

Subapical primary growth responds to distant sensing throughout the stem volume Vs. The internal signal propagated axially along the whole stem and controlling the response of primary growth Sᵢ,₁ can then be written as the sum total of all the

local signal outputs from the strained cells across the stem volume Vs (Coutand and Moulia, 2000):

$$S_{i,1}(\epsilon) = \iiint_V k_{o(\zeta,y,z)} \cdot (\epsilon_{(\zeta,y,z)} - \epsilon_0) \cdot dV \quad (14)$$

where ζ is the distance from the apex and (y, z) describes the position of the tissue elements across the cross-section of the stem and the triple sign means the sum along the three dimensions of the stem volume.

By analogy the thigmomorphogenetical signal controlling secondary growth Sᵢ,₂ can be computed as the sum of the elementary signals dS₀ on a one-cell thick cross-section (see Figure 10).

As can be seen in Equation (14) the mechanosensitive structure of the plant (at a particular time) is given by the spatial distribution of mechanosensitivity k₀(ζ,y,z) and threshold ε₀(ζ,y,z). However, comparative tests of the S³m have shown that the most
determinant factor was the geometry of the stem. For simplicity, more recent studies took mechanosensitivity to be homogeneous over all tissues (e.g., Coutand et al., 2009). If $k_o(\zeta,y,z)$ and $\varepsilon_o(\zeta,y,z)$ are constant, then they can be factorized in the spatial integrals, so that the model for the control of primary growth becomes

$$S_{i,1(\zeta)} = k_o \left( \iint\int_{V_s} \varepsilon_{(\zeta,y,z)} \cdot dV \right) - k_o \varepsilon_o \cdot V_s = k_o S_{1\text{strains}} - \Sigma_0 \quad (15)$$

where $S_{1\text{strains}}$ is the integrated stimulus summing all the strains of the cells and $\Sigma_0$ is the integrated minimal effective strain. The same expression can be used for secondary growth.

This model thus predicts that the integrated signals reaching the two meristems are linearly dependent on integrals of the strain field over the domains of mechanosensitive integration for primary and secondary growth ($S_{1\text{strains}}$ and $S_{2\text{strains}}$, respectively).

This prediction was tested quantitatively against the corresponding growth responses described earlier. Compared with what we studied for the local mechanosensitive gene expression (Equation 11 and Figure 7B), the prediction of $S^3$ m is no longer that the growth response should be linear with $S_{i,\text{strains}}$ but that the growth responds in a dose-dependent manner to the sum of strains $S_{i,\text{strains}}$. Another way of thinking about this is that collecting more cells in the strained tissues for the analysis necessarily entails adding RNA to the sample (linearity), but the biological thigmomorphogenetic response of growing tissues to the supposed integrated signal $S_i$ may not be additive.

As shown in Figure 11A, a tight logarithmic relation was found between the primary growth response (recovery time $\tau_r$) and $S_{1\text{strains}}$ which explains 75% of the 1:10 variation in the response (Coutand and Moulia, 2000).

$$\tau_{\text{recovery}} = a_1 \cdot \ln \left( \frac{S_{1\text{strains}}}{S_{0\text{1strains}}} \right) \quad \text{for } S_{1\text{strains}} > S_{0\text{1strains}} \quad (16)$$

where $a_1$ is the global thigmomorphogenetic sensitivity of a plant (including both the initial sensitivity of the mechanoperceptive structure of the plant $k_o(\zeta,y,z)$ and the responsiveness of the meristem to a long distance signal $S_{i,1}$) and $S_{0\text{1strains}}$ is the threshold over which meristematic cells perceive the global systemic signal reaching the growth zone. (Note that $S_{0\text{1strains}}$ is not the same as the integrated minimal effective strain threshold $\Sigma_0$ presented earlier).

Similarly, a relationship was found between $S_{2\text{strains}}$ and the radial growth response, which again explains 75% of the 1:5 variation generated by varying stem bending with different stem sizes (Coutand et al., 2009). An initial experiment on poplar suggested that a linear relationship between the radial growth response and $S_{2\text{strains}}$ was statistically slightly more significant than a logarithmic relationship. However, analysis of a set of dicot tree species (Coutand et al., 2010) showed that the logarithmic relationship is more generic (also see Telewski, 2006).

It should be noted that using $S^3$ m the global thigmomorphogenetic sensitivity of a plant can be described quantitatively using just two parameters for the primary growth response ($a_1, S_{0\text{1strains}}$) (Coutand and Moulia, 2000) and two for the secondary response

($a_2, S_{0\text{2strains}}$) (Coutand et al., 2010). Varying the load and/or plant size and anatomy affects the $S_{1\text{strains}}$ value along the x-axis in Figure 9, and thus the value of the response, but the relationship expressed in Equation (16) (and the corresponding log response curve) are invariant. This relationship and the $a_1$ parameter in Equation (16) are thus independent of both load intensity and plant size/structure.

Equation (15) involves an explicit integration of the effect of the mechanical and perceptive structures of the plant through the $S^3$ m model, a model that has been validated experimentally. This is not the same as a purely correlative dose-response curves with an “arbitrarily chosen” measure of the stimulus (e.g., force Jaffe, 1980a).

Finally, and very interestingly, we modified the local sensing equation so that local stress was the sensed variable instead of local strain. In this Sum of Stress-Sensing model, the 1:10 variation in the response was no longer explained, clearly disproving the stress-mechanosensing hypothesis for the control of growth by external mechanical loads (Figure 11B).

## INTEGRATIVE STRESS-SENSING MODEL IN SAM SUBJECT TO INTERNAL LOADS

The SFm of a cell has been assembled into a cellular network encompassing a realistic SAM shape to simulate the entire dynamics of meristematic growth and morphogenesis, including primordial bulging, phyllotactic patterning, and distribution of CMT orientation. This integrative SAM Stress Feedback model (SAM SFm) was derived from an existing model relating auxin transport to phyllotactic dynamics (Jonsson et al., 2006). To account for a realistic distribution of growth rate across the meristem an ad-hoc dependency of the growth extensibility parameter $k_g$ on distance to the center of the SAM was included (thereby introducing an 8th parameter, while artificially constraining the possible dynamics of the model).

The overall dynamics of the model was qualitatively satisfactory. Moreover, as shown in Figure 12 the predicted CMT orientation was found to match experimental observation at the primordial boundaries, at the side of a growing primordium and at the stem-side base of the SAM. The robustness of the stabilization of CMT orientation during the formation of the primordial boundary and crease could be also assessed. Last but not least, an alternative Growth-Strain Feedback model (GSFm) for cellular mechanosensing was implemented and simulated within the same framework, and found to produce erroneous qualitative predictions (see replies to Hamant et al., 2008) by Schopfer and Meyerowitz in Science e-letters).

However, in these studies, the prediction by the SAM SFm could not be quantitatively assessed vs. the observed CMT orientation. Growth-induced strains (and strain rates) were not measured concurrently, so the validity of the competing GSFm could not be really assessed. In other words, the GSFm might generate the wrong SAM dynamics not because its sensing hypothesis is wrong but because the modeling of growth and growth-induced auto-stresses is problematic. The elaborate task of testing the models was undertaken by Burian et al. (2013) who combined confocal measurement of CMT orientation dynamics with the sequential-cast (replica) technique, and careful comparative
<figure><img src="image_15.png" /><figcaption>**FIGURE 11** | Experimental assessment of the S³m model.
Dose-response curve of the recovery time of the primary growth response after bending plotted against the candidate internal signal (S₁, strains) predicted by the S³m model (adapted from Coutand and Moulia, 2000 Journal of Experimental Botany, by permission of the Society for</figcaption></figure>

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

Experimental Biology). (A) A logarithmic relationship is obtained under the hypothesis that the mechanosensored variable is the strain and which explains 72% of the overall response. (B) No relationship is obtained under the hypothesis that the mechanosensored variable is the stress. —, log fit; ♦, experimental results.

<figure><img src="image_17.png" /><figcaption>**FIGURE 12** | Experimental assessment of the 2D Stress Feedback model (SFm) of the entire morphogenetic dynamics of the SAM from Hamant et al. (2008), reprinted with permission from AAAS. (A)</figcaption></figure>

marking cortical microtubules (green) and cell shape (red) at the surface of a meristem generating a young primordium (P). Cortical microtubule marking is obtained using the expression of a fusion protein involving the Green Fluorescent Protein and the Microtubule Binding Domain (GFP-MBD) under the control of the constitutive promoter 35S (*p35S::GFP-MBD*) Scale bar, 20 mm. (B) Microtubule orientation (red bars) in cells in the 2D SFm (extracted from confocal data). Note the alignment of the virtual microtubule orientations in the boundary zone and compare to (A). (C) Simulation of an auxin-induced primordium. The 2D SFm results in orthod亊al alignment of microtubules around the growing primordium. (D) Tip-growing simulation with the stress-feedback model generating a growing stem. Microtubules align mainly orthod亊ally in the stem, which has a regular shape.

mapping of CMT orientation statistics, local surface curvature and strain-rates. First Burian et al. (2013) clearly disproved the GSFm and its central hypothesis that CMTs are always oriented perpendicular to the maximal growth. However, the analysis of the relationship between wall stress pattern and CMT orientation proved more complex. A priori two models were available to test for inferred stresses: the simple PVm (Figure 5), and the more complex numerical SAM SFm (Figures 9, 12). However,

as seen before, using the numerical cellular SFm requires the estimation of eight parameters, as well as exact knowledge of the initial geometry of cells. This requires very complex, destructive and time-consuming measurements. Therefore, the authors relied on the much simpler PVm (Hamant et al., 2008) which only requires mapping of the curvatures of the L1 layer. This was achieved through the sequential replica method followed by three-dimensional (3D) reconstruction and differential geometry. Burian et al. (2013) found that the presumed geometry-derived stress distribution is not sufficient to predict CMT orientation throughout the SAM (i.e., other than at the primordium boundary), thereby rejecting the simple stress-feedback hypothesis. They argued that a better, qualitative match between estimated developmental changes in stress and CMT were found when mechanical auto-stresses derived from differential and heterogeneous growth were considered. However, a full assessment of this new hypothesis requires an extension of the numerical cellular SFm to include differences in pressure in the L1 layer, so that in-plane stressing between neighboring cells can be considered. The modeling of possible differential behavior of the inner and outer sides of the L1 layer may also be considered (as has been done in drosophila embryo models Supatto et al., 2005). This is likely to require substantial changes in the model such as a full coupling with water flows, and therefore many additional parameters entailing more experimental assessment of models.

## MAJOR INSIGHTS FROM THE TWO EXAMPLES OF MODELING

These two sets of work on the mechanosensitive control of growth and morphogenesis by mechanical cues have a lot in common. In both cases it was necessary to combine integrative models with guided experiments and to progress iteratively in a model↔experiment loop. The role of the models was to allow for the predictions of mechanistic hypotheses to be assessed experimentally both at the cellular and whole-organ scales. In both cases, ISM integrative biomechanical models were coupled with a mechanobiological model to study the load distribution, the mechanosensing, and the building of the overall growth and morphogenetic responses. This has allowed the clear definition of the three “superimposed” structures of the plant: the mechanical
structure bearing and distributing the load, the mechanosensitive structure, and the responsive morphogenetic structure. The SAM could be considered a borderline case for this theoretical framework because there is very little visible tissue differentiation, so the three structures are merged (even though invisible dynamic patterning of cell fate is underway). On the contrary, the growing dicot stem offers a much clearer spatial and functional distinction between the three structures. The mechanical structure mostly involves the stiff tissues (although parenchyma does act as a filler and stabilizer, Gibson, 2013), the mechanosensitive structure is mostly parenchyma (and phloem), and the morphogenetic structures are the primary and secondary meristems. Therefore, changes in the geometry of the stem, and on the balance between the three structures may change the overall responsiveness of the plant to a given mechanical load (even if the intrinsic sensitivities of the mechanosensing cells and meristematic cells remain unchanged).

In all cases, the major insight is how the geometry of the organ influences both the load distribution and the amount of mechanosensing. Insightful dose-response curves cannot be obtained when only the overall mechanical load is considered (e.g., the force applied on the organ). In some experiments apparently good correlative dose-response curves with the overall load can be produced (e.g., Jaffe, 1980b), but this is because the three structures of the plant displayed very little variability, while there was a wide range of loads. Global dose-response curves risk causing confusion as they do not consider the mechanical and mechanobiological structures and processes involved in the interaction between the organ and its load and may be misleading when analyzing genetic variability of natural populations or mutant or stage variability.

The second insight brought by these two approaches is that the sensing is distributed, so that the relevant variable is the distributed change in mechanical state, i.e., stress or strain. It is now time to revisit briefly the strain-sensing vs. stress-sensing controversy (a more complete discussion can be found in Moulia et al., 2011). The conclusions of the two sets of studies outlined in this review are contradictory. However, these findings can be reconciled by postulating that different sensing pathways may be involved, each having a different mechanical relationship with the cell wall. Mechanosensing of external loads is thought to involve mechanosensitive ionic channels that sit in the soft cellular membranes and are gated by membrane tensional stresses (Haswell et al., 2008). They thus cannot sense wall stresses. But the membrane is more compliant than the cell wall by orders of magnitude and it is attached to and pressed against the cell wall, so it is forced to follow the straining of the wall. As a consequence the intrinsic stretch-activated-channels are stretched according to cell wall strains. On the contrary, cytoskeleton elements link to the cell wall at direct adhesion domains (Baluska et al., 2003). These linkers are partially embedded in the cell wall, so any cell wall stresses are directly transmitted to them. Moreover, unlike the surrounding polysaccharide wall, they are likely to behave elastically. So the change in their configuration, the strain, that ultimately modulates their biological activity (Monshausen and Haswell, 2013) is directly linked with elastic stresses within the surrounding wall. As these proteins are minute inclusions within the cell wall, their

very local stress-strain pattern cannot be resolved from that of the surrounding cell wall, so the changes in their configuration are better predicted by the stresses they receive from the surrounding cell wall than by its global strain. This illustrates how important it is to precisely consider the cellular (and macromolecular) structures involved, and how strain and stress patterns propagate, not only through tissues and cells, but ultimately in the apoplasm and symplasm. This is a new domain in which modeling will prove a very valuable tool.

## BIOMECHANICAL AND MECHANOBIOLOGICAL MODELS FOR PLANT MECHANOSENSING: FROM THE PLANT IN ITS ENVIRONMENT TO GENES AND BACK

### BIOMECHANICAL MODELS AS A TOOL TO DESIGN AND ANALYZE EXPERIMENTS

We now review more systematically some of the many ways in which biomechanical and mechanobiological models have been used to gain insights into the mechanosensitivity of plants in their natural environments.

### Mechanical models as state observers estimating stress and strain distributions

Mechanical models are almost indispensable tools for tracking how the mechanical loads are distributed throughout the organ or plant, scaling down to the mechanosensitive cells themselves. Indeed the load acts through the whole organ or plant structure, with possible mechanical focusing and amplification at specific places (related to lever arms, holes and curvatures). The stress distribution in the SAM (Figure 5) is primarily prescribed by the geometry of the dome and its connection with primordia, the saddle-shape of the crease having a huge effect. If one stem has a diameter that is twice that of another stem, when they are both submitted to the same bending moment M (Figure 4), the strains will be 8-fold higher in the more slender stem. If the two stems have same diameter the Young's modulus of one is twice that of the other, the strains in the stem with the lower Young's modulus will be only twice as high. Without a proper understanding of these effects a mechanobiologist can easily misinterpret experimental data. The case of stresses is particularly compelling and indeed, as already emphasized, stress is not an observable or measurable quantity. Engineers and mathematicians create models that act as "state observers" for the very purpose of inferring stresses in a structure (Villaverde and Banga, 2014). Interestingly though, a model to estimate stress distribution can be rather simple, as long as the equilibrium configuration of the organ and its mechanical load can be measured with sufficient accuracy and the organ is not rheologically too heterogeneous. In this case, the stress field can be calculated by applying the law of mechanical equilibrium. This was done in the PVm by Hamant et al. (2008). However, as soon as tissue patterns and/or growth-induced stresses among cells [also known as residual stresses or auto-equilibrated stresses, see (Moulia and Fournier, 2009; Moulia, 2013) for definitions] are involved, the picture becomes far more complex (Burian et al., 2013). Models then also need to consider strains and growth. Both elastic and growth-induced strains are observable quantities that can be measured even with non-contact techniques as long as the organ is accessible (e.g., Silk
and Erickson, 1980; Barbacci et al., 2013, 2014). However, only models can give insights into the full distribution of strains (and strain rates) across the tissues and their dependency on the structure of the plant and the load. These models acting as strain-state observers (such as the CBmS) are more advanced than stress-state observers, as they handle the changes in configuration of the organ structure under the load as well as the reference length of its constitutive elements, their changes through growth (e.g., Barbacci et al., 2013), and other visco-plastic effects.

### A guide for experimental strategies

We have seen that stem geometry has a major influence on the strains resulting from bending One problem that biologists face is that plants are never uniform. When studying thigmomorphogenesis, the classic experimental approach has been to apply a range of forces to “standardized” sets of plants. However, we now know that force is not the perceived variable. Therefore, a more appropriate strategy is to apply a range of strains in a controlled manner. The natural heterogeneity of stem diameters is a natural source of continuous variation in applied strains. We can therefore, take advantage of the two sources of variation in the amount of strain, the load, and the size of the plant. This strategy was used to disentangle the “strain-sensing vs. force-sensing” issue (Coutand and Moulia, 2000; Coutand et al., 2000). It is important to use appropriate controls as the size variations may confer other biological properties that produce confounding effects.

This experimental strategy can even be extended to study plants growing in windy conditions outdoors. In natural conditions, wind cannot be controlled easily. If the aim is to study a range of strains due to wind, using standardized plants will make the experimental design completely dependent on wind velocity so it can take a long time to accumulate the data on the desired range of strains. However, if plants with different geometries are chosen, the natural variability in geometry will bring a range of vibration frequencies and a range of strains for any given wind velocity. To choose suitable plant geometries, a mechanical model is required (Rodriguez et al., 2008; Sellier et al., 2008). The model would also give insights into the strain field, so the positioning of detectors measuring growth responses can be optimized.

Finally the wind-plant interaction model can be coupled with S³m (instead of the static bending mechanical model CBmS in Figures 4, 10) to test thigmomorphogenetic hypotheses on plant development outdoors.

### A TOOLBOX FOR MECHANISTIC SYSTEMS BIOLOGY

While mechanical and biomechanical models are instrumental in mechanobiology, hypothesis-driven plant mechanobiology can also benefit from well-based models.

### A tool for handling scale changes in integrative biology

Throughout this review, we have emphasized that models are necessary tools whenever effects of changes in organizational levels and scales are involved. Considering that the cell is the structural unit of life, it is intuitive to ask whether the cell with its constitutive molecules should be the usual scale for modeling. A central insight from biomechanical and biological modeling is that there is no absolute need to burrow down to the macromolecular or even cellular scale, even when cellular mechanobiological

responses are the object of the study (e.g., Hamant et al., 2008). This is certainly true for mechanical models. For example, both the CBmS and the PVm do not specify cells as finite elements, but they can still be used to understand the distribution of stress and strain in the apoplast, mechanosensitive cytoskeletal remodeling, and even gene expression. The SFm does not really consider the biological cell as the structural units as the inner mechanical structure of the cell, the cortical cytoskeleton, the transverse actin stress fibers, the vacuole, the nucleus or the endoplasmic reticulum are not considered. The stress and strain distribution through the cell ultrastructure is unknown. The model is simply a honeycomb-like apoplastic structure. Models of the inner structure of the cell are only necessary when the localization of the intracellular mechanosensitive responses are being studied, namely the precise molecular mechanism of mechanoreception (e.g., mechanosensitive channel opening or anisotropic cytoskeleton disassembly). The internal mechanical structure of the cell protoplast, and its links with the cell wall, is beginning to be revealed. Some models of animal cell ultrastructure have been built (Nick, 2011; Asnacios and Hamant, 2012), but they rely on simplified representation of the cellular structures (e.g., microtubule are modeled as cables). Detailed 3D molecular mechanical models are only feasible when studying the molecular dynamics of isolated macromolecules or oligomolecular complexes for which the crystallographic structure is known (e.g., a mechanosensitive channel in a patch of lipid bilayer, Sotomayor and Shulten, 2014). Another viewpoint is that plasmodesmata form a cytoplasmic continuum in tissues for mRNA trafficking between cells, so the cell is not the intuitive unit for mechanotransduction in plants. Finally, quantitative functional genomics has shown that master transcription factors may be less concentrated than their respective promoter sites in a given cell. In this context, transcriptional regulation at the single cell level may have a highly stochastic component that only becomes deterministic at the scale of a tissue element containing many interacting cells (Elowitz et al., 2002; Gandrillon et al., 2012).

### A Tool for Genetic Dissection

Thigmomorphogenetic responses have been described in many plant species and it is interesting to compare responses between species. The responses may be different because (i) the level of applied strain is different, for example, due to the geometry of the different species, (ii) the species sensitivity to mechanical loading is different, or (iii) both. The sensitivity must be an intrinsic character of the species and it must be independent of plant size. From the S³m model we know

$$ \frac{dR}{dt_{\max}} = a_2 \cdot \ln \left( \frac{S_{2\text{strains}}}{S_{0_{2\text{strains}}}} \right) = a_2 \cdot \ln(S_{2\text{strains}}) - a_2 \cdot \ln(S_{0_{2\text{strains}}}) \quad \text{for } S_{2\text{strains}} > S_{0_{2\text{strains}}} \quad (17) $$

Therefore, by applying similar ranges of S₂₃ₐ₅₄ and measuring the daily increment in diameter of different species, we can plot the growth response against the logarithm of the applied S₂₃ₐ₅₄, and fit a regression line in which the slope (a₂) is the intrinsic mechanosensing sensitivity of the species and the intercept can be used to estimate the mechanical signal threshold (S₀₂₃ₐ₅₄).
The values of parameters a₂ and S₀₂ can be compared to see if there is any variability in mechanosensing between species. This conceptual framework has been used to study the mechanosensing variability between five sympatric tropical tree species (Coutand et al., 2010). No variability in mechanosensing sensitivity was found but differences in signal thresholds were found (Figure 13). Note that this conceptual framework could also be used to compare mechanosensing sensitivity between organs for example.

### Enabling model-assisted phenotyping

The previous genetic study can be seen as an example of model-assisted phenotyping in which the S³m combined with a fairly simple experiment was used to extract two intrinsic parameters describing the mechanosensitivity of a species. The same approach can be used to phenotype collections of mutants or varieties. Defining intrinsic characteristics and enabling model-assisted high-throughput phenotyping is becoming a major challenge of systems biology (Bastien et al., 2013).

### Identifying rewarding molecular studies

Another interesting heuristic property of integrative models is that they can help to identify (i) which modules are more influential than others in the global mechanobiological response, and (ii) which are insufficiently understood. This can help to guide molecular studies of gene regulation, for example, so that the most rewarding projects are designed. Genes targeted in this way are likely to trigger particularly significant phenotypes, and identifying key regulators significantly would improve our understanding of the whole system response. An example is our recent study of thigmomorphogenesis in wind acclimation of plants, which clearly pinpointed a current deadlock in

advancing our understanding of mechanosensitivity accommodation in response to successive bending (see Leblanc-Fournier et al., 2014).

### UNDERSTANDING PLANT RESPONSES TO MECHANICAL LOADS IN THEIR NATURAL ENVIRONMENT

Ideally the ultimate goal of systems mechanobiology would be to understand how plants acclimate in their natural environment. For example, for wind resistance, models are available to describe the mechanics of lodging, wind-throw and wind-break (see De Langre, 2008; Gardiner et al., 2008) for reviews), but existing growth models disregard thigmomorphogenesis so they do not deal with wind acclimation (Moulia et al., 2006). A major asset of the S³m is that it can be coupled with mechanical models to analyze the effects of the static and dynamic strains produced by wind-induced vibrations in plants (e.g., Gardiner et al., 2008; Rodriguez et al., 2008; Sellier et al., 2008). On the other hand, as S³m can also handle growth responses to wind, it can be coupled with structure-function growth models (see Moulia et al., 2006; Fourcaud et al., 2008 for general discussion). The overall wind resistance of plants due to genetic variation in intrinsic mechanosensitivity could be assessed in various (present, past or forecasted) climatic scenarios during *in silico* numerical experiments. There is obviously a lot of work to accomplish, but integrative mechanosensing models are surely a key breakthrough paving the way to a better understanding of the ecological and economic relevance of thigmomorphogenetic acclimation, for example, by exploring the consequences of global climate changes on stand growth and resistance to wind hazards.

### CONCLUSION: THE CHALLENGES OF SYSTEMS MECHANOBIOLOGY

Plants respond to internal and external mechanical loads at many scales. The mechanical growth response can theoretically be broken down into four processes.

* Bearing the load, how the mechanical load is distributed by a plant structure.

* Sensing the effect of the load distribution by mechanosensing of local mechanical states.

* Transducing the signal through mechanoreceptor pathways to alter the expression of a specific set of transcription factors.

* Responding by the global retuning of growth rate and direction.

* However, none of these processes are disconnected from one another. The challenge of the integrative biology of mechanosensing is to operate across scales and processes and understand outputs in terms of the overall syndrome of mechanosensitive growth responses and their adaptive relevance.

As this loop involves several organizational levels and scales plus a host of interactions, it cannot be handled without models, placing this challenge within the realm of systems biology (Tardieu, 2003; Moulia and Fournier, 2009; Traas and Moneger, 2010). The bottom line here is that systems biology modeling and

<figure><img src="image_18.png" /><figcaption>FIGURE 13 | S³m-assisted dissection of natural genetic diversity in mechanosensing. Relationship between the standardized diametric growth responses of five neotropical forest species and the log of the candidate internal signal (S₂, strain) predicted by the S³m model. S₂, strain = sum of strains over the cross-section during the bending treatment from Coutand et al. (2010), Annals of Botany, by permission of Oxford University Press.</figcaption></figure>
cross-comparisons against data produced through suitable exper-
imentation makes it possible to test hypotheses. If a combination
of hypotheses cannot be worked out without calculus, modeling
becomes an extension of the experimental method (Legay, 1997).
Therefore, the mechanical models need to be adjusted whenever
necessary, such as to cope with a different experimental set-up or
to natural conditions.

These mechanistic models have to be evaluated on their
outputs, on the performance of their component mechanistic
modules (e.g., Coutand et al., 2009), and their capacity to reflect
natural genetic variation (e.g., Sierra-De-Grado et al., 2008;
Almeras et al., 2009; Coutand et al., 2010). To be useful how-
ever, models need to be kept simple and amenable enough to
use. Multiplying the number of elements and degrees of freedom
makes models more difficult to analyze and assess experimen-
tally. This difficulty may be partly circumvented by implement-
ing a clear modular design and explicitly setting up organiza-
tional levels. Additionally, the development of non-destructive
bio-imaging and micromechanical techniques (see references in
Milani et al., 2013; Moulia, 2013) are ways to gather more infor-
mation to precisely assess the modules. But there are still some
intrinsic limits to model complexity if the model is to be used to
assess mechanobiological hypotheses and/or to capture very non-
linear behavior. Here mechanobiology diverges from standard
mechanical engineering in which the model is mostly an attempt
to assemble well-established physical laws into a given structure,
not a way to assess a biological hypothesis. From this perspective,
it is instructive that the SAM SFm with its eight parameters has yet
to be assessed experimentally and that experimental progress has
relied on the PVm. Hence, an effort to simplify models is called
for (as is common in physics through dimensional analysis; see
Rodriguez et al., 2008; Bastien et al., 2013, 2014) at the same time
addressing them to specific hypotheses that can be experimentally
falsified or upheld.

Overall, this strategy of combining heuristic models with
experiments has provided many insights into the mechanosensi-
tive control of growth. However, the complexity of the structural,
dynamic, and regulatory aspects requires intense interdisciplinary
work (Moulia, 2013). In particular, biologists need to become
familiar with the key concepts of mechanics to collaborate pro-
ductively with physicists and modelers in designing, criticizing,
and experimentally assessing biomechanical and mechanobi-
ological models. It is hoped that our review has convinced
readers of the usefulness and reward of this interdisciplinary
approach and made it more accessible. This interactive approach,
including using modeling as a tool to extend the experimen-
tal method, contrasts with and complements alternative set-ups,
for example, where experimentalists intensively collect data while
bioinformaticians set up models and data mining programs.

The chosen examples considered in this review deal only
with dicot growth, but the approach can be extended easily to
monocots (with obvious adaptions to omit secondary growth
responses), and even to most *Viridiplantae*. Again the examples
all involved the shoot system, ranging from the cell to the organ
level, but the methodology is in no way specific to these sys-
tems. Indeed there is a very active research making extensive
use of integrative modeling to discover how root growth and

morphogenesis is controlled by external and internal mechani-
cal cues (see for example Ditengou et al., 2008; Bengough et al.,
2011; Band et al., 2012; Jin et al., 2013). By the same token, works
at the scale of the intracellular structures would be extremely use-
ful. The way stresses and strains act inside the cell through the
cell wall-membrane-cytoskeleton-nucleus continuum is yet to be
fully grasped. The tensegrity model of the structural mechanics of
the cytoskeleton has started to reveal unexpectedly how mechan-
ical amplifications can take place through the cytoskeleton net-
work (Nick, 2011), but the way it interacts dynamically with wall
stressing and deformation remains to be elucidated (Hamant,
2013). Similarly, the responses of cell trafficking and plasma
membrane straining and turnover upon cell stretching are very
promising fields (Asnacios and Hamant, 2012; Nakayama et al.,
2012). A more complete elucidation of the mechanosensitive
gene networks controlling the cellular response is also a priority
including how the cell division-cell elongation complex is modu-
lated by mechanical cues. Finally this review has mostly focused
on spatial integration, but the challenge of time integration is as
important. When successive signals are applied to a living system,
relative or absolute refractory periods of the mechanosensors and
gene-regulated accommodation of the mechanosensitivity make
it complicated to study the relationship between the signals and
the system's responses (see Leblanc-Fournier et al., 2014). The
analysis of the time aspect will also require two-way interaction
between models and experiments.

The first successes of integrative mechanobiology of growth
control have thus opened up a large set of questions for inter-
disciplinary research, which hopefully will elucidate many more
aspects of the way plants have evolved to bear loads and remain
stable when responding to internal and external mechanical
challenges.

ACKNOWLEDGMENTS

This work was supported by the grant AN5-09-BLAN-0245-01 (ANR, France, project "Senzo").

SUPPLEMENTARY MATERIAL

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

REFERENCES

Almeras, T., Derycke, M., Jaouen, G., Beauchene, J., and Fournier, M. (2009). Functional diversity in gravitropic reaction among tropical seedlings in relation to ecological and developmental traits. *J. Exp. Bot.* 60, 4397–4410. doi: 10.1093/jxb/erp276

Asnacios, A., and Hamant, O. (2012). The mechanics behind cell polarity. *Trends Cell Biol.* 22, 584–591. doi: 10.1016/j.tcb.2012.08.005

Baluska, F., Samaj, J., Wojtaszek, P., Volkmann, D., and Menzel, D. (2003). Cytoskeleton-plasma membrane-cell wall continuum in plants. Emerging links revisited. *Plant Physiol.* 133, 482–491. doi: 10.1104/pp.103.027250

Band, L. R., Fozard, J. A., Godin, C., Jensen, O. E., Pridmore, T., Bennett, M. J., et al. (2012). Multiscale systems analysis of root growth and development: modeling beyond the network and cellular scales. *Plant Cell* 24, 3892–3906. doi: 10.1105/tpc.112.101550

Barbacci, A., Diener, J., Hemon, P., Adam, B., Dones, N., Reveret, L., et al. (2014). A robust videogrammetric method for the velocimetry of wind-induced motion in trees. *Agri. Forest Meteorol.* 184, 220–229. doi: 10.1016/j.agrformet.2013.10.003
Barbacci, A., Lahaye, M., and Magnenet, V. (2013). Another brick in the cell wall: biosynthesis dependent growth model. *PLoS ONE* 8:e74400. doi: 10.1371/journal.pone.0074400

Baskin, T. I., and Jensen, O. E. (2013). On the role of stress anisotropy in the growth of stems. *J. Exp. Bot.* 64, 4697–4707. doi: 10.1093/jxb/ert176

Bastien, R., Bohr, T., Moulia, B., and Douady, S. (2013). Unifying model of shoot gravitropism reveals proprioception as a central feature of posture control in plants. *Proc. Natl. Acad. Sci. U.S.A.* 110, 755–760. doi: 10.1073/pnas.1214301109

Bastien, R., Douady, S., and Moulia, B. (2014). A unifying modeling of plant shoot gravitropism with an explicit account of the effects of growth. *Front. Plant Sci.* 5:136. doi: 10.3389/Fpls.2014.00136

Bengough, A. G., Mckenzie, B. M., Hallett, P. D., and Valentine, T. A. (2011). Root elongation, water stress, and mechanical impedance: a review of limiting stresses and beneficial root tip traits. *J. Exp. Bot.* 62, 59–68. doi: 10.1093/jxb/erq350

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

Braam, J. (2005). In touch: plant responses to mechanical stimuli. *N. Phytol.* 165, 373–389. doi: 10.1111/j.1469-8137.2004.01263.x

Brenner, E. D., Stahlberg, R., Mancuso, S., Vivanco, J., Baluska, F., and Van Volkenburgh, E. (2006). Plant neurobiology: an integrated view of plant signaling. *Trends Plant Sci.* 11, 413–419. doi: 10.1016/j.tplants.2006.06.009

Burian, A., Ludynia, M., Uyttewaal, M., Traas, J., Boudaoud, A., Hamant, O., et al. (2013). A correlative microscopy approach relates microtubule behaviour, local organ geometry, and cell growth at the *Arabidopsis* shoot apical meristem. *J. Exp. Bot.* 64, 5753–5767. doi: 10.1093/jxb/ert352

Coutand, C. (2010). Mechanosensing and thigmomorphogenesis, a physiological and biomechanical point of view. *Plant Sci.* 179, 168–182. doi: 10.1016/j.plantsci.2010.05.001

Coutand, C., and Moulia, B. (2000). Biomechanical study of the effect of a controlled bending on tomato stem elongation: local strain sensing and spatial integration of the signal. *J. Exp. Bot.* 51, 1825–1842. doi: 10.1093/jexbot/51.352.1825

Coutand, C., Chevolot, M., Lacointe, A., Rowe, N., and Scotti, I. (2010). Mechanosensing of stem bending and its interspecific variability in five neotropical rainforest species. *Ann. Bot.* 105, 341–347. doi: 10.1093/aob/mcp286

Coutand, C., Julien, J. L., Moulia, B., Mauget, J. C., and Guitard, D. (2000). Biomechanical study of the effect of a controlled bending on tomato stem elongation: global mechanical analysis. *J. Exp. Bot.* 51, 1813–1824. doi: 10.1093/jexbot/51.352.1813

Coutand, C., Martin, L., Leblanc-Fournier, N., Decourteix, M., Julien, J. L., and Moulia, B. (2009). Strain mechanosensing quantitatively controls diameter growth and PtaZFP2 gene expression in poplar. *Plant Physiol.* 151, 223–232. doi: 10.1104/pp.109.138164

Coutand, C., Mathias, J. D., Jeronimidis, G., and Destrebecq, J. F. (2011). TWIG: a model to simulate the gravitropic response of a tree axis in the frame of elasticity and viscoelasticity, at intra-annual time scale. *J. Theor. Biol.* 273, 115–129. doi: 10.1016/j.jtbi.2010.12.027

Couturier, E., Brunel, N., Douady, S., and Nakayama, N. (2012). Abaxial growth and steric constraints guide leaf folding and shape in acer pseudoplatanus. *Am. J. Bot.* 99, 1289–1299. doi: 10.3732/ajb.1100325

De Langre, E. (2008). Effects of wind on plants. *Annu. Rev. Fluid Mech.* 40, 141–168. doi: 10.1146/annurev.fluid.40.111406.102135

Der Loughian, C., Tadrist, L., Allain, J.-M., Diener, J., Moulia, B., and De Langre, E. (2014). Measuring local and global vibration modes in model plants. *C. R. Mécanique* 342, 1–7. doi: 10.1016/j.crme.2013.10.010

Ding, J. P., and Pickard, B. G. (1993). Modulation of mechanosensitive calcium-selective cation channels by temperature. *Plant J.* 3, 713–720. doi: 10.1111/j.1365-313X.1993.00713.x

Ditengou, F. A., Tealea, W. D., Kochersperger, P., Flittner, K. A., Kneuper, I., Van Der Graaff, E., et al. (2008). Mechanical induction of lateral root initiation in *Arabidopsis thaliana*. *Proc. Natl. Acad. Sci. U.S.A.* 105, 18818–18823. doi: 10.1073/pnas.0807814105

Dumais, J. (2013). Modes of deformation of walled cells. *J. Exp. Bot.* 64, 4681–4695. doi: 10.1093/jxb/ert268

Elowitz, M. B., Levine, A. J., Siggia, E. D., and Swain, P. S. (2002). Stochastic gene expression in a single cell. *Science* 297, 1183–1186. doi: 10.1126/science.1070919

Fourcaud, T., Zhang, X., Stokes, A., Lambers, H., and Korner, C. (2008). Plant growth modelling and applications: the increasing importance of plant architecture in growth models. *Ann. Bot.* 101, 1053–1063. doi: 10.1093/aob/mcn050

Fournier, M., Baillères, H., and Chanson, B. (1994). Tree biomechanics: growth, cumulative prestresses, and reorientations. *Biomimetics* 2, 229–251.

Gandrillon, O., Kolesnik-Antoine, D., Kupiec, J. J., and Beslon, G. (2012). Special Issue: chance at the heart of the cell. *Prog. Biophys. Mol. Biol.* 110, 1–4. doi: 10.1016/j.pbiomolbio.2012.05.006

Gardiner, B., Byrne, K., Hale, S., Kamimura, K., Mitchell, S. J., Peltola, H., et al. (2008). A review of mechanistic modelling of wind damage risk to forests. *Forestry* 81, 447–463. doi: 10.1093/forestry/cpn022

Gibson, L. J. (2013). The hierarchical stucture and mechanics of plant materials. *J. R. Soc. Interface* 9, 2749–2766. doi: 10.1098/rsif.2012.0341

Gibson, L. J., Ashby, M. F., and Easterling, K. E. (1988). Structure and mechanics of the iris leaf. *J. Mater. Sci.* 23, 3041–3048. doi: 10.1007/BF00551271

Hamant, O. (2013). Widespread mechanosensing controls the structure behind the architecture in plants. *Curr. Opin. Plant Biol.* 16, 654–660. doi: 10.1016/j.pbi.2013.06.006

Hamant, O., Heisler, M. G., Jonsson, 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

Haswell, E. S., Peyronnet, R., Barbier-Brygoo, H., Meyerowitz, E. M., and Frachisse, J. M. (2008). Two MscS homologs provide mechanosensitive channel activities in the *Arabidopsis* root. *Curr. Biol.* 18, 730–734. doi: 10.1016/j.cub.2008.04.039

Heisler, L. K., and Lam, D. D. (2010). 5-HT2C receptor ligands in the treatment of obesity and type 2 diabetes. *Eur. Neuropsychopharmacol.* 20, S182–S182. doi: 10.1016/S0924-977X(10)70171-1

Jaffe, M. J. (1980a). Morphogenetic responses of plants to mechanical stimuli or stress. *Bioscience* 30, 239–243.

Jaffe, M. J. (1980b). Morphogenetic responses of plants to mechanical stimuli or stress. *Bioscience* 30, 239–243. doi: 10.2307/1307878

Jin, K., Shen, J. B., Ashton, R. W., Dodd, I. C., Parry, M. A. J., and Whalley, W. R. (2013). How do roots elongate in a structured soil? *J. Exp. Bot.* 64, 4761–4777. doi: 10.1093/jxb/Ert286

Jonsson, H., Heisler, M. G., Shapiro, B. E., Meyerowitz, E. M., and Mjolsness, E. (2006). An auxin-driven polarized transport model for phyllotaxis. *Proc. Natl. Acad. Sci. U.S.A.* 103, 1633–1638. doi: 10.1073/pnas.0509839103

Kwiatkowska, D. (2008). Flowering and apical meristem growth dynamics. *J. Exp. Bot.* 59, 187–201. doi: 10.1093/jxb/erm290

Landrein, B., and Hamant, O. (2013). How mechanical stress controls microtubule behavior and morphogenesis in plants: history, experiments and revisited theories. *Plant J.* 75, 324–338. doi: 10.1111/tpj.12188

Leblanc-Fournier, N., Coutand, C., Crouzet, J., Brunel, N., Lenne, C., Moulia, B., et al. (2008). Jr-ZFP2, encoding a Cys2/His2-type transcription factor, is involved in the early stages of the mechano-perception pathway and specifically expressed in mechanically stimulated tissues in woody plants. *Plant Cell Environ.* 31, 715–726. doi: 10.1111/j.1365-3040.2008.01785.x

Leblanc-Fournier, N., Martin, L., Lenne, C., and Decourteix, M. (2014). To respond or not to respond, the recurring question in plant mechanosensitivity. *Front. Plant Sci.* 5:401. doi: 10.3389/fpls.2014.00401

Legay, J. M. (1997). *L'expérience et Le Modèle, Un Dicours Sur La Méthode*. Paris: Quae Editions.

Lopez-Rodriguez, A. M., Venkataraman, G., and Holmgren, M. (2014). Shaker Kv Channel's sugar remotion in real-time. *Biophys. J.* 106, 537a. doi: 10.1016/j.bpj.2013.11.2996

Martin, L., Leblanc-Fournier, N., Azri, W., Lenne, C., Henry, C., Coutand, C., et al. (2009). Characterization and expression analysis under bending and other abiotic factors of PtaZFP2, a poplar gene encoding a Cys2/His2 zinc finger protein. *Tree Physiol.* 29, 125–136. doi: 10.1093/treephys/tpn011

Martin, L., Leblanc-Fournier, N., Julien, J. L., Moulia, B., and Coutand, C. (2010). Acclimation kinetics of physiological and molecular responses of plants to multiple mechanical loadings. *J. Exp. Bot.* 61, 2403–2412. doi: 10.1093/jxb/erq069

Mattheck, C., and Bethge, K. (1998). The mechanical survival strategy of trees. *Arboric.* J. 22, 369–386. doi: 10.1080/03071375.1998.9747222

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/ert169

Monshausen, G. B., and Haswell, E. S. (2013). A force of nature: molecular mechanisms of mechanoperception in plants. *J. Exp. Bot.* 64, 4663–4680. doi: 10.1093/jxb/ert204
Moulia, B. (2000). Leaves as shell structures: double curvature, auto-stresses, and minimal mechanical energy constraints on leaf rolling in grasses. *J. Plant Growth Regul.* 19, 19–30. doi: 10.1007/s003440000004

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

Moulia, B., Fournier, M., and Guitard, D. (1994). Mechanics and form of the maize leaf: *in vivo* qualification of flexural behavior. *J. Mater. Sci.* 29, 2359–2366. doi: 10.1007/Bf00363427

Moulia, B., and Fournier, M. (1997). Mechanics of the maize leaf: a composite beam model of the midrib. *J. Mater. Sci.* 32, 2771–2780. doi: 10.1023/A:1018604012754

Moulia, B., and Fournier, M. (2009). The power and control of gravitropic movements in plants: a biomechanical and systems biology view. *J. Exp. Bot.* 60, 461–486. doi: 10.1093/jxb/ern341

Moulia, B., Coutand, C., and Lenne, C. (2006). Posture control and skeletal mechanical acclimation in terrestrial plants: implications for mechanical modeling of plant architecture. *Am. J. Bot.* 93, 1477–1489. doi: 10.3732/ajb.93.10.1477

Moulia, B., Der Loughian, C., Bastien, R., Martin, L., Rodriguez, M., Gourcilleau, D., et al. (2011). “Integrative mechanobiology of growth and architectural development in changing mechanical environments,” in *Mechanical Integrative of Plant Cells and Plants* ed P. Wojtaszek (Berlin: Springer-Verlag), 269–302.

Nakayama, N., Smith, R. S., Mandel, T., Robinson, S., Kimura, S., Boudaoud, A., et al. (2012). Mechanical regulation of auxin-mediated growth. *Curr. Biol.* 22, 1468–1476. doi: 10.1016/j.cub.2012.06.050editors

Nick, P. (2011). “Mechanics of the cytoskeleton,” in *Mechanical Integration of Plant Cells and Plants*, ed P. Wojtaszek (Berlin; Heidelberg: Springer), 53–90.

Niklas, K. J., and Spatz, H. C. (2012). Mechanical properties of wood disproportionately increase with increasing density. *Am. J. Bot.* 99, 169–170. doi: 10.3732/ajb.1100567

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

Pivato, D., Dupont, S., and Brunet, Y. (2014). A simple tree swaying model for forest motion in windstorm conditions. *Trees* 28, 281–293. doi: 10.1007/s00468-013-0948-z

Pot, G., Toussaint, E., Coutand, C., Le Cam, J. B., and Saudreau, M. (2014). A model to simulate the gravitropic response and internal stresses in trees, considering the progressive maturation of wood. *Trees* 28, 1235–1248. doi: 10.1007/s00468-014-1033-y

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

Rodriguez, M., De Langre, E., and Moulia, B. (2008). A scaling law for the effects of architecture and allometry on tree vibration modes suggests a biological tuning to modal compartmentalization. *Am. J. Bot.* 95, 1523–1537. doi: 10.3732/ajb.0800161

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

Sato, K., Adachi, T., Matsuo, M., and Tomita, Y. (2005). Quantitative evaluation of threshold fiber strain that induces reorganization of cytoskeletal actin fiber structure in osteoblastic cells. *J. Biomech.* 38, 1895–1901. doi: 10.1016/j.jbiomech.2004.08.012

Schriefer, J., Warden, S., Saxon, L., Robling, A., and Turner, C. (2005). Cellular accommodation and the response of bone to mechanical loading. *J. Biomech.* 38, 1838–1845. doi: 10.1016/j.jbiomech.2004.08.017

Sellier, D., Brunet, Y., and Fourcaud, T. (2008). A numerical model of tree aerodynamic response to a turbulent airflow. *Forestry* 81, 279–297. doi: 10.1093/forestry/cpn024

Sierra-De-Grado, R., Pando, V., Martinez-Zurimendi, P., Penalvo, A., Bascones, E., and Moulia, B. (2008). Biomechanical differences in the stem straightening process among Pinus pinaster provenances. A new approach for early selection of stem straightness. *Tree Physiol.* 28, 835–846. doi: 10.1093/treephys/28.6.835

Silk, W. K., and Erickson, R. O. (1980). Local biosynthesis rates of cytoplasmic constituents in growing tissue. *J. Theor. Biol.* 83, 701–703. doi: 10.1016/0022-5193(80)90197-6

Sotomayor, M., and Shulten, K. (2014). Molecular dynamics study of gating in the mechanosensitive channel of small conductance MscS. *Biophys. J.* 87, 3050–3065. doi: 10.1529/biophysj.104.046045

Supatto, W., Debarre, D., Moulia, B., Brouzes, E., Martin, J. L., Farge, E., et al. (2005). *In vivo* modulation of morphogenetic movements in *Drosophila* embryos with femtosecond laser pulses. *Proc. Natl. Acad. Sci. U.S.A.* 102, 1047–1052. doi: 10.1073/pnas.0405316102

Tardieu, F. (2003). Virtual plants: modelling as a tool for the genomics of tolerance to water deficit. *Trends Plant Sci.* 8, 9–14. doi: 10.1016/S1360-1385(02)00008-0

Telewski, F. W. (2006). A unified hypothesis of mechanopception in plants. *Am. J. Bot.* 93, 1466–1476. doi: 10.3732/ajb.93.10.1466

Telewski, F. W., and Pruyn, M. L. (1998). Thigmomorphogenesis: a dose response to flexing in Ulmus americana seedlings. *Tree Physiol.* 18, 65–68. doi: 10.1093/treephys/18.1.65

Tixier, A., Badel, E., Franchel, J., Lakhal, W., Leblanc-Fournier, N., Moulia, B., et al. (2014). Growth and molecular responses to long-distance stimuli in poplars: bending vs flame wounding. *Physiol. Plant.* 150, 225–237. doi: 10.1111/ppl.12089

Traas, J., and Hamant, O. (2009). From genes to shape: understanding the control of morphogenesis at the shoot meristem in higher plants using systems biology. *C. R. Biol.* 332, 974–985. doi: 10.1016/j.crvl.2009.09.008

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

Villaverde, A. F., and Banga, J. R. (2014). Reverse engineering and identification in systems biology: strategies, perspectives and challenges. *J. R. Soc. Interface* 11:20130505. doi: 10.1098Rsif.2013.0505

**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.

Received: 14 August 2014; accepted: 20 January 2015; published online: 23 February 2015.

Citation: Moulia B, Coutand C and Julien J-L (2015) Mechanosensitive control of plant growth: bearing the load, sensing, transducing, and responding. *Front. Plant Sci.* 6:52. doi: 10.3389/fpls.2015.00052

This article was submitted to Plant Physiology, a section of the journal Frontiers in Plant Science.

Copyright © 2015 Moulia, Coutand and Julien. 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.