Biological systems are organized into well-ordered structures and can evolve new patterns when perturbed. To identify principles underlying biological order, we turned to C. elegans for its simple anatomy and powerful genetics. We developed a method to quantify the arrangement of three dendrites in the main sensory nerve bundle, and found that they exhibit a stereotyped arrangement throughout larval growth. Dendrite order does not require prominent features including sensory cilia and glial junctions. In contrast, loss of the cell adhesion molecule (CAM) CDH-4/Fat-like cadherin causes dendrites to be ordered randomly, despite remaining bundled. Loss of the CAMs PTP-3/LAR or SAX-7/L1CAM causes dendrites to adopt an altered order, which becomes increasingly random as animals grow. Misexpression of SAX-7 leads to subtle but reproducible changes in dendrite order. Our results suggest that combinations of CAMs allow dendrites to self-organize into a stereotyped arrangement and can produce altered patterns when perturbed.https://doi.org/10.7554/eLife.35825.001
Biological systems reflect a remarkable hierarchy of structural organization: proteins assemble into nanometer-scale machines, that in turn build cells, that are then organized and patterned to form structures as complex as the human brain. How do these pieces come together in predictable ways? From first principles, one can imagine a top-down deterministic approach in which the shape, position, and connections of each biological part are specified, similar to an electrical wiring diagram or an architectural blueprint. While this approach is reasonable for simple systems, it breaks down quickly in the face of the complexity we encounter in biology. For example, the human brain contains over 100 billion neurons and glial cells (Herculano-Houzel, 2009), all of which are precisely connected to form a functioning organ. It is hard to imagine a top-down program that can reliably produce such well-organized brains across a population of individuals.
A powerful alternative strategy for organization is a bottom-up, rules-based approach which has been widely considered in studies of morphogenesis and pattern formation. This type of approach has been called ‘emergent’ organization, because well-ordered patterns emerge from a series of local interactions rather than being specified by a blueprint. For example, Alan Turing’s reaction-diffusion mechanism, which consists of an activator molecule that can make more of itself, an inhibitor molecule that inhibits production of the activator, and a mechanism for diffusing these two molecules (Turing, 1952), has been suggested to create the diverse array of biological patterns found in fish stripes (Kondo and Asal, 1995), in seashells (Meinhardt, 1995), and in fur and skin coloration such as those of giraffes and lizards (Manukyan et al., 2017; Walter et al., 1998). A second example of a rule that can generate biological patterns is differential adhesion, in which biological parts sort themselves based on adhesivity to create elaborate structures and patterns. This concept was first demonstrated in 1955 by Townes and Holtfreter, who dissociated different germ layers of embryonic amphibian cells, mixed them together, and found that these cells re-aggregated and sorted out into layers according to their cell type (Townes and Holtfreter, 1955). Differential adhesion has also been shown to organize nerve bundles in the Drosophila visual system (Schwabe et al., 2014). From an evolution standpoint, these emergent, bottom-up strategies may be more advantageous than following a global blueprint, as changes to the starting conditions do not necessarily lead to a disorganized jumble, but instead can give rise to novel well-ordered patterns. For example, we recently showed how altering cell number in the C. elegans nervous system gives rise to an emergent, well-ordered pattern of dendrite arbors (Yip and Heiman, 2016).
To identify rules that help to organize cells into other kinds of ordered structures, we decided to study a simple, stereotyped structure in the C. elegans nervous system called the amphid sense organ. The amphid contains 12 sensory neurons and two glial cells. Each neuron extends a single unbranched dendrite that terminates at the nose tip in a sensory cilium that senses environmental stimuli (Ward et al., 1975). Together, these dendrites and their associated glial processes fasciculate to form one of four bilaterally symmetric nerve bundles that constitute the sensory structures of the head. These structures were the first portion of the C. elegans nervous system to be reconstructed by electron microscopy (Ward et al., 1975). Based on analysis of four animals, Ward and colleagues noted that amphid dendrites appeared to be ordered consistently within the bundle relative to one another, commenting that ‘Although individual worms are not precise replicas of each other down to the finest details, they are remarkably exact copies’ (Ward et al., 1975). These observations are also consistent with recent EM (Doroquez et al., 2014). Together, the studies show that amphid bundles are remarkably well ordered, yet the problem of how dendrite order arises has not been pursued, largely because of the painstaking methods required to examine it.
In addition to the evidence from these EM studies, we chose to study dendrite order in the amphid for four reasons. First, the question of how axons and dendrites are organized within bundles is relatively unexplored, yet a neuron's choice of neighbors may affect its developmental outgrowth, mature activity, or susceptibility to damage or age-related disease. Second, the amphid bundle is a relatively simple and well-isolated system compared to other bundles of processes in the nervous system, containing only 12 unique and identifiable dendrites with no gap junctions or synapses between them (Doroquez et al., 2014; Ward et al., 1975). Third, the system allows us to easily distinguish defects in fasciculation versus guidance. Amphid dendrites do not grow outwards to the nose from a stationary cell body, but instead form by anchoring to the embryonic nose tip while the cell bodies migrate away together, a phenomenon known as retrograde extension (Heiman and Shaham, 2009). Thus, amphid dendrite fasciculation can be studied independently of outgrowth and guidance. Finally, promoters with single-cell-specificity are readily available for all amphid neurons, allowing us to easily visualize and genetically manipulate any single amphid neuron in live animals.
To study the cellular and molecular mechanisms underlying dendrite order in the amphid, we first sought to develop methods to measure, quantify, and compare dendrite order in populations of animals. Ideally, we would label each of the 12 amphid dendrites with a unique fluorescent marker, examine cross sections of the dendrite bundle, and determine how dendrites are ordered across individuals. However, we lacked the technical ability to generate 12 differently colored cell-specific markers and, even if we could uniquely mark each neuron, quantifying the degree of stereotypy of 12 dendrites is a mathematically complex problem. Thus, we decided to simplify the problem by labeling only three dendrites and examining their relative order as a proxy for overall bundle order. We reasoned that, if the overall bundle were well-ordered, then any three dendrites would be similarly well-ordered. Our method can be broken down into three parts: imaging, quantification, and statistical analysis (Figure 1).
First, we generated a strain to label three amphid neurons in different colors (AWA:YFP, AFD:CFP, ASE:mCherry) by combining two separately integrated transgenes (one with AWA:YFP, the second with AFD:CFP and ASE:mCherry) (Figure 1A). We used this YFP + CFP/mCherry approach because we consistently observed recombination artifacts when CFP and YFP were introduced on the same transgene. Using traditional widefield deconvolution microscopy, we found that dendrites in newly hatched larvae were too close together to be reliably resolved (Figure 2—figure supplement 1A) but dendrites in older larvae and adults were readily discerned despite each dendrite being only about 0.5 μm in diameter, close to our effective resolution limit (Figure 1A and Figure 2—figure supplement 1A) (Doroquez et al., 2014; Ward et al., 1975).
Next, we sought to quantify the order of dendrites within individual animals (Figure 1B). To do this, we used a custom MATLAB script that extracts the x-, y-, and z-coordinates of the three dendrites, calculates a centroid line that runs between the three dendrites, and, for each point along the centroid, finds the point on each of the three dendrites that is closest to the centroid (Figure 1Bi). These three dendrite points define a single cross section (Figure 1Bii). Thus, independent of any curvature or rotation in the head of a given animal, this approach identifies a series of contiguous cross-sections along the length of the dendrite bundle that intersect the three dendrites exactly once. Each cross section contains a triangle consisting of one ASE, AFD, and AWA dendrite point. We calculated the pairwise distances between these points (Figure 1Bii).
We used these three pairwise distances to quantify dendrite order in two ways. First, we used the longest pairwise distance as a proxy for bundle width, allowing us to measure how tightly or loosely the dendrites are bundled. To do this, we simply plotted the value of the longest pairwise distance at every point along the bundle (Figure 1Biii). For example, for the animal shown in Figure 1, the longest pairwise distance at every point along the length of the bundle is less than 2 μm, suggesting that the dendrites in this animal are tightly bundled together (Figure 1Biii). Second, we created a categorical variable that describes relative dendrite order by identifying the dendrite ‘in the middle’ of the others at every point along the bundle. A dendrite is considered to be in the middle if it is opposite the longest side of the triangle, that is, the longest pairwise distance. For example, if the longest pairwise distance is d1, which is the distance between AWA (Figure 1Bii, yellow ‘A’) and ASE (Figure 1Bii, red ‘E’), then the dendrite in the middle at that point along the bundle is AFD (Figure 1Bii, blue ‘D’). Importantly, this method describes the arrangement of dendrites relative to each other within the bundle, and is not affected by rotation or twisting of the bundle or head. Next, to visualize dendrite order along the length of the bundle in a single animal, we simply plotted the color of the middle dendrite at each position along the bundle in a column (Figure 1Biii). For example, for the animal in Figure 1, the dendrite in the middle is AFD (blue) near the nose tip and then switches to ASE (red) closer to the cell bodies. To visualize dendrite order within a population, we generated a ‘population plot’ by arranging these individual columns side-by-side (Figure 1Biii).
Finally, we employed a statistical approach based on permutation tests to quantitatively compare the observed dendrite order in a population to a simulated random order (Figure 1Ci, see ‘Statistics’ in Methods for additional details including our rationale for choosing this approach). We also used this approach to compare dendrite order between wild-type and mutant populations. In this approach, our null hypothesis is that the two populations are drawn from the same distribution and any differences between them merely reflect sampling error. First, the observed and test populations are compared using a test statistic (we used Fisher's exact test, see Methods) to yield a nominal p-value. These populations are then computationally merged and re-sampled to create mock populations, and the test statistic is recalculated. Repeating this permutation process (500 iterations) gives a representative set of p-values for populations that have the same composition as our true samples but, by definition, differ only by sampling error. This approach yields 501 p-values (1 true nominal p-value+500 p-values from resampling) for each point along the length of the dendrite bundle. The percentile rank of the true nominal p-value at each point is plotted using a log-scale color bar (Figure 1Cii). Darker colors indicate a low rank, meaning that the observed dendrite order is significantly different from the test distribution. The rank value is equivalent to a corrected p-value, for example rankings less than 25/501 are equivalent to p<0.05 and rankings less than 5/501 are equivalent to p<0.01.
The pairwise distances, dendrite triangles, and p-value ranks for all positions along each bundle in this study (n = 475 bundles) can be explored with a graphical interface using an ‘interactive bundle browser’ we created (http://heimanlab.com/ibb).
To summarize, we have developed a robust and semi-automated pipeline to detect, quantify, and compare dendrite order, which now allows us to determine dendrite order in a wild-type population and to ask whether it is altered by various manipulations.
To determine wild-type dendrite order, we imaged animals expressing CFP, YFP, and mCherry in AFD, AWA, and ASE respectively. Animals were synchronized at the first larval stage (L1) and then collected and imaged at three time points (24 hr, early larval stage (L2/3); 48 hr, late larval stage (L4); 72 hr, one-day adult). We found that dendrites are tightly fasciculated throughout larval growth (Figure 2A) (average longest pairwise distance ± standard deviation: 24 hr, 0.88 ± 0.09 μm; 48 hr, 1.11 ± 0.22 μm; 72 hr, 1.51 ± 0.25 μm). The approximate doubling in bundle width may reflect an increased diameter of each dendrite and thus the entire bundle, and roughly corresponds to the overall doubling in length of the head during these stages.
Next, we looked at dendrite order across the population (Figure 2B,C). We visualized dendrite order using population plots as described above (Figure 2C) as well as summary plots that represent the fraction of animals with the ASE, AWA, or AFD dendrite in the middle at every point along the bundle (Figure 2B). The summary plot provides a compact way to visualize how well-ordered dendrites are in a population. If all of the lines converge at 0.33 then it indicates that dendrite order is random for that population, whereas if any of the lines approach 1 or 0 then it suggests that dendrites are arranged in a consistent, non-random order.
We found that amphid dendrites are well-ordered, especially in younger animals (Figure 2B,C; left plots). To quantify the degree of order, we compared our observed samples to a simulated random sample using a permutation test (Figure 2B, blue bars) as well as a chi-squared test (Figure 2—figure supplement 1B, see Methods). Interestingly, although dendrite order is stereotyped, it is not uniform along the length of the bundle. Instead, it exhibits a switch point, which itself is stereotyped (Figure 2B,C). Close to the nose tip, AFD (blue) is predominantly in the middle, but following the switch point, ASE (red) takes over as the middle dendrite for the rest of the length of the bundle. AWA (yellow) is rarely in the middle, and only for short segments (Figure 2B,C). This order is maintained throughout larval development, although it becomes increasingly ‘noisy’ with age. We further confirmed these results by showing that a different set of three amphid dendrites (AWA, YFP; AWC, CFP; ASG, mCherry) are also well-ordered and exhibit a similar switch point (Figure 2—figure supplement 1C).
The cause of the switch point is unclear. Similar discontinuities were observed in other bundles by classical EM, and it was noted that ‘abrupt changes in neighbourhood exhibited by some neurons’ might arise from mechanical obstacles in the local environment or from changes in the composition of the bundle (White et al., 1986). Indeed, examination of classical EM sections reveals changes in the overall shape of the bundle along its length, switching from a cylinder to a sheet as it is pressed against the basement membrane of the pharynx closer to the nose (Figure 2—figure supplement 1D,E) (Altun and Hall, 2005). The composition of the bundle also changes along its length. The neuron AUA extends a dendrite in the posterior portion of the bundle that does not reach the nose, while the amphid socket glial cell extends a process that fasciculates with only the anterior portion of the bundle (Ward et al., 1975). We measured the average positions of the AUA dendrite ending and the amphid socket glial cell, and found the position of the switch point does not correlate with the AUA ending, but does correlate with the position where the socket glial cell process joins the bundle (Figure 2—figure supplement 1F).
To summarize, we used our imaging pipeline to show that amphid dendrites are well-ordered in young animals and that this order is maintained – albeit imperfectly – during larval growth, despite an approximately two-fold increase in the length and width of the bundle and despite ongoing mechanical perturbations caused by locomotion and pharyngeal pumping of the animal.
Our data show that amphid dendrites exhibit the most consistent order in the distal region of the dendrite bundle, closest to the nose. This region is especially interesting because it is rich in cell biological structures, including dendritic cilia used to detect signals from the environment and cell-cell junctions between each dendrite and the amphid sheath glial cell (Figure 3A). By contrast, the remainder of the dendrites are comparatively featureless, lacking gap junctions, synapses, or any other obvious cell biological specializations. We therefore wanted to test whether cilia or interactions with the sheath glial cell contribute to amphid dendrite order.
We first looked at amphid cilia as a potential source of dendrite order. Early EM reconstructions showed that the order of amphid cilia is also stereotyped, albeit different from that of amphid dendrites (Ward et al., 1975). To test the hypothesis that amphid cilia are required for dendrite order, we crossed our markers into a mutant lacking the RFX transcription factor DAF-19, which is required to activate the genetic program for ciliogenesis (Figure 3B) (Swoboda et al., 2000). Because daf-19 mutants constitutively enter a non-reproductive developmental state called dauer, we also introduced a mutation in daf-12, which encodes a receptor required for dauer entry (Antebi et al., 2000). daf-19; daf-12 animals lack cilia but do not enter dauer, allowing the strain to be maintained as a homozygous stock (Senti and Swoboda, 2008). We found that lack of cilia had no effect on fasciculation or dendrite order (Figure 3B, Figure 3—figure supplement 1). To statistically quantify this observation, we performed permutation tests comparing daf-19; daf-12 with simulated random distributions (Figure 3B, blue column) or with our observed wild-type population (Figure 3B, red column). Darker blues indicate that daf-19 dendrite order is non-random, while the absence of darker reds indicates that daf-19 dendrite order does not differ from wild type. These tests are consistent with the qualitative impression that dendrite order is unaffected. As cilia are required for normal neuronal activity, this result also implies that dendrite order does not depend on normal patterns of activity.
To further test the role of dendrite endings in establishing dendrite order, we examined mutants lacking the splicing factor MEC-8 (Lundquist et al., 1996). Whereas wild-type dendrite endings enter the amphid sheath glial cell in a stereotyped order and converge into a central channel, in mec-8 mutants the amphid dendrite endings are disorganized and diverge into separate, isolated channels in the sheath glial cell (Figure 3C) (Perkins et al., 1986). Despite this disorganization, we found that overall dendrite fasciculation and order along the length of the amphid bundle remain normal (Figure 3C, Figure 3—figure supplement 1). Altogether, our results show that proper arrangement of amphid cilia is not required for dendrite order.
We next tested whether interactions with the amphid glial cell are required for dendrite order. Amphid dendrites form cell junctions with the sheath glial cell, and the sheath glial cell secretes factors that promote the normal development and function of the dendrite sensory endings (Bacaj et al., 2008). To test the possibility that the sheath glial cell might impose order on the dendrites, we genetically ablated sheath glia using diphtheria toxin expressed under the control of a late embryonic-stage amphid sheath promoter (Figure 3D) (Bacaj et al., 2008). We were not able to ablate the sheath glia in earlier embryonic stages, as that causes amphid dendrites to fail to extend (Singhal and Shaham, 2017). Our strain also carried a fluorescent marker for the amphid sheath glial cell, allowing us to identify and exclude rare animals in which ablation failed. In glia-ablated animals, we found dendrites to remain well-ordered (blue color bar) and mostly unchanged from wild type (red color bar), suggesting that sheath glia are not required to maintain dendrite order after the amphid bundle initially develops (Figure 3D, Figure 3—figure supplement 1). Together, our data suggest that the conspicuous cell biological features of amphid dendrites – their cilia and cell-cell junctions – are not required for dendrite fasciculation or order.
Next, we considered the possibility that cell adhesion molecules (CAMs) along the lengths of the dendrites might be required for fasciculation and dendrite order. Previous studies in C. elegans and other organisms have shown that cell adhesion molecules such as SAX-7/L1CAM and the large extracellular matrix molecule DIG-1 are required for dendrite and axon fasciculation (Burket et al., 2006; Bénard et al., 2006; Sasakura et al., 2005). However, few studies in any system have explored how CAMs affect the ordering of axons or dendrites within a bundle (Lin et al., 1994; Schwabe et al., 2014). We wanted to identify CAMs that affect amphid dendrite fasciculation and order. Since it is technically unfeasible to conduct a large-scale forward genetic screen for dendrite order defects using our imaging pipeline, we decided to focus first on a small number of candidate CAMs.
To this end, we began by generating a list of 20 CAMs that are known to play roles in axon fasciculation, that have highly enriched expression in amphid neurons, or that interact with CAMs known to affect fasciculation in other areas of the C. elegans nervous system or in other organisms (Table 1). To further prioritize this list, we reasoned that disrupting CAMs that are important for amphid dendrite fasciculation might cause large-scale dendrite disorganization at some frequency, and that we would be able to detect even infrequent defects by using a dye-filling method to label the amphid neurons. Briefly, animals are soaked in 2 mg/ml of the lipophilic fluorescent dye DiI for 45 min, which for unknown reasons leads to bright and highly specific labeling of six amphid neurons, thus providing a fast and marker-independent method to visualize the overall structure of the amphid bundle. Altered amphid structure can reflect loss of adhesion between dendrites or pleiotropic defects in morphogenesis. We used dye-filling to screen mutations in the 20 candidate CAMs using late larval stage (L4) animals. As a positive control, we found that dig-1 mutants exhibited readily apparent defasciculation defects in this assay (24% of animals, Table 1).
Using this approach, we identified seven additional CAM mutants that cause defasciculation defects (ptp-3, sax-7, vab-1, sax-3, cdh-4, unc-40 and nrx-1; Table 1). Interestingly, even the most pronounced defects among these are weakly penetrant (<25%) suggesting that a redundant system of CAMs contributes to dendrite fasciculation. Loss of sax-7 was previously observed to cause ‘loosening’ of amphid dendrites in some animals, along with other defects, but this aspect of its phenotype has not been pursued (Sasakura et al., 2005). Because two mutants (vab-1, sax-3) exhibited gross head morphology defects and two mutants (unc-40, nrx-1) exhibited only very rare defasciculation defects (1 in 50 animals), we chose to focus our initial characterization on the remaining three mutants – cdh-4, ptp-3/LAR, and sax-7/L1CAM.
From our candidate screen, we found cdh-4 to be required for dendrite fasciculation. cdh-4 encodes a Fat-like cadherin characterized by the large number of cadherin repeats in its extracellular domain (Figure 4A). In C. elegans, cdh-4 has been implicated in axon fasciculation in the dorsal and ventral nerve cords (Schmitz et al., 2007; Schmitz et al., 2008). We used the rh310 allele, which introduces a premature stop codon in the extracellular domain (Figure 4A)(Schmitz et al., 2008).
Since only a small number of cdh-4 mutant animals had dendrite fasciculation defects (6%, Table 1), we wanted to ask whether dendrites are still ordered in cdh-4 individuals that had normally fasciculated dendrite bundles. To do this, we crossed our markers into cdh-4 animals. Consistent with our dye-filling results, we observed grossly defasciculated amphid dendrites in only 1/21 animals examined (Figure 4—figure supplement 1A). Quantification of dendrite distances identified four additional animals in which a portion of the bundle significantly exceeded the average width of wild-type bundles (z-score >3.5, see Methods) and which were therefore also classified as defasciculated. These defasciculated bundles are shown in pink in Figure 4C and were excluded from further analysis. In the remaining 16/21 animals, we found that dendrite order was lost despite the bundle remaining intact and well-fasciculated (Figure 4B), with the ASE, AFD, and AWA dendrites occupying the middle position with approximately equal frequency along the entire length of the bundle (Figure 4D). Permutation tests confirmed that dendrite order in this population is not significantly different from random (light shading in blue color bar) and is significantly different from wild type (dark shading in red color bar).
We found that dendrite order defects are also present in younger animals (24 hr, L2/L3 stage; Figure 4—figure supplement 1C). Because amphid dendrites grow by retrograde extension, we considered the possibility that these early defects in dendrite order might reflect mispositioning of neuronal cell bodies relative to each other during embryonic development. However, cell body positioning appeared normal in newly hatched larvae (Figure 4—figure supplement 1B).
Our results show that, although cdh-4 mutants exhibit only partially penetrant defects in fasciculation, they show nearly complete loss of dendrite order within the bundle, suggesting CDH-4 plays an essential role in specifying dendrite order despite having a more redundant role in overall fasciculation.
We next examined how ptp-3 and sax-7 contribute to amphid dendrite fasciculation and order. ptp-3 encodes a receptor-like protein tyrosine phosphatase that is part of the leukocyte antigen related (LAR) family of proteins that play important roles in nervous system development, including axon guidance and fasciculation (Ackley et al., 2005; Clandinin et al., 2001; Dunah et al., 2005; Harrington et al., 2002; Krueger et al., 1996; Maurel-Zaffran et al., 2001; Wills et al., 1999). sax-7 encodes a homolog of the L1 cell adhesion molecule that is also widely involved in neurodevelopment, including axon guidance and fasciculation, and is disrupted in a human neurological disorder called L1 syndrome (Bénard et al., 2012; Chen and Zhou, 2010; Kim and Emmons, 2017; Sakurai, 2012; Sasakura et al., 2005).
To examine dendrite defects, we crossed our markers to ptp-3 and sax-7 mutants. There are short and long isoforms of ptp-3 (ptp-3a and ptp-3b; Figure 5A), and we made use of ptp-3(mu256), which disrupts both isoforms, as well as ptp-3(ok244), which disrupts only the longer ptp-3a isoform (Figure 5A). There are also short and long isoforms of sax-7 (sax-7S and sax-7L respectively (Bénard et al., 2012; Pocock et al., 2008; Sasakura et al., 2005; Wang et al., 2005), (Figure 5B). We made use of two alleles that disrupt both isoforms – sax-7(ky146) was used for most of our analyses, and we confirmed our results using sax-7(eq1) (Wang et al., 2005; Zallen et al., 1999) – as well as sax-7(nj53) which prevents expression of sax-7L (Sasakura et al., 2005).
As expected from our dye-filling assays, ptp-3(mu256) and sax-7(ky146) exhibit amphid defasciculation defects with low penetrance (pink traces in Figure 5C,F and representative images in Figure 5—figure supplement 1A,D). Fasciculation defects were not observed in ptp-3(ok244), suggesting that the shorter PTP-3B isoform is sufficient for fasciculation (Figure 5—figure supplement 1C). In both ptp-3(mu256) and sax-7(ky146), fasciculation defects become progressively more severe throughout larval development (Figure 5C,F), consistent with previous observations that SAX-7 is required for the maintenance of nervous system architecture (Bénard et al., 2012; Sasakura et al., 2005).
To assess dendrite order, we examined only dendrite bundles that remained tightly fasciculated. In older animals (72 hr, adult stage), ptp-3(mu256) and sax-7 exhibited random dendrite order, similar to what we observed in cdh-4 (Figure 5D,G). These defects are specific to the PTP-3B isoform and to SAX-7, as dendrite order is not affected by loss of PTP-3A alone (ptp-3(ok244), Figure 5—figure supplement 1C) or by disruption of adhesion molecules that interact genetically or physically with SAX-7 or its homologs in other contexts (DMA-1, CLR-1, IGCM-1) (Figure 5—figure supplement 2) (Dong et al., 2013; Islam et al., 2003; Liu et al., 2016; Salzberg et al., 2013). While sax-7 has been reported to exhibit cell body positioning defects in ~10% of newly hatched animals (Sasakura et al., 2005), we found that cell bodies remained ordered normally relative to each other in newly hatched animals for both ptp-3 and sax-7 (Figure 5—figure supplement 1B,E). As a further specificity control, we examined a second allele of sax-7, eq1, and observed similar defects to sax-7(ky146) (Figure 5—figure supplement 1F). Interestingly, sax-7(nj53) exhibited defasciculation defects (10/26 animals) but mild or no defects in dendrite order (Figure 5—figure supplement 1G).
In contrast to adults, ptp-3(mu256) and sax-7 young larvae (24 hr, L2/L3 stage) exhibited a dendrite order that was neither random nor wild type (Figure 5D,G). In ptp-3(mu256) animals at the 24 hr (L2/L3) and 48 hr (L4) time points, the bundle is arranged in a stereotyped order that resembles wild type but has the switch point shifted towards the nose, such that ASE (red) occupies the middle position throughout nearly the entire bundle (Figure 5D). In sax-7 animals at the 24 hr (L2/L3) time point, the middle position is occupied with roughly equal frequency by ASE (red) and AWC (yellow) but almost never by AFD (blue), indicating that dendrite order is not random. As animals reach the 48 hr (L4) time point, the order of dendrites appears increasingly random, and by 72 hr (adult) it is indistinguishable from the random order seen in cdh-4 (Figure 5G). Loss of PTP-3 and SAX-7 together (ptp-3(mu256); sax-7(ky146) double mutant) leads to even greater defasciculation and increased randomness of dendrite order at the 48 hr (L4) time point (Figure 5—figure supplement 1H,I).
Together, our results suggest that PTP-3 and SAX-7 are required, first, to establish the wild-type dendrite order during early development and, second, to maintain dendrite order throughout larval growth. The observation that loss of a single CAM can generate an order that is neither wild type nor random suggests that dendrites might normally ‘choose’ from among several potential neighbors, and that these preferences are altered when a given CAM is removed.
Recent work has shown that cell-specific differences in the expression level of a single CAM helps to determine neighbor preferences among Drosophila retinal axons (Schwabe et al., 2014). Therefore, we examined the expression patterns of ptp-3b and sax-7S by fusing promoter regions upstream of their coding regions (7 kb, ptp-3bpro; 3 kb, sax-7Spro) to a nuclear-localized mCherry (Figure 5E,H). We imaged early larval stage (L1) animals and used the six dye-filling neurons, stained in green, as landmarks to assist in cell identification. We found that ptp-3bpro is expressed at low but detectable levels in many cells, including several amphid neurons, with highest expression in AWB and ASE (Figure 5E). We found that sax-7Spro is expressed throughout the nervous system and, within the amphid, is consistently expressed in AWC, ASE, and ASK (Figure 5H). The differential expression we observe using these synthetic reporter constructs is consistent with a model in which these and other CAMs determine dendrite neighbor preferences through differential adhesion, differences in signaling, or both.
To summarize, we found that loss of either sax-7 or ptp-3 causes amphid dendrites to initially take on a non-random arrangement that is distinct from wild type. In both sax-7 and ptp-3 mutants, this order becomes increasingly random over time, and correlates with increased defasciculation. We also found that sax-7 and ptp-3 were expressed in a subset of amphid neurons, suggesting that these adhesion molecules are not uniformly expressed. Taken together, our data suggest that sax-7 and ptp-3 are required for ordering dendrites within the amphid bundle, such that loss of either CAM causes a weakly-penetrant defasciculation phenotype and a highly-penetrant change in dendrite order.
We considered two models to explain how CAMs promote dendrite order. In the first model, dendrite order is determined during development through a CAM-independent mechanism, for example cell lineage. CAMs would then act like concrete poured over this pre-existing structure to secure it in place such that, without appropriate CAMs, dendrite order would be labile and deteriorate over time. In the second model, dendrite order is determined by differential expression of the CAMs themselves acting as adhesion molecules, signaling receptors, or both. In this model, experimentally misexpressing a single CAM might alter dendrite order (see concept model, Figure 6A).
To test this idea, we misexpressed SAX-7 using promoters that are expressed in one, several, or many additional amphid neurons (odr-10pro, ptp-3pro, osm-6pro respectively). We used SAX-7 because, in our hands, PTP-3 appears to undergo post-transcriptional regulation that makes it difficult to manipulate. Misexpression of SAX-7 did not change the dendrite order drastically, but misexpression in more cells and at higher levels led to subtle yet reproducible changes in the wild-type order, shifting the switch point more posteriorly (Figure 6D). These results do not exclude that additional CAM-independent mechanisms might contribute to dendrite order, but suggest that SAX-7 may play an instructive role in specifying the arrangement of amphid dendrites, possibly by contributing to patterns of differential adhesion.
In this study, we took advantage of the simple nervous system of C. elegans to quantitatively assess the arrangement of individual dendrites within a nerve bundle. We found that amphid dendrites are ordered within the bundle, and this order is maintained over time. Dendrite order does not seem to depend on the conspicuous cell biological structures at the dendrite endings (sensory cilia and dendrite-glia junctions) but is instead imposed by multiple CAMs expressed by the neurons. Loss of the Fat-like cadherin CDH-4 leads to randomized arrangement of dendrites, whereas loss of PTP-3/LAR or SAX-7/L1CAM, or misexpression of SAX-7, can cause dendrites to take on an altered non-random order. Taken together, our results suggest that the arrangement of dendrites within the bundle is determined by multiple CAMs, and altering their relative expression can lead to novel arrangements.
Our results provide insight into the organization of nerve bundles, a defining feature of the nervous system that has previously received relatively little attention due to the technical challenges of studying it. Bundles are known to sometimes correspond to functional groupings of axons, for example each of the cranial nerves is composed of axons with shared functions such as smell, vision, or oculomotor control. Similarly, in the periphery, ‘fast’ and ‘slow’ motor axons are segregated into distinct bundles during development (Milner et al., 1998). In the mammalian cortex, dendrites arising from neurons in different cortical layers come together to form microcolumns (Fleischhauer et al., 1972; Fleischhauer, 1974; Peters and Walsh, 1972). Axons projecting from the retina, as well as the olfactory bulb, have provided evidence that neurites can exhibit stereotyped arrangements within their respective bundles, for example reflecting their topographic or dorsal-ventral origins, or the chronological order of their outgrowth (Chan and Chung, 1999; Walsh and Guillery, 1984). Recent studies have shown that these arrangements arise during development and can affect axon targeting (Cioni et al., 2018; Imai et al., 2009; Zhou et al., 2013). However, the question of how nerve bundles are organized has been difficult to study due to two major technical problems. First, fasciculation defects often lead to axon guidance defects, and vice versa, making it challenging to distinguish one from the other except in rare cases (Luxey et al., 2013). Our approach circumvents this problem by taking advantage of the fact that amphid dendrites grow collectively via retrograde extension (Heiman and Shaham, 2009). Second, most studies have only been able to assess large populations of axons, making it impossible to ask how single neurons or neuron classes are arranged within a bundle. We and others have overcome the latter limitation by turning to invertebrate model organisms with simpler and more stereotyped anatomy that allows the trajectories of defined neurons to be examined.
Two previous sets of studies using invertebrate models are especially helpful to compare with our results. In classical studies using grasshoppers, Goodman and colleagues showed that a certain growing (G) axon selectively adheres to a defined partner (P) axon in a specific bundle. If the P axon is ablated, the G axon does not adhere to other axons in this bundle or in any neighboring bundles, suggesting a remarkable degree of selective adhesion between G and P (Bastiani et al., 1984; Lin et al., 1994; Raper et al., 1984). In a more recent study, Schwabe et al. studied the cartridges of the Drosophila visual system, which exhibit an invariant organization in which neurites from lamina cells (L) are surrounded by photoreceptor (R) axons. They found that this organization reflects differential expression levels of N-cadherin, with the higher-expressing L neurites forming an adhesive core surrounded by the lower-expressing R axons (Schwabe et al., 2014).
These studies led to two major principles, which are further supported by the work described here: first, that neurites can exhibit stereotyped neighbor relationships within a bundle and, second, that these relationships can arise from differential adhesion. Importantly, our study introduces a third principle – that multiple adhesion molecules can act in concert to specify distinct relationships among many neurites, in a pattern that is more complex than the 1:1 pairing of G and P axons or the inside:outside relationship of L and R neurites.
Why might nerve bundles be ordered? One trivial explanation would be that bundle order reflects the developmental order of neurite outgrowth, with the oldest neurite – typically the pioneer axon or dendrite – in the center, and the youngest neurites on the periphery. However, this is not the case in C. elegans, the only organism in which the birth order and neurite position of every neuron is known. As an example, the VD motor neurons are born post-embryonically and extend axons along the ventral nerve cord after the nerve cord is established, yet their axons insinuate deeply into the ventral nerve cord to adhere closely to the DD motor neurons, a functionally similar set of motoneurons that are born much earlier and lie at the center of the ventral nerve cord (White et al., 1976; White et al., 1986). This and similar examples suggest that neurites ‘choose’ their neighbors within a bundle, possibly with functional consequences.
One intriguing hypothesis is that one's immediate neighbors within a bundle influence neuronal activity. For example, the dendrites of gonadotropin-releasing hormone neurons are intertwined in bundles, and this arrangement has been proposed to help synchronize hormone release (Campbell et al., 2009). In a possibly related example, neuropeptide-releasing neurons in Drosophila that serve as circadian pacemakers undergo daily changes in axon fasciculation, driven by cycling expression levels of the NCAM adhesion molecule FAS2 (Fernández et al., 2008; Sivachenko et al., 2013). Interestingly, genetically manipulating FAS2 expression to force constitutive defasciculation leads to changes in circadian behavior, suggesting that fasciculation may indeed affect neuronal activity (Sivachenko et al., 2013).
There are at least three ways in which neighboring neurites could influence each other's activity. First, the presence of chemical synapses or electrical synapses (gap junctions) would couple the activity of immediate neighbors. Notably, no such synapses are present along the amphid dendrites. Second, neuropeptides or other small secreted molecules could locally signal to one's immediate neighbors. Finally, neighboring axons or dendrites could affect each other through passive electrical properties, a concept referred to as ephaptic coupling (Anastassiou et al., 2011; Arvanitaki, 1942). This idea is exciting because it is highly generalizable – anywhere neurites are apposed to one another, such ephaptic effects could be manifested. Consistent with this notion, Ferenczi et al. used an engineering approach to optogenetically manipulate ion flow in hippocampal and cortical axons and observed a ‘bystander’ effect on the membrane current of neighboring axons, independent of synaptic transmission (Ferenczi et al., 2016). Remarkably, in wild-type Drosophila, Carlson and colleagues demonstrated similar ephaptic effects between bundled olfactory dendrites (Su et al., 2012): a short pulse of odor-induced activity in one of these dendrites inhibited tonic firing of its bundled neighbor, independently of chemical or electrical synapses. Importantly, they showed that this ephaptic effect could even alter behavioral responses.
In light of these observations of ephaptic coupling, our finding that dendrites assume stereotyped neighbor relationships could have important implications for how sensory information is processed. Recently developed methods for whole-brain calcium imaging in C. elegans make it an ideal system for asking whether the neighbor relationships we describe here are mirrored by patterns of correlated neuronal activity and, if so, whether such correlations are disrupted by the mutants we found to alter dendrite order (Nguyen et al., 2016; Prevedel et al., 2014; Schrödel et al., 2013).
As Warren Lewis, a pioneer in cell biology, put it in 1922: 'Were the various types of cells to lose their stickiness for one another and for the supporting extracellular white fibers, reticuli, etc., our bodies would at once disintegrate and flow off into the ground in a mixed stream of ectodermal, muscle, mesenchyme, endothelial, liver, pancreatic, and many other types of cells’ (Lewis, 1922). Yet, it is not enough just to stick. Cells care about their neighbors, and arrange themselves into elaborate patterns as they assemble complex organs like the brain. Although it is important to note that CAMs also play major roles in signaling, our work points to the possibility that differential adhesion between cells might contribute to the emergence of biological order. In this model, the arrangement of cells in a tissue is optimized for adhesive strength between them – cells would thermodynamically ‘fold’ into their final arrangement in a manner comparable to protein folding. Importantly, if a single adhesion molecule is absent or misexpressed, the system would reorder itself to maximize adhesion, thus producing a new stable order.
Evolutionarily, this ‘optimize adhesion’ rule provides a general strategy to create a diverse set of well-ordered structures using only small changes in adhesion molecule expression. Interestingly, amphid sensilla in other nematode species such as Acrobeles complexus, Strongyloides stercoralis, and Haemonchus contortus also exhibit well-ordered bundles, albeit with configurations that are different from each other and from C. elegans (Bumbarger et al., 2009). It is interesting to speculate whether this diversity of structures may have arisen partly through altered expression of a small number of CAMs.
Strains were constructed in the N2 background and cultured under standard conditions (Brenner, 1974; Stiernagle, 2006). In addition to the wild-type strain N2, the transgenes and strains used in this study are described in Supplementary files 1–3. Unless otherwise specified, all animals were imaged in the L4 stage, corresponding with ~48 hr after bleach synchronization (see ‘Time point analyses’ below).
To examine dendrite order in cdh-4 mutant animals, we crossed cdh-4(rh310) animals into a strain expressing a three-neuron marker (CHB2646). Because cdh-4 is located on the same chromosome as one of our integrated fluorescent markers (hmnIs23, AWA:YFP), we created a strain (CHB2770) carrying AWA:YFP on an extrachromosomal array (hmnEx1486) and crossed that array into a strain containing hmnIs17 (AFD:CFP, ASE:mCherry) for all cdh-4 analyses. For the permutation test comparing cdh-4 to wild-type dendrite order, we used CHB2646 instead of CHB1963 as the wild-type control.
Animals were bleach-synchronized (20% bleach, 250 mM NaOH in dH2O for 5 min, then hatched overnight in M9 medium), plated onto agar plates containing food, and cultured under standard conditions. We imaged animals at three different time points: second larval stage (L2 stage; 24 hr after plating), fourth larval stage (L4 stage; 48 hr after plating), and 1 day adult (72 hr after plating).
Image stacks were collected on a DeltaVision Core imaging system (Applied Precision) with a UApo 40×/1.35 NA oil-immersion objective (72 hr animals) or a PlanApo 60×/1.42 NA oil-immersion objective (24 hr and 48 hr animals) and a Photometrics CoolSnap HQ2 camera (Roper Scientific). Animals were mounted on an agarose pad with 20–40 mM sodium azide and imaged in yellow (excitation [EX] 513 nm/emission [EM] 559 nm), red (EX 575 nm/EM 632 nm), blue (EX 438 nm/EM 470 nm), and/or green (EX 475 nm/EM 525 nm) channels. To avoid possible complications due to left-right asymmetry, animals were selected for dendrite bundle imaging such that their right-hand side faced the coverslip and only this bundle was imaged.
Deconvolution and analysis of images were performed with Softworx (Applied Precision) and ImageJ (NIH, Bethesda, MD). Maximum-intensity projections were obtained using contiguous optical sections.
Projections were adjusted for brightness, contrast, and were pseudo-colored in Photoshop (Adobe). Merged color images were assembled using the Screen layer mode in Photoshop. Figures were assembled using Photoshop CS5.1 and Illustrator CS5.1.
Image analysis for each animal was done in three parts. First, we imaged three amphid neurons using our three-neuron marker (ASE:mCherry, AWA:YFP, AFD:CFP). We then input this image stack into a custom-made Matlab script that detected the 3D coordinates of each dendrite and generated files containing the distances between each pair of dendrites (AWA-ASE, AWA-AFD, and AFD-ASE) at every point along the dendrite bundle (see Section I: generating 3D coordinates and pairwise distances). Second, we manually inspected and corrected the computer-generated dendrite traces (see Section II: manual inspection and selection of dendrite traces). Finally, we pooled animals belonging to the same population and wrote scripts in Python to generate figures to visualize the data as well as conduct statistical tests to compare populations (see Section III: data visualization, resampling methods, and statistical analysis). Scripts for all the analysis in this paper are available for download at http://github.com/zcandiceyip/dendrite_fasciculation (Yip, 2018; copy archived at https://github.com/elifesciences/dendrite_fasciculation).
The first step in our image-processing pipeline was to obtain 3D coordinates for the ASE, AWA, and AFD dendrites in each animal. To do this, we wrote a script in Matlab that takes an image stack containing three amphid dendrites imaged in three different channels as input and returns several files. One file contains pairwise distances between the three dendrites along the length of the dendrite bundle starting at the dendrite tip. Another set of files contains the digitized dendrite traces superimposed on maximum-intensity projections for manual inspection of the trace accuracy. First, we will describe how the pairwise distances between dendrites are generated.
As mentioned above, our script takes as input an image stack containing three amphid neurons imaged in three different channels. For each amphid neuron, we manually selected the start and end points of the dendrite by clicking on 2D projections of that image generated in Matlab. We consistently chose the starting point to be the ciliated ending at the tip of the dendrite and the ending point to be the cell body. Next, for each dendrite, Dijkstra’s algorithm was used to find the brightest path between the start and end points. Dijkstra’s algorithm finds the shortest path between two nodes, where the objective function is to minimize the distance between the two nodes. In this case, the two nodes are the user-defined start and end points, and distance is defined by the inverse of the intensity of each pixel between the start and end points. Thus, the brightest pixels between the start and end points represent the shortest path. Since the brightest pixels in our image stack between the start and end points correspond to the dendrite itself, Dijkstra’s algorithm yields the 3D coordinates of the dendrite, thus generating a computerized trace of that dendrite. As described below (‘Manual inspection and selection of dendrite traces’), dendrite traces were cropped at a later step to exclude the distal dendrite ending with its complex cilium as well as the proximal portion that crosses the nerve ring (see Figure 1A, arrowheads).
After applying Dijkstra’s algorithm to obtain digital traces of the three amphid dendrites, the next step is to generate a centroid line that runs in the middle of the three dendrites. To do this, we first cropped the start and end points of the three dendrites so that the three dendrite traces were of similar lengths. We did this by averaging the coordinates of the three dendrite tips on each end, selecting the dendrite that gave the shortest distance between the average dendrite tip and the opposite end of that dendrite, and finding dendrite points on the other two dendrites that were closest to the starting point of the shortest dendrite. Those two dendrite points became the new starting points for the two longer dendrites. This process was repeated for the other side of the dendrite bundle. Next, we generated new start and end points of the dendrite bundle by averaging the coordinates of the three cropped dendrites on each end, and found a centroid line by averaging the coordinates of all the points along the three dendrite traces. Finally, we used that centroid line to determine a series of planes that intersect each of the three dendrite traces once by walking along the average trace pixel by pixel, and used Dijkstra’s algorithm to find the shortest path between the centroid line to each of the three traces.
These planes, or cross-sections, yield triangles where each vertex of the triangle is a point along one of the three dendrite traces and the sides of the triangle give the pairwise distances between amphid dendrites in a single cross-section. We then calculated each of these pairwise distances by finding the Euclidean distance between two vertices, and saved them in a comma-delimited (‘.csv’) spreadsheet. In addition, we obtained 2D projections of the computerized traces of each dendrite and superimposed them onto maximum-intensity projections of each image stack to manually confirm that the dendrite trace followed the actual dendrite in the image.
To confirm that the computer-generated dendrite traces and resulting pairwise distances for each animal were accurate, we manually inspected maximum-intensity projections of each image stack projected in XY- and YZ-planes and superimposed the computer-generated dendrite trace onto the projections. If the computer-generated dendrite trace did not follow the dendrite in the maximum-intensity projections, we excluded that animal in further analyses.
In general, the computer-generated dendrite traces followed the dendrites with high fidelity. However, the start and end points of the computer-generated traces were usually inaccurate, as the chosen start and end points were selected a few microns outside the actual start and end points of the dendrite. To correct these inaccuracies and standardize the start and end points, we used the XY-projections to count the number of pixels that were traced inaccurately at the proximal and distal dendrite ends, and deleted that number of rows in the corresponding ‘.csv’ spreadsheet. This technique ensured that we were looking at dendrite tracings that, across all animals, begin at the dendrite tip (excluding the cilium) and end near the nerve ring (see Figure 1A).
To correct for variations in animal size within a population, we measured the length of an anatomical feature of the head of each animal (distance from the anterior bulb of the pharynx to the nose tip) and used that distance to normalize the lengths of the dendrite traces. Because we are making point-by-point comparisons along the dendrite bundle across animals, we further segmented each dendrite bundle into 100 equally-sized bins and took the mean pairwise distance within each bin for further analysis. For example, if an animal had pairwise distance measurements for 200 positions along the dendrite bundle, then we would segment the bundle into 100 bins, with each bin containing pairwise distance measurements from two adjacent positions, and we would assign the mean of those two pairwise distance measurements to that bin. The purpose of binning is to normalize dendrite bundle positions relative to head size; thus, in each individual, bin 50 is at the same position relative to anatomical landmarks but at a different physical distance from the nose (in µm) depending on head size.
For each animal, we plotted the maximum pairwise distance (bundle width) at each position along the dendrite bundle. We defined a dendrite bundle to be defasciculated by comparing its bundle width at each position to the distribution of bundle widths observed at that position in age-matched wild-type animals. If the bundle width exceeded the wild-type average by more than 3.5 standard deviations for at least 10 consecutive positions, then that bundle was considered defasciculated. For reference, in a normally distributed dataset, approximately 1 in 2000 data points exceeds the mean by at least 3.5 standard deviations, and our bundle width dataset consists of 526 dendrite bundles with their widths measured at 100 positions each (52,600 data points). Independently, each bundle was also inspected visually and subjectively classified as fasciculated or defaciculated; these manual calls differed from the automated classification in <10% of cases (46 of 526 bundles). Animals with defasciculated dendrite bundles are overlaid in pink on the bundle width plots and were excluded from further dendrite order analysis.
Each column in the population plot represents the middle dendrite (ASE in red, AWA in yellow, AFD in blue) at each of the 100 bins of the bundle of a single animal. To create summary plots, we calculated and plotted three fractions – the counts of ASE, AFD, and AWA as a fraction of the total for each bin along the dendrite bundle. For populations with ordered dendrite bundles, one of those three fractions should be close to 1 while the other two fractions should be close to 0, whereas populations with highly disordered bundles have all three fractions closer to 0.33.
We used two approaches to test whether the dendrite order for a population is significantly different from random. First, we used a chi-squared test to test for the independence of two populations (genotype vs. random). We chose the chi-squared test because our data is categorical (at each position, the middle dendrite is either AWA, ASE, or AFD) and, for our typical sample sizes (~20), counts of >5 are expected in each category for a random distribution. We calculated the chi-squared values (chi-squared test statistic: ) and associated p-values at each point along the dendrite bundle, where the observed values are the middle dendrite counts of ASE, AWA, and AFD for a given genotype and the expected values are n/3, where n is the number of dendrite bundles analyzed for the given genotype. For example, if n = 21 then we would expect a random distribution to yield 7 counts each of ASE in the middle, AWA in the middle, and AFD in the middle, and these would be compared to our observed values using the formula above.
When comparing a mutant genotype to wild type, we could not use the chi-squared test because the expected value at some positions would be <5 (for example, at the nose tip in 48h L4 animals, both ASE and AWA values are zero, leading to a division by zero error when attempting to use the chi-squared formula; the chi-squared test is also not suitable for comparisons where expected values are <5 in any category). Therefore, we adopted a second approach based on a permutation test using Fisher’s exact test (3 × 2) as the test statistic. In this case, the null hypothesis is that the two populations (for example, mutant genotype vs. wild type) are drawn from the same distribution and differ only by sampling error. We first calculated a nominal p-value at each position along the bundle by comparing the counts of AWA, ASE, and AFD of each population using Fisher’s exact test. Then, for 500 iterations, we merged the two populations of middle dendrite values, randomly split the mixed populations into two equally-sized groups, and calculated p-values comparing these resampled mock populations to each other using Fisher’s exact test. This approach yields 501 p-values (1 true p-value+500 p-values from resampling) for each point along the length of the dendrite bundle. Finally, we determined the percentile rank of the true p-value and plotted that rank using a red log-scale color bar. Darker reds on the color bar indicate that the true p-value is much lower than would be expected if the samples were drawn from the same population; that is, the populations are more different than one would expect from sampling error alone.
We also used this approach to compare each genotype to a random population (Figures 2–6, Figure 2—figure supplement 1, Figure 4—figure supplement 1, Figure 5—figure supplements 1–2 blue bars). In all figures, rather than presenting a mixture of chi-squared and permutation tests, permutation test results are shown for comparison to random (Figures 2–6, Figure 2—figure supplement 1, Figure 4—figure supplement 1, Figure 5—figure supplements 1–2 blue bars) and comparison to wild type (Figures 3–6, Figure 4—figure supplement 1, Figure 5—figure supplements 1–2, red bars). A comparison of chi-squared and permutation tests is shown in Figure 2—figure supplement 1B (green vs. blue bars respectively).
To create this swarmplot (Figure 2—figure supplement 1F), we measured three different lengths in wild-type L4 animals. First, we measured the distance from the switch point to the nose tip. To do this, we defined the switch point to be the point where the middle dendrite changes from one dendrite to another for more than one bin (for definition of a bin, see Section II under ‘Image analysis’ methods section). Second, we measured the distance from the nose tip to the dendrite tip of AUA. Finally, we measured the length of the amphid socket that is fasciculated with the amphid; we used the amphid dendrite AFD (blue) as a proxy for the amphid bundle. We created the figure using the swarmplot function in the Python Seaborn package. To determine the difference between populations, we used a two-sample Kolmogorov-Smirnov test because the measurements for the switch point were not normally distributed.
For each mutant, we used DiO (Sigma, D4292) to dye-fill six amphid neurons and scored 50 L4 animals for amphid defasciculation phenotypes using a fluorescence dissecting microscope (Altun et al., 2002).
To determine the identity of the amphid neurons that express sax-7 or ptp-3 we used DiI (Sigma, 468495) to dye-fill L1 animals expressing a nuclear-localized mCherry (NLS-mCherry-NLS) under the control of sax-7 or ptp-3 promoters (CHB1687 and CHB1840 respectively, see Supp. File 1). We then collected image stacks for these animals in the red and green channels as L1-stage animals as well as positions of all nuclei, visualized using Nomarski optics. Because the identity and positions of the dye-filled amphid neuron cell bodies are known for L1 animals (Sulston et al., 1983), we used the position of the dye-filled amphid cell bodies to infer the identity of the other amphid neurons based on their relative nuclei positioning.
Slice No. 18 and 50SW-Worm Viewer. Accessed October 5, 2017.
DiI and DiO staining in C. elegansWormatlas. Accessed October 5, 2017.
Effects evoked in an axon by the activity of a contiguous oneJournal of Neurophysiology 5:89–108.https://doi.org/10.1152/jn.1922.214.171.124
Three-dimensional reconstruction of the amphid sensilla in the microbial feeding nematode, Acrobeles complexus (Nematoda: rhabditida)The Journal of Comparative Neurology 512:271–281.https://doi.org/10.1002/cne.21882
large-scale screening for targeted knockouts in the Caenorhabditis elegans genomeG3: Genes|Genomes|Genetics 2:1415–1425.
Changes in axon arrangement in the retinofugal [correction of retinofungal] pathway of mouse embryos: confocal microscopy study using single- and double-dye labelThe Journal of Comparative Neurology 406:251–262.https://doi.org/10.1002/(SICI)1096-9861(19990405)406:2<251::AID-CNE8>3.0.CO;2-E
"CRASH"ing with the worm: insights into L1CAM functions and mechanismsDevelopmental Dynamics : An Official Publication of the American Association of Anatomists 239:1490–1501.https://doi.org/10.1002/dvdy.22269
On different patterns of dendritic bundling in the cerebral cortex of the catZeitschrift für Anatomie und Entwicklungsgeschichte 143:115–126.https://doi.org/10.1007/BF00525764
The neurexin superfamily of Caenorhabditis elegansGene Expression Patterns 11:144–150.https://doi.org/10.1016/j.gep.2010.10.008
The human brain in numbers: a linearly scaled-up primate brainFrontiers in Human Neuroscience 3:31.https://doi.org/10.3389/neuro.09.031.2009
Eph:ephrin-B1 forward signaling controls fasciculation of sensory and motor axonsDevelopmental Biology 383:264–274.https://doi.org/10.1016/j.ydbio.2013.09.010
The Algorithmic Beauty of Sea ShellsBerlin: Springer-Verlag.
Caenorhabditis elegans flamingo cadherin fmi-1 regulates GABAergic neuronal developmentJournal of Neuroscience 32:4196–4211.https://doi.org/10.1523/JNEUROSCI.3094-11.2012
Mutant sensory cilia in the nematode Caenorhabditis elegansDevelopmental Biology 117:456–487.https://doi.org/10.1016/0012-1606(86)90314-3
A study of the organization of apical dendrites in the somatic sensory cortex of the ratThe Journal of Comparative Neurology 144:253–268.https://doi.org/10.1002/cne.901440302
Functional dissection of the C. elegans cell adhesion molecule SAX-7, a homologue of human L1Molecular and Cellular Neuroscience 37:56–68.https://doi.org/10.1016/j.mcn.2007.08.014
The role of NrCAM in neural development and disorders--beyond a simple glue in the brainMolecular and Cellular Neuroscience 49:351–363.https://doi.org/10.1016/j.mcn.2011.12.002
Distinct isoforms of the RFX transcription factor DAF-19 regulate ciliogenesis and maintenance of synaptic activityMolecular Biology of the Cell 19:5517–5528.https://doi.org/10.1091/mbc.e08-04-0416
The embryonic cell lineage of the nematode Caenorhabditis elegansDevelopmental Biology 100:64–119.https://doi.org/10.1016/0012-1606(83)90201-4
Directed movements and selective adhesion of embryonic amphibian cellsJournal of Experimental Zoology 128:53–120.https://doi.org/10.1002/jez.1401280105
Fibre order in the pathways from the eye to the brainTrends in Neurosciences 7:208–211.https://doi.org/10.1016/S0166-2236(84)80012-0
Clonal mosaic model for the synthesis of mammalian coat patternsProceedings of Graphics Interface '98.https://doi.org/10.20380/GI1998.11
A role for the C. elegans L1CAM homologue lad-1/sax-7 in maintaining tissue attachmentDevelopmental Biology 284:273–291.https://doi.org/10.1016/j.ydbio.2005.05.020
Electron microscopical reconstruction of the anterior sensory anatomy of the nematode Caenorhabditis elegans.?2uuThe Journal of Comparative Neurology 160:313–337.https://doi.org/10.1002/cne.901600305
The structure of the ventral nerve cord of Caenorhabditis elegansPhilosophical Transactions of the Royal Society B: Biological Sciences 275:327–348.https://doi.org/10.1098/rstb.1976.0086
The structure of the nervous system of the nematode Caenorhabditis elegansPhilosophical Transactions of the Royal Society B: Biological Sciences 314:1–340.https://doi.org/10.1098/rstb.1986.0056
Github, version bf5b296dendrite_fasciculation.
Oliver HobertReviewing Editor; Howard Hughes Medical Institute, Columbia University, United States
Didier YR StainierSenior Editor; Max Planck Institute for Heart and Lung Research, Germany
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Ordered arrangement of dendrites within a C. elegans sensory nerve bundle" for consideration by eLife. Your article has been reviewed by Didier Stainier as the Senior Editor, a Reviewing Editor, and two reviewers. The reviewers have opted to remain anonymous. The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Generally, the overall relevance of the problem addressed is much appreciated by the reviewers, so is the quality of the data and the potential interest of the conclusions. However, there are a number of issues raised by the reviewers (listed below in the individual reviewers’ comments) that both reviewers and the editor feel are very important to address. The issues that particularly stand out are the questions relating (a) to a developmental or a maintenance role of the studied adhesion molecules and (b) the specificity of sax-7 function (is it a generic glue or does it have more specific, instructive roles). These issues shall not simply be addressed by changes in the text but require some very straight-forward experimentation, as suggested by the reviewers.
In this paper Yip and Heiman address how dendrites are organized within nerves. Understanding of how nervous systems assemble and how nerves are organized in bundles is an important problem, which remains poorly understood. Through a simple, but clever method they use three dendritic processes in the nose of C. elegans as a proxy to describe a stereotypic organization of the amphid sensory bundle, which consists of a total of 11 dendrites in C. elegans. The authors show that the stereotypic organization of bundled amphid dendrites is independent of their ciliated endings or associated glia. A candidate genetic screen using cell adhesion molecules (CAMs) allowed them to identify specific CAMs, which are required for dendrite bundle organization and may act redundantly. Gain of function (i.e. overexpression) experiments for one of the CAMs, SAX-7, affect the stereotypic pattern observed in the dendrite bundle. Overall, the data is nicely presented and clearly described, although the text could be somewhat tightened up in places. The authors provide a detailed and thorough statistical analysis to describe and quantify the fasciculation of the amphid dendrite bundle, which may be useful in many other instances and for people in other fields.
A concern is, however, that amphid dendrites develop in a rather unusual way, by retrograde extension; i.e. the ciliated endings of the dendrites attach to the tip of the nose and extend through the retrograde migration of the neuronal cell bodies towards their destination near the posterior bulb of the pharynx. This retrograde cell migration could profoundly influence bundle organization. An analysis of cell body position for the cells and genes being studied at different time points may be required to rule out defects resulting of abnormal migration patterns. Conspicuously absent from the analysis is a time course for the cdh-4 cadherin and an analysis of long isoform sax-7L for which specific alleles exist. Thus, it is not clear, which of the observed bundling defects in adults are actually maintenance defects and which are developmental defects. There are also concerns regarding the overexpression experiments of sax-7, which was overexpressed in all amphid dendrites. These experiments produced rather subtle defects, which are difficult to interpret. Certainly, something was disturbed, but whether that is sufficient to proof the differential adhesion hypothesis or instructiveness remains less clear. I think an opportunity was missed here. More targeted expression of SAX-7 or CDH-4 in specific cells of the amphid bundle could have allowed more concrete predictions and hypotheses to be tested or, at a minimum provide stronger evidence for an instructive function. In fact, more conclusive data for sax-7 misexpression in the same amphid neurons but in regard to the effects of sax-7 on cell body position (and adhesion between cell bodies) was previously described (Sasakura et al., 2005). Finally, the specificity of this SAX-7 misexpression is unclear. Maybe anything expressed (extracellularly or intracellularly) in all amphid dendrites results in subtle defects? Appropriate controls such as for example deletion constructs of SAX-7 or other Ig/FNIII containing proteins would probably have to be tested to make the claim of instructiveness.
In this manuscript, the authors provide some new insights into a largely unexplored area of research, i.e. how do neuronal processes acquire and maintain reproducible positions within a neurite bundle? The authors use C. elegans, more specifically the dendrites of amphid sensory neurons, to identify mechanisms that specify dendrite position. First, they developed a workflow to semiautomatically assess the relative position of three dendrites from light-microscopic images, where dendrites were labelled with different fluorescent proteins. They show that in wildtype the dendrite order is conserved but not completely invariant from animal to animal and is largely maintained during larval growth of an individual animal (thus confirming and expanding the results of previous EM studies). Interestingly the order of the dendrites changes from the tip of the nose (the end of the dendrite) to the cell body with a well-defined switch point close to the nose. The authors then show that cilia and amphid sheath cells (somewhat unexpectedly) are not required to establish the order of the dendrites. Finally, the authors screen a list of candidate cell adhesion molecules for defects in dendrite ordering. They found that mutations in the cadherin cdh-4 randomizes the order of dendrites. Mutations in ptp-3 and sax-7 lead to a dendrite order that is distinct from wild type in early larvae and random later, suggesting a role for these adhesion molecules and receptors in establishing the dendrite order.
The manuscript is well written and easy to follow. The experiments largely support the conclusions.
A major claim stated in the Abstract and the Discussion section is: "Our results suggest that differential expression of CAMs allows dendrites to self-organize …" (quoted from the Abstract). However, the authors provide no evidence that the relevant genes (cdh-4, ptp-3 and sax-7) act as adhesion molecules and not as receptors. Furthermore, only one experiment addresses the "differential expression" aspect (ectopic expression of sax-7). The results of this experiment, essentially no effect on the order of the dendrites (but a shift in the switch point), seems to support the "glue hypothesis" rather than the differential expression hypothesis. I suggest the authors interpret their results more cautiously.
Figure 2: In the 24h plot (Figure 2B), the order of dendrites close to the nose seems random (light blue colors). At 72h, the order close to the nose is non-random (dark color). This is unexpected, since the dendrites attach to the nose during embryogenesis and stay attached. So, dendrite order at the nose should not change during larval development. The authors should discuss this unexpected observation. One possible explanation is that the distances between dendrites at the nose in young animals are at or below the resolution limit of the microscope used. It would be helpful to have more information about the actual distances between dendrites at various stages and positions. In fact, the authors seem to ignore the resolution limit completely. The processing algorithm will always give a position and distance information. However, this is meaningless information, if the distances are smaller than the resolution limit. In general, the authors should exclude such data from the analysis. I consider this a serious issue potentially compromising some of the data.
Resampling for statistical tests (subsection “Statistical analyses”): The authors used permutation tests to compare populations of wildtype and mutant animals. They present the resulting (501) ranked p-values with a color code (e.g. see Figure 3). It is not clear how the ranking of p-values should be interpreted. Darker colors indicate differences that are non-random, but when does this become meaningful? E.g mec-8 (Figure 3C) looks more random and less like wt in the middle section of the dendrite (light blue and dark red in the statistical analysis). Is this meaningful? The authors do not explain how they used the ranked p-values to declare whether mutant and wildtype show "significant" differences.https://doi.org/10.7554/eLife.35825.019
A concern is, however, that amphid dendrites develop in a rather unusual way, by retrograde extension; i.e. the ciliated endings of the dendrites attach to the tip of the nose and extend through the retrograde migration of the neuronal cell bodies towards their destination near the posterior bulb of the pharynx. This retrograde cell migration could profoundly influence bundle organization. An analysis of cell body position for the cells and genes being studied at different time points may be required to rule out defects resulting of abnormal migration patterns.
Altered cell body positioning during retrograde extension would be an interesting mechanism to explain changes in dendrite order, although it would not account for the progressive loss of dendrite order during larval growth. Indeed, sax-7 has been reported to exhibit mispositioned cell bodies in L1 animals, albeit only with 10% penetrance (Sasakura et al., 2005) that would probably not explain the highly penetrant defects in dendrite order that we observe. Cell body positioning for wild-type, cdh-4, sax-7, and ptp-3 animals has been added to Figure 4—figure supplement 1B and Figure 5—figure supplement 1B and E. Consistent with previous reports, we observed occasional cell body mispositioning in newly hatched sax-7 animals, however the relative order of the cell bodies remained wild-type (they did not flip positions like the dendrites).
Conspicuously absent from the analysis is a time course for the cdh-4 cadherin and an analysis of long isoform sax-7L for which specific alleles exist. Thus, it is not clear, which of the observed bundling defects in adults are actually maintenance defects and which are developmental defects.
We have added the 24h time point data for cdh-4 to Figure 4—figure supplement 1C. Dendrites are already randomly ordered at this earlier timepoint. We have also added the sax-7L-specific allele, sax-7(nj53), to Figure 5—figure supplement 1G. This allele is somewhat complicated to interpret, as it prevents SAX-7L expression but may also increase SAX-7S expression (Sasakura et al. 2005), but we agree it makes sense to include it. Dendrite order is mostly unaffected.
There are also concerns regarding the overexpression experiments of sax-7, which was overexpressed in all amphid dendrites. These experiments produced rather subtle defects, which are difficult to interpret. Certainly, something was disturbed, but whether that is sufficient to proof the differential adhesion hypothesis or instructiveness remains less clear. I think an opportunity was missed here. More targeted expression of SAX-7 or CDH-4 in specific cells of the amphid bundle could have allowed more concrete predictions and hypotheses to be tested or, at a minimum provide stronger evidence for an instructive function. In fact, more conclusive data for sax-7 misexpression in the same amphid neurons but in regard to the effects of sax-7 on cell body position (and adhesion between cell bodies) was previously described (Sasakura et al., 2005).
We have added additional sax-7 misexpression experiments, using a promoter strongly expressed in a single neuron (odr-10pro, AWA); a promoter expressed broadly in the amphid at low levels (ptp-3bpro); as well as the promoter we already included that is expressed broadly in the amphid at high levels (osm-6pro). These data are in Figure 6C and suggest a dose-response relationship. However, we appreciate the cautionary note that these experiments do not prove differential adhesion or instructiveness (and do not disprove alternative explanations) and we have revised the text to try to capture this.
Finally, the specificity of this SAX-7 misexpression is unclear. Maybe anything expressed (extracellularly or intracellularly) in all amphid dendrites results in subtle defects? Appropriate controls such as for example deletion constructs of SAX-7 or other Ig/FNIII containing proteins would probably have to be tested to make the claim of instructiveness.
Indeed, we would predict that any construct that modifies cell-surface adhesion could alter dendrite order, including deletion mutants of SAX-7 or other CAMs. We have revised the text to reflect this.
A major claim stated in the Abstract and the Discussion section is: "Our results suggest that differential expression of CAMs allows dendrites to self-organize …" (quoted from the Abstract). However, the authors provide no evidence that the relevant genes (cdh-4, ptp-3 and sax-7) act as adhesion molecules and not as receptors. Furthermore, only one experiment addresses the "differential expression" aspect (ectopic expression of sax-7). The results of this experiment, essentially no effect on the order of the dendrites (but a shift in the switch point), seems to support the "glue hypothesis" rather than the differential expression hypothesis. I suggest the authors interpret their results more cautiously.
These points are all well taken. We have revised the Abstract to highlight "combinations of CAMs" (rather than their differential expression) and "altered order" (rather than a new order) and have made similar changes in the Discussion section. We have added additional data to Figure 6 to support the idea that misexpression of sax-7 produces a subtle but reproducible alteration in dendrite order. We have also reworked the Results section for Figure 6 and parts of the Discussion section to emphasize (i) that our data does not exclude CAM-independent mechanisms that might also contribute to dendrite order (the 'glue' hypothesis) and (ii) that we have not distinguished adhesion vs signaling roles for CAMs.
Figure 2: In the 24h plot (Figure 2B), the order of dendrites close to the nose seems random (light blue colors). At 72h, the order close to the nose is non-random (dark color). This is unexpected, since the dendrites attach to the nose during embryogenesis and stay attached. So, dendrite order at the nose should not change during larval development. The authors should discuss this unexpected observation. One possible explanation is that the distances between dendrites at the nose in young animals are at or below the resolution limit of the microscope used.
This is a good catch; there are several nuggets like this throughout the data that are interesting to speculate about. We hope the online interactive browser will encourage readers to explore them. In this case, the differences are being driven by the 4-6 yellow bins at the nose of the 24h animal. Examining these bundles, we find that the actual pairwise distances are reasonable (~0.4-0.5 µm) so it is unlikely to be a resolution issue, although this was a good guess. We speculate that it may relate to the geometry of the bundle at the nose tip, with younger animals having bundles that are rounder at the nose such that dendrites can engage multiple neighbors at once, and older animals having the bundle squeezed flatter at the nose thus limiting each dendrite to its highest affinity neighbors. But, because it is not a key point and we can only speculate about it, we have chosen not to elaborate on this or other nuggets and to keep our focus on the major conclusions.
It would be helpful to have more information about the actual distances between dendrites at various stages and positions. In fact, the authors seem to ignore the resolution limit completely. The processing algorithm will always give a position and distance information. However, this is meaningless information, if the distances are smaller than the resolution limit. In general, the authors should exclude such data from the analysis. I consider this a serious issue potentially compromising some of the data.
In fact, the resolution limit prevented us from analyzing younger animals (0h, L1 stage) for exactly these reasons. We have added this to Figure 2—figure supplement 1 with a summary of the wild-type average pairwise distances at each developmental stage. Pairwise distances for all genotypes are also now available in the online browser. We had originally considered adding a fourth color (e.g., black) for ambiguous datapoints (using some arbitrary cut-off) but, in practice, sections of the bundle where pairwise distances are essentially equal to each other reveal themselves in the population plots as flickering colors within a column, and they are averaged out in the summary plots and statistics. Indeed, one of the strengths of using our quantitative population-based approach is that it is robust to this kind of noise, and we think its success is demonstrated in the way that the calculated dendrite order is consistent across developmental time points, across a range of genetic perturbations, and agrees well with previous EM analysis.
Resampling for statistical tests (subsection “Statistical analyses”): The authors used permutation tests to compare populations of wildtype and mutant animals. They present the resulting (501) ranked p-values with a color code (e.g. see Figure 3). It is not clear how the ranking of p-values should be interpreted. Darker colors indicate differences that are non-random, but when does this become meaningful? E.g mec-8 (Figure 3C) looks more random and less like wt in the middle section of the dendrite (light blue and dark red in the statistical analysis). Is this meaningful? The authors do not explain how they used the ranked p-values to declare whether mutant and wildtype show "significant" differences.
We have added sentences to the Results section and the Figure 1 legend to clarify that the ranked p-values are equivalent to a corrected p-value (e.g., 25/501 is equivalent to p=0.05). This use of a ranked value statistic is comparable to other commonly used non-parametric ranked tests, for example the Mann-Whitney U-test. There is no cut-off for meaningfulness or significance, but the value of any statistical test is that it allows us to quantitatively compare differences between samples (e.g., to say that the difference between mec-8 and wild type is less than the difference between cdh-4 and wild type). We wanted to show the actual values rather than impose an arbitrary cut-off, especially because our data involve differences that vary in their extent along the length of the bundle (for example, compare the red bars in Figure 6D). However, we appreciate that people perceive color shades differently and some readers may find the color bars challenging to interpret; we hope that the online browser will provide an additional resource to explore and compare the p-value rankings directly (e.g., mec-8 vs wt hovers around p=0.05 for two stretches of about 5 bins each).https://doi.org/10.7554/eLife.35825.020
- Maxwell G Heiman
- Zhiqi Candice Yip
- Maxwell G Heiman
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Tiao Xie and Hunter Elliott of the Harvard Image and Data Analysis Core for writing code for image analysis, and Josh Sanes, Elizabeth Engle, and members of the Heiman lab for comments on the manuscript. We thank Hiroyuki Sasakura and Ikue Mori for the SAX-7S cDNA. We thank the CGC, which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440), and WormBase. We thank David Hall for permission to use Slideable Worm images. This work was supported in part by an NSF Graduate Research Fellowship and Alfred J Ryan Foundation Fellowship (ZCY) and NIH R01GM108754 and the Harvard Milton Fund (MGH). The authors declare that they have no competing interests.
- Didier YR Stainier, Max Planck Institute for Heart and Lung Research, Germany
- Oliver Hobert, Howard Hughes Medical Institute, Columbia University, United States
© 2018, Yip et al.
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.