# The self-organization of plant microtubules inside the cell volume yields their cortical localization, stable alignment, and sensitivity to external cues

Vincent Mirabet, Pawel Krupinski, Olivier O. Hamant, Elliot M. Meyerowitz, Henrik Jönsson, Arezki Boudaoud

► To cite this version:

Vincent Mirabet, Pawel Krupinski, Olivier O. Hamant, Elliot M. Meyerowitz, Henrik Jönsson, et al.. The self-organization of plant microtubules inside the cell volume yields their cortical localization, stable alignment, and sensitivity to external cues. PLoS Computational Biology, 2018, 14 (2), pp.e1006011. 10.1371/journal.pcbi.1006011. hal-02624527

HAL Id: hal-02624527

https://hal.inrae.fr/hal-02624527

Submitted on 26 May 2020

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

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

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

RESEARCH ARTICLE

# The self-organization of plant microtubules inside the cell volume yields their cortical localization, stable alignment, and sensitivity to external cues

Vincent Mirabet¹,²*, Pawel Krupinski³, Olivier Hamant¹, Elliot M. Meyerowitz²,⁴,⁵, Henrik Jönsson²,³,⁶, Arezki Boudaoud¹*

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

¹ Reproduction et Développement des Plantes, Univ. Lyon, ENS de Lyon, UCB Lyon 1, CNRS, INRA, F-69364 Lyon, France, ² Sainsbury Laboratory, University of Cambridge, Cambridge, United Kingdom, ³ Computational Biology and Biological Physics Group, Department of Theoretical Physics, Lund University, Lund, Sweden, ⁴ Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, California, United States of America, ⁵ Howard Hughes Medical Institute, California Institute of Technology, Pasadena, California, United States of America, ⁶ Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge, United Kingdom

* vincent.mirabet@gmail.com (VM); arezki.boudaoud@ens-lyon.fr (AB)

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

## OPEN ACCESS

**Citation:** Mirabet V, Krupinski P, Hamant O, Meyerowitz EM, Jönsson H, Boudaoud A (2018) The self-organization of plant microtubules inside the cell volume yields their cortical localization, stable alignment, and sensitivity to external cues. PLoS Comput Biol 14(2): e1006011. https://doi.org/10.1371/journal.pcbi.1006011

**Editor:** Eric N Cytrynbaum, University of British Columbia, CANADA

**Received:** July 25, 2017

**Accepted:** January 31, 2018

**Published:** February 20, 2018

**Copyright:** © 2018 Mirabet et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by Agence Nationale de la Recherche (http://www.agence-nationale-recherche.fr, grants ANR-10-BLAN-1516 to OH and 12-BSV2-0023-02 to AB), by the European Research Council (https://erc.europa.eu, grants #307387 to AB and #615739 to OH), and by the Gatsby Charitable Foundation (http://www.

## Abstract

Many cell functions rely on the ability of microtubules to self-organize as complex networks. In plants, cortical microtubules are essential to determine cell shape as they guide the deposition of cellulose microfibrils, and thus control mechanical anisotropy of the cell wall. Here we analyze how, in turn, cell shape may influence microtubule behavior. Building upon previous models that confined microtubules to the cell surface, we introduce an agent model of microtubules enclosed in a three-dimensional volume. We show that the microtubule network has spontaneous aligned configurations that could explain many experimental observations without resorting to specific regulation. In particular, we find that the preferred cortical localization of microtubules emerges from directional persistence of the microtubules, and their interactions with each other and with the stiff wall. We also identify microtubule parameters that seem relatively insensitive to cell shape, such as length or number. In contrast, microtubule array anisotropy depends on local curvature of the cell surface and global orientation follows robustly the longest axis of the cell. Lastly, we find that geometric cues may be overcome, as the network is capable of reorienting toward weak external directional cues. Altogether our simulations show that the microtubule network is a good transducer of weak external polarity, while at the same time, easily reaching stable global configurations.

## Author summary

Plants exhibit an astonishing diversity in architecture and morphology. A key to such diversity is the ability of their cells to coordinate and grow to reach a broad spectrum of
gatsby.org.uk, grants GAT3395/DAA to EMM and GAT3395-PR4 to HJ). EMM also acknowledges support from the Howard Hughes Medical Institute. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

sizes and shapes. Cell growth in plants is guided by the microtubule cytoskeleton. Here, we seek to understand how microtubules self-organize close to the cell surface. We build upon previous two-dimensional models and we consider microtubules as lines growing in three dimensions, accounting for interactions between microtubules or between microtubules and the cell surface. We show that microtubule arrays are able to adapt to various cell shapes and to reorient in response to external signals. Altogether, our results help to understand how the microtubule cytoskeleton contributes to the diversity of plant shapes and to how these shapes adapt to environment.

## Introduction

Despite their amazing diversity in shapes, biological organisms share some common structural components at the cellular level. Among those, one of the best conserved proteins across eukaryotes, tubulin, assembles into protofilaments, which in turn form 25 nm nanotubes known as microtubules, usually made of 13 protofilaments. The network of microtubules is highly labile and can reshape itself in a matter of minutes. In plants, microtubules form superstructures before (the preprophase band), during (the spindle) and after (the phragmoplast) cell division. Plant microtubules also form dense and organized arrays at the periphery of the cell during interphase [1] and these arrays are known as cortical microtubules (CMTs). The behavior of CMTs has been studied extensively at the molecular level [2]. One of their main functions is to guide the trajectory of the transmembrane cellulose synthase complex and thus to bias the orientation of cellulose microfibrils in the wall. This in turn impacts the mechanical anisotropy of the cell wall and controls growth direction [3–6]. This function explains why most mutants impaired in microtubule-associated proteins exhibit strong morphological defects [7].

Whereas this provides a clear picture on how microtubules impact cell shape, in turn how cell shape impacts microtubule behavior has been less explored. There is evidence that microtubule orientation depends on cell shape [8–10], with microtubules being mostly transverse to the longest axis, but this might require specific regulation because microtubules orient along the longer axis of the confining domain in vitro [11]. There is also evidence that shape-derived mechanical stress can bias cellulose deposition, possibly through microtubule orientation towards the direction of maximal tension, both at the tissue and single cell scales [8, 12–21]. The molecular mechanism behind remains largely unknown. Finally, cortical microtubules orientation may change in response to signals such as blue light or hormones, see for instance [17, 22] and may be oriented by the hydrodynamic forces due to cytoplasmic streaming [23, 24]. Here we use modeling to explore the relative contributions of cell geometry and external directional cues in the final microtubule organization.

The molecular basis for microtubule dynamics is rather well established. Consistent with the absence of centrosome in land plants, microtubule nucleation is dispersed in plant cells, as it occurs at the cell cortex [25], along existing microtubules during branching events [25–27], and at the nuclear envelope [28]. As they grow, microtubules form stiff and polar structures. They can alternate growth, pause and shrinking at the so-called plus end [29], whereas they mainly shrink or pause at the minus end [30]. The combination of an average shrinkage at the minus end and dynamic instability at the plus end leads to an overall displacement of the microtubule, also called hybrid treadmilling [30], with a dominant contribution of short treadmilling microtubules in the final microtubule organization [31]. The growth of microtubules in persistent directions is the main cause for microtubule encounters.
When one microtubule encounters another microtubule, different outcomes can be observed [9, 32]: if the encounter angle is shallow, zippering can occur, i.e. the growing microtubule bends and continues its polymerization along the encountered microtubule, which leads to the creation of microtubule bundles; if the encounter angle is steep, crossover can occur, i.e. a microtubule polymerizing without deviating its trajectory and crossing over the encountered microtubule; or alternatively catastrophe is triggered, i.e. a rapid plus end shrinkage after contact with the encountered microtubule. Such selective pruning of microtubules may explain how microtubules can form parallel arrays from initially random orientations, and conversely change the net orientation of their arrays over time, through a phase of randomisation [33, 34]. Selective pruning has indeed an essential ingredient of most models for microtubule dynamics [35–44]. The presence or absence of different microtubules associated proteins (MAPs) can modulate the stability of microtubules or their capacity to form bundles and to self-organize. For instance, the microtubule severing protein Katanin accounts for most of the pruning events at crossover sites [45].

The microtubule network is a typical example of a self-organizing system, where properties of individual elements and their interactions induce specific and sometimes counter-intuitive global properties. To predict how regulation at the level of each microtubule can give rise to specific global outcomes, one can resort to computational models. Modeling approaches have been developed, simplifying microtubule interactions by restricting them to the plasma membrane, i.e. a simpler 2D space [40, 46]. In those agent-based models, several microtubule properties were coded and interactions between CMTs, based on these properties, were simulated. The outcome is an emergent network, whose characteristics can be analyzed. For instance, increased microtubule severing was predicted to generate a larger number of free microtubules, more amenable to bundle into aligned arrays [9, 42] and this was observed in experiments [47].

So far, most of the microtubule models have been implemented in a 2D space or with microtubules confined to the surface of the cell. A major outcome of such models was to demonstrate that global orientations of the network can spontaneously emerge from the interactions between microtubules [48]. Many combinations of parameters and behaviors have been studied: instability at the plus end [35, 37], role of zippering [9, 36, 38, 39], nucleation modes [39, 42], and severing [44]. Beyond their differences, a global orientation emerges in most of these combinations suggesting that converging toward a global orientation is a robust feature of microtubule networks. Conversely, the diverse combinations of microtubule properties provide different scenarios for the fine-tuning of the network structure and stability of this emerging behavior.

Some aspects of cell geometry were related to microtubule behavior in certain simulations. Simulations showed how different directional biases in nucleation can induce an ordering of the array toward directions that are correlated to cell geometry [35, 43]. Further, branching nucleation rules can elicit handedness of the global direction of microtubule arrays, provided that the branching is biased toward one direction [41]. Other studies used a simulation space where borders, analog to cell edges, induce more or fewer catastrophe events or are more or less permissive toward microtubule growth [10, 33, 35, 41]. Most studies concluded that a global orientation of microtubules can be correlated to cell face orientations.

The contribution of the third dimension to microtubule behavior has started to be investigated. Computational models for animal systems have focused on 3D considerations but the nucleation hypotheses are too different from that in plants to be transposed directly [49, 50]. Fully 3D models suited for plants are still lacking: almost all existing studies have confined microtubules to surfaces embedded in 3D [10, 11, 41, 51]. In [35], a 2D model was extended into a full 3D model but it did not include cell boundaries, which yielded microtubules
distributed over the whole simulated domain, in contrast with the cortical localization of
microtubules in planta.

In this paper, we explore the influence of 3D cell shape on the basic properties of a dynamic microtubule network. We do not assume that microtubules are confined to the cell surface; rather, we simulate a closed volume where microtubules are more or less free to grow in all directions. Anchoring to the membrane is not imposed by the model and instead becomes a variable in the model. Using this framework, we investigate to what extent microtubule interaction with the membrane can influence microtubule dynamics. Our study also addresses the relative contributions of cell shape, microtubule interactions, and external directional cues in network organization.

## Results

### A model for interphase microtubules growing at the plasma membrane or within the cytoplasm

Following previous studies, we modelled microtubules as a set of line segments that nucleate, grow, shrink, and interact with each other and with the cell surface represented as a triangular mesh (see Methods for details). Nucleation of the minus end occurs randomly at the surface. Growth occurs from the plus end with a small directional noise that is related to the persistence length of microtubules (Fig 1A). Shrinkage starts randomly at the nucleation site (minus end) and then continues at constant velocity.

A microtubule that encounters a preexisting microtubule either changes direction to that of the preexisting microtubule if the encounter angle is shallow (a process known as zippering, see Fig 1A), and otherwise starts shrinkage from the plus end (“head-on” collision). We considered two types of interaction with the cell surface: strong anchoring, whereby microtubules remain on the surface, as in previous studies [10, 11, 41, 51], and weak anchoring, whereby microtubules are prevented from leaving the cell interior. More specifically, in the case of weak anchoring, the interaction between a growing microtubule and the nearby surface is similar to the interaction between two microtubules: the microtubule encountering the surface at a steep angle starts shrinking, and otherwise starts growing tangentially to the surface (Fig 1A); a microtubule may leave the surface because of the directional noise.

We considered three base shapes: cube, “square” (flattened cube), and “long” (elongated small cube), of dimensions in the order of 10 µm, typical of plant cells. These shapes were smoothed so that the maximal curvature corresponded to a radius of either ~1 µm (“sharp”) or ~5 µm (“smooth”), which corresponds to typical radii of curvature measured in root epidermis [10], see Fig 1B. We also considered an ellipsoidal shape when investigating the effects of an external cue. Although these shapes are not fully realistic, they make it easier to disentangle the geometrical parameters influencing microtubule dynamics.

A typical simulation with weak anchoring is given in S1 Video. Corresponding snapshots are shown in Fig 1C and 1D with various 3D and 2D views. A few first observations can be made: the microtubules tend to bundle; a well-defined local orientation appears; microtubules appear to be mostly close to the cortex.

### Microtubules become cortical in a 3D space because of their directional persistence and growth mode

We first considered the effect of the anchoring of the microtubules to the membrane. There are proteins that have been shown to be associated with both microtubules and a plasma membrane component in plants [52]. For instance, CELLULOSE SYNTHASE INTERACTIVE
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Shapes of simulated cells. (A) Top left: Two time steps of a simulation when the microtubule grows freely inside the cell and simultaneously shrinks from the minus end; at the plus side (arrowhead), the new vector is calculated by adding a small random deviation to the previous vector. Bottom left: Example of a zippering process—When the new microtubule encounters an existing one at a shallow angle (if the angle $a$ is smaller than threshold $\alpha$), the new direction of the vectors follows that of the leading microtubule; a steep angle would lead to plus end shrinkage. Top right: Strong anchoring to the plasma membrane—</figcaption></figure>
The microtubule grows tangentially to the local plane. Bottom right: Weak anchoring to the membrane—When the microtubule encounters the membrane at a shallow angle (a < α), it continues tangentially until the random deviation leads it to leave the membrane. (B) Shapes of the envelope of the cell: cube, “square”, and “elongated”. First row: smooth shapes. Second row: sharp shapes. Scales are not respected. (C,D) Example of a simulation in a sharp square with weak anchoring and default parameters (see Methods); the length of the simulation is 50000 timesteps, corresponding to approximately 100 minutes, also see S2 Video. (C) Snapshot showing all microtubules (left), microtubules located in a central slice that is ~2.4µm-thick (center), and the projection of this slice in 2D (right). (D) Snapshots taken after 10000 timesteps intervals (~20 min); the pictures show a confocal microscopy-like z-projection (see Methods).

https://doi.org/10.1371/journal.pcbi.1006011.g001

PROTEIN 1 (CSI1) interacts with both CMTs and the cellulose synthase (CESA) complex [6, 53] and CSI1 was also proposed to stabilize the microtubule network [54]. Yet, the influence of CMT-CESA interactions on the microtubule network is still poorly understood. Thus, we took advantage of the 3D nature of our model to study the impact of the anchoring rule to the membrane on the global parameters of the microtubule network.

We investigated the microtubule dynamics when the anchoring to the membrane is weak. As a reference case, we also considered strong anchoring, whereby microtubule are constrained to grow on the membrane, as if putative anchoring proteins were highly concentrated. As expected, strong interactions led to a surface-localized cortical zone with microtubules trajectories embedded in the plane parallel to the mesh (Fig 2A).

In the case of weak anchoring, microtubules grow in all directions, but occasionally, as they encounter the membrane, their direction may be transiently tangent to the membrane. The typical distance between microtubules and membrane ranges from 50 to 250 nm (Fig 2A) according to the shape. Even if weak anchoring allows microtubules to grow through all the cell volume, we find that such weak interaction with the membrane is enough to elicit the existence of cortical microtubules. Therefore, the three-dimensional nature of our model helps us demonstrate that strong anchoring is not required for the presence of large populations of cortical microtubules in plant cells: the directional persistence, together microtubule growth mode, can cause such sub-cellular localization.

As the microtubules tend to stay at close to the membrane, they also bundle, with the proportion of tubulin in bundles varying from 30% to 75% (Fig 2B). The ability for the microtubule network to generate a spontaneous bundled structure is consistent with previous models constraining microtubules to the surface. Strikingly, this effect is also present in the case of weak anchoring.

## Strong anchoring decreases microtubule number and length, and increases microtubule bundling

Independently of the encounter rule, weak anchoring strength increases the total number and the size of microtubules (by about 20%, in length and in number) when compared to strong anchoring. Weak anchoring also yields less bundling (reduction of about 20%). A likely explanation is that a weak anchoring to the membrane allows microtubules to “escape” inside the cell, thus diminishing the encounter probability. Consequently, microtubules weakly bound to the membrane have more space to grow and are less subject to shrinkage-induced collisions or bundling with other microtubules (Fig 2B).

## The number and length of microtubules, as well as the proportion of microtubules in bundles, are relatively independent of cell shape

Next, we used our model to determine the consequences of changing cell shape on global properties of the network. We simulated the network in three main sharp shapes represented
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Influence of the anchoring to the membrane. (A) Mean distance of microtubules to the membrane according to cell shape and to weak/strong anchoring to the membrane. (B) Effect of anchoring on three properties of the microtubule network (data for all shapes pooled together). Left: total number of microtubules. Middle: mean microtubule length. Right: proportion of tubulin in bundles. Model parameters have the default values (see Methods).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006011.g002

on Fig 1A. The results from the simulations indicate that these shapes only have a small effect on the number of microtubules (less than 10% difference), on the length of microtubules (5 to 15% difference) and on the proportion of bundles (less than 15% difference). Overall, elon-gated cells have more microtubules, more bundles, and longer microtubules, while cubes show the lower values (Fig 3). A likely explanation is that microtubules are more likely to follow the long axis in the long shape (see below) so that they are less affected by cell boundaries.
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Influence of global shape on three properties of the network. Data for sharp shapes: cube, square, and long (Fig 1). From left to right: number of microtubules, mean microtubule length, and proportion of tubulin in bundles. Top row: weak anchoring; bottom row: strong anchoring. Model parameters have the default values (see Methods).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006011.g003

Microtubule array anisotropy is influenced by cell curvature but not by global shape

We also analyzed the effect of cell shape on the microtubule array anisotropy, averaged over the cell (see Methods); this is quantified with an order parameter with values between 0 and 1.
As microtubule array anisotropy is skewed towards low values (found inside the cell) we used a non-parametric test based on ranks for statistical comparisons. Consistently with its effect on bundling (see above), we found that anchoring strength slightly affects the anisotropy of the microtubule arrays. Weak anchoring decreases the anisotropy of the network by about 10% (Fig 4C).

<figure><img src="image_8.png" /><figcaption>**Fig 4.** Anisotropy of the network. (A) Effect of global shape: Anisotropy for the three shapes (sharp/smooth and weak/strong anchoring were pooled together). (B) Effect of curvature: Anisotropy for sharp shapes vs. smooth shapes with strong anchoring. (C) Effect of curvature: Anisotropy for sharp shapes vs. smooth shapes with strong anchoring. (D) Effect of anchoring strength: Anisotropy for all shapes (pooled). In all panels, the data are pooled except for the parameter that is varied in the subfigure; model parameters have the default values (see Methods).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006011.g004
The type of cell shape did not appear to influence the global anisotropy of the microtubule network (Fig 4A). This is interesting as it suggests that different cell types with various shapes do not require differential regulation of the network in order to maintain the anisotropic properties of their CMT arrays.

However, smooth and sharp shapes differ significantly in anisotropy of the microtubule network (by 10%, *p* < 0.001); more curved cell edges in sharp shapes correspond to a lower anisotropy (Fig 4B). Larger surface curvature induce more “heads-on” collisions when microtubules grow nearby, leading to a smaller microtubule density, as they become more distant from the surface (Fig 2A); this would lead to more spatial variations in the orientation of microtubules and hence lower anisotropy.

## The average microtubule orientation is strongly influenced by cell shape

In order to determine how the shape of the cell influences the global orientation of the network, we measured the distribution of orientations along two opposite faces of the cells. We chose two arbitrary faces for cube shape, the largest faces for square shape, and two of the largest faces for the long shape; accordingly, these two opposite faces are among the two largest faces in area for the shape of interest. An angle of 0° corresponds to the long axis in the case of the elongated cell. Angles of 90° or -90° correspond to directions perpendicular to that axis.

First, we observed that in the case of square shaped faces, most of the microtubules align along the cell face diagonal, i.e. the longest path (Fig 5). This occurred whatever the anchoring strength and the shape of faces at the side. Second, we observed a strong correlation between the axis of the cell and microtubule orientation, and this correlation is the highest for the elongated cell: The microtubule distribution is always either maximal or minimal at an angle of 0. This effect is higher in the case of strong anchoring, whereas in the case of weak anchoring, secondary peaks show that the diagonals are also overrepresented. These results indicate that the microtubule network is able to read the longest axis of the cell and orient toward that axis, by default. Interestingly, cortical microtubules become longitudinal in hypocotyl cells when growth stops [5, 55, 56], suggesting that they may adopt their configuration by default in that situation.

## Small external directional bias is sufficient to affect the orientation of the whole network

Next, we investigated the robustness of microtubule network. Microtubule arrays entirely reorganize during cell division [57]. Light and hormones can also completely reorient the microtubule network within minutes [33, 34, 58], suggesting that the constraints on microtubules array must not be too strong to allow such rapid reorganization in vivo. We tested whether our model provides such adaptability, using the case of external, directional cues (Fig 6 and S2 Video).

We investigated the effect of an external cue, assuming that, as microtubules grow in the cell, their growth direction is biased toward the direction of cue, with a specific weight. We used an ellipsoidal cell with a circumferential cue, which is orthoradial to the long axis of the cell. Such a cue could, for instance, be related to cell polarity (vector along the long axis) or could be a proxy of mechanical tension [19]. In the absence of the cue, microtubules are on average parallel to the main axis of the cell. Strikingly, our simulations indicate that even a very low bias (0.1%) could disrupt the main orientation of the network. When the weight reaches 1% or more, microtubules massively reorient toward the direction of the cue. Interestingly, the transition between the longitudinal and circumferential orientation occurs when the bias is comparable to the directional noise that is used to account for the persistence length of
<figure><img src="image_9.png" /><figcaption>**Fig 5. Orientation of the network.** Distribution of the orientation of unit vectors (local orientation of microtubules) located close to the top and bottom faces. Top and bottom faces refer to the largest faces of the square shape represented in C, and to arbitrary opposite large faces of the cube or the long shapes. The reference angle 0 corresponds to the longest axis of the long cell and is parallel to one of the edges of the top face in the two other cases. (A,B) Distribution of orientations for sharp cube, sharp square, and sharp long shapes, as labeled in Fig 1. (A) Weak</figcaption></figure>
anchoring and (B) strong anchoring. (C) Arrow showing the reference for angle measurement for the long shape (top) and for the square and cube shapes (bottom). The shapes are not at the same scale. Model parameters have the default values (see Methods).

https://doi.org/10.1371/journal.pcbi.1006011.g005

microtubules. This reinforces the idea that fluctuations of microtubules and the interactions between them lead to alignment with the long axis of the cell and that the small cue is sufficient to counteract these fluctuations and influence the orientation of microtubules. Accordingly, we found that anisotropy increases from no cue to a weight of about 1% and then seems to saturate (S1 Fig).

These results indicate that, despite an apparent robust organization, the microtubule network remains extremely sensitive to directional cues. As such, it is capable of reading slight directional cues and generating an ordered array in that direction. This ability to read directional cues is probably linked to the self-enhancement of microtubule orientation through their interactions: As more microtubules orient toward a direction, they prevent growth of microtubules in the perpendicular direction.

Our simulations indicate that the network should exhibit two behaviors concerning its orientation. When no external directional cue is present, the network orients toward the main

<figure><img src="image_10.png" /><figcaption>**Fig 6.** Influence of directional cues on the global orientation of the network. Simulations in an ellipsoidal cell with a circumferential cue. (A) Histogram of the scalar product between the unit vectors (defining microtubules locally) and the long axis of the ellipsoid. Black: no directional cue. Colors: directional cue with influence on the direction. The strength, $b_d$, of the signal is: 0.1% (yellow), 0.2% (green), 1% (red), and 2% (blue). (B) Simulation result for a weight of 2%, also see S2 Video. (C) Simulation result with no directional cue. Model parameters: $n_p = 2.4 \cdot 10^{-7}$, $n_s = 10^{-3}$, $\alpha = 40^\circ$, weak anchoring.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006011.g006
axis of the cell, and generates a polarity that is a direct reading of the global shape. When a
directional cue is present, the network reorients so as to emphasize the direction of the cue.

Discussion

A model for microtubules growing in the cell volume

The impact of external cues on the microtubule network has been extensively characterized in experiments. However, because real cells and tissues are never really devoid of external cues, the behavior of the microtubule network by default in a plant cell remains an open question. Our work provides some clues to address this question, by taking into account the 3D shape of the cell and by using a minimal set of parameters for microtubule behavior. For instance, the model with weak anchoring does not require a specific rule for the crossing of microtubules (unlike when microtubules are confined to the cell surface), because microtubules naturally cross each other if they are farther than their diameter. We found that microtubule directional persistence largely determines the subcellular localization and orientation of the microtubule network in various cell shapes. We also identified parameters that seemed relatively insensitive to cell shape, such as microtubule length or number. Last, we found that microtubule dynamics yields at the same time stable orientation and high sensitivity to directional cues, even when such cues go against the default orientation. Altogether, this provides a conceptual framework to dissect the exact contribution of microtubule regulators to the microtubule network organization, in relation to 3D cell shape.

Cortical localization is an emerging property of microtubule directional persistence and dynamics in convex cells

Based on TEM images where microtubules are often seen very close to the plasma membrane, it is assumed that anchoring of microtubules to the plasma membrane is relatively strong. However it remains unknown how such anchoring would be mediated [52]. Many biochemi- cal studies have been performed to extract proteins that would link the plasma membrane to cortical microtubules, and so far the only published candidate is a phospholipase [59], for which no follow-up results have been obtained to the best of our knowledge. Other links have been put forward, such as CLASP [10] or CELLULOSE SYNTHASE INTERACTIVE PRO- TEIN 1 (CSI1) [53, 60], but they might represent rather indirect regulators of the link between microtubules and plasma membrane. Does this mean that microtubule could be cortical with- out any anchoring module? Our results suggest that in most of the measured cell shapes, microtubules do not need to be strongly connected to the membrane to remain cortical. This prediction implies that modulating anchoring would affect self-organizing properties tangen- tially to the cell surface, rather than modulating the density of microtubules inside the cell. This would allow molecular regulators to modify the microtubule organization without directly affecting the rate of cellulose deposition.

A high curvature increases anisotropy

In our simulations, we observed that anisotropy varied according to curvature of face-face con-
tacts of the cell, which is an emerging property of the weak anchoring mode. Cells that have
sharper edges exhibit decreased array anisotropy compared to smoother cells. This result is
interesting, as, for instance, epidermal cells possess edges with different curvatures [10]. Simi-
larly, in the L1 layer of the shoot apical meristem, the outermost wall has a stronger curvature
than all the walls separating the cell from its neighbors. Pavement cells also display a broad
range of curvature values. Based on our results, higher anisotropy would be produced in L1 cells without the need for specific regulation.

Such results are difficult to account for in models with microtubules constrained to surfaces embedded in 3D, in which for instance additional hypotheses linking curvature and catastrophe rate are implemented [10]. Nevertheless, our results do not preclude microtubule associated proteins from additionally modulating microtubule dynamics according to curvature or to other membrane-localized cues [52]; we only account for the default behavior of microtubule networks according to surface curvature.

Similarly, changes in the curvature of the epidermal wall can occur through changes in internal pressure, which in return influence the anisotropy of the network [61]. An increase in the microtubule organization could be a result of an increased pressure, that would increase the curvature of the epidermal cell. Further experimental work is required to investigate the role of curvature on microtubule behavior and its relation to mechanical stress in the epidermis.

## Influence of shape on the global directionality of the network

Our simulations show that the cell aspect ratio has an important impact on the global orientation of the network. The predicted default behavior of microtubules is their alignment parallel to the long axis of a cell, due to the directional persistence of microtubules. This is in agreement with previous models with microtubules confined to surfaces embedded in 3D [11, 51], where the default orientation is longitudinal for long cylinders.

This default state was observed with microtubules polymerizing in vitro inside elongated 3D chambers [11]. In slowly growing cells of the hypocotyl, microtubules are oriented along the long axis of the cell, whereas microtubules are circumferential in rapidly elongating cells [55, 56]. Our model suggests that directional cues are needed to avoid this default orientation is growing cells. At the boundary between the shoot apical meristem and the primordia, cell division leads to cell shapes that are elongated along the axis of the boundary; our model predicts that microtubules will be oriented along the same direction by default, amplifying their response to mechanical stress [19].

## A weak directional bias is sufficient to change the orientation of the microtubule network

In this study, we show that the microtubule network is oriented by default along the longest axis of the cell. However, microtubules in plants often show supracellular orientation, independently of cell shape, a behavior that has been ascribed to tissue-level signals, notably mechanical stress [18–20]. Moreover, It has been demonstrated that inside the cell, microtubules orientation is coupled to polarity markers such as proteins from the PIN FORMED and RHO OF PLANTS families [62, 63]. Simulations have assessed how localized membrane heterogeneity could result in a biased orientation of the microtubule network [10, 41]. In this study, we show that a weak directional cue influencing microtubule growth rapidly modifies the orientation of the network towards the direction of the cue. This cue could be due to mechanical stress, hormone gradients [22], or to cytoplasmic streaming [23, 24], for instance. As such, the network behaves as a sensor translating an external directional information into a structural polarity inside the cell. The coexistence of a default orientation and a strong ability to reorient could shed a new light on orientation changes in cells. Changes in microtubule orientation need not be always related to specific regulation but may also be related to the arrest of signals and the return to the default state. This could be occurring in the shift from transverse to longitudinal orientation in hypocotyls responding to light or hormones [33, 34, 58].
## Perspectives

The shape of the cell has little influence on mean length, number of microtubules or bundles proportions. In addition, anisotropy of the network is not highly correlated to changes in global cell shape. This prediction of a robust network suggests that plant cells do not need specific regulations to compensate for their great variations in shapes. Accordingly, the microtubule network appears as a good polarity system, with a default orientation and a high sensitivity to directional cues. It was recently shown that, for global polarity to emerge in a tissue, an important requirement is the existence of internal cellular polarity [64]. In this work we show that the microtubule network is suited for such a requirement.

Overall, microtubules and associated proteins form a complex self-organizing system that is difficult to comprehend without resorting to models. The results obtained here demonstrate that our three-dimensional model provides a framework to test hypotheses on the regulation of the microtubule cytoskeleton in plant cells.

The model given here is only the beginning to a more complete analysis. We have not yet incorporated microtubule severing [44, 58], microtubule branching [26, 27], nucleation at the nuclear envelope [28], or the possible effects of connections between cortical microtubules and cellulose fibrils outside of the cell as mediated by the cellulose synthase complex [65]. Severing in particular has been shown to be key to microtubule reorientation [44] following mechanical signals [66]. We also have not included limiting levels of tubulin [51], which could affect overall microtubule number. Altogether, we expect our model to help progress in understanding how microtubule self-organization integrates directional cues with three-dimensional cell shape and how microtubule-associated proteins modulate this integration.

## Methods

The dynamical microtubule model was implemented in C++. The code is available from https://gitlab.com/slcu/teamHJ/vincent/microtubule_simulations. Simulations were performed on Intel/AMD desktop computers running Debian and Ubuntu operating systems.

### The microtubule network

Microtubules were coded as 3D multi-segment vectors of constant length. A ring of tubulin of the length of a dimer is represented as a unit vector in the simulation. Microtubule growth in the model occurs by adding one vector element at the plus end of the microtubule, at the position of the end of the last vector. In plants, microtubules are considered to be mostly static, their growth and shrinkage are the result of treadmilling processes. To code for microtubule directional persistence (which relates physically to bending stiffness), the direction at which a new vector is added to an existing microtubule changes by a small random amount. Microtubule shrinkage occurs at the minus end by removing the first vector from the list. The main model parameters are shown in Table 1.

<table><caption>Table 1. Main parameters for microtubule dynamics.</caption><thead><tr><th>Parameter</th><th>Nucleation frequency</th><th>Shrinkage probability</th><th>Random direction</th><th>Encounter angle</th><th>Anchoring strength</th><th>Bias strength</th></tr></thead><tbody><tr><td>Label</td><td>n<sub>p</sub></td><td>n<sub>s</sub></td><td>r<sub>d</sub></td><td>α</td><td></td><td>b<sub>d</sub></td></tr><tr><td>Default value</td><td>2.4 · 10<sup>-7</sup>, 4.7 · 10<sup>-7</sup>, and 7.1 · 10<sup>-7</sup> per unit surface per time step</td><td>0.001 per nucleation site per time step</td><td>0.025</td><td>40°</td><td>strong or weak</td><td>0</td></tr></tbody></table>

https://doi.org/10.1371/journal.pcbi.1006011.t001
## Time and space units

Typically, a cell has a width of several micrometers, and we take the unit of length as 8nm, the height of a ring of tubulin. Considering a measured speed of growth at plus ends of 3-5 µm/min [67], a simulation time step is approximately 0.1 to 0.2 s. A typical simulation of 10000 time steps thus represents 15 to 30 min of real time.

## Nucleation and minus end behavior

The microtubules are nucleated on the cell surface [25–27] at a constant rate. The default value is $n_p = 4.7 \cdot 10^{-7}$ per time step per unit surface, corresponding to 1 to 2 nucleation events per time step, or 5 to 20 per second. Once nucleated, the microtubules do not immediately shrink. At each of the time steps that follow nucleation, a microtubule has a probability $n_s$ to begin shrinkage. Once a microtubule has started to shrink, one vector is removed from the minus end at each time step.

## Plus-end growth

**Default growth.** At each time step, a new vector is added to the plus end. In the absence of interactions with the membrane or with neighbouring microtubules, the new direction is that of the previous vector (tubulin), modified by a random direction. The new unit vector is computed as

$$ \mathbf{r}_{n+1} = \frac{(1 - r_d)\mathbf{r}_n + r_d\mathbf{u}}{\|(1 - r_d)\mathbf{r}_n + r_d\mathbf{u}\|}, \quad (1) $$

where $\mathbf{r}_n$ is the previous unit vector, $\mathbf{u}$ is a random unit vector, and $\|\|$ is the standard Euclidean norm.

At lower order in $r_d$ (which is a small parameter), the correlation between two consecutive vectors take the form $\langle \mathbf{r}_{n+1} \cdot \mathbf{r}_n \rangle = 1 - r_d^2/4$ and should be equal to $1 - 1/2\ell/p$ where $\ell$ is the length of one vector and $p$ is the persistence length [68]. Therefore, the default value of the random direction parameter $r_d = 0.025$ corresponds to a persistence length $p = 2\ell/r_d^2$ of 26 µm. This is in the range of values for persistence length measured in vivo, 20 µm [69] to 30 µm [70], which is much smaller than those found in vitro [71, 72], likely due to the action of proteins interacting with microtubules.

**Directional cue.** In the presence of directional cue, the new vector unit vector is computed as the weighted average of its default direction and of the direction of the cue (unit vector $\mathbf{b}$) weighted by $b_d$:

$$ \mathbf{r}_{n+1} = \frac{(1 - r_d - b_d)\mathbf{r}_n + r_d\mathbf{u} + b_d\mathbf{b}}{\|(1 - r_d - b_d)\mathbf{r}_n + r_d\mathbf{u} + b_d\mathbf{b}\|} \quad (2) $$

The directional cue was only implemented for the ellipsoidal shape, where it was assumed to be circumferential (orthoradial) to the axis of revolution (long axis).

**Interaction with a pre-existing microtubule.** When a plus end is closer than 25 nm to another (pre-existing) microtubule, it follows one of two behaviors according to the angle between the two microtubules:

* if the angle is smaller than a threshold $\alpha$, the microtubule begins zipping, and its direction becomes permanently that of the pre-existing microtubule,

* if the angle is larger than $\alpha$, the microtubule begins to shrink at the plus end (catastrophe event). It will continue shrinking by one vector at each of the following timesteps until it disappears.
$α$ was set to 40° following experimental observations [9].

**Growth near the surface.** When a microtubule reaches the plasma membrane, it follows one of two behaviors depending on the interaction strength with the membrane.

* Strong anchoring—The microtubule growth is limited to the local plane tangential to the surface, which is calculated at each time step. The new direction will only have two degrees of freedom, corresponding to the projection of the unit vector (Eqs 1 and 2) on the tangential plane. For this reason, microtubules remain at the surface of the cell. This is equivalent to previous models constraining microtubules on the cell surface [10, 11, 41, 51].

* Weak anchoring—Microtubules may grow inside the cell. If a microtubule reaches a distance of 10 nm from the membrane, the angle between the microtubule and the local tangential plane is calculated. We use the same rule as for the interaction with another microtubule. If the angle is greater than $α$, then the microtubule starts shrinking at the plus end; it will continue shrinking by one vector at each of the following timesteps until it disappears. If the angle is smaller than $α$, then the new vector is computed according to Eqs (1) and (2). If the vector points outside the cell then it is re-computed (with new random choices of **u**) until it points inside the cell or tangentially. Thus, the microtubules remain encased in the cell at all times.

## Cell shape and interaction with the membrane

The cell contour is described with a triangular mesh. Each vertex is endowed with the information of the vector normal to the surface, which is used during the simulation to calculate a local approximation of the tangential plane. It is possible to add other informations at the vertex level that can be read during the simulation and serve as extrinsic input. Inputs can be scalars or tensors. The distance from any point in space to the membrane is calculated using the nearest point on the surface. At this point, the membrane is approximated by the plane perpendicular to the normal of the mesh. The distance between a tubulin ring (unit vector) and the membrane is calculated as the shortest distance between the endpoint of the vector and this plane. A collection of standard cell shapes was generated for our simulations (Figs 1 and 6). The main shapes were constructed starting from square parallelepipeds of dimensions 8.8µm×8.8µm×8.8µm (cube), 9µm×9µm×4.7µm ("square"), and 4.8µm×4.8µm×15.6µm ("long"), which are comparable to typical plant cell dimensions. Then the meshes were smoothed so that the minimal radius of curvature was 1.3µm and 4.7µm for sharp and smooth shapes, respectively, which roughly spans the range 0.5-5µm measured in plant roots [10]. The ellipsoid shape is an ellipsoid of revolution around the long axis; the short and long axis have dimensions of 10.3 µm and 16.8 µm, respectively.

## The simulation process

The simulation progress is made through discrete timesteps. At each timestep new vectors are added and vectors are removed from the simulation space according to the rules specified in the previous subsections (Nucleation and minus end behavior; Plus end growth). Collision tests are performed so as to implement the different growth or shrinkage rules. In order to increase the simulation speed, space is divided into subelements and the vectors are identified according to the subelement to which they belong, which diminishes the number of particles involved in the collision test (locality-sensitive hashing [73]).
## Visualization

To visualize microtubule density and orientations an image is created using a matrix of resolution ($r_x, r_y, r_z$). $r_z$ is larger than $r_x = r_y$ by typically a ratio of 10 to 1, which mimics the anisotropic resolution of a confocal microscope. The simulation space is then screened. When a vector is located inside a cube of the matrix, the value of this cube is incremented by one, and the immediate neighbouring cubes are incremented by a lower number (typically 0.3). At the end of the process, a stack is formed where microtubules appear as blurred intensity signals. One can either visualize each sub-image from the stack by moving along the z axis, or create a projection that sums the matrix along the z axis.

## Quantifying anisotropy of microtubule arrays

We used the standard nematic order parameter to quantify anisotropy. The space is subdivided into cubes of arbitrary size, typically segmenting the structure into circa 27 cubes (S2 Fig). Segmenting the structure into 216 smaller pieces gives similar trends, with globally higher anisotropy values. All tubulin ring directions are extracted as a 3 × N matrix, D. We then compute the square symmetric (3 × 3) matrix $M = D^T \cdot D/N$, where T stands for the transpose. M is diagonalised, yielding three eigenvalues $\lambda_i, i \in 1, 2, 3$. Local anisotropy, A, of the microtubule array in each cube of space is defined as

$$A = \sqrt{\frac{3}{2} \frac{\sqrt{(\lambda_1 - \lambda_m)^2 + (\lambda_2 - \lambda_m)^2 + (\lambda_3 - \lambda_m)^2}}{\sqrt{\lambda_1^2 + \lambda_2^2 + \lambda_3^2}}}; \quad \lambda_m = \frac{\lambda_1 + \lambda_2 + \lambda_3}{3}. \quad (3)$$

The value of A is such that $0 \le A \le 1$. Anisotropy value is then computed as the average of the local value A over the whole cell.

## Statistics

The simulations were run 5 times for each parameter value or shape considered; to avoid artificial correlations between data, only the last snapshot was considered for further statistical analysis. As data for the three default values of $n_p$ were statistically identical, they were pooled together for Figs 1–5. The plots were produced with the boxplot function of R: the boxes extend between the first and third quartiles, the segments in the box indicate the medians, and the whiskers are representative of extreme values.

## Supporting information

**S1 Fig. Influence of directional cues on the average anisotropy of the network.** Simulations in an ellipsoidal cell with a circumferential cue. (A) Anisotropy as a function of the strength of the signal, $b_d$. (B-E) Snapshots showing the microtubules in half the cell (the half to the back was removed for clarity). The strength, $b_d$, of the signal is: 0.1% (B), 0.2% (C), 1% (D), and 2% (E), respectively. Model parameters have the default values (see Methods) except for the variable $b_d$ and for weak anchoring.

(PDF)

**S2 Fig. Calculation of anisotropy.** (A) Smooth “square” cell shape. (B) Microtubules. (C) Unit vectors are shown as small white arrows plotted at the center of each cube. (D) The matrix M corresponding to each cube (see Methods) is diagonalised; the red vectors point in the direction of the eigenvectors and the length of the red vectors is proportional to the corresponding eigenvalues.

(PDF)
**S1 Video. Typical simulation.** Smooth “square” cell shape with the default parameter for microtubule dynamics and weak anchoring. The duration of the movie is 5000 timesteps.
(AVI)

**S2 Video. Typical simulation with external cue.** Ellipsoid cell shape with the following parameters for microtubule dynamics: $n_p = 2.4 \cdot 10^{-7}$, $n_s = 10^{-3}$, $\alpha = 40.1^{\circ}$, weak anchoring, $b_d = 0.02$; the duration of the movie is 5000 time steps. The microtubules are shown in half the cell (the half to the back was removed for clarity).
(AVI)

## Acknowledgments

We thank the members of RDP and Sainsbury laboratories for discussions and support.

## Author Contributions

**Conceptualization:** Vincent Mirabet, Olivier Hamant, Elliot M. Meyerowitz, Henrik Jönsson, Arezki Boudaoud.

**Formal analysis:** Vincent Mirabet.

**Funding acquisition:** Olivier Hamant, Elliot M. Meyerowitz, Henrik Jönsson, Arezki Boudaoud.

**Investigation:** Vincent Mirabet.

**Methodology:** Vincent Mirabet.

**Resources:** Pawel Krupinski.

**Software:** Vincent Mirabet.

**Supervision:** Arezki Boudaoud.

**Validation:** Vincent Mirabet.

**Visualization:** Vincent Mirabet.

**Writing – original draft:** Vincent Mirabet, Olivier Hamant, Arezki Boudaoud.

**Writing – review & editing:** Vincent Mirabet, Olivier Hamant, Elliot M. Meyerowitz, Henrik Jönsson, Arezki Boudaoud.

## References

1. Ledbetter MC. A "microtubule" in plant cell fine structure. The Journal of Cell Biology. 1963; 19(1): 239-250. https://doi.org/10.1083/jcb.19.1.239 PMID: 19866635

2. Brandizzi F, Wasteneys GO. Cytoskeleton-dependent endomembrane organization in plant cells: an emerging role for microtubules. The Plant Journal. 2013; 75(2):339-49. https://doi.org/10.1111/tpj. 12227 PMID: 23647215

3. Baskin, Meekes, Liang, Sharp. Regulation of growth anisotropy in well-watered and water-stressed maize roots. II. Role of cortical microtubules and cellulose microfibrils. Plant Physiology. 1999; 119(2): 681-92. https://doi.org/10.1104/pp.119.2.681 PMID: 9952465

4. Paredez AR, Somerville CR, Ehrhardt DW. Visualization of cellulose synthase demonstrates functional association with microtubules. Science. 2006; 312(5779):1491-5. https://doi.org/10.1126/science. 1126551 PMID: 16627697

5. Baskin TI. Anisotropic expansion of the plant cell wall. Annual Review of Cell and Developmental Biol- ogy. 2005; 21(1):203-22. https://doi.org/10.1146/annurev.cellbio.20.082503.103053 PMID: 16212493
6. Bringmann M, Landrein B, Schudoma C, Hamant O, Hauser MT, Persson S. Cracking the elusive alignment hypothesis: the microtubule-cellulose synthase nexus unraveled. Trends in Plant Science. 2012; 17(11):666-74. https://doi.org/10.1016/j.tplants.2012.06.003 PMID: 22784824

7. Buschmann H, Lloyd CW. Arabidopsis mutants and the network of microtubule-associated functions. Molecular Plant. 2008; 1(6):888-898. https://doi.org/10.1093/mp/ssn060 PMID: 19825590

8. Williamson R. Alignment of cortical microtubules by anisotropic wall stresses. Australian Journal Of Plant Physiology. 1990; 17(6):601-613. https://doi.org/10.1071/PP9900601

9. Dixit R, Cyr R. Encounters between dynamic cortical microtubules promote ordering of the cortical array through angle-dependent modifications of microtubule behavior. The Plant Cell. 2004; 16(12):3274-84. https://doi.org/10.1105/tpc.104.026930 PMID: 15539470

10. Ambrose C, Allard JF, Cytrynbaum EN, Wasteneys GO. A CLASP-modulated cell edge barrier mechanism drives cell-wide cortical microtubule organization in Arabidopsis. Nature Communications. 2011; 2(May):430. https://doi.org/10.1038/ncomms1444 PMID: 21847104

11. Cosentino Lagomarsino M, Tanase C, Vos JW, Emons AMC, Mulder BM, Dogterom M. Microtubule organization in three-dimensional confined geometries: evaluating the role of elasticity through a combined in vitro and modeling approach. Biophysical Journal. 2007; 92(3):1046-57. https://doi.org/10.1529/biophysj.105.076893 PMID: 17098802

12. Castle E. Membrane tension and orientation of structure in the plant cell wall. Journal Of Cellular And Comparative Physiology. 1937; 10(1):113-121. https://doi.org/10.1002/jcp.1030100110

13. Green P, King A. A mechanism for origin of specifically oriented textures in development with special reference to Nitella wall texture. Australian Journal Of Biological Sciences. 1966; 19(3):421-37. https://doi.org/10.1071/BI9660421

14. Preston RD. Cellulose-microfibril-orienting mechanisms in plant cells walls. Planta. 1988; 174(1): 67-74. https://doi.org/10.1007/BF00394875 PMID: 24221419

15. Wasteneys GO, Williamson RE. Microtubule orientation in developing internodal cells of Nitella—a quantitative analysis. European Journal Of Cell Biology. 1987; 43:14-22.

16. Wasteneys GO, Williamson RE. Reassembly of microtubules in Nitella tasmanica: assembly of cortical microtubules in branching clusters and its relevance to steady-state microtubule assembly. Journal of Cell Science. 1989; 93(4):705-714.

17. Fischer K, Schopfer P. Interaction of auxin, light, and mechanical stress in orienting microtubules in relation to tropic curvature in the epidermis of maize coleoptiles. Protoplasma. 1997; 196(1-2):108-116. https://doi.org/10.1007/BF01281064

18. Hejnowicz Z, Rusin A, Rusin T. Tensile Tissue Stress Affects the Orientation of Cortical Microtubules in the Epidermis of Sunflower Hypocotyl. Journal Of Plant Growth Regulation. 2000; 19(1):31-44. https://doi.org/10.1007/s003440000005 PMID: 11010990

19. Hamant O, Heisler MG, Jönsson H, Krupinski P, Uyttewaal M, Bokov P, et al. Developmental patterning by mechanical signals in Arabidopsis. Science. 2008; 322(5908):1650-5. https://doi.org/10.1126/science.1165594 PMID: 19074340

20. Jacques E, Verbelen JP, Vissenberg K. Mechanical stress in Arabidopsis leaves orients microtubules in a'continuous' supracellular pattern. BMC Plant Biology. 2013; 13(1):163. https://doi.org/10.1186/1471-2229-13-163 PMID: 24138025

21. Sampathkumar A, Yan A, Krupinski P, Meyerowitz EM. Physical forces regulate plant development and morphogenesis. Current Biology. 2014; 24(10):R475-R483. https://doi.org/10.1016/j.cub.2014.03.014 PMID: 24845680

22. Bouquin T, Mattsson O, Naested H, Foster R, Mundy J. The Arabidopsis lue1 mutant defines a katanin p60 ortholog involved in hormonal control of microtubule orientation during cell growth. Journal of Cell Science. 2003; 116(Pt 5):791-801. https://doi.org/10.1242/jcs.00274 PMID: 12571277

23. Foissner I, Wasteneys GO. Microtubules at wound sites of Nitella internodal cells passively co-align with actin bundles when exposed to hydrodynamic forces generated by cytoplasmic streaming. Planta. 1999; 208(4):480-490. https://doi.org/10.1007/s004250050585

24. Ambrose JC, Wasteneys GO. CLASP modulates microtubule-cortex interaction during self-organization of acentrosomal microtubules. Molecular Biology of the Cell. 2008; 19(11):4730-4737. https://doi.org/10.1091/mbc.E08-06-0665 PMID: 18716054

25. Fishel EA, Dixit R. Role of nucleation in cortical microtubule array organization: variations on a theme. The Plant Journal. 2013; 75(2):270-277. https://doi.org/10.1111/tpj.12166 PMID: 23464654

26. Chan J, Sambade A, Calder G, Lloyd C. Arabidopsis cortical microtubules are initiated along, as well as branching from, existing microtubules. The Plant Cell. 2009; 21(8):2298-2306. https://doi.org/10.1105/tpc.109.069716 PMID: 19706794
27. Nakamura M, Ehrhardt DW, Hashimoto T. Microtubule and katanin-dependent dynamics of microtubule nucleation complexes in the acentrosomal Arabidopsis cortical array. Nature Cell Biology. 2010; 12(11): 1064-1070. https://doi.org/10.1038/ncb2110 PMID: 20935636

28. Ambrose C, Wasteneys GO. Microtubule initiation from the nuclear surface controls cortical microtubule growth polarity and orientation in Arabidopsis thaliana. Plant & Cell Physiology. 2014; 5(0):1-10.

29. Mitchison T, Kirschner M. Dynamic instability of microtubule growth. Nature. 1984; 312(5991):237-42. https://doi.org/10.1038/312237a0 PMID: 6504138

30. Shaw SL, Kamyar R, Ehrhardt DW. Sustained microtubule treadmilling in Arabidopsis cortical arrays. Science. 2003; 300(5626):1715-8. https://doi.org/10.1126/science.1083529 PMID: 12714675

31. Chomicki G, Wightman R, Turner SR. A specific class of short treadmilling microtubules enhances cortical microtubule alignment. Molecular Plant. 2016; 9(8):1214-1216. https://doi.org/10.1016/j.molp.2016.05.008 PMID: 27235548

32. Wightman R, Turner SR. Severing at sites of microtubule crossover contributes to microtubule alignment in cortical arrays. The Plant Journal. 2007; 52(4):742-751. https://doi.org/10.1111/j.1365-313X.2007.03271.x PMID: 17877711

33. Sambade A, Pratap A, Buschmann H, Morris RJ, Lloyd C. The influence of light on microtubule dynamics and alignment in the Arabidopsis hypocotyl. The Plant Cell. 2012; 24(1):192-201. https://doi.org/10.1105/tpc.111.093849 PMID: 22294618

34. Vineyard L, Elliott A, Dhingra S, Lucas JR, Shaw SL. Progressive transverse microtubule array organization in hormone-induced Arabidopsis hypocotyl cells. The Plant Cell. 2013; 25(2):662-676. https://doi.org/10.1105/tpc.112.107326 PMID: 23444330

35. Baulin VA, Marques CM, Thalmann F. Collision induced spatial organization of microtubules. Biophysical Chemistry. 2007; 128(2-3):231-44. https://doi.org/10.1016/j.bpc.2007.04.009 PMID: 17512654

36. Hawkins RJ, Tindemans SH, Mulder BM. Model for the orientational ordering of the plant microtubule cortical array. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics. 2010; 82(1 Pt 1): 011911. https://doi.org/10.1103/PhysRevE.82.011911 PMID: 20866652

37. Shi XqXq, Ma Yq. Understanding phase behavior of plant cell cortex microtubule organization. Proceedings of the National Academy of Sciences of the United States of America. 2010; 107(26): 11709-11714. https://doi.org/10.1073/pnas.1007138107 PMID: 20547876

38. Tindemans SH, Hawkins RJ, Mulder BM. Survival of the aligned: Ordering of the plant cortical microtubule array. Physical Review Letters. 2010; 104(5):058103. https://doi.org/10.1103/PhysRevLett.104.058103 PMID: 20366797

39. Allard JF, Wasteneys GO, Cytrynbaum EN. Mechanisms of self-organization of cortical microtubules in plants revealed by computational simulations. Molecular Biology of the Cell. 2010; 21(2):278-286. https://doi.org/10.1091/mbc.E09-07-0579 PMID: 19910489

40. Allard JF, Ambrose JC, Wasteneys GO, Cytrynbaum EN. A mechanochemical model explains interactions between cortical microtubules in plants. Biophysical Journal. 2010; 99(4):1082-1090. https://doi.org/10.1016/j.bpj.2010.05.037 PMID: 20712991

41. Eren EC, Dixit R, Gautam N. A three-dimensional computer simulation model reveals the mechanisms for self-organization of plant cortical microtubules into oblique arrays. Molecular Biology of the Cell. 2010; 21(15):2674-84. https://doi.org/10.1091/mbc.E10-02-0136 PMID: 20519434

42. Deinum EE, Tindemans SH, Mulder BM. Taking directions: the role of microtubule-bound nucleation in the self-organization of the plant cortical array. Physical Biology. 2011; 8(5):056002. https://doi.org/10.1088/1478-3975/8/5/056002 PMID: 21791726

43. Lindeboom JJ, Lioutas A, Deinum EE, Tindemans SH, Ehrhardt DW, Emons AMC, et al. Cortical microtubule arrays are initiated from a nonrandom prepattern driven by atypical microtubule initiation. Plant Physiology. 2013; 161(3):1189-201. https://doi.org/10.1104/pp.112.204057 PMID: 23300168

44. Deinum EE, Tindemans SH, Lindeboom JJ, Mulder BM. How selective severing by katanin promotes order in the plant cortical microtubule array. Proceedings of the National Academy of Sciences of the United States of America. 2017; 114(27):6942-6947. https://doi.org/10.1073/pnas.1702650114 PMID: 28630321

45. Zhang J, Guo WH, Wang YL. Microtubules stabilize cell polarity by localizing rear signals. Proceedings of the National Academy of Sciences of the United States of America. 2014; 111(46):16383-8. https://doi.org/10.1073/pnas.1410533111 PMID: 25368191

46. Wasteneys GO, Ambrose JC. Spatial organization of plant cortical microtubules: close encounters of the 2D kind. Trends in Cell Biology. 2009; 19(2):62-71. https://doi.org/10.1016/j.tcb.2008.11.004 PMID: 19144522
47. Stoppin-Mellet V, Gaillard J, Vantard M. Katanin's severing activity favors bundling of cortical microtubules in plants. Plant Journal. 2006; 46(6):1009-1017. https://doi.org/10.1111/j.1365-313X.2006.02761.x PMID: 16805733

48. Eren EC, Gautam N, Dixit R. Computer simulation and mathematical models of the noncentrosomal plant cortical microtubule cytoskeleton. Cytoskeleton. 2012; 69(3):144-154. https://doi.org/10.1002/cm.21009 PMID: 22266809

49. Buxton GA, Siedlak SL, Perry G, Smith MA. Mathematical modeling of microtubule dynamics: insights into physiology and disease. Progress in Neurobiology. 2010; 92(4):478-83. https://doi.org/10.1016/j.pneurobio.2010.08.003 PMID: 20713128

50. Nedelec F, Foethke D. Collective Langevin dynamics of flexible cytoskeletal fibers. New Journal of Physics. 2007; 9. https://doi.org/10.1088/1367-2630/9/11/427

51. Tindemans SH, Deinum EE, Lindeboom JJ, Mulder BM. Efficient event-driven simulations shed new light on microtubule organization in the plant cortical array. Frontiers in Physics. 2014; 2(April):1-15.

52. Oda Y. Emerging roles of cortical microtubule-membrane interactions. Journal of Plant Research. 2017; 5:72.

53. Li S, Lei L, Somerville CR, Gu Y. Cellulose synthase interactive protein 1 (CSI1) links microtubules and cellulose synthase complexes. Proceedings of the National Academy of Sciences of the United States of America. 2012; 109(1):185-90. https://doi.org/10.1073/pnas.1118560109 PMID: 22190487

54. Mei Y, Gao HB, Yuan M, Xue HW. The Arabidopsis ARCP protein, CSI1, which is required for microtubule stability, is necessary for root and anther development. The Plant Cell. 2012; 24(3):1066-80. https://doi.org/10.1105/tpc.111.095059 PMID: 22427339

55. Crowell EF, Timpano H, Desprez T, Franssen-Verheijen T, Emons AM, Höfte H, et al. Differential regulation of cellulose orientation at the inner and outer face of epidermal cells in the Arabidopsis hypocotyl. The Plant Cell. 2011; 23(7):2592-605. https://doi.org/10.1105/tpc.111.087338 PMID: 21742992

56. Chan J, Eder M, Crowell EF, Hampson J, Calder G, Lloyd C. Microtubules and CESA tracks at the inner epidermal wall align independently of those on the outer wall of light-grown Arabidopsis hypocotyls. Journal of Cell Science. 2011; 124(7):1088-1094. https://doi.org/10.1242/jcs.086702 PMID: 21363888

57. Müller S. Universal rules for division plane selection in plants. Protoplasma. 2012; 249(2):239-53. https://doi.org/10.1007/s00709-011-0289-y PMID: 21611883

58. Lindeboom JJ, Nakamura M, Hibbel A, Shundyk K, Gutierrez R, Ketelaar T, et al. A mechanism for reorientation of cortical microtubule arrays driven by microtubule severing. Science. 2013; 342(6163): 1245533. https://doi.org/10.1126/science.1245533 PMID: 24200811

59. Gardiner JC, Harper JD, Weerakoon ND, Collings Da, Ritchie S, Gilroy S, et al. A 90-kD phospholipase D from tobacco binds to microtubules and the plasma membrane. The Plant Cell. 2001; 13(9):2143-58.

60. Bringmann M, Li E, Sampathkumar A, Kocabek T, Hauser MT, Persson S. POM-POM2/CELLULOSE SYNTHASE INTERACTING1 is essential for the functional association of cellulose synthase and microtubules in Arabidopsis. The Plant Cell. 2012; 24(1):163-177. https://doi.org/10.1105/tpc.111.093575 PMID: 22294619

61. Liu Q, Qiao F, Ismail A, Chang X, Nick P. The plant cytoskeleton controls regulatory volume increase. Biochimica et Biophysica Acta (BBA)—Biomembranes. 2013; 1828(9):2111-2120. https://doi.org/10.1016/j.bbamem.2013.04.027

62. Heisler MG, Hamant O, Krupinski P, Uyttewaal M, Ohno C, Jönsson H, et al. Alignment between PIN1 polarity and microtubule orientation in the shoot apical meristem reveals a tight coupling between morphogenesis and auxin transport. PLoS Biology. 2010; 8(10):e1000516. https://doi.org/10.1371/journal.pbio.1000516 PMID: 20976043

63. Fu Y, Xu T, Zhu L, Wen M, Yang Z. A ROP GTPase signaling pathway controls cortical microtubule ordering and cell expansion in Arabidopsis. Current Biology. 2009; 19(21):1827-1832. https://doi.org/10.1016/j.cub.2009.08.052 PMID: 19818614

64. Abley K, De Reuille PB, Strutt D, Bangham A, Prusinkiewicz P, Maree AFM, et al. An intracellular partitioning-based framework for tissue cell polarity in plants and animals. Development. 2013; 140(10):2061-2074. https://doi.org/10.1242/dev.062984 PMID: 23633507

65. Paradez A, Wright A, Ehrhardt DW. Microtubule cortical array organization and plant cell morphogenesis. Current Opinion in Plant Biology. 2006; 9(6):571-578. https://doi.org/10.1016/j.pbi.2006.09.005 PMID: 17010658

66. Uyttewaal M, Burian A, Alim K, Landrein B, Borowska-Wykręt D, Dedieu A, et al. Mechanical stress acts via katanin to amplify differences in growth rate between adjacent cells in Arabidopsis. Cell. 2012; 149(2):439-51. https://doi.org/10.1016/j.cell.2012.02.048 PMID: 22500806
67. Ehrhardt DW, Shaw SL. Microtubule dynamics and organization in the plant cortical array. Annual Review of Plant Biology. 2006; 57(1):859-875. https://doi.org/10.1146/annurev.arplant.57.032905.105329 PMID: 16669785

68. Bicek AD, Tüzel E, Kroll DM, Odde DJ. Analysis of microtubule curvature. Methods in Cell Biology. 2007; 83:237-268. https://doi.org/10.1016/S0091-679X(07)83010-X PMID: 17613311

69. Pallavicini C, Levi V, Wetzler DE, Angiolini JF, Benseñor L, Despósito MA, et al. Lateral motion and bending of microtubules studied with a new single-filament tracking routine in living cells. Biophysical Journal. 2014; 106(12):2625-2635. https://doi.org/10.1016/j.bpj.2014.04.046 PMID: 24940780

70. Brangwynne CP, MacKintosh FC, Weitz DA. Force fluctuations and polymerization dynamics of intra-cellular microtubules. Proceedings of the National Academy of Sciences of the United States of America. 2007; 104(41):16128-16133. https://doi.org/10.1073/pnas.0703094104 PMID: 17911265

71. Bicek AD, Tüzel E, Demtchouk A, Uppalapati M, Hancock WO, Kroll DM, et al. Anterograde microtubule transport drives microtubule bending in LLC-PK1 epithelial cells. Molecular Biology of the Cell. 2009; 20(12):2943-2953. https://doi.org/10.1091/mbc.E08-09-0909 PMID: 19403700

72. Hawkins T, Mirigian M, Selcuk Yasar M, Ross JL. Mechanics of microtubules. Journal of Biomechanics. 2010; 43(1):23-30. https://doi.org/10.1016/j.jbiomech.2009.09.005 PMID: 19815217

73. Gionis A, Indyk P, Motwani R. Similarity search in high dimensions via hashing. VLDB'99 Proceedings of the 25th International Conference on Very Large Data Bases. 1999;99(1):518-529.