Tracking transcription factor mobility and interaction in Arabidopsis roots with fluorescence correlation spectroscopy
Abstract
To understand complex regulatory processes in multicellular organisms, it is critical to be able to quantitatively analyze protein movement and protein-protein interactions in time and space. During Arabidopsis development, the intercellular movement of SHORTROOT (SHR) and subsequent interaction with its downstream target SCARECROW (SCR) control root patterning and cell fate specification. However, quantitative information about the spatio-temporal dynamics of SHR movement and SHR-SCR interaction is currently unavailable. Here, we quantify parameters including SHR mobility, oligomeric state, and association with SCR using a combination of Fluorescent Correlation Spectroscopy (FCS) techniques. We then incorporate these parameters into a mathematical model of SHR and SCR, which shows that SHR reaches a steady state in minutes, while SCR and the SHR-SCR complex reach a steady-state between 18 and 24 hr. Our model reveals the timing of SHR and SCR dynamics and allows us to understand how protein movement and protein-protein stoichiometry contribute to development.
https://doi.org/10.7554/eLife.14770.001eLife digest
Stem cells are a specific type of cell found in both plants and animals. These cells can divide to produce daughter cells that can take on the role of any of the different tissues and organs within the plant or animal. A plant known as Arabidopsis is often used as a model in scientific research. In Arabidopsis, two proteins called SHORTROOT and SCARECROW are known to control the ability of stem cells in the roots to divide.
SHORTROOT is made in cells at the center of the root known as the vasculature. From there, it moves to the next cell layer (called the endodermis) where it interacts with SCARECROW to form a protein complex. Here, Clark et al. investigated how quickly SHORTROOT moves between cells, the direction it moves in, and how it interacts with SCARECROW.
The experiments used a new imaging technique called scanning fluorescence correlation spectroscopy to track the movements of SHORTROOT molecules in the root. This technique relies on the protein of interest (in this case, SHORTROOT) being attached to a fluorescent protein so that it is visible when the cells are examined. In plants that had lower levels of SCARECROW, SHORTROOT moves between cells more quickly and in an unrestricted manner. This suggests that SCARECROW forms a complex with SHORTROOT to restrict its movement in the endodermis. The experiments also show that SHORTROOT is only able to leave the endodermis to return to the vasculature when SCARECROW levels are lower than normal.
Clark et al. developed a model to map the behavior of SHORTROOT and SCARECROW in the root and predict how the levels of these proteins change over time. One of the next steps following on from this work would be to test whether other proteins restrict the movement of SHORTROOT, perhaps by studying mutant plants in which SHORTROOT is less able to move.
https://doi.org/10.7554/eLife.14770.002Introduction
During development, multicellular organisms must coordinate patterning, maintenance, and growth of different cell types. This dynamic coordination is achieved through complex spatio-temporal signaling mechanisms (Voas and Rebay, 2004; Scheres, 2007; Han et al., 2014; Sozzani and Iyer-Pascuzzi, 2014; Sablowski, 2015; Tuazon and Mullins, 2015; Fisher and Sozzani, 2016). Many of the signaling mechanisms regulating development utilize mobile signals that cross cell boundaries (Koizumi et al., 2011; Xu et al., 2011; Gallagher et al., 2004). In plants, the intercellular movement of transcription factors and the spatio-temporal control of protein complex formation regulate many processes including cell fate specification.
The key to a systems-level understanding of development in multicellular organisms is the ability to obtain quantitative information about various signaling factors. In the Arabidopsis root, the SHORTROOT-SCARECROW (SHR-SCR) mediated gene regulatory network (GRN) is a well-characterized developmental pathway that controls ground tissue patterning and endodermal cell fate specification (Levesque et al., 2006; Cui et al., 2007; Sozzani et al., 2010; Cruz-Ramırez et al., 2012; Long et al., 2015; Moreno-Risueno et al., 2015). SHR is transcribed in the vasculature (Helariutta et al., 2000; Nakajima et al., 2001), and then the protein moves to the adjacent cell layer where it is retained in the nuclei of the Quiescent Center (QC), Cortex/Endodermis Initials (CEI) and endodermis (Nakajima et al., 2001; Gallagher et al., 2004; Gallagher and Benfey, 2009). In these cells, SHR activates the expression of the downstream transcription factor SCR (Helariutta et al., 2000; Levesque et al., 2006; Cui et al., 2007; Sozzani et al., 2010), which, as shown by yeast two hybrid experiments, interacts with SHR and prevents further SHR movement (Heidstra et al., 2004; Cui et al., 2007; Long et al., 2015). Although there have been important advances in identifying the essential features that govern the SHR-SCR GRN (Gallagher et al., 2004; Sena et al., 2004; Cui et al., 2007; Gallagher and Benfey, 2009; Sozzani et al., 2010; Cruz-Ramırez et al., 2012), the ability to measure key network parameters that may contribute to patterning and cell fate specification remains a fundamental bottleneck.
New imaging tools that enable parameter quantification and acquisition of in vivo kinetic parameters could provide quantitative information that describes temporal and spatial dynamics of proteins in multicellular organisms. Thus, we explored the possibility of combining scanning Fluorescence Correlation Spectroscopy (scanning FCS) techniques. Unlike more common time correlation FCS techniques, which only use temporal information, scanning FCS techniques utilizes both the spatial and temporal information present in a confocal raster scan to measure protein movement, protein-protein interactions, and the stoichiometry of protein complexes. (Petrasek and Schwille, 2008; Digman and Gratton, 2011). Previously, these techniques have only been used to quantify protein mobility and the dynamics of protein association in cell cultures (Digman et al., 2005a; 2005b; Digman and Gratton, 2009a; Digman et al., 2009a; Jameson et al., 2009; Rossow et al., 2010; Hinde et al., 2010; 2011; Vetri et al., 2011). We combined the techniques of Raster Image Correlation Spectroscopy (RICS), Pair Correlation Function (pCF) and Number and Brightness (N&B) to analyze SHR and SCR mobility and interaction at high spatio-temporal resolution. By using RICS and pCF (Digman et al., 2005a; 2005b; Brown et al., 2008; Digman and Gratton, 2009a; 2009b; Jameson et al., 2009; Rossow et al., 2010; Hinde et al., 2010; 2011; Vetri et al., 2011; Digman and Gratton, 2011), we quantified the rate and directionality of SHR movement. Specifically, we used RICS and a 3D Gaussian diffusion model (Digman et al., 2005a; 2005b) to calculate the diffusion coefficient of SHR in different root cell types. We also acquired line scans and performed pCF analyses (Hinde et al., 2010) to assess the directionality of SHR movement in these different cell types. Moreover, we used N&B and cross-correlation analyses (Digman et al., 2008; 2009a; 2009b) to characterize the oligomeric state of SHR and the stoichiometry of the SHR-SCR complex, respectively. Finally, we incorporated the diffusion coefficient of SHR and the stoichiometry of the SHR-SCR complex into a mathematical model of SHR and SCR dynamics. Our results demonstrate that these parameters can be used in predictive mathematical models, allowing us to understand how protein movement and stoichiometry of protein complexes contribute to developmental processes. Further, our study highlights how these non-invasive scanning FCS techniques can be used to experimentally measure protein movement and protein-protein interactions within multicellular organisms.
Results
Protein movement quantified using raster image correlation spectroscopy (RICS)
A key parameter in biological models is molecular diffusion, which is frequently estimated based on published measurements from single cell organisms (Spiller et al., 2010). In order to measure protein movement in a multicellular organism such as Arabidopsis, we used Raster Image Correlation Spectroscopy (RICS), which returns an autocorrelation function (ACF) by correlating fluorescence intensity fluctuations in pixels in an image over time and space. The diffusion coefficient (DC) is then calculated by fitting the ACF with a diffusion model (Figure 1). Since the fit of the diffusion model depends on the choice of imaging parameters, such as pixel size and pixel dwell time, we first set these parameters by performing RICS analysis on Green Fluorescent Protein (GFP) driven by the CaMV 35S constitutive promoter (35S:GFP) (Table 1, see Materials and methods). In the Arabidopsis root, the resulting DC for free GFP was 6.33 ± 0.37 μm2/s (n = 34 for different cells, including vascular, endodermal, and QC cells) (Figure 1—figure supplement 1). We obtained similar diffusion coefficients using two different confocal microscopes (Zeiss 780 and Zeiss 710).
-
Figure 1—source data 1
- https://doi.org/10.7554/eLife.14770.004
-
Figure 1—source data 2
- https://doi.org/10.7554/eLife.14770.005
-
Figure 1—source data 3
- https://doi.org/10.7554/eLife.14770.006
-
Figure 1—source data 4
- https://doi.org/10.7554/eLife.14770.007
We next used RICS to determine if the DC of SHR differs between the vasculature, where it is produced, and the endodermis and QC, where it moves and interacts with SCR (Heidstra et al., 2004; Cui et al., 2007; Long et al., 2015). To this end, we used a SHR:SHR-GFP translational fusion, which complements the shr2 mutant phenotype (SHR:SHR-GFP in shr2) (Nakajima et al., 2001). Our RICS analysis showed that SHR moves at a rate of 2.45 ± 0.26 μm2/s (n = 40) and 1.73 ± 0.15 µm2/s (n = 16) within a population of vascular and endodermal cells and vascular and QC cells, respectively (Figure 1). Notably, in a population of only endodermal cells, SHR moves at a significantly slower rate (1.29 ± 0.14 μm2/s n = 19, Wilcoxon with Steel-Dwass, p = 0.0303) (Figure 1). To investigate whether the reduction in SHR movement in the endodermis may be due to its interaction SCR (Cui et al., 2007; Levesque et al., 2006; Sozzani et al., 2010; Moreno-Risueno et al., 2015), we measured the DC of SHR in a SCR RNAi (SCRi) line (SHR:SHR-GFP in SCRi) (Cui et al., 2007). In this line, in which the levels of SCR are reduced, SHR diffusion in populations of only endodermal cells was similar to that in populations of vascular and endodermal cells (2.40 ± 0.29 μm2/s, n = 14, Wilcoxon with Steel-Dwass, p = 0.9337) (Figure 1). Overall, our RICS data are in agreement with genetic and molecular data that show SHR movement is restricted by SCR (Heidstra et al., 2004; Cui et al., 2007; Long et al., 2015). Most importantly, our RICS data provide a precise quantification of how SHR diffusion is affected by the presence of SCR.
The direction of protein movement determined using pair correlation function (pCF) analysis
An open question is whether SHR moves unidirectionally, only from the inner to the outer cell layers, or in a bidirectional fashion. The former would be consistent with an active transport mechanism from the vasculature into the endodermis. To detect the route of intercellular movement of SHR between the vasculature and endodermis as well as between the endodermis and cortex, we used pair correlation function (pCF) analysis on line scans (Figure 2). We used pCF to measure the directionality of movement by correlating pixels that are separated by a specific pixel distance (Hinde and Cardarelli, 2011). To account for differences in cell size, as well as cell wall orientation within our images, we used three different distances: 5 pixels, 7 pixels, and 9 pixels. Generally, the pCF analysis returns a carpet, or heatmap, that shows the fluorescence correlation over time (y-axis) and space (x-axis). If proteins move across the cell wall, there is an arch that represents the delay in movement. If proteins are unable to cross the barrier, the arch is absent (Hinde et al., 2010) (Figure 2) (see Materials and methods). Thus, we performed a binary analysis on each carpet to determine the movement of SHR between cells by looking at the presence, or absence, of these arches. Specifically, we recorded a 1 if the carpet showed an arch and a 0 if no arch was present (Figure 2). We took the average of these values from the different pixel distances (5, 7, and 9) to represent one biological replicate. We then calculated the protein Movement Index (MI), which is the average of all biological replicates (Figure 2). As a positive control, we acquired pCF data for 35SGFP, which had a MI = 0.71 ± 0.07 (n = 15). As a negative control we used 3xGFP, which was shown to restrict free GFP movement in roots (Kim et al., 2005). To this end, we drove the 3xGFP using a root vasculature promoter (TMO5:3xGFP, Schlereth et al., 2010) which had a MI = 0.26 ± 0.05 (n = 19) (Figure 2—figure supplement 1).
-
Figure 2—source data 1
- https://doi.org/10.7554/eLife.14770.011
-
Figure 2—source data 2
- https://doi.org/10.7554/eLife.14770.012
We then used pCF to measure SHR movement across the vasculature, endodermis, and cortex. We found that SHR moves from the vasculature to the endodermis (MI = 0.58 ± 0.03, n = 20, Wilcoxon with Steel-Dwass, p = 0.6356), consistent with previous experimental data (Helariutta et al., 2000; Nakajima et al., 2001) (Figure 2). However, from the endodermis back to the vasculature, the MI of SHR was significantly lower than the MI of free GFP (MI = 0.36 ± 0.07, n = 20, Wilcoxon with Steel-Dwass, p = 0.0368) and not significantly different from the MI of the 3xGFP (Wilcoxon rank-sum test, p = 0.5974), suggesting that SHR is unable to move back to the vasculature (Figure 2). In addition, we found that SHR moves bidirectionally between the endodermis and cortex, as the MI is not significantly different from that of free GFP (MI = 0.61 ± 0.06 from endodermis to cortex, MI = 0.53 ± 0.06 from cortex to endodermis, n = 20, Wilcoxon with Steel-Dwass, p = 0.9524 and p = 0.3909 respectively).
Given that our RICS analysis showed that the presence of SCR reduces SHR mobility in the endodermis, we measured SHR intercellular movement with pCF in the absence of SCR. Accordingly, we performed a pCF analysis using SHR:SHR-GFP in the SCRi line. Line scan measurements were taken across vascular, endodermal, and cortical cells as well as across the extra layer between the endodermis and cortex which is a direct consequence of the RNAi reduction of SCR (Cui et al., 2007). We not only observed SHR movement from the vasculature to the endodermis with a MI = 0.66 ± 0.07 (n = 14) but also from the endodermis back to the vasculature with a MI = 0.78 ± 0.09 (n = 14) (Figure 2). Since the MI for SHR:SHR-GFP in SCRi in both directions is not significantly different from the MI for 35S:GFP (Wilcoxon rank-sum test, p = 0.9553 for endo to vasc, p = 0.9999 for vasc to endo) it suggests that movement is unrestricted in both directions (Figure 2). Similarly, we found that SHR is still able to move bidirectionally between the endodermis and cortex (MI = 0.61 ± 0.05 from endodermis to cortex, MI = 0.59 ± 0.07 from cortex to endodermis, n = 17, Wilcoxon with Steel-Dwass, p = 0.9997 and p = 0.5510 respectively). These results suggest that the presence of SCR prevents SHR movement specifically from the endodermis back to the vasculature while it does not affect the movement between the endodermis and cortex. To determine if SCR can move with SHR from the endodermis to the cortex, we used pCF on the SCR:SCR-GFP. Our results show that SCR does not move between the endodermis and the cortex, as the MI in either direction is significantly lower than that of 35S:GFP (MI = 0.31 ± 0.08 from endodermis to cortex, MI = 0.31 ± 0.07 from cortex to endodermis, n = 14, Wilcoxon with Steel-Dwass, p = 0.0236 and p = 0.0054 respectively). Taken together, our pCF results provide information about the directionality of SHR movement and indicate that SCR restricts SHR movement from the endodermis to the vasculature.
Protein oligomeric state determined by number and brightness (N&B) analysis
Stoichiometry is an important feature of protein complexes, as some transcription factors must form higher order complexes in order to function (Nakashima et al., 2012; Sornaraj et al., 2016). We used the Number and Brightness technique (N and B), which relies on the RICS image acquisition, to investigate the oligomeric state of the SHR protein in different root cell types. We used the average fluorescence intensity, and the variance in fluorescence, to determine the brightness of particles and their number in an image (Digman et al., 2008; 2009a) (Figure 3, see Materials and methods). In order to use N&B to measure SHR oligomeric state, we first obtained the brightness of the autofluorescence (immobile fraction) and of monomeric GFP protein. We used roots expressing 35SGFP to calculate the S-factor, an imaging parameter that shifts the brightness of the image, such that the immobile fraction has a brightness value of 1 (Digman et al., 2008). In addition, we used the 35S:GFP line to measure the brightness of monomeric GFP protein (Figure 3 and Table 2) (see Materials and methods).
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.14770.015
-
Figure 3—source data 2
- https://doi.org/10.7554/eLife.14770.016
-
Figure 3—source data 3
- https://doi.org/10.7554/eLife.14770.017
-
Table 2—source data 1
- https://doi.org/10.7554/eLife.14770.019
-
Table 2—source data 2
- https://doi.org/10.7554/eLife.14770.020
-
Table 2—source data 3
- https://doi.org/10.7554/eLife.14770.021
After calibration, we used N&B on SHR:SHR-GFP to determine the oligomeric state of SHR in different cell types. Our results indicated that SHR exists primarily as a monomer in the vasculature, given that SHR homodimer formation was only 2.6% ± 0.3% in these cells (Figure 3). Notably, SHR monomeric proteins are observed in both the nuclei and in the cytoplasm of vasculature cells. By contrast, in the nuclei of the endodermis, there are significantly more SHR homodimers present (7.5% ± 0.9% homodimer, Wilcoxon with Steel-Dwass, p<0.0001) (Figure 3). Moreover, in the endodermis SHR is present at very low levels outside of the nuclei (Figure 3). To further understand the differences in SHR oligomeric state, we determined if the homodimer of SHR present in the endodermis could be affected by the presence of SCR. Therefore, we performed N&B on SHR:SHR-GFP in the SCRi line and observed that SHR is mainly found as a monomer in the endodermis (1.9% ± 0.4% homodimer, Wilcoxon with Steel-Dwass, p = 0.2546), indicating that SCR influences the oligomeric composition of SHR (Figure 3). These N&B results provide a quantitative assessment of SHR’s oligomeric states and their distribution in the vasculature and endodermis.
Stoichiometry of the SHR-SCR complex
The N&B analysis revealed that SHR exists both as a monomer and as a homodimer in the endodermis. We next asked if both the monomer and homodimer are able to form a complex with SCR. To test this hypothesis, we performed cross-N&B, which can determine the stoichiometry of the SHR-SCR complex. Cross N&B requires that each protein be tagged with a different fluorophore. The analysis then determines which proteins are in a complex by calculating the cross-correlation between the two channels at each pixel (Figure 4, see Materials and methods). Accordingly, we generated a transgenic line containing both SHR and SCR tagged with different fluorophores (SHR:SHR-GFP & SCR:SCR-mCherry) (Figure 4—figure supplement 1). We reasoned that SCR may also exist in higher oligomeric states, which would increase the possible binding ratio of the SHR-SCR complex; therefore we first used N&B on SCR:SCR-mCherry to determine the oligomeric state of SCR. Initially, we determined the S-factor and brightness of monomeric mCherry protein in the root (UBQ10:mCherry) as we did for the monomeric 35SGFP protein (Figure 4—figure supplement 2 and Table 2). In the SCR:SCR-mCherry line, we detected mostly monomers with 4.7% ± 0.5% of homodimers (Figure 4). We tested the oligomeric state of SCR using a SCR:SCR-GFP fusion protein (Figure 4—figure supplement 2), as it was shown that the type of fluorescent tag can change the behavior of a protein and its aggregation (Brown et al., 2008). We found that the SCR:SCR-GFP line had 5.4% ± 0.7% of homodimers, which is similar to the SCR:SCR-mCherry line (Figure 4—source data 1). Therefore, when performing the cross-N&B analysis, we considered the possibility that the homodimer of SCR could be part of the SHR-SCR complex.
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.14770.023
-
Figure 4—source data 2
- https://doi.org/10.7554/eLife.14770.024
After analyzing the SHR:SHR-GFP and SCR:SCR-mCherry separately, we performed cross-N&B on the SHR:SHR-GFP/SCR:SCR-mCherry line to determine the stoichiometry of the complex. The cross-N&B analysis returns a stoichiometry diagram that represents the proportion of different complexes (Figure 4). We found that both the monomer and the homodimer of SHR bind the monomer of SCR, suggesting that SHR and SCR bind with a 1:1 and 2:1 stoichiometry (Figure 4). We were not able to detect any complexes that contain the homodimer of SCR. Additionally, we determined that 84.8% ± 1.6% of the SHR-SCR complexes have 1:1 stoichiometry while 15.2% ± 1.6% have 2:1 stoichiometry (Figure 4). Our cross N&B results reveal that both the monomer and homodimer of SHR are able to bind the monomer of SCR, while the homodimer of SCR does not seem to be part of this complex.
A mathematical model of SHR and SCR dynamics
SHR and SCR dynamics have been previously modeled in the endodermis, but this model did not take into account either the intercellular movement of SHR or the stoichiometry of the SHR-SCR complex (Cruz-Ramırez et al., 2012). After experimentally determining the rate of SHR movement, SHR oligomeric state, and the binding ratio of the SHR-SCR complex, we sought to incorporate this information into a mathematical model of SHR and SCR, with the goal of determining how protein movement and stoichiometry affect SHR and SCR dynamics. Therefore, we constructed a compartmental model that only measures SHR, SCR, and the SHR-SCR complex in both the vasculature and the endodermis (see Materials and methods).
After constructing the model, we performed a sensitivity analysis to determine the most influential parameters. We chose to use Sobol indices to measure how sensitive the model is to each parameter (Sobol, 2001; see Materials and methods). Briefly, the total Sobol effects index measures how much the model outcome varies as the parameters are changed. If small changes in the parameter values cause large changes in the model outcome, then that parameter is more influential. Using this measure, we found that the rate of movement of SHR from the vasculature to the endodermis (a1) and from the endodermis to the vasculature (a2) are both highly influential parameters (Figure 5—figure supplement 1). Thus, this suggests that SHR movement is a key component of the model that greatly influences the dynamics of SHR in the endodermis.
Next, we sought to use the model, in conjunction with our experimentally determined diffusion coefficients, to (i) simulate SHR and SCR dynamics in the endodermis, and (ii) estimate values for the other parameters. Most of the parameter values were chosen based on the previous mathematical model (Cruz-Ramírez et al., 2012). However, since this model did not account for different oligomeric states of SHR we estimated some of the parameters using the N&B data (Supplementary file 3; see Materials and methods). Our mathematical model predicts that SHR reaches a steady state in the vasculature and endodermis in a matter of minutes. The levels of SCR increase greatly in the first 3 hr, which is supported by data that show that SCR expression in a SHR inducible system is significant after 3 hr (Sozzani et al, 2010). While the 1:1 SHR-SCR complex increases greatly in the first 3 hr, the 2:1 complex does not form until after 9 hr. This is because the SHR homodimer does not form until about 9 hr into the simulation (Figure 5). We reasoned that this is a plausible scenario because SCR should exceed 60% of steady-state levels to trigger homodimer formation (see Materials and methods). Finally, the entire system reaches a steady-state between 18 and 24 hr (Figure 5). This suggests that cell division occurs once SCR and the SHR-SCR complexes reach their steady state values.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.14770.028
-
Figure 5—source data 2
- https://doi.org/10.7554/eLife.14770.029
In addition to simulating what happens in wild type, we also wanted to observe how decreasing SCR affects the model dynamics and how it reflects our pCF and N&B data. Accordingly, in the SCRi model, SHR movement is bidirectional and no SHR homodimer forms (see Materials and methods). In this model, we assume that SCR levels are maintained at below 60% of wild type levels as in agreement with previous experimental data (Cui et al, 2007). At these lower levels of SCR, we observe that the steady-state value of the 1:1 SHR-SCR complex is reduced by half (Figure 5). In addition, we observe that SHR homodimers do not form as seen in the N&B analysis (Figure 3). Since SHR homodimers do not form, the 2:1 SHR-SCR complex does not exist (Figure 5). These results suggest that the extra layer in the SCRi line could form due to a reduction in the levels of the SHR-SCR complex. However, there are still other unknown factors that could contribute to the formation of the extra layer.
Discussion
The methodologies we describe are rapid, direct, and convenient approaches for characterizing the quantitative and qualitative behavior of proteins in vivo. Combining RICS and pCF provided detailed understanding of SHR movement both within and between cells of the Arabidopsis root. Trafficking along the radial axis is normally unidirectional - SHR moves from the vasculature to the endodermis, and not in the opposite direction (Figure 2). In the absence of SCR, SHR trafficking is bidirectional. In agreement with the pCF results, our RICS data showed that the DC of SHR in the endodermis increases significantly in the absence of SCR as compared to wild type, indicating that removal of SCR affects the diffusive behavior of SHR. Taken together, these results provide further evidence that SCR, which is present only in the endodermis, spatially restricts SHR movement. The mechanisms underlying this restricted movement remain unclear, and future experiments should focus on uncovering other factors, in addition to SCR, that can bind SHR (Long et al., 2015; Moreno-Risueno et al., 2015). In this study, combining RICS and pCF provided quantitative information about the speed and direction of protein movement.
Using RICS and pCF, which allow both intra- and inter-cellular analysis of protein movement, we were able to determine the in vivo flow of molecules in a multicellular organism. The average value of the diffusion coefficient obtained with RICS could be obtained with other techniques such as single-point FCS or Fluorescence Recovery After Photobleaching (FRAP), but these techniques have limitations, such as returning only temporal information (Miyawaki, 2011). On the other hand, both RICS and pCF are applied to frames or lines, respectively, which increase the spatial resolution and the statistical power of the analysis by simultaneously measuring the movement of many individual molecules in time and space. Moreover, RICS and pCF introduce a spatial component that can determine whether slow movement is due to a binding interaction or slow diffusion, which FRAP and single-point FCS cannot do. Thus, these scanning FCS methods are robust against cell movement and other artifacts that could bias measurements taken with other single-point methods.
The DC measurements of SHR in both the shr2 complemented line and SCRi are significantly lower than those of freely-moving 35S:GFP (Wilcoxon rank-sum test, p<0.0001) (Figure 1). The slower movement of SHR could be due to its size alone. However, since the DC of SHR is three-fold smaller than the DC of 35S:GFP, and the DC is inversely proportional to the cube root of the molecular weight of the protein (Young et al., 1980), the change in DC is not attributable to the molecular weight of SHR (~60 kDa). If it were, then the molecular weight of SHR would have to be 27-fold higher than that of GFP, or approximately 729 kDa (Prendergast and Manni, 1978). Further, the decrease in the DC of SHR relative to 35S:GFP, in addition to the pCF analysis which shows unidirectional movement, suggests an active regulation of SHR movement, which is in agreement with previous data (Gallagher et al., 2004; Sena et al., 2004; Gallagher and Benfey, 2009). Future work could expand our understanding of SHR movement by examining it in callose synthase gain-of-function mutants, which block transport via the plasmodesmata (Vaten et al., 2011), or in Shortroot interacting embryonic lethal (SIEL) mutant alleles, which have reduced SHR movement (Koizumi et al, 2011).
The N&B results provided a quantitative assessment of SHR’s oligomeric states and their distribution across cell types. SHR exists in the vasculature primarily as a monomer and in the endodermis as a monomer and a homodimer (Figure 3). The presence of the SHR homodimer in the nucleus and loss of the homodimer in the SCRi line indicate that SCR is required to maintain SHR in its higher oligomeric forms. Taken together, these results provide insights into the molecular mechanisms by which SCR regulates SHR movement. They also raise new questions, such as how SCR maintains SHR in multimeric forms, and how formation of higher oligomeric complexes helps restrict SHR movement.
By using our experimentally determined parameters in a mathematical model, we were able to determine that SHR reaches a steady state in a matter of minutes, while SCR and the SHR-SCR complex stabilize within 24 hr. This suggests that future experiments that aim to understand the details of the SHR-SCR GRN should focus on a time scale of hours to measure its dynamics. In addition, the sensitivity analysis revealed that diffusion is one of the most important parameters in our model, motivating the need to experimentally measure the diffusion coefficient. In the simulation of SCRi, our model showed that the 2:1 complex does not form and the 1:1 complex is reduced to half of normal levels. This suggests that the mutant layer in the SCRi line is likely due to a reduction in the SHR-SCR complex. These results only scratch the surface of what is likely to be a complex network regulating the spatial localization of SHR as a mobile transcription factor and developmental regulator. Furthermore, they highlight the importance of physical interactions between transcription factors as a regulatory component of transcription factor intercellular trafficking networks. The three scanning FCS methodologies utilized here (pCF, RICS, and N&B) gave consistent results providing strong evidence for their reliability. Our data suggest that the application of in vivo molecule tracking techniques is virtually limitless, which opens exciting new opportunities in all fields of biology.
Materials and methods
Plant material and growth conditions
Request a detailed protocolPrior to plating, Arabidopsis seeds were dry sterilized using 50% bleach and 1.5 ml of HCl for at least 1 hr, imbibed with 500–700 μL of sterile water, and vernalized for 2 days at 4°C in complete dark. After vernalization, seeds were plated on 1x MS (Murashige and Skooge) media supplemented with 1% sucrose and grown vertically at 22°C in long day conditions (16 hr light/8 hr dark). Seedlings were 5 days old when imaged. The 35SGFP, SHR:SHR-GFP in shr2, SCR:SCR-GFP in scr4, SHR:SHR-GFP in SCRi, and TMO5:3xGFP lines are described in (Ruiz et al., 1998; Nakajima et al., 2001; Cui et al., 2007; Sabatini et al., 2003; Schlereth et al, 2010). For the SCR:SCR-mCherry construct, the coding sequence (CDS) of the mCherry fluorescent protein (Goedhart et al., 2007) was amplified using primers with flanking attB sites: mCherry-R2R3 F: 5’-GGGGACAGCTTTCTTGTACAAAGTGGCTATGGTGAGCAAGGGCGAGGAG-3’ and mCherry -R2R3 R: 5’-GGGGACAACTTTGTATAATAAAGTTGCTTACTCACTTGTACAGCTCGTCCATGCC-3 and recombined into pGEMTeasyR2R3 vector by Gateway BP reaction. SCR coding sequence in pDONR221-derived entry clones was previously described (Welch et al., 2007). The root expression vector was created using endogenous SCR promoter in pH7m34GW binary vectors by multiple Gateway LR reactions as described (Long et al., 2015).
Statistical analysis
Request a detailed protocolWe performed all statistical analyses using the Wilcoxon rank-sum test with Steel-Dwass for multiple comparisons at significance level α = 0.05. We chose the Wilcoxon test as not all of our data are normally distributed (Shapiro-Wilk Goodness of Fit test, p<0.0001, Supplementary file 1). In addition, the Wilcoxon test adjusts for our sample sizes, which were unequal between groups. The Steel-Dwass test is used after the Wilcoxon test to perform pairwise comparisons. All statistical analyses were performed using JMP software (http://www.jmp.com).
Point spread function (PSF) measurement
Request a detailed protocolThe Point Spread Function (PSF) measures the radius of the laser beam and is experimentally measured in order to perform Raster Image Correlation Spectroscopy (RICS) analysis (Rigler et al., 1993). The PSF was measured for each objective and each confocal microscope used for image acquisition. To calibrate the PSF for our objective lens, we performed RICS on free EGFP in an aqueous solution. While we determined the PSF using free EGFP in solution, it is possible that the PSF varies depending on the cell or tissue we are imaging. However, the exact value of the PSF is not as important for these scanning FCS techniques as it is for single-point FCS since the scanning techniques measure the time a molecule takes to go from one pixel to the next (Digman et al., 2005b; Petrasek and Schwille, 2008). Thus we made the assumption that the PSF determined using the EGFP solution is an accurate estimation of the true PSF. Accordingly, raster images of a solution of 0.6 µM EGFP were acquired using commercial CLSMs, including the Zeiss 780 and 710 (Zeiss Inc, Germany). We determined the PSF by fitting the autocorrelation function (ACF) to the intensity fluctuations of free EGFP in aqueous solution, obtained from the raster images, while fixing the known EGFP diffusion coefficient of 78 μm2/s (Chen et al., 2002). The ACF then returned the experimental PSF beam waist (Supplementary file 2). We found that the PSF measurement was different for different objectives (Supplementary file 2).
Raster image correlation spectroscopy (RICS)
Request a detailed protocolRICS is a technique that has been developed to determine the rate of movement (either passive diffusion or regulated movement) of fluorescently-labeled particles in a small volume (i.e femtoliter volume of the PSF). In particular, RICS can be used to determine the rate of movement of a GFP-labeled protein, such as SHR:SHR-GFP. During a RICS measurement, the laser performs a raster scan across a selected region of interest. The raster scan involves scanning from left to right across a set number of pixels and then retracing, without backtracking across the already imaged pixels, to perform the next line below. The scanning is repeated for multiple consecutive lines until a 256x256 pixel frame is created, and then starts over again to obtain between 50 to 100 time points. Because each pixel is collected at a different time, and we know where each pixel is located, there is spatial and temporal information included in each individual image (Figure 1). This spatio-temporal information can be used to extract the diffusion coefficient of a population of molecules using the RICS algorithm, which has two steps: (1) background subtraction and (2) image correlation. The background subtraction removes stationary and slow-moving objects so that the image autocorrelation only detects the dynamics of diffusing species (Digman et al., 2005a) while the image correlation compares each pixel to its adjacent pixel in both the x and y direction.
For the RICS analysis the pixel size, pixel dwell time, line scan time, number of frames, and image size were set prior to acquisition. The parameters were set so that each pixel is sampled at a slightly faster rate than the particles move in the solution (Digman et al., 2005a). The laser intensity was set so that the signal to noise ratio is high, but within a range that does not cause photobleaching (Table 2). Since all of these parameters affect the ACF fit and, therefore, the estimation of the diffusion coefficient, we determined the optimal parameters for the A. thaliana root system by testing different values while imaging the 35S:GFP genotype. We chose parameter values that resulted in an ACF with a good fit to the data and low residuals (Table 1). We found that these parameter values remained the same between different confocal microscopes. Moreover, since sample movement can create artifacts in the ACF, which results in an erroneous diffusion coefficient (Digman et al., 2005a), no data were analyzed that had significant sample movement (i.e. sample movement that shifts the region of interest outside of the imaging frame).
The RICS-ACF is decomposed into two correlation functions that depend on ξ (the spatial lag in x) and ψ (the spatial lag in y). The first correlation function, S(ξ,ψ), calculates the spatio-temporal correlation due to the scanning of the microscope. The second correlation function, G(ξ,ψ), calculates the spatio-temporal correlation due to particles diffusing in the medium. The ACF, GS(ξ,ψ), takes both of these correlations into account by multiplying them: GS(ξ,ψ)= S(ξ,ψ)* G(ξ,ψ). The functions are constructed assuming that the distribution of fluorescence intensities follows a 3D Gaussian distribution. The decomposition of the ACF into two parts allows RICS to distinguish random, Brownian motion from diffusing particles in the medium (Digman et al., 2005a).
It is possible for artifacts to occur in the ACF due to slow, mobile structures within the volume. To eliminate these artifacts from our analysis, a moving average with a 10 frame time window was applied to the time series images to subtract the immobile fraction (Digman et al., 2005a). The first and last 5 frames of the time series were ignored since in this time window there is not enough information to calculate the moving average. For the remainder of the frames, the average intensity of the 5 preceding and 5 following frames was subtracted. For example, the average of frames 10–20 was subtracted from frame 15.
The ACF was calculated using the PSF of the beam (Supplementary file 2) and the other parameters that were set before imaging (Table 1). The only unknown parameters in the ACF function were DC, the diffusion coefficient of the particles, and G(0), which is inversely proportional to the number of molecules present. The diffusion coefficient returned by the software was the value that best fits the data (see SimFCS Software Analysis).
Pair correlation function
Request a detailed protocolPair correlation function (pCF) is a technique that allows us to measure the movement of a protein along a line. To do this, we scanned a 32 pixel line through a region of interest in the root. The region of interest was chosen such that 2–3 cells are contained within the frame. To visualize cell walls, we used propidium iodide (PI), which causes the walls to fluoresce red when excited with the 488 nm emission line of an argon laser. In our images, the line scans across the middle of the frame, and the sample was positioned such that the line does not scan directly over a cell wall. Additionally, the line was placed towards the middle of the cell so that movement in both the nucleus and cytoplasm can be measured. The imaging area was cropped such that the edges of the line overlap exactly with the outer cell walls. Once the imaging area was cropped, a reference image was taken to check for image movement. The selected 32 pixel line was then rapidly scanned 2 × 105 times at a pixel dwell time of 12.61 µs. Note that the pixel size was not set as a constant, but rather changed for each image depending on the position of the line scan and the size of the cells analyzed. As soon as the line scan finished, another reference image was taken to check if the root had moved during imaging. If the reference images suggested that the root had moved during imaging, then that line scan was not used for analysis.
The pair correlation function (pCF) for two points at a distance δr as a function of the delay time τ is calculated using Equation (1):
where F(t,0) is the fluctuation in fluorescence intensity at pixel 0 and F(t+τ,δr) is the fluctuation in fluorescence intensity at some other pixel position (δr) at different time delays (τ) (Hinde et al., 2010; Hinde et al., 2011; Digman and Gratton, 2009b). The result of the pCF analysis is a carpet, or heatmap, that displays the correlation in fluorescence over time (y-axis) and space (x-axis). Molecules that move across a barrier display a characteristic ‘arch’ pattern in the pCF carpet output, whereas molecules that do not move across the barrier do not (Hinde et al., 2010) (Figure 2). Since the delay time recovered by the pair correlation function analysis is variable, we performed a binary analysis on the pCF carpets to look for movement (presence of an arch pattern) or no movement (no arch) (see SimFCS Software Analysis).
Number and brightness (N&B)
Request a detailed protocolN&B is used to determine the number (N) and brightness (B) of particles in a volume, which allows us to determine the amount of aggregation of particles. This is useful in determining the oligomeric state of GFP- and mCherry-labeled proteins such as SHR:SHR-GFP, SCR:SCR-GFP, and SCR:SCR-mCherry. A time course of raster-scanned images was obtained using the confocal microscope (see RICS). Certain imaging parameters had to be determined so that the pixels were not under or over sampled (see RICS). In addition, lines expressing monomeric forms of GFP (35S:GFP) and mCherry (UBQ10:mCherry) were used to set the background fluorescence and to measure monomer brightness (see SimFCS Software Analysis).
We obtained both the mean and the variance of the intensity distribution at each pixel in order to determine the number (N) and brightness (B) of the particles. The mean, <k>, and the variance, σ2, of the intensity distribution are given by Equations 2 and 3:
where K is the number of time points and ki is the fluorescence counts for time point i. The number and brightness of the particles can be determined from the mean and variance of the intensity distribution alone due to the assumption that the occupation of particles follows a Poisson distribution (Digman et al., 2008). Using moment analysis (Qian and Elson, 2000), the apparent number (N) and apparent brightness (B) of the particles are defined in Equations 4 and 5.
Note that if the average intensity is fixed and the variance increases, B increases but N decreases. This is because fewer, larger particles cause greater intensity fluctuations as the laser scans than many small particles. The true number of particles, n, and the true brightness, ε, can be calculated from N and B respectively. However, the apparent brightness B is used for the software analysis (see SimFCS Software Analysis).
Cross number and brightness (Cross N&B)
Request a detailed protocolCross N&B follows the same theory as N&B but involves two particles that are marked with different fluorescent proteins (Digman et al., 2009b). Thus, Cross N&B is used to determine the binding ratio of a protein-protein interaction. We specifically use Cross N&B to look at the stoichiometry of the SHR-SCR complex using the SHR:SHR-GFP/SCR:SCR-mCherry line. The apparent number of particles, N, and the apparent brightness, B, were calculated for the green and the red channels separately (see N&B). The cross-variance, , is defined in Equation 6:
where Gi and Ri are the pixel intensities in the green and red channels, respectively, at time i, and <G> and <R> are the mean intensities of the green and red channels. When the cross-variance is zero, the fluctuations in the two channels are independent. When the cross-variance is positive or negative, the two channels are correlated or anti-correlated, respectively. The cross-brightness, BCC, and the cross-number, NCC, are defined in Equations 7 and 8.
To determine the stoichiometry of a protein-protein complex, the brightness of each protein was compared to the cross-brightness at each pixel (Digman et al., 2009b). A large, positive cross-brightness indicates that the two proteins bind at that pixel in the image. The brightness of each protein in the complex determines the stoichiometry (see SimFCS Software Analysis).
SimFCS software analysis
Request a detailed protocolThe SimFCS Software (Digman et al., 2005a), developed at the Laboratory for Fluorescence Dynamics (www.lfd.uci.edu), is used to perform RICS, pCF, N&B, and Cross N&B analysis on raster or line scans obtained using a confocal microscope. For the RICS analysis, the software can reduce the region of interest (ROI) from 256x256 pixels to 128 × 128 pixels in order to obtain a more enriched cell population. For example, we used this feature to obtain QC-enriched populations (Figure 1). After selecting the ROI to use, the software uses the moving average (see RICS) to eliminate any artifacts from immobile fractions. Then, the software fits the RICS-ACF using the imaging parameters provided (Table 1) and returns the diffusion coefficient of the protein. The diffusion coefficient returned results in the ACF curve that best fits the data (Figure 1). Goodness of fit was determined by comparing the residuals to the amplitude of the ACF. We only kept images where the maximum value of the ACF was three fold larger than the greatest residual in order for the RICS analysis to be reliable (Figure 1). Images that had residuals larger than this threshold generally had low laser intensity or sample movement were not used for analysis.
For pCF, the line scan file was loaded as a 32 pixel by 2 × 105 pixel image where the x axis represents position along the line and the y axis represents time. When there was photobleaching in the sample, we eliminated some of the frames acquired at later time points. A period average of 800 frames was then applied so that trends in the fluorescence carpet are easier to see. The carpet is displayed as a gradient, with red corresponding to high correlation and blue corresponding to low correlation. Next, the autocorrelation of each of the 32 pixels in the line scan was calculated using a moving average of 200 frames. The resulting image was a 32 column image where each of the columns represented the autocorrelation value of that pixel. At this step, the 32 columns were aligned with the reference image to determine the pixel location of the cell wall. The location of the cell wall was then used as the column distance ( in order to calculate the pCF. The software calculated the pCF by correlating pixels that are apart, moving from left to right. The pCF in the opposite direction was then calculated by correlating pixels from right to left, instead of left to right. To account for the fact that the cell wall is not straight, the pCF was calculated using pixel distances of 5, 7, and 9. Given that there is heterogeneity in the cell size and cell wall orientation in the root, we would not capture differences by using only one pixel distance. Finally, the color scale of the pCF was adjusted such that high correlation is represented as red, low correlation is represented as blue, and no correlation is black.
N&B analysis was applied to the same data set as RICS. First, we calibrated the software using a monomeric form of the fluorescent proteins, namely, 35S:GFP and UBQ10:mCherry. Images were taken of A. thaliana roots expressing 35S:GFP and UBQ10:mCherry that contained and did not contain the background (background refers to a region of the image that does not contain the root) using the same experimental settings as the RICS analysis (Figure 3). RICS analysis was run on the images to ensure a good ACF fit. Then, N&B analysis was first run on the background images. The software plots the brightness versus the intensity of each pixel. Since the intensity is not used in the N&B calculation, the exact values of the intensity for each fluorophore do not matter. Background images have two distinct populations representing the monomer and the background. The background brightness was standardized by setting the S-factor parameter such that the background population is centered at B = 1 (Figure 3). This ensured that the brightness of the monomer was calibrated for the detector output. Since background brightness can vary between microscopes and laser or detector settings, the S-factor was calculated for each microscope and experimental set-up (Table 2).
Once the S-factor was set, N&B analysis was run on the 35S:GFP and UBQ10:mCherry images with no background. All particles with B > 1 were bounded by a rectangle, or cursor, on the brightness graph, and the position and size of the cursor on the B axis were recorded (Table 2, Figure 3). The size of the cursor measures the distribution of the monomer in the image. The position on the B axis represents the brightness of the GFP or mCherry monomer. The x axis represents the intensity and is not used in the N&B analysis: thus, the intensity axis can vary between images. Since these parameters can vary by microscope, the cursor size and monomer brightness were calculated for each microscope (Table 2).
The quality of acquired data for N&B analysis of GFP- and mCherry-labeled proteins was first determined by RICS analysis of the acquired images to ensure a good ACF fit. The S-factor was set to the value determined by the background image analysis (Table 2). A cursor size was selected that took into account the entire distribution of monomer detected in the brightness histogram. Another cursor of the same size was then positioned at the brightness value that corresponds to twice the brightness of the monomer, as any pixels inside the higher rectangle represent a homodimer of the protein (Figure 3). The percentage of monomer and homodimer were then calculated by dividing the number of pixels inside the monomer, or homodimer, cursor by the total number of fluorescent pixels (monomer plus homodimer).
For Cross N&B using images that contain both GFP- and mCherry-labeled protein, the B1-B2 channel was used to determine the pixels in the green channel that were positively correlated with the red channel. The cursor was positioned in the area of the B1-B2 channel that corresponded to where the GFP and mCherry monomers are located. The cursor was then expanded to include any higher oligomeric states that are present, and that area was set as the correlation area (Figure 4). Once the correlation area was set, the Brightness cross correlation (Bcc) channel was used to determine the brightness of the green and red channels at each of the correlated pixels (Digman et al, 2009b). The software highlighted the mCherry pixels that correlated with the GFP pixels using a green curve. The percentage of each complex stoichiometry (1:1, 2:1, etc) was then calculated by overlaying cursors on GFP monomer, mCherry monomer, GFP homodimer, and mCherry homodimer. As in N&B analysis, the percentage was computed by dividing the pixels in the 1:1 complex by the total number of fluorescent pixels (1:1 plus 2:1). Finally, the software returned a stoichiometry plot that displayed the most likely stoichiometry of the protein-protein complex (Figure 4).
Mathematical model formulation
Request a detailed protocolWe constructed a mathematical model that incorporated our experimentally determined parameters and simulated SHR and SCR dynamics in the root. Our model assumed that transcription and translation happens quickly. Because of this, we modeled transcription and protein movement terms in the same equation. Additionally, we assumed that all proteins have linear degradation terms.
First, we developed a model reflecting wild type conditions (Model 1). We modeled six different variables: SHR in the vasculature (Sv), SHR monomer in the endodermis (Se), SHR homodimer in the endodermis (S2e), SCR in the endodermis (C), 1:1 SHR-SCR complex in the endodermis (SC), and 2:1 SHR-SCR complex in the endodermis (S2C). The model consisted of six ordinary differential equations (ODEs) that measure how each of the variables changes over time.
We assumed that SHR is constantly produced at rate k1 as there is no information on upstream regulators of SHR. Since our pCF analysis showed that SHR only moves from the vasculature to the endodermis, possibly through an active transport mechanism (Gallagher et al., 2004; Sena et al., 2004; Gallagher and Benfey, 2009), we modeled the movement of SHR using an active transport term, where a1is the active transport rate. We defined a1 as the experimentally determined diffusion coefficient (D1) divided by the area of a vasculature cell (A1). We measured the area of vascular cells (n = 19) using ImageJ, and averaged them to determine A1 (Supplementary file 3). Although the diffusion coefficient returned by RICS is from a population of vascular and endodermal cells, we assumed that it is a good approximation of SHR movement between one vascular and one endodermal cell. We included a second active transport term for movement in the reverse direction, from the endodermis to the vasculature, where a2is the active transport rate. However, based on the pCF analysis, a2 = 0 since normally there is no bidirectional movement. Adding linear degradation gave us the equation for the change in SHR in the vasculature over time.
Given that SHR is not produced in the endodermis, there is no production term in the equation. Thus, the only SHR present in the endodermis is the SHR that moves from the vasculature. This leads to the equation for the change in SHR monomer in the endodermis over time.
The SHR homodimer forms from two SHR monomers. Our N&B analysis revealed that SHR homodimer does not form in a SCRi line. This suggested that the homodimer formation rate, k2, should depend on the concentration of SCR. To account for this, we modeled k2 as a logistic function of the concentration of SCR. Once the SCR concentration reaches a critical value C0, k2 will increase at rate k until it reaches a maximum value of L (Figure 5—figure supplement 2). We chose values for C0, k, and L based on the N&B data (see Parameter Estimation). Thus, the equation for the change in SHR homodimer over time is
Unlike the other variables, SCR production is not a linear term but rather a Hill equation. This structure was chosen because it has been shown that the SHR-SCR complex activates SCR expression (Cui et al, 2007). We assumed that both the 1:1 and 2:1 SHR-SCR complexes can activate SCR. In addition, SCR has been shown to autoregulate itself (Cui et al, 2007; Moreno-Risueno et al, 2015). Therefore, the change in SCR over time is expressed as:
Finally, our N&B analysis revealed that both the monomer and homodimer of SHR can bind SCR and form a complex. The final two equations in our model measure the change in these complexes over time.
In addition to the wild type model, we constructed a model that simulates SCRi (Model 2). Our pCF analysis revealed that SHR movement in the SCRi line is bidirectional, so now a2 is defined as the experimentally determined diffusion coefficient (D2) divided by the area of an endodermal cell (A2). We determined the average area of endodermal cells by averaging the area of representative cells (n = 19) as we did for vascular cells (Supplementary file 3). Since it had been shown that a 60% reduction of SCR is required to produce the mutant layer, we assumed that SCR concentrations are maintained below 60% and that the change in SCR over time is zero (Cui et al., 2007).
Sensitivity analysis
Request a detailed protocolThe sensitivity analysis was performed to determine the most influential parameters in our model. Notably, we reasoned that small changes in highly influential parameters could result in large changes in the model outcome. In addition, it has been shown that parameter estimation can become more computationally complex and produce more uncertainty in the parameter values as the number of estimated parameters increases (Smith, 2014). Therefore, we focused on estimating only the most influential parameters. We chose to use Sobol decomposition to measure how sensitive the model is to a particular parameter (Sobol, 2001). Sobol decomposition is a variance-based method, meaning that the sensitivity of the model to a parameter is quantified by calculating the variance in the model outcome. In addition, Sobol decomposition is a global sensitivity method, so we are exploring the entire parameter space in order to determine the most influential parameters (Smith, 2014).
Since the Sobol decomposition allows us to calculate numerous Sobol indices, we chose to use the total effects index to measure sensitivity in our model. Accordingly, the total effects index takes into account how sensitive the model is to a single parameter as well as combinations of more than one parameter. Thus, by using the total-effects index, we take into account any parameter interactions in our model.
To calculate the Sobol total effects index, we rewrote model 1 in the form where the Xi represent the 15 parameters. Y is the model outcome, which in this case must be a scalar. The ODE solution is a set of values over time, so the solution must be numerically integrated to obtain a single, constant value (Smith, 2014). Consequently, the total effects index for parameter i, STi, is defined as
where E( ) denotes the expected value, V( ) denotes the variance, and X~i is the vector of parameters without the ith parameter value. The numerator represents the expected variance in the model if all factors except Xi are fixed. If this term is very large, then that means that Xi contributes greatly to variance in the model: in other words, Xi is an influential parameter. We divided by the variance in the model outcome so that the value of the total index is normalized across different outcomes (Saltelli et al, 2010; Sobol, 2001).
We calculated the index for each parameter using Monte Carlo evaluations and built two matrices A and B, which are 1000 rows by 15 columns. Each row represented a random draw of all 15 parameters from the parameter space (Supplementary file 3). For each parameter i = 1,2,…,15, we constructed a third matrix C that was identical to A except that the ith column of A was replaced with the ith column of B. Model 1 was then evaluated for all 1000 rows. When we evaluated the ODE, we obtained 6 different outcomes: one for each variable we measured. To understand how each parameter affects all of the different variables we calculated the total effects index for each variable separately.
After evaluating the model, we numerically integrated the solution to obtain a constant value. The result is a 1000 by 6 matrix, where each row represents a model evaluation and each column represents one of the variables (SHR in the vasculature, SHR monomer, etc). We then used each column to evaluate (1), giving us the total effects index for parameter i with respect to each of the variables. Since the numerator of (1) cannot be computed exactly, we estimated it using one of the most accurate estimators,
where N is the number of Monte Carlo iterations (in our case, N = 1000) (Saltelli et al., 2010). Finally, we repeated the entire process 10 times to obtain technical replicates (Supplementary file 4). We used the Wilcoxon test with the Steel-Dwass test for multiple comparisons at significance level α = 0.1 to determine parameters that are significantly more influential.
Parameter estimation
Request a detailed protocolFor the initial model simulations, we derived most of the model parameters from the previous model of SHR and SCR in the CEI (Supplementary file 3). It has been shown that the mutant layer in the SCRi line occurs when the concentration of SCR drops to 60% of WT levels (Cui et al., 2007). Therefore, for the homodimer formation term k2, we chose C0 to be 60% of the steady-state value of SCR. By choosing this value, the homodimer will not form until the concentration of SCR is above 60%. Additionally, we chose k = 0.1 so that homodimer formation occurs rapidly after SCR passes the threshold (Supplementary file 3). We expected that our model would have high levels of SHR monomer and 1:1 complex relative to the levels of SHR homodimer and 2:1 complex as shown by our N&B data. However, using these parameters, we were unable to replicate these results.
We reasoned that estimating better values for some of the parameters could improve our model simulation. Thus, we used the results from our sensitivity analysis to identify the most influential parameters to estimate. We chose to estimate d2, k3, K2D, and d4 based on their high sensitivity indices (Figure 5 and Figure 5—figure supplement 1). Although k1 and d1 are influential parameters, we chose not to estimate them as varying their values only changes the overall level of SHR and not the oligomeric dynamics. We also chose to estimate L, which is the maximum value of k2. Although k2 is not as influential as some of the other parameters, its functional structure comes from the N&B analysis, so its maximum value should be derived from the N&B data.
For the parameter estimation, we started parameters at their default values (Supplementary file 3) and varied them one at a time until we were able to replicate the N&B data. To replicate the N&B data, we required that the SHR homodimer accounted for 7.5% of the total SHR proteins, while the 2:1 complex accounted for 15.2% of the total SHR-SCR complexes (Figure 3 and Figure 4). We found that increasing d2, decreasing K2D, and decreasing L achieved this relationship. In addition, decreasing K2D resulted in SCR increasing over time (Figure 5D). We did not change the values of k3 and d4. Using N&B to experimentally determine the oligomeric state and protein-protein stoichiometry of SHR and SCR allowed us to estimate better values for some parameters and, further, improve the reliability of the conclusions we derived from the model.
Data availability
-
Data from: Tracking transcription factor mobility and interaction in Arabidopsis roots with fluorescence correlation spectroscopyAvailable at Dryad Digital Repository under a CC0 Public Domain Dedication.
References
-
Measuring fast dynamics in solutions and cells with a laser scanning microscopeBiophysical Journal 89:1317–1327.https://doi.org/10.1529/biophysj.105.062836
-
Mapping the number of molecules and brightness in the laser scanning microscopeBiophysical Journal 94:2320–2332.https://doi.org/10.1529/biophysj.107.114645
-
Analysis of diffusion and binding in cells using the RICS approachMicroscopy Research and Technique 72:323–332.https://doi.org/10.1002/jemt.20655
-
Imaging barriers to diffusion by pair correlation functionsBiophysical Journal 97:665–673.https://doi.org/10.1016/j.bpj.2009.04.048
-
Lessons in fluctuation correlation spectroscopyAnnual Review of Physical Chemistry 62:645–668.https://doi.org/10.1146/annurev-physchem-032210-103424
-
Stoichiometry of molecular complexes at adhesions in living cellsProceedings of the National Academy of Sciences of the United States of America 106:2170–2175.https://doi.org/10.1073/pnas.0806036106
-
Uncovering the networks involved in stem cell maintenance and asymmetric cell division in the Arabidopsis rootCurrent Opinion in Plant Biology 29:38–43.https://doi.org/10.1016/j.pbi.2015.11.002
-
Mechanisms regulating SHORT-ROOT intercellular movementCurrent Biology 14:1847–1851.https://doi.org/10.1016/j.cub.2004.09.081
-
Transcription factor-mediated cell-to-cell signalling in plantsJournal of Experimental Botany 65:1737–1749.https://doi.org/10.1093/jxb/ert422
-
In vivo pair correlation analysis of EGFP intranuclear diffusion reveals DNA-dependent molecular flowProceedings of the National Academy of Sciences of the United States of America 107:16560–16565.https://doi.org/10.1073/pnas.1006731107
-
Measuring the flow of molecules in cellsBiophysical Reviews 3:119–129.https://doi.org/10.1007/s12551-011-0051-x
-
Subdomains for transport via plasmodesmata corresponding to the apical-basal axis are established during Arabidopsis embryogenesisProceedings of the National Academy of Sciences of the United States of America 102:11945–11950.https://doi.org/10.1073/pnas.0505622102
-
NAC transcription factors in plant abiotic stress responsesBiochimica Et Biophysica Acta (BBA) - Gene Regulatory Mechanisms 1819:97–103.https://doi.org/10.1016/j.bbagrm.2011.10.005
-
Distribution of molecular aggregation by analysis of fluctuation momentsProceedings of the National Academy of Sciences of the United States of America 87:5479–5483.https://doi.org/10.1073/pnas.87.14.5479
-
Fluorescence correlation spectroscopy with high count rate and low background: analysis of translational diffusionEuropean Biophysics Journal 22:169–175.https://doi.org/10.1007/BF00185777
-
Raster image correlation spectroscopy in live cellsNature Protocols 5:1761–1774.https://doi.org/10.1038/nprot.2010.122
-
Control of patterning, growth, and differentiation by floral organ identity genesJournal of Experimental Botany 66:1065–1073.https://doi.org/10.1093/jxb/eru514
-
Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity indexComputer Physics Communications 181:259–270.https://doi.org/10.1016/j.cpc.2009.09.018
-
Stem-cell niches: nursery rhymes across kingdomsNature Reviews. Molecular Cell Biology 8:345–354.https://doi.org/10.1038/nrm2164
-
BookUncertainty Quantification: Theory, Implementation, and ApplicationsPhiladelphia: Siam.
-
Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimatesMathematics and Computers in Simulation 55:271–280.https://doi.org/10.1016/S0378-4754(00)00270-6
-
Postembryonic control of root meristem growth and developmentCurrent Opinion in Plant Biology 17:7–12.https://doi.org/10.1016/j.pbi.2013.10.005
-
Temporally coordinated signals progressively pattern the anteroposterior and dorsoventral body axesSeminars in Cell & Developmental Biology 42:118–133.https://doi.org/10.1016/j.semcdb.2015.06.003
-
Callose biosynthesis regulates symplastic trafficking during root developmentDevelopmental Cell 21:1144–1155.https://doi.org/10.1016/j.devcel.2011.10.006
-
Signal integration during development: insights from the Drosophila eyeDevelopmental Dynamics 229:162–175.https://doi.org/10.1002/dvdy.10449
-
Estimation of diffusion coefficients of proteinsBiotechnology and Bioengineering 22:947–955.https://doi.org/10.1002/bit.260220504
Article and author information
Author details
Funding
National Science Foundation (MCB-1453130)
- Rosangela Sozzani
National Science Foundation (DGE-1252376)
- Natalie M Clark
- Adam P Fisher
National Institutes of Health (GM-106690-01)
- Cara M Winter
National Institutes of Health (P41-GM103540)
- Enrico Gratton
National Institutes of Health (P50-GM076516)
- Enrico Gratton
National Science Foundation (IOS-1021619)
- Philip N Benfey
National Institutes of Health (R01-GM043778)
- Philip N Benfey
Howard Hughes Medical Institute (GBMF3405)
- Philip N Benfey
Gordon and Betty Moore Foundation (GBMF3405)
- Philip N Benfey
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
NMC is supported by a NSF GRF (DGE-1252376). CMW is supported by a NRSA F32 (GM-106690-01). EG is supported by NIH (P41-GM103540 and P50-GM076516). APF is supported by a NSF GRF (DGE-1252376). This work was funded by a grant to PNB from the NSF Arabidopsis 2010 program (IOS-1021619), from the NIH (R01-GM043778), by the Howard Hughes Medical Institute and the Gordon and Betty Moore Foundation (through Grant GBMF3405), and by a NSF CAREER grant (MCB-1453130) to RS.
Copyright
© 2016, Clark et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,078
- views
-
- 1,120
- downloads
-
- 72
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Stem Cells and Regenerative Medicine
Stimulation of pancreatic beta cell regeneration could be a therapeutic lead to treat diabetes. Unlike humans, the zebrafish can efficiently regenerate beta cells, notably from ductal pancreatic progenitors. To gain insight into the molecular pathways involved in this process, we established the transcriptomic profile of the ductal cells after beta cell ablation in the adult zebrafish. These data highlighted the protein phosphatase calcineurin (CaN) as a new potential modulator of beta cell regeneration. We showed that CaN overexpression abolished the regenerative response, leading to glycemia dysregulation. On the opposite, CaN inhibition increased ductal cell proliferation and subsequent beta cell regeneration. Interestingly, the enhanced proliferation of the progenitors was paradoxically coupled with their exhaustion. This suggests that the proliferating progenitors are next entering in differentiation. CaN appears as a guardian which prevents an excessive progenitor proliferation to preserve the pool of progenitors. Altogether, our findings reveal CaN as a key player in the balance between proliferation and differentiation to enable a proper beta cell regeneration.
-
- Developmental Biology
Mutations in the well-known Myostatin (MSTN) produce a ‘double-muscle’ phenotype, which makes it commercially invaluable for improving livestock meat production and providing high-quality protein for humans. However, mutations at different loci of the MSTN often produce a variety of different phenotypes. In the current study, we increased the delivery ratio of Cas9 mRNA to sgRNA from the traditional 1:2 to 1:10, which improves the efficiency of the homozygous mutation of biallelic gene. Here, a MSTNDel73C mutation with FGF5 knockout sheep, in which the MSTN and FGF5 dual-gene biallelic homozygous mutations were produced via the deletion of 3-base pairs of AGC in the third exon of MSTN, resulting in cysteine-depleted at amino acid position 73, and the FGF5 double allele mutation led to inactivation of FGF5 gene. The MSTNDel73C mutation with FGF5 knockout sheep highlights a dominant ‘double-muscle’ phenotype, which can be stably inherited. Both F0 and F1 generation mutants highlight the excellent trait of high-yield meat with a smaller cross-sectional area and higher number of muscle fibers per unit area. Mechanistically, the MSTNDel73C mutation with FGF5 knockout mediated the activation of FOSL1 via the MEK-ERK-FOSL1 axis. The activated FOSL1 promotes skeletal muscle satellite cell proliferation and inhibits myogenic differentiation by inhibiting the expression of MyoD1, and resulting in smaller myotubes. In addition, activated ERK1/2 may inhibit the secondary fusion of myotubes by Ca2+-dependent CaMKII activation pathway, leading to myoblasts fusion to form smaller myotubes.