Conserved neural circuit structure across Drosophila larval development revealed by comparative connectomics

Abstract
Introduction
Results
Discussion
Materials and methods
References
Article and author information
Metrics

Abstract

During postembryonic development, the nervous system must adapt to a growing body. How changes in neuronal structure and connectivity contribute to the maintenance of appropriate circuit function remains unclear. Previously , we measured the cellular neuroanatomy underlying synaptic connectivity in Drosophila (Schneider-Mizell et al., 2016). Here, we examined how neuronal morphology and connectivity change between first instar and third instar larval stages using serial section electron microscopy. We reconstructed nociceptive circuits in a larva of each stage and found consistent topographically arranged connectivity between identified neurons. Five-fold increases in each size, number of terminal dendritic branches, and total number of synaptic inputs were accompanied by cell type-specific connectivity changes that preserved the fraction of total synaptic input associated with each pre-synaptic partner. We propose that precise patterns of structural growth act to conserve the computational function of a circuit, for example determining the location of a dangerous stimulus.

https://doi.org/10.7554/eLife.29089.001

Introduction

As an animal undergoes postembryonic development, its nervous system must continually adapt to a changing body. While developing neural circuits can produce new behaviors, such as the addition of new swimming strategies in zebrafish (Björnfors and El Manira, 2016), in many cases circuit function is conserved as an animal grows. For example, as a Drosophila larva grows from a first instar just out of the egg to a third instar ready to pupariate, its body wall surface area grows by two orders of magnitude (Keshishian et al., 1993). To accommodate this growth, mechanosensory neurons grow their dendrites to maintain receptive fields (Grueber et al., 2002), while larval motor neurons add more synapses at the neuromuscular junction and change firing properties to maintain functional responses in much larger muscles (Keshishian et al., 1993; Guan et al., 1996; Davis and Goodman, 1998; Rasse et al., 2005). Similar functional maintenance has been observed in central circuits as well, from the frequency selectivity of cricket mechanosensory interneurons (Murphey and Chiba, 1990) to motor rhythms in crustacean stomatogastric ganglion (STG) (Bucher et al., 2005).

A neuron’s function arises from the combination of its morphology, synaptic connectivity, and ion channel properties. If morphology and membrane properties co-vary in precise ways, a neuron’s integration properties can be consistent across homologous cells, even between species with very different brain sizes (Cuntz et al., 2013). Homeostatic regulation of functional and structural properties has been proposed as a key principle in neuronal development, allowing consistent output in the presence of both growth and an uncertain or ever-changing environment (Kämper and Murphey, 1994; Bucher et al., 2005; Marder and Goaillard, 2006; Tripodi et al., 2008; Giachello and Baines, 2017).

It remains unclear how neuronal circuits adapt during development by changing their anatomical structure — varying size, adding branches, or producing new synaptic connections — as opposed to adaptation in intrinsic functional properties like ion channel expression and distribution. Studies of circuit variability offer hints, since variability reflects differences in the outcomes of neurons following the same developmental rules. For rhythmic pattern generator circuits, similar temporal dynamics can be produced in many different ways. Simulations of STG have found that numerous different combinations of intrinsic functional parameters and synaptic weights are able to produce extremely similar dynamics (Grashow et al., 2010; O'Leary et al., 2014; Prinz et al., 2004). Correspondingly, the morphological structure of neurons (Otopalik et al., 2017) and their functional connection strengths (Goaillard et al., 2009) have been observed to have high inter-animal variability while still generating similar motor patterns. Observations of inter-animal variability in leach heartbeat networks (Norris et al., 2011; Roffman et al., 2012) suggest that this may be a general principle for rhythm generating circuits.

However, synaptic-resolution electron microscopy (EM) reconstructions from Drosophila sensory systems have found relatively low intra-animal variation in number of synaptic contacts between columnarly repeated neurons in the adult visual system (Takemura et al., 2015) or bilaterally repeated neurons in the mechanosensatory (Ohyama et al., 2015), visual (Larderet et al., 2017), and olfactory (Berck et al., 2016) systems. Comparisons between individuals at this scale have been limited due to incomplete image volumes (Ohyama et al., 2015) or high error rates with early reconstruction methods (Takemura et al., 2013).

Here, we used detailed circuit reconstruction from EM to study the circuitry of identified neurons across postembryonic development in two Drosophila larvae. Despite considerable growth in body size between hatching and pupariation, almost no new functional neurons are added to the larval nervous system (Truman and Bate, 1988) and behavior remains largely unchanged (Almeida-Carvalho et al., 2017). Nonetheless, electrophysiological and light microscopy analysis has shown that central neurons become larger (Zwart et al., 2013) and have more synapses both in total (Zwart et al., 2013) and in specific connections (Couton et al., 2015).

Results

mdIV axon terminals increase in size and number of synapses

We focused on nociception, a somatosensory modality crucial for larvae to avoid wide-ranging sources of damage, such as parasitoid wasp attack (Hwang et al., 2007; Robertson et al., 2013) or intense light (Xiang et al., 2010). Nociceptive stimuli are detected by the three multidendritic class IV sensory neurons (mdIVs) (Hwang et al., 2007) in each hemisegment, with dendrites that tile the body wall (Figure 1A) (Grueber et al., 2002). We began by investigating the structure of the mdIV axon terminals at first and third instar stages. The mdIV terminals in abdominal segment A1 of an early first instar larva EM volume (L1v) were previously reconstructed (Ohyama et al., 2015). We generated a new serial section transmission EM volume of a late third instar larva (L3v) spanning several abdominal segments of the ventral nerve cord (VNC) (Figure 1—figure supplement 1A,B). In the L3v, we manually reconstructed the six mdIV terminals (three per hemisegment) in abdominal segment A3 (Figure 1B,C and Figure 1—figure supplement 1C) using the web-based tool CATMAID (Saalfeld et al., 2009; Schneider-Mizell et al., 2016). Segment A3 was chosen due to its centrality in the L3v and lack of image artifacts or missing sections. In all cases, we reconstructed neurons as skeletons, expressing the 3D topology of neuronal arbors, but not their diameter or volume. The identity of each mdIV terminal was determined based on stereotyped morphological features such as antero-posterior projections, midline crossing, and nerve bundle (Figure 1B,C and Figure 1—figure supplement 2) (Merritt and Whitington, 1995; Grueber et al., 2007; Ohyama et al., 2015).

Figure 1 with 2 supplements see all

Download asset Open asset

Structure of mdIV terminals through postembryonic development.

(A) Cartoon comparison of the dendritic fields of the three nociceptive mdIV sensory neurons from a single hemisegment at first and third instar stages; sagittal view; anterior to left. (B) Dorsal view of EM reconstructions of all mdIV terminals from a single abdominal segment in the L1v (left; segment A1, first instar larva) and L3v (right; segment A3, third instar larva) data. Colors are as in A. The vertical extent of the gray box indicates the width of the neuropil; anterior to left; dashed line indicates midline. (C) Morphology of the terminals of each mdIV subtype, presented as in B. Unbranched primary projections from the nerve are cropped. (D) Dorsal view of a single vdaB terminal from the L1v and L3v, shown with synapses (outputs, red; inputs, cyan). Dashed line indicates midline. (E) Number of synaptic outputs on each mdIV terminal. L1v (solid bars), L3v (empty bars); left/right bar corresponds to left/right neuron. (F) Fold-change in synaptic outputs in the L1v and L3v. For each mdIV subtype, left and right neurons were averaged. (G) A standard polyadic synapse. In this example, taken from the L3v, the single pre-synaptic site (red arrowhead) has four post-synaptic contacts (cyan arrowheads). (H) Normalized histogram of number of post-synaptic contacts per pre-synaptic site on mdIV terminals (No significant difference; p=0.5641, two-sided Kolmogorov-Smirnov test). n.s. not significant; *p $<$ 0.05. **p $<$ 0.01. ***p $<$ 0.001.

https://doi.org/10.7554/eLife.29089.002

The morphology of mdIV axon terminals remained similar across larval stages, with growth in overall size but no change in branching patterns (Figure 1B–D and Figure 1—figure supplement 2). However, the number of synaptic outputs increased by a factor of 4.7, from a mean of 185 synapses per terminal to 872 synapses per terminal (Figure 1E,F). Insect synapses are polyadic, with multiple post-synaptic targets per pre-synaptic site (Figure 1G), thus this increase could arise from either changes in the number of distinct pre-synaptic sites or the number of targets per pre-synaptic site. We found no significant difference between the distribution of number of post-synaptic targets for mdIV pre-synaptic sites in the L1v compared to the L3v (Figure 1H), suggesting the structure of individual polyadic synapses remains unchanged.

Nociceptive interneurons increase in total dendritic cable length and synaptic inputs

The pattern of sensory input onto second-order interneurons is a key component of early sensory processing. To comprehensively identify all second-order mdIV neurons in the L1v, we used all pre- or post-synaptic contacts with mdIV terminals to seed further reconstructions (Figure 2A). We found that there are 13 distinct cell types stereotypically connected to mdIV terminals (Figure 2—figure supplement 1A–C). Five types were local neurons (LNs), with dendrites covering 1–2 segments (Figure 2B); three were regional, with dendrites covering 3–5 segments; four were ascending neurons projecting across the entire VNC; and one was a descending neuron with an axon that projected along the whole VNC (Figure 2—figure supplement 1A). One cell type (A02n) was comprised of two indistinguishable cells per hemisegment, unusual for the larva, making a total of twelve LNs cells per segment. Note that two LNs, A09a and A09c, have been the focus of previous work under the names ’Basin-2’ and ’Basin-4’ (Ohyama et al., 2015; Jovanic et al., 2016). Second-order nociceptive interneurons formed a sparse network (Figure 2—figure supplement 1B,C), without the densely connected local interneurons found in other early processing of other Drosophila sensory modalities like olfaction (Berck et al., 2016; Liu and Wilson, 2013) or mechanosensation (Jovanic et al., 2016; Tuthill and Wilson, 2016).

Figure 2 with 2 supplements see all

Download asset Open asset

Morphology and properties of second-order nociceptive LNs.

(A), Starting from the synapses of the mdIVs in segment A1 of the L1v, we reconstructed all synaptic partners (grays). See Figure 2—figure supplement 1 for details of each cell type. Dorsal view, gray outline indicates CNS boundary; anterior is to left. (B) Examples of the anatomy of all five classes of LNs from the L1v. Posterior view; gray outline indicates neuropile boundary, orange shows mdIV position. (C) Based on the mdIV reconstructions in the L3v (orange), we reconstructed the same populations of all mdIV LNs in segment A3 (grays; 12 LN cells in total). (D) Examples of the anatomy of all five classes of LNs from the L3v, shown as in B. (E) All neurons were split into axonal and dendritic compartments based on well-separated synaptic input and output domains. The example shown is the A02n from D. (F), Total dendritic cable length for all LNs. (G) Number of synaptic inputs onto LN dendrites. (H) Fold-change in dendritic cable length and dendritic synaptic inputs between the L1v and L3v LNs.

https://doi.org/10.7554/eLife.29089.005

To measure how second-order nociceptive interneurons change across larval growth, we reconstructed all twelve third instar LNs in the L3v (Figure 2C). Each LN was morphologically identifiable, despite increases in size, arbor complexity, and synaptic count (Figure 2B,D). For every LN, the spatial segregation of synaptic input and synaptic output made it possible to split neuronal arbors into a separate dendritic domain and axonal domain (Figure 2E). Dendritic cable length, defined as the sum total length of all dendritic neurites, increased by an average factor of $4.69 \pm 0.28$ , consistent with the increase measured from light microscopy in larval motor neurons (Zwart et al., 2013) (Figure 2F,H). The number of synaptic inputs onto LN dendrites increased similarly, by an average factor of $5.28 \pm 0.52$ (Figure 2G,H). Only three LN types had axons fully contained in the L3v (A02n, A08l, and A10a), but our data suggest that axons and dendrites differed in their overall morphological growth (Figure 2—figure supplement 2A–E). In particular, axonal cable length increased by an average factor of only $2.15 \pm 0.33$ , significantly less than the scale-up of dendritic cable (Figure 2—figure supplement 2A,E), and close to the overall 1.7–1.8 times scale-up of neuropile width (L1v: 43 $μ$ m; L3v, 72 $μ$ ) and segment length (L1v, 15 $μ$ m; L3v, 27 $μ$ m). Only those LN types that exhibited dendritic outputs in the L1v also did so in the L3v (Figure 2—figure supplement 2D). For each LN the broad pattern of segregation between synaptic inputs and outputs, which indicates the degree of local dendritic output and axonal pre-synaptic input, was preserved over postembryonic development (Figure 2—figure supplement 2F).

Nociceptive interneurons maintain a topographically-arranged distribution of mdIV synaptic inputs across larval development

Since both mdIV terminals and LN dendrites grow more synapses, we next measured how the synaptic connectivity from mdIVs onto LNs changed across larval development. Every LN in the L3v received synaptic input from mdIVs (Figure 1—figure supplement 1D–F). On average, the total count of synaptic input from mdIVs differed by a factor of $5.77 \pm 1.11$ from the L1v (Figure 3A,C). However, the normalized synaptic input, defined here as the number of synapses in a connection divided by the total number of dendritic input synapses on the post-synaptic neuron, remained strikingly stable, changing by an average factor of $1.09 \pm 0.20$ (Figure 3B,C).

Figure 3 with 1 supplement see all

Download asset Open asset

Connectivity of second-order nociceptive LNs is topographically arranged and consistent across larval development.

(A) Number of synaptic inputs onto LNs from mdIV terminals in the same segment. (B) Normalized dendritic synaptic input from mdIV terminals for each LN. (C) Fold-change in number of synapses and normalized synaptic inputs from mdIVs for each LN type. (D) Heatmap of normalized dendritic input from each mdIV terminal onto each LN for L1v (left) and L3v (right). Note that mdIV terminals are ordered clockwise from ventral left. (E) Normalized dendritic input from mdIVs onto LNs is strongly correlated across animals and developmental time points. Each data point corresponds to average normalized dendritic input from an mdIV type onto an LN type. (Pearson’s r = 0.77, p $<$ 0.001 to be different from zero). (F) Asymmetry between normalized mdIV synaptic input into left and right LNs, measured as coefficient of variation. Asymmetry in the L3v is significantly lower (p=0.006, paired two sided t-test). (G) Cartoon of the larval body wall viewed from posterior. The dendritic receptive field of each mdIV covers approximately 1/6 of the circumference of the animal. (H) Mean body wall orientation of mdIV input into each LN in the L1v (left) and L3v (right), computed as the average of unit vectors pointing at the center of each mdIV dendrite receptive field, weighted by number of synaptic inputs from that neuron. Arrow color corresponds to LN type. n.s. not significant; *p $<$ 0.05. **p $<$ 0.01. ***p $<$ 0.001.

https://doi.org/10.7554/eLife.29089.008

LNs do not receive synaptic input equally from all mdIV subtypes. The normalized synaptic input into each LN from each mdIV axon was highly structured in both L1v and L3v data (Figure 3D). For each mdIV type and LN type, the average normalized synaptic input was significantly correlated between L1v and L3v (Figure 3E). Moreover, the variability between left and right cells of the same type was significantly lower in the L3v than the L1v (Figure 3F). These observations suggest that there is, effectively, a target value for the normalized synaptic input for each connection and this value is achieved more precisely as the nervous system develops postembryonically.

We speculated that the ability to respond according to location of stimuli on the body wall is likely to be an important conserved function of mdIV circuitry. Each segment of the body wall is spanned by six mdIVs whose dendritic fields divide the left and right sides into dorsal, lateral, and ventral thirds (Figure 3G) (Grueber et al., 2002). For each LN, we approximated the mean orientation of its input as the average of unit vectors oriented toward the center of each mdIV dendritic field, weighted by its associated synaptic count (Figure 3H) (see Materials and methods). We found that LN orientations span the body wall, and the orientation of LNs are conserved across development. Further, LN inputs are arranged so that a nociceptive stimulus smaller than a single mdIV’s dendrite, for example a wasp ovipositor (Robertson et al., 2013), is likely to drive different populations of LNs based on its exact location, with the smallest difference being between left and right ventral regions (Figure 3D). Interestingly, only LNs with similar input orientations synaptically connect to one another (Figure 3—figure supplement 1), suggesting that convergent feed-forward motifs are specifically present within sets of neurons likely to be driven at the same time. Conservation of synaptic input through larval development thus preserves the topographical structure of the nociceptive circuit, both in sensory input and interactions between interneurons.

The likelihood of synaptic contact between nearby neurons is stereotyped, cell type-specific, and conserved across larval development

To better understand how synaptic and morphological changes work together to maintain specific patterns of normalized synaptic input, we analyzed the relationship between the spatial location of neuronal arbors and their connectivity. A post-synaptic neuron can only connect to pre-synaptic sites that are nearby in space, or ‘potential synapses’. Numerically strong connections could arise either due to a low probability of connecting to many nearby potential synapses, or to a high probability of connecting to fewer potential synapses. To distinguish these scenarios, we measured ‘filling fraction’ (Stepanyants et al., 2002), defined as the fraction of potential synapses that are actually connected (Figure 4A,B) (see Methods). In both the L1v and L3v, filling fraction ranged from 0.01 to 0.47, indicating that some connections from mdIV types to LNs were realized much more often than others. Filling fraction correlated strongly with the overall count of synapses in a connection, suggesting that numerically strong connections are produced through high connection probability, not only increased potential synapse counts (Figure 4C,D). Moreover, filling fraction was significantly correlated between the L1v and L3v (Figure 4E), suggesting that the local propensity to form stable synapses with a nearby cell type is preserved across development.

Figure 4

Download asset Open asset

Numerically strong connections are associated with stereotypically high filling fraction.

(A) Description of 'filling fraction’ for a connection from Neuron 1 (purple) to Neuron 2 (black). Neurons can only be connected where they are adjacent to one another in space. A pre-synaptic site on Neuron 1 is a potential synapse from Neuron 1 to Neuron 2 if any part of Neuron 2 passes within a given radius (dashed circles). Filling fraction is defined as the number of potential synapses (red and green dashed circles) that are actually connected (green dashed circles only). (B) Dependence of filling fraction on the potential synapse radius for four example connections. For subsequent figures, we chose 2 $μ$ m (filled circles) as a compromise between the typical size of a terminal branch and a shoulder in the filling-fraction versus radius curve. (**C–D**) Mean filling fraction vs. mean number of synapses in the L1v (C) and L3v (D). Each data point represents the average value for connections from mdIV types onto LN types. The high correlation in both (L1v, Pearson $r$ = 0.99, p $<$ 0.001 different from zero; L3v, Pearson $r$ = 0.93, p $<$ 0.001) suggests that increased connection probability, not merely access to differing numbers of pre-synaptic synapses, helps set cell type-specific differences in synaptic counts. (E) Filling fraction of mdIV type to LN type connections in the L1v and L3v are significantly correlated with one another (Pearson $r$ = 0.64, p=0.009).

https://doi.org/10.7554/eLife.29089.010

Post-synaptic connections are not evenly distributed throughout a neuron’s dendrite. Most synaptic input onto a neuron is located on ‘twigs’, spine-like microtubule-free terminal branches hosting a small number of synapses, in contrast to the microtubule-containing ‘backbone’ that spans the soma and all of the main branches of a neuron (Leiss et al., 2009; Schneider-Mizell et al., 2016) (Figure 5A,B). In order to host an increased number of synaptic inputs, a neuron’s twigs would need to change, growing more twigs or hosting more synapses per twig. To measure this, we manually identified all twigs in the twelve LNs (Figure 5—figure supplement 1). We found that the number of twigs increased by an average factor of $2.70 \pm 0.36$ (Figure 5C). The total length of dendritic cable that twigs span increased by a factor of $5.85 \pm 0.31$ , significantly more than dendritic backbone ( $3.31 \pm 0.20$ ) (Figure 5D). Both the fraction of dendritic cable comprised of twigs (Figure 5E,F) and the fraction of dendritic input synapses onto twigs (Figure 5G,H) increased significantly, suggesting that twigs become even more central to dendritic input.

Figure 5 with 1 supplement see all

Download asset Open asset

The structure of terminal dendritic branches across postembryonic development.

(A) Definition of microtubule-containing ‘backbone’ (black) and microtubule-free, spine-like ‘twigs’ (red). (B) Example A02n cell (from the L3v) where all twigs are labeled, posterior view. (C) Number of twigs in each LN in the L1v and L3v. Inset: Fold-change in number of twigs between the L1v and L3v. (D) Fold-change in length of cable comprised of twigs or backbone in the L1v and L3v. Twigs increase more than backbone (two sided t-test). (E) Fraction of dendritic cable comprised of twigs for all LNs. (F) The average fraction of dendritic cable comprised of twigs per cell type was larger in L3v than L1v (two sided, paired t-test). (G) Input synapses that contact twigs as a fraction of all input synapses for all LNs. (H) The fraction of input synapses that are onto twigs increased significantly (two-sided, p=0.003, paired t-test). (I) Cartoon definition of Strahler order. Terminal tips are defined to have Strahler order 1. Where two branches with the same Strahler order converge, the value increments by one. The most core, proximal neurites thus have the highest Strahler order. (J) An example A09a cell from the L1v with branches labeled by Strahler order (Dorsal view). (K) Fraction of dendritic cable for each LN cell by Strahler order. The relative amount of cable with low Strahler order (i.e. distal) is approximately conserved between the L1v and L3v neurons. n.s. not significant; *p $<$ 0.05. **p $<$ 0.01. ***p $<$ 0.001.

https://doi.org/10.7554/eLife.29089.011

An increased number of small twigs host a larger fraction of synaptic input

Measuring twigs requires painstaking visual inspection of EM imagery, so we also looked at a purely topological measure of neuronal arbor structure, Strahler order (Binzegger et al., 2005), that matches intuitive definitions of proximal and distal branches (Figure 5I,J). We found that the fraction of dendritic cable that is last or next-to-last (Strahler order 1 or 2) order is similar not only across development, but also across cell types (Figure 5K). For this observation to be consistent with the relative increase in dendritic twigs for cable, the properties of individual twigs must change so that twigs in the L3v have branches with higher Strahler order than in the L1v. This suggests that in the larva, neurons grow their dendrites by both increasing the number of twigs, while also modestly increasing the length of the backbone neurites from which they sprout.

Twigs remain short and continue to host few synaptic inputs

To get better insight into how twigs changed between the first and third instar, we measured the properties of individual twigs on LNs. Typical dendritic twigs in both the L1v and L3v are small. They were short in both total length and maximum depth from twig root, had few branch points, and few post-synaptic sites (Figure 6A). However, twigs in the L3v were slightly longer than their L1v counterparts and had significantly more branch points (Figure 6A). The median distance between adjacent twigs along neuronal backbone remains similar (L1v: 0.83 $μ m$ ; L3v: 1.03 $μ m$ ), suggesting that the density of twigs on branches remains similar even as neurons grow. In a few cases, we also found that there were quantitative differences between the twig properties of different cell types (e.g. A02n twigs were significantly longer and had greater maximum depth than those of other LNs Figure 6—figure supplement 1), suggesting that individual cell types can deviate from the typical case.

Figure 6 with 1 supplement see all

Download asset Open asset

Twig properties across postembryonic development.

(A) Distribution of morphological and synaptic properties of distinct twigs in the L1v and L3v LNs. Total twig length and maximum twig depth are in $μ$ m, branch points and synapses are integer. Boxes are interquartile intervals, blue dashes are median values, and whiskers correspond to 5/95 percentiles. Wilcoxon rank sum test with Bonferroni correction. (B) Number of distinct twigs in a connection versus number of synapses in the same connection. The blue line indicates a linear fit to connections with five or more synapses (slope shown). The gray region corresponds to the disallowed situation of more twigs connected than synapses. (C) Histogram of number of synapses per twig in each mdIV $\to$ LN connection. (D) The fraction of synapses in the connection from a v’ada to an A09a in the L1v (9 synapses onto twigs) and L3v (56 synapses onto twigs) recovered after simulated random omission of twigs (N = 5000 instances), as a function of omission probability. Thick lines show median value, shaded region the 5/95 percentile value. The dashed horizontal line indicates the 25% of synapses recovered. (E) Maximum error rate permitting the recovery of 25% of synapses with probability $>$ 0.05 for each observed mdIV $\to$ LN connection (i.e. where the horizontal dashed line crosses into the shaded area in D). Each data point is a single mdIV axon’s synapses with a single LN. The vertical line indicates the error rate for twigs found previously for manual annotation of motor neurons (Schneider-Mizell et al., 2016). n.s. not significant; *p $<$ 0.05. **p $<$ 0.01. ***p $<$ 0.001.

https://doi.org/10.7554/eLife.29089.013

The number of distinct twigs involved in a connection increases with the number of synaptic contacts

We next asked how the input from a pre-synaptic sensory neuron is distributed across the twigs on an LN’s dendrite for each mdIV $\to$ LN connection. Consistent with previous work in the first instar motor system (Schneider-Mizell et al., 2016), mdIV $\to$ LN connections with many synaptic contacts were distributed across many twigs in both the L1v and L3v — approximately one twig for every 2.5 synapses in a connection in the L3v (Figure 6B). Within a single mdIV $\to$ LN connection, the vast majority of twigs (L1v: 92.5%; L3v: 81%) hosted only 1 or 2 of the many possible synaptic contacts (Figure 6C).

A practical consequence of numerically strong but anatomically distributed synaptic connectivity is that EM reconstruction becomes robust to random errors. The vast majority of manual errors in previous larval reconstructions was the omission of single dendritic twigs (Schneider-Mizell et al., 2016). To measure how twig omission rate would affect accuracy in measuring mdIV input into LNs, we simulated the effect of removing random twigs from LN reconstructions. For each mdIV $\to$ LN connection, we simulated removing twigs from our anatomical reconstructions with a given omission probability between 0 and 1 (N = 5000 simulations per value), and measured the fraction of synapses that would remain observed in the resulting arbor (Figure 6D,E). For concreteness in comparing connections, we found the maximum error rate for which the probability of detecting fewer than 25% of the observed synapses was $\leq$ 5% (Figure 6E). The anatomical and numerical redundancy of synapses on LN dendrites resulted in connections that would be detectable with at least 25% of the actual number of synapses, even if twigs were missed at the same rate (12%) as observed in previous work with the same reconstruction method (Figure 6E). Numerically strong connections in the L3v were particularly robust, and would still be detectable with a 50% false negative rate (Figure 6E). For large neurons, a strategy of incomplete sampling could thus quickly identify numerically strong synaptic partners at the cost of precise measurement of synaptic count.

Discussion

We have shown how in Drosophila neuronal arbor morphology changes across postembryonic development while circuit connectivity properties remain largely unchanged. Our findings establish a quantitative foundation for the previous observation that numerically strong connectivity in the L1v predicted the presence of functional connectivity in third instar larvae tested experimentally (Ohyama et al., 2015; Zwart et al., 2016; Fushiki et al., 2016; Heckscher et al., 2015; Jovanic et al., 2016). In all neurons measured, the basic anatomical elements of connectivity — polyadic synapses and small post-synaptic twigs — remained similar, while neurons grew five-fold in total synaptic input and cable length. For the highly stereotyped, numerically strong mdIV $\to$ LN connections, the number of synaptic contacts scaled almost identically to the total number of inputs, suggesting the fraction of total inputs per connection is a developmentally conserved value. Interestingly, although cell types ranged considerably in size at any given time point, the fold-increase in total cable length and synapse count was nearly constant across cell types. We note that the sensory connections we focused on here are excitatory (Ohyama et al., 2015). An interesting avenue of future work would be to examine if inhibitory connections follow similar developmental rules.

Compensatory changes in synaptic connectivity and the maintenance of circuit function

The tight control of normalized synaptic input is likely to be in the service of circuit function. Our data suggest that, as neurons grow, there is a consistent compensatory growth in synaptic inputs from sensory neurons. This observation suggests that central neurons adapt structurally to compensate for increasing volume with concurrent increases in excitatory synaptic currents by adding synaptic contacts, as seen at the neuromuscular junction (Rasse et al., 2005). It is possible that such structural changes are also accompanied by functional changes, for example in neurotransmitter receptor or release properties.

Neuronal computations depend on how dendrites integrate synaptic inputs. In visual system interneurons in the adult fly, dendritic geometry and membrane properties work together so that, near the spike initiation zone, the functional weight of a synaptic input does not depend strongly on its location on the dendrite (Cuntz et al., 2013). Similarly, simulations based on adult Drosophila olfactory projection neurons reconstructed from EM found that the functional responses were simply proportional to the number of synapses activated, even after shuffling input locations (Tobin et al., 2017). Taken together, this suggests linear dendritic integration of excitatory input may be common, at least in early sensory processing. In our data, each mdIV input into LNs typically increased by a common factor, irrespective of specific pre-synaptic cell type. Linear integration would thus imply that the relative functional weight of each mdIV type is preserved across development. The higher scaling of synaptic count in the numerically weakest connections (e.g. mdIV $\to$ A09a) could potentially reflect small deviations from linear integration for low numbers of synaptic input.

The same developmental rules that allow neurons to maintain circuit function as the body grows would also be well-suited to handle natural variability, for example from reduced growth due to food restriction (Mirth and Riddiford, 2007). Indeed, it is possible that the use of consistent homeostatic rules for cell type-specific connectivity and integration could allow circuits to remain functionally or computationally similar over large evolutionary changes in neuron size. Such homology has been observed in Calliphora and Drosophila visual system neurons, which were found to differ in scale by a factor of four in each spatial dimension but have retained similar electrotonic structures (Cuntz et al., 2013).

Stringent structural stereotypy observed here stands in contrast to rhythm-generating circuits in other invertebrates, in which large variability can be found in morphological and functional properties (Goaillard et al., 2009; Norris et al., 2011; Roffman et al., 2012; Otopalik et al., 2017). One possibility is that the computation of certain features from sensory input imposes tighter constraints on circuit structure than the production of periodic activity. The ability to combine detailed measurements of structure with cell type-specific genetic reagents (Pfeiffer et al., 2008; Pfeiffer et al., 2010) will allow this hypothesis to be tested across different circuits in the fly and to better elucidate the detailed mechanisms underlying their structural development.

Materials and methods

Sample preparation and electron microscopy

Request a detailed protocol

The L1v is fully described in (Ohyama et al., 2015). In brief, the central nervous system from a 6 hr old [iso] Canton S G1 x [iso] w $^{1118}$ 5905 female larva were dissected and, after chemical fixation, stained en bloc with 1% uranyl acetate, dehydrated, and embedded in Epon resin. Serial 50 nm sections were cut and stained with uranyl acetate and Sato’s lead (Sato, 1968). Sections were imaged at 3.8 $\times$ 3.8 nm using Leginon (Suloway et al., 2005) on an FEI Spirit TEM (Hillsboro). Images were montaged in TrakEM2 (Saalfeld et al., 2010; Cardona et al., 2012) and aligned using elastic registration (Saalfeld et al., 2012).

For the L3v, the central nervous systems from a 96 hr wandering third instar [iso] Canton S G1 x [iso] w $^{1118}$ 5905 larva was dissected in PBS and immediately transferred to 125 $μ$ l of 2% glutaraldehyde in 0.1 M Na cacodylate buffer, pH 7.4 in a 0.5 dram glass vial (Electron Microscopy Sciences, cat. no. 72630–05) on ice. 125 $μ$ l of 2% OsO $_{4}$ in 0.1 M Na-cacodylate buffer, pH 7.4 was then added and briefly mixed immediately before microwave assisted fixation on ice conducted with a Pelco BioWave PRO microwave oven (Ted Pella, Inc.) at 350W, 375W and 400W pulses for 30 s each, separated by 60 s intervals. Samples were rinsed 3 × 30 s at 350W with 0.1 M Na-cacodylate buffer, separated by 60 s intervals, and post-fixed with 1% OsO $_{4}$ in 0.1 M Na-cacodylate buffer at 350W, 375W and 400W pulses for 30 s each, separated by 60 s pauses. After rinsing with distilled water 3 × 30 s at 350W with 60 s pauses between pulses, the samples were stained en bloc with 7.5% uranyl acetate in water overnight at 4 $^{\circ}$ C. Samples were then rinsed 3 $\times$ 5 min with distilled water, dehydrated in an ethanol series followed by propylene oxide, infiltrated and finally embedded in Epon resin. Serial 50 nm sections were cut using a Diatome diamond knife and a Leica UC6 ultramicrotome, and picked up on Pioloform support films with 2 nm C on Synaptek slot grids. Sections were stained with uranyl acetate followed by Sato’s lead (Sato, 1968) prior to imaging. An FEI Spirit TEM operated at 80kV was used to image the serial sections at 2.3 × 2.3 nm pixel resolution using Leginon (Suloway et al., 2005).

L3v image volume registration

Request a detailed protocol

The L3v consisted of $\approx$ 300,000 4k $\times$ 4k image tiles, which were montaged and aligned using linear and nonlinear methods (Saalfeld et al., 2012) in TrakEM2 (Cardona et al., 2012). Filters for brightness and contrast correction were applied before montaging (Default min and max, normalized local contrast, enhance contrast). Images were first montaged in a section with two passes of linear montaging, first targeting only a translation transformation, and in the second pass targeting an affine transformation. This was followed by an elastic, non-linear montaging pass. For alignment between sections, parameter exploration was performed on a scaled down substack (scale factor 10) of 5 sections, targeting extraction of approximately 2000 features, 100 correspondences and an average displacement of 10 pixels. Linear alignment was applied to all sections using an affine transformation model. The ”Test Block Matching Parameters’ tool (http://imagej.net/Test_Block_Matching_Parameters) was used on five adjacent sections to find optimal parameters for the elastic registration pass. Elastic alignment was applied with local smoothness filter approximating an affine local transformation. The resulting aligned image stack was exported to an image tile pyramid with six scale levels for browsing and circuit reconstruction in CATMAID (Saalfeld et al., 2009; Schneider-Mizell et al., 2016). The L3v image stack is available at https://neurodata.io/.

Neuron reconstruction

Request a detailed protocol

For the annotation of mdIV targets in the L1v, we manually reconstructed all neurons pre- and post-synaptic to the previously-described mdIV terminals in segment A1 (Ohyama et al., 2015). Circuit reconstruction in both datasets was performed in CATMAID following annotation and review procedures described previously (Schneider-Mizell et al., 2016). For 1004/1096 post-synaptic connections and 85/85 pre-synaptic connections, we were able to reconstruct an identifiable neuron. This included 173 neurites spanning a total of 30.2 mm in cable length, 13,824 synaptic inputs, and 18,624 synaptic outputs. For each cell type that exhibited more than 3 synapses of input from or output onto mdIV terminals on both left and right sides of the body, we fully reconstructed and comprehensively reviewed a left and right pair of neurons. No unpaired medial neurons were found. For segmentally repeated cell types that exhibited multiple segments of connection, we chose to review examples from the segment with the most synapses from mdIVs, typically segment A1. Reconstructions here were performed by CMSM (30.2%, 104,335/345,917 nodes), IA (28.0%, 96,999/345,917 nodes), Javier Valdes Aleman (8.7%, 29,942/345,917 nodes), Laura Herren (8.1%, 28,165/345,917 nodes), Waleed Osman (7.1%, 24,727/345,917 nodes), and 3% or less each from several other contributors. Comprehensive reviews of arbors and synapses in the L1v were performed by AC and CMSM.

For annotation of the new L3v, we specifically targeted mdIV axons in segment A3 using characteristic anatomical features, particularly entry nerves and the ventromedial location of pre-synaptic boutons. This segment was selected for its centrality in the EM volume and lack of section gaps. Interneurons were identified based on cell body location, neuropil entry point of the primary neurite and characteristic branching structures. To identify target cells from imagery, the principal branches of candidates were reconstructed until they could be conclusively identified from characteristic features. The reconstruction of 6 mdIV terminals and 12 specific LNs spanned 15.3 mm, 10035 synaptic inputs, and 13499 outputs. Reconstructions were performed by IA (50.8%, 91,836/180,753 nodes), SG (21.9%, 39,654/180,753 nodes), CMSM (19.3%, 34,927/180,753 nodes), AC (6.4%, 11,611/180,753 nodes), and Waleed Osman (1.5%, 2,725/180,753 nodes). Comprehensive reviews of arbors and synapses in the L3v were performed by SG, AC and CMSM.

Analysis

Neurons were exported from CATMAID (Saalfeld et al., 2009; Schneider-Mizell et al., 2016) through custom python scripts and imported into python or MATLAB (The Mathworks, Inc.) environments for analysis. Analysis was performed with custom MATLAB scripts with statistics performed using SciPy and R. Morphology and connectivity data were exported from CATMAID and imported into Matlab as a custom neuron data structure to ease analysis. Neuron data structures contained the spatial and topological information for every skeleton node in reconstructions, as well as their polyadic synapses, and annotations such as the location of twig roots and cell bodies. The data structures permitted network-based analysis and visualization using custom scripts and the Brain Connectivity Toolbox (Rubinov and Sporns, 2010). Analysis scripts and files describing neuronal morphology, synapse locations and connectivity can be found at https://github.com/ceesem/Larva_development_structure_2017 (Schneider-Mizell, 2017). A copy is archived at https://github.com/elifesciences-publications/Larva_development_structure_2017.

Neurons were split into axonal and dendritic compartments to maximize spatial segregation along the arbor between synaptic inputs and outputs using previously describes algorithms (Schneider-Mizell et al., 2016). The synaptic segregation index ( $S$ ) was defined as before (Schneider-Mizell et al., 2016):

S = - \frac{1}{S_{0} (N_{a x} + N_{d e n})} \sum_{i = a x, d e n} N_{i} (l o g (p_{i}) + l o g (1 - p_{i}))

where $N_{i}$ is the number of synaptic contacts in compartment $i$ (either axon or dendrite), $p_{i}$ is the fraction of synaptic contacts that are inputs, and $S_{0} = - (l o g (p) + l o g (1 - p))$ for $p$ being the fraction of all synaptic contacts that are inputs (Schneider-Mizell et al., 2016). $S_{0}$ is the maximum possible value of $S$ for a fully unsegregated neuron with the same numbers of synaptic inputs and outputs.

For the receptive field orientation analysis, we defined six unit vectors in a 2D plane $\hat{u_{j}} = \cos (θ_{j}) \hat{x} + \sin (θ_{j}) \hat{y}$ , with the angle $θ_{j}$ corresponding to the approximate center of each of the mdIV terminals ( $θ_{j} = j π / 3$ , with v’ada R corresponding to $j = 0$ and the mdIVs ordered counterclockwise). The mean orientation of interneuron $i$ , ${\vec{r}}_{i}$ , was computed as

\vec{r_{i}} = \frac{1}{\sum_{j = 0}^{5} A_{i j}} \sum_{j = 0}^{5} A_{i j} \hat{u_{j}}

where $A_{i j}$ is the number of synapses from mdIV neuron $j$ to LN $i$ and the sums are over all six mdIVs.

For the filling fraction analysis, we computed potential synapses for a connection from an mdIV terminal onto an LN by computationally removing all terminal branches (Strahler order 1) from the LN dendrites and measuring the number of pre-synaptic sites that were within a distance ( $d = 2 μ m$ unless specified) of the arbor. This approximates a distance that could feasibly be spanned by typical twig growth without overestimating a neuron’s spatial extent.

For the random twig omission errors for a given mdIV $\to$ LN connection, we assumed that each twig could be omitted with an independent probability $p$ . We generated 5000 random instances for each LN and value of $p$ . The synaptic counts were computed by considering only synaptic connections on remaining, non-removed twigs.

References

1. Almeida-Carvalho MJ
2. Berh D
3. Braun A
4. Chen YC
5. Eichler K
6. Eschbach C
7. Fritsch PMJ
8. Gerber B
9. Hoyer N
10. Jiang X
11. Kleber J
12. Klämbt C
13. König C
14. Louis M
15. Michels B
16. Miroschnikow A
17. Mirth C
18. Miura D
19. Niewalda T
20. Otto N
21. Paisios E
22. Pankratz MJ
23. Petersen M
24. Ramsperger N
25. Randel N
26. Risse B
27. Saumweber T
28. Schlegel P
29. Schleyer M
30. Soba P
31. Sprecher SG
32. Tanimura T
33. Thum AS
34. Toshima N
35. Truman JW
36. Yarali A
37. Zlatic M
(2017) The Ol1mpiad: concordance of behavioural faculties of stage 1 and stage 3 Drosophila larvae
The Journal of Experimental Biology 220:2452–2475.

https://doi.org/10.1242/jeb.156646
- PubMed
- Google Scholar
1. Berck ME
2. Khandelwal A
3. Claus L
4. Hernandez-Nunez L
5. Si G
6. Tabone CJ
7. Li F
8. Truman JW
9. Fetter RD
10. Louis M
11. Samuel AD
12. Cardona A
(2016) The wiring diagram of a glomerular olfactory system
eLife 5:e14859.

https://doi.org/10.7554/eLife.14859
- PubMed
- Google Scholar
(2005) Axons in cat visual cortex are topologically self-similar
Cerebral Cortex 15:152–165.

https://doi.org/10.1093/cercor/bhh118
- PubMed
- Google Scholar
1. Björnfors ER
2. El Manira A
(2016) Functional diversity of excitatory commissural interneurons in adult zebrafish
eLife 5:e18579.

https://doi.org/10.7554/eLife.18579
- PubMed
- Google Scholar
(2005) Animal-to-animal variability in motor pattern production in adults and during growth
Journal of Neuroscience 25:1611–1619.

https://doi.org/10.1523/JNEUROSCI.3679-04.2005
- PubMed
- Google Scholar
(2012) TrakEM2 software for neural circuit reconstruction
PLoS ONE 7:e38011.

https://doi.org/10.1371/journal.pone.0038011
- PubMed
- Google Scholar
1. Couton L
2. Mauss AS
3. Yunusov T
4. Diegelmann S
5. Evers JF
6. Landgraf M
(2015) Development of connectivity in a motoneuronal network in Drosophila larvae
Current Biology 25:568–576.

https://doi.org/10.1016/j.cub.2014.12.056
- PubMed
- Google Scholar
1. Cuntz H
2. Forstner F
3. Schnell B
4. Ammer G
5. Raghu SV
6. Borst A
(2013) Preserving neural function under extreme scaling
PLoS ONE 8:e71540.

https://doi.org/10.1371/journal.pone.0071540
- PubMed
- Google Scholar
1. Davis GW
2. Goodman CS
(1998) Synapse-specific control of synaptic efficacy at the terminals of a single neuron
Nature 392:82–86.

https://doi.org/10.1038/32176
- PubMed
- Google Scholar
1. Fushiki A
2. Zwart MF
3. Kohsaka H
4. Fetter RD
5. Cardona A
6. Nose A
(2016) A circuit mechanism for the propagation of waves of muscle contraction in Drosophila
eLife 5:e13253.

https://doi.org/10.7554/eLife.13253
- PubMed
- Google Scholar
1. Giachello CN
2. Baines RA
(2017) Regulation of motoneuron excitability and the setting of homeostatic limits
Current Opinion in Neurobiology 43:1–6.

https://doi.org/10.1016/j.conb.2016.09.014
- PubMed
- Google Scholar
(2009) Functional consequences of animal-to-animal variation in circuit parameters
Nature Neuroscience 12:1424–1430.

https://doi.org/10.1038/nn.2404
- PubMed
- Google Scholar
(2010) Compensation for variable intrinsic neuronal excitability by circuit-synaptic interactions
Journal of Neuroscience 30:9145–9156.

https://doi.org/10.1523/JNEUROSCI.0980-10.2010
- PubMed
- Google Scholar
1. Grueber WB
2. Jan LY
3. Jan YN
(2002)
Tiling of the Drosophila epidermis by multidendritic sensory neurons

Development 129:2867–2878.
- PubMed
- Google Scholar
1. Grueber WB
2. Ye B
3. Yang CH
4. Younger S
5. Borden K
6. Jan LY
7. Jan YN
(2007) Projections of Drosophila multidendritic neurons in the central nervous system: links with peripheral dendrite morphology
Development 134:55–64.

https://doi.org/10.1242/dev.02666
- PubMed
- Google Scholar
1. Guan B
2. Hartmann B
3. Kho YH
4. Gorczyca M
5. Budnik V
(1996) The Drosophila tumor suppressor gene, dlg, is involved in structural plasticity at a glutamatergic synapse
Current Biology 6:695–706.

https://doi.org/10.1016/S0960-9822(09)00451-5
- PubMed
- Google Scholar
1. Heckscher ES
2. Zarin AA
3. Faumont S
4. Clark MQ
5. Manning L
6. Fushiki A
7. Schneider-Mizell CM
8. Fetter RD
9. Truman JW
10. Zwart MF
11. Landgraf M
12. Cardona A
13. Lockery SR
14. Doe CQ
(2015) Even-Skipped(+) Interneurons Are Core Components of a Sensorimotor Circuit that Maintains Left-Right Symmetric Muscle Contraction Amplitude
Neuron 88:314–329.

https://doi.org/10.1016/j.neuron.2015.09.009
- PubMed
- Google Scholar
1. Hwang RY
2. Zhong L
3. Xu Y
4. Johnson T
5. Zhang F
6. Deisseroth K
7. Tracey WD
(2007) Nociceptive neurons protect Drosophila larvae from parasitoid wasps
Current Biology 17:2105–2116.

https://doi.org/10.1016/j.cub.2007.11.029
- PubMed
- Google Scholar
1. Jovanic T
2. Schneider-Mizell CM
3. Shao M
4. Masson JB
5. Denisov G
6. Fetter RD
7. Mensh BD
8. Truman JW
9. Cardona A
10. Zlatic M
(2016) Competitive Disinhibition Mediates Behavioral Choice and Sequences in Drosophila
Cell 167:858–870.

https://doi.org/10.1016/j.cell.2016.09.009
- PubMed
- Google Scholar
1. Keshishian H
2. Chiba A
3. Chang TN
4. Halfon MS
5. Harkins EW
6. Jarecki J
7. Wang L
8. Anderson M
9. Cash S
10. Halpern ME
(1993) Cellular mechanisms governing synaptic development in Drosophila melanogaster
Journal of Neurobiology 24:757–787.

https://doi.org/10.1002/neu.480240606
- PubMed
- Google Scholar
1. Kämper G
2. Murphey RK
(1994) Maturation of an insect nervous system: Constancy in the face of change
Comparative Biochemistry and Physiology Part A: Physiology 109:23–32.

https://doi.org/10.1016/0300-9629(94)90308-5
- Google Scholar
(2017) Organization of the Drosophila larval visual circuit
eLife 6:e28387.

https://doi.org/10.7554/eLife.28387
- Google Scholar
1. Leiss F
2. Koper E
3. Hein I
4. Fouquet W
5. Lindner J
6. Sigrist S
7. Tavosanis G
(2009) Characterization of dendritic spines in the Drosophila central nervous system
Developmental Neurobiology 69:221–234.

https://doi.org/10.1002/dneu.20699
- PubMed
- Google Scholar
1. Liu WW
2. Wilson RI
(2013) Glutamate is an inhibitory neurotransmitter in the Drosophila olfactory system
PNAS 110:10294–10299.

https://doi.org/10.1073/pnas.1220560110
- PubMed
- Google Scholar
1. Marder E
2. Goaillard JM
(2006) Variability, compensation and homeostasis in neuron and network function
Nature Reviews Neuroscience 7:563–574.

https://doi.org/10.1038/nrn1949
- PubMed
- Google Scholar
1. Merritt DJ
2. Whitington PM
(1995)
Central projections of sensory neurons in the Drosophila embryo correlate with sensory modality, soma position, and proneural gene function

The Journal of Neuroscience 15:1755–1767.
- PubMed
- Google Scholar
1. Mirth CK
2. Riddiford LM
(2007) Size assessment and growth control: how adult size is determined in insects
BioEssays 29:344–355.

https://doi.org/10.1002/bies.20552
- PubMed
- Google Scholar
1. Murphey RK
2. Chiba A
(1990) Assembly of the cricket cercal sensory system: genetic and epigenetic control
Journal of Neurobiology 21:120–137.

https://doi.org/10.1002/neu.480210109
- PubMed
- Google Scholar
(2011) Constancy and variability in the output of a central pattern generator
Journal of Neuroscience 31:4663–4674.

https://doi.org/10.1523/JNEUROSCI.5072-10.2011
- PubMed
- Google Scholar
(2015) A multilevel multimodal circuit enhances action selection in Drosophila
Nature 520:633–639.

https://doi.org/10.1038/nature14297
- PubMed
- Google Scholar
(2017) Sloppy morphological tuning in identified neurons of the crustacean stomatogastric ganglion
eLife 6:e22352.

https://doi.org/10.7554/eLife.22352
- PubMed
- Google Scholar
(2014) Cell types, network homeostasis, and pathological compensation from a biologically plausible ion channel expression model
Neuron 82:809–821.

https://doi.org/10.1016/j.neuron.2014.04.002
- PubMed
- Google Scholar
1. Pfeiffer BD
2. Jenett A
3. Hammonds AS
4. Ngo TT
5. Misra S
6. Murphy C
7. Scully A
8. Carlson JW
9. Wan KH
10. Laverty TR
11. Mungall C
12. Svirskas R
13. Kadonaga JT
14. Doe CQ
15. Eisen MB
16. Celniker SE
17. Rubin GM
(2008) Tools for neuroanatomy and neurogenetics in Drosophila
PNAS 105:9715–9720.

https://doi.org/10.1073/pnas.0803697105
- PubMed
- Google Scholar
1. Pfeiffer BD
2. Ngo TT
3. Hibbard KL
4. Murphy C
5. Jenett A
6. Truman JW
7. Rubin GM
(2010) Refinement of tools for targeted gene expression in Drosophila
Genetics 186:735–755.

https://doi.org/10.1534/genetics.110.119917
- PubMed
- Google Scholar
(2004) Similar network activity from disparate circuit parameters
Nature Neuroscience 7:1345–1352.

https://doi.org/10.1038/nn1352
- PubMed
- Google Scholar
1. Rasse TM
2. Fouquet W
3. Schmid A
4. Kittel RJ
5. Mertel S
6. Sigrist CB
7. Schmidt M
8. Guzman A
9. Merino C
10. Qin G
11. Quentin C
12. Madeo FF
13. Heckmann M
14. Sigrist SJ
(2005) Glutamate receptor dynamics organizing synapse formation in vivo
Nature Neuroscience 8:898–905.

https://doi.org/10.1038/nn1484
- PubMed
- Google Scholar
(2013) Larval defense against attack from parasitoid wasps requires nociceptive neurons
PLoS ONE 8:e78704.

https://doi.org/10.1371/journal.pone.0078704
- PubMed
- Google Scholar
(2012) Animal-to-animal variability of connection strength in the leech heartbeat central pattern generator
Journal of Neurophysiology 107:1681–1693.

https://doi.org/10.1152/jn.00903.2011
- PubMed
- Google Scholar
1. Rubinov M
2. Sporns O
(2010) Complex network measures of brain connectivity: uses and interpretations
NeuroImage 52:1059–1069.

https://doi.org/10.1016/j.neuroimage.2009.10.003
- PubMed
- Google Scholar
(2009) CATMAID: collaborative annotation toolkit for massive amounts of image data
Bioinformatics 25:1984–1986.

https://doi.org/10.1093/bioinformatics/btp266
- PubMed
- Google Scholar
(2010) As-rigid-as-possible mosaicking and serial section registration of large ssTEM datasets
Bioinformatics 26:i57–i63.

https://doi.org/10.1093/bioinformatics/btq219
- PubMed
- Google Scholar
(2012) Elastic volume reconstruction from series of ultra-thin microscopy sections
Nature Methods 9:717–720.

https://doi.org/10.1038/nmeth.2072
- PubMed
- Google Scholar
1. Sato T
(1968)
A modified method for lead staining of thin sections

Journal of Electron Microscopy 17:158–159.
- PubMed
- Google Scholar
1. Schneider-Mizell CM
2. Gerhard S
3. Longair M
4. Kazimiers T
5. Li F
6. Zwart MF
7. Champion A
8. Midgley FM
9. Fetter RD
10. Saalfeld S
11. Cardona A
(2016) Quantitative neuroanatomy for connectomics in Drosophila
eLife 5:e12059.

https://doi.org/10.7554/eLife.12059
- PubMed
- Google Scholar
Software
1. Schneider-Mizell CM
(2017) larva_development_structure_2017, version 35c485c
Github.

https://github.com/ceesem/larva_development_structure_2017
(2002) Geometry and structural plasticity of synaptic connectivity
Neuron 34:275–288.

https://doi.org/10.1016/S0896-6273(02)00652-9
- PubMed
- Google Scholar
1. Suloway C
2. Pulokas J
3. Fellmann D
4. Cheng A
5. Guerra F
6. Quispe J
7. Stagg S
8. Potter CS
9. Carragher B
(2005) Automated molecular microscopy: the new Leginon system
Journal of Structural Biology 151:41–60.

https://doi.org/10.1016/j.jsb.2005.03.010
- PubMed
- Google Scholar
1. Takemura SY
2. Bharioke A
3. Lu Z
4. Nern A
5. Vitaladevuni S
6. Rivlin PK
7. Katz WT
8. Olbris DJ
9. Plaza SM
10. Winston P
11. Zhao T
12. Horne JA
13. Fetter RD
14. Takemura S
15. Blazek K
16. Chang LA
17. Ogundeyi O
18. Saunders MA
19. Shapiro V
20. Sigmund C
21. Rubin GM
22. Scheffer LK
23. Meinertzhagen IA
24. Chklovskii DB
(2013) A visual motion detection circuit suggested by Drosophila connectomics
Nature 500:175–181.

https://doi.org/10.1038/nature12450
- PubMed
- Google Scholar
1. Takemura SY
2. Xu CS
3. Lu Z
4. Rivlin PK
5. Parag T
6. Olbris DJ
7. Plaza S
8. Zhao T
9. Katz WT
10. Umayam L
11. Weaver C
12. Hess HF
13. Horne JA
14. Nunez-Iglesias J
15. Aniceto R
16. Chang LA
17. Lauchie S
18. Nasca A
19. Ogundeyi O
20. Sigmund C
21. Takemura S
22. Tran J
23. Langille C
24. Le Lacheur K
25. McLin S
26. Shinomiya A
27. Chklovskii DB
28. Meinertzhagen IA
29. Scheffer LK
(2015) Synaptic circuits and their variations within different columns in the visual system of Drosophila
PNAS 112:13711–13716.

https://doi.org/10.1073/pnas.1509820112
- PubMed
- Google Scholar
(2017) Wiring variations that enable and constrain neural computation in a sensory microcircuit
eLife 6:e24838.

https://doi.org/10.7554/eLife.24838
- PubMed
- Google Scholar
1. Tripodi M
2. Evers JF
3. Mauss A
4. Bate M
5. Landgraf M
(2008) Structural homeostasis: compensatory adjustments of dendritic arbor geometry in response to variations of synaptic input
PLoS Biology 6:e260.

https://doi.org/10.1371/journal.pbio.0060260
- PubMed
- Google Scholar
1. Truman JW
2. Bate M
(1988) Spatial and temporal patterns of neurogenesis in the central nervous system of Drosophila melanogaster
Developmental Biology 125:145–157.

https://doi.org/10.1016/0012-1606(88)90067-X
- PubMed
- Google Scholar
1. Tuthill JC
2. Wilson RI
(2016) Parallel Transformation of Tactile Signals in Central Circuits of Drosophila
Cell 164:1046–1059.

https://doi.org/10.1016/j.cell.2016.01.014
- PubMed
- Google Scholar
1. Xiang Y
2. Yuan Q
3. Vogt N
4. Looger LL
5. Jan LY
6. Jan YN
(2010) Light-avoidance-mediating photoreceptors tile the Drosophila larval body wall
Nature 468:921–926.

https://doi.org/10.1038/nature09576
- PubMed
- Google Scholar
(2013) Dendritic growth gated by a steroid hormone receptor underlies increases in activity in the developing Drosophila locomotor system
PNAS 110:E3878–E3887.

https://doi.org/10.1073/pnas.1311711110
- PubMed
- Google Scholar
1. Zwart MF
2. Pulver SR
3. Truman JW
4. Fushiki A
5. Fetter RD
6. Cardona A
7. Landgraf M
(2016) Selective inhibition mediates the sequential recruitment of motor pools
Neuron 91:615–628.

https://doi.org/10.1016/j.neuron.2016.06.031
- PubMed
- Google Scholar

Article and author information

Author details

Stephan Gerhard

Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, United States

Present address
Friedrich Miescher Institute, Basel, Switzerland

Contribution
Conceptualization, Data curation, Investigation, Visualization, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4454-6171
Ingrid Andrade

Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, United States

Present address
Department of Molecular, Cell and Developmental Biology, University of California Los Angeles, Los Angeles, United States

Contribution
Investigation

Competing interests
No competing interests declared
Richard D Fetter

Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, United States

Present address
Department of Biochemistry and Biophysics, University of California San Fransisco, San Fransisco, United States

Contribution
Resources, Investigation, Writing—review and editing

Competing interests
No competing interests declared
Albert Cardona
1. Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, United States
2. Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge, United Kingdom
Contribution
Conceptualization, Supervision, Investigation, Methodology, Writing—review and editing

For correspondence
cardonaa@janelia.hhmi.org

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4941-6536
Casey M Schneider-Mizell

Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, United States

Contribution
Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Writing—original draft

For correspondence
schneidermizellc@janelia.hhmi.org

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9477-3853

Funding

Howard Hughes Medical Institute

Stephan Gerhard
Ingrid Andrade
Richard D Fetter
Albert Cardona
Casey M Schneider-Mizell

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Stephan Saalfeld for assistance registering the L3v, Maarten Zwart, Matthias Landgraf, Chris Q Doe, and Marta Zlatic for helpful comments, and James Truman for generously sharing light microscopy data to help identify neurons. This work was funded by the Howard Hughes Medical Institute.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.