Toric Spines at a Site of Learning

<jats:title>Abstract</jats:title><jats:p>We discovered a new type of dendritic spine. It is found on space-specific neurons in the barn owl inferior colliculus, a site of experience-dependent plasticity. Connectomic analysis revealed dendritic protrusions of unusual morphology including topological holes, hence termed “toric” spines (<jats:italic>n</jats:italic>= 76). More significantly, presynaptic terminals converging onto individual toric spines displayed numerous active zones (up to 49) derived from multiple axons (up to 11) with incoming trajectories distributed widely throughout 3D space. This arrangement is suited to integrate input sources. Dense reconstruction of two toric spines revealed that they were unconnected with the majority (∼84%) of intertwined axons, implying a high capacity for information storage. We developed an<jats:italic>ex vivo</jats:italic>slice preparation and provide the first published data on space-specific neuron intrinsic properties, including cellular subtypes with and without toric-like spines. We propose that toric spines are a cellular locus of sensory integration and behavioral learning.</jats:p>

section_content: Neuronal computations are shaped by postsynaptic morphology (Rall, 1974). Across circuits and species, dendritic spines receive the majority of excitatory synaptic input (Hering and Sheng, 2001). Most spines have narrow necks that serve to compartmentalize, electrically (Yuste, 2013; Kwon et al., 2017) and biochemically (Müller and Connor, 1991; Sabatini et al., 2001), a small number of synaptic inputs -usually one. In this arrangement the dendritic branch performs the first stage of integration of inputs coming from multiple sources (Poirazi et al., 2003). Here, we identify an alternative motif involving first stage integration via a new type of dendritic spine. 

Space-specific neurons in the barn owl inferior colliculus (Knudsen and Konishi, 1978) compute sound source direction via integration of binaural cues, interaural time difference (ITD) and interaural level difference (ILD), and the elimination of phase ambiguity by convergence across frequency channels (Wagner et al., 1987; Olsen et al., 1989; Brainard et al., 1992; Konishi, 2003; Singheiser et al., 2012), computations known to involve multiplication and non-linear thresholding, respectively (Peña and Konishi, 2000, 2001, 2002, 2004). In addition, the responses of space-specific neurons can be modulated by visual input (Bergan and Knudsen, 2009) and attention (Winkowski and Knudsen, 2006; Mysore and Knudsen, 2014). Thus, space-specific neurons act as pattern detectors via the integration of complex sensory inputs, and make response adjustments based on behavioral context. 

Space-specific neurons provide an opportunity to understand how neuronal computations are shaped by microanatomical patterns of synaptic convergence, but their ultrastructure is unknown. To fill this gap we used serial block-face scanning electron microscopy (SBEM; Denk and Horstmann, 2004) and stimulated emission depletion (STED) microscopy (Hell, 2003). Space-specific neurons in young adult owls were labeled in vivo by focal injection of retrograde tracer at functionally defined map locations (Rodriguez-Contreras et al., 2005). STED imaging revealed two cellular subtypes, one studded with typical spines and the other with a previously unreported type of spine with unusually complex and variable morphology, termed "toric" spines. SBEM imaging was used to reconstruct all 76 toric spines found on the soma and proximal dendrites of one space-specific neuron, along with their associated presynaptic terminals. In length, cytoplasmic volume and innervation density, toric spines resemble thorny excrescences found on pyramidal cells in mammalian hippocampus and amygdala (Amaral and Dent, 1981). Thorny excrescences receive numerous active zones (Chicurel and Harris, 1992) coming from one or at most two input sources (Wilke et al., 2013; e. g., individual mossy fibers in the hippocampus), forming the so-called detonator synapse (Urban et al., 2001). We found that toric spines receive a comparable number of active zones, but in contrast, these come from up to 11 different inputs. The latter arrangement is suited for integration, not detonation nor compartmentalization. 

To date, knowledge of space-specific neuron function derives from extracellular recordings (numerous studies) and three reports using intracellular sharp electrodes (Peña and Konishi, 2001, 2002, 2004). Neither technique provides information on intrinsic electrical properties or can resolve individual synaptic events. We performed ex vivo patch-clamp recordings from putative space-specific neurons in brain slices from juvenile owls. Many of these cells exhibited large, atypical spines including toric-like structures, the apparent developmental precursors to mature toric spines. In total, these findings and newly established methods provide a path to investigate the causal connections between microanatomical structure and the neuronal computations that underlie high-level pattern detection.
section_name: Introduction

section_name: Materials and Methods

section_content: A total of 12 barn owls (Tyto alba) of either sex were used in this study. Adult animals (Ͼ250 d old) were grouphoused in large flight aviaries while juveniles (35-40 d) remained in isolated nest boxes with parents and siblings. All husbandry and experimental methods were approved by the University of California Davis Institutional Animal Care and Use Committee.
section_name: Animals

section_content: Four adult owls were used for in vivo retrograde labeling. Owls were anesthetized using 2% isofluorane, a mixture of nitrous oxide and oxygen (1:1) and wrapped in soft cloth. Craniotomies were opened above the tectal lobes and the owl was secured in a stereotax located inside a soundproof recording room. Insulated tungsten recording electrodes ϳ6-M⍀ impedance (FHC) were lowered through the forebrain to the optic tectum (OT) based initially on stereotaxic coordinates (17 mm rostral to the neck muscle insertion, 6 mm lateral to the midline) and then refined on the basis of response criteria. Visual receptive fields were measured by computer-controlled projection of static or moving dots onto a calibrated screen and the responses were used to navigate to the deep layers of OT (dOT; in 5/6 injections) based on characteristic latencies, preference for negative contrast stimuli, and absence of bursting behavior. Sites representing frontal space near 0°azimuth were targeted. ITD tuning was determined via dichotic stimulation delivered through speakers (ED-21913-000, Knowles Electronics) positioned 5 mm from the tympanic membrane. Tungsten electrodes were removed and replaced with glass electrodes in the same location (1. 5 mm borosilicate glass, 10 -20 m tip), containing 10% biotinylated dextranconjugated tetramethylrhodamine 3000 MW (microruby, #D-7162; Thermo Fisher) or 10% lysine fixable dextran 3000 MW (Texas-red, #D3328; Thermo Fisher) in 1% potassium chloride solution. Chloridized silver wires were used to record physiologic responses to confirm injection location. Ionophoretic injections were performed using 7-s on/off cycles of 3 A for 15 min. Glass electrodes were withdrawn and the procedure was repeated on the opposite tectal lobe. In 1/6 injections, the external nucleus of the inferior colliculus (ICX) was directly targeted (Table 1 ). Following injections, craniotomies were disinfected with 1% chloroptic ointment. After 3-to 7-d survival time, owls were anesthetized using 5% isofluorane ϩ 1:1 nitrous oxide/oxygen. Heparin (300 U) was injected into the left cardiac ventricle and the owl was perfused with 0. 1 M phosphate buffer (PO4, pH 7. 4), followed by 4% paraformaldehyde in PO4 buffer (for STED imaging; see below for tissue preparation for SBEM imaging). Brains were removed and placed in 30% sucrose and 0. 1 M PB for 24 h. Tectal lobes were isolated, sectioned at 50 m on a vibratome (Leica), and mounted with an antifade agent (Prolong gold, P36930; Thermo Fisher).
section_name: Microelectrode recording and retrograde labeling in vivo

section_content: Imaging was performed on a Leica TCS SP8 STED 3ϫ confocal microscope (Leica Microsystems), using a 100ϫ oil-immersion objective (HC PL APO CS2, 1. 4 NA; Leica Microsystems). Parameters for STED imaging were white light laser, line 551 nm for stimulating, and 660-nm depletion laser, both between 20% and 40%. Image frame size was 4800 ϫ 4800. 

Step size for z-stacks was 17. 44 m, resulting in pixel/voxel dimensions of 0. 024, 0. 024, and 0. 18 m (x, y, and z, respectively). Huygens Professional SVI software was used for deconvolution of all images, using a theoretical point spread function (PSF). All well-labeled neurons in the rostrolateral aspects of IC were imaged (n ϭ 39 cells in perfusion-fixed tissue from three adult owls; n ϭ 16 cells in immersion-fixed tissue from five juvenile owls as described later). 

STED images were analyzed using Imaris 8. 1 (Bitplane AG). Surface creation was used to measure surface area and volume of the soma and primary dendrites. To isolate the primary dendrite, the slice tool was used to separate the base of the dendrite from the soma at the point where the dendritic caliber first became uniform (i. e., avoiding sharp curvature at the interface). The cross-sectional surface of the ensuing dendritic face was used to measure thickness of primary dendrites. 

Cell classification was accomplished using Imaris' 3D view. For adult tissue, cells containing toric spines were classified as Type I, while cells lacking toric spines were classified as Type II. For juvenile tissue the cell population was more heterogeneous. Cells were not classified as Type I or Type II but scored for the presence or absence of toric-like spines.
section_name: STED imaging and morphometric analysis

section_content: In one owl, microelectrode recordings were made as described above followed by multiple tracer injections in deep OT to increase the likelihood of generating a completely filled neuron. The owl was perfused as before but with 1% glutaraldehyde in cacodylate buffer (Wilke et al., 2013) ; 150-m vibratome sections were collected and shipped to University of California San Diego for further processing and imaging. Sections were freeze-fractured, reacted en bloc with streptavidin-HRP (which binds to the biotin moiety on microruby) and detected using diaminobenzidine (DAB). Light microscope image (LM) of the section at this point was obtained. Following osmication and staining with heavy metals the tissue was extremely opaque to visible light, obscuring internal structure. To target the labeled cell for SBEM volume imaging, a Zeiss Versa 510 x-ray microscope (XRM) was used to generate 3D tomographic volumes relative to tissue landmarks with 0. 7-m pixel resolution. Based on the XRM volume, a 45 ϫ 75 ϫ 75 m tissue block centered on the labeled cell was targeted and imaged with a Gatan 3View unit on a Zeiss Merlin SEM at 2. 5 kV, with a magnification of 2800ϫ, 15 ϫ 15 k raster size, and 1-s pixel dwell time. The final pixel size was 5 nm and the Z step size was 70 nm. The imaging parameters yield sufficient resolution to identify synaptic contacts while maximizing the rate of data collection. The SBEM image volume was postprocessed for proper alignment at University of California, San Diego and annotated at University of California, Davis.
section_name: SBEM

section_content: A team of eight annotators used the open source platform IMOD (https://bio3d. colorado. edu/imod/) to analyze the SBEM image volume. It contained 54 cell nuclei representing six distinct cell types (neuronal and glial somatic profiles), four putative space-specific neurons, one of which was well labeled throughout entire cytosolic compartment including spines, an estimated 100s of myelinated axons, 1000s of unmyelinated axons and 10,000s of synapses. The labeled space-specific neuron was traced through all sections, meshed and rendered in 3D. This reconstruction includes the entire soma and proximal dendrites but not distal dendrites as they extended outside of the volume (space-specific neuron dendritic fields can be ϳ200 m in diameter). Each toric spine was numbered and isolated for further analysis. The surface areas and contour volumes were determined using the imodinfo command, which calculates surface area by adding the areas of the triangles making up the mesh and the volume by summing the areas of each contour times the Z distance connecting the contours. The number of holes was assessed by visual inspection of the 3D model. Because the cytosolic label was dense, in most cases it was not clear whether these tori were cytoplasmically contiguous or interrupted by tight junctions, although one clear example of each was found. Twenty-seven toric spines representing the full spectrum of complexities were selected for further analysis. All unlabeled axons innervating these 27 spines were traced outwards from their points of synaptic contact until the axonal process became ambiguous or reached the end of the volume.
section_name: SBEM analysis

section_content: Synapses were identified by the presence and distribution of presynaptic vesicles. Active zones were identified using the following criteria: (1) a cluster of five or more presynaptic vesicles, (2) vesicle cluster apposed to the presynaptic membrane, (3) parallel presynaptic and postsynaptic membranes (no gaps), (4) these features present in at least two consecutive sections. In some cases, the postsynaptic density (PSD) could be discerned but the presence of the dark DAB label made this an unworkable requirement.
section_name: Synapse identification in SBEM

section_content: Seven spines whose axons had been extended for longer distances without ambiguity were selected for trajectory analysis. The length, number of boutons and outgoing trajectories at each axonal end were measured. Length was measured by manual skeletonization. Boutons were defined as axonal swellings containing presynaptic vesicles. Trajectories were determined by marking the xyz-coordinates of the skeleton endpoint as well as a second internal skeletal point within 2 m of the skeleton endpoint. A small number of axons branched and were split into segments. If the axon had a clear main segment with a smaller side branch it was classified as a "T" branch and the axon was split into two segments; the main segment and side branch. The length and number of boutons were measured separately for each segment. For trajectory analysis, only the two ends of the main segment of these axons were included. If the axon branched into two pieces that were comparable and a "main segment" could not be distinguished, it was classified as a "Y"branch and split into three segments at the branch point. End 1 of each segment was defined as the open end, and end 2 was the branch point. Only one analyzed axon was classified as a Y-branch. For this axon, the trajectory was measured at all three open ends. 

Axonal trajectories were analyzed by custom scripts written in MATLAB (MathWorks). The two points at each end were converted to a 3D vector with the outside point as the vector head and the inside point as the vector tail. All z-coordinates were scaled by a z ratio of 13. 78 pixels in x-y/z section. Each end vector was scaled to a length of 1 and plotted in 3D with its tail at the origin. End vectors were grouped together by spine, with the two vectors from the same axon paired together. To determine whether input trajectories were clustered, the distance between vectors was measured as the minimum angle of arc drawn from the head of one vector to the head of the other. Vector distance () was calculated from the dot product of the Euclidean vectors:

Distances were calculated pairwise between each end vector and every other end vector on the same spine, except the other end of the same axon. The distance between the two end vectors from the same axon was called "angle of separation" to differentiate from other vector distances. This is a proxy measure for curvature of the axon. A perfectly straight axon yields an angle of separation of 180°. The single Y-branched axon on TS24 was included in the vector plots but excluded from the distance calculations for simplicity.
section_name: Axonal trajectory analysis

section_content: Potential connectivity for TS1 and TS7 was assessed by partial reconstruction of all unmyelinated axons that passed within 2 m of the spine. CF was calculated as (the number of axons that synapse with the spine/the total number of axons within 2 m). To determine the probability envelope of observing the actual distribution of active zones across each potential axonal input, a bootstrap analysis was constructed based on the simple assumption of equal access. The simulation was run 10,000 times for each spine to determine mean values and Ϯ2 SD.
section_name: Connection fraction (CF)

section_content: Six juvenile owls (35-40 d old) were used for current clamp experiments to determine intrinsic properties and two juveniles (40 -50 d old) were used for voltage-clamp experiments to measure EPSCs. Juvenile owls were anesthetized with isoflurane and perfused transcardially with heparin-containing (300 units/l) artificial CSF (aCSF; 124 mM NaCl, 1. 3 mM KCl, 1 mM NaH 2 PO 4, 26. 2 mM NaHCO 3, 1. 3 mM MgCl 2, 2. 5 mM CaCl 2, and 11 mM glucose). Brains were rapidly removed, blocked, and placed in ice-cold modified aCSF (248 mM sucrose, 11 mM glucose, 1. 3 mM KCl, 1 mM NaH 2 PO 4, 26. 2 mM NaHCO 3, 1. 3 mM MgCl 2, and 2. 5 mM CaCl 2 ). Horizontal midbrain sections (250 M) containing ICX and lateral shell of the central nucleus of the inferior colliculus (ICCls) were cut on a vibratome (Leica VT1200S) and transferred to an incubation chamber containing aCSF at 32°C for 30 min before moving to aCSF at room temperature until used for recordings. Recordings were made in a submersion chamber perfused with aCSF (2 ml/min) at room temperature. All solutions were bubbled with 95% O 2 -5% CO 2 continuously.
section_name: In vitro slice preparation

section_content: Whole-cell patch-clamp recordings were made from visually identified cells in the ICX region using borosilicate glass pipettes (2-4 M⍀). For current-clamp experiments, pipettes were filled with K-methanesulfonate-based internal solution (135 mM CH3KO3S, 10 mM HEPES, 1 mM EGTA, 4 mM Mg-ATP, 0. 4 mM Na-GTP, and 10 mM phosphocreatine; 310 mOsm) containing 0. 5% biocytin (Thermo Fisher). For voltage-clamp experiments, a Csmethanesulfonate-based internal (120 mM CH 3 CsO 3 S, 15 mM CsCl, 8 mM NaCl, 10 mM TEA-Cl, 10 mM HEPES, 0. 5 mM EGTA, 2 mM QX314, 4 mM Mg-ATP, and 0. 3 mM Na-GTP; 310 mOsm) was used. All chemicals were from Sigma. Recordings were acquired using a Multiclamp 700B amplifier (Molecular Devices), digitized at 20 kHz with a Digidata 1550 digitizer (Molecular Devices), and low-pass filtered at 8 kHz. Evoked EPSCs were elicited in cells voltage-clamped at -70 mV in response to a brief electrical stimulus (0. 2 ms, 20 -50 A) delivered through a concentric bipolar stimulating electrode (FHC). The stimulating electrode was placed locally in the ICCls. Synaptic blockers (bicuculine: 20 M; NBQX: 10 M; both from Sigma) were washed in during recordings, as indicated. Slices were immersion-fixed in 4% paraformaldehyde and stained with SA-488 for STED imaging.
section_name: Patch-clamp recordings

section_content: Resting membrane potential (in mV) was measured while injecting no current (I ϭ 0) immediately after breaking into a cell (average of first 500 ms). Membrane capacitance was estimated using a single-exponential fit to the first 200 ms of the response to a 0. 5-s injection of hyperpolarizing current. Membrane resistance was measured from the linear portion of a current clamp input-output curve as the slope of the line fitted to the voltages recorded at increasing current injections. Action potentials were elicited by increasing depolarizing steps (0. 1-3 nA, 0. 5 s) delivered every 15 s. Spikes were identified using a 20 mV/ms threshold of the first-order derivative of membrane potential. Results were identical when a more conservative 50 mV/ms detection threshold was applied. Spike onset was determined using the 2nd-order derivative (maximum acceleration) of membrane potential within a 2-ms time window before and after threshold crossing. Spike peak voltage (V P ) was defined as the maximal V within a 2-ms time window from threshold crossing and was used to calculate spike amplitude (peak V minus V at spike onset). First spike latency was defined as the time between stimulus onset and first spike onset, and inter-spike interval as the time between consecutive spike onsets. All data were analyzed with custom scripts written in MATLAB.
section_name: Data analysis

section_content: The strategy for labeling space-specific neurons is shown in Figure 1. Tungsten and glass electrodes were used to record multiunit activity in the tectal lobes, which include the inferior colliculus (IC) and OT. The lobes are shown in a magnetic resonance image of an intact owl brain at mid-transverse plane (Fig. 1A ). Space-specific neurons found in the ICX are tuned for distinct values of ITD, as illustrated in Figure 1B. Their auditory spatial receptive fields are constructed from convergent inputs which originate in the ICCls and other structures not shown (see Discussion). Space-specific neurons were labeled by injection of the tracer microruby or Texas-red in the dOT (5/6 injections), which receive monosynaptic input from space-specific neurons. Injection sites were confined to the deep layers with minimal rostrocaudal spread (Fig. 1C, D ). ITD and visual spatial tuning at the injection sites (Fig. 1E ) represented frontal space; as expected, the retrogradely labeled neurons were found within the region of the inferior colliculus representing frontal space, shown by the white circle in Figure 1C. This region includes the full mediolateral extent of ICX as well as lateral aspects of ICCls. Because no cytoarchitectonic marker that delineates the border between ICCls and ICX (Takahashi et al., 1987; Wagner et al., 2003) was compatible with the three tissue preparation methods used in this study (glutaraldehyde perfusion, 150-m section; paraformaldehyde perfusion, 50-m section; in vitro slice, 250-m section), neurons were not assigned to one or the other structure. In one case the injection was directly targeted to ICX; labeled neurons resulting from this injection were indistinguishable in both location and morphology from those labeled by the five dOT injections and therefore included in the analysis below. No neurons outside of the ICCls-ICX microcircuit were analyzed. 

The tracer often produced complete cell-fills, revealing fine aspects of dendritic structure, similar to previous reports using in vivo labeling (Peña and Konishi, 2001) or immunostaining for CaMKII (Rodriguez-Contreras et al., 2005; Niederleitner and Luksch, 2012). From these studies, it was known that space-specific neurons are large multipolar cells with dendritic protrusions which often included typical spines. Unlike previous studies, we used high-resolution methods to reveal, for the first time, their ultrastructure. STED microscopy of 39 well-labeled space-specific neurons derived from five injections in three owls revealed two cell classes within a spectrum of overlapping morphologic features. The most striking difference was the presence of highly atypical dendritic spines on a subset of cells, defined here as Type I spacespecific neurons. These large structures exhibited tubular build, tended to lack spine heads and were often riddled with holes (Fig. 2A ). The morphologic diversity of toric spines was so large that no single feature distinguishing them from typical spines was present in all. Type I space- specific neurons were studded with toric spines and a paucity of typical spines, sometimes none (Fig. 2A ). Type II space-specific neurons, defined here, were studded with a high density of typical spines and devoid of toric spines (Fig. 2B ). In addition, Type I space-specific neurons tended to have larger somas (Fig. 2C, F, G ) and thicker primary dendrites (Fig. 2D, H ). The number of primary dendrites was not significantly different (Fig. 2E, I ). 

While the resolution of STED was sufficient to identify toric spines when they presented in the right orientation, high anisotropy (ϳ40 nm lateral, ϳ600 nm axial resolution) introduced unacceptable distortions to the apparent 3D structure of spines presenting in other orientations (the majority of spines). Combined with high morphologic diversity and narrow tubular structure, volumetric measurement (toric vs typical) were not informative for the population of spines as a whole. We therefore turned to SBEM. Space-specific neurons were labeled as before, this time with multiple injections in the deep OT to increase the chance of detection, the brain perfused with glutaraldehyde and reacted en bloc with DAB to produce a durable reaction product that survived all subsequent processing. One well-labeled space-specific neuron was found on the ICCls-ICX border (Fig. 3, LM panel). After osmication the labeled space-specific neuron was isolated using x-ray microscopy (Fig. 3, XRM panel) and a 75 ϫ 75 ϫ 45 m volume (253,125 m 3 ) imaged at 5 nm lateral resolution. 

The DAB reaction product completely filled the neuron, which received profuse synaptic contacts from many unlabeled afferents (Fig. 3, SBEM panel). Volumetric reconstruction revealed a large multipolar neuron (Fig. 4 ) resembling the Type I space-specific neurons identified by STED. In particular, the soma and its six thick primary dendrites were studded with 76 toric spines (TS1 -TS76) and devoid of typical spines. Comparison of representative "zoos" of toric spines identified by SBEM or STED (Fig. 5 ) indicated that these populations are qualitatively indistinguishable. A thin process resembling an axon emerged from near the base of one primary dendrite. No somatic axon was found. 

To map synaptic input, all active zones were identified by a team of annotators using criteria defined in Materials and Methods. Active zones occurring on the soma, dendrites and toric spines appeared similar as a relatively homogenous population of chemical synapses with clusters of synaptic vesicles apposed to the postsynaptic membrane (Fig. 6 ). Bifurcated synapses with two or more active zones belonging to the same axonal bouton were occasionally observed on all compartments. 622 active zones were found on toric spines and 754 on the soma and dendrites, indicating that toric spines are a major input hub (Fig. 7A-C ). The innervation density (# active zones per surface area) of toric spines was on average higher than that observed on the six primary dendrites, soma, or thin process (Fig. 7D ). 

Quantification of the ultrastructural features of toric spines is shown in Figure 8. Volumes (mean ϭ 1. 99 ϫ 10 9 nm 3 ) and surface areas (mean ϭ 8. 20 ϫ 10 7 nm 2 ) varied over a wide range and were strongly correlated, the latter consistent with structures built largely of narrow tubes (Fig. 8A ). Holes (up to seven on a single spine) were found on 37% of spines and were more prevalent in larger spines (Fig. 8B ). Some holes were formed by tight junctions between two membrane-bound filopodial extensions while others appeared as cytoplasmic continuities (Fig. 9 ). The number of active zones per spine ranged from 1 to 49 and was strongly correlated with size (Fig. 8C ) and number of holes (Fig. 8D ). These data demonstrate large morphologic diversity within the toric spine population. 

To determine patterns of synaptic convergence, 27 spines representing the morphologic spectrum were se-lected for reconstruction of axonal inputs. Three examples are shown in Figure 10. In each case, the active zones derived from multiple axons: 10 axons provided 25 AZs to TS1; three axons provided six AZ to TS2; and 10 A, Volume versus surface area. Data were marginally better fit (dashed line) by a power law ( 2 ϭ 8. 15 ϫ 10 14 ) than linear regression ( 2 ϭ 9. 28 ϫ 10 14 ; r 2 ϭ 0. 97, p Ͻ 0. 00001. B, Holes versus surface area: marginally better fit (dashed line) by power law ( 2 ϭ 90. 5) than linear regression ( 2 91. 1; r 2 ϭ 0. 42, p Ͻ 0. 00001). C, Active zones versus surface area: marginally better fit (dashed line) by linear regression ( 2 ϭ 1254; r 2 ϭ 0. 73, p Ͻ 0. 00001) than power law ( 2 ϭ 1282). D, Active zones versus holes: best fit by linear regression (r 2 ϭ 0. 50, p Ͻ 0. 00001). There was a weak tendency for spines located on the soma to be larger than those located on dendrites, in both volume (mean values 3. 13 vs 1. 62 ϫ 10 9 nm 3 ; p ϭ 0. 026 Mann-Whitney U test) and surface area (mean values 2. 75 vs 1. 50 ϫ 10 7 nm 2 ; p ϭ 0. 016 Mann-Whitney U test). Spines on the thin process tended to be smaller though there were not enough data points to make a statistical comparison. axons provided 23 AZ to TS7. In fact, the majority of toric spines (all except three of the 27 representative spines) received more than one axonal input. The number of axonal inputs strongly correlated with the number of AZs, and this was true for spines located on the soma or dendrites; spines on the soma tended to receive more inputs than those on the dendrites or thin process (Fig. 11A ), consistent with their tendency to be larger. 

To determine the degree to which inputs to a single spine were balanced or unbalanced in anatomical strength we measured the input fraction, defined as the number of active zones made by a single input divided by the total number of active zones received (metric range ϭ 0 -1). A low fraction indicates the input has a weak contribution to the overall innervation of the spine, while a high fraction indicates the input is a dominant driver. The distribution of input fractions was skewed low (mean ϭ 0. 20, SD ϭ 0. 19, n ϭ 138) indicating the most common arrangement was integration of many individually weak inputs, a situation of anatomical balance (Fig. 11B ). Large input fractions were mostly limited to spines that received few inputs (Fig. 11C ) and active zones (Fig. 11D ), both of which set a high lower limit on the possible values of input fraction. Overall, the common motif was anatomically balanced integration. 

To determine whether multiple inputs derived from the same parent axon, seven spines whose axons had been extended for considerable distance without ambiguity were selected for trajectory analysis (mean length ϭ 20. 6 m; see Materials and Methods; Fig. 12A, B ). None of these 40 axons joined together. Further, the population of trajectories was distributed near randomly in 3D (Fig. 12C ), and this divergent profile (Figs. 12D, 13 ) was evident at individual spines (with one exception) including all of the most complex spines. The axonal segments tended to exhibit minimal curvature between incoming and outgoing ends (Fig. 14 ). In total, wide divergence of axons with linear trajectories strongly suggests that the axons originate from diverse input sources (see Discussion). 

To investigate the potential for toric spines to form synapses with novel inputs we investigated the CF of two complex spines, TS1 and TS7. Prior work defined "filling fraction" as the number of synapses formed on a postsynaptic target (in that case, a dendrite or neuron) divided by the total number of synapses available within a local volume of surrounding neuropil extending 2 m in all directions away from the postsynaptic target (Stepanyants et al., 2002). This is the typical range of filopodial extensions that occur across circuits and species in healthy tissue undergoing experience-dependent synaptic plasticity (Grutzendler et al., 2002; Trachtenberg et al., 2002). We define a closely related concept, CF, which is the number of axons connected to the spine divided by the total number of axons that could be connected by virtue of proximity and lack of intervening myelin (Mishchenko et al., 2010). To determine CF, all unmyelinated axons within two m of the target spine were reconstructed (Fig. 15A ). Infiltrating TS1 were 61 unmyelinated axons, 10 of which were synaptically connected to it (CF ϭ 0. 164). For TS7 the numbers were 64 axons, 10 connected, CF ϭ 0. 156. These low CFs are on par with prior reports of filling fractions for circuits in mammalian A, The number of axonal inputs strongly correlated with the number of active zones (linear regression: r 2 ϭ 0. 813, p Յ 0. 00001, 2 ϭ 49. 4) and the data were better fit by an exponential function ( 2 ϭ 14. 59; dotted line). This was also true for both subpopulations of toric spines on the dendrites (linear regression: r 2 ϭ 0. 895, p Յ 0. 00001, 2 ϭ 8. 88; exponential fit 2 ϭ 6. 177) and on the soma (r ϭ 0. 837, r 2 ϭ 0. 700, p ϭ 0. 00125, 2 ϭ 25. 46; exponential fit 2 ϭ 6. 65). There was a weak tendency for spines located on the soma to have more inputs (mean values 7. 1 vs 4. tissue, e. g., 0. 18 and 0. 19 for two densely reconstructed hippocampal dendrites (Mishchenko et al., 2010). Synaptic connectivity of the infiltrating axons could in principle be random or structured. To test this, we constructed a bootstrap analysis based on the assumption that each axon has equal ability to form an active zone with the target spine during a hypothesized epoch of microstructural plasticity. We determined the probability envelope for distributing the total number of AZs received (25 for TS1, 23 for TS7) across each potential axonal input (61 for TS1, 64 for TS7). The actual data fell well outside the probability envelope based on random access (Fig. 15C ). This indicates that connectivity is structured such that certain axons have, or develop, preferential access. Such non-random connectivity is also found in mammalian cortex (Kasthuri et al., 2015) and thalamus (Morgan et al., 2016), suggesting a common principle for circuit development. 

To open a path for linking microanatomical patterns of synaptic convergence with functional measurement of neu- ronal computation, we developed an ex vivo slice preparation to record from space-specific neurons. Recording electrodes were targeted to large somata in ICX under visual guidance (Fig. 16A ) and used to record membrane potential responses using a current-clamp protocol (Fig. 16B ). Neurons (n ϭ 20 from five owlets) exhibited a spec-Figure 13. Angular distances between the input trajectories onto toric spines. Histograms of the angular distance between each trajectory vector and every other vector for each spine, excluding the vector from the other end of same axon. Distances were calculated with the dot product (see Materials and Methods). Mean across all spines ϭ 89. 7°; SD ϭ 40; N ϭ 78 vectors, 976 pairings; bin width ϭ 5). A population of random trajectories predict a Gaussian distribution of distances with a mean of 90°, as observed from 6/7 spines and the overall population, while a population of bundled trajectories predicts a bimodal distribution with a mean of 90°, as observed for TS76. To assess potential connectivity of TS1 and TS7, all unmyelinated axons that passed within 2 m of the spine were reconstructed. A, SBEM section showing four annotated axons (shaded in color) synapsing onto TS1 (dark purple). Arrows point to nearby axons which do not make any synapses onto TS1. B, All 61 unmyelinated axons within 2 m of TS1 (solid gray). There were also five myelinated long-range axons that came within 2 m of TS1 (data not shown and excluded from this analysis; see text for explanation). In all, 10 unmyelinated axons synapsed onto TS1 (solid axons) and 51 did not (meshed axons). Thus, the CF for TS1 is 0. 164. The CF of TS7 (data not shown) is 0. 156 (10 connected axons, 54 unconnected axons). C, Bootstrap analysis of potential connectivity of TS1 (left panel) and TS7 (right panel). Histograms indicate the relative frequency of observing axons connected with the specified number of active zones to TS1 in the real data (solid circles) and 10,000 simulations (gray line; light gray envelope indicates Ϯ2 SD). In both spines, four of the 10 data points lie outside the envelope of expectation based on random connectivity. Insets display the same data over a magnified range. trum of intrinsic membrane properties (Fig. 16C, D ; Table 2 ). Some neurons (n ϭ 16) were filled with biocytin for retrospective imaging to determine cell type. In this population of neurons in juvenile owls, large atypical spines were prevalent including toric-like structures of narrow tubular build with multiple arms whose ends often infolded on themselves, but which infrequently presented as closed holes (Fig. 17D ). The morphologic diversity was greater than in the adult space-specific neuron population, with many juvenile cells expressing some mixture of atypical and typical spines (Fig. 17A, B ), while a minority (n ϭ 5) lacked a high density of spines (Fig. 14C ). Therefore, juvenile cells were not classified as Type I or Type II, but instead scored for the presence (n ϭ 7) or absence (n ϭ 9) of toric-like spines. There was a trend for juvenile cells with toric-like spines to have larger soma and thicker primary dendrites than those lacking them (Fig. 17E-G ), although these differences were not statistically significant, and both juvenile and adult cells had on average the same number of primary dendrites (mean ϭ 4. 5 for both), a comparative profile similar to the differences between Type I and Type II adult cells (Fig. 2 ). In total, the juvenile cells appeared as developmental precursors to toric space-specific neurons. Finally, in separate voltage-clamp experiments, electrical stimulation of the known input source to ICX, the ICCls, evoked typical AMPA-mediated EPSCs (Fig. 16E, F ; n ϭ 2), consistent with all prior studies of this monosynaptic connection.
section_name: Results

section_content: We describe the morphology and synaptic innervation patterns of a novel brain structure, toric spines. Across circuits and species, toric spine morphology appears unique. In contrast, while their high integrative capacity is unparalleled in degree, it is evocative of reports of complex spine types found in mammalian thalamus, hippocampus and cerebellum. This suggests that toric spines may represent a standout example of a processing motif common across the animal kingdom. The discussion below critically analyzes toric spine morphology, integrative capacity, implications for cellular computation, potential role in plasticity, and concludes by outlining testable hypotheses to probe more deeply into their structure and function. 

Both STED and SBEM imaging revealed a broad spectrum of spine morphologies. Features that varied mark-edly across the population include size, number of holes and degree of branching. Features that tend to be shared include narrow tubular build, paucity of bulbous heads, and wide necks at the process interface. This profile distinguishes toric spines from all previously documented spine types including typical thin, stubby, mushroom (Bourne and Harris, 2007; Chirillo et al., 2019) or branched spines (Bourne and Harris, 2008), and atypical dendritic structures including thorny excrescences, dendritic appendages or glomerular excrescences (Jones and Powell, 1969; Spacek and Lieberman, 1974; Rapisardi and Miles, 1984; Eris ¸ir et al., 1997; Morgan et al., 2016). Perhaps most striking, toric spines received up to 49 active zones contributed by up to 11 axonal inputs -a stark contrast to typical spines (e. g. Arellano et al., 2007) that compartmentalize single inputs or thorny excrescences which detonate their postsynaptic target (Table 3 ). Given their size and integrative capacity (discussed below), toric spines could be classified as atypical dendrites. We call them spines because, like typical spines, they are substantially smaller in caliber than the parent process, a size differential that exceeds the difference between higher and lower order dendritic branches. 

The narrow tubular build of toric spines should produce a short length constant due to high axial resistance, and therefore predicts that single synaptic inputs made onto one spine are isolated from those onto another (Rall, 1974). At the same time neighboring inputs onto the same spine should be subject to summation. Compartmental simulations of the reconstructed space-specific neuron to test this intuition have so far been unsuccessful due computational limitations of available modeling platforms. In comparison, simulations of neocortical pyramidal cells reconstructed from light microscopy (i. e., lacking ultrastructural detail and therefore amenable to compartmen- tal simulation) have shown that neuronal computation can be represented by a two-stage model in which individual dendritic branches perform the first-stage processing of summation of EPSPs (Poirazi et al., 2003). Active properties (voltage-gated ion channels) modify the rules of dendritic integration providing a range of operations from supra-linear to sub-linear. The second stage is linear summation of dendritic branch output at the soma. The multipolar structure of space-specific neurons is well suited to emulate this processing scheme with the significant addition of an earlier-stage operation at toric spines. 

Direct tests of these computational strategies require whole-cell recording in combination with targeted synaptic stimulation. We demonstrate, for the first time to our knowledge, the ability to make high-quality patch-clamp recordings from ICX neurons. Development of this ex vivo preparation opens a path for understanding the casual connections between microanatomical convergence and cellular computation. 

All large toric spines with analyzed convergence patterns (n ϭ 24) were found to receive synaptic input from multiple axons, indicating that they are structured to act as integrators. Synaptic convergence onto branched spines and other atypical dendritic protrusions has been suggested, though not firmly established, by earlier studies. Pioneering work using serial EM to reconstruct rat cerebellar Purkinje cells demonstrated a complex spine type with up to five branches. Independent axonal boutons could be found on distal branches of the same spine, several microns apart, raising the possibility that branched spines receive inputs from multiple axons, though this remains untested (Harris and Stevens, 1988). The prevalence of branched spines on Purkinje cells has been estimated between 2% and ϳ6% (4/64; Harris and Stevens, 1988; Lee et al., 2004 ) and ϳ15% on adult-generated granule cells in the dentate gyrus, where they are also polyinnervated (Bosch et al., 2015), suggesting an important role in neuronal computation. Synaptic glomeruli in mammalian thalamus receive convergent input from anatomically and neurochemically diverse sources (Jones and Powell, 1969; Spacek and Lieberman, 1974; Eris ¸ir et al., 1997). In a recent connectomics study, two densely reconstructed glomeruli received input from 15 different thalamocortical afferents (Morgan et al., 2016), a level of convergence similar to that observed in our dataset. The central postsynaptic compartment of glomeruli are dendritic shafts from which multiple tiny excrescences emerge, distinct electrotonically from toric spines. Nonetheless, the above results in combination with our findings raise the intriguing possibility that synaptic convergence onto complex spine types is more common than currently appreciated. 

The conclusion that toric spines are innervated by multiple input sources is supported by analysis of axonal convergence onto 27 spines and trajectory analysis for seven. Not one shared branchpoint was found (40 axons, mean length ϭ 20. 6 m). Until the axons are traced back to their parent cell bodies, a single source explanation cannot be logically eliminated. We think it is highly unlikely because it would require parent axons to branch, diverge at least 10-100s of micrometer, then re-converge and terminate within an extremely small volume of tissue. There is no evidence for such an arrangement. Owls are known for unusual axonal specializations, delay lines (Carr and Konishi, 1988 ), yet these are found in the projection from nucleus magnocellularis to nucleus laminaris (Carr and Konishi, 1990), distant brainstem structures, not ICX. In fact, ICCls axons, which provide the bulk of input to ICX, exhibit typical terminal fields in which axonal branches diverge, not converge (DeBello et al., 2001; Rodriguez-Contreras et al., 2005; McBride et al., 2008; McBride and DeBello, 2015). In principle, short branchlets with terminal boutons could be within range of a toric spine, but because those branchpoints are mere microns away from the parent axon, such arbors could not have contributed to the divergent trajectories of much longer axonal segments observed in our dataset. 

The fact that individual ICCls axons cannot account for the convergence results does not rule out ICCls axons as a potential input source to toric spines. In fact, spacespecific neurons are known to receive input from multiple ICCls neurons encoding different sound frequencies. Integration across frequency is required to eliminate phase ambiguity and ultimately provide the owl with an accurate estimation of sound source azimuth. At the level of spiking output, frequency convergence is a non-linear operation. Thus, one hypothesis is that axons originating from multiple ICCls neurons encoding the same ITD but different frequencies converge onto individual toric spines where their EPSPs are combined and thresholded (Peña and Konishi, 2000, 2002). These predictions are testable using dual anterograde injections placed at different dorsoventral laminae, and whole-cell recordings, respectively. 

A second possibility is that toric spines integrate ITD and ILD, a multiplicative operation. Binaural integration first occurs in the ICCls and is largely inherited by spacespecific neurons in ICX. Yet the precise anatomical and functional boundaries between ICCls and ICX are unclear. Several cytoarchitectonic markers delineate the lateral border of ICX, and separately, the border of the core of the central nucleus. But on close inspection of commonly used markers calbindin, calretinin and CaMKII, the ICCls-ICX border appears as a gradient. Microelectrode surveys also support the idea of functional heterogeneity along this border. Moreover, while cytoarchitectonic markers were not routinely employed in our study (see Materials and Methods), estimates based on comparison to similar sections from other owls indicate that at least some Type I space-specific neurons, including the one in the SBEM volume, were located in ICCls. This is a minor revision to the notion that IC-OT projection neurons are strictly confined to ICX. More importantly, it motivates the second hypothesis, that toric spines in ICCls receive convergent input encoding behaviorally relevant ITD-ILD pairs, from contralateral ICC (for ITD) and contralateral dorsal lateral lemniscal nucleus pars posterior (for ILD), and perform a multiplicative operation (Peña and Konishi, 2001, 2004 ) via a different complement of ion channels. A third possibility is that both hypotheses are true: toric spines are multipurpose integrators whose computation can be customized to cellular need. 

Other potential sources of input to toric spines include the arcopallial gaze field (Cohen et al., 1998), which provides top-down attentional modulation and the intermediate/dOT, which provides a visually-based instructive signal to guide calibration of the auditory map (Hyde and Knudsen, 2000; Luksch et al., 2000). Both are anatomically sparse and unlikely to be dominant contributors. space-specific neurons do receive lateral connections from GABAergic neurons within ICX. The inputs to toric spines we observed, however, more closely resemble asymmetric glutamatergic synapses. 

To our knowledge, this is the first report of a toric postsynaptic structure. Presynaptic perforations described extensively in the literature are cytoplasmically contiguous (Sorra et al., 1998). In contrast, at least some tori in our SBEM volume were interrupted by plasma membranes that appeared as sites of contact between two "arms" of the same spine (Fig. 9A, B ). The simplest interpretation is these represent enduring adhesions between filopodial extensions emanating from the same base structure. These could be maintained or later converted to cytoplasmic continuities via anastomosis. Regardless of the developmental mechanics, the prevalence and functional properties of contact sites versus cytoplasmic continuities is expected to significantly impact current flow during synaptic stimulation. 

The morphologies of ICX neurons from 35-to 40-d old juveniles labeled ex vivo and adult space-specific neurons labeled in vivo were qualitatively similar with a few clear differences: (1) juvenile cells often expressed a mix of atypical (branched and toric-like) and typical spines whereas adult cells were dominated by either toric (Type I) or typical (Type II) spines, but never both; (2) juvenile cells appeared smaller; and (3) juvenile cells sometimes appeared ciliated. We consider three interpretations: ex vivo artifact, unrelated cell type or developmental age. 

Preparation of acute brain slices can induce changes in dendritic morphology in comparison to the in vivo condition as reflected in perfusion-fixed tissue (Fiala et al., 2003). In this scenario, toric-like spines would represent the sliceinduced disassembly of fully-formed toric spines. We think this is highly unlikely due to the number of neuropil elements involved in such a disassembly, specifically the high density of intertwining axons. Whole cell recordings in juvenile slices were targeted under visual guidance to the largest somatic profiles in ICX, the exact location of retrogradely labeled toric space-specific neurons in adults. That the measured size of juvenile cells was smaller is expected from the differential impact of immersion versus perfusion fixation. Specifically, measurements of three neuron types in primate hippocampus indicate that apparent soma volume is 1. 36 -2. 7 times larger in perfusion-fixed than immersion-fixed tissue (Lavenex et al., 2009). On average, our adult cells (perfused brain) were 2. 42 (volume) and 1. 81 (surface area) times as large as our juvenile cells (immersed slices). Further, the one morphologic attribute not subject to differential shrinkage, number of primary dendrites, was the same in juveniles (4. 5) and adults (4. 7). In summary, the most straightforward interpretation is that the observed juvenile cell types represent developmental precursors of adult cell types. 

Lesion reconstruction (Brainard and Knudsen, 1993), anatomic (Feldman and Knudsen, 1997; DeBello et al., 2001), and pharmacological (Feldman et al., 1996; Zheng and Knudsen, 1999) experiments pinpoint space-specific neurons as the cellular locus of learning (Knudsen, 2002). This raises the question of whether toric spines have a role in experience-dependent synaptic plasticity (DeBello et al., 2014). We found that toric spines are conspicuously well structured to be hubs for microstructural change: of all unmyelinated axons passing within 2 m of TS1 and TS7, only ϳ16% were synaptically connected. Consequently, the number of different wiring diagrams that could be achieved by biologically plausible, small-scale dendritic extensions is very large (Stepanyants et al., 2002; Stepanyants and Chklovskii, 2005). We propose a model in which toric spines actively search for axonal inputs in their local environment, make "tester" synapses, then stabilize or eliminate these on the basis of efficacy in firing the postsynaptic cell (DeBello, 2008). This would provide a mechanism for pattern detection within a complex array of input sources, and the means to adapt connectivity to changing experiential needs. The high degree of cytoskeletal dynamics implicit in this model could also explain the broad spectrum of observed morphologies, with each spine stabilized not by genetic programs but by chance encounters of its extensions within a local, yet information diverse, axonal milieu.
section_name: Discussion

section_content: This work was supported by National Institutes of Health Grants R01DC05640 (from NIDCD ; to W. M. D. ), BRAIN-STIM grant from the University of California Davis Office of Research and Behavioral Health Center for Excellence (W. M. D. ), R01NS104911 (from the BRAIN Initiative ; to W. M. D. and M. H. E. ), P41GM103412 (for support of National Center for Molecular Imaging Research ; to M. H. E. ), 2018 NARSAD Young Investigator grant (E. G. A. ), NSF1754831 and MH114178 (to D. F. ), and a UCDMC Interdepartmental Grant (D. F. and W. M. D.