Analysis of axonal trafficking via a novel live-imaging technique reveals distinct hedgehog transport kinetics

ABSTRACT The Drosophila melanogaster (Dmel) eye is an ideal model to study development, intracellular signaling, behavior, and neurodegenerative disease. Interestingly, dynamic data are not commonly employed to investigate eye-specific disease models. Using axonal transport of the morphogen Hedgehog (Hh), which is integral to Dmel eye-brain development and implicated in stem cell maintenance and neoplastic disease, we demonstrate the ability to comprehensively quantify and characterize its trafficking in various neuron types and a neurodegeneration model in live early third-instar larval Drosophila. We find that neuronal Hh, whose kinetics have not been reported previously, favors fast anterograde transport and varies in speed and flux with respect to axonal position. This suggests distinct trafficking pathways along the axon. Lastly, we report abnormal transport of Hh in an accepted model of photoreceptor neurodegeneration. As a technical complement to existing eye-specific disease models, we demonstrate the ability to directly visualize transport in real time in intact and live animals and track secreted cargoes from the axon to their release points. Particle dynamics can now be precisely calculated and we posit that this method could be conveniently applied to characterizing disease pathogenesis and genetic screening in other established models of neurodegeneration.


INTRODUCTION
Since the discovery of the first eye-specific mutant (white) in 1910, the Drosophila melanogaster (Dmel) visual system has been the focus of innumerable genetic screens, ranging from research in development, intracellular signaling, behavior, and in particular, neurodegenerative disease. In addition to the practical implications, Dmel disease models are popular because many (>75%) genes associated with human disease are conserved between fly and human and the basic fundamentals of cell biology are similar between these two species. In particular, the developing visual system has been a model of choice since the fly eye is amenable to genetic disruption and is dispensable for the organism's survival (Lenz et al., 2013;Sang and Jackson, 2005).
Although many eye-specific models of neurodegenerative disease have been developed (Alzheimer's disease, Parkinson's disease and trinucleotide expansion disorders such as Huntington's disease), these studies have mainly used static readouts of phenotype (e.g. biochemistry, immunohistochemistry, pathology, rescue of 'rough-eye') (Lenz et al., 2013;Sang and Jackson, 2005). Only a select few of these models have focused on dynamic data such as the velocity, flux, and positional distribution of pathogenically relevant proteins (e.g. PINK1, polyglutaminecontaining proteins, Lis1) or organelles in these mutated tissues (e.g. mitochondria) (Liu et al., 2000, 1;Stowers et al., 2002;Wang et al., 2011;Wyan-Ching Mimi Lee et al., 2004). This disparity between static and dynamic data is significant in the field of development especially with respect to the morphogen Hedgehog (Hh), which is trafficked down Dmel photoreceptor axons and is integral to eye-brain development. Amongst Hh-expressing tissues (Dmel embryo and wing disc, planarian CNS regeneration, and vertebrate central nervous system (CNS), limb, organs, and foregut) developing photoreceptor axons offer a unique and accessible means to study Hh intracellular trafficking (Jiang and Hui, 2008;Yazawa et al., 2009).
Hedgehog (Hh) is a highly conserved secreted morphogen capable of patterning many developing tissues. More recently, Hh signaling has been implicated in both Dmel and mouse adult stem cell maintenance and human neoplastic disease (Bale, 2002;Ihrie et al., 2011;Jiang and Hui, 2008;Rojas-Ríos et al., 2012;Scales and de Sauvage, 2009;Varjosalo and Taipale, 2008). While a number of papers have investigated Hh extracellular movement and downstream signaling, little is known about the mechanisms regulating its intracellular transport and secretion. To investigate the nature and factors governing Hh transport we developed a novel method to directly visualize and characterize this process in the photoreceptor axons of live and intact Drosophila larvae.
Using visualization of Hh in Dmel photoreceptors, we demonstrate the ability to comprehensively quantify the movement and directionality of important cargoes during eye and brain development in a practical and accessible model system. To our knowledge, this is the first instance of dynamic Hh data from neurons that has comprehensively analyzed its transport.
Unlike most model systems of axonal transport, which involve cultured neurons, our data are collected from intact organisms in which the normal environment and interactions with other cell types are preserved. While this Hh system could be applied to many cell biological questions (e.g. protein and/or organelle transport throughout disease progression) we focus on the apical/basal transport of Hh in developing photoreceptors in healthy and disease models, as this is particularly difficult to observe and quantify in other established models of Hh in development (Jiang and Hui, 2008;Yazawa et al., 2009). More broadly, this method could be adapted to additional neurodegenerative disease models (e.g. in the Dmel eye) and should enable future research to characterize the role of pathogenic proteins and organelles, candidate drug treatment, apical/basal cell polarity, and axon transport in development in any practical, transparent, and genetically tractable organism.

RESULTS
In fly larval photoreceptor neurons the developmental signal Hh is guided to two receptive fields; the apical (retina) and basal (growth cone, GC) ends where secretion of the morphogen is an inductive factor in photoreceptor differentiation and establishment of eyebrain neural connections ( Fig. 1A-C) (Huang and Kunes, 1996;Roignant and Treisman, 2009). Hh released apically induces ommatidial development (Fig. 1C,D, left), while Hh transported down the photoreceptor's axon (Fig. 1D, middle) and released in the brain induces the development and specification of postsynaptic target neurons (Fig. 1D, right). This phenomenon has also been described in retinal ganglion cells during rat visual system development (Beug et al., 2011;Dakubo et al., 2003Dakubo et al., , 2008Soukkarieh et al., 2007;Wallace and Raff, 1999;Wang et al., 2005). Thus, a balance of apical/basal transport to two receptive fields is integral to the precise timing and specification of the Drosophila eye and brain.
Hh undergoes autocleavage from its full-length form (HhU; 46 kDa) to become two products, a cholesterol-modified N-terminal signaling molecule (HhN; 19 kDa) and a C-terminal molecule (HhC; 24 kDa), the full role of which is still unknown (Lee et al., 1994). HhN contains all the information necessary for downstream signaling but when the C-terminal domain (HhC) is mutated or deleted from the full-length protein it is unable to undergo axonal Photoreceptor neurons differentiate temporally with the posterior-to-anterior progression (top, right to left) of the morphogenetic furrow across the eye disc. These neurons project their axons (R1-R8) into the brain through the optic stalk, where they spread to dorsal and ventral retinotopic positions (dorsal is up). Micrograph shows the Dmel eye from the lateral perspective. Image is an immunostaining against the eyespecific protein Chaoptin (24B10). Os, optic stalk; gc, growth cone (axon termini). (C) The Dmel visual system shown from the horizontal perspective. The R1-R6 axons terminate in the lamina (lam), while R7 and R8 axons terminate in the deeper medulla ganglion. Hh secreted from developing photoreceptor neurons is required for both eye and lamina (brain) development (top). Micrograph shows the Dmel eye from the horizontal perspective. As in (B), image is an immunostaining against Chaoptin (24B10). (D) Representative still frames from time-lapse movies of Dmel expressing ss::GFP::HhC (HhC) in developing photoreceptors (GMR-Gal4>UAS-ss::GFP::HhC, see transport. Under these circumstances all HhN is released at the apical (retina) side. The C-terminal domain, however, contains a conserved motif that designates the protein for axonal transport and release at the growth cone. Indeed, HhC (which travels with HhN) is necessary for the proper axonal transport of the N-terminal signaling domain (Chu et al., 2006).
To determine whether HhC undergoes axonal transport, we observed trafficking via real-time imaging of axons in live Drosophila larvae ( Fig. 2 and see Materials and Methods for stepby-step directions regarding data acquisition and analysis). Fluorescently tagged HhC (ss::GFP::HhC) was first expressed with CCAP-GAL4, which drives expression in a small population of motor neurons that extend axons along segmental nerves ( Fig. 3A) (Luan et al., 2006;Vömel and Wegener, 2007). The directionality of transport can be unambiguously determined because sensory axons sharing the nerve are not labeled (Wang et al., 2011). To graphically represent transport behavior we generated kymographs, which denote distance translocated on the x-axis versus time on the y-axis (Fig. 3B,C) (Miller and Sheetz, 2006). The slope of this trace yields velocity of transport. Such an analysis revealed that velocity of Hh transport was 0.31 µm s −1 in the anterograde direction, and 0.22 µm s −1 in the retrograde direction ( Fig. 3D, grey bars and Movie 1). A similar analysis on HhC expressed in photoreceptor neurons (GMR-Gal4; Fig. 3A,D, white bars and Movie 2; Tare et al., 2011) showed a larger disparity between anterograde and retrograde transport, with the velocities calculated at 0.50 and 0.17 µm s −1 , respectively. The ranges of velocities strongly suggest that Hh is trafficked within motor neuron and photoreceptor axons via microtubule-based fast axonal transport and that kinesin and dynein are its primary molecular motors (Goldstein and Yang, 2000). These depictions revealed significant anterograde bias of HhC particles, with 93.5% of total particles moving toward the motor neuron axon terminal and 86.6% of total particles moving anterogradely in photoreceptors ( Fig. 3E and Movies 1 and 2, respectively).
Visibly more Hh was transported in an anterograde direction than retrogradely (Fig. 3). Quantification analysis of photoreceptor neurons corroborates this, with anterograde flux being several-fold higher than retrograde flux (0.54 and 0.07 particles s −1 , respectively). Flux is defined as the average number of punctae that cross a certain position in the axon at any given time, and can be measured by counting the number of traces in any one direction that run across a line drawn through the kymograph (e.g. see Fig. 3C).
With this stark disparity in velocity and flux of anterograde versus retrograde particles we wondered whether HhC particles behaved differently depending on whether they were proximal (to the cell body), medial (in the optic stalk) or distal (near the growth cone) (Fig. 4A). While the percentage of anterograde particles was similar in the proximal and medial axon, this number decreases in the distal axon ( Fig. 4B and Movies 3 and 4). Indeed, the anterograde flux is most different between the proximal and distal axon (Fig. 4C, Movies 3 and 4). Focusing on just the proximal and distal axon, we found that anterograde HhC particles moved quickly in the proximal axon but then appear to slow down near the growth cone tip ( Fig. 4D and Movies 2-4). A plot of the velocities of these particles revealed a 'two-peak' distribution in the proximal axon (Fig. 4E, orange bars) with peaks in the ∼0.21-0.25 µm s −1 and ∼0.46-0.55 µm s −1 ranges, while the distal axon showed just one peak in the ∼0.26-0.30 µm s −1 range (Fig. 4E, white bars). These results suggest that there may be two potential modes of anterograde transport for HhC in the proximal axon which appear to resolve into one mode at the distal axon.
Finally, we wondered how Hh might behave in a mutant system. We chose to observe Hh in a recently developed model (in Drosophila and mice) of photoreceptor (PR) neurodegeneration, caused by mutations in fatty acid transport protein ( fatp) (Dourlen et al., 2012(Dourlen et al., , 2015. This protein is especially important to Hh processing since it controls the intracellular levels of palmitic acid, a lipid moiety which, in addition to cholesterol, is covalently linked to the mature Hh ligand (Buglino and Resh, 2012;Luiken et al., 1999;Seeßle et al., 2015). Although it is believed that accumulation of either Rhodopsin-1 protein or intracellular ceramide causes the acute adult-onset PR toxicity observed in fatp mutants, we wondered if altered Hh transport might also be contributing to the decline in PR health. Importantly, we also chose this model because PR neurodegeneration occurs in the adult eye, and thus the developing PRs are presumed healthy and should be free of any obvious toxic traits that might alter transport (e.g. apoptosis, loss of MT polymerization, damaged mitochondria).
Indeed, Hh flux was significantly different in the PRs of animals driving fatp RNAi with 28% of particles shifting to the retrograde direction of transport (ant/ret for: WT, 87%/13%; fatp RNAi, 59%/ 41%) (Fig. 5A). Although the retrograde velocity of these particles remained unchanged (WT, 0.17 µm s −1 ; fatp RNAi, 0.21 µm s −1 ), the anterograde velocity was dramatically different in the mutant condition (WT, 0.50 µm s −1 ; fatp RNAi, 0.23 µm s −1 ) (Fig. 5B). Possibly the most dramatic difference we observed was the nearly 4× increase in the number of particles moving in the retrograde direction (WT, 0.07 particles s −1 ; fatp RNAi, 0.26 particles s −1 ), while anterograde flux was not significantly different (Fig. 5C). Notably, no significant differences in Hh localization were observed with conventional immunohistochemistry (Fig. S1), which demonstrates the unique benefits of dynamic versus static data. Thus, despite PR death being reported exclusively in adult fatp mutants (Dourlen et al., 2012), we describe here, for the first time, that altered transport of free fatty acids (FFAs) could also alter Hh transport in the early stages of PR development.

DISCUSSION
The Drosophila melanogaster (Dmel) eye is a powerful system to study mechanisms of development, intracellular signaling, behavior, and neurodegenerative disease. Few studies have used dynamic data in the Dmel eye to characterize the pathology of neurodegeneration, in particular. By applying live imaging methods in photoreceptors and motor neurons, we were able to characterize transport of the morphogen Hedgehog (Hh), a protein integral to eye-brain development and implicated in stem cell maintenance and neoplastic disease. As an example, we demonstrate the facility of the Dmel visual system to derive unprecedented positional and temporal detail of complex phenotypes in real time. We posit that this method could be used to characterize the behavior of pathologically relevant cargoes or organelles in existing models of neurodegeneration or in screening for new phenotypes in disease models that are already characterized (e.g. PINK1, polyglutamine-containing proteins, Lis1) (Liu et al., 2000;Stowers et al., 2002;Wang et al., 2011;Wyan-Ching Mimi Lee et al., 2004).
In fly larval photoreceptor neurons the morphogen Hh is simultaneously released apically to activate ommatidial development and basally, down photoreceptor axons, to the brain to induce the specification of postsynaptic target neurons (Fig. 1). This balance of apical/basal transport to two receptive fields is integral to the precise timing and specification of the Dmel eye and brain, though little is known about the mechanisms governing this process. Unlike the other models of Hh in development Hh transport can be comprehensively quantified in various Dmel neuron types (Fig. 2). Our dynamic data of neuronal Hh suggests distinct transport kinetics and the possibility of two modes of Hh transport in the proximal axon, which appear to resolve distally (Figs 3 and 4). While little is known about the intracellular movement of Hh in cells that produce the ligand, our bimodal velocity analysis in the proximal axon may be indicative of differential trafficking and secretion of 'long-range' and 'short-range' Hh complexes, or perhaps of 'early' and 'late' stages of Hh maturation (Guerrero and Kornberg, 2014). Unlike previous studies of photoreceptor axon transport (e.g. behavioral readouts and retina/growth cone ratios in the brains of live or fixed animals), we have demonstrated the ability to directly visualize transport in real time in intact and live animals, and the option to convert 3D time-lapse data into a 2D kymograph (Liu et al., 2000;Stowers et al., 2002;Wyan-Ching Mimi Lee et al., 2004). Kymographs, which are routinely applied to mitochondrial axon transport, can be a great source of information as they reveal subtle nuances in the transport behavior of cargo (Baqri et al., 2009). Precise analytical details such as direction, volume and velocity that are difficult to derive from movies can be conveniently quantified from kymographs. The direction of traces reflects directionality of transport. For example, Hh-positive puncta that move in the anterograde direction trace a forward trajectory and vice-versa for retrograde transport. Therefore, even the most abrupt transient reversals in transport are easily detected on kymographs. Such fine analysis of transport is inconceivable using techniques such as immunohistochemistry and western blotting that have been traditionally applied to Hedgehog transport studies.
Our method also offers a unique opportunity to characterize how disruption of cell biological components (endoplasmic reticulum, Golgi, endocytosis) or mutation of specific residues in Hh (e.g. mutations in holoprosencephaly) affect transport. Since developing Dmel photoreceptors resemble neurons (Liu et al., 2000;Wyan-Ching Mimi Lee et al., 2004), it is also possible to use this technique to complement a screen for regulators of Hh axonal transport to see which larvae have enhanced transport of Hh to the brain or the eye (Stowers et al., 2002). More broadly, our method can be applied to characterize the movement and directionality of other biologically important cargoes (e.g. presynaptic markers, labeled Golgi or mitochondria) implicated in development or in various disease models.
We also characterized Hh transport in a developed model of photoreceptor (PR) neurodegeneration. Gene knockdown of fatp was chosen because of Fatp's role in the transport of palmitic acid, a fatty acid that is covalently attached (in addition to cholesterol) to create mature Hh ligand. In our characterization (Fig. 5) we observed dramatically reduced anterograde transport of Hh, accompanied by a ∼2.1× decrease in anterograde velocity. Retrograde transport also increased ∼3.1× relative to wild type with a concomitant increase in retrograde flux (∼3.9×). These results suggest that less 'mature', lipidated Hh is being sent down axons and that a large proportion of Hh that reaches growth cones is being sent back, since fatp RNAi larvae exhibited slower-moving, lower flux, anterograde Hh. Notably, we do not believe that defective secretion is causing the aberrant transport of Hh since another post-translationally lipidated, membrane-targeted protein, Rhodopsin-1, appears to properly localize and oligomerize in fatp mutant Drosophila (Dourlen et al., 2012;Kock et al., 2009;O'Brien and Zatz, 1984). Thus, we have identified a novel mode of PR toxicity which is perceptible long before PR death occurs. This finding would not have been apparent if only static data had been used.
Finally, it should be noted that a number of Dmel neurodegenerative disease models exist in tissues other than the eye (e.g. motor neurons, mushroom bodies), but the dynamic characterization of axon transport in these models involves ex vivo culturing of dissected brains (Medioni et al., 2015;Millecamps and Julien, 2013;Rabinovich et al., 2015;Sinadinos et al., 2009;Vos et al., 2008;Williamson and Hiesinger, 2010). While these techniques are excellent and will likely yield many new insights, they are labor-intensive and thus, less amenable to use in genetic screens. Our method, by contrast, requires minimal preparation and allows the complement of adult eye phenotypes and live imaging of intact larvae.

CONCLUSIONS
We demonstrate the ability to comprehensively quantify the movement and directionality of a developmentally important cargo in neurons during eye-brain development in a living organism in healthy and mutant conditions. This enables us to ask questions regarding the directionality, velocity, and flux of particles and investigate their behavior relative to axonal position. We posit that this method could be adapted to neurodegenerative disease models (e.g. in the Dmel eye) and allow future research to characterize the role of important cargoes (e.g. secreted ligands, organelles) and candidate drugs in the fields of cell polarity, axonal transport, and eye-brain development in any practical, transparent, and genetically tractable model system.

Time-lapse imaging
Crawling third-instar larvae expressing fluorescently labeled HhC (ss:: GFP::HhC) (Chu et al., 2006) were anaesthetized in halocarbon oil with 15% chloroform, sufficient to inhibit muscular contraction. Larvae were mounted between a slide and coverslip, and stretched under pressure to bring nerves closer to the cuticle ( Fig. 2A,B). Images were acquired with a 63× PlanApo oil objective, NA 1.4, on a Zeiss 7 Live upright scanning confocal with a CCD detector. Images were captured at 0.5 Hz, for a total elapsed time of ∼3-8 min at ∼25°C. All methodology and statistics (including choice of sample size, exclusion criteria, double blind test, randomization, and choice of statistical test) were performed as previously described and according to standard procedures for this type of Drosophila live-imaging data (Baqri et al., 2009;Miller and Sheetz, 2006). A detailed description of this preparation is included below.

Analysis
To generate kymographs, time-lapse frames were aligned horizontally with cell body on the left and axon terminal on the right, re-sliced and z-projected. Kymographs in Fig. 3B,C were color inverted in Adobe Photoshop (Adobe, San Jose, CA, USA) to facilitate visibility of transport events. For quantification of anterograde and retrograde transport, total numbers of transport events in each direction were counted at two different positions in each axon and averaged. For velocity of transport, lines were hand-traced over the path of moving punctae. The slope was calculated for each case.
For Hh distribution analysis in Fig. S1, specimens were viewed on a Zeiss LSM700 Inverted confocal microscope 40× PlanNeofluar oil objective, NA 1.3, with constant acquisition settings when comparing specimens within a given experiment. Quantification of growth cone, optic stalk and eye disc fluorescence was performed with Image J as described previously (Chu et al., 2006). Growth cone/optic stalk, eye disc/optic stalk, and growth cone/ eye disc average fluorescence ratios were computed as in (Chu et al., 2006) as well. A detailed description of this preparation is also included at the end of the methods section.
Step-by-step method for data acquisition and analysis Reagents