Anatomy and activity patterns in a multifunctional motor neuron and its surrounding circuits

  1. Mária Ashaber
  2. Yusuke Tomina
  3. Pegah Kassraian
  4. Eric A Bushong
  5. William B Kristan
  6. Mark H Ellisman
  7. Daniel A Wagenaar  Is a corresponding author
  1. Division of Biology and Biological Engineering, California Institute of Technology, United States
  2. Division of Biological Sciences, University of California, San Diego, United States
  3. National Center for Microscopy and Imaging Research, University of California, San Diego, United States
  4. Department of Neurosciences, UCSD School of Medicine, United States
14 figures, 1 video, 7 tables and 1 additional file

Figures

Approach.

Several fictive behaviors were induced in the isolated nervous system of a medicinal leech while one segmental ganglion was imaged using a voltage-sensitive dye. After fixation and resin embedding, the ganglion was X-ray-imaged to verify that the geometry of somata was preserved. Finally, the neuropil was imaged at nanometer resolution with serial blockface electron microscopy and selected neurons were manually traced. Scale bars: 1 cm (leech photograph); 100 μm (all others).

Fictive swimming imaged using VSD.

(a) Images of the dorsal (left) and ventral (right) aspects of a leech ganglion simultaneously obtained using a double-sided microscope. ‘R’ indicates the right side of the ganglion (i.e., the animal’s right when dorsal side up). (b) Selected VSD traces during fictive swimming. From the dorsal surface: dorsal and ventral inhibitory and excitatory motor neurons DI-1, VI-2, DE-3, and VE-4; from the ventral surface: the Retzius cells (neuromodulatory interneurons) and cell 153L (an interneuron). All those cells are known to be rhythmically active during swimming. Also shown is CVL, an excitor of ventrolateral circular muscles that was intracellularly recorded during the trial as a control to verify that fluorescence signals reflect membrane potential changes as they should. Scale bars: 0.2% relative fluorescence change, 5 mV membrane potential change. (c) Magnitude (radial axis from 0 to 1) and phase (angular coordinate) of the coherence of activity in individual neurons with the swim rhythm in motor neuron DE-3R. Error bars indicate confidence intervals based on a multi-taper estimate. (d) Coherence maps of the VSD signals of all cells on the dorsal (left) and ventral (right) surfaces of the ganglion. Colors of cell bodies indicate coherence relative to DE-3R. The yellow needle on CVL indicates a sharp electrode for intracellular recording. Color scale applies to all panels.

Fictive crawling imaged using VSD.

(a) Selected VSD traces during fictive crawling. From the dorsal surface: dorsal and ventral inhibitory and excitatory motor neurons DI-1, VI-2, DE-3, and VE-4; from the ventral surface: the Retzius cells and CV cells. All those cells are known to be rhythmically active during crawling. Below the VSD traces, a simultaneously recorded intracellular trace of the annulus erector motor neuron AER is displayed (in both columns). Scale bars: 0.2%, 10 mV. Gray lines mark hyperpolarized phase of AER. (b) Magnitude and phase of the coherence of activity in individual neurons with the crawl rhythm in motor neuron DE-3R. (c) Coherence maps of the VSD signals of all cells on the dorsal (left) and ventral (right) surfaces of the ganglion. Colors of cell bodies indicate coherence relative to DE-3R. The yellow needle on AER indicates a sharp electrode for intracellular recording. Color scale applies to all panels.

Fictive local bending imaged using VSD.

(a) Selected VSD traces during fictive local bending. From the dorsal surface: dorsal and ventral inhibitory and excitatory motor neurons DI-1, VI-2, DE-3, and VE-4; from the ventral surface: the Retzius cells, ‘Anterior Pagoda’ cells (‘AP’; well-known postsynaptic partners of the P cells with unknown function). Below the traces, a simultaneously recorded intracellular trace of the PVL cell is displayed with injected current trains (in both columns). Scale bars: 0.2% relative fluorescence change, 100 mV membrane potential change, 2 nA injected current. (b) Magnitude and phase of the coherence of activity in individual neurons with the local bend rhythm in DE-3R. (c) Coherence maps of the VSD signals of all cells on the dorsal (left) and ventral (right) surfaces of the ganglion. Colors of cell bodies indicate coherence relative to DE-3R. A yellow needle on PVL indicates a sharp electrode for electrical stimulation. Note that PVL was only weakly stained by the VSD because it was left partially covered with sheath tissue to preserve its health.

Mapping between functional and anatomical images.

(a) Light micrograph (left) and X-ray image (right) of the ganglion in which we recorded neuronal activity using a VSD. Matching colors label the same cells on the two images. A: Anterior, P: Posterior, L: Left, R: Right. (b) Transverse section from X-ray tomographic image stack of the ganglion (top). Only the shaded area was imaged with SBEM (bottom). Arrowheads mark several somata that are (white) or are not (red) included in the SBEM volume. Circles mark neurites that facilitated complete mapping between the two imaging modalities. D: Dorsal, V: Ventral. Note that the X-ray image in (b) is shown in reverse video relative to the x-ray image in (a) for easier visual comparison with the SBEM image.

Figure 6 with 3 supplements
Electron microscopic tracing: neurites and synapses of motor neuron DE-3R.

(a) The principal neurite of DE-3R near its entrance to the neuropil (dashed yellow outline). (b) Two branches of the neurite of DE-3R (dashed outlines). (c) A synaptic connection onto DE-3R from an inhibitory motor neuron (DI-1R). Arrowheads: synapses, Pre: presynaptic terminal, v: vesicles. (d) A synapse onto DE-3R from an interneuron (cell 24 on the canonical map [Wagenaar, 2017]).

Figure 6—figure supplement 1
Comparison of SEM with TEM for interpreting synapses.

(a) Our interpretation of a small section of our SBEM image: a process of cell DE-3 (yellow) and a presynaptic partner (purple). (b) Same area without overlay. Arrowhead: synapse. Stars: Mitochondria. (c) A similar area imaged with TEM (from a ganglion from another leech). (d) Same area as (c), computationally blurred (radius: 15 μm) to simulate the lower resolving power of SEM.

Figure 6—figure supplement 2
Basis for estimating true resolution of SEM images: Spectral power in the images.

Average based on four representive 2048 × 2048 pixel regions from the volume.

Figure 6—figure supplement 3
Workflow for identifying neurons in our data with the canonical map of the ganglion.
Traced neuronal arborizations.

(a) Fully reconstructed arborization of DE-3R overlaid on a dorsal projection and a transverse section of the micro-CT data. (b) A selection of presynaptic partners with previously known identities. (c) A selection of presynaptic partners not matched to neurons previously described in the literature. (Top row: cell bodies on dorsal aspect; bottom row: cell bodies on ventral aspect of the ganglion.) Other cells in this category are 20R, 25R, 154R, 156L. Numbers in parentheses are number of synapses between each cell and DE-3R.

Full tracing of DE-3R (thick black line, soma location marked ‘DE-3’) and backtracings of all synaptic partners.

Partners that we could identify with cells in the VSD recording are marked with (arbitrary) colors. Small gray disks indicate partner neurons that could not be cross-identified between EM and VSD image stacks.

Correlation between anatomical and functional measures of synaptic strength.

(a) Scatter plot of synapse count vs coherence magnitude during swimming for presynaptic partners with at least two synapses, with linear regression line (data from swim trial #1, R = 0.26, p = 0.20, n.s.). (b) Same for crawling (data from crawl trial #1, R = 0.38, p = 0.05, n.s.). (c) Correlation coefficients for all trials (left) and when raw synapse count was replaced by proximity weight (right; see Materials and methods).

Figure 10 with 1 supplement
Spatial distribution of synapses onto DE-3.

(a) Distribution of excitatory and inhibitory synapses. (b) Histogram of the length of the paths between those synapses and the soma. (c) Distribution of synapses more strongly associated with a certain behavior. (d) Histogram of the length of the paths between those synapses and the soma (LB: Local bend).

Figure 10—figure supplement 1
Path lengths between synapses and the trunk of DE-3R.

(a) Path length for excitatory vs inhibitory synapses. (b) Path lengths for synapses from neurons more strongly associated with a certain behavior.

Figure 11 with 3 supplements
Synaptic clustering.

(a) Tracing of DE-3R with synapses (arbitrarily) colored by presynaptic partner. (b) Clusters (elliptic areas) associated with synchronized synapses for a swimming trial (top) and a crawling trial (bottom). Within the same connectome, clusters of synchronized synapses differ with respect to their spatial extent for the two behaviors: During swimming, synchronization extends over larger areas along the neurite than during crawling. Shown are synaptic clusters obtained with parameter values (dNN, dext) = (5 μm, 65 μm) for swimming and (dNN, dext) = (7.5 μm, 15 μm) for crawling, respectively. Synapses are colored by the coherence between the activity of their presynaptic partner and DE-3R during the behaviors (as in Figure 2) and clusters are colored by the average coherence of their constituent presynaptic partners. Inset: Explanation of clustering parameters (see Materials and methods). (c) All clustering results for swim and crawl trials. Color indicates the degree of correspondence between spatial clusters and functional grouping expressed as an F-ratio from complex ANOVA (see Materials and methods) as a function of clustering parameters.

Figure 11—source data 1

Results of the ANOVA analysis of synaptic clusters in all trials.

Included in the table are all parameter values that resulted in a significant link between anatomy and activity (p < 0.05). Parameter values are expressed in μm. Only the 10 most significant parameter value pairs are included for each trial.

https://cdn.elifesciences.org/articles/61881/elife-61881-fig11-data1-v2.csv
Figure 11—figure supplement 1
Clustering results for the local bend trials.

Out of the four trials, the first three involved stimulation of the left PV cell; the final trial involved stimulation of the right PV cell. In contrast to the swim and crawl trials (Figure 10b), the clustering results for local bend trials are inconsistent.

Figure 11—figure supplement 2
Peak F-ratios for connecting spatial clusters to functional activity.

These graphs represent slices through the data shown in Figure 11c, vertically through the highest peak in the parameter space. Lines are Gaussians fitted through the data; see Materials and methods.

Figure 11—figure supplement 3
Demonstration of the F-ratio method using synthetic data.

(a) Synthetic data where the coherence of a neuron’s activity (location of each dot) is unrelated to its membership in a spatially defined cluster (color), hence F = 1. (b) Synthetic data where such a relationship does exist, hence F > 1. See Materials and methods.

Wiring diagram of DE3.

The diagram shows synaptic connections onto DE3, including newly found presynaptic neurons in this study (yellow) and previously known ones (green if also found in this study; white otherwise). Only confidently identified neurons shown in Figure 7 and Table 1 are represented as ‘found.’ Other presynaptic partners found but not positively identified in this study are not shown here. This diagram does not depict whether connections are ipsi- or contralateral.

Author response image 1
Author response image 2

Videos

Video 1
Full tracing of motor neuron DE-3R with all of its input synapses, visualized growing out from the soma to the distal branches (0:00 to 0:05).

Back tracings of all DE-3R’s presynaptic partners to their somata (0:05 to 0:12). Rotational view of the completed tracing (0:12 to 0:15).

Tables

Table 1
Identified partner neurons of DE-3R.

These synaptic partners could be confidently assigned as previously described neurons.

CellSynapse countKnown function
DI-1L4Inhibitor of dorsomedial longitudinal muscles
DI-1R52Inhibitor of dorsomedial longitudinal muscles
VI-2R2Inhibitor of ventral longitudinal muscles
8R13Excitor of ventral longitudinal muscles
101R2Inhibitor of dorsoventral muscles
102R6Inhibitor of dorsal longitudinal muscles
107R12Excitor of dorsomedial longitudinal muscles
109R7Excitor of lateral dorsoventral muscles
117R16Excitor of medial dorsoventral muscles
LR3Excitor of dorsal and ventral longitudinal muscles
Table 2
Frequency of clusters of different sizes for clustering parameters (dNN, dext) = (5 µm, 65 µm), the parameters that gave the largest F-ratio for the swim trials.

Synapse count: Number of synapses in a cluster. Frequency: Number of clusters with the respective synapse count. Number of presynaptic partners: Number of unique presynaptic neurons contributing to the synapses in the respective clusters.

Synapse countFrequencyNumber of presynaptic partners
291 (6x), 2 (3x)
441 (1x), 2 (3x)
511 (1x)
622 (2x)
712 (1x)
822 (1x), 3 (1x)
933 (1x), 4 (2x)
1113 (1x)
1513 (1x)
2316 (1x)
Table 3
Frequency of clusters of different sizes for clustering parameters (dNN, dext) = (7.5 µm, 15 µm), the parameters that gave the largest F-ratio for the crawl trials.

Columns as in Table 2.

Synapse countFrequencyNumber of presynaptic partners
2281 (23x), 2 (5x)
3111 (7x), 2 (4x)
491 (4x), 2 (5x)
521 (2x)
622 (2x)
731 (1x), 2 (1x), 3 (1x)
1213 (1x)
Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Hirudo verbana)Wild-type backgroundNiagara Leeches
Other (voltage- sensitive dye)VF2.1(OMe).HWoodford et al., 2015Courtesy Miller lab
Software, algorithmSBEMAlignThis paperhttps://github.com/wagenadl/sbemalign; Wagenaar, 2021a;copy archived at swh:1:rev:d76dcc55e7dad3e7bca91de24d20d201696a5339
Software, algorithmSBEMViewerThis paperhttps://github.com/wagenadl/sbemviewer; Wagenaar, 2021b;copy archived at swh:1:rev:8f8d3d2bcae39e165993d9e11ffe173640b940db
Software, algorithmGVoxThis paperhttps://github.com/wagenadl/gvox; Wagenaar, 2021c;copy archived at swh:1:rev:5e7ccd2273caed49bac3e09ca39de68a0b182fc5
Table 4
List of raw data trials and how they are referred to in the paper.
FigureBehaviorTrial no. in paperTrial no. in raw data
Figure 2Swim16
Swim28
Figure 3Crawl115
Crawl217
Figure 4Local bend19
Local bend210
Local bend311
Local bend412
Author response table 1
Correlation between synapse count and coherence magnitude.
TrialPearsonp-valueSpearmanp-value
Swim 10.1900.3420.0570.78
Swim 20.1390.4890.1130.576
Crawl 10.1730.387-0.0450.823
Crawl 20.1540.4430.1390.49
Local Bend 10.1520.450-0.0140.944
Local Bend 20.2020.3110.1170.562
Local Bend 3-0.2790.158-0.1930.336
Local Bend 40.0720.721-0.0470.816
Author response table 2
Correlation between weighted synapse count and coherence magnitude.
TrialPearsonp-valueSpearmanp-value
Swim 10.2510.2060.1600.425
Swim 20.1580.4310.1180.558
Crawl 10.1340.504-0.1530.447
Crawl 20.0670.740-0.0680.737
Local Bend 10.1490.459-0.0280.891
Local Bend 20.2680.1770.1940.333
Local Bend 3-0.2400.227-0.0890.661
Local Bend 40.1130.5750.0940.641

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Mária Ashaber
  2. Yusuke Tomina
  3. Pegah Kassraian
  4. Eric A Bushong
  5. William B Kristan
  6. Mark H Ellisman
  7. Daniel A Wagenaar
(2021)
Anatomy and activity patterns in a multifunctional motor neuron and its surrounding circuits
eLife 10:e61881.
https://doi.org/10.7554/eLife.61881