Introduction {#sec001}
============

Plants are constantly moving to reach the light and to maintain the architecture and posture of their aerial organs: stems, branches, leaves... These movements are active, generally powered by differential growth and are controlled by environmental signals \[[@pcbi.1004037.ref001]\]. Tropisms are the class of movements that are oriented by a vectorial environmental factor. Light is the main source of energy for plants and is a major cue for tropic movement. In phototropism, shoots grow in the direction of the light source. Under natural conditions on Earth, gravity is unavoidable and drives gravitropism, in which shoots usually grow against the direction of gravity. Proprioception, the ability of plants to perceive their own deformations, has recently been identified as a major factor in tropic movement; it stimulates the straightening of a curved organ, a response that can be thought of as autotropism \[[@pcbi.1004037.ref002]--[@pcbi.1004037.ref004]\]. Such classification of control mechanisms into three different tropic drivers is however merely a conceptual convention: it is likely that, in natural conditions, the three processes interact constantly. Plant shoots generally exhibit negative gravitropism, positive phototropism, and negative autotropism in the same organs at the same time.

The way phototropism and gravitropism interact to control the movement of plant shoots is acknowledged to be important in determining plant habit in natural conditions but is not well understood \[[@pcbi.1004037.ref005]\]. Advanced genetic studies in the model species *Arabidopsis thaliana*, which used mutants with impaired gravitropism or phototropism \[[@pcbi.1004037.ref006], [@pcbi.1004037.ref007]\], have shown that the two processes interact in a complex manner to control movement in wild-type plants \[[@pcbi.1004037.ref008]\]. However, part of this complexity lies in the variety of phototropic responses observed in response to different qualities and intensities of the light stimulus. For example, the effects of low-fluence pulses differ from those of continuous light \[[@pcbi.1004037.ref007]\]; but only the latter are directly relevant to plant growth under natural conditions, and will therefore be considered here.

The focus of the current study is the regulation and the control of organ movements during the interaction of tropisms. To do so, it is proposed to extend the recent dynamical model for gravitropism called the *AC* model \[[@pcbi.1004037.ref002]\] to phototropism. In this model, active tropic bending is controlled by the additive (but opposing) effects of graviception and proprioception, expressed by $$\frac{\partial C\left( s,t \right)}{\partial t} = - \beta A\left( s,t \right) - \gamma C\left( s,t \right)\text{for} s > L - L_{gz}\text{and} 0\text{otherwise}$$ where *s* is the position along the organ, *L* ~*gz*~ is the length of the growth zone, *L* is the length of the entire organ, *t* is time, *A*(*s*, *t*) is the local angle of the organ to the vertical, *C*(*s*, *t*) is the local curvature (i.e., the spatial rate of change of *A* along *s*) and the parameters *β* and *γ* are, respectively, the gravi- and proprio-ceptive sensitivities. The *AC* model was shown to explain the complex kinematics of gravitropism in eleven species covering a broad taxonomical range of angiosperms, major growth habits and organ types. \[[@pcbi.1004037.ref003], [@pcbi.1004037.ref004]\]

Although graviception and proprioception act additively, their control over the dynamics of tropic movement and the steady-state final shape actually depends only on the ratio between gravisentitivity and propriosentitivity, scaled to the size of the growth zone. In the model, this is formalized through the definition of the dimensionless graviproprioceptive bending number *B* = *β L* ~*gz*~/*γ. B* fully defines both the time to reach the steady state and the final shape at steady state, and can be measured in simple morphometric experiments. As *B* is dimensionless, it can be used to make quantitative comparisons between experiments involving very different sizes and growth velocities. This enables universal behaviors and control mechanisms to be identified.

To understand and describe properly the interaction between gravitropism and phototropism, successive steps are undertaken: First, the hypotheses that lead the construction of the model are discussed, then general specifications of the geometry of the organ and of the gravity and light fields are established. In the following part, the construction of model is given. Starting with the simple case of phototropism in isolation which allows to dissect the interaction between photosensitivity and propriosensitivity. Finally, the dynamics of the interaction among photoception, graviception and proprioception are explored, including the influence of light intensity. All these models are investigated for different localization of photoception, at the tip or distributed along the organ. Quantitative experimental data on photogravitropic equilibrium are used to test the validity of the model. These results and the implications for plant biology are discussed.

Materials and Methods {#sec002}
=====================

Before presenting the construction of a dynamic model, it is important to state properly what are the hypotheses that lead our work, as well as the simplified geometries of the organ and of the source fields that are considered. Once these bases are well defined, we will discuss the construction of successive model variants of increasing complexity, first with the simple phototropic model, then with the model describing the interaction between gravitropism and phototropism. For each model two sub-models are considered, depending on the distribution of the photoception along the organ.

Hypotheses {#sec002a}
----------

In order to formulate hypotheses regarding the extension of the AC model to include phototropism, it is now useful to review our knowledge of the distribution of sensing mechanisms and differential growth responses along the plant organ, the mechanisms of photo- gravi- and proprio-ception, as well as their possible interactions.

The localization of gravisensitivity has been established in more detail than the localization of photosensitivity. The plant perception of gravity is related to the presence of statoliths within specialized cells called statocysts. In aerial organs, statocysts and statoliths are found throughout the growth zone, and both perception and the bending response are local \[[@pcbi.1004037.ref009]\]. Our knowledge of the mechanisms of proprioception is limited, but it is known to be exclusively local and likely to involve cytoskeleton remodeling \[[@pcbi.1004037.ref001], [@pcbi.1004037.ref002]\].

In phototropism, blue-light-sensitive phototropins sense the direction of incoming light, but they interact with other photoreceptors in a network that remains to be fully elucidated \[[@pcbi.1004037.ref010]--[@pcbi.1004037.ref013]\]. Photoreceptor localization studies were pioneered by Darwin \[[@pcbi.1004037.ref014]\]. By masking different parts of coleoptiles, Darwin found that the perception of light occurred at the coleoptile apex. As tropic bending occurred all along the growth zone of the organ, a secondary basipetal signal was postulated to be involved, now elucidated as a lateral redistribution of the polar transport of the plant hormone auxin \[[@pcbi.1004037.ref010]\]. Very few studies have further systematically analyzed the localization of photoreceptors or even phototropic responses along plant organs. It has been found, however, that apical phototropic sensing is not universal. For example, the hypocotyl of *Arabidopsis* exhibits photoreceptors all along the growth zone, allowing for distributed local photoception and local phototropic growth responses \[[@pcbi.1004037.ref011], [@pcbi.1004037.ref015]\]. The functional significance of such differences in the localization of photosensitivity among organs has not been investigated.

A few quantitative studies of the effects of continuous light \[[@pcbi.1004037.ref016]--[@pcbi.1004037.ref018]\] have shown that phototropism and gravitropism seem to act in an additive manner. When a plant organ is lit from a direction that differs from that of gravity, the plant undergoes an active movement until a steady-state shape called the photogravitropic equilibrium is reached \[[@pcbi.1004037.ref017], [@pcbi.1004037.ref019]\]. The photogravitropic equilibrium angle (*PGEA*) of the organ tip was found to follow the phenomenological equation $$PGEA = k\mspace{720mu}\log\left( \frac{I}{I_{0}} \right) - g\mspace{720mu}\sin A_{0}$$ where *I* is the fluence rate of illumination, *I* ~0~ is the light-sensing threshold, *g* is the gravitational force, and *A* ~0~ is the initial angle of inclination towards gravity. The term *g*sin*A* ~0~ is the sine law of gravitropism \[[@pcbi.1004037.ref004]\]. This sine law was also used when defining the gravisensing term in the *AC* model for gravitropism1 but in this case it has been reduced to a linear term using the approximation sin(*A*(*s*, *t*)) ≈ *A*(*s*, *t*) + *O*(*A* ^3^). The term *k* *log*(*I*/*I* ~0~) reflects the phototropic stimuli. Both phototropic and gravitropic sensing have been shown to control the relocalization of polar auxin transporters, controlling the formation of lateral gradients in auxin concentration and hence differential growth. These molecular dynamics might explain why the effects of the two stimuli on PGEA are additive \[[@pcbi.1004037.ref007], [@pcbi.1004037.ref020], [@pcbi.1004037.ref021]\]. Thus far, the dynamics of the processes leading to this photogravitropic equilibrium have not been analyzed or modeled in detail.

Building on this current knowledge, we propose herein to extend the graviproprioceptive *AC* modeling approach to study phototropism and its interactions with gravitropism. We focus on investigating dynamic control of photosensitivity, gravisensitivity and propriosensitivity throughout the whole movement.

More precisely, our working hypotheses are as follows:

(H1) The action of the tropic motor is fully driven by the perception-regulation process and results in a change in the local curvature.

(H2) The angles formed between the axis of the growing organ and the respective gravity and light fields influence graviception and photoception, and hence influence gravity- and light components driving the tropic response.

(H3) Proprioception takes place regardless of whether external tropic signals are present. The occurrence of proprioception has been shown for gravitropism \[[@pcbi.1004037.ref002], [@pcbi.1004037.ref022]\]. Concerning phototropism, plants growing in a clinostat (an apparatus that suppresses graviception) were also shown to straighten after a transient light pulse \[[@pcbi.1004037.ref023]\]. Thus, we assume that each constituent element of the organ perceives its own local deformation during active bending, namely the curvature, and responds in order to restore local straightness \[[@pcbi.1004037.ref002]\].

(H4) The different types of perception act additively within the model, meaning that their respective contributions have equivalent roles in driving the movement. It is therefore expected that much of the apparent complexity of the motion is due to spatio-temporal integration of several responses over the changing geometry of the organ with respect to the light and gravity fields.

In order to formulate testable and refutable predictions regarding the dynamic control of tropic movement, we combine these working hypotheses into successive variants of a dynamic model, extending the graviproprioceptive *AC* modeling approach. These models are expected to provide relevant dimensionless numbers that control the dynamics, as well as a means of estimating their values on the basis of experimental measurements. This will enable the model and hypotheses to be evaluated experimentally, and eventually open the way to obtaining accurate phenotypes of tropisms.

We have chosen to neglect the effects related to the elongation of the organ. Recent work has shown that elongation tends to destabilize tropic movements and should increase the oscillations that take place during these movements \[[@pcbi.1004037.ref024]\]. We have found, however, that the proprioceptive sensitivities of modern angiosperms have evolved such that these effects are fully controlled by the plants. Therefore, these growth-related effects can be neglected in the description of tropic movement \[[@pcbi.1004037.ref024]\]. Similarly, even though plants are expected to bend under their own weight, we propose that the effects of a plant's weight on its shape can also be neglected \[[@pcbi.1004037.ref025]\]. This is based on the fact that the *AC* model, which takes into account the plant's perception as a sole driving process, has proved to be sufficient to explain the movement of a wide variety of organs and species of many orders of magnitude of size.

To finish, it is important to note that the *AC* model accounts for orthogravitropism, wherein the organ aligns with the direction of gravity. Some organs align in a direction different from that of gravity, called the gravitropic set-point angle (GSA) \[[@pcbi.1004037.ref026], [@pcbi.1004037.ref027]\]. It is possible to adapt the *AC* model for such cases by modifying the model's graviceptive term, −*β A*(*s*, *t*)→−*β* (*A*(*s*, *t*)−*GSA*) such that the organ aligns with the GSA. However, the *AC* model has not been tested for non-orthogravitropic organs. And this study concentrates on the case of ortho-gravi and photo-tropisms (a very common case in primary shoots).

As stated before, the present study involve successive variants of the model. To unease their handling we named the successive models presented below according to i) the sensory processes that are involved and their locations (biological specification)---e.g., photoproprioceptive; and ii) the driving variables under investigation (mathematical specification) (in accordance with the naming convention in \[[@pcbi.1004037.ref002]\], *e.g.* model AC is driven par the inclination angle A versus the vertical and the curvature C all along the axis, whereas model A^a^C is driven by the apical angle of inclination and by the curvature all along the axis). The reader may find a list of name of the variables and parameters in [Table 1](#pcbi.1004037.t001){ref-type="table"} and of the model in [Table 2](#pcbi.1004037.t002){ref-type="table"}.

10.1371/journal.pcbi.1004037.t001

###### Variables and Parameters.

![](pcbi.1004037.t001){#pcbi.1004037.t001g}

  ---------------------------------- -------------------------------------------------------------------------------------------
  *s*                                Curvilinear abscissa from the base to the apex
  *t*                                Time
  *L*                                Length of the organ
  *L* ~*c*~                          Length of the curved zone
  *A*(*s*, *t*)                      Local angle
  *C*(*s*, *t*)                      Local curvature
  *β*                                Graviceptive sensitivity
  *γ*                                Proprioceptive sensitivity
  *ν*                                Photoceptive sensitivity
  $B = \frac{\beta L}{\gamma}$       Graviproprioceptive number
  $D = \frac{\nu L}{\gamma}$         Photoproprioceptive number
  $M = \frac{\nu}{\beta}$            Photograviceptive number
  *A* ~0~                            Initial angle of the organ
  *A* ~*P*~                          Orientation angle of the light
  $A_{R} = \frac{A_{P}}{1 + M}$      Resultant photogravitropic orientation angle
  *A* ^*a*^ = *A*(*L*, *t*)          apical angle of the organ
  $A_{R}^{a} = A(L,t) - A_{R}$       Resultant photogravitropic orientation apical angle
  *I*                                Intensity of the light
  Φ~*S*~(*I*) = *aI* ^*b*^           relation between the intensity of the light and the observed response. Stevens law.
  Φ~*F*~(*I*) = *c* + *d* log(*I*)   relation between the intensity of the light and the observed response. Weber-Fechner law.
  *GSA*                              Gravitropic Setpoint Angle
  *PGSA* = *A* ~*R*~                 PhotoGravitropic Setpoint Angle
  *PGEA*                             PhotoGravitropic Equilibrium Angle
  ---------------------------------- -------------------------------------------------------------------------------------------

10.1371/journal.pcbi.1004037.t002

###### Models.

![](pcbi.1004037.t002){#pcbi.1004037.t002g}

  *name*                              *notation*      *equation*
  ----------------------------------- --------------- ---------------------------------------------------------------------------------------------
  graviproprioceptive                 *AC*            $\frac{\partial C(s,t)}{\partial t} = - \beta A(s,t) - \gamma C(s,t)$
  photoproprioceptive apical          *A* ^*a*^ *C*   $\frac{\partial C(s,t)}{\partial t} = - \nu A(L,t) - \gamma C(s,t)$
  photograviproprioceptive            *A* ~*R*~ *C*   $\frac{\partial C(s,t)}{\partial t} = - \beta A(s,t) - \nu(A(s,t) - A_{P}) - \gamma C(s,t)$
  photo-gravi-proprioceptive apical   $A_{R}^{a}C$    $\frac{\partial C(s,t)}{\partial t} = - \beta A(s,t) - \nu(A(L,t) - A_{P}) - \gamma C(s,t)$

The different models used in this paper, with their full names and their respective notations.

Geometrical specifications of the organ and of the source fields {#sec002b}
----------------------------------------------------------------

One drawback of most studies on phototropism is that they have overlooked the true geometry of the organ, measuring only the orientation of the apical tip of the organ \[[@pcbi.1004037.ref028]--[@pcbi.1004037.ref031]\]. One of the greatest insight of the *AC* model was to show that the geometry of the organ and of the source field is central to describe the movement properly. It is then important to have a clear description of the studied geometry.

Like the *AC* model, the models developed herein describe the shape of the organ in terms of its median line, i.e., its central axis ([Fig. 1](#pcbi.1004037.g001){ref-type="fig"}). We parameterize the position along this median line according to the curvilinear abscissa *s* going from the base *s* = 0 to the apex *s* = *L*. The angle *A*(*s*, *t*) then describes the local orientation of the median with respect to the vertical pointing up (and hence to the direction of the gravity vector **g**) at time *t*. The angle of the apex at the tip of the organ, *A* ^*a*^, is then equal to *A*(*L*, *t*). The local curvature *C*(*s*, *t*) is the spatial rate of change of *A*(*s*, *t*) along *s*, and from differential geometry we know that $$C\left( s,t \right) = \frac{\partial A\left( s,t \right)}{\partial s}\mspace{720mu}\mspace{720mu}{or}\mspace{720mu}\mspace{720mu} A\left( s,t \right) = A_{0} + \int_{0}^{s}C\left( l,t \right)dl$$

![Geometric description of organ shape.\
The median line of an organ of total length *L* is in a plane defined by coordinates *x*, *y*. The arc length *s* is defined along the median line with *s* = 0 referring to the base and *s* = *L* referring to the apex. *A*(*s*) is the local orientation of the organ with respect to the vertical, and *C*(*s*) the local curvature. The orientation of the gravity field vector *g* is parallel to the *y* axis. The orientation of the light field vector *l* forms an angle *A* ~*P*~ with the *y* axis. The resulting angle of the two tropisms is called *A* ~*R*~.](pcbi.1004037.g001){#pcbi.1004037.g001}

At the scale of the plant, gravity acceleration **g** is homogeneous, uniform, and invariant to translation and to rotation around the main direction of the gravity field ([Fig. 1](#pcbi.1004037.g001){ref-type="fig"}). The case of light perception seems to be less straightforward. The geometry of light fields in nature can be quite diverse. At the scale of the plant, the light field produced by the sun on the earth is supposed to be homogeneous, uniform, and invariant to translation and to rotation around the main direction of the light field, but its orientation relative to the gravity field varies over the course of the day. In the lab, a point light source can emit a spherical field whose intensity (irradiance) decreases with the square of the distance. Furthermore, multiple point sources might be used, thereby increasing the complexity of the analysis. Thus far, most controlled experiments have used a distant punctual source with a collimated beam \[[@pcbi.1004037.ref016]--[@pcbi.1004037.ref018], [@pcbi.1004037.ref032]\]. For simplicity, we assume that the light field **l** has the same symmetry as the gravity field. This assumption holds true if the source of light is far from the organ. To further reduce the complexity of the analysis, we do not consider curvature outside the plane; that is, we explore the case in which the main direction of the plant's organ, the gravity field **g** and the light field **l** are all in the same plane *P* ~*gl*~ ([Fig. 1](#pcbi.1004037.g001){ref-type="fig"}). This again was the case in most controlled experimental studies (e.g. \[[@pcbi.1004037.ref017]\]). Therefore, the direction of the lighting in the *P* ~*gl*~ plane is defined by the angle *A* ~*p*~ and by a radiance (or light intensity) *I* ~0~ ([Fig. 1](#pcbi.1004037.g001){ref-type="fig"}).

The simplifications we adopt enable us to assess the model by comparing it with experimental data, while still providing insights regarding the regulation of the movement. Moreover, the assumption that the initial direction of the organ, the direction of gravity and that of the light are in the same plane is fulfilled for organs that are already aligned with the direction of gravity.

Optical aspects of photoception {#sec002c}
-------------------------------

Coleoptiles and other plant shoots appear to sense light direction via the light gradient across the organ \[[@pcbi.1004037.ref033]\]. Treatments that change the steepness of this gradient (e.g., infiltration with dyes or changes in the amounts of natural pigments) alter the phototropic response: the steeper the gradient, the more the coleoptile bends towards the light. This gradient depends on internal light transfer through plant tissues but also depends linearly on the irradiance of the light impinging on one side of the surface of an organ. This irradiance is given by Lambert's cosine law of irradiance of geometrical optics (e.g. \[[@pcbi.1004037.ref017]\]), which states that the irradiance falling on any surface varies in proportion to the cosine of the incident angle. $$I\left( s,A_{P} \right) = I_{0}\mspace{720mu}\cos\left( {\frac{\pi}{2} - \left( A\left( s,t \right) - A_{P} \right)} \right) = - I_{0}\mspace{720mu}\sin\left( {A\left( s,t \right) - A_{P}} \right)$$

As the *AC* model is a first-order model, and in the limit of small angles *A*(*s*, *t*)−*A* ~*P*~, we can use the approximation sin(*A*(*s*, *t*)−*A* ~*P*~) ≈ *A*(*s*, *t*)−*A* ~*P*~ + *O*(*A* ^3^), so that the photosensitivity is proportional to the angle between the organ and the light direction i.e., *A*(*s*, *t*)−*A* ~*P*~ (the influence of the intensity of the light irradiance will be considered in a later section). Although this approximation is only valid for small angles, such an approximation gives a good and efficient understanding of the dynamics of the system. The influence of each term on the tropic movement can be easily understood and discussed.

Phototropism {#sec002d}
------------

We may now start with the simplest case of phototropism without gravitropism. Using [equation 4](#pcbi.1004037.e014){ref-type="disp-formula"}, a phototropic model can be constructed. Two sub-models are considered, when the perception is apical and when the perception is local.

### Apical photoception: *A^a^C* model {#sec002da}

In cases where perception of light is found to be purely apical, the previous dependency, [equation 4](#pcbi.1004037.e014){ref-type="disp-formula"}, is expressed as *A*(*L*, *t*)−*A* ~*P*~. The photosensitivity corresponding to the inclination of the apex in relation to the light beam is translated into a secondary signal, propagated basipetally along the organ through asymmetric polar transport of the plant hormone auxin \[[@pcbi.1004037.ref010]\]. As auxin transport is faster than the characteristic time of the movement, this long-distance signaling along the organ can be approximated as instantaneous, so that, at any position *s*, a small segment of the organ is submitted to a signal proportional to *A*(*L*, *t*)−*A* ~*P*~ ([SI](#sec005){ref-type="sec"}). Indeed, it can be shown that the propagation of the signal has only slight effects on the dynamics and the steady state of the organ, as long as the characteristic time of propagation (the time for the signal to go from the apex to the base) is smaller than the characteristic time of movement (given by the elongation rate \[[@pcbi.1004037.ref025]\]) ([SI](#sec005){ref-type="sec"}). The perceived signal is thus the value at the apex, *A*(*L*, *t*), and this is the only value accessible along the organ.

Let us now assume for a moment that the light is zenithal, *A* ~*P*~ = 0, so that the phototropic reaction tends to bring the tip to the vertical.

The conditions of symmetry described in \[[@pcbi.1004037.ref002]\] are then fulfilled. The behavior remains the same, according to the change *A*(*s*, *t*) → −*A*(*s*, *t*) and then *C*(*s*, *t*) → −*C*(*s*, *t*).

According to hypotheses H1--4, in the case of zero gravity or zero graviception (agravitropic mutants), the most basic linear equation for an apical perception can then be formulated as $$\frac{\partial C\left( s,t \right)}{\partial t} = - \nu A\left( L,t \right) - \gamma C\left( s,t \right)$$ where *νA*(*L*, *t*) is the apical photoceptive term, readily transmitted all along the organ, and *γ C*(*s*, *t*) is the proprioceptive term. This will be referred to as the apical photo-proprioceptive model, or the *A^a^C* model, with ^*a*^ expressing the apical perception of the stimulus (see [Table 1](#pcbi.1004037.t001){ref-type="table"}). As noted, this model assumes that the transmission is instantaneous. Every part of the organ is able at any point to perceive the exact orientation of the apex. It is also assumed that the transmission is homogenous along the coleoptile, meaning that *ν* is not dependent on the position *s* (see [SI](#sec005){ref-type="sec"}).

Now let us assume the following initial conditions: a straight but tilted organ and the boundary conditions of perfect basal clamping: $$A\left( s,0 \right) = A\left( 0,t \right) = A_{0}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu} C\left( 0,t \right) = 0$$ It is easy to see that the *A^a^C* model is in fact independent of space coordinate *s*. First, the organ is initially straight, so the initial curvature is the same everywhere. Then as the perception of the light is purely apical, at any given time all parts of the organ receive the same signal, proportional to the apical orientation at that time. The variation of curvature, and hence the curvature, are the same at any point of the organ. The solution of the *A^a^* model is then given by $$A(s,t) = A(t)\mspace{180mu} = \mspace{180mu} A_{0}\left( {1 - \frac{s}{L}\frac{1}{1 + \frac{\gamma}{\nu L}}\left( {1 - e^{- {({\nu L + \gamma})}t}} \right)} \right)$$ $$C(s,\mspace{180mu} t) = C(t)\mspace{180mu} = \mspace{180mu} - \frac{A_{0}}{L}\frac{1}{1 + \frac{\gamma}{\nu L}}\left( {1 - e^{- {({\nu L + \gamma})}t}} \right)$$

It is possible to define a dimensionless number *D*, the photoproprioceptive number, that expresses the ratio between photoception and proprioception: $$D = \frac{\nu L}{\gamma}$$

### Local photoception: *AC* model {#sec002db}

If the perception of light is set to be local (as observed, for example, in *Arabidopsis* hypocotyls), [equation 5](#pcbi.1004037.e015){ref-type="disp-formula"} should be modified to the local photoceptive equation $$\frac{\partial C\left( s,t \right)}{\partial t} = - \nu A\left( s,t \right) - \gamma C\left( s,t \right)$$ with the same initial conditions as described in [equation 6](#pcbi.1004037.e016){ref-type="disp-formula"}. [Equation 11](#pcbi.1004037.e021){ref-type="disp-formula"}, which describes the photoceptive model, is now strictly equivalent to the graviproprioceptive [equation 1](#pcbi.1004037.e001){ref-type="disp-formula"}, i.e., the *AC* model \[[@pcbi.1004037.ref002]\].

Photo-gravitropic interaction {#sec002e}
-----------------------------

The previous model can now be easily extended to take into account the perceptions of both light and gravity. As described in [Fig. 1.A](#pcbi.1004037.g001){ref-type="fig"}, the direction of the light stimulus is now considered to form an angle *A* ~*P*~ with the vertical, the direction of the gravity. The orientation angle *A*(*s*, *t*) along the organ is still measured with respect to the vertical, and hence to the direction of gravity.

If *A* ~*P*~ ≠ 0 the symmetry of the system is broken. When only one of the two tropisms influences the system, rotation around the axis defined by the direction of the corresponding field should have no effect on the system. When the two tropisms act simultaneously, however, the angle between the two fields prevents such global symmetry. It should also be noted that if the gravitropic term or the phototropic term tends to 0, the system converges to the photo-proprioceptive models outlined above.

As for the phototropic models, effects due to local and apical perception are considered and discussed.

### Local perception: *A~R~C* model {#sec002ea}

When the organ lies in the plane defined by the direction of gravity and of light, the local equation of purely local photogravitropism, called the local photo-gravi-proprioceptive equation, is given by $$\frac{\partial C\left( s,t \right)}{\partial t} = - \nu\left( {A\left( s,t \right) - A_{P}} \right) - \beta A\left( s,t \right) - \gamma C\left( s,t \right)$$ The two external tropisms act in the same way through perception of the local angle, except for the constant angle *A* ~*P*~, indicating the direction of phototropism.

Simple variable substitution can be performed to simplify [equation 12](#pcbi.1004037.e022){ref-type="disp-formula"}: $$A^{\prime}\left( s,t \right) = A\left( s,t \right) - A_{P}\frac{\nu}{\nu + \beta} = A\left( s,t \right) - A_{R}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu} C^{\prime}\left( s,t \right) = \frac{dA^{\prime}\left( s,t \right)}{dt} = C\left( s,t \right)$$ where *A* ~*R*~ is a Resultant Angle linked to the balance of the two tropic factors through combined gravi- and photo-ception $$A_{R} = A_{P}\frac{\nu}{\nu + \beta}$$ [Equation 12](#pcbi.1004037.e022){ref-type="disp-formula"} can be rewritten in a more compact form $$\frac{\partial C\left( s,t \right)}{\partial t} = - \left( \nu + \beta \right)A^{\prime}\left( s,t \right) - \gamma C\left( s,t \right)$$ with no constant term. Starting with the set of initial conditions defined in [equation 6](#pcbi.1004037.e016){ref-type="disp-formula"}, the initial conditions to be considered are then $$A^{\prime}\left( s,0 \right) = A^{\prime}\left( s,t \right) = A_{0}^{\prime} = A_{0} - A_{R}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu}\mspace{720mu} C^{\prime}\left( 0,t \right) = C\left( 0,t \right) = 0$$

The variable substitutions and calculations shown above are equivalent to the case in which the plant has been rotated at an angle *A* ~*R*~, such that the two tropisms then act together in the same direction. The local photo-gravi-proprioceptive [equation 15](#pcbi.1004037.e025){ref-type="disp-formula"} can thus be named the *A~R~C* model.

Another dimensionless number *M*, called the photograviceptive number, can then be defined as the ratio between the gravisensitivity and the photosensitivity: $$M = \frac{\beta}{\nu}$$ Note that *M* is not independent from the graviproprioceptive number *B* and the photoproprioceptive number *D* but instead can be expressed as $$M = \frac{B}{D}$$

### Apical photoception: $A_{R}^{a}C$ model {#sec002eb}

In the case in which the perception of light is apical, the equation that describes this apical-photo/local-gravi-proprioception-driven movement, called the $A_{R}^{a}C$ model, is given by $$\frac{\partial C\left( s,t \right)}{\partial t} = - \nu\left( {A\left( L,t \right) - A_{P}} \right) - \beta A\left( s,t \right) - \gamma C\left( s,t \right)$$ Using the variable substitutions shown in [equation 13](#pcbi.1004037.e023){ref-type="disp-formula"}, the dynamical [equation 19](#pcbi.1004037.e031){ref-type="disp-formula"} can be simplified to $$\frac{\partial C\left( s,t \right)}{\partial t} = - \nu A^{\prime}\left( L,t \right) - \beta A^{\prime}\left( s,t \right) - \gamma C\left( s,t \right)$$ This equation converges to a steady-state shape that is described as $$A^{\prime}\left( s,t\rightarrow\infty \right) = A_{0}^{\prime}\left( {e^{- Bs/L} - e^{- B}\frac{1 - e^{- Bs/L}}{M + \left( {1 - e^{- B}} \right)}} \right)$$ $$C\left( s,t\rightarrow\infty \right) = - A_{0}^{\prime}\frac{B}{L}e^{- Bs/L}\left( {1 + \frac{e^{- B}}{M + \left( {1 - e^{- B}} \right)}} \right)$$

### Apical photoception and light intensity in the *A~R~C* model {#sec002ec}

The sensitivity of the outputs of the photo-gravi-proprioceptive models *A~R~C* and $A_{R}^{a}C$ can vary as a function of the photoceptive, graviceptive or proprioceptive sensitivities (captured in the dimensionless ratios *D*, *B*, and *M*). Variations in these sensitivities could be due, for example, to natural or artificial genetic variations within or across species, such as those reported by \[[@pcbi.1004037.ref002]\] for *B*. However, the photograviceptive number, *M*, may also depend on a current environmental factor, namely, the intensity of the light *I* (in addition to the angle with the incident light direction *A* ~*p*~). Many studies have indeed reported that the light fluence rate (irradiance) *I* influences the rate and strength of the phototropic reaction (e.g \[[@pcbi.1004037.ref017], [@pcbi.1004037.ref034]\]). It is thus expected that the photoceptive term *ν* in our model, and therefore the photograviceptive number *M*, should be a function of *I*, the fluence rate of the incident light beam. For example, \[[@pcbi.1004037.ref034]\] proposed a biochemical model of photoception in which the rate constant of the primary photosensitive reaction is proportional to the fluence rate *I* reaching the photosensitive tissues. The situation may be even more complex, as there have been reports of a "tonic" effect of light intensity on gravisensitivity (acting through other photoreceptors besides phototropins) \[[@pcbi.1004037.ref008]\]. Therefore, in cases in which the intensity of the incident light is subject to change, it seems that at least the photograviceptive number *M* should depend on the intensity of the incident light. Clearly, however, this dependence has yet to be quantitatively and experimentally elucidated. Several forms for the intensity dependency has been proposed in the literature but there are no clear theoretical or empirical reasons for choosing one specific formulation over another. The photobiology literature has classically assumed that this relation is given by the linear relation between the stimulus and the plant's response, the Bunsen-Roscoe reciprocity law \[[@pcbi.1004037.ref035]--[@pcbi.1004037.ref037]\]. However, this law seems to be valid for only a small range of light dosage levels. The Schwarzschild law (a power-law generalization of the reciprocity law) appears to more accurately describe the relationship between the response and the intensity of the light stimuli, for a wide range of materials and systems \[[@pcbi.1004037.ref038]\] and has been found to be a fairly generic approximation, where the Bunsen-Roscoe reciprocity law can be seen as a special case of the power law. Galland \[[@pcbi.1004037.ref017]\] used a logarithm relation to express the relation between the observed *PGEA* and the intensity of the light, but \[[@pcbi.1004037.ref017]\] does not fully assess or discuss whether the assumption of a logarithmic relation is indeed preferable to the use of a power law.

The lack of an unequivocal formulation of the dependency on the intensity of the signal is also reflected in the broader domain of sensory biology, where two competing general models of the relation between sensation and stimuli from sensory biology are considered. The first is the Weber-Fechner Law \[[@pcbi.1004037.ref039]\], which postulates that the intensity of the response is proportional to the logarithm of the stimulus. The second is Stevens' law, which postulates a power law relationship between the stimulus and the sensation \[[@pcbi.1004037.ref040]\]. Although these relations have been assessed primarily in the domain of human behavior, quantitative similarities in sensing relations between animals and plants have been identified \[[@pcbi.1004037.ref022]\]. The question of which model better describes the effect of the intensity of the stimulus (in our case light) thus remains open, and the adoption of a given law seems to be mostly dependent on the context and on the conventions of the discipline. Therefore, in what follows, we will consider the two types of formulations (log or power law) in the following and will subsequently compare them using experimental data.

We may now extend the photograviproprioceptive model by incorporating a new working hypothesis stating that

(H5) the photograviceptive number *M* is a function of the fluence rate of the incident light, i.e., $$M \sim \Phi\left( I \right)$$

This yields light-intensity-dependent photo-gravi-proprioceptive models *A~R~*(*I*)*C* and $A_{R}^{a}(I)C$. The function Φ can be expressed either by a power law $$\Phi_{S}\left( I \right) = aI^{b}$$ or by logarithmic relation $$\Phi_{F}\left( I \right) = c + d\mspace{180mu}\log\left( I \right)$$

Experiments {#sec002f}
-----------

A series of experiments were conducted or reprocessed from the bibliography in order to assess the Photo-Gravitropic models against experimental data in continuous light and 1g conditions (i.e.reminiscent of natural outdoor conditions on Earth): an experiment in the typical set-up used for phototropic studies, but tracking the entire kinematics of the movement (LLE experiment) and a set of experiments conducted by \[[@pcbi.1004037.ref017]\] and focusing only on the PhotoGravitropic Equilibrium tip Angle (PGEA), but at various tilting angles and light intensities (referred to as "Galland's experiments" in the following).

### Lateral light experiments with Wheat coleoptile (LLE experiment) {#sec002fa}

In order to asses the effect of lateral light on vertical organs, experiments were conducted on 12 etiolated wheat coleoptiles (Recital). Wheat seeds were grown in in vermiculite in the dark at 23°*C* until their size reached at least 2*cm* (approximatively 72h). In order to avoid effects of already curved coleoptiles, the straightest coleoptiles were selected for the experiment. The coleoptiles were then displayed in front of a camera (Nikon DS5200) with a source of blue light (LED VAOL-5701SBY4-ND collimated with a parabolic mirror RCphotonics) perpendicular to the vertical. Timelapse were taken every five minutes for 24h.

### Galland experiments {#sec002fb}

Galland carried out a detailed quantitative study of the relationship between the *PGEA* ([equation 2](#pcbi.1004037.e002){ref-type="disp-formula"}) and the fluence rate *I* in light-grown *Avena* coleoptiles under continuous light \[[@pcbi.1004037.ref017]\], hence contrary to our experiment the detailed kinematics were not reported. In these experiments, the coleoptiles were tilted at different angles from the vertical while the tip of each coleoptile was illuminated orthogonally and in the plane of bending at different fluence rates. The experiments comprised 15 replicates. Data have been reprocessed from \[[@pcbi.1004037.ref017]\]. Two experimental protocols were used, called **PROT1** and **PROT2** here.

**PROT1** The coleoptiles were tilted at different initial angles *A* ~0~, while the light was always maintained perpendicular to the initial orientation of the organ, *A* ~*P*~ = *A* ~0~ + *π*/2. The dependency of the *PGEA* on the fluence rate of the incident light *I* was then assessed by varying *I* at each value of the initial angle *A* ~0~. *PGEA* can be plotted as a function of *I* for the different values of *A* ~0~ (and the corresponding *A* ~*P*~ = *A* ~0~ + *π*/2).

**PROT2** The position of the *PGEA* was also measured for different values of *A* ~0~, where *A* ~*P*~ = *A* ~0~ + *π*/2. However, the fluence rate of the incident light *I* was tuned experimentally so that the angle of the straight organ was immediately steadied at *A* ~*R*~ = *A* ~0~ stated in terms of the *A~R~C* model,. This provides an estimate of the value of *I*(*A* ~*R*~ = *A* ~0~) at photogravitropic equilibrium. The results of both protocols can be used to independently assess the $A_{R}^{a}(I)C$ model.

The data produced in \[[@pcbi.1004037.ref017]\] according to the two protocols have been reprocessed. Linear orthogonal fit was carried out in a log-log plot and in a semi-log plot. A straight line in a log-log plot (as in [Fig. 2.A](#pcbi.1004037.g002){ref-type="fig"}) would reflect a power law, whereas a linear fit in a semi-log plot (as in [Fig. 2.B](#pcbi.1004037.g002){ref-type="fig"}) would reflect a logarithmic relation.

![A. −1 + *A* ~*P*~/*A* ~*R*~ as a function of the fluence rate of the light.\
Data are reprocessed from Figs. [4](#pcbi.1004037.g004){ref-type="fig"}, [5](#pcbi.1004037.g005){ref-type="fig"} and [6](#pcbi.1004037.g006){ref-type="fig"} in \[[@pcbi.1004037.ref017]\] Solid lines correspond to the fit $log( - 1 + \frac{A_{P}}{A_{R}}) = a^{\prime} + b^{\prime}\log(I)$. An etiolated *Avena* coleoptile is tilted from the vertical for different values of the angle *A* ~0~ (0, 10, 30, 90 and 120) while a light beam is collimated on the tip of the coleoptile. For different fluence rates, controlled by a neutral density filter, the apical angle is measured after variable duration. This is said to be the equilibrium angle. B. −1 + *A* ~*P*~/*A* ~*R*~ as a function of the fluence rate of the light; data are reprocessed from Figs. [4](#pcbi.1004037.g004){ref-type="fig"}, [5](#pcbi.1004037.g005){ref-type="fig"} and [6](#pcbi.1004037.g006){ref-type="fig"} of \[[@pcbi.1004037.ref017]\]. Solid lines correspond to the fit $- 1 + \frac{A_{p}}{A_{R}} = a^{\prime} + b^{\prime}\log(I)$. C. −1 + *A* ~*p*~/*A* ~*R*~ as a function of the fluence rate of the light; data are reprocessed from Figure 8 in in \[[@pcbi.1004037.ref017]\]. The solid lines correspond to the fits from A. The empty black symbols (circles and squares) correspond to the measured fluence rates that precisely compensate the gravitropic reaction. The equilibrium state then corresponds to the case in which no movement is observed and *A* ~*R*~ = *A* ~0~. Empty circles: *A* ~0~\<90°; empty squares: *A* ~0~\>90°. D. *A* ~*p*~/*A* ~*R*~ as a function of the fluence rate of the light. Solid lines correspond to the fits from A.](pcbi.1004037.g002){#pcbi.1004037.g002}

Results {#sec003}
=======

First, a simple experiment on wheat coleoptile with lateral light (LLE) is presented to get a rough estimate of the expected behavior during combined photo- and gravi-tropism movements.

Then to get a better estimate of the validity of the photogravitropic model, measurements of the apical angle of coleoptiles \[[@pcbi.1004037.ref017]\]) are compared with the prediction of the model. Besides providing a good agreement with the *A* ~*R*~ *C* model, this agreement is only made for one specific organ. The different models are providing a more richer diversity that needs to be discussed.

The effect of each parameter in the model needs to be discerned as well as the respective influence of local and apical perception. To this end, the dynamics of the model of increasing complexity are presented.

Finally to obtain a full understanding of the model $A_{R}^{a}C$, it is proposed to study the different limits of the model. When one parameter is negligible or is dominating the dynamic, the specific effect of each parameter is accentuated. A deep, qualitative understanding of a complex system can then easily be reached.

The steady states of the different models are summarized in [Table 3](#pcbi.1004037.t003){ref-type="table"}.

10.1371/journal.pcbi.1004037.t003

###### Steady states of the different models.

![](pcbi.1004037.t003){#pcbi.1004037.t003g}

  *notations*     *steadystate*
  --------------- --------------------------------------------------------------------------------------------------------------------------------------------------------------
  *AC*            *A*(*s*) = *A* ~0~ *e* ^−*βs*/*γ*^
  *A* ^*a*^ *C*   $A(s) = A_{0}\left( {1 - \frac{s}{L}\frac{1}{1 + \frac{\gamma}{\nu L}}\left( {1 - e^{- {({\nu L\mspace{180mu} + \mspace{180mu}\gamma})}t}} \right)} \right)$
  *A* ~*R*~ *C*   $A(s) = (A_{0} - A_{R})e^{- (\beta\mspace{180mu} + \mspace{180mu}\nu)s/\gamma} + A_{R}$
  $A_{R}^{a}C$    $A(s) = (A_{0} - A_{R})\left( {e^{- Bs/L} - e^{- B}\frac{1 - e^{- Bs/L}}{M + \left( {1 - e^{- B}} \right)}} \right) + A_{R}$

The steady state, $\frac{\partial C}{\partial t} = 0$, of the different models.

Experimental assessment of the *A~R~C* model {#sec003a}
--------------------------------------------

A major prediction of the of *A* ~*R*~-based models is that even when photoception and graviception act in different directions (i.e., *A* ~*P*~ ≠ 0), they seem to act together to bring the organ towards a single direction defined by the angle *A* ~*R*~. Additionally the *A~R~C* model predicts that the whole coleoptile curves in the direction *A* ~*R*~ which is first achieved at the tip and the curvature then concentrates near the base.These two predictions are independent of the position at which light is perceived, i.e., at the tip of the organ or along the entire organ, and is thus common to all the variants of the *A* ~*R*~ *C* model, (*i.e.* the local photoception model *A* ~*R*~ *C* and apical photoception model $A_{R}^{a}C$ and their versions including explicit light intensity effects *A~R~*(*I*)*C* and $A_{R}^{a}(I)C$). For the sake of simplicity, whenever this common core is addressed, the notation *A* ~*R*~ *C* will be used to refer to all of them.

When illuminated by a source of light from a direction that differed from that of gravity (LLE experiment), the coleoptile curved and finally reached a steady state. The apical part aligned with an angle that was an intermediary between the direction of gravity (*A* = 0) and the direction of the light (*A* = *π*/2) ([Fig. 3.A and B](#pcbi.1004037.g003){ref-type="fig"}). This shows qualitative agreement with the *A~R~C* model. Additionnally the transient kinematic pattern in which i) the whole coleoptile curves until the direction *A* ~*R*~ is first achieved at the tip and ii) the curvature then concentrates near the base. is qualitatively observed in ([Fig. 3.B and C](#pcbi.1004037.g003){ref-type="fig"}). However, this kinematic experiment also revealed some limits of grass coleoptile as a "model system" for tropic studies. Indeed, when a coleoptile is illuminated, photomorphogenetic effects take place, slowing down the expansion of the coleoptile and increasing that of the inner leaf. After approximately 10*h*, the leaf inside the coleoptile pierces the coleoptile. The dynamics of the coleoptile then ceases, and there is no further elongation of coleoptile. The dynamics of the movement is then fully determined by the leaf that was inside and this dynamics differs from that of the coleoptile. Moreover the influence of the inner leaf may explain why the rather steady tip angle (between 3 and 6 h in the example [Fig. 3.A and B](#pcbi.1004037.g003){ref-type="fig"}) is subsequently varied again (between 7 and 10h), and it is thus difficult to determine whether the coleoptile has had time to reach a steady state before been disrupted. However, the apical tip is predicted in the model to converge to the direction defined by *A* ~*R*~ before the shape reaches a steady state, and we saw that this is consistent with what is observed experimentally. We may thus assume that in Galland \[[@pcbi.1004037.ref017]\] estimate of *A* ~*R*~ could be measured even when the full steady state of the rest of the coleoptile could not be well defined.

![A. Timelapse of a coleoptile during phototropic movement, 50 minutes between each picture.\
The source light is collimated on the right of the organ and is perpendicular to the vertical. The white bar is 1*cm*. B. The kinematics of the angle *A*(*s*, *t*) plotted with respect to time *t* and curvilinear abscissa *s* (from the base to the apex of the organ). C. The kinematics of the curvature *C*(*s*, *t*) plotted with respect to time *t* and curvilinear abscissa *s* (from the base to the apex of the organ). At first, the whole organ curves; then the curvature concentrates near the base.](pcbi.1004037.g003){#pcbi.1004037.g003}

We therefore tested further the validity of the functional form of *A* ~*R*~-based models through the analysis of the apical tip angle at PhotoGravitropic Equilibrium (*PGEA*) in the work by Galland \[[@pcbi.1004037.ref017]\]. If the *A* ~*R*~ *C* core-model is correct, then the experimental *PGEA* should be equal to *A* ~*R*~, and according to [equation 14](#pcbi.1004037.e024){ref-type="disp-formula"}, it should depend on the angle between the direction of gravity and direction of the light *A* ~*P*~, and on the ratio between the gravisensitivity and the photosensitivity levels.

Indeed *A* ~*R*~ can be expressed in the model as a function the photograviceptive number *M*: $$A_{R} = A_{P}\frac{1}{1 + M}$$ It follows directly that when photoception dominates, the organ bends in the direction of the light, *M* \< \< 1, *A* ~*R*~ = *A* ~*P*~. However when graviception dominates, *M* \> \> 1, *A* ~*R*~ = 0 ([Fig. 4](#pcbi.1004037.g004){ref-type="fig"}). It is then possible to express *M* directly as a function of *A* ~*R*~/*A* ~*P*~. $$\begin{array}{r}
{M = \frac{A_{P}}{A_{R}} - 1} \\
\end{array}$$

![Variation of the reference orientation *A* ~*R*~ as a function of *M*.\
When photoception dominates, *M* \<\< 1, *A* ~*R*~ = *A* ~*P*~. When graviception dominates, *M* \>\> 1, *A* ~*R*~ = 0.](pcbi.1004037.g004){#pcbi.1004037.g004}

A direct quantitative assessment of the functionnal form of the *A~R~*(*I*)*C* and $A_{R}^{a}(I)C$ models---and of the competing *M* ∼ Φ(*I*) sub-models---vs. experimental data can then be conducted through the analysis of the experimental relation between the steady-state tip angle *PGEA* and the fluence rate of the incident light *I*. Indeed *A* ~*P*~ is known, and we have seen that we may assume *A* ~*R*~ = *PGEA*. Then using [equation 27](#pcbi.1004037.e055){ref-type="disp-formula"} an estimate of M can be obtained. And combining [equation 27](#pcbi.1004037.e055){ref-type="disp-formula"} with [equation 23](#pcbi.1004037.e037){ref-type="disp-formula"}, a prediction of the model is then that a plot of *M* = *A* ~*P*~/*A* ~*R*~ − 1 as a function of *I* should display either a power-law dependancy, [equation 24](#pcbi.1004037.e038){ref-type="disp-formula"}, or a log-dependancy, [equation 25](#pcbi.1004037.e050){ref-type="disp-formula"}. Moreover, these relations should fit to a single curve *M* ∼ Φ(*I*), independent of the initial tilting angle *A* ~0~. Note however that due to the linearization of the sine terms of the angle dependency of gravi- and photoception (equations [1](#pcbi.1004037.e001){ref-type="disp-formula"} and [4](#pcbi.1004037.e014){ref-type="disp-formula"}) this should be only valid for inclination angles *A* ~0~ \< 90°. [Fig. 2](#pcbi.1004037.g002){ref-type="fig"} shows the compiled data from the first experimental protocol (**PROT1**) \[[@pcbi.1004037.ref017]\], reprocessed for comparison with the predictions of [equation 27](#pcbi.1004037.e055){ref-type="disp-formula"}. The results of fitting equations [24](#pcbi.1004037.e038){ref-type="disp-formula"} and [25](#pcbi.1004037.e050){ref-type="disp-formula"} to the experimental data are also shown in [Fig. 2 A and B](#pcbi.1004037.g002){ref-type="fig"}. And [Table 4](#pcbi.1004037.t004){ref-type="table"} provides the values of the coefficients *a* and *b* of the curves fitting the power law ([equation 24](#pcbi.1004037.e038){ref-type="disp-formula"}) for each initial angle *A* ~0~. In the log-log plot in [Fig. 2.A](#pcbi.1004037.g002){ref-type="fig"}, the results of each individual experiment at a given *A* ~0~ fit well to a straight line *M* ∼ *I* ^−*b*^ with values of *b* between 0.36 and 0.44 ([Table 4](#pcbi.1004037.t004){ref-type="table"}). Except in cases in which *A* ~0~\>90°, the coleoptile response is independent of the initial angle *A* ~0~, as predicted by the *A~R~C* model (in [Fig. 2.A](#pcbi.1004037.g002){ref-type="fig"}), and these data collapse on a single master curve. On the semi-log plot, however, no convincing fit can be identified. The power law relationship therefore seems to better describe the results than the logarithmic relationship does. Under the power law dependency of *M* on light intensity *I*, the coefficient of determination was *R*2 ∼ 0.91, such that our model captured 91% of the total changes in the steady-state tip angle attributable to changes in the initial angle *A* ~0~ and fluence rate *I*.

10.1371/journal.pcbi.1004037.t004

###### Parameters of the fit $log(1 - \frac{A_{P}}{A_{R}}) = a^{\prime} + b^{\prime}\log(I)$.

![](pcbi.1004037.t004){#pcbi.1004037.t004g}

  *A* ~0~   *a* ^′^   *b* ^′^   *R* ^2^
  --------- --------- --------- ---------
  0         18        0.36      0.88
  10        47        0.40      0.95
  30        43        0.41      0.94
  90        10        0.44      0.72
  120       28        0.41      0.89

Value of the parameters *a* ^′^ and *b* ^′^ given by the fit $log(1 - \frac{A_{P}}{A_{R}}) = a^{\prime} + b^{\prime}\log(I)$ for different values of the initial angle *A* ~0~ of the plants ([Fig. 2A](#pcbi.1004037.g002){ref-type="fig"}).

The response for higher angles does not follow the first order relation expressed in the *AC* model anymore, but the sine laws. It is then expected that, for inclination angles *A* ~0~ \> 90°, the experimental data are adrift from this fit and vary non-monotonously, as *A* varies in a way that can no longer be accounted for by the *A~R~C* model. This is indeed observed in [Fig. 2](#pcbi.1004037.g002){ref-type="fig"}.

The second protocol (**PROT2**) provides a different way to measure the same parameter than **PROT1**, even if the protocol is different. The plot of *A* ~*R*~/*A* ~*P*~ = *A* ~0~/*A* ~*P*~ as a function of the experimental values of *I* should then collapse on the single master curve defined with **PROT1** in [Fig. 2.C](#pcbi.1004037.g002){ref-type="fig"}. Again, the model predictions were fair at small angles, *A* ~*R*~/*A* ~*P*~ ≤ 0.4, but for larger angles the experimental data diverged non-monotonously from the model. It is interesting to note that the results fit well to both the logarithmic law and the power law \[[@pcbi.1004037.ref017]\].

Finally a last testable prediction of the *A~R~C* model is that the plot of *A* ~*P*~/*A* ~*R*~ as a function of *M*, and hence of the light fluence rate *I*, should be S-shaped with two asymptotes: i) when the light intensity is low, gravity dominates and *A* ~*R*~ = 0 and ii) when the light intensity is high. Moreover the equilibrium *A* ~*R*~ is shifted toward the orientation of incident light such that *A* ~*R*~ = *A* ~*P*~ ([Fig. 4](#pcbi.1004037.g004){ref-type="fig"}). These predictions are also consistent with the experimental data shown in [Fig. 2.D](#pcbi.1004037.g002){ref-type="fig"}.

Now that experimental agreement has been found with features of the photogravitropic models, it is interesting to discuss more in details the properties of the models and the insights it provides on the photo- and gravi-tropic control. This will be conducted through considering variant models of increasing complexity.

Insights on phototropism brought by the models: Photoception and proprioception {#sec003b}
-------------------------------------------------------------------------------

### Apical photoception: *A^a^C* model {#sec003ba}

When photoception is only apical (the *A^a^C* model; [equation 7](#pcbi.1004037.e017){ref-type="disp-formula"}), the curvature is constant along the organ (no dependency on *s*), so the shape of the organ is an arc of a circle ([Fig. 5.A](#pcbi.1004037.g005){ref-type="fig"}). This can be understood as the organ trying to bring the apical part towards the vertical with an identical driving signal everywhere. Finally, at the steady state, the curvature is also constant along the organ so that the orientation *A*(*s*) changes linearly from the base to the tip, where it reaches the vertical ([Fig. 5.B](#pcbi.1004037.g005){ref-type="fig"}). The solution of [equation 7](#pcbi.1004037.e017){ref-type="disp-formula"} depends only on the photoproprioceptive number *D* and on the propriosensitivity *γ*, but in a contrasted way.

![A. Straightening dynamics of the *A~a~C* model (*D* = 4).\
At each time the shape is an arc of circle. B. Steady-state shape of the *A^a^C* model. The shape is an arc of circle. As photoception dominates over proprioception, *D* increases (from blue to red); the apical angle reorients in the direction of the light field vector.](pcbi.1004037.g005){#pcbi.1004037.g005}

The apical angle of the steady state of the *A^a^* model is given by $$A\left( L,t\rightarrow\infty \right) = A_{0}\frac{1}{1 + D}$$ and the curvature of the steady state is simply defined by $$C\left( t\rightarrow\infty \right) = A_{0}L^{- 1}\frac{1}{1 + D^{- 1}}.$$

The apical angle and the curvature at the steady state are thus direct indicators of the ratio *D* between photoception and proprioception. When the photoceptive term dominates over the proprioceptive term (*D* → ∞), the apical angle tends to approach the direction of the light ([Fig. 5.B](#pcbi.1004037.g005){ref-type="fig"}). In sharp contrast to what was found for gravitropism \[[@pcbi.1004037.ref002]\], when apical perception is involved, proprioception is not necessary to reach a steady state. If *γ* = 0, the movement stops when the apical part is parallel to the direction of the light. Consequently (and concurrently), all parts of the organ near the base also cease to curve. Proprioception only modifies the maximal curvature that the system can reach. When proprioception increases (i.e., *D* decreases) the apical angle at steady state no longer reaches the vertical. In any case, however, the apical perception shapes the organ as an arc of a circle.

Considering now the transient dynamics before convergence to the steady state, the characteristic time at which the apical part of the organ reaches the vertical is given by $$T_{c} = \gamma^{- 1}\frac{1}{1 + D} = \frac{1}{\gamma + \nu}$$ Note that this characteristic time depends on both the proprioceptive and the photoceptive terms. This differs from the case in the *AC* model \[[@pcbi.1004037.ref002]\], in which, due to the purely local control, the characteristic time depends only on the proprioceptive term. The inclusion of a non-local photoceptive term thus makes the organ converge to its steady state faster than in the case of purely local perception. Also in contrast to the case of gravitropic movement \[[@pcbi.1004037.ref002]\], the dynamics of phototropic movement, governed by photoception and proprioception, cannot not be fully described by a unique dimensionless number, such as in this case *D*.

### Local photoception: Back to the *AC* model {#sec003bb}

In the case in which light perception is local, the photoproprioceptive model takes the form of an *AC* model ([equation 11](#pcbi.1004037.e021){ref-type="disp-formula"}). The dynamics and the steady- state shape are thus fully described by the dimensionless number, in this case *D*, alone. The salient features of the movement can then be summarized as follows: i) Initially the whole part of the plant that responds to the stimuli starts curving; ii) then, the curvature concentrates near the base while the apical part straightens; and iii) the number of oscillations around the vertical before convergence to the steady state increases with *D*. The dimensionless number *D* also reflects the ratio between the length of the organ and the length of the curved zone at steady state, in a strict analogy with what was found for the dimensionless number *B* in the *AC* model for gravitropism \[[@pcbi.1004037.ref002]\]. And in this case proprioception is necessary to avoid infinite oscillations and reach a steady state.

Having established a clear understanding of the behavior under the sole influence of either phototropism or gravitropism \[[@pcbi.1004037.ref002]\], we can now discuss the interaction between phototropism and gravitropism.

Insights on photo-gravitropic interaction brought by the models {#sec003c}
---------------------------------------------------------------

In order to reach insights on the dynamic of the interactions between photo -, gravi- (and proprio-)-ceptions on the control of the tropic movement, we will consider the case of a tropic motion at constant incident light intensity *I*, therefore studying the *A~R~C* model (the influence of *I* through changes in *M*, the ratio between photoception and graviception, having been clarifyied previously).

### Local photoception: the *A~R~C* model {#sec003ca}

The dynamics of the *A~R~C* model ([equation 15](#pcbi.1004037.e025){ref-type="disp-formula"}), where graviception and photoception both occur locally, is strictly equivalent to the dynamics of the graviproprioceptive*AC* model \[[@pcbi.1004037.ref002]\]. It is thus controlled by a unique dimensionless number. This dimensionless number drives the shape of the steady state (for example, the length *L* ~*c*~ over which the organ is curved at steady state), as well as the number of transient oscillations before the steady state is reached. But there are some differences between the two models: the dimensionless number, denoted *B* ^′^, that drives the full dynamics of the *A~R~C* model is a composition of the graviproprioceptive number $B = \frac{\beta L}{\gamma}$ and of the photoproprioceptive number *D* $$B^{\prime} = \frac{\left( \beta + \nu \right)L}{\gamma} = B + D$$ As *D* ≥ 0, the dimensionless number *B* ^′^ is larger than the graviproprioceptive number *B*, *B* ^′^ ≥ *B*. The convergence length *L* ~*c*~ is then modified $$L_{c} = \frac{L}{B^{\prime}}$$ This involves quantitative differences in the dynamics of the *A~R~C* model compared to that of the graviproprioceptive *AC* model. As *B* ^′^ is larger than *B* (*D* ≥ 0), the tropic movement driven by the local photograviproprioceptive *A~R~C* model ([equation 12](#pcbi.1004037.e022){ref-type="disp-formula"}) exhibits more transient oscillations than does the graviproprioceptive tropic movement driven by the *AC* model (for the same *B*), as well a smaller convergence length *L* ~*c*~.

In this case we can demonstrate that the tip is converging to the resultant photogravitropic angle *A* ~*R*~. Indeed, according to [equation 15](#pcbi.1004037.e025){ref-type="disp-formula"}, the apical angle of the steady state of the *A~R~C* model is given by $$\left. A^{\prime}(L,\mspace{180mu} t\rightarrow\infty)\mspace{180mu} = \mspace{180mu}{A^{\prime}}_{0}e^{- B^{\prime}} = {A^{\prime}}_{0}e^{- (B + D)} \right.$$ $$\left. A(L,\mspace{180mu} t\rightarrow\infty)\mspace{180mu} = \mspace{180mu}(A_{0} - A_{R})e^{- (B + D)} + A_{R} \right.$$ In the absence of light, *D* = 0, if the apical part of the organ reaches the vertical, *A*(*L*, *t* → ∞) ∼ 0, it is expected that *e* ^−*B*^ ∼ 0. As the length of convergence is smaller when phototropism and gravitropism interact than in the gravitropic case, *B*+*D*\>*B*, the term *e* ^−(*B*\ +\ *D*)^ can therefore be neglected $$A\left( L,t\rightarrow\infty \right) \sim A_{R}$$

Just as for *B* in the gravitropism case \[[@pcbi.1004037.ref002]\], it then becomes possible to design a simple morphometric estimate of *M* (the ratio between photoception and graviception at the current incident light intensity) for a plant organ with local and distributed gravisensitivity and photosensitivity, such as the *Arabidopsis* hypocotyl. The measurement of the apical angle of the steady state provides a simple way to directly measure *A* ~*R*~ and the ratio between photoception and graviception *M* can then be readily estimated using [equation 27](#pcbi.1004037.e055){ref-type="disp-formula"}.

### Apical photoception: $A_{R}^{a}C$ model {#sec003cb}

The $A_{R}^{a}C$ model ([equation 19](#pcbi.1004037.e031){ref-type="disp-formula"}), in which photoception is apical and graviception is local (as in the case of grass coleoptiles), is now explored. The steady-state shape of the $A_{R}^{a}C$ model is given by [equation 21](#pcbi.1004037.e033){ref-type="disp-formula"}. The left term inside the brackets of the steady-state shape defined by *A* ^′^(*s*, *t* → ∞) ([equation 21](#pcbi.1004037.e033){ref-type="disp-formula"}) is the steady-state shape of the graviproprioceptive solution, the *AC* model, alone $$A_{0}^{\prime}e^{- Bs/L}$$ whereas the right term can be considered as a corrective term resulting from the interaction between photoception and graviception. $$A_{0}^{\prime}e^{- B}\frac{1 - e^{- Bs/L}}{M + \left( {1 - e^{- B}} \right)}$$

This $A_{R}^{a}C$ model shares some properties with the *A~R~C* model. The perceptions of light and gravity act together to bring the system in the direction *A* ~*R*~, which results from the combination of the two different tropisms. But unlike the empirical equation for the photogravitropic equilibrium *PGEA* assumed by \[[@pcbi.1004037.ref017]\] ([equation 2](#pcbi.1004037.e002){ref-type="disp-formula"}), the angle *A* ~*R*~ defined through the [equation 26](#pcbi.1004037.e051){ref-type="disp-formula"} does not depend on the addition of the two stimuli but on the ratio between the two sensitivities, the photograviceptive number *M*. The length of the curved zone, *L* ~*c*~, depends only on the graviproprioceptive term *B*, $L_{c} = \frac{\gamma}{\beta} = \frac{L}{B}$. However, in contrast to the case in which graviception is dominated by proprioception, *L* ~*c*~ → ∞, the curvature is not necessarily equal to 0 (see [equation 22](#pcbi.1004037.e035){ref-type="disp-formula"}) and depends on the value of the photograviceptive number *M*.

Finally, it is still possible to design a simple morphometric estimate of *A* ~*R*~, of the ratio between photosensitivity and gravisensitivity *M* and of the ratio between gravisensitivity and propriosensitivity *B*. But this now requires the combination of two experiments, one in the light, and one in the dark. As in the *A~R~C* model, the measurement of the apical angle of the steady state in a photogravitropic experiment provides a simple way to directly measure *A* ~*R*~ (this will be further investigated when studying the limit cases of the $A_{R}^{a}C$ model in the next section). *B* can then be estimated independently using a gravitropic tilting experiment in the dark on the same organ (as in \[[@pcbi.1004037.ref002]\]). Finally the ratio between photosensitivity and gravisensitivity *M* can be calculated using the steady-state ([equation 21](#pcbi.1004037.e033){ref-type="disp-formula"}) at *s* = *L*.

### Limit cases and sensitivity analysis of the $A_{R}^{a}C$ model {#sec003cc}

The behavior of the $A_{R}^{a}C$ model is complex. To obtain a better understanding of it, in this section we explore the limit cases. The model's dynamics can be described using only *B*, the ratio between graviception and proprioception, and *M*, the ratio between photoception and graviception. Therefore only four different limit cases are to be considered: *B* → ∞, *B* → 0, *M* → ∞, *M* → 0. In what follows we characterize each case. Moreover, starting from each given limit case, we gather insights regarding the control dynamics by varying one of the dimensional numbers. This essentially constitutes an analysis of sensitivity of the output of the photo-gravi-proprioceptive model at steady state to changes in its dimensionless control parameters *B* and *M*.

i\) *B* → ∞. Graviception dominates proprioception. This case has already been shown to be a non-physiological case \[[@pcbi.1004037.ref002]\]. Indeed, in this case there is no steady state. The organ cannot converge to the vertical and oscillates indefinitely, displaying no tropism.

ii\) *B* → 0. Proprioception dominates graviception. *A* ~*R*~ is given by the direction of the light *A* ~*R*~ = *A* ~*P*~. The final shape is then given by $$A^{\prime}\left( s \right) = A_{0}^{\prime}\left( {1 - \frac{s}{1 + D^{- 1}}} \right)$$ This is the solution of the *A^a^C* model, [equation 3](#pcbi.1004037.e013){ref-type="disp-formula"}.

Starting from this limit case of pure photoproprioceptive equilibrium, it is now useful to consider the changes of the steady-states of the photograviproprioceptive model that occur when the graviproprioceptive number *B* is increased from 0 while the photoproprioceptive number *D* is maintained at a constant value ([Fig. 6.A](#pcbi.1004037.g006){ref-type="fig"}). This corresponds to increasing the gravisensitivity *β* relative to the photosensitivity *ν*, and hence to increasing the photograviceptive number *M*. As this change first affects the graviproprioceptive part of the $A_{R}^{a}C$ model, the changes in steady state of the graviproprioceptive *AC* model with increasing *B* are also plotted for comparison (dashed lines in [Fig. 6.A](#pcbi.1004037.g006){ref-type="fig"}). Starting from the arc-shaped solution of the photoproprioceptive *A^a^C* model, increasing *B* (and *M*) for a constant *D* brings the photo-gravi-proprioceptive model into steady states in which the apical part orients toward the changing *A* ~*R*~ (see insert in [Fig. 6.A](#pcbi.1004037.g006){ref-type="fig"}), whereas the curved zone becomes more and more concentrated at the base. Note also that the solution for the photo-gravi-proprioceptive $A_{R}^{a}C$ model is no longer equivalent to a rotation of the graviproprioceptive *AC* model toward *A* ~*R*~, but it is striking that the alignment of the apical part toward *A* ~*R*~ is maintained and actually increases with *B*.

![A. Steady-state shape of the $A_{R}^{a}C$ model for D = 4.\
The black line is the solution when *B* = 0. As *B* increases (solid colored line from red to blue), the orientation and the shape of the organ is modified. The dashed line shows the steady-state shape of the graviceptive equation. The expected orientation of the organ *A* ~*R*~ is shown at the bottom. B. Steady-state shape of the $A_{R}^{a}C$ model for *B* = 4. The black line is the solution when *D* = 0. As *D* increases (solid colored line from red to blue), the organ's orientation and shape are modified. The dashed line shows the steady-state shape of the photoceptive equation, *A^a^C* model. The expected orientation of the organ *A* ~*R*~ is shown at the bottom.](pcbi.1004037.g006){#pcbi.1004037.g006}

iii\) *M* → ∞. Graviception dominates photoception. For a given graviproprioceptive ratio *B*, this also means a very low photoproprioceptive number *D*, i.e., *D* → 0. The direction of the expected orientation is the direction of gravity, *A* ~*R*~ = 0. It follows easily that the steady state is given by: $$A\left( s \right) = A_{0}e^{- Bs/L}$$ This is the solution of the graviproprioceptive equation, the *AC* model, [equation 1](#pcbi.1004037.e001){ref-type="disp-formula"}.

We may now study the effect of increasing photosensitivity, i.e., increasing the photo-graviceptive number *D* and, by necessity, decreasing the photograviceptive number *M* ([Fig. 6.B](#pcbi.1004037.g006){ref-type="fig"}) for a constant value of *B*. As this change is first affecting the photoproprioceptive part of the model, the changes in the steady state of the photoproprioceptive *A^a^C* model with increasing *D* are also plotted for comparison (dashed lines in [Fig. 6.B](#pcbi.1004037.g006){ref-type="fig"}). When increasing *D* (and decreasing *M*) at constant *B*, changes in the steady-state shapes of the photo-gravi-proprioceptive $A_{R}^{a}C$ model are driven by changes in *A* ~*R*~ ([Fig. 6.B](#pcbi.1004037.g006){ref-type="fig"}). However, the solution never converges toward the photoproprioceptive solution of the *A* ^*a*^ *C* model. The orientation of the stem changes as *D* changes, but the length of the curved zone remains the same (as *B* is fixed). Even if graviception and photoception contribute equally towards orienting the apical part, they play different roles in the control of the global shape of the organ.

iv\) *M* → 0. Photoception dominates graviception. In this case, the expected orientation is driven simply by the direction of the light, *A* ~*R*~ = *A* ~*P*~. The steady-state shape is then given by $$A^{\prime}(s)\mspace{180mu} = \mspace{180mu} A_{0}\left( {e^{- \frac{Bs}{L}} - \frac{1 - e^{- Bs/L}}{e^{B} - 1}} \right)$$ $$C(s)\mspace{180mu} = \mspace{180mu}\frac{BA_{0}}{L}\frac{e^{- \frac{Bs}{L}}}{1 - e^{- B}}$$ It should be noted that, for any *B* that is not equal to zero, this solution is not the photo-proprioceptive *A^a^C* model but rather is a modification of the graviproprioceptive *AC* model (compare the plain line and the dashed lines in [Fig. 6.B](#pcbi.1004037.g006){ref-type="fig"}, especially for high values of *D*; green and blue-colored curves). The photoception sets the apical angle $$A^{\prime}\left( s \right) = 0$$ but the global shape is just a perturbation of the graviproprioceptive shape of the *AC* model. So even if photoception dominates over graviception, if the graviception does not tending towards *B* = 0, the distribution of the curvature in the steady state is still clearly determined by graviception and proprioception and thus by *B*.

The complete mechanism of control over the steady-state shape in the case of photogravitropism can now be summarized. The orientation of the apical part is determined by the ratio between graviception and photoception, the photograviceptive number *M*, whereas the length of the curved zone is determined only by the ratio between graviception and proprioception, the dimensionless number *B*. This can be illustrated clearly by considering the case in which the two tropisms act in the same direction (*A* ~*P*~ = 0, [Fig. 7.A](#pcbi.1004037.g007){ref-type="fig"}). As in previous examples, the initial condition is *D* = 0 (pure graviproprioceptive steady state) but with a lower value of *B* (so that proprioception prevents the tip from aligning in the direction of gravity). In these conditions, decreasing the value of *M* (by increasing *D* for a constant *B*) modifies the apical angle toward the vertical (phototropism reinforces gravitropism in this case), but the length of the curved zone remains unchanged (being determined solely by the graviproprioceptive process, and thus by the value of *B*).

![A. Steady-state shape of the $A_{R}^{a}C$ model for *B* = 2.\
The solid black line is the solution when *D* = 0; the dashed line is the solution when *D* → ∞ and *B* = 0. As *M* increases (from red to blue), the organ's orientation and shape are modified. B. Straightening dynamics of the $A_{R}^{a}C$ model (*B* = 2 *D* = 20). The solid black line is the solution when *D* = 0; the dashed line is the solution when *D* → ∞ and *B* = 0. During the movement, the organ reaches the steady state of the photoceptive equation but then goes to the steady state of the gravi-photo-proprioceptive equation. C. Space-time mapping of the curvature *C*(*s*, *t*) during straightening.](pcbi.1004037.g007){#pcbi.1004037.g007}

It is possible to use a similar analysis procedure to investigate the dynamics of movement within the $A_{R}^{a}C$ model, before convergence to the steady state (*i.e.* shape of the transients and time to converge to a steady state).. The first cases are straightforward.

i\) *B* → ∞ Graviception dominates proprioception. The organ never reaches the vertical \[[@pcbi.1004037.ref002]\].

ii\) *B* → 0 Proprioception dominates graviception. The dynamics are described by apical phototropism alone, *i.e.* the dynamics of the *A^a^C* model, as described in the corresponding section ([Fig. 5](#pcbi.1004037.g005){ref-type="fig"}).

iii\) *M* → ∞ Graviception dominates photoception. The dynamics of the $A_{R}^{a}C$ model are equivalent to the dynamics of the graviproprioceptive *AC* model as described in \[[@pcbi.1004037.ref002]\]. Briefly, for higher values of *B* (i.e., greater dominance of graviception over proprioception), the transients are more contorted, and the organ oscillates for a longer time around the steady state before converging.

iv\) *M* → 0 Photoception dominates graviception. This case must be discussed more carefully. In this case, the two processes, graviception and photoception, act as dominant drivers at two different time scales. In the photo-gravi-proprioceptive $A_{R}^{a}C$ model, photosensitivity has a faster effect than gravisensitivity does. This can be illustrated clearly by considering again the case in which the two tropisms act in the same direction (*A* ~*P*~ = 0, [Fig. 7.B](#pcbi.1004037.g007){ref-type="fig"}). The model predicts that the organ will bend faster towards the light in order to bring the apical part to the direction of the light, and to reach the steady-state shape of the photoproprioceptive equation. Then, the graviproprioceptive process takes over. The curvature concentrates near the base while the apical part remains oriented in the direction of *A* ~*P*~ ([Fig. 7.B and C](#pcbi.1004037.g007){ref-type="fig"}). This behavior can be compared to the experiments of the [Fig. 3](#pcbi.1004037.g003){ref-type="fig"}. First, the whole coleoptile is curving and the apical tip reaches an angle bigger the angle at steady state. Then the curvature concentrates near the base as the apical straightens. Similarities between the theoritical behavior and the experimental behavior are then observed. However due in part to the specific dynamics of the coleoptile and the inner leaf, it seems difficult to discuss further away the similarities and disparities between experiments and theory.

Discussion {#sec004}
==========

This study shows that a minimal modeling approach used to study shoot gravitropism, which gave rise to the *AC* model \[[@pcbi.1004037.ref002]\], can be readily extended to study the interactions among photoception, graviception and proprioception. The extended model, the *A~R~C* model, is minimal but complete. Simple dimensionless control parameters can be estimated from the *A~R~C* model, and the model is tractable and easy to understand. Furthermore, it is possible and straightforward to include the dependency of photosensitivity on light intensity *I*, thereby producing an augmented *A~R~C* model. Predictions from the latter model were assessed against detailed quantitative experimental results \[[@pcbi.1004037.ref017]\]. The model agreed very well with the data for small initial angles. All the experimental curves at various inclination angles and light intensities collapse to a single curve when using the *A~R~C* model with a power-law dependency between the photograviceptive control number *M* and the intensity of light *I*. Although the assessment is not comprehensive, the evidence presented supports the validity of the model. The hypotheses that the action of the tropic motor is fully driven by perception-regulation processes (H1), that the angles between the gravity and the light fields are first-order variables influencing the tropic dynamical movement (H2), that the two types of perceptions act additively (H4), and that the balance between gravisensing and photosensing depends on the perception of light intensity (H5) are upheld in the model.

The simplified versions and limiting cases of the *A~R~C* model served as useful tools for analyzing the dynamical interactions associated with perception during tropic movement, yielding five major insights.

i\) Despite having different set-points, gravitropism and phototropism act together to align the organ with the direction defined by *A* ~*R*~ ([equation 26](#pcbi.1004037.e051){ref-type="disp-formula"}). *A* ~*R*~ is determined by the dimensionless number *M* and hence by the ratio between gravisensitivity and photosensitivity. As in the *AC* model, the tropic movement is controlled globally through the different types of tropic perception, which together drive the local curving velocity. The final steady-state shape reflects the ratio of the respective sensitivities to light and gravity \[[@pcbi.1004037.ref002]\].

This may seem to be at odds with the phenomenological model for *PGEA* proposed by Galland, summarized in [equation 2](#pcbi.1004037.e002){ref-type="disp-formula"} \[[@pcbi.1004037.ref017]\], which states that the global steady-state shape at photo-gravitropic equilibrium depends additively on the phototropic and gravitropic stimuli, not on their ratio. However this phenomenological model for *PGEA* becomes problematic when considering the limiting cases. According to [equation 2](#pcbi.1004037.e002){ref-type="disp-formula"}, in darkness (i.e., when *I* = 0) the *PGEA* = −*g* sin(*A* ~0~). In other words, the gravitropic equilibrium orientation is a function of the initial angle. However, when there is no light, the set-point becomes simply the gravitropic set-point angle (*GSA*) \[[@pcbi.1004037.ref028]\], and the equilibrium shape is driven by the *GSA* and by the graviproprioceptive number *B* \[[@pcbi.1004037.ref002]\]. For ortho-gravitropic organs (i.e., in which the *GSA* is vertical) in plants in which graviception dominates proprioception (where the value of *B* is large), experimental observations confirm that the equilibrium tip angle is indeed vertical, that is, *PGEA* = 0 for any initial angle *A* ~0~ \[[@pcbi.1004037.ref002], [@pcbi.1004037.ref022]\]. The final gravitropic steady-state shape does not depend on sin(*A* ~0~) at all. Rather, the sine law dependency applies to the transient rate of changes in curvature, not to the final steady-state shape \[[@pcbi.1004037.ref001], [@pcbi.1004037.ref002], [@pcbi.1004037.ref004], [@pcbi.1004037.ref034]\]. In the other limiting case, when photoception dominates gravisensing (e.g., when the light is very bright, in microgravity or in agravitropic mutants), the set-point angle is determined by the direction of the light alone. The steady-state tip angle should therefore tend to *A* ~*P*~. The expression of the *A* ~*R*~ in [equation 26](#pcbi.1004037.e051){ref-type="disp-formula"} accommodates this condition, whereas the phenomenological model for *PGEA* in [equation 2](#pcbi.1004037.e002){ref-type="disp-formula"} predicts a tip angle that would depend on the light intensity with no saturation. The phenomenological model for *PGEA* \[[@pcbi.1004037.ref017]\] is thus not consistent for these two limiting cases.

The expression of the *A* ~*R*~ as the ratio between the sensitivities overcomes these problems directly and is also consistent with a complete dynamical model, the *A~R~C* model. *A* ~*R*~ is thus a better indicator than the previously defined *PGEA* \[[@pcbi.1004037.ref017]\]. In order to to unify the notation, we propose to follow \[[@pcbi.1004037.ref026]\], and to refer to this angle *A* ~*R*~ as the PhotoGravi Set Point Angle (*PGSA*).

ii\) The relation between the intensity of light (*I*) and the tropic response is better described by a power law, *ν* ∼ *I* ^−*b*^, than by a logarithmic law \[[@pcbi.1004037.ref017]\]. The power-law relationship was found to be more robust, fitting all the different experiments over more orders of magnitude of light intensity. For small angles, it was possible to collapse the different experimental results at various inclination angles and various light intensities into a single master curve. As the power law, the log law only relies on two parameters, however the log law also requires the assumption of a threshold in the perception of light. The power law is then more parsimonious. Further studies should be performed to investigate the photochemical basis of such an invariant scaling of the photo-gravisensitivity balance with light intensity, and to confirm the values of the exponent of the power law, estimated here as *a* ∼ 0.4.

iii\) The respective effects of local distributed perception along the organ and of apical perception are clearly defined for the first time. The *A^a^C* model, which deals with purely apical perception, predicts that the shape formed by the responding organ is an arc of a circle. It is reasoned that all the cells located along the organ receive the same secondary signal from the apical sensory apparatus, so every part of the organ curves to the same extent in order to bring the apical part into the direction of light. This prediction could be assessed experimentally on aerial organs but would require agravisensitive experimental conditions in which the $A_{R}^{a}(I)C$ model can be reduced to the *A^a^C* model. Experimental conditions that approximate the absence of gravisensitivity may be achieved using microgravity, clinostat experiments \[[@pcbi.1004037.ref032], [@pcbi.1004037.ref041]\], or mutant genotypes with impaired gravisensing \[[@pcbi.1004037.ref030], [@pcbi.1004037.ref041]\]. The relevant control parameter is the dimensionless photoceptive number *D*, which measures the ratio between photoception and proprioception. *D* can be easily estimated by measuring the apical angle of the steady state and the basal angle *A* ~0~, two variables routinely measured in most experiments on tropism.

iv\) Local perception (gravitropism) and apical perception (in species with apical phototropic sensing) interact so that aerial organs can forage for light without losing mechanical stability. Indeed, according to the $A_{R}^{a}(I)C$ model, the plant strives to align with the direction set by the ratio between photoception and graviception, but the curvature concentrates near the base while the apical part straightens. The length of the curved zone is under graviproprioceptive control only. This could be a functional adaption: Although the shoot must grow toward the light, the posture of the stem also needs to be controlled to ensure long-term stability. This may also explain how the apical part of some plants such as sunflowers can track the sun daily, while the plants continue to maintain control over their posture.

v\) Information from the *A~R~*(*I*)*C* model can be used to design methods for high-throughput phenotyping of the complete tropic control of organs by measuring the dimensionless parameters *B*, *D* and *M* of individual plants. *B* can be determined by measuring a plant's entire shape in simple gravitropic experiments in the darkness \[[@pcbi.1004037.ref002]\]. *D* can technically be measured in non-gravisensing conditions. However, it is simpler and faster to directly estimate *M* from the *PGSA*. Indeed, *M* is the ratio between gravisensitivity and photosensitivity, so the experiment can be conducted in gravisensing conditions. The apical angle at the steady state is the only measurement needed in order to estimate *M*. Once *B* and *M* are known, the value of *D* can be readily calculated, and the entire set of parameters controlling the dynamics of the system is then quantified. The dimensionless numbers *B*, *D* and *M* are real quantitative genetic traits that can be used to phenotype tropic mutants in genetically amenable plant models \[[@pcbi.1004037.ref030], [@pcbi.1004037.ref041]\] and subsequently identify the genetic and molecular mechanisms controlling tropic movements in plants.

The main limitation of the *A~R~*(*I*)*C* model in its present form is that it only deals with a geometrically simple light field, i.e., with uniform light and with the gravity field **g**, where the plant and the light field **l** lie within the same plane. It is possible, for example, to create spherical light fields by setting up point sources of light in proximity to the plant. The model presented here can be extended to include this case by making the light sensitivity parameter *ν* dependent on the spatial position *ν*(**x**, t). It is likely that the geometry of light fields in nature is even more diverse. For example, the direction of the sun's light changes over the course of the day. The direction of the light may also be outside the plane defined by the plants and the gravitational field. The modeling framework developed herein could be extended to deal with more complex light-field geometries, but this would require further mathematical and programming development. Such extensions would enable diverse movements to be simulated and would provide a more complete understanding of the ecological function of tropisms \[[@pcbi.1004037.ref005]\] beyond the core example of tropic control presented in this study.

It should also be noted that the use of etiolated organs (e.g., coleoptiles and hypocotyls) can complicate the analysis of phototropism. These organs are subject to photomorphogenetic effects that modify their behavior, independently of phototropic movements. Hypocotyl geometry is modified by the opening of the hook \[[@pcbi.1004037.ref042], [@pcbi.1004037.ref043]\]; in addition, it has been shown that pretreatment with light affects hypocotyl behavior \[[@pcbi.1004037.ref044]\]. Fully developed organs, which are less susceptible to photomorphogenetic effects, are therefore a powerful tool in the study of phototropism. Thus, kinematic analysis of the inflorescence of *Arabidopsis thaliana* might yield powerful insights regarding the phototropism process \[[@pcbi.1004037.ref019]\].

Finally and more generally, our modeling approach may also be used to study the tropic responses of other organisms, such as fungal stripes, sporangiophores or hyphae \[[@pcbi.1004037.ref018], [@pcbi.1004037.ref045]\], or of other plant organs such as roots. Plant roots might provide an interesting system in which to test the effects of apical perception and to assess the validity of the *A^a^C* model. As graviception in roots is purely apical \[[@pcbi.1004037.ref041]\], the *A^a^C* model can be readily extended to root gravitropism with the caveat that the extent of proprioception in roots remains unknown. This model could then be assessed by measuring the steady-state curvature (predicted to be constant along the organ) in hydroponics experiments to avoid confounding effects of interactions with substrates. If the model is validated, it may be possible to confirm a unified theory of the tropic movements of fixed organisms in natural conditions.

Supporting Information {#sec005}
======================

###### Supplementary information, Analytical results of the models and modeling aspects concerning the effect of the time of signal propagation along the organ.

Figure S1---Solution of the *A~a~C* model for instantaneous propagation, *T* ~*B*~ = 0. Figure S2---Simulation of the *A~a~C* model for different value of the ratio *T* ~*B*~/*T* ~*C*~. As long as the the propagation is faster than the characteristic time of the movement, the solution is similar to the solution with instantaneous propagation. Figure S3---Orientation of the organ in the *AaC* model for different value of the ratio *T* ~*B*~/*T* ~*C*~.

(PDF)

###### 

Click here for additional data file.

The authors thank Catherine Coutand for discussion and support. RB thanks the team of Herman Hofte at IJPB where model development was finalized. The English was edited by Karen Marron.

[^1]: The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: RB SD BM. Performed the experiments: RB. Analyzed the data: RB SD BM. Wrote the paper: RB BM.

[^3]: ⋆ These authors have contributed equally to the work
