Chemical connections were detected manually, curve fit, and then characterized for strength, kinetics, and short-term plasticity in a multi-stage process. Manual connectivity calls were used to train a machine classifier (see Seeman, Campagnola, et al. 2018 for more detail). After training, the machine classifier revealed a set of potential false positives or false negatives that were manually re-checked. Electrical connections were also detected and characterized.
Chemical connections were visually identified from average postsynaptic responses sorted by recording mode (voltage and current clamp) and membrane potential (see QC criteria for holding potential ranges for excitatory and inhibitory connections). Individual responses were time aligned to the presynaptic spike. Responses that failed QC were not included in the average.

The User Latency value was used to initialize automated curve fitting of the postsynaptic response. Parameters from the curve fits were shown for each quadrant. When the user was satisfied with the curve fitting output, the goodness of fit was manually passed or failed. The user would take into consideration gap junctions, electrical artifacts, the shape of the fit or other factors when deciding whether the curve fit accurately reflected the response. Notes were also made by the user and used to update, and subsequently test, new fitting algorithms.

Once a chemical connection was identified, the strength and kinetics of that connection, as well as the dynamics (short-term plasticity (STP) and variability), were characterized utilizing the curve fits generated during connection detection.
The "strength" of a connection is an important characteristic particularly when we want to start comparing connections across cell class. However, we know that strength changes dynamically over time (as is highlighted in our analysis of STP) as well as stochastically from spike to spike. Given that our stimulus set utilizes trains of various frequencies we want to ensure that our metric of strength would be useful in comparisons across connections and not be stimulus dependent. Our resting-state strength is a metric that works well for most connections; however, some facilitating connections have a very small resting-state strength which may be misleading. Rise and decay kinetics are more faithfully preserved even as the amplitude of the connection changes and thus, was more straightforward to characterize.
Resting-state amplitude was determined from responses in which the presynaptic spike follows a period of quiescence. Usually this is the first pulse in each stimulus train. Individual responses were averaged and curve fit; the amplitude output of the fit served as our metric of strength.

Latency, rise and decay kinetics were also determined from curve fits to average synaptic responses. Kinetics were calculated for postsynaptic currents (PSCs) and potentials (PSPs) held at either -70 or -55 mV depending on whether the synapse was excitatory or inhibitory (see QC requirement).
Short-term plasticity was analyzed from data recorded in current clamp utilizing different stimuli depending on the analysis. We measured three main metrics of short-term plasticity: paired-pulse ratio (PPR), train-induced STP, and recovery from train-induced STP.
Short-term plasticity is often calculated as a ratio of response amplitudes for each pulse in the stimulus train. Many connections in our dataset have small responses that are close to the noise level in our recordings making the use of an amplitude ratio unstable. This situation is accentuated for connections that are strongly depressing or facilitating resulting in spurious ratio measurements. To avoid this instability, we use the 90th percentile PSP amplitude as an approximation of the "maximum" amplitude of a synapse and normalize our STP measurements by this value. Below are the calculations for each STP metric.
Paired-pulse STP was measured from the 50Hz stimulus as:

Train-induced STP was measured from the 50 Hz stimulus as:

Recovery from STP was measured from all stimulus frequencies with a 250 ms delay between the last induction pulse (8th) and first recovery pulse (9th):

The amplitude of PSPs for a synaptic connection vary randomly each time neurotransmitter is released. This is often reported using the coefficient of variation (CV). In the mouse cortex, however, typical PSP amplitudes can be much smaller than the background electrical noise in the cell. CV in this regime is thus dominated by noise and tells us little about the physiology of the synaptic connections. In our dataset, PSP amplitude variability is reported using a metric that is adjusted (aCV) to correct for the effect of background noise:

This metric has a value of 1.0 when the standard deviation of the PSP amplitudes (after noise correction) is the same as the median amplitude. The noise correction itself introduces a new source of variance, however, which can sometimes lead to this value being negative. We measure variability in the resting state for each connection as well as in various states of induced short-term plasticity.

We developed a new model of synaptic vesicle release that provides a more comprehensive description of each synaptic connection and also allows to predict the behavior of the synapse in response to arbitrary stimuli. The model includes basic quantal release parameters (release probability, number of release sites, and quantal size) as well as short term plasticity (vesicle depletion, depression, and facilitation with varying recovery time constants) and accounts for recording noise.
Our dataset includes best fit parameters (using a maximum likelihood estimation) of the model for many connections that can be used to simulate different types of connections, or as a basis for comparison between connection types.
Electrical connections, formed by gap junctions, were also detected in our dataset. We used the long-pulse stimuli to characterize the strength of electrical connections.

Learn about Synaptic Physiology Analysis Methods Quality Control with comprehensive guides and examples from Allen Institute for Brain Science.
Ensuring the quality of the Synaptic Physiology Dataset is of critical importance. We applied quality control criteria at several stages of data processing and analysis. During the processing stages, in order for data to proceed into successive stages it had to pass quality control for the prior stage. During analysis of connection properties, each analysis (e.g., strength, kinetics, short-term plasticity) had independent quality control metrics. Below is a workflow of the quality control metrics required for inclusion at each stage of processing and analysis.
Data was first processed to ensure that it was of good enough quality to detect a connection

During characterization of a synapse, different QC criteria were applied based on the metric that was being analyzed.

Learn about Synaptic Physiology Analysis Methods with comprehensive guides and examples from Allen Institute for Brain Science.
Data generated during the experiment were analyzed in stages to identify and characterize connections between distinct cell subclasses. Below is a workflow describing how this dataset was processed and analyzed; more detail can be found in Seeman, Campagnola et al. 2018 and Campagnola, Seeman et al. 2022.
Transgenic cell subclass was identified from the overlap of the fluorescent reporter of the cell and that of the recording pipette. Spiny and aspiny cells were defined from morphological analysis. Cells were annotated with a target layer during the experiment, and later a corrected layer determined from biocytin and DAPI staining.

Connections were identified from manual inspection of the postsynaptic response aligned to the presynaptic spike. Excitatory and inhibitory connections could be observed in both voltage and current clamp at a holding potential that increased the driving force for each connection class. Weak connections were detectable by averaging hundreds of postsynaptic responses.

Resting-state synaptic strength, rise/decay kinetics, and depth of short-term plasticity were calculated from curve fits of postsynaptic responses.
More on connection characterization

Data passed through multiple quality control filters as part of the analysis workflow.

Learn about Synaptic Physiology Experimental Methods Experimental Protocol with comprehensive guides and examples from Allen Institute for Brain Science.
During a single experiment, a stereotyped set of stimuli were delivered to recorded cells to characterize the strength, kinetics, and short-term plasticity of synaptic connections. Each stimulus was delivered to cells in turn while recording the response from all other cells. Each stimulus train consisted of eight pulses to induce short-term plasticity (induction) followed by a delay and four pulses to measure recovery from dynamic effects. The frequency of pulses and delay were varied to measure these effects on synaptic connections. There was a 15 second delay between each stimulus to allow the cells to fully recover from stimulation. A separate set of stimuli were delivered to characterize the intrinsic properties of cells.
A variety of stimulus frequencies were used, each with a 250 ms delay between the eight induction pulses and four recovery pulses.

For the 50 Hz stimulus, a range of recovery delays were interposed between the eight induction pulses and four recovery pulses to evaluate the time course of recovery from short-term effects.

The "mixed frequency" stimulus was composed of eight pulses delivered at 30Hz immediately followed by 30 pulses whose intervals were a random resequencing of 29 exponentially increasing intervals between 5 and 100ms. The intervals were fixed across sweeps and experiments. This stimulus allowed us to explore a wider range of stimulus frequencies that could be used to inform our model.
Each stimulus, one frequency with one delay or mixed frequency, was repeated at least five times with a 15 second rest period between repeated trials. Spikes were recorded in the presynaptic cell along with postsynaptic responses.

Long, sustained stimuli were delivered to each cell to measure intrinsic properties as well as quantify electrical synapses in current clamp. Intrinsic features were extracted from these stimuli using the IPFX package.
Subthreshold, mostly hyperpolarizing, stimuli were used to measure properties such as input resistance and sag. This stimulus set was initiated with a pulse at -20 pA while keeping the neuron at -70 mV. The voltage response to each current step was measured online and successive current steps were titrated to target response voltages of -68, -72, -75, -80, and -85 mV so as to reliably activate Ih when present.

Suprathreshold, depolarizing stimuli were delivered to measure spiking properties. These stimuli started at rheobase and increased 25 pA for 6 intervals.

A 15-second sinusoidal chirp that increased in frequency from 0.2 to 40 Hz. The amplitude was targeted to evoke a response magnitude that measured ~10 mV from peak to trough.
