The MATLAB-based, open source software package, plusTipTracker, can be used to analyze image series of fluorescently-labeled +TIPs to quantify microtubule dynamics.
Microtubule (MT) plus-end-tracking proteins (+TIPs) localize to the growing plus-ends of MTs and regulate MT dynamics1,2. One of the most well-known and widely-utilized +TIPs for analyzing MT dynamics is the End-Binding protein, EB1, which binds all growing MT plus-ends, and thus, is a marker for MT polymerization1. Many studies of EB1 behavior within growth cones have used time-consuming and biased computer-assisted, hand-tracking methods to analyze individual MTs1-3. Our approach is to quantify global parameters of MT dynamics using the software package, plusTipTracker4, following the acquisition of high-resolution, live images of tagged EB1 in cultured embryonic growth cones5. This software is a MATLAB-based, open-source, user-friendly package that combines automated detection, tracking, visualization, and analysis for movies of fluorescently-labeled +TIPs. Here, we present the protocol for using plusTipTracker for the analysis of fluorescently-labeled +TIP comets in cultured Xenopus laevis growth cones. However, this software can also be used to characterize MT dynamics in various cell types6-8.
The goal of this method is to obtain quantitative information regarding microtubule (MT) plus-end-tracking protein (+TIP) dynamics in living growth cones. MT +TIPs are a group of proteins that localize to the plus-ends of MTs9,10. They perform a range of functions to regulate parameters of MT dynamic instability11, including rates of polymerization, catastrophe, and rescue. One well-used method for analyzing MT dynamics is to track the behavior of the +TIP EB1, which binds specifically to growing MT plus-ends1,12. EB1 is known to recruit several other proteins to growing MT plus-ends13,14, and has recently been established as a MT maturation factor15, promoting both MT growth and catastrophe frequency15,16.
Many studies of MT dynamics within growth cones have utilized hand-tracking methods to measure changes in EB1-GFP dynamics over time1-3, as EB1 localization to MT plus-ends can be used as a marker for MT polymerization. A key benefit for examining EB1-GFP comets as a proxy for MT growth is that MT dynamics can be measured even in regions of significant MT overlap. While the method of hand-tracking EB1-GFP comets has provided useful insights into MT behaviors1-3, it is time-consuming and can be biased. Additionally, as aberrant growth cone behaviors are likely the result of minute shifts in cytoskeletal dynamics, analyzing only a small subset of MTs (usually necessary when hand-tracking) may miss significant information.
Thus, we measure global MT dynamics parameters using the software package, plusTipTracker4, after the acquisition of high-resolution, live images of tagged EB1 in cultured embryonic growth cones5. This software, developed in the Danuser Lab, has been used in several studies characterizing MT dynamics in various cell types6-8. It is an open-source, user-friendly, MATLAB-based package that includes automated detection, tracking, visualization, and analysis for movies of fluorescently-labeled +TIPs. A long list of specific parameters of MT dynamics are calculated by this software (see Reference 4 for details), but for analysis of MT dynamics in growth cones, the most useful parameters are MT growth track velocity (in microns/minute), growth track lifetime (in seconds), and growth track length (in microns). The software can be downloaded directly from the Danuser Lab website (under “Software”). While the Danuser Lab currently supports a newer interface for +TIP tracking analysis, which is incorporated into a software package called u-track 2.0, the original, stand-alone software will remain available. The underlying algorithms between the two programs are the same (at least as of 2014), with only a difference of interface and analysis outputs. For the novice user with little MATLAB and/or computational analysis experience, plusTipTracker has more user-friendly features, including automated statistical parameter outputs.
Here, we describe the steps for analyzing images of EB1-GFP dynamics in cultured Xenopus laevis growth cones. This protocol was utilized in a recent paper examining MT dynamics17. See also Lowery et al. 20125 for detailed instructions regarding culturing growth cones expressing EB1-GFP. While this paper primarily focused on examining EB1-GFP dynamics in growth cones, the same protocol can be used for other cell types17. For all cell types, the time interval between frames should be between 0.5-2 sec for optimal +TIP tracking. A time interval of up to 4 seconds between frames is possible, but this increased interval time results in additional tracking errors.
This protocol and video are meant to serve as a companion to the original paper describing the software package in more detail4, as well as the Technical Report that comes with the software download on the Danuser Lab website. Readers are encouraged to review these documents carefully if there are additional questions regarding using the software.
1. Prior to Image Analysis
2. plusTipGetTracks
Note: The first step in image analysis is to detect the EB1-GFP comets, link the comets into tracks, and determine the parameters of microtubule dynamics. This is obtained with the command “plusTipGetTracks”4.
3. plusTipSeeTracks
Note: Now that the microtubule tracks have been defined, the function “plustipSeeTracks” is used for track visualization4. This function can provide multiple outputs for visualization, including spatial MT dynamics maps and speed movies, but here, the focus is solely on using “Track Movies” to display MT tracks superimposed on the growth cone images. While plusTipGetTracks can analyze multiple movies at a time, plusTipSeeTracks can only analyze one movie at a time.
4. plusTipGroupAnalysis
Note: This final function is used to create groups of movies for analysis and comparison of their MT track parameters.
Using this software as described here will provide several files of information that quantify +TIP dynamics in living cells.
The function plusTipGetTracks identifies the tracks (using example settings shown in Figure 1), and then provides parameters regarding the +TIP tracks. To view the information that the software has obtained, go into the roi_X directory that was created in step 2.2. The “feat” folder contains “overlayImages”, which is a series of images showing the detected comets. Examining these images using image analysis software can demonstrate the accuracy of comet detection. The “meta” folder also contains detailed information regarding +TIP comet statistics, including the “projData” file, as well as the “Stats” file. To view the “Stats” file, drag it into an open worksheet of a spreadsheet application. This file contains the calculated microtubule parameters for each movie (Figure 2). As noted above, a long list of specific parameters of MT dynamics are calculated by this software (see Reference 4 for details), including MT growth track velocity (in microns/minute), growth track lifetime (in seconds), and growth track length (in microns).
The function plusTipSeeTracks saves a movie of the tracked comets, which can be reviewed by opening the file “allTracks_X_X_X” (Figure 3).
The function plusTipGroupAnalysis combines multiple individual data sets into groups and creates folders (named perCell or pooledData, depending on which analysis is selected) that contain group parameter data, including histograms, plots, and spreadsheets for comparing groups and individual parameters within each group (Figure 4).
Figure 1. PlusTipGetTracks settings used for EB1-GFP comets in Xenopus laevis growth cones. This figure shows the specific settings for “Detection”, “Tracking”, and “Post-Processing” steps that may be used for the analysis of EB1-GFP comets in Xenopus laevis growth cones. The plusTipParamSweepGUI tool, available within the plusTipTracker package, may be used to optimize tracking settings for other model organisms and/or cell types7.
Figure 2. Screenshot of MT parameters obtained from plusTipGetTracks analysis. The “meta” folder, created by running plusTipGetTracks, contains information regarding +TIP comet statistics. By dragging the “Stats” file into a spreadsheet application, microtubule dynamics parameters can be examined.
Figure 3. Screenshot of MT track movie obtained from plusTipSeeTracks analysis. PlusTipSeeTracks not only allows for microtubule track visualization, but also serves as a verification tool by allowing the user to view the validity of the data acquired from plusTipGetTracks.
Figure 4. Screenshot of MT parameters obtained from plusTipGroupAnalysis. PlusTipGroupAnalysis offers the user a simple method for comparing groups and individual parameters between within each group by combining multiple individual data sets and generating statistical output, which can be examined in a spreadsheet application.
PlusTipTracker provides a straightforward, graphical user interface to quickly and automatically detect virtually all visible EB1-GFP comets in a cell or growth cone, link the comets into tracks, and calculate MT parameters. Other publications have reported the design of similar types of software (for example, Marx et al. also utilized quantitative analysis of tagged EB1 dynamics in growth cones18). But, this software appears to be unique in its ease of access, as it is freely downloadable from the website of the Danuser Lab, which specializes in designing open-source, turn-key software useful to the cell biological community. While access to MATLAB is required, one does not need to be fully familiar with this computer application in order to utilize the software. However, there are a few points that need to be addressed for ease of use.
First of all, one of the most common issues that arise when using the software and computer application for the first time is related to the file path. If this error occurs (with the prompt “Error using cd – Argument must contain a string. Error in formatPath…”), then the easiest solution is to ensure that the plusTipTracker software, as well as the directory with all of the “images” subdirectories, are both in the same “MATLAB” file path. It is best if these are not in the “Program Files” directory, as it has been suggested that the space in the “Program Files” name may be an issue. Related to this, it is important to note that plusTipTracker saves the file path that was utilized when first computing the plusTipGetTracks analysis, and as such, this file path must be maintained when this data is accessed and employed by another component of the software. The functions plusTipGetTracks, plusTipSeeTracks, and plusTipGroupAnalysis all use the original saved file path, and thus, trying to call those functions for a given movie, after moving the files to another path, will result in an error.
Another common error to occur during analysis is when Tracking fails during the plusTipGetTracks step. This will occur if a frame in the image series contains no detectable comets. This will completely halt the analysis and no post-processing will occur. An easy fix to circumvent this issue and allow the analysis to proceed, is to create a mock comet on the image in an area where it will not be incorrectly linked to any actual tracks. This will not impact the final track parameters, as any comet that does not exist in a minimum number of consecutive frames will be filtered out of the final analysis.
One other issue that may arise is faulty comet detection. This can usually be fixed by improving the region-of-interest selection in step 2.3. It is important to draw the region of interest closely around the cell and not to draw a wider region than is necessary. The software uses this region to determine the background used during comet detection. If comet detection is still sub-optimal with the default settings, the settings can be adjusted in the plusTipGetTracks window (during step 2.7).
After any analysis, it is critical to validate the automated track linkages by eye, using plusTipSeeTracks. Tracking settings may need to be modified to reduce the numbers of false positive or false negative comet linkages. See the original plusTipTracker documentation4 as well as the Technical Report PDF that accompanies the software download for details on optimizing the settings. The performance of this software compared to hand-tracking has previously been tested in non-neuronal cells4. Growth cones pose a slightly different challenge, however, as growth cone MTs exhibit frequent translocations in all directions17, in addition to MT growth and shrinkage. One issue that was not found to be a major concern is whether the tightly packed MTs in growth cones pose tracking difficulties17. As only a subset of MTs are in the growing phase, with EB1-GFP on the ends, resolving and tracking individual EB1-GFP comets was not problematic. However, it should be noted that these previous studies used Xenopus laevis growth cones, which were specifically chosen because of their relatively large growth cone size (approximately 10 microns), compared to other vertebrate growth cones. Using these larger growth cones allows for more accurate EB1-GFP comet analysis.
To assess the utility and accuracy of this software for analyzing EB1-GFP tracks in Xenopus laevis growth cones, we compared the experience of using plusTipTracker with hand-tracking of an identical data series (data not shown). The time it took to hand-track EB1-GFP comets in an average growth cone of 39 comet tracks (in a 1 minute time-lapse series, with 2 seconds between each frame) was over two hours, compared to two minutes with the software. The parameters obtained with the two methods were similar for MT growth velocity (7.4 microns per minute for automated tracking versus 7.0 microns per minute for hand-tracking). However, for growth lifetime and length, software analysis leads to significantly shorter tracks (by about half the time and distance). This is due to growth tracks being split by the software if a comet goes in and out of focus over time. While the human eye can easily identify that it is the same comet, the software does not. This issue is not problematic though, if one is using the software in order to compare multiple conditions. Since identical tracking parameters are used for all conditions (and assuming that comets go in and out of focus at the same rate in multiple conditions), then the relative lifetimes and lengths are still quite useful measurements for comparison. As for automated analysis error rates, these depend greatly on the quality of the images. In high signal-to-noise movies, the percent of misjoined or incorrect tracks is in the single digits. Even in lower-quality movies (where individual comets are still clearly visible by eye, but the background noise is greater), the error rates are still low enough (5-15%) that the significant time saved by using the software is worth the costs in error. This is especially the case when analyzing hundreds of growth cones (sixty to eight growth cones per condition were analyzed in a previous study17).
It is important to note that this software was designed for detecting +TIP comets that only bind to growing MT ends, such as EB1-GFP. Given that the linkage and tracking algorithms expect that comets only exist on polymerizing MTs, using this software to analyze the dynamics of a fluorescently-tagged +TIP that binds to shrinking MT ends in addition to growing ends will lead to incorrect information regarding calculated MT growth velocities.
One of the unique features of this software compared to other single-particle-tracking software, is that it takes into account known MT behaviors to calculate not only polymerization parameters, but also shrinkage parameters. It does this by linking an EB1-GFP comet that has disappeared with one that has newly formed directly behind it in the same trajectory (this is called a backgap, or bgap track). While this algorithm works well for some cell types, such as HeLa cells4, it is a less effective feature when analyzing MT dynamics in growth cones. This is because MT tracks frequently follow each other along the exact same paths in growth cones (often following along F-actin bundles), and so it is usually impossible to tell if bgap linkages are correct. For this reason, it is not recommended to utilize the bgap data outputs in growth cones.
Despite these minor caveats and issues which must be taken into consideration when using plusTipTracker (and most any automated image processing program), this software can be a very useful tool for analyzing thousands of EB1-GFP comets in a relatively short amount of time.
The authors have nothing to disclose.
We thank Dr. Gaudenz Danuser and members of his lab for creating the plusTipTracker software and for helpful discussion regarding using the software, in particular Maria Bagonis and Sebastien Besson. We especially thank the Boston College Media Center for their assistance and support in the creation and editing of the video. We also thank members of the Lowery Lab for useful discussions and constructive criticism, and Abigail Antoine for proof-reading the manuscript. This work was funded by an NIH R00 MH095768 award to LAL.
plusTipTracker software | Danuser Lab | http://lccb.hms.harvard.edu/software.html | This software may be hosted by another website in the future. If the listed site does not exist, search "Danuser Lab Software" on a web search engine to find the site. |
Matlab software | Mathworks | http://www.mathworks.com/products/matlab/ |