components of a virtual tissue
play

Components of a virtual tissue Christophe Godin INRIA Project-team - PowerPoint PPT Presentation

Plant Bioinformatics, Systems and Synthetic Biology Summer school Nottingham, 27-31 July 2009 Components of a virtual tissue Christophe Godin INRIA Project-team Virtual Plants Growth areas in plants Shoot apical meristem Root apical


  1. Result of virtual auxin transport on digitized PIN1 maps • � Auxin accumulates at the primordia locations • � Auxin accumulates at the initium location • � Auxin accumulates in the center • � Accumulation patterns do not depend on the location of auxin production

  2. Back to experiment … t = 0 t = 22h NPA 1. The center is not sensitive to auxin with NPA DR5::GFP Addition of auxin 2. Anti-auxin 3. High levels of auxin observed in immunolabelling the CZ of the clv3 mutants

  3. What drives the polarization of PIN pumps ? Integrating dynamics of tissue development

  4. Allocation of PIN to membranes � a i ( t ) � = D � a i ( t ) + � ( P j , i a j ( t ) � i , j a i ( t )) � � a i ( t ) + � P � t j • � Hypothesis 1 : – � Pumps are oriented so that local auxin spots are amplified ( concentration-based hypothesis ) (Jönsson et al. 06, Smith et al., PNAS, 06) P i Available amount of PINs in cell i s i , j Surface between cell i and j (Smith et al., 06)

  5. Simulating tissue growth V = d � � = � � r Velocity field: r P O � dt � � r : relative elementary rate of growth Constant speed Linear speed

  6. Simulating tissue growth � � Velocity Field V = d � � dt = f ( � r r , t ) � � Division rules (Nakielski, …) – � Volume > threshold. – � Location and orientation of the new wall • � Minimal length, • � Right angle between new and old walls.

  7. Concentration-based hypothesis (Smith et al., PNAS, 06)

  8. Candidate hypotheses � a i ( t ) � = D � a i ( t ) + � ( P j , i a j ( t ) � i , j a i ( t )) � � a i ( t ) + � P � t j • � Hypothesis 1 : – � Pumps are oriented so that local auxin spots are amplified ( concentration-based hypothesis ) (Jönsson et al. 06, Smith et al., PNAS, 06) • � Hypothesis 2 : – � Pumps are oriented so that fluxes are amplified ( canalization = flux-based hypothesis ) (Sachs 69, Mitchison 81, Feugier et al. 05, Rolland-Lagand et al. 05) (Runions et al., SIGGRAPH, 05)

  9. Could canalization explain auxin transport in the L1 layer ? Flux-based hypothesis: dP = f ( � i , j ) � � P i , j + � i , j dt f = feedback function quadratic linear (Feugier et weak strong (canalization) al. JTB, 05)

  10. Flux-based polarization allows pumping with or against the auxin gradient Pumping against the Pumping with the gradient gradient (infinite sink strength) (finite sink strength)

  11. Flux-based polarization may create dynamic patterning Decreasing the threshold of primordia initiation

  12. Weak flux-based polarization can create inhibitory fields � P The size of the inhibitory field is a = � � i , j � � P i , j + � i , j � function of the feedback parameter ( ) � t � � � � =1.3 =1.5 =1.7 =2.0

  13. Simulation of auxin fluxes on digitized PIN1 maps - � Auxin is produced and degraded in each cell - � Diffusive and active transport - � Primordia are perfect sinks ? Simulated PIN1 maps Observed PIN1 maps (weak flux-based polarization)

  14. Influence zone of a region Definition : set of cells connected in the map with cells of a given region by an oriented path of pumps Central zone Primordia Observed maps Simulated maps

  15. Role of the central zone Central zone has no distinct Observed map Central zone degrades behaviour auxin

  16. Comparison of the influence zones 15% more pumps are correctly oriented (78% in total) Simulated map with CZ degrading auxin Observed maps Simulated map without CZ

  17. Dynamic simulation of phyllotaxy

  18. Flux-based simulation of phyllotaxy

  19. Simulated divergence angle

  20. Simulation of the generation of provascular tissues

  21. Flux-based simulation of vascularisation

  22. Flux-based polarization makes it possible to pump both with and against the gradient ( Ottenschläger et al. PNAS, 03) DR5::GFP AIA flux all around AIA flux AIA reflux

  23. An alternative dual model (Bayer et al., 2008) Simulated PIN Simulated Auxin

  24. Experimental verification

  25. Summary on transport Concentration-based Flux-based polarization polarization YES YES (weak FBP) Phyllotaxis (Smith et al. 06, Johnson et al. 06) � (Stoma et al. 08) YES (strong FBP) Being investigated/Mixed model Venation patterns (Mitchison 81, Rolland-Lagan 06, Runion 06, Feugier 05) � (Merks et al. 07) � , / (Bayer,08) ? YES (strong FBP) Fountain model (root apex) (Stoma et al. 08) No Molecular interpretation No Assessment (Phyllotaxis): Ok Divergence angles Ok Phyllotactic pattern stability To improve To improve Consistent with observed Fairly consistent / quantitative Partially/qualitative PIN maps if center degrades auxin (role?) Predicted event sequence Maximum is maintained / Maximum / leaks / minimum Pumps pointing upwards initially

  26. Building of a virtual meristem 2 – Transport model 1 – Geometric model 3 – Physical model 4 – Cell model Division and Growth Interaction network + + a h -

  27. Mechanical aspects of growth Cell-cell physical interactions ?

  28. Local/Bottom up specification of growth « The growing Canvas », The art of genes , E. Coen, 1999 « The genetics of geometry », (Coen et al, PNAS, 2004) Shape as an emerging property of region growth …

  29. A general conceptual framework « The genetics of geometry », (Coen et al, PNAS, 2004) Alphabet of elementary geometric transformations : Growth rate Anisotropy Direction Rotation Local information : Global constraints : - Mechanical - � genes activity, - � microtubules, - � fluxes, forces, - hormones - stresses, - … - … - …

  30. Strain description l � l 0 • � Strain in 1D l 0 � = l � l 0 l 0 l • � Strain in 2D � � � xx � xy � � � = � yx � yy � � � � Strain tensor

  31. Elementary transforms in mathematical terms Decomposition of the strain tensor (2D) : � = � � e . s c a l a n i V � V e = � � � 1 + � � 2 . I � = r e f s c a l 2 V r e f i = T � R D R a n T R � scale R D

  32. Development controlled by gene expression « The genetics of geometry », (Coen et al, PNAS, 2004) lobe late- dor- tube ral sal extern intermediate central • � High growth rate: + + + • � High anisotropy: Modeling the growth of a petal shape

  33. Integration of local changes development How to assemble these local changes consistently ?

  34. Deformation constraints L l 0 0 1 2 3 new ? � n l i = L Geometric constraint: i 0 =

  35. Different admissible solutions Different combinations: 0 1 new 2 3 0 1 new 2 3 0 1 new 2 3

  36. Cost of a deformation (Energy) Physical interpretation: Translation W = M g h h Deformation 1 x 2 W = k x 2

  37. Total energy of a transformation new 0 1 2 3 W W W n e w 1 2 W W 3 0 0 1 new 2 3 W = W 0 + W 1 + W 2 + W 3 + W n e w Solution : transformation with minimum energy

  38. Integration • � Set of admissible deformations development cost = W i • � Energy minimization over Use of integration methods : • � mass-spring systems • � finite elements

  39. Mechanics and Differential growth - � Each region grows isotropically - � Geometric anisotropy results from global constraints

  40. Residual stresses Growing “petal” Problem of residual stresses Solution: introduce a feedback of the stress on the growth

  41. Cell wall • � Cell wall : Cosgrove (2001) – � Main determinant of cell shape – � Regularly synthesized by the cell – � Composed of bundles of microfibrils linked together by elastic links • � Mechanical aspects: – � Each microfibril resist axial load P � – � Resistance perpendicular to microfibrils is less important – � Turgor pressure induces cell wall strain

  42. Individual cell growth � = l � l 0 • � Elasticity of a rod : Hook’s law l 0 l 0 l � l 0 � = F s = E � l F • � Cell is elastically deformed by turgor pressure P � > 0 P � = 0 � � � � 1 0 0 P � E � � � 1 = P � = P � I = � � = � � E x � 1 0 � � 0 P � � E � y Stress in the region Elastic strain (Hook’s law)

  43. Individual cell growth • � Cell deformation � > 0 P • � Growth induces plastic deformations (Cosgrove 98,01,03,04)

  44. Taking into account cell growth � � > � • � Cell growth P � > 0 � G = f ( � � ) P � = 0 � G = � � t � Example: � Elastic strain Wall synthesis speed � � 1 0 E � G = � � t P x � 1 0 E y

  45. Mechanical interpretation of growth parameters Growth strain of the reference configuration: � 2 � E 0 y E x + E E x + E � G = � � t P y y � 2 E 2 E x E 0 x y E x + E y � � s c a l e a n i V � V r e f • � Scaling represents the relative variation of volume V r e f • � Anisotropy distributes the growth along the principal axes

  46. Simulation • � Without retroaction • � With retroaction

  47. Role of microtubules in growth Microtubules re-orient according to main stresses (Hamant et al., Sience, 2008)

  48. Cell growth decomposition Cell growth is controled by 2 factors : • � Growth intensity (e.g hormone concentration, gene activity) • � Growth anisotropy (polarization of microtubules)

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend


More recommend