PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Cytometry A. Author manuscript; available in PMC 2010 April 1.
Published in final edited form as:
PMCID: PMC2661008
NIHMSID: NIHMS76355

Automated Analysis of NeuronJ Tracing Data

Abstract

Background

Studies of neuronal differentiation in vitro often involve tracing and analysis of neurites. NeuronJ (1,2) is a program that can be used for semi-automated tracing of individual neurons; when tracing is completed, a text file containing neurite length measurements is generated. Using cultured hippocampal neurons, we have found that to reach statistical significance it is generally necessary to trace about 100 neurons in each treatment group. Post tracing data analysis requires importing each text file into a statistics program. Analysis of distinct parameters, such as effects of a treatment on axonal versus dendritic branching, requires a great deal of time consuming post tracing data manipulation.

Methods

We have developed XL_Calculations, a Java based program that performs batch analysis on NeuronJ measurement files and automatically makes multiple calculations, including the number, length and total output (sum length) of primary, secondary and tertiary neurites on axons and dendrites and writes the calculations into an Excel worksheet.

Results and Conclusion

Batch processing of NeuronJ measurement files dramatically reduces the time required to analyze neuronal morphology. In addition, our program performs more than 45 distinct calculations, enabling detailed determination of treatment effects on neuronal differentiation. Using this program to analyze NeuronJ tracing data, we demonstrate that continuous exposure of differentiating hippocampal neurons to Netrin 1 increases the number of secondary branches on both axons and dendrites, without significantly altering the length of the axon, dendrites or branches. Similar results were obtained whether neurons were grown on poly-d-lysine or laminin.

Keywords: neurite tracing, NeuronJ, ImageJ, axon guidance, growth cone, Netrin, laminin

Background

Studies of neuronal development and regeneration often use primary neuronal cultures or neuronal cell lines to test the effect of a given molecule on the formation, elongation and branching of axons and dendrites. Typically, cultures are analyzed by immunofluorescence and the number and lengths of the neurites are quantified by tracing of the fluorescent images. Tracing can be done manually using common image analysis software (e. g. ImageJ (3), Adobe Photoshop™), but this is very labor intensive. Fully automated programs are available (e.g. Neurolucida (4) or HCA-Vision (5,6)) and are much less labor intensive, but are often cost prohibitive. An alternative is to use the semi-automated tracing program NeuronJ.

NeuronJ (1) is distributed as a plugin for ImageJ, free image analysis software distributed by the National Institutes of Health (3). Semi-automated tracing with NeuronJ has been shown to be as accurate as both fully manual (1) and fully automated tracing (5). Using NeuronJ, fluorescently labeled neurons are manually traced, but an algorithm compares the pixel intensity on the neurite with its adjacent pixel neighborhoods, automatically updating the cursor to follow an estimated path of the neurite and increasing the accuracy and speed of tracing. After the neurites are traced, a text file is generated that contains measurements of the lengths of all the neurites.

In a typical experiment to determine the effects of a given treatment on neurite outgrowth, we find it is necessary to trace about 100 neurons in each treatment group to obtain statistical significance. Thus, for a simple experiment with a single treatment and control, at least 200 measurement files would be generated. To analyze the morphology of individual neurons, the user must import each measurement text file into a database program such as Excel for statistical analyses. Because of all the manipulations involved, this post-measurement analysis can be tedious and error prone.

Here we report development of a software that performs batch analysis on NeuronJ measurement files and makes multiple calculations for a typical cell culture experiment, including the number, length and total output (sum length) of primary, secondary and tertiary neurites. These data are subsequently exported to a single Excel worksheet. Use of this program dramatically reduces the amount of time required to analyze NeuronJ tracing data and facilitates analysis of multiple morphometric parameters.

We have used the XL_Calculations program to analyze the effects of Netrin 1 on the differentiation of embryonic hippocampal neurons and report that treatment with Netrin 1 increases the number of axonal and dendritic branches, without significantly changing total neurite output. In addition, similar results were obtained whether neurons were grown on poly-d-lysine or laminin, indicating that unlike the growth cone guidance effects, the effects of Netrin 1 on branching are not substrate dependent.

Methods

Hippocampal cultures and immunofluorescence

Hippocampal cultures were prepared from embryonic day 16 mice as previously described (7) and plated on coverslips coated with either 1.0 mg/ml poly-d-lysine (PDL) or a mixture of 0.1 mg/ml PDL and 4 µg/ml laminin. Netrin 1 (R&D Systems) was applied at 250 ng/ml and left on continuously. For cultures analyzed in Fig. 1, replication deficient recombinant adenoviruses were produced using the AdEasy system (8) and were used to express enhanced green fluorescent protein (EGFP) in a subset of the neurons. After 72 hrs, cultures were fixed for 30 min. at either 37°C or 4°C in PPS (4% paraformaldehyde in PHEM buffer (60 mM PIPES pH 7.0/25 mM HEPES pH 7.0/10 mM EGTA/2 mM MgCl2) with 0.12M sucrose); the 4°C fix was used for staining with Tau and MAP2, the 37° fix was used for all other cultures. After rinsing in PBS, coverslips were incubated in 3% fatty acid free bovine serum albumin (BSA) in PBS for 30 min, permeablized for 10 min. in 0.2% triton/PBS, rinsed, and re-blocked in 3% BSA/PBS for 30 min. Primary antibodies included anti-Tau (Chemicon), anti-MAP2 (Covance) and anti-Tubulin Beta3 (Promega). Incubations with primary and secondary antisera were done in the presence of 1% BSA/PBS and coverslips were mounted with 2.5% 1,4-Diazabicyclo-[2.2.2]Octane/150 mM Tris pH 8.0/30% glycerol to reduce photo bleaching. Cells were imaged on a Zeiss Axiovert 200. Neurites were identified by staining with EGFP and anti-Tubulin Beta3 (Fig. 1A-E) or with a combination of anti-MAP2 and anti-Tau (Fig 1F, 2A-D).

Figure 1
Analysis of NeuronJ tracing data reveals altered neuronal morphologies
Figure 2
Effects of Netrin 1 on axonal and dendritic branching in hippocampal neurons grown on Poly-D-lysine (PDL) or laminin

NeuronJ tracing and labeling

Neurites were traced using NeuronJ following the developers' instructions. After the neurites are traced, they were labeled as primary (emanating directly from the soma), secondary (branching from a primary) or tertiary (branching from a secondary). The neurites could be labeled as axon or dendrite; however, assignment of the axon/dendrite labels requires the use of fluorescent markers that are specific to axonal and dendritic regions within the cell. Since many users are tracing neurons that are either labeled with an antibody to neuronal specific Tubulin Beta3 (which labels all neurites) or expressing a fluorescent protein such as EGFP, it is not possible to accurately assign axon and dendrite labels based on staining. Therefore, our program use the primary/secondary/tertiary labels to perform calculations.

Once labels were assigned, the neurite tracings appeared color coded by type (Fig. 2A-E). Neurites were then assigned to clusters, with each cluster comprised of a primary neurite and all its associated branches. After tracing was completed and clusters had been assigned, a text file containing neurite length measurement data was generated for each neuron traced and a snapshot of the tracings overlaid on the neuron was saved as a TIFF file (e.g. Fig. 1A-E and 2A-D). For the analysis Table 2, 124-145 neurons from each treatment group were traced. Because some of these were stage 2 neurons, the analysis of stage 3 neurons (Fig. 2) was performed on 95-112 neurons.

Table 2
Percent of Neurons at Stage 2 vs. Stage 3

Use of XL_Calculations for morphometric analysis

The XL_Calculations program was installed by copying a folder containing the Jar program and an Excel template file onto the computer. The program was used to perform arithmetic calculations on the NeuronJ measurement data and write the calculations directly into an Excel worksheet, with measurements for each neuron in a different row. Detailed step-by-step instructions for using the XL_Calculations program is provided in Supplemental Fig. 1 and Supplemental methods.

The program was created with Java Development Kit version 1.5.0_13 and requires installation of Java version 1.5 (http://java.sun.com/) or higher on the native operating system and the third party driver Jxl for Excel (http://jexcelapi.sourceforge.net/). This driver is required for Java Database Connectivity (JDBC) import and is included in our software package. The XL_Calculations program also requires Excel from Microsoft Office version 2003 (Windows) or 2004 (Mac) or later and has been tested on both the Mac and Windows versions of Excel. The XL_Calculations program, source code, manual, and parameter customization are available upon request.

Results and Discussion

The steps required to run the XL_Calculations program are outlined in Supplemental Fig. 1 and detailed in the Supplemental Methods. The XL_Calculations program performs more than 45 different calculations pertaining to the length and branching of neurites (Table 1). The calculations performed by our program are designed to distinguish between neurons at different stages of differentiation; stage 1, filopodia rich lamellapodium with no identifiable neurites; stage 2, neuron has neurites, but no identifiable axon; stage 3, neuron has a rapidly elongating axon; stage 4, dendrites elongate; stage 5, synapse formation (9). Stage 2 and 3 neurons can be distinguished biochemically and morphologically (9). The XL_Calculations program uses established morphometric criteria (10) to determine if a neuron is stage 2 or 3; a neuron is defined as stage 3 if its longest neurite is at least twice as long as the second longest neurite (Fig. 1B), otherwise the neuron is classified as stage 2 (Fig. 1A). Once the stage of the neuron is determined, stage-specific calculations are performed. For stage 2 neurons, the values are reported for all neurites. For stage 3 neurons, separate calculations are made for axons and dendrites.

Table 1
Calculations made by the XL_Calculations Program

The arithmetic calculations made by our program could be performed manually (e.g. average axon length or average number of dendrites) and the values obtained by manual calculation and the XL_Calculations program are identical (data not shown); however, XL_Calculations performs multiple calculations on a batch of files, greatly reducing the processing time. For example, in a timed trial, it took approximately 15 minutes to manually calculate the average total neurite output from 25 NeuronJ measurement files; this included time to open each file, copy the data into an Excel worksheet, determine the total neurite output for each neuron and finally determine the average total neurite output for the population. Importantly, this manual calculation did not distinguish stage 2 vs. 3 neurons. In contrast, the XL_Calculations program made all 45 calculation in under 40 seconds; determining the average total neurite output for stage 2 vs. 3 neurons was a simple matter of sorting and averaging using the appropriate Excel functions; this entire procedure took less than 2 minutes.

By separating the calculations for stage 2 and stage 3 neurons, it is possible to determine if a given treatment differentially affects axonal and dendritic outgrowth. In addition, treatments that affect neuronal polarization (i.e. the transition from stage 2 to stage 3) are revealed because they alter the percentage of cells at each stage. A treatment that blocks or delays axon formation will increase the percentage of stage 2 cells. Sometimes, a treatment results in the formation of multiple axons; such a situation increases the frequency of stage 2 categorizations because the longest neurites are not twice as long as the second longest (Fig. 1C); however, comparison of the lengths of the second longest neurites in treatments and the control reveals that both are significantly longer than normal dendrites, suggesting abnormal polarity.

XL_Calculations also determines the number and lengths of secondary and tertiary neurites (i.e. branches) and the total neurite output. When analyzing total neurite output, it is informative to know the developmental stage and the contributions from axons, dendrites and branches. For example, one treatment may enhance axon elongation (Fig. 1D), while another treatment reduces elongation of the primary neurites, but increases branching (Fig. 1E); these likely reflect different underlying biological mechanisms of action. Thus, the calculations made by our program reveal nuances in neuronal morphology not evident from the ‘total neurite output’ calculation alone.

While our program is optimized for analysis of stage 2 and 3 neurons in the first few days of differentiation (generally 2-5 days in culture), it can also be used for analysis of stage 4 neurons, in which dendrite elongation has begun. When analyzing more mature cultures, the user should empirically establish an average dendrite length for stage 3 neurons, then set a threshold for dendrite length that distinguishes stage 3 and stage 4 neurons. It is important to note that results based on morphometric analyses should always be confirmed by immunohistochemical methods. For example, staining with antibodies to the microtubule binding proteins Tau and MAP2 can be used to identify axons and dendrites, respectively (Fig 1F).

To validate the usefulness of the XL_Calculations program, we have analyzed the effect of Netrin 1 on the differentiation of cultured embryonic hippocampal neurons. Netrin 1 was first identified as a guidance molecule that attracts commissural axons toward the floor plate of the spinal cord (11). Subsequently, cultured Xenopus retinal growth cones were shown to turn toward a gradient of Netrin 1 (12). Interestingly, the attraction to Netrin 1 was converted to repulsion when Xenopus spinal cord neurons were grown on laminin (13). This observation lead to the speculation that laminin might reverse other effects of Netrin 1.

Laminin is an extracellular matrix molecule that is found in many regions of the developing CNS, including the hippocampus and spinal cord, but which is not present in most layers of the cortex (14,15). Laminin has been shown to selectively enhance axonal growth and accelerate the development of polarity by several types of CNS neurons, including hippocampal and cervical ganglion neurons in culture (16,17). When cortical neurons are cultured on polylysine alone (in the absence of laminin), exposure to Netrin 1 increases the number and length of axonal branches, without significantly changing the length of the primary axon (18). Under similar culture conditions, short-term exposure of hippocampal neurons to Netrin 1 induces formation of axonal shaft filopodia (18,19).

Here we investigate whether Netrin 1 enhances branching in hippocampal neurons and determine if this effect is altered by growth on laminin. As previously shown (16), growth on laminin significantly increased the development of polarity, as reflected by the increased percent of stage 3 neurons on laminin compared to PDL alone (Table 2). Unexpectedly, treatment with Netrin-1 appeared to slightly delay development of polarity for neurons grown on PDL alone, but not for neurons grown on laminin (Table 2). Because of the relatively small number of stage 2 neurons at this time point (72 hrs in vitro), further analysis focused on stage 3 neurons with identifiable axons and dendrites. For stage 3 neurons, growth on laminin increased the number of dendrites (Fig. 2G), axon length (Fig. 2N) and the number and length of dendritic and axonal branches (Fig. 2E-F, 2I-M), resulting in an increase in total neurite output compared to neurons grown on PDL alone (Fig. 2O).

When neurons grown on PDL were treated with Netrin 1, there was a significant increase in the number of secondary branches on dendrites (Fig. 2E) and axons (Fig. 2I) and a slight, but not statistically significant, increase in the length of secondary branches. Treatment of neurons grown on laminin with Netrin I also resulted in increased branching; however, on laminin only the increase in axonal branches was statistically significant (Fig. 2I, L; tertiary dendritic branches are not reported because there were too few (<0.05 3° branches/dendrite) for statistical analysis). On both PDL and laminin, Netrin 1, increased the frequency of neurons with multiple secondary axonal branches (Fig. 2K).

Surprisingly, on both PDL and laminin, treatment of neurons with Netrin 1 resulted in a slight decrease in the length of primary dendrites (Fig. 2H) and axons (Fig. 2N) and decreased the length of secondary branches on laminin (Fig. 2F, J). Although the decrease in neurite lengths was not always statistically significant, it apparently balanced the increase in the number of branches, resulting in no net increase in total neurite output in Netrin 1 treated neurons compared to controls on either PDL or laminin (Fig. 2O).

Together, our results indicate that, contrary to predication, laminin does not block the ability of Netrin 1 to induce branching of hippocampal neurons. In addition, we demonstrate for the first time that treatment with Netrin 1 increases dendrite branching. The observation that the effects of Netrin 1 were more significant of PDL than on laminin may indicate that, on laminin, the neurons are already close to maximizing their neurite output. It is interesting to note that in hippocampal neurons, Netrin 1 does not appear to significantly alter branch length or total neurite output, where as in cortical neurons, Netrin 1 increase both branch length and total neurite output (18). These observations suggest that Netrin 1 may differentially regulate arborization of cortical and hippocampal neurons in vivo.

In conclusion, use of the XL_Calculations program with NeuronJ greatly reduces the time required to analyze tracing data; analysis of a batch of 50-100 measurement files takes only a few minutes. These calculations can distinguish among treatments that differentially alter neuronal polarity, growth of primary neurites (i.e. axons and dendrites) and branching of axons and dendrites. We have used these detailed calculations to reveal the differential effects of exposure to Netrin 1 on hippocampal axonal and dendritic outgrowth and branching. It is our hope that by facilitating these calculations, our program will make it easier for other investigators to perform a sophisticated analysis of multiple morphometric parameters pertaining to neuronal differentiation.

Acknowledgements

We would like to thank Teresa Nick for critical reading of this manuscript. This research was supported by NIH RO1NS049178 to L.M.L.

Supplementary Material

Supp Fig 1

Supp Fig Leg

Supp Methods

References

1. Meijering E, Jacob M, Sarria JC, Steiner P, Hirling H, Unser M. Design and validation of a tool for neurite tracing and analysis in fluorescence microscopy images. Cytometry A. 2004;58(2):167–76. [PubMed]
5. Vallotton P, Lagerstrom R, Sun C, Buckley M, Wang D, De Silva M, Tan SS, Gunnersen JM. Automated analysis of neurite branching in cultured cortical neurons using HCA-Vision. Cytometry A. 2007;71(10):889–95. [PubMed]
7. Lanier LM, Gates MA, Witke W, Menzies AS, Wehman AM, Macklis JD, Kwiatkowski D, Soriano P, Gertler FB. Mena is required for neurulation and commissure formation. Neuron. 1999;22(2):313–25. [PubMed]
8. He TC, Zhou S, da Costa LT, Yu J, Kinzler KW, Vogelstein B. A simplified system for generating recombinant adenoviruses. Proc Natl Acad Sci U S A. 1998;95(5):2509–14. [PubMed]
9. Dotti CG, Sullivan CA, Banker GA. The establishment of polarity by hippocampal neurons in culture. J Neurosci. 1988;8(4):1454–68. [PubMed]
10. Goslin K, Banker G. Rat hippocampal neurons in low density culture. In: Banker G, Goslin K, editors. Culturing Nerve Cells. 2nd edition ed. MIT Press; Cambridge, MA: 1998. pp. 339–370.
11. Kennedy TE, Serafini T, de la Torre JR, Tessier-Lavigne M. Netrins are diffusible chemotropic factors for commissural axons in the embryonic spinal cord. Cell. 1994;78(3):425–35. [PubMed]
12. de la Torre JR, Hopker VH, Ming GL, Poo MM, Tessier-Lavigne M, Hemmati-Brivanlou A, Holt CE. Turning of retinal growth cones in a netrin-1 gradient mediated by the netrin receptor DCC. Neuron. 1997;19(6):1211–24. [PubMed]
13. Hopker VH, Shewan D, Tessier-Lavigne M, Poo M, Holt C. Growth-cone attraction to netrin-1 is converted to repulsion by laminin-1. Nature. 1999;401(6748):69–73. [PubMed]
14. Dent EW, Kwiatkowski AV, Mebane LM, Philippar U, Barzik M, Rubinson DA, Gupton S, Van Veen JE, Furman C, Zhang J. Filopodia are required for cortical neurite initiation. Nat Cell Biol. 2007;9(12):1347–59. [PubMed]
15. Zhou FC. Four patterns of laminin-immunoreactive structure in developing rat brain. Brain Res Dev Brain Res. 1990;55(2):191–201. [PubMed]
16. Lein PJ, Banker GA, Higgins D. Laminin selectively enhances axonal growth and accelerates the development of polarity by hippocampal neurons in culture. Brain Res Dev Brain Res. 1992;69(2):191–7. [PubMed]
17. Lein PJ, Higgins D. Laminin and a basement membrane extract have different effects on axonal and dendritic outgrowth from embryonic rat sympathetic neurons in vitro. Dev Biol. 1989;136(2):330–45. [PubMed]
18. Dent EW, Barnes AM, Tang F, Kalil K. Netrin-1 and semaphorin 3A promote or inhibit cortical axon branching, respectively, by reorganization of the cytoskeleton. J Neurosci. 2004;24(12):3002–12. [PubMed]
19. Strasser GA, Rahim NA, VanderWaal KE, Gertler FB, Lanier LM. Arp2/3 is a negative regulator of growth cone translocation. Neuron. 2004;43(1):81–94. [PubMed]