Trichome patterning in Arabidopsis serves as a model system for de novo pattern formation in plants. It is thought to typify the theoretical activator–inhibitor mechanism, although this hypothesis has never been challenged by a combined experimental and theoretical approach. By integrating the key genetic and molecular data of the trichome patterning system, we developed a new theoretical model that allows the direct testing of the effect of experimental interventions and in the prediction of patterning phenotypes. We show experimentally that the trichome inhibitor TRIPTYCHON is transcriptionally activated by the known positive regulators GLABRA1 and GLABRA3. Further, we demonstrate by particle bombardment of protein fusions with GFP that TRIPTYCHON and CAPRICE but not GLABRA1 and GLABRA3 can move between cells. Finally, theoretical considerations suggest promoter swapping and basal overexpression experiments by means of which we are able to discriminate three biologically meaningful variants of the trichome patterning model. Our study demonstrates that the mutual interplay between theory and experiment can reveal a new level of understanding of how biochemical mechanisms can drive biological patterning processes.
Arabidopsis trichomes are leaf hairs formed by single epidermal cells and their patterning serves as a model system for de novo pattern formation in plants. Trichomes are initiated in a regular spacing pattern in a growing field of rapidly dividing cells at the base of the young leaf. All current data indicate that the relative position of newly emerging trichomes is governed by cellular interactions between initially equivalent cells. Genetic and molecular experiments have identified genes either by promoting or inhibiting trichome fate. The trichome‐promoting genes encode the R2R3 MYB transcription factors GLABRA1 (GL1) and MYB23 (Oppenheimer et al, 1991; Kirik et al, 2001, 2005), the bHLH factors GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) (Payne et al, 2000; Zhang et al, 2003; Bernhardt et al, 2005) and the WD40‐repeat protein TRANSPARENT TESTA GLABRA1 (TTG1) (Galway et al, 1994; Walker et al, 1999a). The trichome inhibitors are represented by a class of six single‐repeat MYB‐related transcription factors (Wada et al, 1997a; Schellmann et al, 2002; Kirik et al, 2004b; Wang et al, 2007; Tominaga et al, 2008), most importantly TRIPTYCHON (TRY) (Schellmann et al, 2002) and CAPRICE (CPC) (Wada et al, 1997a). Both groups of genes are initially expressed in all cells. It is thought that differences are generated by an activator–inhibitor‐like mechanism such that the activators form a R2R3 MYB/bHLH/WD40 trimeric complex, which is counteracted by the mobile inhibitors through competition with the R2R3 MYB for binding to bHLH proteins. By integrating the key genetic and molecular data of the trichome patterning system, we develop a new theoretical model that allows the direct testing of the effect of experimental interventions and in the prediction of patterning phenotypes. We test several key assumptions of our model experimentally. We show experimentally that the trichome inhibitor TRIPTYCHON is transcriptionally activated by the known positive regulators GLABRA1 and GLABRA3. By particle bombardment of protein fusions with GFP, we show that TRIPTYCHON and CAPRICE but not GLABRA1 and GLABRA3 can move between cells. Further, we confirm the yeast two‐hybrid interaction between TRY and bHLH by bimolecular fluorescence complementation assays and pull‐down experiments. With both techniques, we discover a new interaction between TRY and GL1. To address the potential of this interaction, we formulate three variants of our model (cf. Figure 1) and challenge them experimentally. In the single competitive inhibition model, TRY binds to free GL3 and prevents GL1 GL3 dimerization. In the double competitive inhibition model, TRY binds to free GL3 or GL1, thereby preventing the interaction between GL1 and GL3. In the uncompetitive inhibition model, TRY binds to the GL1 GL3 dimer and represses its function. All models show similar patterning, yet for different parameter ranges and sensitivities. To determine the biological relevance of the three cases, we perform overexpression experiments. We use the 35S promoter to express TRY and GL3 in all cells and the GL2 promoter to express them trichome specific. All GL2:TRY plants are completely glabrous in two independent lines. 35S:GL3 lines show a higher trichome density in the patterning zone and GL2:GL3 lines exhibit a normal trichome density (cf. Figure 2C and D). We simulate the corresponding overexpression of TRY and GL3 using the three inhibition models (cf. Figure 2E–G). Because the parameters of these models are unknown, we randomly sample the parameter space. By exploiting our experimental results, we formulate a specific set of criteria and test each sample whether these criteria are matched. On the basis of our simulations, we conclude that the uncompetitive scenario is not consistent with our data. Moreover, we find that the single competitive inhibition is sufficient to explain all experimental findings. Our study demonstrates that the mutual interplay between theory and experiment can reveal a new level of understanding of how biochemical mechanisms can drive biological patterning processes.
A new theoretical model of trichome patterning is developed using a combined approach of mathematical modeling and experimental validation.
Several assumptions of the model such as transcriptional regulation of the inhibitors, the mobility of the inhibitors, and molecular interactions are confirmed on a genetic and biochemical level.
Three competitive inhibition scenarios can be distinguished that all have the property to create de novo patterns.
A combination of promoter swapping and basal over‐expression experiments with mathematical modeling identifies the most relevant scenario for trichome patterning.
Most theoretical models in biology aiming to explain de novo creation of regular spacing patterns are based on either of two principles: activator–inhibitor or substrate depletion (Gierer and Meinhardt, 1972). In both cases, pattern formation relies on a dynamic instability first suggested by Turing (1952) in the context of morphogenesis. The initiation of Arabidopsis trichomes on the leaf blade serves as an excellent model system to study the molecular mechanism underlying de novo patterning. Trichomes are leaf hairs derived from epidermal cells that are formed in a regular spacing pattern in a rapidly growing cell layer at the leaf base. New trichomes are formed at a minimal distance of three or four cells from already existing ones and their position is not correlated with any other recognizable positional landmark (Hulskamp et al, 1994). A mechanism by which the spacing pattern is achieved by a conserved cell division pattern was excluded by clonal analysis (Larkin et al, 1996; Schnittger et al, 1999). Thus, all data indicate that a de novo patterning mechanism is operating. Several independent mutational screens have identified genes that appear to act as positive or negative regulators of trichome initiation. The corresponding mutants of the positive regulators have fewer or no trichomes. They include the R2R3 MYB transcription factors GLABRA1 (GL1) and MYB23 (Oppenheimer et al, 1991; Kirik et al, 2001, 2005), the bHLH factors GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) (Payne et al, 2000; Zhang et al, 2003; Bernhardt et al, 2005) and the WD40‐repeat protein TRANSPARENT TESTA GLABRA1 (TTG1) (Galway et al, 1994; Walker et al, 1999a). Yeast two‐hybrid data suggest the formation of a trichome‐promoting trimeric complex, due to binding of one R2R3 MYB factor and TTG1 to a bHLH factor (Payne et al, 2000). The negative regulators TRIPTYCHON (TRY) and CAPRICE (CPC) were initially identified by mutants showing trichome clusters (Schellmann et al, 2002) and a higher trichome density (Wada et al, 1997a), respectively. Both genes encode homologous single‐repeat MYB‐related transcription factors (Wada et al, 1997a; Schellmann et al, 2002). Later, four further homologues were found that act in a partially redundant manner, namely ENHANCER OF TRY AND CPC1 (ETC1) (Kirik et al, 2004a), ETC2 (Kirik et al, 2004b), TRICHOMELESS1 (Wang et al, 2007) and CAPRICE‐LIKE MYB3 (CPL3) (Tominaga et al, 2008). These inhibitors can compete with the R2R3 MYB factor for binding to GL3/EGL3 in yeast three‐hybrid assays (Esch et al, 2003). Therefore, current models assume that these inhibitors counteract the trimeric active complex by this competition mechanism under the assumption that the single‐repeat MYB factor‐containing complex is inactive (Larkin et al, 1996; Scheres, 2002; Marks and Esch, 2003; Pesch and Hulskamp, 2004; Ishida et al, 2008). It is assumed that the generation of the actual spacing pattern is mediated by the movement of the inhibitors. One downstream target of this machinery is the homeobox transcription factor GLABRA2 (GL2), which is thought to trigger the actual process of trichome formation (Rerie et al, 1994; Cristina et al, 1996).
Virtually the same gene cassette is operating during epidermal root hair determination. Root hairs are arranged in cell files and only epidermal cells overlying a cleft between two cortex cells develop into root hairs (Dolan et al, 1994; Berger et al, 1998). Although this suggests that cell fate is determined by position; root hair fate depends largely on the R2R3 MYB factor WEREWOLF, TTG1, GL3/EGL3 and the inhibitors CPC, TRY and ETC1. Taken together, these genes activate the expression of GL2 in non‐root hair cells, where GL2 represses root hair development (Larkin et al, 1996; Scheres, 2002; Marks and Esch, 2003; Pesch and Hulskamp, 2004; Ishida et al, 2008). Thus, in contrast to trichome development on leaves, the default fate in the root epidermis is the formation of root hairs.
In this study, we create a theoretical model based on the current knowledge of the trichome patterning system. We verify experimentally the previous assumption that the expression of the inhibitor TRY is induced by the positive regulators. Further, we find (Payne et al, 2000; Zhang et al, 2003) that the inhibitors TRY and CPC are mobile, whereas GL1 and GL3 are cell autonomous. Protein interaction assays suggest three alternative scenarios by which the positive regulators are inhibited. A combination of overexpression experiments and theoretical modelling allows in the identification of the most relevant inhibition scenario for trichome patterning.
Mathematical modelling and simulation of trichome patterning
We develop our initial model from the following assumptions, most of which are based on already published data in the context of trichome or root hair patterning: (1) the GL1 GL3 complex is considered to be the transcriptionally active complex (Morohashi et al, 2007); (2) similar to GL3 and WEREWOLF in the root hair system, it is assumed that GL1 GL3 activate TRY (Bernhardt et al, 2005); (3) TRY counteracts the GL1 GL3 activity by competing with GL1 for binding to GL3 (Esch et al, 2003); (4) similar to that shown for CPC in the root, it is assumed that TRY protein can move between cells (Kurata et al, 2005). Further, we consider GL3 to be mobile because GL3 can travel between cells in the root (Bernhardt et al, 2005) and GL1 to be immobile (Hulskamp et al, 1994); (5) as observed for the GL3 homologue TT8 and the GL1 homologue TT2, it is assumed that the active complex, GL1 GL3, activates the expression of GL3 (Baudry et al, 2006); (6) TTG1 is not taken into account as it is not essential for trichome formation as indicated by the fact that GL1 and GL3 overexpression or GL1 overexpression in try mutants rescues the ttg1 phenotype (Schnittger et al, 1998); (7) GL2 is considered as a downstream target gene of this machinery (Rerie et al, 1994; Cristina et al, 1996); (8) GL1 and GL3 are expressed ubiquitously in the patterning zone (Larkin et al, 1993; Zhang et al, 2003); (9) in addition, it is assumed that GL1 is activated by the active complex. This is included as a prerequisite for the model; (10) it is further assumed that all proteins are linearly degraded. The resulting interactions are shown in Figure 1 and the corresponding mathematical model is presented in the Materials and methods section.
Simulations of the corresponding differential equations reveal a regular spacing pattern for biologically reasonable parameter ranges. To further validate our model, we directly tested several of our key assumptions that have so far only been based on indirect genetic experiments or made by analogy to the root hair system.
TRY is transcriptionally activated by the active complex and suppresses its own transcription
The transcriptional regulation of TRY was analysed using a TRY:GUS construct containing a promoter previously shown to truly reflect the endogenous expression (Schellmann et al, 2002). Trichome initiation takes place only at the base of young leaves (Figure 2, indicated by a red square). In this region, TRY is expressed in all epidermal cells. Expression in trichomes is stronger than in the surrounding epidermis (Schellmann et al, 2002). To corroborate the assumption that the active complex promotes the inhibitor, only those aspects of TRY expression are relevant that occur before or at the time of trichome initiation, i.e. the ubiquitous epidermal expression. This ubiquitous expression is absent in gl1, gl3 and gl3 egl3 mutants, indicating that TRY is transcriptionally induced by the active complex (Figure 2). According to our model, the inhibitors should repress the active complex and therefore indirectly themselves. This is confirmed in plants overexpressing TRY or CPC. No TRY expression is observed in a 35S:TRY or 35S:CPC background (data not shown).
Subcellular localization and ability for intercellular movement of TRY, CPC, GL1 and GL3
In our model, we assume that TRY and GL3 can move into neighbouring cells similar to that shown for CPC (Kurata et al, 2005) and GL3 (Bernhardt et al, 2005) in the root hair system and that GL1 acts cell autonomously (Hulskamp et al, 1994). To test this in leaves, we fused the coding sequence of TRY, CPC, GL1 and GL3 to GFP and placed them under the control of the constitutive CaMV 35S promoter. 35S:GFP and 35S:GFP:YFP constructs served as controls and a 35S:YFP:peroxisome marker was used to label the targeted cells. The functionality of the TRY and CPC fusion proteins was demonstrated by showing that their overexpression results in a loss of trichomes. The GL1 and GL3 fusion proteins have previously been shown to be functional by Esch et al (2003). Arabidopsis cotyledon and leaf epidermal cells were transformed by micro‐projectile bombardment and analysed after 6–10 h. The TRY and CPC fusion proteins were localized in the targeted cell and in approximately one‐third of the neighbouring cells. This demonstrates that these fusion proteins can move from the originally transformed cell to its neighbours (Figure 3; Tables I and II). GFP‐GL1 and surprisingly also GFP‐GL3 proteins did not move in this assay (Figure 3; Table I). The finding that GL3 is cell autonomous is incorporated into our model.
TRY physically interacts with GL1 and GL3
By assessing the interaction of TRY with other patterning genes using the bimolecular fluorescence complementation (BiFC) system and pull‐down experiments, we confirmed the interaction between TRY and GL3 and EGL3 (Kirik et al, 2004b; Bernhardt et al, 2005) previously shown using yeast two‐hybrid experiments (Figure 4). To our surprise, we also found an interaction of TRY with GL1 in both assays. This suggests an additional inhibitory interaction where TRY competes with GL3 for binding to GL1.
Challenging the system by overexpression experiments
The finding that GL1 can bind to TRY as well as the fact that TRY can bind to GL3 raises the possibility for three different inhibition scenarios (Figure 1B). (1) Single competitive inhibition: TRY binds to free GL3 and prevents GL1 GL3 dimerization; (2) Double competitive inhibition: TRY binds to free GL3 or GL1, thereby preventing the interaction between GL1 and GL3; (3) Uncompetitive inhibition: TRY binds to the GL1 GL3 dimer and represses its function. Although the double competitive model contains the single competitive model as an extreme case, the uncompetitive model is independent of the other two inhibition scenarios. In all three cases, we found regular spacing patterns for biologically meaningful parameters, although with different parameter ranges and sensitivities (Supplementary Figure 4). This suggests manipulation of the system experimentally to discriminate the biological relevance of the three scenarios. We decided to consider overexpression of TRY and GL3 either constitutively under the 35S promoter or under the common downstream promoter of GL2. We created the respective plant lines to investigate the results of these experimental interventions. All GL2:TRY plants are completely glabrous in two independent lines. 35S:GL3 lines show a higher trichome density in the patterning zone (compare Figure 5B and C). In the most basal part of the leaf spanning the region containing only unbranched trichomes, a significant difference is found between wild‐type Ler (9.2±2.7, n=30) and 35S:GL3 (14.7±4.3, n=30) by the Student's t‐test (P<0.01). However, in GL2:GL3 lines, we observe no significant difference in trichome density in the patterning zone (wild type: 6.3±1.7, n=30; GL2:GL3: 7±2.77, n=30; compare Figure 5A and D). It is noteworthy that in both GL3 overexpression lines, the mature part of the leaf displays additional young trichomes regularly scattered between fully mature trichomes. This situation is rarely found in wild‐type Columbia (7.3%±7.3, n=20) and Ler (7.8%±8.4, n=20). In GL2:GL3 and 35S:GL3, they represent a significant fraction of all trichomes (47%±7.3, n=21) and (44.1%±11.5, n=22), respectively. This indicates that these young trichomes are formed later in leaf development, suggesting that the time window allowing trichome formation is expanded in GL2:GL3 and 35S:GL3 lines.
We simulated the corresponding overexpression experiments of TRY and GL3 using the three inhibition scenarios. Because the model parameters are unknown, we cannot directly test the accordance of the different simulation results with our experimental findings. However, we can confine the model parameters to biologically meaningful ranges and evaluate the accordance within these ranges (see Table III). We randomly sampled parameters inside the reduced parameter space and tested for each sample whether the following five criteria are fulfilled: (1) pattern formation from homogeneous initial conditions is possible, i.e. the parameter set is inside the Turing space; (2) 35S:TRY overexpression leads to a loss of patterning; (3) GL2:TRY overexpression leads to a loss of patterning; (4) 35S:GL3 overexpression yields a higher trichome density compared to wild type; (5) GL2:GL3 overexpression results in a similar trichome density compared to wild type. For details of the simulations, see Materials and methods.
The frequency of matched criteria for the three competition scenarios is presented on a logarithmic scale in Figure 6. The single competitive inhibition scenario meets each criterion most often, followed by the double competitive scenario. The uncompetitive scenario matches each criterion at least one order of magnitude less compared to the two other cases. In particular, the frequency of randomly hitting the Turing space is only 2%, whereas it is more than 30% in the two other scenarios. Interestingly, we found parameter samples that match all criteria simultaneously only for the single and double competitive scenarios but not for the uncompetitive case. As we could not reproduce all the experimental overexpression phenotypes simultaneously with the uncompetitive inhibition scenario, we conclude that it is not consistent with our data.
Single competitive inhibition is a submodel of double competitive inhibition as it contains only the GL3 TRY interaction, whereas the latter additionally contains the GL1 TRY interaction. Hence, we investigate the relation between the parameter samples matching all criteria of these two models simultaneously. In particular, we ask whether the double competitive inhibition is only able to meet all criteria in the extreme case of single competitive inhibition. We calculate the pairwise correlations between all parameters of the double competitive inhibition model and find a strong negative correlation (ρk3, k7= −0.55, P<0.01) between the binding rates of TRY to GL1 and GL3. Parameter sets with comparable binding affinities are found far less frequent than these two extreme cases (see Supplementary Figure 3). Thus, the double competitive inhibition explains in most cases all experimental data either by a strong binding affinity of TRY to GL3 and a weak binding affinity to GL1 or vice versa. Therefore, we conclude that single competitive inhibition is sufficient to explain all experimental findings and double competitive inhibition is not necessary to describe trichome patterning in the context of our data.
A qualitative comparison between the simulated and experimental overexpression of GL3 in the single competitive inhibition scenario is presented in Figure 5C–G. Note that the irregularities of the simulated trichome pattern are due to the fact that only TRY acts non‐cell autonomously and can also be observed experimentally. In Figure 5H and I, the effect of the GL3 and TRY overexpression is illustrated by a two‐dimensional section of the parameter space. Parameter regions that give rise to a Turing instability, as determined from a linear stability analysis, are indicated in white, whereas regions outside the Turing space are given in grey. For details of the analysis, see the Supplementary information. This local visualization of the Turing space provides an explanation why the overexpression of GL3 and TRY yield different phenotypes. Because of the local convexity of the Turing space, an increased TRY expression yields a parameter set outside the Turing space and consequently leads to a loss of pattern. In contrast, an increased GL3 expression preserves the pattern‐forming ability of the model. Note that we only consider the relative increase in expression rates, as the absolute values of the expression rates in wild type and mutant are unknown. The increased sensitivity towards a manipulation of the TRY level in comparison to the GL3 level holds for constitutive as well as AC‐driven overexpression. It is noteworthy that this increased sensitivity is not a local feature of the Turing space at this particular parameter set given in Figure 5, but hold for very different parameter sets as indicated by the results presented in Figure 6 and Supplementary Figure 2.
The use of phenomenological descriptions in the study of development has a long tradition (Hofmeister, 1868; Thompson, 1942). Essential principles of development were formulated as mathematical models that describe the observed developmental pattern on a phenomenological basis (Gierer and Meinhardt, 1972; Cooke and Zeeman, 1976; Othmer, 1977; Mitchison, 1977; Erneux, 1978; Mitchison, 1980). The overwhelming advances in the molecular understanding of developmental processes in the last decades demand an adaptation of this successful approach to the details of current knowledge and experimental possibilities. Mechanistic in contrast to phenomenological descriptions can fill this gap. Mechanistic models can predict the effect of specific genetic or biochemical interventions and can suggest the design of specific experiments. In this study, we present a mechanistic model for trichome patterning that can be used to address several experimentally feasible interventions of the system.
Our study reveals several new properties of trichome patterning in Arabidopsis leaves. First, we show that GL1 and GL3 are actually positive regulators for the expression of TRY and CPC, which has been assumed previously. This fact is a necessary requirement of theoretical models and has not been experimentally confirmed before. Second, we demonstrate by micro‐projectile bombardment that the inhibitors TRY and CPC can move into neighbouring cells. In contrast to the root system, we find that GL3 is cell autonomous. Both results are in agreement with the requirement of the activator–inhibitor model for a farther transport of the inhibitors compared to the activators. As predicted by theoretical models including only the mobility of the inhibitor, the resulting pattern is less regular than in a model including the mobility of both substances. This is also consistent with the less orderly pattern observed in the initiation zone of the young leaf. Further, we confirm the previously shown interaction of GL3 and TRY and find an additional interaction between TRY and GL1. This opens the possibility for different ways how TRY can exert its inhibitory function. We construct three alternative models that reflect different inhibition scenarios: (1) single competitive inhibition in which TRY binds to free GL3 and blocks the formation of the active complex, (2) double competitive inhibition that additionally includes the binding of TRY to free GL1 and (3) uncompetitive inhibition whereby TRY binds only to the activator complex and renders it inactive. Analysis of the three models suggests experimental interventions in the form of different overexpressions to identify the most relevant of these scenarios. We use the 35S promoter and the promoter of the downstream target GL2 to overexpress GL3 and TRY either constitutively or under the control of the active complex, respectively. For each of the three scenarios, we compare simulation results with the overexpression phenotypes. This comparison depends on the particular choice of parameter values. As these are unknown, we employ a global sampling strategy that allows us to evaluate the different scenarios and identify constraints on the parameters. Using this approach, we can predict that the interaction between TRY and GL3 is the most relevant for trichome patterning. In particular, the uncompetitive inhibition scenario cannot reproduce our experimental findings, whereas the double competition scenario mostly does so in the extreme case of single competition. Interestingly, the single competitive inhibition scenario tolerates larger parameter variations than both other scenarios while still forming the correct patterns (cf. Supplementary Figure 4). In this sense, the model can explain why the trichome patterning system is robust and can yield a similar qualitative output despite naturally occurring perturbations, e.g. in the form of genetic and environmental variability.
Although we have considered only few components of the patterning system, our combined experimental and theoretical effort reveals specific system properties far beyond any intuitive judgement. Yet the model is not sufficient to explain many additional observations, e.g. the single and double mutant phenotypes of the inhibitors. Additional experiments on the specific properties of the inhibitors are needed. With joint endeavour of theory and experiment, it should be possible to uncover the principles of trichome pattern formation.
Materials and methods
From the interactions depicted in Figure 1, we derive the following system of coupled ordinary differential equations describing the time evolution of GL1, GL3, TRY and the active complex (AC).
The model does not include an equation describing the time evolution of any inactive complex as these do not feed back into the system. Note that the square in equation (3) is needed for the stability of the system. For biological relevance, all variables and parameters have to be non‐negative. The model is formulated on a hexagonal grid with coordinates j=(y, x) as depicted in Supplementary Figure 1. Here, 1⩽x⩽xmax and 1⩽y⩽ymax, and xmax and ymax are the number of cells in x and y directions, respectively. We assume periodic boundary conditions of the grid. To describe the coupling between neighbouring cells, we define the passive transport of a variable [C]i as
where C stands for GL3 or TRY. The parameters in our model are rates for basal expression (σi), regulated expression (αi), degradation (ρi), complex formation (βi) and transport (γi) of the corresponding species. Note that the basal expression of TRY is only incorporated to enable simulations of 35S:TRY overexpression experiments (σ3>0i); in all other cases, σ3=0. The complex formation rates, which together determine the type and strength of the inhibition, are adapted to reflect three different scenarios. These are single competitive (β2=0, β3>0, β4=0), double competitive (β2>0, β3>0, β4=0) and uncompetitive (β2=0, β3=0, β4>0) inhibition. As all model parameters are unknown, a rescaling of time and concentration is applied. This approach allows us to reduce the total number of model parameters and to confine the resulting dimensionless parameters to biologically reasonable ranges, while retaining a mathematically equivalent set of equations. All concentrations are rescaled by the factor β1/ρ1. Time is expressed in units of half‐life of GL1, i.e. τ=ρ1t. The resulting dimensionless model has the following form:
Table III lists the dimensionless parameters together with their corresponding functional relation to the dimensional parameters. Note that due to our experimental finding that GL3 is cell autonomous, we set k8=0 in all numerical simulations.
All simulations are performed with MATLAB from Math Works Inc. The ODE system (6)–(9) is integrated using the ode15s function of MATLAB, which is designed to solve stiff ODEs. As we could not obtain an analytical expression for the homogeneous steady state of equations (6), (7), (8) and (9), we determined the steady state by numerical integration of the single‐cell model starting with zero initial conditions for all protein concentrations. This steady state is used with additional 1% random variation per cell as initial conditions for the simulation of the grid of cells, including spatial coupling.
For each inhibition model, 106 random parameter samples are drawn from an exponential distribution to sample the orders of magnitude of each parameter uniformly. Each sample is confined to the ranges specified in Table III. For a given parameter sample, the condition for Turing instability is checked by a linear stability analysis as described in detail in the Supplementary information. Basal GL3 overexpression is simulated by a five‐fold increase in the expression rate k4. As a basal expression of TRY is absent in planta, 35S:TRY overexpression is simulated by setting the basal expression rate of TRY, i.e. k9, equal to the value of k4 in the simulated 35S:GL3 overexpression. The overexpression of GL3 and TRY under the GL2 promoter is simulated by a five‐fold increased rate of active complex‐dependent expression, i.e. rate k5 and k10, respectively, as the active complex binds to the promoter of GL2. The effect of each of these four parameter perturbations is checked by a linear stability analysis (for details of the analysis, see Supplementary information). Additionally, numerical simulation of the 35S:GL2 and GL2:GL3 plants is used to determine the trichome density of the corresponding plants. The simulated 35S:GL3 plant is said to fulfil the experimental observation of an increased trichome density if it is at least 1.5 times higher than the simulated WT density. The density of the simulated GL2:GL3 plants DGL2:GL3 is said to match the WT density DWT if (DGL2:GL3−DWT)2⩽0.0004 holds.
Parameters used for simulation presented in Figure 5
Wild‐type parameter values: k1=8.2707, k2=3.4869, k3=0, k4=15.0952, k5=1.3488, k6=0.4503, k7=7.9509, k8=0, k9=0, k10=0.4117, k11=0.9565, k12=0, k13=10 and k14=0.2703. Parameter values of the overexpressed mutants are chosen as described above.
Plant material, growth condition and genetic methods.
Plants were grown at 22°C for 16 h white light a day. The Arabidopsis lines TRY:GUS, gl1‐1, gl3, gl3 egl3, 35S:TRY and 35S:CPC have been described previously (Oppenheimer et al, 1991; Larkin et al, 1993; Hulskamp et al, 1994; Lloyd et al, 1994; Wada et al, 1997b; Szymanski and Marks, 1998; Walker et al, 1999b; Schellmann et al, 2002). Crosses of the TRY:GUS lines to recessive mutants were performed by using the respective GUS line as male parent. The F2 generation was screened for BASTA resistance.
For the swapping experiment, GL2:GL3 lines have been described previously (Kirik et al, 2005). 35S:GL3 was generated by recombination of full‐length GL3 cDNA into pAMPAT vector. 35S:TRY and GL2:TRY were generated by cloning of TRY cDNA in pCAMBIA 1300 (Mathur et al, 2003). The constructs were transformed into the Agrobacterium strain GV3101 by electroporation (Bio‐Rad gene pulser) and Landsberg erecta plants were transformed using the floral dip method (Clough and Bent, 1998). Transformants were selected in the T1 generation using MS plates containing 50 mg/l kanamycin or on soil with a 0.1% BASTA solution.
Constructs, biolistic transformation and BiFC.
The constructs used for the transient assay in leaf epidermis were created by fusing GFP (Clontech) to the N terminus of the coding sequence of GL1, GL3, TRY and CPC and placed under the control of the CaMV 35S promoter. Details are available on request. Transient expression analysis was carried out by using the particle bombardment method as described previously (Mathur et al, 2003).
TRY, GL1 and GL3 were recombined into BiFC vectors (pBatTL) and transfected in Arabidopsis protoplast (Uhrig et al, 2007). The transfected cells were incubated at 23°C for 16–20 h in the dark before microscopic observation.
Histology and microscopy.
GUS staining was performed as described by Malamy and Benfey (1997). Plant specimens were analysed using the LEICA‐DMRE microscope equipped with a high‐resolution KY‐F70 3‐CCD JVC camera and frame‐grabbing software (DISKUS; Technisches Büro, Königswinter). Photos were edited with Adobe Photoshop 6.
Expression and purification of recombinant fusion proteins.
GL1, GL3, EGL3 and TRY cDNA were recombined into pGEX2TMGW and pDEST 17 vectors by Gateway Cloning (Invitrogen) to create GST and His fusion proteins, respectively. Recombinant proteins were expressed in Escherichia coli BL21DE3RIL cells (Stratagene) by the induction of GST fusion proteins with a final IPTG concentration of 0.1 mM at 37°C for 4 h and of His fusion proteins with a final concentration of IPTG of 1 mM at 37°C for 3 h. Bacteria expressing GST fusion proteins were lysed as described by Frangioni and Neel (1993) and proteins were purified through the glutathione sepharose (GE Healthcare cat. no: 17‐5279‐01) by batch method as described by Sambrook and Russell (2001). The buffer of the purified proteins was changed against 50 mM Tris–HCl pH 7.8, 1 mM EDTA, 150 mM NaCl, 0.1% Nonidet P‐40 and 4% BSA by the Amicon Ultra 10k. Bacteria expressing His fusion proteins were pelleted after induction, and lysis was achieved by sonicating five times using 10‐s pulses in the lysis buffer (100 mM NaCl, 50 mM Tris pH 7.9, 2 mM EDTA, 1% Triton X‐100). His fusion proteins were purified through Ni‐NTA resin (Qiagen), according to the product instructions. Final elutions were done with phosphate‐buffered saline buffer containing a final concentration of 20 mM EDTA. Purified proteins were dialysed against 50 mM Tris–HCl pH 7.8, 1 mM EDTA, 150 mM NaCl. All of the buffers used contained Roche Complete Protease Inhibitor Cocktail™.
GST and His fusion purified proteins (0.5 μg) were incubated together for 2 h at 4°C on a rocking platform. Then 50 μl of glutathione resin (Mathur et al, 2003; GE Healthcare cat. no: 17‐5279‐01) was added and the mixture was incubated further for 2 h at 4°C on a rocking platform. After incubation, the mixture was washed for five times with buffer (50 mM Tris–HCl pH 7.8, 1 mM EDTA, 150 mM NaCl, 0.1% Nonidet P‐40). Pull‐down assays and the mixture was washed either in the presence or absence of BSA. SDS gel extraction buffer was added to the samples. Samples were run on SDS–PAGE and blotted against PVDF membranes. Anti‐His antibodies were used for the detection of the His‐tagged pulled down proteins.
We appreciate the support of Amanda Walker for the completion of some work in her laboratory. This study was supported by the Sonderforschungsbereich 572 of the Deutsche Forschungsgemeinschaft to MH. SD was supported by the Graduate School for Biological Sciences and BD by the European ADOPT program. BG was supported by the FP6 COSBICS Project (512060), FG by BMBF NGFN II 101 35 05 201 and CF by BMBF FRISYS 0313921. pGEX‐2TM‐GW vector is kindly provided by Dr Bekir Uelker and Imre Somssich, Max‐Planck Institute for Plant Breeding Research, Department of Plant Microbe Interactions, Cologne, Germany.
Supplementary Figures 1–4 [msb200854-sup-0001.pdf]
Supplementary Information [msb200854-sup-0002.zip]
This is an open‐access article distributed under the terms of the Creative Commons Attribution License, which permits distribution, and reproduction in any medium, provided the original author and source are credited. This license does not permit commercial exploitation without specific permission.
- Copyright © 2008 EMBO and Nature Publishing Group