Abstract
Understanding binaural perception requires detailed analyses of the neural circuitry responsible for the computation of interaural time differences (ITDs). In the avian brainstem, this circuit consists of internal axonal delay lines innervating an array of coincidence detector neurons that encode external ITDs. Nucleus magnocellularis (NM) neurons project to the dorsal dendritic field of the ipsilateral nucleus laminaris (NL) and to the ventral field of the contralateral NL. Contralateral-projecting axons form a delay line system along a band of NL neurons. Binaural acoustic signals in the form of phase-locked action potentials from NM cells arrive at NL and establish a topographic map of sound source location along the azimuth. These pathways are assumed to represent a circuit similar to the Jeffress model of sound localization, establishing a place code along an isofrequency contour of NL. Three-dimensional measurements of axon lengths reveal major discrepancies with the current model; the temporal offset based on conduction length alone makes encoding of physiological ITDs impossible. However, axon diameter and distances between Nodes of Ranvier also influence signal propagation times along an axon. Our measurements of these parameters reveal that diameter and internode distance can compensate for the temporal offset inferred from axon lengths alone. Together with other recent studies, these unexpected results should inspire new thinking on the cellular biology, evolution, and plasticity of the circuitry underlying low-frequency sound localization in both birds and mammals.
Introduction
Binaural processing of acoustic signals is essential for localizing sound and extracting signals in a noisy environment. Low-frequency sounds are localized by interaural differences in the arrival times of sound created by a spatial separation of the ears [interaural time differences (ITDs)]. ITDs are thought to be computed by a neural mechanism similar to that proposed by Jeffress (1948), where external differences of sound arrival times at the ears are represented along internal delay lines.
In the avian auditory brainstem, a modified Jeffress model is represented by nucleus magnocellularis (NM) and nucleus laminaris (NL) (see Fig. 1A). Neurons in NM receive monaural input from the ipsilateral acoustic sensory epithelium of birds via the auditory nerve (AN). NM axons bifurcate and send bilateral projections to neurons in NL (see Fig. 1A). The ipsilateral NM axon provides a simultaneous input to the dorsal dendrites of an isofrequency array of NL neurons (Hyson et al., 1994) and by traveling an extended looped trajectory putatively equalizes conduction times with the contralateral pathway (Jhaveri and Morest, 1982; Young and Rubel, 1983) (see Fig. 1B). Coincidence detection for 0 ITD in the NL has recently been shown in vivo (Köppl and Carr, 2008). The contralateral axons cross the midline and travel on the ventral side along an isofrequency band of NL neurons. Terminal arbors branch off the main axon in serial order, thereby creating a physical delay line (Young and Rubel, 1983).
NL neurons function as an array of coincidence detectors (Carr and Konishi, 1990; Joseph and Hyson, 1993; Hyson et al., 1994; Reyes et al., 1996; Kuba et al., 2002). NL neurons distributed along the delay line respond maximally to particular ITDs, thus establishing a neural map of sound source location in the azimuth. When a stimulus comes from straight ahead (ITD = 0), the specific location on the map receiving equivalent delay from each ear responds maximally. Sound source locations generating ITDs different from 0 are encoded by an array of neurons along internal delay lines that compensate for the different external arrival times (see Fig. 1B).
We analyzed the anatomy of the NM-NL circuit to determine what mechanisms are used to compensate for external delays and its fit for the model described above. Our three-dimensional (3-D) measurements of axon length do not comply with the proposed model in its simple form (Overholt et al., 1992). The contralateral axon is on average >1600 μm longer than its ipsilateral counterpart, rendering encoding of physiological ITDs based on conduction length alone impossible; coincident binaural inputs at NL would only be created by ITDs larger than 180 μs, hence outside the physiological range (Hyson et al., 1994; Hyson, 2005). We also measured axon diameter and distances between Nodes of Ranvier in different segments of the NM axon. These parameters vary across the NM axonal arbor and act to slow transmission of the acoustically evoked signal in the ipsilateral branch of the NM axon compared to the longer contralateral branch. It appears that these determinates of conduction velocity are adjusted to compensate for the axon length offset. Our results provide implications for the ITD detection mechanism not only in birds, but also in mammals, as recent studies have challenged the status quo of its proposed mechanism, i.e., its fit for the Jeffress model (McAlpine et al., 2001; Brand et al., 2002).
Materials and Methods
Animals.
Over 50 white leghorn chickens (Gallus gallus) of varied postnatal ages were studied. The day of hatching was considered to be postnatal day 0 (P0). The University of Washington Institutional Animal Care and Use Committee approved all procedures.
3-D reconstruction of auditory brainstem nuclei.
Chickens of different ages (P0, P1, P14, P15, P27) were deeply anesthetized (urethane 1 mg/kg, ketamine 20 mg/kg) and perfused intracardially with 0.9% saline plus 0.4% heparin, followed by 4% paraformaldehyde in phosphate buffer (PFA) plus 1% methylene blue. Adequate perfusion was achieved when the beak and feet turned dark blue. The entire head was further fixed by immersion in 4% PFA for 24 h at room temperature (RT; 21°C) and then decalcified in RDO (Apex Engineering Products Cooperation) for 3–4 h. Specimens were passed through a series of sucrose gradients over 48 h for cryoprotection and then blocked in the desired plane of sectioning. Then entire heads were put into an embedding medium (Tissue-Tek O.C.T. Compound, Sakura Finetek) under vacuum for 2–3 h to remove air from the middle ear cavities.
Digitized image stacks of aligned and stained serial sections undistorted by histological manipulation were made with a “photomicrotome” (Harris et al., 1990), a freezing microtome with a digital camera (Nikon D200, 10.5 Megapixels, Nikon) attached, focused on the cutting plane. The embedded chicken head was positioned on the stage of the microtome, covered with embedding medium and packed in dry ice. Frozen sections were serially cut at 22 μm. Before each frozen section was cut, an image of the surface of the block was captured and saved online. In this way, the whole brain was imaged in situ along with the skull. The histological sections containing the nucleus angularis (NA), AN, NM, and NL were then stained with cresyl violet to aid in identifying nuclear boundaries and structures during subsequent digital segmentation.
The acquired image stacks were contrast-enhanced in Photoshop (Adobe) and then loaded into Amira software (Mercury Software) for reconstruction of the three-dimensional anatomy. The volume and length of the resulting Amira 3-D models were measured using ImageJ (National Institute of Mental Health, Bethesda, MD).
3-D tracing of individual NM-axons.
Thick slices (1000 μm) containing the auditory brainstem of 12 newly hatched chickens (P0/1) were prepared in ice-cold artificial CSF (ACSF) (Sorensen and Rubel, 2006) using a Vibratome (Vibratome Series 1000, The Vibratome Company) and a sapphire knife (Leica Microsystems). One slice per animal was harvested and immediately incubated in oxygenated ACSF at RT for 45 min. With the caudal side of the slice facing upward, NM was visualized and small clusters of NM cells were filled with Alexa Fluor568 dextran (20 mm in sterile 0.9% saline; Invitrogen) by electroporation through a glass pipette positioned in NM [tip diameter 1–2 μm, adapted from Haas et al. (2001), Burger et al. (2005), and Sorensen and Rubel (2006)]. After electroporation, the slice was incubated for another 4–5 h at RT in oxygenated ACSF. The slices were subsequently fixed in 4% PFA for 10 min, rinsed in phosphate buffered saline (PBS), dehydrated in a series of ethanol steps and put into a clearing solution [3:5 mixture of Benzyl benzoate and Methyl salicylate (MacDonald and Rubel, 2008)]. Filled cells and axons could be imaged in their entirety in these 1000 μm tissue slabs using confocal microscopy (Fluoview 1000, Olympus). Brain slices containing filled cells and axon extensions to both ipsi- and contralateral NL were imaged by recording several Z-stacks of tiff-images. The resulting Z-stacks were loaded into Neurolucida (MBF Bioscience, Williston, VT) and the cell and its projections traced in 3-D over every Z-stack containing labeled ipsi- and contralateral projections to NL. The resulting 3-D trace was analyzed and measured in Neurolucida Explorer (MBF Bioscience).
Measurement of individual dorsal and ventral NM terminal fields on NL.
Individual dorsal and ventral NM terminal fields were measured in the 3-D tracings we described in the previous paragraph. On the ipsilateral side, we measured the axon length from the second branch point of the ipsilateral axon (BP-2i) (Fig. 1B) to the respective axon ending. On the contralateral side, we measured axon length from the point where the axon segment branched off the main axon to the axonal ending, i.e., we did not include the delay line in the measurements. All measurements were done in Neurolucida Explorer.
Additionally, previously published material of Young and Rubel (1983, 1986) was used to measure the axonal terminal trees of NM neurons onto NL cells. Briefly, horseradish peroxidase (HRP) tracer injections were made in the crossed dorsal cochlear tract (XDCT) of chick embryos at embryonic day 17 (E17). NM cells and fibers were filled with HRP and both ipsi- and contralateral terminal trees of NM axons onto NL were identified. Camera lucida drawings of HRP-labeled terminal fields of single NM axons were scanned and saved as digital files. All measurements were calibrated by scale bars depicted in respective images. Individual segments of axonal terminal trees to the ventral and dorsal of NL were traced and the lengths measured in ImageJ (NIH, Bethesda, Maryland). When fibers extended over several histological sections (50 μm section thickness), the Pythagorean theorem was applied to adjust the measured length to compensate for foreshortening.
In Figure 8A a schematic of the terminal tree demonstrates the pattern of terminal field branching of NM terminates onto NL dendrites. Branch Point 1 (BP-1) is the first bifurcation point. BP-2i is the second ipsilateral bifurcation point; BP-2c is the second contralateral bifurcation point (Fig. 1B).
Plastic embedding for determination of axon diameter.
Electron microscopy of ultrathin (0.1 μm) and light microscopy of semithin (1–2 μm) sagittal sections of chicken brainstem were used to measure axon diameter at the midline of the XDCT and through areas lateral of the midline that contain NM axons as they approach the ipsi- and contralateral NL. Chickens (P0/1) were deeply anesthetized with an intraperitoneal injection of pentobarbital and prepared for cardiac perfusion. Initial washout was made with 0.1 m cacodylate buffer (CB, pH 7.4) containing 0.001% CaCl2 and 4U/ml heparin. This was followed by perfusion with 2% PFA, 2.5% glutaraldehyde in CB. Chicks were decapitated, skulls opened and heads immersed in the same fixative for 1 h on a platform rotator and then overnight in fresh fixative at 4°C. Following 3 washes (10 min each) in CB, brains were dissected out and 2 mm slabs containing NM and NL were obtained by free-hand sectioning with a razor blade. Tissue was postfixed in 1% OsO4 in CB for 30 min, followed by 3 washes in CB, and embedded in Eponate. Semithin sections (1–2 μm) were stained with 1% toluidine blue in 1% sodium borate and trimmed to area of interest before ultrathin sectioning at 90 nm.
Ultrathin sections were mounted on 200-mesh thin bar Athene grids, stained with uranyl acetate and lead citrate and examined in a JEOL JEM 1200 EX II transmission electron microscope (JEOL). Digital images were obtained and axon diameter measured with an Olympus Morada camera (Olympus) and iTEM software (iTEM Software).
Semithin sections in the sagittal plane at the XDCT and NL were also imaged with light microscopy at 100× with a Marianas imaging system (Intelligent Imaging Innovations), incorporating a Zeiss Axiovert 200M microscope with an X–Z motorized stage (ASI) and a CoolSnap HQ digital camera (Princeton Instruments). We obtained high-power images of either NL and NM or XDCT using a 100× oil objective [numerical aperture (N.A.) 1.4] and Slidebook software with autofocus function (Intelligent Imaging Innovations). The resulting images were photomontaged in Photoshop (Adobe). The axon diameter and distance from the NL cell line were measured in ImageJ.
To avoid axon diameter measurements of fibers not cut exactly orthogonal to their long axis and to be conservative, we chose only fibers appearing round and measured diameter at the smallest value when the axon was slightly oval.
Internode distance.
In 8 chickens (P1 and P3) the distances between Nodes of Ranvier along individual NM axons were determined using an antibody against Caspr/paranodin/neurexin IV (NeuroMab clone K65/35, NeuroMab, University of California, Davis, Davis CA) to label paranodal proteins (Einheber et al., 1997). NM axons were labeled as described in the section above (“3-D tracing of individual NM axons and terminal trees”). The slices were then fixed in 4% PFA for 30 min at RT. After rinsing in PBS, the slices were resectioned at 50 μm on a Vibratome (Leica VT 1000S, Leica Microsystems). After further rinsing in PBS, the sections were immunostained with 1:500 NeuroMab clone K65/35 antibody in PBS with a standard block solution [0.3% Triton X-100 (Sigma), 5% normal goat serum] for at least 6 h at RT or overnight at 4°C. The sections were then rinsed in PBS (3×) and then exposed to the secondary antibody (1:200, Alexa Fluor 488 goat anti-mouse, Invitrogen) overnight at 4°C using the same block solution as with the primary antibody. After a final rinse in PBS (3×), the sections were coverslipped with Glycergel (Dako). Sections containing both labeled NM axons and paranodal staining were imaged using a confocal microscope (Fluoview 1000, Olympus) and the resulting 3-D images loaded into Neurolucida. Nodes of Ranvier colocalizing with dye-filled axons were identified and the internode distances measured in 3-D images. If necessary, the 3-D images were contrast-enhanced in ImageJ before the analysis in Neurolucida. The image in Figure 10B was surface-rendered with Huygens Essential (Scientific Volume Imaging).
Conduction time difference calculations.
To estimate actual conduction times between NM neurons and NL on each side of the brain, we used the average axon segment length data of each axon segment from our 3-D reconstruction measurements and divided them by sample conduction velocities (CV). Conduction times were then calculated as a function of each CV. Conduction time differences (NM to NL-ipsilateral vs NM to NL-contralateral) were established for XDCT conduction velocities of 2–50 m/s.
For the unadjusted conduction times, we assumed equivalent CVs in all NM axon segments. For the adjusted conduction times, CV was normalized to XDCT and CV in every other segment was calculated as a function of CV in XDCT multiplied with a CV correction factor. The correction factor was determined by the functions relating axon diameter and internode distance to CV (Rushton, 1951; Brill et al., 1977) of each segment relative to XDCT. All the segment conduction times for either the ipsilateral or contralateral axon length were added together and the difference in conduction times between the ipsilateral and contralateral pathway was determined.
Statistical analysis.
All our statistical analyses were done with Prism 5 software (GraphPad Software). Tests for statistical significance were made by t test and all error bars shown in figures represent SDs.
Results
Three major findings are presented below. First, the ipsilateral and contralateral axons from NM neurons to NL are not equal in length. This result does not support the common assumptions derived from the modified Jeffress model, that axons from the ipsi- and contralateral NM match length at the medial edge of NL. Second, measurements of axon diameter and internode distance at strategic positions along the NM axon reveal large and systematic differences in these parameters influencing conduction velocity. Third, calculations of conduction velocity (CV) using axon diameter and internode distance measurements demonstrate that the disparity in axon length can be compensated for by the parameter values that were measured.
3-D reconstruction of chicken brainstem nuclei and pathways
3-D reconstructions of NM and NL nuclei are illustrated in Figure 2(see also movie in supplemental data, available at www.jneurosci.org). NM is situated between the concave curvature of NL and the floor of the IVth ventricle. The ventral aspect of NM is separated from dorsal NL by ∼200–300 μm of neuropil. The entire NM lies slightly caudal relative to NL; the pair of nuclei is angled in the brainstem ventrocaudally to dorsorostrally (Fig. 2C) and its rostral ends are tilted medially. For most of its extent NL is a monolayer of bipolar cells that is curved, folded and positioned at an oblique angle. The NL resembles a saddle-shaped twisted plane (Fig. 2C). There is a genu close to the rostral end of NL with a bend of ∼70–90° that varies across animals (n = 3). A 3-D reconstruction of a P7 chicken NL was measured and the dimensions of NL are illustrated on a flattened schematic in Figure 3. The whole area of one NL nucleus is ∼1000 × 2000 μm; in the example shown the surface area is 1.89 mm2. The cell center to cell center spacing between NL cells was measured in 40 stained histological serial sections in two P15 animals. The average inter-cell distance (measured from cell center to cell center) was 27.4 μm (SD ± 8.7 μm) with no apparent rostrocaudal gradient in packing density (data not shown).
Integral tracing of NM neurons
To make an integral assessment of the axonal length of single NM neurons, fluorescent markers were electroporated into small clusters of NM cell bodies in thick slices, thereby labeling their axons. Figure 4A shows a montage image of multiple maximum intensity projection images, depicting a slice with an injection into the left NM. The labeled ipsilateral and contralateral branches of single cells were subsequently imaged, traced and measured. Three representative tracing examples are shown in Figure 4B. For each tracing, the cell body, the ipsilateral loop and terminals, as well as the contralateral axon and terminals can be identified. The numbers indicate the axon branch length (in μm) from the soma to the respective ipsilateral and contralateral terminal. In all cases the contralateral distance was measured to its most medial contralateral NL neuronal termination.
The main result of axon tracing measurements is presented in Figure 5A. Data are derived from 22 cells where axon length was determined from the NM soma to the contralateral NL. In 7 of these cells the length of the axon branch from the soma to ipsilateral axon endings was measured. Again, data for the contralateral axon is given as the distance to the most medial terminal identified and the length of the ipsilateral axon is given as the average across all terminals identified.
The distance from the soma to the point where the axons first bifurcates was determined as well (BP-1) (see Materials and Methods and Fig. 1B) and the distance from BP-1 to the most medial terminal identified in 15 cells (Fig. 5B). In 7 of those cells we additionally measured the axon segment length from BP-1 to the ipsilateral terminal (Fig. 5B).
The average length (± SD) for the ipsilateral NM axonal projection was 1480 ± 257 μm and the mean of the contralateral projection was 3116 ± 258 μm (Fig. 5C). Thus, the shortest (most medial) contralateral NM-NL projection averages >1600 μm longer than the corresponding ipsilateral projection. The distance from the soma to BP-1 has a mean of 378 ± 139 μm. Consequently, the mean value for the distance from BP-1 to the ipsilateral terminal is 1049 ± 166 μm and 2766 ± 200 μm to the most medial contralateral terminal (Fig. 5C). There is a profound and statistically significant length difference between the axon branch extending ipsilaterally and the branch extending contralaterally (Fig. 5C, p < 0.0001).
The 3-D reconstructions also showed that NM axons project overall slightly rostrally (Fig. 6). The ipsilateral loop originates in an NM neuron that is located more caudally in the brainstem than its terminal field, consistent to what the tonotopic organization of both nuclei would dictate (Rubel and Parks, 1975). The contralateral branch projects across the midline in the coronal plane. When it reaches the medial edge of NL it angles rostrally, to travel along an isofrequency contour of NL, confirming previous findings by Young and Rubel (1983).
Analysis of terminal trees
Axon terminal trees from 3-D reconstructions of NM axons in P0/1 animals were measured and compared with archival material published previously (Young and Rubel, 1983, 1986). Figure 7A shows the distribution of axon length in the ipsilateral terminal tree (BP-2i to terminal, mean 356 ± 118 μm, n = 37) and of axon length in the contralateral terminal tree, measured from the main axon branch coursing ventral to NL to the terminal (mean 168 ± 113 μm, n = 40) (Fig. 7B). The same analysis was obtained from the archival E17 data (Fig. 7A). Here the mean for the ipsilateral terminal tree is 337 μm (± 74 μm, n = 38) and for the contralateral terminal tree 147 μm (± 49 μm, n = 38) (Fig. 7B). Therefore, for both age groups the ipsilateral terminal trees are almost 200 μm longer than the contralateral terminal trees (p < 0.0001). There is no significant difference of terminal tree length between the two age groups (ipsilateral: p = 0.4142, contralateral: p = 0.4815).
The example shown in Figure 8A depicts a graphical schematic of the ipsilateral and contralateral terminal trees of an E17 animal in detail. Distances were measured from soma to BP-1, from BP-1 to BP-2i, from BP-2i to the ipsilateral terminals (orange) and from BP-2c to the contralateral terminals (white: delay line, main axon; green: contralateral terminal tree). The distance from soma to BP-1 was 354 μm. The distances from BP-2i to the terminal endings are plotted in Figure 8B and average 262 μm (± 36), and are very similar across the nucleus from the most medial to the most lateral NL neuron. The complete lengths of the ipsilateral pathways from BP-1 to the ipsilateral terminals are 848 μm on average. The ventral terminal arborization line (delay line), i.e., from BP-2c to the most lateral terminal, measures 885 μm (Fig. 8C). The projections from branching points along the delay line running parallel to the NL ventrally to the terminal endings average 136 μm (± 35) (Fig. 8C, green segments). As expected (Young and Rubel, 1986) the distance from BP-2c to terminal endings on the ventral side increases from medial to lateral.
Axon diameter
Since axon diameter influences conduction velocity and differences in diameter could compensate for differences in the NM axon path lengths to the ipsilateral and contralateral NL, NM axon diameter was measured in sagittal semithin and ultrathin sections in different segments of the NM axon (Fig. 9A). NM axon diameter was evaluated in the XDCT near the midline (Fig. 9B), within the ipsilateral terminal tree branching dorsal onto NL (Fig. 9C), and within the contralateral terminal tree branching ventral onto NL (Fig. 9D). Axon diameter measurements in the NL neuropil were also made as a function of distance from the NL cell line (Fig. 9G). Experiments in which we electroporated dye into the XDCT confirmed the homogenous population of axons, as only fibers extending to NM and NL were labeled (data not shown). We refrained from measuring axon diameter in the ipsilateral loop, i.e., in the segment from BP-1 to BP-2i. The mixed content of fibers in this area, containing auditory nerve axons as well as axons to the lateral lemniscal nuclei and midbrain, would make it impossible to get a clean database of NM to NL axons for this segment.
Three of the data groups shown in Figure 9 are from measurements made from semithin, osmium fixed, plastic embedded sections at maximum light microscopic resolution [100× (N.A. 1.4) objective]. To validate our measurement method using light microscopy, we measured axon diameter in electron microscopic images and light microscopic images at the midline (in the XDCT) from the same tissue (Fig. 9E,F; see Materials and Methods). We found no significant difference between the two groups (p = 0.7280) with similar values of means and SDs (EM: 1.917 ± 0.943 μm; light microscopy: 1.889 ± 0.933 μm).
Figure 9, E and F, presents the results of measurements of axon diameters in the XDCD and the terminal arbors approaching the dorsal and ventral neuropil regions of NL. Average axon diameters dorsal (1.104 ± 0.417 μm) and ventral of NL (1.470 ± 0.646 μm) are both significantly smaller than the mean axon diameter along the midline (p < 0.001). Dorsal axon diameters are also significantly smaller than diameters in ventral terminal tree (p < 0.001).
Axon diameter measurements in the NL neuropil as a function of their distance from the NL cell line (Fig. 9G) show a decrease in diameter as they approach NL (Fig. 9G). The two distribution correlations are significant (p < 0.001) and the regression line slopes are significantly different from each other (p < 0.05). However, if similar distance ranges (≤120 μm distance from NL cell line) are considered, i.e., datapoints of the ventral group with >120 μm distance from the NL cell line are excluded, the regression line slopes are not significantly different from each other (p = 0.19), while the distributions are (p < 0.05). The r2 value for the linear regression of ventral axon diameters within 120 μm from the NL cell line is 0.2491.
Internode distances
Another important factor regulating conduction velocity of axons is the distance between the nodes of Ranvier; the farther the distance between nodes, the faster the conduction velocity (Brill et al., 1977). Nodes in the NM axon were evaluated by filling axons with a fluorescent dye and counterstaining the resectioned tissue containing labeled axons with an antibody against a paranodal protein. Distances between the nodes were measured in 3-D using confocal microscopy along 5 different segments of the NM axon (Fig. 10A). Figure 10B shows an axon segment of the ipsilateral loop (axon: red, paranodal protein: green, white arrowheads indicate Nodes of Ranvier).
The measurements shown in Figure 10, C and D, provide the distributions of internode distance (C) and the means with SDs (D) (proximal segment 92.69 ± 32.42, ipsi loop 81.29 ± 41.54, ipsi terminal tree 61.59 ± 24.58, XDCT 157.8 ± 36.83, contra terminal tree 54.69 ± 28.29, all μm). These data reveal that the internode distances along the XDCT are significantly larger than all other segments (p < 0.0001), again indicating the fastest conduction velocity in the XDCT.
Conduction time difference
Next, we evaluated the data presented above for the difference in conduction time between the ipsilateral and contralateral axon. To accomplish this, axon diameter and internode distance were normalized to the values found in the XDCT and the adjusted conduction times calculated for all segments along the NM axon (see Material and Methods).
Figure 11A–D shows the calculated conduction times for each axon segment as a function of conduction velocities along XDCT. The conduction velocity for the contralateral axon along the ventral side of NL was measured in a chicken brainstem slice preparation by Overholt et al. (1992). Based on these data we assume that the conduction velocity will be between 3.0–8.8 m/s.
Figure 11, A and B, shows conduction times for an axon with uniform conduction velocity. Each axon segment contributes to the overall conduction time purely by its length. Figure 11, C and D, shows conduction times for an axon in which conduction velocities vary and are adjusted depending on axon diameter and internode distance. The graphs in Figure 11, A and C, show values for the ipsilateral branch and Figure 11, B and D, show values for the contralateral branch. Figure 11E is derived from this analysis and shows the difference in arrival time of action potentials between the ipsilateral and contralateral NL based on whether velocities are assumed to be equivalent in all axon segments of the NM axonal trajectory (blue line, unadjusted, data from Fig. 11A,B) or velocities are compensated by differences in axon diameter and internodal distances (red line, adjusted, data from Fig. 11C,D). The dashed horizontal green line shows the maximum range of physiological ITDs in chicks at ±180 μs, defined by the speed of sound, the inter ear distance and the interaural canal frequency-dependent effects (Calford and Piddington, 1988; Overholt et al., 1992; Hyson et al., 1994; Hyson, 2005).
Assuming equal conduction velocity throughout the entire length of the NM to NL axon branches, the difference between ipsilateral and contralateral conduction time is too large to be within the physiological relevant range. For conduction velocity values between 3.0 and 8.8 ms (Overholt et al., 1992) (Fig. 11, dashed vertical lines), the calculated difference in conduction time is entirely outside the physiological range of ITDs. With the adjusted conduction times from Figure 11, C and D, the difference in conduction time between ipsilateral and contralateral is substantially decreased. Here, the conduction velocities in the XDCT and along the NL neuropil differ. Therefore, based on measured conduction velocities in the NL neuropil (Overholt et al., 1992), as well as the differences in axon diameter and internode distance, it is assumed that the axons in the XDCT have a conduction velocity between 6.5 and 18.8 m/s (mean 11.9 m/s, solid vertical lines). Given these values, the difference between ipsilateral and contralateral conduction time falls into the physiological range (Fig. 11E). Hence, the variations in axon diameter and internode distance can provide a means to compensate for the axon length disparity.
Discussion
In this study we measured axon length in the chicken brainstem sound localization circuit responsible for the processing of ITDs, as well as parameters responsible for conduction velocity, specifically axon diameter and internode distance. Surprisingly, axon length dimensions do not to comply with the proposed modified Jeffress model (Overholt et al., 1992), as the contralateral axon is >1600 μm longer than the ipsilateral one, making coincidence detection based on travel distance impossible. However, the data indicate that variations in both axon diameter and internode distances may counterbalance the offset in axon length. This axon length disparity and anatomical evidence for temporal compensation provokes a reassessment of current thinking about the organization of the avian and mammalian sound localization circuits responsible for ITD coding.
In accord with the hitherto proposed model embodied by the NM–NL circuit (Overholt et al., 1992), our results illustrate similar lengths of axons from NM onto the dorsal NL dendrites and confirm that each of the contralateral axons forms successive series of short collaterals stemming from a parent axon coursing immediately ventral to NL neurons orthogonal to the tonotopic axis. These collaterals form a “delay line” of axon terminals on NL neurons from the contralateral NM. The most interesting and unforeseen observation of this study is a large and reliable difference in the lengths of contralateral and ipsilateral pathways from NM to NL, which does not comply with the proposed modified Jeffress model. In their seminal paper about the barn owl sound localization circuit, Carr and Konishi (1990) also noted a difference in axon lengths from NM to the ipsilateral and contralateral NL. The differences (∼1 mm, and a much larger head) were not as great as those reported here, but that may be due to the fact that three-dimensional reconstructions were not attempted. In this paper, Carr and Konishi (1990) also foreshadowed our findings of variations in axon caliber and internodal differences, noting that these two parameters varied along the pathways and may modulate conduction velocities in meaningful ways. Physiological data demonstrate the existence of neurons responding best to 0 ITD (Köppl and Carr, 2008). Indeed, rough calculations of the conduction velocity for the ipsilateral and contralateral axon using response delay data from Köppl and Carr (2008), their Figure 4 and from Köppl (1997), her Figure 10 yield values that average at 3.8 and 8 m/s, respectively. Not only are the resulting values for the contralateral axon twice as large, the conduction velocity for the contralateral axon is also consistent with the values proposed for the XDCT in this paper (Fig. 11E). Given the difference in length between ipsi- and contralateral axon collaterals, the encoding of ITDs in the physiological range for the chicken is impossible unless other morphological or physiological features regulate the timing in these circuits. In myelinated axons, their diameter is linearly proportional to conduction velocity (Gasser and Grundfest, 1939; Rushton, 1951; Hutchinson et al., 1970). Our data show significantly smaller axon diameters in the terminal trees compared to the major branches of the axons as they traverse the XDCT. Additionally, the ipsilateral terminal tree has an average axon diameter that is only three quarters of the average diameter in the ventral terminal tree, further delaying the arrival of ipsilateral action potentials at NL relative to the contralateral signal. This effect is amplified by the fact that the average axon length in the ipsilateral terminal tree is more than double the length of the mean axon in the ventral terminal tree. We did not measure axon diameter in the ipsilateral loop, but variations in this segment might enable the system to further alter conduction velocities and to slow down the signal on the ipsilateral side.
Brill and colleagues showed that conduction velocity increases linearly with internode distance to up to 2000 μm (Brill et al., 1977). Moreover, when the ratio of internode distance (L) and axon diameter (d) is small (<150), conduction velocity is quite sensitive to variations in L/d (Brill et al., 1977). Carr and Konishi (1990) compared internode distance of axons as they traverse the cell body layers within the owl NL with axons outside NL, and also report shorter internode distances of the axons within NL. At that time, methods to label individual axons were not readily available; hence the origin of the axons they measured could not be confirmed. We were able to measure internode distances in different segments of labeled NM axons and found significant variations along the NM axon. In particular, internode distances along XDCT is significantly larger compared to other parts of the NM axon by at least 50%. It also appears that internode distances within the terminal arbors on both the ipsilateral/dorsal axons and the contralateral/ventral axons are smaller than the main branches leading up to these terminal branches. Again, this disparity of terminal tree length enhances the effect of short internode distances along the terminal endings and contributes to delaying the propagating signal ipsilateral relative to the contralateral side of the brain. Taking into account the axon diameter variations and the resulting sensitivity to L/d ratio, there seems to be ample possibility for variations in signal speed.
In summary, our results show that at least two parameters influencing conduction velocity, axon diameter and internode distances, are regulated at different sites within individual axons of NM neurons to specifically adjust the conduction velocities of the ipsilateral and contralateral signals, presumably to optimize coincidence detection. Slowing down the propagating signal in the ipsilateral part of the axons allows the contralateral signal to “catch up” and enables encoding of ITDs in the physiological range in NL neurons. The regulation of these axonal parameters within individual axons seems quite remarkable from a cell biological point of view, but it is not unprecedented. Measurements in other neuronal systems have been shown to provide temporal compensation for conduction distance disparities. For example, in the electromotor system in fish internode distance is adjusted so that conduction velocity compensates for different axon lengths (Bennett, 1970). In rat Purkinje cells, conduction velocity is varied to enable isochronicity of different length climbing fibers (Sugihara et al., 1993). Other parameters, such as myelin sheet thickness (Rushton, 1951), may also influence signal propagation in the NM-NL circuit. Moreover, one might inquire if there is a period of adjustment for one or several of these morphological features during development, wherein these circuits are temporally sharpened and compensated for changes in head size to optimize discrimination of binaural timing differences.
We show that axon length alone cannot be solely responsible to compensate for external ITDs in the chicken sound localization circuit. This unexpected result has implications for the mammalian brainstem binaural system as well. It is noteworthy that the mechanism responsible for low-frequency sound localization in mammals has recently come under dispute. There is only marginal anatomical evidence for delay lines in mammals (Smith et al., 1993; Beckius et al., 1999) and a proposed alternative mechanism, incorporating the glycinergic inputs to MSO, requires a very fast inhibitory input (τdecay = 0.1 ms) (Brand et al., 2002) that has not been confirmed by physiological measurements (Magnusson et al., 2005). Variations of parameters such as axon diameter, internode distance, and others (Pecka et al., 2008) in the mammalian brainstem might be responsible for precise adjustments of physiological delays, thereby creating the framework and adjustments of the ITD detection circuit.
Footnotes
-
This work was supported by National Institutes of Health–National Institute on Deafness and Other Communication Disorders Grants DC003829, DC004661, and DC008042. We thank Glen McDonald and Dale Cunningham for technical assistance. Dale Cunningham was pivotal in preparing the semithin and ultrathin sections. We also thank Rock Levinson for suggesting the caspr antibody. We are grateful to Yuan Wang, Kathryn Tabor, Melissa Caras, Yong Lu, R. Michael Burger, Jason Tait Sanchez, and Daniel T. Kashima for discussion and criticism on the manuscript. We thank the anonymous reviewers for useful comments and constructive criticism. Parts of this study were presented previously (Harris et al., 2008).
- Correspondence should be addressed to Edwin W Rubel, Virginia Merrill Bloedel Hearing Research Center, University of Washington, Box 357923, Seattle, WA 98195-7923. rubel{at}uw.edu