Skip to main content

Main menu

  • For Authors
    • Submit a Manuscript
    • Instructions for Authors
  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
    • Focus Collections
    • Classics Collection
    • Upcoming Focus Issues
  • Advertisers
  • About
    • About the Journal
    • Editorial Board and Staff
  • Subscribers
  • Librarians
  • More
    • Alerts
    • Contact Us
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Plant Cell Teaching Tools
    • ASPB
    • Plantae

User menu

  • My alerts
  • Log in

Search

  • Advanced search
Plant Physiology
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Plant Cell Teaching Tools
    • ASPB
    • Plantae
  • My alerts
  • Log in
Plant Physiology

Advanced Search

  • For Authors
    • Submit a Manuscript
    • Instructions for Authors
  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
    • Focus Collections
    • Classics Collection
    • Upcoming Focus Issues
  • Advertisers
  • About
    • About the Journal
    • Editorial Board and Staff
  • Subscribers
  • Librarians
  • More
    • Alerts
    • Contact Us
  • Follow plantphysiol on Twitter
  • Visit plantphysiol on Facebook
  • Visit Plantae
Research ArticleBREAKTHROUGH TECHNOLOGIES
Open Access

VirtualLeaf: An Open-Source Framework for Cell-Based Modeling of Plant Tissue Growth and Development

Roeland M.H. Merks, Michael Guravage, Dirk Inzé, Gerrit T.S. Beemster
Roeland M.H. Merks
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: roeland.merks@sysbio.nl
Michael Guravage
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dirk Inzé
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gerrit T.S. Beemster
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site

Published February 2011. DOI: https://doi.org/10.1104/pp.110.167619

  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading
  • © 2011 American Society of Plant Biologists

Abstract

Plant organs, including leaves and roots, develop by means of a multilevel cross talk between gene regulation, patterned cell division and cell expansion, and tissue mechanics. The multilevel regulatory mechanisms complicate classic molecular genetics or functional genomics approaches to biological development, because these methodologies implicitly assume a direct relation between genes and traits at the level of the whole plant or organ. Instead, understanding gene function requires insight into the roles of gene products in regulatory networks, the conditions of gene expression, etc. This interplay is impossible to understand intuitively. Mathematical and computer modeling allows researchers to design new hypotheses and produce experimentally testable insights. However, the required mathematics and programming experience makes modeling poorly accessible to experimental biologists. Problem-solving environments provide biologically intuitive in silico objects (“cells”, “regulation networks”) required for setting up a simulation and present those to the user in terms of familiar, biological terminology. Here, we introduce the cell-based computer modeling framework VirtualLeaf for plant tissue morphogenesis. The current version defines a set of biologically intuitive C++ objects, including cells, cell walls, and diffusing and reacting chemicals, that provide useful abstractions for building biological simulations of developmental processes. We present a step-by-step introduction to building models with VirtualLeaf, providing basic example models of leaf venation and meristem development. VirtualLeaf-based models provide a means for plant researchers to analyze the function of developmental genes in the context of the biophysics of growth and patterning. VirtualLeaf is an ongoing open-source software project (http://virtualleaf.googlecode.com) that runs on Windows, Mac, and Linux.

Plant scientists have gathered a wealth of knowledge about plant growth and plant genetics. Yet, it is still impossible to predict by which mechanisms genetic changes affect plant growth and development. To unravel how genetic information determines the morphogenesis of plants, analyzing and reconstructing the dynamics of the genetic regulatory networks (Gonzalez et al., 2009) is only the first step (Merks and Glazier, 2005). To mechanistically predict morphogenesis from gene function, it is essential to have a detailed understanding not only of genetics but also of (1) how genetic networks regulate cell behaviors, (2) how cell behaviors lead to tissue growth and patterning, and (3) how the tissue-level phenomena, including spatial patterns (Swarup et al., 2005; Jönsson et al., 2006; Bayer et al., 2009) and strains and stresses (Green, 1999; Hamant et al., 2008), induce responses at the molecular level. How these multilevel feedbacks between levels of organization produce biological function and form is perhaps the most central question in systems biology (Noble, 2006).

Traditionally, biologists are familiar with presenting and developing hypotheses using “block-and-arrow” diagrams, symbolically representing identified and hypothetical gene interactions. However, despite their frequent use, it is not possible to understand the behavior of gene networks just by reading such diagrams, except for very straightforward cases. Mathematical models and computational models can bring the gene networks represented by block-and-arrow diagrams to “life” and can explain potential dynamic behaviors of the network even where precise quantitative information is unavailable. In this way, the biologist may discover that the hypothetical biochemical network cannot exhibit the expected behavior or that it occurs only for biologically unfeasible values of the reaction rates or binding constants. The researcher can then propose new interactions or system components and add them to the model to check if the revised network better reproduces the observations the original diagram was attempting to explain. Subsequent experiments will test the validity of the new elements in the hypothesis and can reveal a set of new observations that will be used to refine the mathematical model. This iterative approach of alternating experimental steps and model refinements is known as the “systems biology cycle” (Kitano, 2002). Mathematical and computational modeling supports experimental biologists in testing alternative hypotheses for “internal consistency” and helps them fine-tune their ideas.

A typical systems biology approach is taken in analyzing the de novo patterning of root hairs and in analyzing the regular arrangement of trichomes in the leaf epidermis of Arabidopsis (Arabidopsis thaliana). These two processes are controlled by the same molecular regulators (Schnittger et al., 1999; Pesch and Hülskamp, 2004; Kang et al., 2009). Based on a series of knockout experiments, an initial block-and-arrow model for the interaction of the transcription factors TTG, GL1, and TRY was proposed (Schnittger et al., 1999). By applying a local activation and lateral inhibition model for pattern formation, it was predicted that TRY is activated by GL1 and/or TTG and inhibits GL1 and/or TTG, two predictions that subsequently could be confirmed experimentally. Newly identified components have been added to the initial model (Pesch and Hülskamp, 2004), which was fine-tuned in a series of computer simulation studies (Benítez et al., 2007; Bouyer et al., 2008; Dupuy et al., 2008) that were validated experimentally (Kang et al., 2009). Ongoing cycles of computational modeling and experimental validation continue to fine-tune an awareness of epidermal patterning in Arabidopsis (Bouyer et al., 2008).

The analysis of epidermal patterning is one of several specific questions in plant morphogenesis that could be addressed with this approach, including gravitropism (Swarup et al., 2005), phyllotaxis (de Reuille et al., 2006; Jönsson et al., 2006; Smith et al., 2006), root tip specification (Grieneisen et al., 2007), leaf venation (Merks et al., 2007), mechanical signaling at the shoot apex (Hamant et al., 2008), lateral root initiation (Laskowski et al., 2008; Dun et al., 2009; Prusinkiewicz et al., 2009), and patterning of the cell wall (Tindemans et al., 2010). These projects demonstrate how plant biologists benefit from computational modeling by developing new ideas, insights, and hypotheses and how it helps predict missing system components. To date, modeling typically requires close collaboration between wet-lab scientists and computational biologists or mathematicians, even for testing relatively straightforward ideas or variants of existing models for which no new modeling and simulation codes and techniques need to be developed. Developments in bioinformatics demonstrate that, presumably, this only represents a transient state in the development of modeling. Initially, only computer scientists or computational biologists developed and applied sequence-based algorithms. Since then, experimental biologists have integrated several now-standard bioinformatics tools in their research, while computer scientists and mathematicians continue to develop new algorithms and methods. For example, to perform a BLAST search, experimental biologists no longer require collaborations with an expert bioinformatician. We hope that experimental biologists will soon also have a standard set of tools available for testing simple dynamic hypotheses on morphogenesis (Merks et al., 2006b). Apart from developing their own dynamic models and testing those in the laboratory, it will also help experimental biologists develop an intuition of how computational models can help them advance their research and assist them in framing a set of specific questions to be answered using more advanced models in collaboration with computational biologists.

A number of simulation systems exist for modeling multicellular plant systems. A popular and mature system is L-studio (Karwowski and Prusinkiewicz, 2004). As one of its features, L-studio implements vv-systems, a two-dimensional rewriting grammar to model cell division (Smith, 2006) derived from L-systems (Lindenmayer, 1968a, 1968b, 1975). L-studio and vv-systems have been applied in a number of recent studies on plant development (Smith et al., 2006; Bayer et al., 2009; Prusinkiewicz et al., 2009). L-studio has an intuitive user interface, and users can specify models using L-system or vv-system rules. In vv-systems, a morphological transformation of the tissue as a whole is specified by the user or by a model variable (Smith et al., 2006). The cell division algorithm then partitions the resulting space. This is a powerful approach if one is primarily interested in how growth can drive pattern formation (i.e. by increasing the domain in which pattern formation takes place; Crampin et al., 2002), such as leaf development (Runions et al., 2005), phyllotaxis (Jönsson et al., 2006; Smith et al., 2006), or trichome initiation (Dupuy et al., 2008) or if one is interested in how stable patterns can drive morphogenesis entirely downstream (see the emerging primordia in the phyllotaxis model of Smith et al., 2006). In contrast to such “morphostatic” development (Salazar-Ciudad et al., 2003), in many cases morphogenesis is a closed feedback loop of pattern formation, growth, and morphogenesis (Holloway, 2010); here, cell divisions guided by spatiotemporal patterns lead to domain growth that, in turn, changes pattern formation. In such “morphodynamic” systems (Salazar-Ciudad et al., 2003), the morphology emerges from a feedback loop between pattern formation and simultaneous morphogenesis, while chemical patterns drive morphogenesis by locally affecting cell behavior.

Other popular multicellular modeling techniques include the cell-based models (Merks and Glazier, 2005). Cell-based models derive morphogenesis from the genetically regulated behavior of individual cells that can be modeled as individual particles or as collections of particles. Chaste (for “cancer, heart, and soft-tissue environment”; Pitt-Francis et al., 2009) provides a set of libraries for developing biological simulations of animal tissues, including a cell-based simulation system for modeling tissue growth. It represents cells by its centers and connects cells with virtual springs. A disadvantage of such single-particle methods is that they cannot represent any morphological, subcellular detail, including cell polarity or cell shape. An example of a multiparticle method is the cellular Potts model (CPM; Graner and Glazier, 1992; Glazier and Graner, 1993). The CPM was developed originally for simulating animal development (for review, see Merks and Glazier, 2005) and was recently applied to plants (Grieneisen et al., 2007). In their model of root development, Grieneisen and coworkers (2007) propose that the auxin maximum at the root tip determines the location of the quiescent center in the root apical meristem. The auxin maximum results from a dynamic equilibrium between central auxin transport toward the root tip (acropetal transport), lateral auxin transport toward the shoot (basipetal transport), and auxin reflux from the lateral epidermis toward the stele at the root basal meristem. The same polar auxin transport mechanism also forms an auxin gradient along the root meristem, which is thought to regulate cell division, cell elongation, and cell maturation in the growing root via the PLT1 to PLT4 genes (Galinha et al., 2007). Although most of the work focuses on polar auxin transport on static cell lattices or cellular configurations derived from microscopic images, the CPM shows that the mechanism can produce realistic root division patterns while maintaining a dynamically stable auxin maximum at the root apex.

Because Chaste and the CPM in their present forms lack an accurate description of the physical properties of the cell wall, adjacent cells slide along one another, making their use in the modeling of growing plant tissues relatively limited. Therefore, most developmental phenomena simulated with Chaste and the CPM are examples of plastic morphogenesis, as it occurs in most animal tissues where the cells can migrate and move relative to one another. Because plant cells cannot migrate and do not slide along each other (an exception is the intercalary growth of phloem fibers [Ageeva et al., 2005]), plant morphogenesis depends exclusively on patterned cell division and cell expansion. This mode of development is often called symplastic growth (Priestley, 1930; Erickson, 1986).

Apart from “gluing” plant cells to each other and maintaining the shape of plant cells, cell walls determine the mechanical properties of plant tissues. The directed orientation of cell wall microfibrils determines the differential directionality of tensile strength in cell walls, allowing for anisotropic cell expansion that underlies organ morphogenesis (Green, 1980; Baskin et al., 2004). Moreover, cells respond to forces in the tissue, such as by setting their division plane parallel to the strain direction in the shoot apical meristem (Hamant et al., 2008). Thus, both the mechanical properties of the cell wall and the tight intercellular connections are crucial for understanding the morphogenesis of plant tissues.

A number of cell-based modeling methodologies consider detailed cell wall mechanics. CellModeller (http://www.archiroot.org.uk/doku.php/navigation/cellmodeller; Rudge and Haseloff, 2005; Dupuy et al., 2008) is a two-dimensional environment for plant tissue simulation. It has explicit representations of the cells, the cell walls, and the relative interactions between the cells and includes differential equation models of the dynamics of cellular properties, including biochemical networks. CellModeller describes cell walls as viscoelastic rods, which are maintained tensed by the turgor pressure. Cell wall dynamics is simulated with differential equations, while architectural, geometric models describe cell division and cell death. The model of Corson et al. (2009) also represents cell walls as viscoelastic beams while using an energy minimization approach to simulate cell wall dynamics and cell growth.

A shortcoming in most of the current modeling approaches is that they require a high level of mathematical and computer programming skills, which often keeps biologists from using them. Therefore, we introduce VirtualLeaf, a modeling environment that allows biologists to turn their hypotheses of multicellular development and patterning into a dynamic simulation model with little programming effort. For simplicity and computational efficiency, the underlying computational method describes cells in relatively little detail. Nevertheless, the model explicitly considers the balance between turgor pressure and cell wall tension. Our method relates to CellModeller and the method of Corson et al. (2009) in its representation of the cell walls and its consideration of turgor pressure. We illustrate the use of our simulation environment with a number of example simulations, including simple growing tissues and systems in which chemical pattern formation interacts with local growth.

MODEL DESCRIPTION

Cellular Mechanics

Our computational method represents cell walls in a tissue in two dimensions using polygonal finite elements (Fig. 1), an approach similar to that taken in other plant tissue models (Nagai and Honda, 2001; Rudge and Haseloff, 2005; Dupuy et al., 2008; Corson et al., 2009). A cell is bound by a set of elastic cell wall elements. To restrict relative cell movement, adjacent cells share a cell wall, where each cell wall consist of multiple cell wall elements, so that cell walls are flexible. A Monte Carlo-based energy minimization algorithm straightforwardly describes cell behaviors, including expansion, division, and active shape change. This approach conceptually derives from the CPM (Graner and Glazier, 1992; Glazier and Graner, 1993), which is most suitable for plastic, animal development (for review, see Merks and Glazier, 2005). We allow cell walls to move according to rules derived from the classic equations for turgor-driven cell expansion (Lockhart, 1965): (1) the intracellular turgor pressure exerts a uniform force on cell walls, attempting to enlarge cells; (2) the elastic cell walls counteract the turgor pressures; and (3) walls expand irreversibly if stretched over a threshold. Wall relaxation occurs under the influence of enzymes secreted into the wall that break specific bonds that cross-link cellulose microfibrils (Cosgrove, 2000) and is also called cell wall yielding (Lockhart, 1965; Cosgrove, 1985).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Cell-based modeling of plant tissue growth. Polygons represent cells, surrounded by cell walls, which consist of viscoelastic wall elements linked by nodes. A wall furthermore represents the membranes and transporter proteins on either side of the cell wall and the apoplast (i.e. the intercellular space). The algorithm iteratively displaces the nodes to balance the turgor pressure with the wall tension and possible extra constraints (e.g. a cell shape constraint). A, Overview of the simulation algorithm. The algorithm attempts to displace each node in random order; one such cycle is a Monte Carlo Step (MCS). MCSs are repeated until the energy fails to drop significantly. This iteration of MCSs is called a relaxation cycle. Once the system has relaxed, we simulate the biochemical networks and other cell behavior rules. These potentially affect cell properties (e.g. turgor pressure, wall extensibility, etc.), and a new relaxation cycle is required. B, Visualization of the simulation algorithm. Step 1, This move of a cell wall node relaxes local turgor pressure and minimally strains the cell walls; the move will be kept. Step 2, This move compresses the right cell and extends the left cell. The top wall element is overstrained, and the bottom wall element is compressed. Because this potential move would be energetically unfavorable, it will typically be rejected. Step 3, Each cell potentially contains a biochemical network described by a set of differential equations; the chemical species often regulate cell and wall attributes, including turgor pressure and wall extensibility. Step 4, A molecular species (shown in purple) channels back and forth between the cytoplasm and the cell membrane, as described by sets of differential equations. Step 5, Transport of the blue molecular species from cell to cell; diffusion passively transports molecules from high to low concentrations, and molecules at the cell membrane (e.g. PINs) actively transport molecules upstream or downstream.

Based on this set of rules, our algorithm defines for each cell i a resting area AT(i) at which the cell’s turgor pressure balances the ambient pressure and a resting length LT(j) for each wall element j; this is the length a section of cell wall would assume in the absence of turgor pressure. The actual cell area is a(i), and actual wall-element lengths are denoted l(j). We can describe the balance between turgor pressure and cell wall resistance in terms of a generalized potential energy or Hamiltonian (H):Embedded Imagewhere indices i and j sum over all cells and polygon edges, respectively, λA is a parameter setting the cells’ resistance to compression or expansion, and λM is a spring constant.

The simulation uses Metropolis dynamics (Metropolis et al., 1953) to minimize the potential energy H: The algorithm iteratively selects a random node and attempts to move it in a random direction Embedded Image, where Embedded Image is a random vector chosen uniformly within the unit circle (i.e. ρ ∈[0, 1] and θ ∈[0, 2π]), and ξ is the step size that is fixed for a given simulation. The algorithm calculates the energy difference ΔH associated with the attempt and always accepts the move if it leads to an energy drop. To prevent the system from getting stuck in local energy minima, we also accept moves producing a small energy increase according to the Boltzmann probability function, P(ΔH) = e−(ΔH/T). T is a parameter setting the amount of “noise” added in this way. To model wall relaxation (cell wall yielding), we introduce new nodes whenever a polygon edge’s length is stretched to the extent that it exceeds a threshold value ϑyield, typically 4 times the original length of the wall element. After the algorithm has attempted to move each of the nodes in the simulation in random order (a Monte Carlo step), it checks if the forces described in the Hamiltonian have balanced out sufficiently. Assuming that shape relaxation is fast compared with biological changes, the algorithm repeats this procedure until turgor pressures, cell expansion, and the resulting cell displacements have equilibrated sufficiently, and the energy difference drops below a relaxation threshold (ΔH < θH). After this relaxation cycle has completed, the algorithm solves the biologically motivated processes, including turgor pressure increases, cell divisions, active cell shape changes, and transport of chemical signals.

Thanks to the energy minimization formalism used for the model, we can include additional cellular properties in a straightforward way by describing them as energy constraints (e.g. body forces like gravity). For example, the length of cells (Merks et al., 2006a) can be constrained using:Embedded Imagewhere H′ is the extended Hamiltonian, e(j) is the actual length of cell j, E(j) is the resting length of cell j, and λE gives the strength of the length constraint and can be interpreted as a deformability parameter. Assuming that cell shapes are ellipsoid, we estimate the cell lengths from the cells’ inertia tensors as Embedded Image with λb(j), the largest eigenvalue of cell j’s inertia tensor:Embedded Imagewhere ∫∫D is the area integral over the cell and {xc, yc} is the cell's center of mass. Using Green’s theorem, the inertia tensor derives from the vertices of the polygon describing the cell. The algorithm reconstructs the inertia tensor from the raw moments of the cells’ coordinates, quantities that can be efficiently updated using local information after each move of the cell boundary. Thus, this algorithm produces a computationally efficient, instantaneous estimate of the cell lengths.

For cell division, we define a division axis and build a new cell wall (consisting of several wall elements) shared by the two daughter cells. The nodes of the parent cell are assigned to the two daughter cells. Typically, cell division is over the shortest axis, which is derived from the eigenvectors of the cell’s inertia tensor (Eq. 3).

Transport of Chemical Signals and Metabolites

Sets of ordinary differential equations in each cell describe the biochemical networks and genetic regulatory networks. Additional differential equations model diffusive transport across the two cell membranes and the cell wall separating adjacent cells. To model active transport by transport proteins (e.g. PINs), we define an extra object “wall” (Fig. 1B), which typically spans several “wall elements.” It stores the amounts of transport proteins present in the cell membranes at either side of the walls. Thus, the export rate of chemical A from cell 1 to an adjacent cell 2 would depend on the concentration of A in cell 1 and the amount of A’s export protein in cell 1’s membrane adjacent to cell 2. Although diffusion between apoplast compartments is not implemented yet, the wall elements can also describe an extracellular space (apoplast). For simplicity, we do not consider apoplasts in the examples described here. After each relaxation cycle (Fig. 1A), we dynamically construct a set of differential equations describing the intracellular reactions and transport of chemicals (with one equation for each chemical per cell and one equation for each chemical per cell wall) and numerically integrate them over a period Δt. This setup implies that the set of differential equations increases in size after each cell division, because division adds new cells and new cell walls. As a numerical integrator, we use a fifth-order, adaptive-step-size Runge-Kutta algorithm (Press et al., 1992).

Modeling Framework

VirtualLeaf separates the conceptual biological model from its implementation in terms of the algorithms for cellular mechanics and chemical transport as well as the visualization and steering of the simulations. The model definition itself is contained in a model plugin, a software module that is developed independently from the main framework. It is defined by combining software building blocks representing objects familiar to the experimental plant biologist, like cells, cell walls, molecules, transporter proteins, or signaling molecules. New model plugins are developed by adapting a small section of C++ code (Supplemental Text S1). The model framework runs fast (minutes rather than hours, depending on model complexity), enabling quick interactive desktop model evaluation. The models can also be run without a user interface, enabling high-throughput parameter screening on a computer cluster.

RESULTS

Based on the developed computational VirtualLeaf framework, it is possible to investigate various aspects of tissue growth. We will illustrate this with a number of example models that investigate emergent behavior of biological systems based on the characteristics attributed to individual component cells.

Example 1: Growth, Margin Stiffness, and Cell Shape

As a first study, we test the behavior of a simple model of plant tissue growth (Fig. 2A; Supplemental Video S1). The simulation was initiated with a single cell, whose target area A [and hence its “turgor pressure” (Ai − ai)2] gradually increases at a constant rate. Cells divide over the short axis once their area a has doubled (for implementation details, see Supplemental Text S1, Tutorial 1). This basic simulation rule creates a rounded “callus” of plant tissue. In these simulations, we assumed that cell walls at the leaf perimeter were thicker and stiffer than the internal walls [λm(perimeter) = ζλm(internal)], with ζ = 2, the stiffness ratio between internal and perimeter walls. Perimeter walls are marked with blue nodes in Figures 1 and 2 and Supplemental Videos S1 to S3. This assumption agrees with the observation that margin walls in plant leaves are thick relative to walls of internal cells (Reinhardt et al., 2007) and that epidermal tissues are stiff relative to internal tissues in meristems. With these assumptions, our model rules produced a rounded aggregate of plant tissue. Interestingly, after reducing the stiffness of perimeter cell walls [λM(perimeter) = λM(internal)], a more “tumorous” morphology develops (Fig. 2B; Supplemental Video S2). In the absence of a neighboring cell or a thick cell wall resisting turgor, the perimeter cells “bulge” outward and expand faster than the internal tissue. Indeed, in plant leaves, the margin cells play a crucial role in morphogenesis (Reinhardt et al., 2007), which was proposed to be due to a biophysical function of the margin cells. In this simple model of plant tissue growth, we have not included the elongated shape of margin cells and we neglect many other mechanical features of plant cells (e.g. microtubules) and of the whole organ (e.g. the venature). Nevertheless, already in this simple model, the mechanical properties of the margin can regulate the shape of the tissue as a whole by constraining the expansion of the interior tissue. Figure 2C and Supplemental Video S3 illustrate the effect of a cell shape constraint (Eq. 2). In this example, the cells elongate into a direction that depends both on the initial shape of the cells and on the forces that the adjacent cells exert on them. Illustrating the shape constraint, this simulation is not intended to correspond to a real plant organ.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Symplastic cell-based plant tissue modeling. Cells adhere tightly to each other and expand at fixed rate. Cells divide once their area ai has doubled. Simulation states are shown after 50, 100, 150, 200, 250, and 300 relaxation cycles (from left to right). A, Walls at the leaf perimeter have double stiffness relative to internal cell walls [λm(perimeter) = 2 λm(internal)]. B, If perimeter walls have “single stiffness” [λm(perimeter) = λm(internal)], perimeter cells bulge out, generating a tumorous morphology. C, Growing tissue with elongating cells. [See online article for color version of this figure.]

Example 2: Signaling Molecules Regulating Growth and Differentiation

Plant morphogenesis results from regulated cell divisions and pattern formation via polar or diffusive transport of chemical signals. Using the VirtualLeaf environment, it is possible to investigate the interaction between a diffusive morphogen and growth (Fig. 3; Supplemental Video S4). One of the cells produces a morphogen (e.g. cytokinin or auxin; concentrations indicated in shades of gray), which transports diffusively from cell to cell. If the concentration c of the morphogen is higher than a threshold, c > θ1, where θ1 is a concentration threshold, cells expand and divide once the area has doubled. At lower concentrations, θ1 > c > θ2, where θ2 < θ1 is a second concentration threshold, cells only expand. At concentrations c < θ2, cells neither expand nor divide (i.e. they enter a maturation state). If the source cell divides, the source is assigned to one of the daughter cells at random. Note that in this illustration, we do not intend to model a biological system in particular; but even in this simple model, the cells naturally progress from a division phase to an expansion phase and finally enter a maturation phase, analogous to what happens in a root tip, for example (Beemster and Baskin, 1998). The model does not consider gravitropy; hence, the growth tip does not correct for the accidental sharp bend in Figure 3D. A local growth stagnation causes the bend: the morphogen concentration in the large expanding remains too high for the cell to enter division. This example simulation could form the basis for more complicated, realistic models focusing on the mechanisms responsible for morphogen gradient formation (Grieneisen et al., 2007) and on modeling the relation between cell division rate, cell expansion rate, and the size of the meristem and expansion zone (Chavarría-Krauser and Schurr, 2004; Chavarría-Krauser et al., 2005). Supplemental Text S1, Tutorial 2, shows how to construct such plant growth models involving signaling molecules with VirtualLeaf.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Diffusive morphogen stimulates growth. One tip cell produces a morphogen. At concentrations c > θ1, cells expand and divide once the area has doubled, for θ1 > c > θ2, cells only expand, while for c < θ2, cells neither expand nor divide (i.e. they go into a maturation state). Simulation state is shown after 100 (A), 500 (B), 1,000 (C), 1,500 (D), and 2,000 (E) relaxation cycles.[See online article for color version of this figure.]

A more realistic and complicated simulation involving chemical reactions that illustrates the use of the VirtualLeaf framework for developing new hypotheses of biological development involves studies of vascular patterning during Arabidopsis leaf development (Merks et al., 2007). The model predicts the existence of a traveling wave of auxin moving from the leaf tip to the leaf base. Auxin is concentrated into a peak by PIN1, which polarizes toward neighboring cells at a rate proportional to their concentration of auxin, a mechanism introduced by Jönsson et al. (2006) and Smith et al. (2006) for phyllotaxis. In our leaf venation model, auxin stimulates the production of auxin pumps. The pumps polarize toward the next cells and pump it onward. In this way, traveling waves form, which leave behind trails of polarized cells that we hypothesize to differentiate into the vascular system (Fig. 4; Supplemental Video S5). The long-term goal of this continuing project is to reconstruct in the VirtualLeaf environment how leaves develop in the interplay between patterned cell division and cell expansion, the transport of auxin and other signaling molecules between cells, and the physical characteristics of various cell types.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

A traveling wave of auxin forming a polar auxin channel. Auxin (shown in green) flows in from the leaf tip at the purple cell walls. PIN (red) localizes preferentially at cell walls adjacent to cells with the highest auxin concentration, as in recent phyllotaxis models. A, At 3 h, auxin accumulates at the leaf edge, and some auxin leaks into the lower cell layers. PIN1 production rises. B, At 13 h, auxin converges into a small peak, toward which the adjacent cells polarize. C, At 20 h, the auxin channel is completed, traced out by the traveling wave of auxin. See also Merks et al. (2007).

Example 3: Interplay between Growth, Cell Division, and Patterning

The next modeling example illustrates how tissue growth can interact with auxin-driven patterning. Here, each cell implements the upstream auxin transport model (Heisler and Jonsson, 2006; Jönsson et al., 2006; Sahlin et al., 2009) and lets local auxin concentrations drive cell expansion. For simplicity, we assume that auxin is produced along the tissue’s perimeter. Cells divide along their shortest axis when they reach twice their original size. After a small cell cluster has formed, auxin accumulates in the cluster’s center and induces further growth. New auxin maxima appear at the tissue perimeter once sufficient tissue has formed, inducing further, localized growth. This feedback between growth and pattern formation eventually forms a bulbous morphology (Fig. 5; Supplemental Video S6). Supplemental Text S1, Tutorials 3, 4, and 5, show in detail how such a model can be constructed in VirtualLeaf.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Interaction between auxin accumulation and growth. Auxin accumulates according to the model of Jönsson et al. (2006). Auxin induces cell expansion; cells divide over the shortest axis after doubling in size. Simulation state is shown after 10 (A), 30 (B), 50 (C), 70 (D), 90 (E), and 110 (F) simulation cycles.

In the previous sections, we illustrated the use of VirtualLeaf to model the biomechanics of plant tissue growth and its use in modeling pattern formation and morphogen-regulated growth. These relatively simple models give us more insight into interactions that exist between pattern formation, tissue growth, and tissue mechanics that drive the formation of plant growth. Supplemental Text S1 provides a step-by-step introduction to building models in VirtualLeaf.

DISCUSSION

In this paper, we introduce VirtualLeaf, a cell-based modeling framework for the simulation of symplastic growth, the typical mode of development for plants and some animal tissues in which adjacent cells cohere tightly and cannot slide relative to one another. The example simulations, ranging from simple tissue growth to models where growth and biochemical patterning interact dynamically, illustrate how cell-based modeling frameworks like VirtualLeaf can help developmental plant biologists turn a hypothesis into a dynamic model and illustrate which tissue-level patterns follow from the observed or hypothetical behavior of the individual cells (Merks and Glazier, 2005). Thus, modeling can be used to check hypotheses for internal consistency by showing whether a putative mechanism indeed produces the shapes and patterns it was meant to explain.

One of our motivations for developing VirtualLeaf was the observation that plant (symplastic) development seems fundamentally different from animal development, because the tight coherence of cell walls prohibits plant cells from moving relative to each other. Would modeling symplastic development be impossible with a “plastic” cell-based simulation method like the CPM or Chaste? Balter et al. (2007) propose that connecting cellular Potts cells with elastic springs, in an approach much like our “cell wall elements,” would prevent relative cell movement. However, doing so would (1) suppress the most powerful feature of the CPM (relative cell movement) at additional computational cost, while (2) the model would not include the interaction between turgor pressure and cell wall tension, one of the main features of plant development.

A further advantage of using an off-lattice method for modeling plant development is in the descriptions of chemical concentrations. In VirtualLeaf, these are mathematically part of the cells and move along with them; in lattice-based formalisms, including the CPM from which VirtualLeaf was derived, the chemical concentrations are described as parts of the lattice instead. Hence, because cells move relative to the lattice, in lattice-based simulation formalisms, the morphogen concentrations may “lag behind” if cells move fast relative to reaction or transport rates.

Currently, developing a new model for VirtualLeaf requires some C++ programming. Although models are defined using biological concepts and it is possible to modify existing models with only basic experience in programming, the required coding and compilation still are a hurdle for computationally untrained users. We are currently defining and implementing a provisional domain-specific language based on the draft Cell Behavior Ontology (http://bioportal.bioontology.org/ontologies/39336), a well-defined set of terms for describing the behavior of animal, plant, or bacterial cells. A biological modeling language derived from Cell Behavior Ontology would make it possible to define the model entirely in a conceptual language familiar to biologists.

Apart from a new modeling language, future versions of VirtualLeaf will require several improvements. Currently, cells only contain a single set of chemical concentrations, while many simulation models of developmental processes depend on intracellular, nonhomogeneous spatial distributions of molecular signals (Grieneisen et al., 2007; Laskowski et al., 2008). Intracellular, triangular grids will allow us to consider intracellular chemical distributions. We also hope to develop a three-dimensional VirtualLeaf environment to allow the simulation of tissues consisting of multiple cell layers; the required algorithms are relatively straightforward, because entities of our two-dimensional framework (polygons, line pieces, etc.) have one-on-one three-dimensional counterparts (polyhedra, triangles, etc.). Further extensions of VirtualLeaf would allow modeling of the physiological interactions of the growing plant organ with the rest of the plant and its environment. Functional-structural plant models (Godin and Sinoquet, 2005) focus on the growth and functioning of the whole plant in interaction with its environment. A number of functional-structural plant modeling frameworks, including OpenAlea (Pradal et al., 2008) and CrossTalk (Draye and Pagès, 2006), integrate structural whole-plant models, environmental models, and tissue growth models to single multiscale plant models. Interfacing VirtualLeaf with OpenAlea or CrossTalk would involve (1) adding appropriate boundary conditions to simulate the influx and efflux of materials from the rest of the virtual plant and (2) adding the control handles to synchronize the VirtualLeaf simulations with the functional-structural plant model.

We hope that future cell-based modeling frameworks will integrate off-lattice and lattice-based cell-based methodology and interface whole-plant and single-organ approaches, thus enabling researchers to use the most suitable approach for the problem or even making it possible to combine approaches. Crucial steps toward this goal include the development of advanced model interfacing algorithms, generic, standard modeling languages for cell-based modeling, and multiscale model integration platforms like OpenAlea or CrossTalk.

In its present form, VirtualLeaf provides an accessible problem-solving environment (Merks et al., 2006b; Cickovski et al., 2007) that allows biologists with minimal programming skills to turn their hypotheses on the mechanisms of plant development and symplastic tissue growth into dynamic simulations. VirtualLeaf is available as an open-source project at http://virtualleaf.googlecode.com and as Supplemental Code S1. Precompiled executables for Windows (Supplemental Code S2) and MacOSX platforms (Supplemental Code S3) demonstrate VirtualLeaf based on a series of example models. A tutorial (Supplemental Text S1) provides instructions for developing new models. We hope that VirtualLeaf’s open-source model will encourage users to use VirtualLeaf for their modeling projects (see e.g. Merks et al., 2007 and Wabnik et al., 2010) and to share new models and framework extensions with other VirtualLeaf users.

Supplemental Data

The following materials are available in the online version of this article.

  • Supplemental Text S1. Modeling plant development with VirtualLeaf, a tutorial.

  • Supplemental Video S1. Symplastic cell-based plant tissue simulation, as in Figure 2A.

  • Supplemental Video S2. Symplastic cell-based plant tissue simulation, as in Figure 2B.

  • Supplemental Video S3. Symplastic cell-based plant tissue simulation, as in Figure 2C.

  • Supplemental Video S4. Diffusive morphogen stimulates growth, as in Figure 3.

  • Supplemental Video S5. A traveling wave of auxin forming a polar auxin channel, as in Figure 4.

  • Supplemental Video S6. Interaction between auxin accumulation and growth, as in Figure 5.

  • Supplemental Code S1. Source codes of VirtualLeaf version 1.0, written in C++ with Qt 4.6.

  • Supplemental Code S2. Installer for Windows.

  • Supplemental Code S3. Installer for MacOSX.

Footnotes

  • ↵1 This work was supported by the Netherlands Consortium for Systems Biology, which is part of the Netherlands Genomics Initiative/Netherlands Organisation for Scientific Research; by grants from Ghent University (Bijzonder Onderzoeksfonds Methusalem project no. BOF08/01M00408 and Geconcerteerde Onderzoeksacties no. 12051403); by the Interuniversity Attraction Poles Programme (grant no. VI/33), initiated by the Belgian State, Science Policy Office; and by Marie Curie Intra-European Fellowship number MEIF–CT–2005–025084 and Marie Curie European Reintegration Grant number PERG03–GA–2008–230974 to R.M.H.M.

  • ↵2 These authors contributed equally to the article.

  • The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Roeland M.H. Merks (roeland.merks{at}sysbio.nl).

  • ↵[C] Some figures in this article are displayed in color online but in black and white in the print edition.

  • ↵[W] The online version of this article contains Web-only data.

  • ↵[OA] Open Access articles can be viewed online without a subscription.

  • www.plantphysiol.org/cgi/doi/10.1104/pp.110.167619

REFERENCES

  1. ↵
    1. Ageeva MV,
    2. Petrovská B,
    3. Kieft H,
    4. Sal’nikov VV,
    5. Snegireva AV,
    6. van Dam JEG,
    7. van Veenendaal WLH,
    8. Emons AMC,
    9. Gorshkova TA,
    10. van Lammeren AAM
    (2005) Intrusive growth of flax phloem fibers is of intercalary type. Planta 222: 565–574
    OpenUrlCrossRefPubMed
  2. ↵
    1. Anderson ARA,
    2. Chaplain MAJ,
    3. Rejniak KA
    1. Balter A,
    2. Merks R,
    3. Popławski N,
    4. Swat M,
    5. Glazier J,
    6. Rejniak K,
    7. Anderson A,
    8. Chaplain M
    (2007) The Glazier-Graner-Hogeweg model: extensions, future directions, and opportunities for further study. In Anderson ARA, Chaplain MAJ, Rejniak KA, eds, Single-Cell-Based Models in Biology and Medicine. Birkhaüser, Basel, pp 151–167
  3. ↵
    1. Baskin TI,
    2. Beemster GTS,
    3. Judy-March JE,
    4. Marga F
    (2004) Disorganization of cortical microtubules stimulates tangential expansion and reduces the uniformity of cellulose microfibril alignment among cells in the root of Arabidopsis. Plant Physiol 135: 2279–2290
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Bayer EM,
    2. Smith RS,
    3. Mandel T,
    4. Nakayama N,
    5. Sauer M,
    6. Prusinkiewicz P,
    7. Kuhlemeier C
    (2009) Integration of transport-based models for phyllotaxis and midvein formation. Genes Dev 23: 373–384
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Beemster GTS,
    2. Baskin TI
    (1998) Analysis of cell division and elongation underlying the developmental acceleration of root growth in Arabidopsis thaliana. Plant Physiol 116: 1515–1526
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Benítez M,
    2. Espinosa-Soto C,
    3. Padilla-Longoria P,
    4. Díaz J,
    5. Alvarez-Buylla ER
    (2007) Equivalent genetic regulatory networks in different contexts recover contrasting spatial cell patterns that resemble those in Arabidopsis root and leaf epidermis: a dynamic model. Int J Dev Biol 51: 139–155
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bouyer D,
    2. Geier F,
    3. Kragler F,
    4. Schnittger A,
    5. Pesch M,
    6. Wester K,
    7. Balkunde R,
    8. Timmer J,
    9. Fleck C,
    10. Hülskamp M
    (2008) Two-dimensional patterning by a trapping/depletion mechanism: the role of TTG1 and GL3 in Arabidopsis trichome formation. PLoS Biol 6: e141
    OpenUrlCrossRefPubMed
  8. ↵
    1. Chavarría-Krauser A,
    2. Jäger W,
    3. Schurr U
    (2005) Primary root growth: a biophysical model of auxin-related control. Funct Plant Biol 32: 849–862
    OpenUrlCrossRef
  9. ↵
    1. Chavarría-Krauser A,
    2. Schurr U
    (2004) A cellular growth model for root tips. J Theor Biol 230: 21–32
    OpenUrlCrossRefPubMed
  10. ↵
    1. Cickovski T,
    2. Aras K,
    3. Alber MS,
    4. Izaguirre JA,
    5. Swat M,
    6. Glazier JA,
    7. Merks RMH,
    8. Glimm T,
    9. Hentschel HGE,
    10. Newman SA
    (2007) From genes to organisms via the cell: a problem-solving environment for multicellular development. Comput Sci Eng 9: 50–60
    OpenUrlPubMed
  11. ↵
    1. Corson F,
    2. Hamant O,
    3. Bohn S,
    4. Traas J,
    5. Boudaoud A,
    6. Couder Y
    (2009) Turning a plant tissue into a living cell froth through isotropic growth. Proc Natl Acad Sci USA 106: 8453–8458
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Cosgrove DJ
    (1985) Cell wall yield properties of growing tissue: evaluation by in vivo stress relaxation. Plant Physiol 78: 347–356
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Cosgrove DJ
    (2000) Loosening of plant cell walls by expansins. Nature 407: 321–326
    OpenUrlCrossRefPubMed
  14. ↵
    1. Crampin EJ,
    2. Gaffney EA,
    3. Maini PK
    (2002) Mode-doubling and tripling in reaction-diffusion patterns on growing domains: a piecewise linear model. J Math Biol 44: 107–128
    OpenUrlCrossRefPubMed
  15. ↵
    1. de Reuille PB,
    2. Bohn-Courseau I,
    3. Ljung K,
    4. Morin H,
    5. Carraro N,
    6. Godin C,
    7. Traas J
    (2006) Computer simulations reveal properties of the cell-cell signaling network at the shoot apex in Arabidopsis. Proc Natl Acad Sci USA 103: 1627–1632
    OpenUrlAbstract/FREE Full Text
    1. Draye X,
    2. Pagès L
    (2006) CrossTalk: a simulation platform for the linking of existing soil, plant and atmosphere models. In PMA ’06: Second International Symposium on Plant Growth Modeling and Applications. IEEE, Los Alamitos, CA, pp 93–100
  16. ↵
    1. Dun EA,
    2. Hanan J,
    3. Beveridge CA
    (2009) Computational modeling and molecular physiology experiments reveal new insights into shoot branching in pea. Plant Cell 21: 3459–3472
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Dupuy L,
    2. Mackenzie J,
    3. Rudge T,
    4. Haseloff J
    (2008) A system for modelling cell-cell interactions during plant morphogenesis. Ann Bot (Lond) 101: 1255–1265
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Erickson R
    (1986) Symplastic growth and symplasmic transport. Plant Physiol 82: 1153
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Galinha C,
    2. Hofhuis H,
    3. Luijten M,
    4. Willemsen V,
    5. Blilou I,
    6. Heidstra R,
    7. Scheres B
    (2007) PLETHORA proteins as dose-dependent master regulators of Arabidopsis root development. Nature 449: 1053–1057
    OpenUrlCrossRefPubMed
  20. ↵
    1. Glazier JA,
    2. Graner F
    (1993) Simulation of the differential adhesion driven rearrangement of biological cells. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 47: 2128–2154
    OpenUrlCrossRefPubMed
  21. ↵
    1. Godin C,
    2. Sinoquet H
    (2005) Functional-structural plant modelling. New Phytol 166: 705–708
    OpenUrlCrossRefPubMed
  22. ↵
    1. Gonzalez N,
    2. Beemster GTS,
    3. Inzé D
    (2009) David and Goliath: what can the tiny weed Arabidopsis teach us to improve biomass production in crops? Curr Opin Plant Biol 12: 157–164
    OpenUrlCrossRefPubMed
  23. ↵
    1. Graner F,
    2. Glazier JA
    (1992) Simulation of biological cell sorting using a two-dimensional extended Potts model. Phys Rev Lett 69: 2013–2016
    OpenUrlCrossRefPubMed
  24. ↵
    1. Green P
    (1980) Organogenesis: a biophysical view. Annu Rev Plant Physiol 31: 51–82
    OpenUrlCrossRef
  25. ↵
    1. Green P
    (1999) Expression of pattern in plants: combining molecular and calculus-based biophysical paradigms. Am J Bot 86: 1059–1076
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Grieneisen VA,
    2. Xu J,
    3. Marée AFM,
    4. Hogeweg P,
    5. Scheres B
    (2007) Auxin transport is sufficient to generate a maximum and gradient guiding root growth. Nature 449: 1008–1013
    OpenUrlCrossRefPubMed
  27. ↵
    1. Hamant O,
    2. Heisler MG,
    3. Jönsson H,
    4. Krupinski P,
    5. Uyttewaal M,
    6. Bokov P,
    7. Corson F,
    8. Sahlin P,
    9. Boudaoud A,
    10. Meyerowitz EM
    et al (2008) Developmental patterning by mechanical signals in Arabidopsis. Science 322: 1650–1655
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Heisler MG,
    2. Jonsson H
    (2006) Modeling auxin transport and plant development. J Plant Growth Regul 25: 302–312
    OpenUrlCrossRef
  29. ↵
    1. Holloway DM
    (2010) The role of chemical dynamics in plant morphogenesis (1). Biochem Soc Trans 38: 645–650
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Jönsson H,
    2. Heisler MG,
    3. Shapiro BE,
    4. Meyerowitz EM,
    5. Mjolsness E
    (2006) An auxin-driven polarized transport model for phyllotaxis. Proc Natl Acad Sci USA 103: 1633–1638
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Kang YH,
    2. Kirik V,
    3. Hulskamp M,
    4. Nam KH,
    5. Hagely K,
    6. Lee MM,
    7. Schiefelbein J
    (2009) The MYB23 gene provides a positive feedback loop for cell fate specification in the Arabidopsis root epidermis. Plant Cell 21: 1080–1094
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Godin C,
    2. Hanan J,
    3. Kurth W,
    4. Lacointe A,
    5. Takenaka A,
    6. Prusinkiewicz P,
    7. DeJong T,
    8. Beveridge C,
    9. Andrieu B
    1. Karwowski R,
    2. Prusinkiewicz P
    (2004) The L-system-based plant-modeling environment L-studio 4.0. In Godin C, Hanan J, Kurth W, Lacointe A, Takenaka A, Prusinkiewicz P, DeJong T, Beveridge C, Andrieu B, eds, Proceedings of the 4th International Workshop on Functional-Structural Plant Models. UMR AMAP, Montpellier, France, pp 403–405
  33. ↵
    1. Kitano H
    (2002) Systems biology: a brief overview. Science 295: 1662–1664
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Laskowski M,
    2. Grieneisen VA,
    3. Hofhuis H,
    4. Hove CA,
    5. Hogeweg P,
    6. Marée AFM,
    7. Scheres B
    (2008) Root system architecture from coupling cell shape to auxin transport. PLoS Biol 6: e307
    OpenUrlCrossRefPubMed
  35. ↵
    1. Lindenmayer A
    (1968a) Mathematical models for cellular interactions in development. I. Filaments with one-sided inputs. J Theor Biol 18: 280–299
    OpenUrlCrossRefPubMed
  36. ↵
    1. Lindenmayer A
    (1968b) Mathematical models for cellular interactions in development. II. Simple and branching filaments with two-sided inputs. J Theor Biol 18: 300–315
    OpenUrlCrossRefPubMed
  37. ↵
    1. Lindenmayer A
    (1975) Developmental algorithms for multicellular organisms: a survey of L-systems. J Theor Biol 54: 3–22
    OpenUrlCrossRefPubMed
  38. ↵
    1. Lockhart JA
    (1965) An analysis of irreversible plant cell elongation. J Theor Biol 8: 264–275
    OpenUrlCrossRefPubMed
  39. ↵
    1. Merks RMH,
    2. Brodsky SV,
    3. Goligorksy MS,
    4. Newman SA,
    5. Glazier JA
    (2006a) Cell elongation is key to in silico replication of in vitro vasculogenesis and subsequent remodeling. Dev Biol 289: 44–54
    OpenUrlCrossRefPubMed
  40. ↵
    1. Merks RMH,
    2. Glazier JA
    (2005) A cell-centered approach to developmental biology. Physica A 352: 113–130
    OpenUrlCrossRef
  41. ↵
    1. Merks RMH,
    2. Hoekstra AG,
    3. Kaandorp JA,
    4. Sloot PMA,
    5. Hogeweg P
    (2006b) Problem-solving environments for biological morphogenesis. Comput Sci Eng 8: 61–72
    OpenUrlCrossRef
  42. ↵
    1. Merks RMH,
    2. Van de Peer Y,
    3. Inzé D,
    4. Beemster GTS
    (2007) Canalization without flux sensors: a traveling-wave hypothesis. Trends Plant Sci 12: 384–390
    OpenUrlCrossRefPubMed
  43. ↵
    1. Metropolis NA,
    2. Rosenbluth AW,
    3. Rosenbluth MN,
    4. Teller AH
    (1953) Equation of state calculation by fast computing machines. J Chem Phys 21: 1087–1092
    OpenUrlCrossRef
  44. ↵
    1. Nagai T,
    2. Honda H
    (2001) A dynamic cell model for the formation of epithelial tissues. Philos Mag 81: 699–719
    OpenUrl
  45. ↵
    1. Noble D
    (2006) The Music of Life. Oxford University Press, Oxford
  46. ↵
    1. Pesch M,
    2. Hülskamp M
    (2004) Creating a two-dimensional pattern de novo during Arabidopsis trichome and root hair initiation. Curr Opin Genet Dev 14: 422–427
    OpenUrlCrossRefPubMed
  47. ↵
    1. Pitt-Francis J,
    2. Pathmanathan P,
    3. Bernabeu MO,
    4. Bordas R,
    5. Cooper J,
    6. Fletcher AG,
    7. Mirams GR,
    8. Murray P,
    9. Osborne JM,
    10. Walter A,
    11. et al.
    (2009) Chaste: a test-driven approach to software development for biological modelling. Comput Phys Commun 180: 2452–2471
    OpenUrlCrossRef
  48. ↵
    1. Pradal P,
    2. Dufour-Kowalski S,
    3. Boudon F,
    4. Fournier C,
    5. Godin C
    (2008) OpenAlea: a visual programming and component-based software platform for plant modeling. Funct Plant Biol 35: 751–760
    OpenUrlCrossRef
  49. ↵
    1. Press WH,
    2. Flann BP,
    3. Teukolsky SA,
    4. Vetterling WT
    (1992) Numerical Recipes in C: The Art of Scientific Computing, Ed 2. Cambridge University Press, Cambridge, UK
  50. ↵
    1. Priestley JH
    (1930) Studies in the physiology of cambial activity. II. The concept of sliding growth. New Physiol 29: 96–140
    OpenUrlCrossRef
  51. ↵
    1. Prusinkiewicz P,
    2. Crawford S,
    3. Smith RS,
    4. Ljung K,
    5. Bennett T,
    6. Ongaro V,
    7. Leyser O
    (2009) Control of bud activation by an auxin transport switch. Proc Natl Acad Sci USA 106: 17431–17436
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Reinhardt B,
    2. Hänggi E,
    3. Müller S,
    4. Bauch M,
    5. Wyrzykowska J,
    6. Kerstetter R,
    7. Poethig S,
    8. Fleming AJ
    (2007) Restoration of DWF4 expression to the leaf margin of a dwf4 mutant is sufficient to restore leaf shape but not size: the role of the margin in leaf development. Plant J 52: 1094–1104
    OpenUrlCrossRefPubMed
  53. ↵
    1. Rudge T,
    2. Haseloff J
    (2005) A computational model of cellular morphogenesis in plants. Lect Notes Comput Sci 3630: 78–87
    OpenUrlCrossRef
  54. ↵
    1. Runions A,
    2. Fuhrer M,
    3. Lane B,
    4. Federl P,
    5. Rolland-Lagan AG,
    6. Prusinkiewicz P
    (2005) Modeling and visualization of leaf venation patterns. ACM Transactions on Graphics 24: 702–711
    OpenUrlCrossRef
  55. ↵
    1. Sahlin P,
    2. Söderberg B,
    3. Jönsson H
    (2009) Regulated transport as a mechanism for pattern generation: capabilities for phyllotaxis and beyond. J Theor Biol 258: 60–70
    OpenUrlCrossRefPubMed
  56. ↵
    1. Salazar-Ciudad I,
    2. Jernvall J,
    3. Newman SA
    (2003) Mechanisms of pattern formation in development and evolution. Development 130: 2027–2037
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Schnittger A,
    2. Folkers U,
    3. Schwab B,
    4. Jürgens G,
    5. Hülskamp M
    (1999) Generation of a spacing pattern: the role of triptychon in trichome patterning in Arabidopsis. Plant Cell 11: 1105–1116
    OpenUrlAbstract/FREE Full Text
  58. ↵
    1. Smith C
    (2006) On vertex-vertex systems and their use in geometric and biological modelling. PhD thesis. University of Calgary, Calgary, Canada
  59. ↵
    1. Smith RS,
    2. Guyomarc’h S,
    3. Mandel T,
    4. Reinhardt D,
    5. Kuhlemeier C,
    6. Prusinkiewicz P
    (2006) A plausible model of phyllotaxis. Proc Natl Acad Sci USA 103: 1301–1306
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Swarup R,
    2. Kramer EM,
    3. Perry P,
    4. Knox K,
    5. Leyser HM,
    6. Haseloff J,
    7. Beemster GTS,
    8. Bhalerao R,
    9. Bennett MJ
    (2005) Root gravitropism requires lateral root cap and epidermal cells for transport and response to a mobile auxin signal. Nat Cell Biol 7: 1057–1065
    OpenUrlCrossRefPubMed
  61. ↵
    1. Tindemans SH,
    2. Hawkins RJ,
    3. Mulder BM
    (2010) Survival of the aligned: ordering of the plant cortical microtubule array. Phys Rev Lett 104: 058103
    OpenUrlCrossRefPubMed
  62. ↵
    1. Wabnik K,
    2. Kleine-Vehn J,
    3. Balla J,
    4. Sauer M,
    5. Naramoto S,
    6. Reinöhl V,
    7. Merks RM,
    8. Govaerts W,
    9. Friml J
    (2010) Emergence of tissue polarization from synergy of intracellular and extracellular auxin signaling. Mol Syst Biol 6: 447
    OpenUrlCrossRefPubMed
View Abstract
PreviousNext
Back to top

Table of Contents

Print
Download PDF
Email Article

Thank you for your interest in spreading the word on Plant Physiology.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
VirtualLeaf: An Open-Source Framework for Cell-Based Modeling of Plant Tissue Growth and Development
(Your Name) has sent you a message from Plant Physiology
(Your Name) thought you would like to see the Plant Physiology web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
VirtualLeaf: An Open-Source Framework for Cell-Based Modeling of Plant Tissue Growth and Development
Roeland M.H. Merks, Michael Guravage, Dirk Inzé, Gerrit T.S. Beemster
Plant Physiology Feb 2011, 155 (2) 656-666; DOI: 10.1104/pp.110.167619

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
VirtualLeaf: An Open-Source Framework for Cell-Based Modeling of Plant Tissue Growth and Development
Roeland M.H. Merks, Michael Guravage, Dirk Inzé, Gerrit T.S. Beemster
Plant Physiology Feb 2011, 155 (2) 656-666; DOI: 10.1104/pp.110.167619
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • MODEL DESCRIPTION
    • RESULTS
    • DISCUSSION
    • Footnotes
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

In this issue

Plant Physiology: 155 (2)
Plant Physiology
Vol. 155, Issue 2
Feb 2011
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
  • Front Matter (PDF)
View this article with LENS

More in this TOC Section

  • Rapid Affinity Purification of Tagged Plant Mitochondria (Mito-AP) for Metabolome and Proteome Analyses
  • GRANAR, a Computational Tool to Better Understand the Functional Importance of Monocotyledon Root Anatomy
  • An Online Database for Exploring Over 2,000 Arabidopsis Small RNA Libraries
Show more BREAKTHROUGH TECHNOLOGIES

Similar Articles

Our Content

  • Home
  • Current Issue
  • Plant Physiology Preview
  • Archive
  • Focus Collections
  • Classic Collections
  • The Plant Cell
  • Plant Direct
  • Plantae
  • ASPB

For Authors

  • Instructions
  • Submit a Manuscript
  • Editorial Board and Staff
  • Policies
  • Recognizing our Authors

For Reviewers

  • Instructions
  • Journal Miles
  • Policies

Other Services

  • Permissions
  • Librarian resources
  • Advertise in our journals
  • Alerts
  • RSS Feeds

Copyright © 2021 by The American Society of Plant Biologists

Powered by HighWire