A new kymogram-based method reveals unexpected effects of marker protein expression and spatial anisotropy of cytoskeletal dynamics in plant cell cortex
© The Author(s) 2017
Received: 30 August 2016
Accepted: 22 March 2017
Published: 27 March 2017
Cytoskeleton can be observed in live plant cells in situ with high spatial and temporal resolution using a combination of specific fluorescent protein tag expression and advanced microscopy methods such as spinning disc confocal microscopy (SDCM) or variable angle epifluorescence microscopy (VAEM). Existing methods for quantifying cytoskeletal dynamics are often either based on laborious manual structure tracking, or depend on costly commercial software. Current automated methods also do not readily allow separate measurements of structure lifetime, lateral mobility, and spatial anisotropy of these parameters.
We developed a new freeware-based, operational system-independent semi-manual technique for analyzing VAEM or SDCM data, QuACK (Quantitative Analysis of Cytoskeletal Kymograms), and validated it on data from Arabidopsis thaliana fh1 formin mutants, previously shown by conventional methods to exhibit altered actin and microtubule dynamics compared to the wild type. Besides of confirming the published mutant phenotype, QuACK was used to characterize surprising differential effects of various fluorescent protein tags fused to the Lifeact actin probe on actin dynamics in A. thaliana cotyledon epidermis. In particular, Lifeact-YFP slowed down actin dynamics compared to Lifeact-GFP at marker expression levels causing no macroscopically noticeable phenotypic alterations, although the two fluorophores are nearly identical. We could also demonstrate the expected, but previously undocumented, anisotropy of cytoskeletal dynamics in elongated epidermal cells of A. thaliana petioles and hypocotyls.
Our new method for evaluating plant cytoskeletal dynamics has several advantages over existing techniques. It is intuitive, rapid compared to fully manual approaches, based on the free ImageJ software (including macros we provide here for download), and allows measurement of multiple parameters. Our approach was already used to document unexpected differences in actin mobility in transgenic A. thaliana expressing Lifeact fusion proteins with different fluorophores, highlighting the need for cautious interpretation of experimental results, as well as to reveal hitherto uncharacterized anisotropy of cytoskeletal mobility in elongated plant cells.
KeywordsActin Microtubules Lifeact Variable angle fluorescence microscopy Spinning disc confocal microscopy Kymogram Structure stability Lateral mobility Anisotropy FH1 (At3g25500)
Advanced microscopy methods such as spinning disc confocal microscopy (SDCM) or variable angle epifluorescence microscopy (VAEM) are nowadays becoming available to a wide community of researchers [1–4]. At the same time, new fluorescent protein markers are being developed, enabling “minimally invasive” detection of intracellular structures in vivo in practically any type of cells. Both main cytoskeletal systems—microtubules and microfilaments—are nowadays being routinely visualized in living cells of plant epidermis in situ using fluorescent protein marker expression in stable transgenic plants or transient transformants. Effects of biologically relevant conditions such as e.g. developmental changes, external influences or mutations on cytoskeletal organization and dynamics can be thus studied, leading to better understanding of the cellular mechanisms underlying plant ontogeny and responses towards environmental challenges.
Multiple markers for labeling either cytoskeletal system have been employed in plant studies. Widely used are fusion proteins based on fimbrin [5–7], albeit they can affect actin-dependent processes including organelle motility  and cause actin bundling when overexpressed . At present, a 17 amino acid peptide called Lifeact is probably the least artifact-prone microfilament tag available, with negligible effects on in vitro actin polymerization kinetics over a wide concentration range , although its overexpression may stabilize microfilaments in vivo . Expression levels and fluorophore properties can also influence actin structure and dynamics. For instance, red fluorescent proteins such as dsRed2 or mOrange either cause increased microfilament bundling or preferentially decorate existing bundles, while GFP and its derivatives fused to the same actin-binding moiety labeled finer filament networks. Overexpression of fimbrin derivatives tagged with mOrange in Arabidopsis seedlings caused growth inhibition, and even mild expression of a Lifeact-based marker slightly reduced root growth .
Several markers are available also for in vivo microtubule visualization [12, 13], and their use is plagued by similar problems. Frequently used are derivatives of the microtubule-associated protein MAP4 [14–16], which, however, induce organ twisting in Arabidopsis at high expression levels (e.g. ), indicating interference with cytoskeletal structure and dynamics, with subsequent effects on growth. Low-level expression of fluorescent protein-tagged tubulin subunits [18–20] can be used to label microtubules, since fluorescent protein-tagged tubulin, unlike tagged actin, can co-polymerize with native monomers. However, studies in a non-plant system (Dictyostelium) suggest that labeled tubulin may affect microtubule dynamics .
Only major changes to actin or microtubule organization and dynamics are obvious upon visual examination of microscopic images or video sequences. Additional insight can be obtained by inspection of kymograms—two-dimensional projections of the video sequence produced by plotting in the horizontal dimension pixels chosen from individual frames along a straight (linear) or curved transect, while the temporal dimension generates the vertical axis. Kymograms are sometimes used to qualitatively demonstrate plant cytoskeletal dynamics (e.g. [22–26]).
Quantitative image analysis is usually required to extract biologically relevant information from images and video sequences produced by advanced microscopy methods. Besides proprietary software (distributed often by the manufacturers of the microscopy equipment), freeware and shareware programs are being developed by the scientific community. In particular, ImageJ and its Fiji distribution  has recently become a standard part of the biologists´ toolkit, including utilities for generating kymograms.
For evaluation of plant cytoskeletal structure (i.e. static distribution of filaments), a number of particularly ImageJ-based protocols has been employed (e.g. [28–33]). Nevertheless, quantitative analysis of cytoskeletal dynamics remains challenging, especially in the case of rapidly changing and often weakly labeled actin structures, and techniques involving manual tracking or counting of individual filaments or filament bundles thus remain the method of choice (e.g. [23, 25, 26, 33, 34]). Such a manual approach is usually very laborious and may be prone to observer bias (which can be eliminated by a blinded study design involving an evaluator who does not know which sample is which). Development of partially or fully automated methods is thus desirable. To our knowledge, only one such technique—image-to-image correlation analysis—has found wider application in studies of plant cortical cytoskeleton dynamics [35, 36]. While this method produces informative and reproducible results, it is far from intuitive from the point of view of mathematically often somewhat naive biologists. It also can only measure a single composite parameter that reflects overall structure mobility, namely the rate at which the correlation coefficient between two subsequent images departs from the value of 1, characterizing a static system. Another drawback of the image-to-image correlation method is its dependence on the costly commercial MATLAB software package, as no version for free data analysis packages such as SciLab is available to date. Recently published ImageJ-compatible or standalone tools that enable semi-automated or automated extraction of quantitative information from kymograms [37, 38] might provide an interesting high-throughput alternative to manual single filament tracking, but their suitability for plant cells remains to be tested.
Here we present a new semi-manual, intuitive method for quantitative analysis of SDCM or VAEM recordings of fluorescence protein-labeled plant cytoskeleton. Our method, QuACK (Quantitative Analysis of Cytoskeletal Kymograms), is based on Fiji and can provide separate measurements of actin filaments or microtubules lifetime and lateral mobility. QuACK has been developed and validated using image data from wild type Arabidopsis thaliana plants and mutants defective in the actin-nucleating protein FH1 (At3g25500), a member of the formin (FH2 protein) family. Mutants lacking FH1 were previously shown by traditional methods to exhibit increased bundling and stability of microfilaments, as well as enhanced turnover of microtubule ends and generally faster microtubule dynamics [25, 26]. Besides of reproducing these phenotypes, QuACK revealed surprising differences in the dynamics of microtubules in plants expressing Lifeact with various fluorescent protein tags. We also successfully applied QuACK to characterize the spatial anisotropy of actin and microtubule dynamics of fully differentiated elongated epidermal cells of Arabidopsis hypocotyls and petioles. The results suggest, to our knowledge for the first time, reduced dynamics of actin bundles parallel to the longitudinal cell axis and decreased lateral mobility of microtubules aligned with the predominant direction of the cortical microtubule array.
To implement this concept, we developed a protocol consisting of several steps, some of them manual, others automated by means of ImageJ macros (Fig. 1). Prior to analysis, an initial selection of raw, microscope-generated video files, which must be all of the same lenght (i.e. timespan and number of frames) is performed based on visual inspection, discarding any obviously incomplete (shorter), out of focus or very noisy recordings, as well as those suffering by visible lateral drift. A sufficient number of biological replicates (individual plants/organs) and technical replicates (repeated recordings from different places of the same organ) should be used. Based on preliminary experiments in Arabidopsis seedlings, we recommend analyzing at least 10 recordings from at least five plants, with no more than 2 recordings originating from the same individual, per sample or treatment, and measuring at least 4 transects per recording, as a standard design. This setup was used to generate all the data shown in this paper unless stated otherwise. The number of transects may need to be increased to 6 or 8 to ensure statistical robustness for samples exhibiting only minor differences, or if comparing multiple samples. It is recommended to record at least 20% more video sequences than the above-suggested minimum to ensure that there is enough data left after discarding problematic files. If there is any doubt about the recording quality, it should be kept until the next (contrast enhancement) step, where visual selection is easier to make.
If working with isotropic or non-oriented samples (such as e.g. randomly positioned cotyledons), transects for generating kymograms can be located parallel to the horizontal or vertical edge of the video frame. For anisotropic samples such as e.g. root, hypocotyl or petiole epidermis, recordings must be aligned with a specific axis (e.g. organ axis, or the predominant direction of the cortical microtubule array), positioned either horizontally or vertically, prior to further processing. This can be done using built-in Fiji commands.
Movies are then converted to grayscale and contrast is enhanced by histogram stretching for each frame. This eliminates effects of photobleaching, i.e. gradual loss of fluorescence during recording. A randomly positioned grid is then overlaid over each movie. At this stage, final selection of recordings to be analyzed should be made, discarding videos with excessive amount of random noise (see below for criteria).
Transects of a pre-defined length are then positioned manually along the grid across a well-focused portion of the movie, as equally spaced as possible. (Optical sections from SDCM or VAEM usually have only a part of their area in focus. Thus, automated transect positioning, which would in theory eliminate a source of subjective error, would lead to unacceptable data loss.). Kymograms are then generated across the transects, converted to binary images and skeletonized. All resulting skeletons are manually measured to determine the length of the longest trace in the time (vertical) dimension and the width of the widest trace in the space (horizontal) dimension. Resulting numerical measurements for each condition (or set of experimental parameters) are then processed using a spreadsheet program such as Libre Office Calc or Microsoft Excel and statistically evaluated.
While it is tempting to use statistical approaches such as the t test or ANOVA to assess significance of between-group differences, this is not recommended, because long-living events (or large movements) are not completely recorded and thus the data may not fulfill the condition of normal distribution. Other methods, such as the Wilcoxon–Mann–Whitney test or the Kruskal–Wallis test, should thus be used. For the same reason, results should be presented as box plots with median and quartile/range rather than mean values ± SE or SD.
All the ImageJ macros required for performing the above-described procedure are provided, together with detailed step by step instructions, in Additional file 1.
Validation and optimization on actin data
Besides of missing relatively rare events (such as the stable actin bundles in the fh1 mutant), there is a theoretical possibility that frequent presence of stable structures might result in all or nearly all kymograms containing events that reach the maximum available value (i.e. the whole duration of the recording in case of lifetime values). While this was never a problem with actin dynamics, we had to address this issue when applying QuACK to microtubules (see below).
To further examine the effects of noise on QuACK analysis results, we compared actin structure lifetime obtained from a particular wt and fh1 dataset using the standard protocol with results derived from the same videos previously subjected to controlled introduction of random noise. On these noisy data, the difference between mutant and wt actin dynamics appears to be diminishing (Fig. 3b), obviously due to fragmentation of long kymogram traces (Fig. 3c), which were mainly responsible for the increased lifetime in the mutant. Thus, to avoid artifacts due to trace fragmentation, video recordings with a similar level of noise should be compared, and very noisy files should be discarded prior to kymogram construction. The level of noise can be best assessed visually after the contrast adjustment and grid insertion step (Fig. 3c). If comparison between videos of varying noise levels cannot be avoided, interpretation should reflect the tendency towards loss of long traces in noisy data.
Validation and optimization on microtubule data
Microtubules in plant cells are generally easier to visualize with an excellent signal to noise ratio, and their dynamics, slower than in case of microfilaments, is also relatively easy to quantify. Nevertheless, the general limitations of conventional methods, as outlined in the Introduction, apply to them as well. We thus examined the applicability of the QuACK approach to microtubules using existing VAEM recordings of cortical microtubule dynamics in cotyledon pavement cells of wt and fh1 mutant Arabidopsis seedlings expressing GFP-MAP4, produced for our recently published study . Individual microtubule end tracking has previously revealed an increased microtubule end turnover and decreased lateral microtubule end movement in the mutant compared to the wt, while an overall increase of microtubule mobility in the mutant was determined by image to image correlation measurements .
We have next employed QuACK to address two model questions, namely (i) how do various fluorophores affect the dynamics of microfilaments tagged by constructs utilizing the same actin-binding moiety, and (ii) how anisotropic are the lifetime and lateral mobility parameters of cortical microfilaments and microtubules in differentiated elongated plant cells.
Case study I: effects of fluorescent protein markers on actin dynamics
We used the same data for quantification of actin dynamics using QuACK, with inclusion of an additional Lifeact-YFP-expressing line which exhibited very high fluorescence intensity and some reduction in growth and morphological abnormalities (the Lifeact-YFP overexpression, or Lifeact-YFP ox, line). The qualitative differences previously detected by the image to image correlation technique were confirmed, and an even more obvious decrease in actin dynamics was apparent in the Lifeact-YFP ox line (Fig. 8b). In addition, QuACK analysis revealed that the observed differences are due to significant differences both in actin structures lifetime and their lateral mobility.
Case study II: anisotropy of cytoskeletal dynamics in elongated cells
Spatial anisotropy of cytoskeletal organization in elongating or elongated cells is a phenomenon is well-documented especially for microtubules [41, 42], but reported also for actin ; however, much less is known about the anisotropy of cytoskeletal dynamics in such cells.
Here we present a new method, QuACK (Quantitative Analysis of Cytoskeletal Kymograms) for evaluating the in vivo dynamics of plant cortical cytoskeleton from SDCM or VAEM video recordings obtained from cells tagged by expression of actin- or microtubule-specific fluorescent protein markers. Our approach is based on quantitative analysis of kymograms, projections of the video recordings to two dimensions with a transect across the image field providing the x dimension and time contributing the y dimension. Kymograms are alternatively termed “kymographs” in the literature (e.g. [37, 44]), albeit this word probably ought to be reserved to hardware or software devices producing these graphical representations of motion, at least per analogy with similar terms such as “telegram” versus “telegraph” . While kymogram-based techniques have been introduced into quantitative cytoskeletal dynamics studies previously, they were usually used to follow the movement of individual filaments, bundles or filament ends along linear or curved trajectories with subsequent analyses of population behavior (e.g. [37, 44, 46]) rather than for capturing parameters reflecting the overall dynamics of the actin or microtubule networks. Moreover, these methods are often based on sophisticated theoretical background that may make them less than user-friendly from the point of view of experimental biologists. This is also the case of the (to our knowledge only) plant study employing kymogram analysis to study global microtubule dynamics .
QuACK presents several advantages over commonly used techniques for estimating cytoskeletal dynamics, such as manual tracing of individual filaments or their ends (e.g. [3, 23, 25, 26, 33, 34]) or automated image to image correlation [35, 36]. It is faster than manual tracing, relatively simple and intuitive. Unlike the correlation method, it independently estimates the contribution of dynamic filament end turnover (reflected in the maximum event duration parameter) and lateral mobility (expressed as the maximum lateral displacement parameter) to overall dynamics of the structures under study. QuACK analysis can be performed using exclusively free software on any commonly used operation system, since both ImageJ/Fiji and freely available spreadsheet programs exist in versions compatible with Microsoft Windows, Linux and Mac-OS. While the method involves some manual steps that introduce a theoretical possibility of observer bias, we found this not to present a practical problem during testing (as even experienced observers could not predict the appearance of the skeletonized kymogram from the source video). Results of our method appear to be somewhat sensitive to the amount of noise in input video data. However, we provide suggestions on post hoc qualitative evaluation of input data based on the frequency of fragmented traces and meshwork patterns.
To validate our approach, we re-examined raw video data from our published study examining the effects of functional impairment of the Arabidopsis FH1 formin (encoded by the At3g25500 locus) on cortical actin and microtubule dynamics . Using static image analysis and manual tracking of filaments labeled with several different marker proteins, we have previously shown that fh1 mutants exhibit increased microfilament bundling and thus increased actin structure lifetime both in rhizodermis and in cotyledon pavement cells [25, 26], possibly as a consequence of weakening the trans-plasmalemma attachment of cortical microfilaments to the cell wall  that may allow for easier lateral movement and thus bundling of individual microfilaments . QuACK analyses convincingly confirmed increased lifetime of actin structures in fh1 mutants. However, lateral movement of recorded actin structures did not significantly differ between the mutant and the wt, most likely because increased lateral mobility of individual microfilaments was compensated for by decreased mobility of the frequently occurring thick bundles in the mutant. Indeed, such extremely stable bundles are frequently found in mutant plants [25, 26]. At present, QuACK does not distinguish between thick and thin bundles or even single filaments. While development in such a direction is conceivable, utilizing, e.g., thresholding functions of ImageJ, it is beyond the scope of the present work.
Unlike microfilaments, microtubules are generally more dynamic in fh1 mutants than in wt plants, as determined by individual filament end tracking, which showed increased dynamic instability of microtubule ends [25, 26], and by the image to image correlation method that documented an overall increase in microtubule rearrangements . Using QuACK, we have demonstrated both a decrease in microtubule lifetime and an increase in overall lateral mobility in the fh1 mutants. The later is somewhat surprising, since decreased frequency of microtubule end translocations (“end waving”) was found in the mutant. Since QuACK cannot distinguish between lateral movements of whole microtubules and their ends, our observations suggest that the later may represent only a minor fraction of overall microtubule mobility. This is also consistent with the change in microtubule end behavior being one of the very few aspects of the fh1 mutant phenotype that could not be replicated by pharmacological inhibition of formin activity , and with the noticeable, though not statistically significant, trend towards increased lateral microtubule mobility detected by individual microtubule tracking in the present study.
After validating and optimizing the QuACK method on the formin mutant data, we employed this approach to characterize microfilament dynamics in A. thaliana lines expressing the actin marker Lifeact fused to various fluorophores, expressed under the control of the UBQ10 promoter, previously proposed as a promoter of choice for live cytoskeletal labeling due to its ability to achieve relatively low expression levels, minimizing thus overexpression artifacts . While transgenes driven by the 35S promoter often exhibit silencing when crossed into the genetic backgrounds associated with commonly used T-DNA lines, including the SALFK collection , no such silencing occurred in our lines after crosses into several SALK-derived lines which previously exhibited silencing of 35S promoter constructs past the early seedling stage . Thus, actin can be now followed in all developmental stages of Arabidopsis mutants.
We found surprising differences in both microfilament lifetime and lateral mobility depending on the fluorophore used, even at comparable marker expression levels. Lifeact-YFP-expressing plants had less dynamic actin structures compared to those carrying Lifeact-GFP, while microfilaments labeled by Lifeact-mRFP were more dynamic than Lifeact-GFP-tagged ones. Fluorophore-dependent differences in actin organization have been reported for fimbrin- based markers, with mOrange labeling aggregated rather than filamentous structures compared to GFP-, YFP- or CFP-based tags . While mOrange is a derivative of dsRed, a protein that, unlike GFP, must oligomerize in order to become fluorescent, it is a descendant of a non-dimerizing mutant ; the reported effects are thus unlikely to be due solely to fluorophore dimerization or oligomerization. However, differences in fluorophore sequence (and thus structure) might affect the tagged proteins´ interactions with microfilaments or with proteins modulating actin organization and dynamics, leading to selective tagging of distinct microfilament structures or to changes in tagged structure behavior. In case of Lifeact, a recent study demonstrated that neither the position nor the identity of the fluorophore (mEGFP or mCherry) affected its affinity to actin in vitro. Interestingly, at low concentrations, Lifeact-mCherry even promoted actin nucleation and filament elongation, while at high doses it exhibited a clear inhibitory effect . The latter is also consistent with inhibitory effects of high levels of Lifeact expression on growth observed by others  and in our present study. However, the differences between Lifeact-GFP and Lifeact-YFP behavior are rather surprising, given the 98% sequence identity between these proteins, which differ only in several amino acids located inside the folded tertiary structure . Major changes in protein surface properties are thus unlikely, albeit subtle changes in conformation may affect protein–protein interactions. Nevertheless, YFP fusions of the prokaryotic actin homolog MreB were reported to exhibit distinctly filamentous organization not reproduced by homologous GFP-tagged constructs in another bacterial species , hinting at possible hitherto overlooked effects of subtle fluorophore differences. In any case, our data show that there can never be too much caution when interpreting results obtained using fluorescent protein tags expression.
We have next employed QuACK for initial characterization of spatial anisotropy of cytoskeletal dynamics in two markedly polar cell types, namely in differentiated, elongated epidermal cells of true leaf petioles (where actin was observed) and light-grown hypocotyls (where we followed microtubule behavior).
While much attention has been devoted to the signaling and regulatory aspects of actin rearrangements , very few studies focus on the anisotropy of microfilament organization and dynamics, usually in the context of pathological situations [55, 56], cell division  or tip growth  rather than during or following diffuse anisotropic cell expansion. Our analyses in mature elongated epidermal cells document decreased turnover and lateral mobility of actin filaments or bundles parallel to the longitudinal cell axis compared to those crossing this axis. This supports the established model (derived in part from work in Chara and in tip-growing cells) where longitudinal actin bundles provide “tracks” for cytoplasmic streaming, which can be expected to be rather stable compared to the rest of the highly dynamic actin meshwork [59–61]. Actin organization and dynamics is also known to be intimately interlinked with that of microtubules (e.g. [30, 62]); it thus remains to be clarified what is the role of the microtubule cytoskeleton in controlling the anisotropic behavior of microfilaments.
Compared to microfilaments, the role of cortical microtubules in polarized diffuse cell growth is well characterized [41, 42]. In longitudinally expanding cells, microtubules undergo constant dynamic rearrangements  resulting in gradual rotation of highly organized cortical arrays, whose orientation correlates with that of nascent cellulose microfibril deposition . Cortical microtubules with reduced lateral mobility tend to form arrays perpendicular to the direction of maximal growth in expanding shoot apical meristem cells . Our observations in fully expanded Arabidopsis hypocotyl cells suggest a reciprocal conclusion—namely that those microtubules which are collinear with the array tend to be less mobile. The cortical array can therefore be understood as an “attractor” stabilizing the position of microtubules once they become part of it.
Indirect evidence suggests that the stability of cortical microtubules may reflect their orientation, since in growing maize coleoptiles transversal microtubules participating in the cortical array are enriched in de-tyrosinylated tubulin that inhibits microtubule turnover . While we did not observe increased lifetime of microtubules aligned to the cortical array, predicted by this hypothesis, we did find that these microtubules are less laterally mobile. It is, of course, possible that observations from growing maize coleoptiles cannot be generalized to fully expanded epidermal cells of Arabidopsis hypocotyls (and vice versa). The differences in lateral microtubule mobility observed in our experiments might still reflect different post-translational modifications of aligned versus non-aligned microtubules, resulting in altered affinity towards microtubule-associated proteins that control lateral mobility, possibly via linking microtubules to the cell wall matrix through transmembrane proteins such as e.g. certain formins [48, 66] or components of the cell wall biosynthesis machinery .
In any case, we hope that QuACK can be used as a versatile tool in future studies aimed at closer characterization of cytoskeletal dynamics anisotropy in a variety of cell types, developmental situations and environmental contexts.
We provide a simple, freeware-based, operational system-independent method for estimating the lifetime and lateral mobility of fluorescent protein-tagged cytoskeletal structures in plant cells. The technique, which is suitable for analyzing the dynamics of both actin and microtubules from SDCM or VAEM recordings, was validated by confirming the previously documented increase in microfilament assembly stability and microtubule turnover in Arabidopsis fh1 mutants and then applied to characterize the unexpected effects of various fluorescent marker fusions on the dynamics of the actin marker Lifeact, as well as to obtain initial insights into the anisotropy of the turnover and lateral mobility of cortical actin and microtubules in elongated plant cells.
A. thaliana ecotype Columbia lines carrying the Lifeact-mRFP  and GFP- MAP4  markers in wt or fh1-1 mutant background have been described previously [25, 26]. Lifeact-GFP and Lifeact-YFP marker lines were prepared using following constructs: Lifeact in pDONR207  was transferred to the pUBC-GFP-dest and pUBC-YFP-dest vectors  by the LR Clonase II (Invitrogen). Both constructs were then transformed into A. thaliana Col-0 plants using the floral dip method and transformants selected using BASTA. Seeds produced by T2 generation plants have been used to grow seedlings for the experiments reported here.
For propagation and transformation, plants were grown in soil or peat pellets (Jiffy). Seedlings for imaging were grown in vitro at 22 °C with a 16 h light/8 h dark cycle on vertical Murashige and Skoog (MS) plates from surface-sterilized seeds stratified by several days of post-imbibition storage at 4 °C to synchronize germination. 5 day old seedlings were observed unless stated otherwise.
Image data acquisition
Image series (videos) were recorded for 2 min unless stated otherwise from seedlings mounted in water on chambered slides as described previously [3, 25, 26]. SDCM recordings were performed using an inverted spinning disc confocal microscope (Yokogawa CSU-X1 on a Nikon Ti-E platform, laser box Agilent MLC400, camera Andor Ixon) with plan apochromat × 100 oil (NA = 1.45) lens, laser lines set at 488 and 561 nm, and image interval 1 s. VAEM recordings were generated using the Leica AF6000 LX fluorescence platform with integrated TIRF module and a Leica DFC350FXR2 digital camera, with plan apochromat × 100 oil (NA = 1.46) lens, 400 nm peak excitation, and 210 ms exposure time, allowing for image interval 0.5 s.
Image data analysis
The Fiji distribution of ImageJ  was used for all image analysis steps except frame-to-frame autocorrelation measurements. Microfilament pause duration and fraction of pausing microtubules was quantified by visual end tracking (aided by manual labeling of tracked structures in ImageJ) as described previously . Lateral motility of microfilament structures was assessed after prior uniform contrast enhancement of video recordings by scoring the presence or absence of movement for randomly selected points located within a structure (using the point tool in ImageJ) over the span of 10 frames. Lateral motility of microtubules was determined in an analogous manner, following the fate of selected points over the span of 2 min. In all manual analyses, at least 120 filament ends or points in at least 12 cells from at least 7 plants per genotype were evaluated.
A summary of the QuACK method is described in Results and a detailed protocol provided in Additional file 1. Before QuACK measurements in anisotropic samples, a straight marker line was manually placed along a visually detected landmark, i.e. the organ axis or several parallel thick bundles of the cortical microtubule array. In the later case, marker line direction was chosen to minimize intersections with prominent cortical microtubules or bundles. Subsequently, the movie frame was rotated to align the marker line with the horizontal axis (for detailed instructions see Additional file 1).
To determine overall microtubule or actin dynamics by the image-to-image correlation method, correlation coefficients of pairs of video frames separated by varying intervals were calculated according to , using at least 10 movies from at least five plants per sample.
Online utilities [69–71] or a freely available MS Excel template  have been used for statistical evaluation using the Kruskal–Wallis, Wilcoxon–Mann–Whitney and Chi square tests. The various tools were used alternately after verification that they provide identical results on a subset of data. Box plots were generated using the BoxPlotR tool .
Formin homolog 1
Green fluorescent protein
Microtubule-associated protein 4
Monomeric red fluorescent protein
Quantitative analysis of cytoskeletal kymograms
Spinning disc confocal microscopy
Variable angle epifluorescence microscopy
Yellow fluorescent protein
FC designed the method, wrote the macros, performed some of the QuACK analyses and drafted the manuscript. DO constructed the transgenic plant lines, recorded the SDCM video sequences, conducted the image correlation analyses, performed some of the QuACK measurements and wrote parts of the text. Both authors have read and approved the final manuscript.
We thank Amparo Rosero for recording the VAEM data from the previous published studies which were re-analyzed in this paper, Ben Staiger for sharing MATLAB scripts for image to image correlation measurements, Ondrej Horváth for VAEM microscopy support, Jaroslav Flegr for advice in statistics, Jana Albrechtová for helpful discussion, and Kateřina Schwarzerová for critical reading of the manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Consent for publication
Re-use of our own data previously published as Fig. 8D of  in Fig. 4c of the present paper is permitted by the copyright policy of the original publisher (Oxford University Press), who has been duly notified and acknowledged the notification.
Ethics approval and consent to participate
The experiments with genetically modified bacteria and A. thaliana were performed under closed laboratory conditions in a facility approved according to the valid legislation of the Czech Republic, which does not require ethics approval for this kind of work.
This work has been supported by the Czech Science Foundation (GACR) 15-02610S and Ministry of Education, Youth and Sports (MSMT) LO14017 Grants.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gräf R, Rietdorf J, Zimmermann T. Live cell spinning disk microscopy. Adv Biochem Eng Biotechnol. 2005;95:57–75.PubMedGoogle Scholar
- Konopka CA, Bednarek SY. Variable-angle epifluorescence microscopy: a new way to look at protein dynamics in the plant cell cortex. Plant J. 2008;53:186–96.View ArticlePubMedGoogle Scholar
- Rosero A, Žárský V, Cvrčková F. Visualizing and quantifying the in vivo structure and dynamics of the Arabidopsis cortical cytoskeleton using CLSM and VAEM. Methods Mol Biol. 2014;1080:87–97.View ArticlePubMedGoogle Scholar
- Higaki T. Real-time imaging of plant cell surface dynamics with variable-angle epifluorescence microscopy. JoVE. 2015;106:e53437.Google Scholar
- Wang YS, Motes CM, Mohamalawari DR, Blancaflor EB. Green fluorescent protein fusions to Arabidopsis fimbrin 1 for spatio-temporal imaging of F-actin dynamics in roots. Cell Motil Cytoskeleton. 2004;59:79–93.View ArticlePubMedGoogle Scholar
- Sheahan MB, Staiger CJ, Rose RJ, McCurdy DW. A green fluorescent protein fusion to actin-binding domain 2 of Arabidopsis fimbrin highlights new features of a dynamic actin cytoskeleton in live plant cells. Plant Physiol. 2004;136:3968–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Voigt B, Timmers AC, Šamaj J, Müller J, Baluška F, Menzel D. GFP-FABD2 fusion construct allows in vivo visualization of the dynamic actin cytoskeleton in all cells of Arabidopsis seedlings. Eur J Cell Biol. 2005;84:595–608.View ArticlePubMedGoogle Scholar
- Holweg CL. Living markers for actin block myosin-dependent motility of plant organelles and auxin. Cell Motil Cytoskeleton. 2007;64:69–81.View ArticlePubMedGoogle Scholar
- Dyachok J, Sparks JA, Liao F, Wang YS, Blancaflor EB. Fluorescent protein-based reporters of the actin cytoskeleton in living plant cells: fluorophore variant, actin binding domain, and promoter considerations. Cytoskeleton. 2014;71:311–27.View ArticlePubMedGoogle Scholar
- Riedl J, Crevenna AH, Kessenbrock K, Yu JH, Neukirchen D, Bista M, Bradke F, Jenne D, Holak TA, Werb Z, Sixt M, Wedlich-Soldner R. Lifeact: a versatile marker to visualize F-actin. Nat Methods. 2008;5:605–7.View ArticlePubMedPubMed CentralGoogle Scholar
- van der Honing HS, van Bezouwen LS, Emons AM, Ketelaar T. High expression of Lifeact in Arabidopsis thaliana reduces dynamic reorganization of actin filaments but does not affect plant development. Cytoskeleton. 2011;68:578–87.View ArticlePubMedGoogle Scholar
- Buschmann H, Sambade A, Pesquet E, Calder G, Lloyd CW. Microtubule dynamics in plant cells. Methods Cell Biol. 2010;97:373–400.View ArticlePubMedGoogle Scholar
- Celler K, Fujita M, Kawamura E, Ambrose C, Herburger K, Holzinger A, Wasteneys GO. Microtubules in plant cells: strategies and methods for immunofluorescence, transmission electron microscopy, and live cell imaging. Methods Mol Biol. 2016;1365:155–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Marc J, Granger CL, Brincat J, Fisher DD, Kao Th, McCubbin AG, Cyr RJ. A GFP-MAP4 reporter gene for visualizing cortical microtubule rearrangements in living epidermal cells. Plant Cell. 1998;10:1927–40.PubMedPubMed CentralGoogle Scholar
- Granger CL, Cyr RJ. Microtubule reorganization in tobacco BY-2 cells stably expressing GFP-MBD. Planta. 2000;210:502–9.View ArticlePubMedGoogle Scholar
- Camilleri C, Azimzadeh J, Pastuglia M, Bellini C, Grandjean O, Bouchez D. The Arabidopsis TONNEAU2 gene encodes a putative novel protein phosphatase 2A regulatory subunit essential for the control of the cortical cytoskeleton. Plant Cell. 2002;14:833–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Hashimoto T. Molecular genetic analysis of left-right handedness in plants. Philos Trans R Soc Lond B Biol Sci. 2002;357:799–808.View ArticlePubMedPubMed CentralGoogle Scholar
- Ueda K, Matsuyama T, Hashimoto T. Visualization of microtubules in living cells of transgenic Arabidopsis thaliana. Protoplasma. 1999;206:201–6.View ArticleGoogle Scholar
- Kumagai F, Yoneda A, Tomida T, Sano T, Nagata T, Hasezawa S. Fate of nascent microtubules organized at the M/G1 interface, as visualized by synchronized tobacco BY-2 cells stably expressing GFP-tubulin: time-sequence observations of the reorganization of cortical microtubules in living plant cells. Plant Cell Physiol. 2001;42:723–32.View ArticlePubMedGoogle Scholar
- Nakamura M, Naoi K, Shoji T, Hashimoto T. Low concentrations of propyzamide and oryzalin alter microtubule dynamics in Arabidopsis epidermal cells. Plant Cell Physiol. 2004;45:1330–4.View ArticlePubMedGoogle Scholar
- Kimble M, Kuzmiak C, McGovern KN, de Hostos EL. Microtubule organization and the effects of GFP-tubulin expression in Dictyostelium discoideum. Cell Motil Cytoskeleton. 2000;47:48–62.View ArticlePubMedGoogle Scholar
- Shaw SL, Kamyar R, Ehrhardt DW. Sustained microtubule treadmilling in Arabidopsis cortical arrays. Science. 2003;300:1715–8.View ArticlePubMedGoogle Scholar
- Ambrose C, Wasteneys GO. Cell edges accumulate gamma tubulin complex components and nucleate microtubules following cytokinesis in Arabidopsis thaliana. PLoS ONE. 2011;6:e27423.View ArticlePubMedPubMed CentralGoogle Scholar
- Sambade A, Pratap A, Buschmann H, Morris RJ, Lloyd C. The influence of light on microtubule dynamics and alignment in the Arabidopsis hypocotyl. Plant Cell. 2012;24:192–201.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosero A, Žárský V, Cvrčková F. AtFH1 formin mutation affects actin filament and microtubule dynamics in Arabidopsis thaliana. J Exp Bot. 2013;64:585–97.View ArticlePubMedGoogle Scholar
- Rosero A, Oulehlová D, Stillerová L, Schiebertová P, Grunt M, Žárský V, Cvrčková F. Arabidopsis FH1 formin affects cotyledon pavement cell shape by modulating cytoskeleton dynamics. Plant Cell Physiol. 2016;57:488–504.View ArticlePubMedGoogle Scholar
- Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Schmid B, Tinevez JY, White DJ, Hartenstein V, Eliceiri K, Tomancak P, Cardona A. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9:676–82.View ArticlePubMedGoogle Scholar
- Higaki T, Kutsuna N, Sano T, Hasezawa S. Quantitative analysis of changes in actin microfilament contribution to cell plate development in plant cytokinesis. BMC Plant Biol. 2008;8:80.View ArticlePubMedPubMed CentralGoogle Scholar
- Higaki T, Kutsuna N, Sano T, Kondo N, Hasezawa S. Quantification and cluster analysis of actin cytoskeletal structures in plant cells: role of actin bundling in stomatal movement during diurnal cycles in Arabidopsis guard cells. Plant J. 2010;61:156–65.View ArticlePubMedGoogle Scholar
- Sampathkumar A, Lindeboom JJ, Debolt S, Gutierrez R, Ehrhardt DW, Ketelaar T, Persson S. Live cell imaging reveals structural associations between the actin and microtubule cytoskeleton in Arabidopsis. Plant Cell. 2011;23:2302–13.View ArticlePubMedPubMed CentralGoogle Scholar
- van der Honing HS, Kieft H, Emons AM, Ketelaar T. Arabidopsis VILLIN2 and VILLIN3 are required for the generation of thick actin filament bundles and for directional organ growth. Plant Physiol. 2012;158:1426–38.View ArticlePubMedGoogle Scholar
- Akita K, Higaki T, Kutsuna N, Hasezawa S. Quantitative analysis of microtubule orientation in interdigitated leaf pavement cells. Plant Signal Behav. 2015;10:e1024396.View ArticlePubMedPubMed CentralGoogle Scholar
- Chi Z, Ambrose C. Microtubule encounter-based catastrophe in Arabidopsis cortical microtubule arrays. BMC Plant Biol. 2016;16:18.View ArticlePubMedPubMed CentralGoogle Scholar
- Dixit R, Cyr R. Encounters between dynamic cortical microtubules promote ordering of the cortical array through angle-dependent modifications of microtubule behavior. Plant Cell. 2004;16:3274–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Vidali L, Burkart GM, Augustine RC, Kerdavid E, Tüzel E, Bezanilla M. Myosin XI is essential for tip growth in Physcomitrella patens. Plant Cell. 2010;22:1868–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Li J, Staiger BH, Henty-Ridilla JL, Abu-Abied M, Sadot E, Blanchoin L, Staiger CJ. The availability of filament ends modulates actin stochastic dynamics in live plant cells. Mol Biol Cell. 2014;25:1263–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Mangeol P, Prevo B, Peterman EJ. KymographClear and KymographDirect: two tools for the automated quantitative analysis of molecular and cellular dynamics using kymographs. Mol Biol Cell. 2016;27:1948–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Chaphalkar AR, Jain K, Gangan MS, Athale CA. Automated multi-peak tracking kymography (AMTraK): a tool to quantify sub-cellular dynamics with sub-pixel accuracy. PLoS ONE. 2016;11:e0167620.View ArticlePubMedPubMed CentralGoogle Scholar
- Daxinger L, Hunter B, Sheikh M, Jauvion V, Gasciolli V, Vaucheret H, Matzke M, Furner I. Unexpected silencing effects from T-DNA tags in Arabidopsis. Trends Plant Sci. 2008;13:4–6.View ArticlePubMedGoogle Scholar
- Fendrych M, Synek L, Pečenková T, Janková Drdová E, Sekereš J, de Rycke R, Nowack MK, Žárský V. Visualization of the exocyst complex dynamics at the plasma membrane of Arabidopsis thaliana. Mol Biol Cell. 2013;24:510–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Crowell EF, Gonneau M, Vernhettes S, Höfte H. Regulation of anisotropic cell expansion in higher plants. C R Biol. 2010;333:320–4.View ArticlePubMedGoogle Scholar
- Landrein B, Hamant O. How mechanical stress controls microtubule behavior and morphogenesis in plants: history, experiments and revisited theories. Plant J. 2013;75:324–38.View ArticlePubMedGoogle Scholar
- Jacques E, Lewandowski M, Buytaert J, Fierens Y, Verbelen JP, Vissenberg K. MicroFilament Analyzer identifies actin network organizations in epidermal cells of Arabidopsis thaliana roots. Plant Signal Behav. 2013;8:e24821.View ArticlePubMedPubMed CentralGoogle Scholar
- Smal I, Grigoriev I, Akhmanova A, Niessen WJ, Meijering E. Microtubule dynamics analysis using kymographs and variable-rate particle filters. IEEE Trans Image Process. 2010;19:1861–76.View ArticlePubMedGoogle Scholar
- Oxford Dictionaries: Telegram. http://www.oxforddictionaries.com/definition/english/telegram. Accessed 22 Aug 2016.
- Hamilton RS, Parton RM, Oliveira RA, Vendra G, Ball G, Nasmyth K, Davis I. ParticleStats: open source software for the analysis of particle motility and cytoskeletal polarity. Nucleic Acids Res. 2010;38:W641–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Dhonukshe P, Vischer N, Gadella TW Jr. Contribution of microtubule growth polarity and flux to spindle assembly and functioning in plant cells. J Cell Sci. 2006;119:3193–205.View ArticlePubMedGoogle Scholar
- Martinière A, Gayral P, Hawes C, Runions J. Building bridges: formin1 of Arabidopsis forms a connection between the cell wall and the actin cytoskeleton. Plant J. 2011;66:354–65.View ArticlePubMedGoogle Scholar
- Cvrčková F, Oulehlová D, Žárský V. On growth and formins. Plant Signal Behav. 2016;11:e1155017.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaner NC, Steinbach PA, Tsien RY. A guide to choosing fluorescent proteins. Nat Methods. 2005;2:905–9.View ArticlePubMedGoogle Scholar
- Courtemanche N, Pollard TD, Chen Q. Avoiding artefacts when counting polymerized actin in live cells with LifeAct fused to fluorescent proteins. Nat Cell Biol. 2016;18:676–83.View ArticlePubMedGoogle Scholar
- Wachter RM, Elsliger MA, Kallio K, Hanson GT, Remington SJ. Structural basis of spectral shifts in the yellow-emission variants of green fluorescent protein. Structure. 1998;6:1267–77.View ArticlePubMedGoogle Scholar
- Margolin W. The price of tags in protein localization studies. J Bacteriol. 2012;194:6369–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Li J, Blanchoin L, Staiger CJ. Signaling to actin stochastic dynamics. Annu Rev Plant Biol. 2015;66:415–40.View ArticlePubMedGoogle Scholar
- Goodbody KC, Lloyd CW. Actin filaments line up across Tradescantia epidermal cells, anticipating wound-induced divison planes. Protoplasma. 1990;157:92–101.View ArticleGoogle Scholar
- Tilsner J, Linnik O, Wright KM, Bell K, Roberts AG, Lacomme C, Santa Cruz S, Oparka KJ. The TGB1 movement protein of Potato virus X reorganizes actin and endomembranes into the X-body, a viral replication factory. Plant Physiol. 2012;158:1359–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Panteris E. Cortical actin filaments at the division site of mitotic plant cells: a reconsideration of the ‘actin-depleted zone’. New Phytol. 2008;179:334–41.View ArticlePubMedGoogle Scholar
- Qu X, Jiang Y, Chang M, Liu X, Zhang R, Huang S. Organization and regulation of the actin cytoskeleton in the pollen tube. Front Plant Sci. 2015;5:786.View ArticlePubMedPubMed CentralGoogle Scholar
- Shimmen T, Hamatani M, Saito S, Yokota E, Mimura T, Fusetani N, Karaki H. Roles of actin filaments in cytoplasmic streaming and organization of transvacuolar strands in root hair cells of Hydrocharis. Protoplasma. 1995;185:188–93.View ArticleGoogle Scholar
- Verchot-Lubicz J, Goldstein RE. Cytoplasmic streaming enables the distribution of molecules and vesicles in large plant cells. Protoplasma. 2010;240:99–107.View ArticlePubMedGoogle Scholar
- Woodhouse FG, Goldstein RE. Cytoplasmic streaming in plant cells emerges naturally by microfilament self-organization. Proc Natl Acad Sci USA. 2013;110:14132–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Havelková L, Nanda G, Martinek J, Bellinvia E, Sikorová L, Šlajcherová K, Seifertová D, Fischer L, Fišerová J, Petrášek J, Schwarzerová K. Arp2/3 complex subunit ARPC2 binds to microtubules. Plant Sci. 2015;241:96–108.View ArticlePubMedGoogle Scholar
- Chan J, Calder G, Fox S, Lloyd C. Cortical microtubule arrays undergo rotary movements in Arabidopsis hypocotyl epidermal cells. Nat Cell Biol. 2007;9:171–5.View ArticlePubMedGoogle Scholar
- Burian A, Ludynia M, Uyttewaal M, Traas J, Boudaoud A, Hamant O, Kwiatkowska D. A correlative microscopy approach relates microtubule behaviour, local organ geometry, and cell growth at the Arabidopsis shoot apical meristem. J Exp Bot. 2013;64:5753–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Wiesler B, Wang QY, Nick P. The stability of cortical microtubules depends on their orientation. Plant J. 2002;32:1023–32.View ArticlePubMedGoogle Scholar
- Deeks MJ, Fendrych M, Smertenko A, Bell KS, Oparka K, Cvrčková F, Žársky V, Hussey PJ. The plant formin AtFH4 interacts with both actin and microtubules, and contains a newly identified microtubule-binding domain. J Cell Sci. 2010;123:1209–15.View ArticlePubMedGoogle Scholar
- Liu Z, Schneider R, Kesten C, Zhang Y, Somssich M, Zhang Y, Fernie AR, Persson S. Cellulose-microtubule uncoupling proteins prevent lateral displacement of microtubules during cellulose synthesis in Arabidopsis. Dev Cell. 2016;38:305–15.View ArticlePubMedGoogle Scholar
- Grefen C, Donald N, Hashimoto K, Kudla J, Schumacher K, Blatt MR. A ubiquitin-10 promoter-based vector set for fluorescent protein tagging facilitates temporal stability and native protein distribution in transient and stable expression studies. Plant J. 2010;64:355–65.View ArticlePubMedGoogle Scholar
- Lowry R. VassarStats: Website for statistical computation. http://vassarstats.net. Accessed 18 Aug 2016.
- Stangroom J. Social science statistics. http://www.socscistatistics.com. Accessed 18 Aug 2016.
- Free math help–MathCracker. http://www.mathcracker.com. Accessed 18 Aug 2016.
- McDonald JH. Handbook of biological statistics: Kruskal–Wallis test. http://www.biostathandbook.com/kruskalwallis.html. Accessed 18 Aug 2016.
- Spitzer M, Wildenhain J. BoxPlotR: a web-tool for generation of box plots. http://boxplot.tyerslab.com. Accessed 18 Aug 2016.