Synergistic plasticity of intrinsic conductance and electrical coupling restores synchrony in an intact motor network

Motor neurons of the crustacean cardiac ganglion generate virtually identical, synchronized output despite the fact that each neuron uses distinct conductance magnitudes. As a result of this variability, manipulations that target ionic conductances have distinct effects on neurons within the same ganglion, disrupting synchronized motor neuron output that is necessary for proper cardiac function. We hypothesized that robustness in network output is accomplished via plasticity that counters such destabilizing influences. By blocking high-threshold K+ conductances in motor neurons within the ongoing cardiac network, we discovered that compensation both resynchronized the network and helped restore excitability. Using model findings to guide experimentation, we determined that compensatory increases of both GA and electrical coupling restored function in the network. This is one of the first direct demonstrations of the physiological regulation of coupling conductance in a compensatory context, and of synergistic plasticity across cell- and network-level mechanisms in the restoration of output. DOI: http://dx.doi.org/10.7554/eLife.16879.001


Introduction
The hallmarks of robust central pattern generator (CPG) output are appropriately tuned excitability of individual neurons combined with circuit-level interactions that maintain appropriate temporal coordination (i.e., phasing) of these neurons. Through both developmental and ongoing tuning processes, CPGs can maintain reliable network output for decades across the lifespan of an individual, despite constant feedback from the changing nature of both the organismal and natural environment. Yet underlying this constant reliability of network output exists a surprising amount of variability in the individual parameters necessary for producing activity. For instance, despite having nearly identical output across animals, networks can exhibit a five-fold or more range in intrinsic and synaptic conductance values Schulz et al., 2006;Marder, 2011;Roffman et al., 2012). Such variability in intrinsic conductances is not limited to CPGs, but has been documented in several cell types of the mammalian brain, including cerebellar Purkinje cells (Swensen and Bean, 2005), and globus pallidus neurons (Günay et al., 2008). Additionally, the synaptic strengths between mammalian central neurons have been shown to vary in several brain regions (Nelson and Turrigiano, 2008;Turrigiano, 2008;Maffei et al., 2012).
The origins and implications of this variability are still an intense area of investigation (Krubitzer and Kahn, 2003;Turrigiano and Nelson, 2004;Ciarleglio et al., 2015). We hypothesize that this variability might be a result of ongoing compensatory changes required to maintain reliable output over time. This compensation, termed homeostatic plasticity, has been well documented both for plasticity of intrinsic excitability via changing ionic conductances (Turrigiano et al., 1994), as well as for changes in chemical synaptic strength (Desai, 2004;Turrigiano, 2012). Variability in conductances may also be an adaptive trait in and of itself: variable solutions that produce convergent circuit output may provide a selective advantage, or perhaps be a substrate for adaptation and evolution Grashow et al., 2009).
Regardless of whether such variability is the result of homeostatic compensation, differential tuning across networks, or a combination of these and other heretofore undiscovered causes, a potential cost to such variability has recently been identified. In the cardiac ganglion (CG) of the Jonah crab (Cancer borealis), five Large Cell motor neurons (LCs) generate completely synchronous output, as a result of pacemaker inputs within the network, to drive simultaneous heart muscle contraction in the crab (Tazaki, 1972). Despite completely uniform and synchronous activity within the network, LCs show highly variable underlying maximal conductances (Ransdell et al., 2012). These variable conductances render the neurons susceptible to perturbations that target a subset of ionic conductances: when high-threshold K + currents were blocked with tetraethylammonium (TEA), the motor neurons lost coordinated output and became divergent in their patterns of firing (Ransdell et al., 2013a). These CG neurons compensate for this change in excitability, presumably to homeostatically maintain a target level of excitability (Ransdell et al., 2012). However, none of the network level impacts of this perturbation and plasticity have been investigated. Indeed, it is difficult to study homeostatic plasticity in intact networks and to simultaneously take into account both properties of individual cells as well as their network interactions. In the present study, we discovered that LC variability makes the network vulnerable to desynchronization as a result of TEA exposure, but that compensation resynchronizes the network within 30-60 min via both intrinsic cellular and circuit-level physiological mechanisms. To examine the underlying mechanisms, we developed a biophysical computational model of the entire cardiac network. The network model enabled a comprehensive search of the conductance space for potential compensatory mechanisms that preserved network synchrony, and we used these findings to guide further experimentation. Our study revealed cooperative homeostatic plasticity among intrinsic conductances and electrical coupling across multiple cells in the cardiac network. We interpret this as a novel homeostatic compensatory mechanism contributing to the overall robustness of CPG output. eLife digest Neurons can communicate with each other by releasing chemicals called neurotransmitters, or by forming direct connections with each other known as gap junctions. These direct connections allow electrical impulses to flow from one neuron to another via pores in the membranes between the cells. Unlike communication via neurotransmitters, gap junctions are usually thought to be hard-wired and unchanging over the life of the animal.
Lane et al. recorded electrical activity in a network of neurons that generates rhythmic heart contractions in the Jonah crab. Neurons in this network usually all fire an electrical impulse at the same time, which is crucial to make sure that the whole heart contracts at the same time. The experiments show that drugs that block potassium channel pores in the membrane cause the neurons to fire too much and at different times to each other.
However, the network of neurons soon adapted to the changes caused by the drugs and returned to working as normal. Mimicking these changes in a computer model of the neuron network, together with experimental data, showed that changes to the gap junctions play a major role in restoring normal activity to the network. The next step following on from this research is to understand how a network of neurons 'senses' that it is not working normally and changes its electrical activity.

Results
Exposure to TEA desynchronizes LC burst waveforms and increases excitability Ransdell et al. (2013a) were able to repeatably reduce the magnitude of high-threshold K + currents (I Kd + I BKKCA ) by~92% in isolated LCs with 25 mM TEA. We used this experimental manipulation on the 3 anterior LCs in intact CGs ( Figure 1A) causing LCs in the same ganglion to change from identical ( Figure 1B) to divergent, asynchronous output ( Figure 2A,I,II) after exposure to TEA. After application of TEA, motor neurons became noticeably more depolarized during burst potentials and LC spiking seen on the extracellular recordings increased substantially ( Figure 2A). Additionally, comparison of intracellular voltage waveforms revealed a loss of conserved output (Figure 2A). These results are consistent with the hypothesis that variable underlying conductances of the LCs makes them vulnerable to a uniform perturbation of a subset of conductances such as the TEA blockade that targets high-threshold K + conductances. Because our previous results demonstrated that the change in excitability that accompanies TEA exposure in LCs is accompanied by an increase in A-type K + current, we hypothesized that compensation may also occur at the network level to restore synchrony among LCs subsequent to the TEA block.

Compensation restores both excitability and synchrony following TEA exposure
To determine whether compensatory responses can restore both excitability and synchrony of LC output following TEA block, anterior LCs were superfused with TEA for at least 1 hr while the activity of the individual LCs and the network were continuously recorded. Our data demonstrate that both synchrony and excitability are restored towards baseline levels over a period of 30-60 min following TEA exposure. Figure 2A and B illustrate a typical progression through the loss and subsequent restoration of synchrony among LCs during one hour of continuous exposure to TEA. Acutely after application of TEA, we saw a significant reduction in synchrony as measured by R 2 (see Materials and methods; Ransdell et al., 2013a) across LC voltage waveforms ( Figure 2B I II  III  IV   I  II  III  IV   - I II III IV  I II III IV  I II III IV  I II III IV   10 Figure 2 continued on next page baseline levels, often by 30 min, and were no longer statistically different from baseline by 1 hr of treatment with TEA ( Figure 2C, right). To monitor compensatory changes in excitability, five different measures of excitability were calculated using both extracellular and intracellular recordings ( Figure 2D; see Experimental Methods). The spikes per burst, average spike frequency, total burst depolarization and burst amplitude all were significantly increased immediately after exposure to TEA ('acute'), while the small cell pacemaker (SC)-to-LC phase delay was significantly decreased ( Figure 2D). All 5 measures also then showed a significant change back towards their baseline levels between the 30 min and 60 min time points. While all measures showed clear shifts towards restoration of baseline excitability, the number of spikes per burst, spike frequency within each burst, burst amplitude, and total burst depolarization were not completely restored to control levels; the exception was the SC-LC phase delay which was fully restored ( Figure 2D). Preparations exposed to TEA for 2-3 hr showed no further change in excitability (data not shown). For this reason, time scales longer than 1 hr were not included in our analyses.

Modelling predicts compensation based on intrinsic conductances
TEA exposure reduces LC synchrony and induces hyperexcitability. Our model development and selection criteria resulted in a population of 27 model CG networks with variable underlying conductances of the constituent neurons that successfully recapitulated the biological data observed in TEA (see Methods, Supplemental Information). Our previous results identified an approximate 2.2 ± 0.8fold change in I A in LCs as a result of 60 min of TEA exposure (Ransdell et al., 2012). Therefore, we used the model networks to explore potential mechanisms of compensation by first increasing and decreasing each individual maximal membrane conductance by a similar factor of 2. We searched for changes that would increase LC spike synchrony while countering the hyperexcitability induced by TEA. To easily visualize the trends, each network was normalized to its initial value for spike synchrony. These data are shown for all conductances in Figure 3. Our initial goal with the model was to determine whether changes in single conductances were sufficient to elicit compensatory changes in output that help restore both excitability and synchrony. While it is not difficult to conceive of a change in multiple aspects of the parameter set that could achieve restoration of output, it is perhaps not as intuitive -but presumably the most parsimonious solution -for a single conductance to have such an impact. True to this expectation, while various manipulations of G max values improved either excitability or synchrony, very few conductance changes improved both. The optimal solution of significantly improving spike synchrony and also significantly decreasing the total number of action potentials was achieved in only one case: two-fold increase in G A resulted in a mean synchrony score that was significantly different from the TEA case (p<0.05, paired t-test) but not significantly different from control (p=0.157). No other change in a given conductance resulted in this combination of statistical outcomes. Not every model cell or networked improved uniformly with this conductance change. Therefore, while these results do not rule out a contribution for other conductances, they do suggest that an increase in G A , as seen in previous experimental studies on isolated LCs (Ransdell et al., 2012), may be the most likely candidate for a change in intrinsic conductance promoting synchrony at the network level. These data suggest that while a single conductance change (increased G A ) can help restore both excitability and synchrony, variations in a single voltage-dependent conductance may not be sufficient to account for correlation analyses of LC voltage waveforms for pairs of LCs from N = 11 preparations. Lines within boxes mark the median, box boundaries represent 25th and 75th percentiles, whiskers represent 5th and 95th percentiles, and points represent outlying observations. Groups with significant differences in median synchrony (p<0.05; Wilcoxon signed rank tests) are denoted with different letters. (D) Excitability of LCs was quantified by five measurements (mean ± SD). Analysis of each preparation used the average of 10 consecutive bursts at each time point (N = 8 preparations). Number of spikes per burst, spike frequency within each burst, and the latency between pacemaker firing and first motor neuron spike (SC-LC Phase Delay) were calculated from extracellular traces. Total depolarization and amplitude of each burst are based on intracellular recordings. Significant differences across groups (p<0.05; paired t-tests) are denoted with different letters, such that any two bars with a letter in common are not significantly different.   Somata contain 9 conductances: G CaS , G CaT , G LEAK , G CAN , G A , G BKKCa , G SKKCa , G Kd , and G NaP . Paired LCs (1+2, 4+5) have stronger local coupling (black resistor symbols), and all 5 LCs are reciprocally electrically coupled via weaker gap junctions (gray resistor symbols). An example of LC3 and LC4 model output within a network burst activity is shown in the red and blue traces under both control and TEA (90% reduction in both G Kd and G BKKCa ) conditions. Graphical representations of spike synchrony (raster plots) and waveform synchrony (scatterplots; as in Figure 2) are shown for the model neurons, demonstrating that both measures reflect the loss of LC synchrony as a result of TEA. (B) Measurements of both output variables (# of spikes and spike synchrony) were made under three model conditions: control, TEA, and TEA + either a 2x increase (G") or 2x decrease (G#) in a given conductance. N = 27 distinct model networks. All output measurements are normalized to their initial (control) conditions. Red lines represent the mean for a given group. Dashed line represents the 1.0 value (baseline) for a given measure. Compensatory responses that restore excitability and synchrony will tend to move the mean towards baseline. DOI: 10.7554/eLife.16879.005 the full compensation response. In addition to perturbing only individual conductances, we also varied current kinetics and activation parameters (half-activation voltage V 1/2 , ± 10 mV, and slope factor k, by 0.5x and 2x (Ballo et al., 2010) and time constant by ± 10 ms) for all the cell currents individually, and found that no changes in parameters for a single current could simultaneously restore excitability at the single cell level, and synchrony at the network level (data not shown). While the analysis has focused only on the parameters of a single current, simultaneous changes in parameters of multiple currents could also potentially provide similar compensation, and that remains to be explored. However, our analysis does reveal the substantial contribution of changing a single parameter -G Aon multiple aspects of network compensation, to an extent that is beyond simple intuition. Importantly, the model also extends the biological data by demonstrating that waveform synchrony can translate into spike synchrony. Because of the electrotonic distance between the somata and axons of LCs, we cannot measure spike synchrony directly in this preparation. The model allows us to infer that waveform synchrony (and loss of synchrony) can indeed translate to the level of the most proximal cellular output -spiking.

Intrinsic compensation contributes to restoration of synchrony
Model runs predicted that increases in I A help restore LC excitability and synchrony. To test this experimentally, we silenced pacemaker activity with isotonic sucrose solution (see Methods, Supplemental; Figure 1A), and tested the similarity of responses of each individual LC to a biologically realistic current stimulus (Ransdell et al., 2013a). We compared LC3 and LC5 to the same current injection at three time points: control, 5 min post-TEA, and 1 hr post-TEA. Between current injections, pacemaker activity was restarted by removal of the sucrose block. This allows us to test each cell in isolation, but compensation occurs in the intact network. The initial voltage responses to our stimulation protocol in LC3 and LC5 in control conditions are highly similar to one another, and their level of waveform synchrony was not significantly different from the synchronous activity across these LCs during intact control network activity ( Figure 4A,B). Immediately following TEA application, LC3 and LC5 show disparate output when driven with a common stimulus protocol ( Figure 4A). Finally, our data show significant increases in R 2 of voltage activity within 1 hr across isolated LCs ( Figure 4B), demonstrating that intrinsic compensation does improve network synchrony. However, after 1 hr the synchrony values were significantly lower than control values ( Figure 4B,C), suggesting that intrinsic compensation alone is insufficient to restore synchrony. To determine whether compensatory changes in I A occur in the intact network, we measured I A with two-electrode voltage clamp in LCs before and after 1 hr of TEA exposure. Measurements were made while the network activity was temporarily halted with sucrose block, and compensation occurred with ongoing network activity. In all cases peak I A current increased ( Figure 4D), with a mean increase of 56% (p<0.05, n = 6, Wilcoxon signed rank test). These data are consistent with the hypothesis that a compensatory increase in I A can help promote synchrony in these networks. Although the waveforms of LC3 and LC5 were not different from one another after 1 hr of compensation, anterior LC burst potentials did not reproduce their original waveform after compensation (see Figure 2; Panels I and IV). We also used the LCs from data shown in Figure 4A to compare the voltage responses of individual LCs to a fixed stimulus before and after compensation. Repeatable voltage responses under control conditions indicate that trial-to-trial variability is negligible (mean R 2 = 0.997; Figure 4E, control). However, the voltage response after compensation was significantly different from the control voltage response (p<0.01, n = 8, Wilcoxon signed rank test), indicating that intrinsic compensation does not restore the original cellular output ( Figure 4E; control vs. TEA [1 Hour]).

Increased model electrical coupling conductance helps restore synchrony
If intrinsic compensation does not fully restore synchrony, another mechanism must be present to explain the results observed during network compensation. LCs receive common excitatory inputs from the pacemakers and one hypothesis is that changing the strength of these chemical synapses might help to restore LC firing to appropriate levels. LCs in the network are also electrically coupled to one another via gap junctions which presumably promotes synchrony, although clearly the native coupling is not able to maintain LC synchrony in TEA (Hagiwara et al., 1959; Figure 4. Intrinsic compensation involving G A partially restores synchrony after TEA block. A reversible sucrose block was used to temporarily stop network activity and a current stimulus protocol was delivered to individual LCs at three time points: control, after 5 min of TEA perfusion (acute), and after 1 hr of TEA perfusion (after compensation). (A) Representative traces from 2 different preparations compare the responses of LC3 (orange) and LC5 (blue) to the same current injection (Stimulus Protocol) at each of the three time points. (B) R 2 values for N = 6 preparations at each time point are plotted with the same preparation connected across time points. All 3 conditions were significantly different from one another (p<0.01; paired t-tests). (C) When isolated cell output was compared with the output in the network, there was no difference in R 2 at the control state, but following 1 hr of compensation in the network, there was a significant (p<0.002t-test) difference in synchrony scores (mean ± SD) between cells when isolated vs. when they are in the intact network. (D) I A was measured by two-electrode voltage clamp before and after compensation in N = 6 LCs in the intact CG. Voltage clamp data were obtained by temporarily silencing network activity with 750 mM sucrose. There was a significant increase in I A after 1 hr in TEA (mean 56 ± 65% increase, p<0.05 -Wilcoxon signed rank test). (E) Similarity of waveform in the same neurons before and after 1 hr TEA exposure. R 2 values were calculated for the output of the same cells before and after TEA exposure, and are shown as before-and-after values in the same cell connected by a line. Box plots show distributions of R 2 values from cross-correlation analysis of LC voltage waveforms N = 8 cells. Lines within boxes mark the median, box boundaries represent 25th and 75th percentiles. 'Control' is the comparison of voltage waveforms to 2 separate rounds of current injection in the absence of TEA. Although there was improvement in similarity of waveforms after 1 hr in TEA (panel A), the newly compensated output after 1 hr of TEA exposure does Figure 4 continued on next page 1983; Cooke, 2002). A second hypothesis is that increased electrotonic coupling between LCs could buffer against disparate output and help to restore synchrony.
Using our set of model networks, we increased and decreased the strength of chemical synapses in 10% increments to test the effects on excitability and synchrony. We then did the same with model electrical coupling conductance. We found that increasing the strength of either chemical synapses or electrical coupling increased both synchrony and excitability ( Figure 5 left). However, increasing the chemical synaptic conductance in conjunction with TEA blockade also increases spiking of the LCs~25-30% in contrast to the biological decrease in excitability relative to the acute TEA exposure seen with compensation. Conversely, only a small change in LC spiking occurs with an increase in electrical coupling (~9%, Figure 5 left). Reducing the strength of either chemical synapses or electrical coupling decreased overall spike synchrony ( Figure 5 right), violating the assumptions of compensation based on the biological data. Reducing chemical synaptic strength eventually ceased LC firing altogether (data not shown).

Increased electrical coupling helps restore synchrony across LCs in experiments
Increased electrical coupling restored synchrony in model LCs with only a modest effect on excitability. To investigate this relationship in experiments, we isolated LCs and used a dynamic clamp to   add an artificial coupling conductance. Pairs of LCs from the same network were physically isolated by thread ligature, exposed to TEA, and simultaneously received the same stimulus protocol, while the dynamic clamp added a non-rectifying artificial coupling conductance (from 0 to 0.2 mS) between the cells. The driving force was equivalent to the voltage difference in membrane potential between the coupled cells. Increasing the artificial coupling conductance significantly increased the correlation coefficient of the waveform between the two cells ( Figure 6A), with a synaptic conductance value of 0.2 mS able to rescue synchrony of LCs to levels observed in intact networks ( Figure 6A). These results provided proof of principle that increasing electrical coupling could be responsible for resynchronization in the network. With acute exposure to TEA, isolated LCs produce disparate output in response to an identical stimulus. The stimulus protocol consists of current injections that mimic biological synaptic currents and back propagating action potentials (see Ransdell et al., 2013a) from four consecutive bursts of network activity. The current was injected simultaneously into isolated cells, while the dynamic clamp was used to provide an artificial coupling conductance. Representative traces of the same two cells shown in TEA with different levels of synaptic current applied. Only the final burst of the fourburst input stimulus is shown for clarity. N = 7 different preparations (bottom panel) show an increase in synchrony with increasing coupling conductance. (B) Biological coupling increases during network compensation. Hyperpolarizing current injections were used to measure coupling coefficients between LC3 and LC5 during control conditions, with acute TEA exposure (5 min) and after 1 hr compensation in TEA. The input resistance of LCs was measured and showed no significant change at any time point. No significant differences were observed for coupling coefficient or coupling conductance between control and acute TEA conditions. Coupling coefficient significantly increased after 1 hr in TEA (Mean increase 85 ± 82% from control, N = 11, p<0.01, Wilcoxon signed rank test). Measurements from LC4 and LC5 show that coupling conductance increased significantly as a result of 1 hr TEA exposure (mean increase 49.5 ± 36% from control, N = 13, p<0.001, Wilcoxon signed rank test). Significant differences across groups are denoted by different letters. Plots show mean ± SD. (C) Representative traces from two different preparations of changes in coupling observed before (Control) and after 1 hr of TEA exposure. Top traces are between LC3-LC5 and the bottom traces are between LC4-LC5. Measurements were made in two-electrode current clamp, and the current was injected into LC5. Because current injections were manually timed to occur between bursts of network activity, slightly different durations of current pulses occurred in the two recordings in the top recordings. Recordings from LC4-LC5 were used for coupling conductance measurements seen in panel B, as their close proximity allows for much less influence of electrotonic distance on calculations of conductance. DOI: 10.7554/eLife.16879.008 We then measured coupling coefficients between LCs during compensation in the intact network. Coupling coefficients between cells increased significantly after 1 hr in TEA (mean increase 85%, p<0.01; Figure 6B,C). The coupling coefficient is a useful description of the functional relationship in coupling, but does not identify the electrophysiological mechanism. Plasticity of coupling properties can ultimately be influenced by two fundamentally different mechanisms: altered resistance of the non-junctional membrane of the coupled cells, or modification of gap junctional conductance. Using two-electrode current clamp, we saw no significant differences in the apparent input resistance of LCs in control physiological saline, after acute TEA exposure, or after compensation (1 hr TEA exposure; Figure 6B,C). These results indicate that changes in passive membrane conductance (G Leak ) are not responsible for increased coupling coefficients.
To directly test whether coupling conductance between LCs increases as a result of TEA-induced compensation, we focused on LC4 and LC5 pairs within the same network. In the crab CG, anterior LCs exhibit strong local electrical and dye coupling (Tazaki andCooke, 1979, 1983). The branch containing LC4 and LC5 somata can be separated and electrotonically 'sealed' from the network by thread ligature to create ideal conditions for measuring coupling conductance. With two electrodes in each cell, we used hyperpolarizing current injections to measure resistance and calculate junctional conductance independent of membrane resistance (as in Bennett, 1966; see Materials and methods). Coupling conductance (G C ) between LC4 and LC5 significantly increased during 1 hr of TEA exposure (mean increase 49.5%, p<0.01, N = 8; Figure 6B,C).

Interaction of intrinsic and electrical synaptic compensation
Taken together, our experimental and modelling results suggest that an increase in G A is able to counter the increase in excitability of LCs in TEA in a compensatory fashion, as well as promote restoration of synchrony, but was insufficient to restore synchrony fully. Additionally, our results suggest that an increase in coupling among LCs can greatly promote synchrony with only a modest effect on excitability. Therefore, we next used our model networks to investigate how G A and G C might interact to promote synchrony by calculating synchrony scores as conductances of all 27 model networks were adjusted. First, we increased G A alone in 10% increments up to a 100% increase (Figure 7). Increasing G A up to +40% promoted greater synchrony after TEA blockade, but was unable to fully restore synchrony even with increasing conductance levels, consistent with our biological data (Figure 4). Increasing G A beyond +40% did not further improve synchrony (Figure 7), and ultimately caused LCs to cease firing altogether. We also increased model G C incrementally (from +10% to +150%), and found that electrical coupling alone was capable of restoring synchrony fully, but this required a 140% increase in its value (Figure 7). Finally, we increased both G A and G C together in 10% increments, revealing a potentially synergistic relationship: a smaller increase of 70% in each conductance was able to produce spike synchrony that was indistinguishable from control ( Figure 7).

Discussion
Homeostatic or compensatory plasticity in the nervous system has been the subject of intense interrogation, with studies focusing on both homeostatic synaptic scaling (Turrigiano, 2012;Lee et al., 2014) and tuning of ionic conductance relationships to maintain a target level of excitability (Turrigiano et al., 1994;Desai et al., 1999). However, few studies have sought to integrate multiple mechanisms to directly address emergent network stability from compensatory processes acting at the level of single neurons. Moreover, only recently has there been an appreciation that variability in underlying cellular parameters such as conductance magnitudes (Schulz et al., 2006;Pratt and Aizenman, 2007;Wilhelm et al., 2009), activation properties of channels (Olypher et al., 2006;Amendola et al., 2012), and synaptic strengths (Olypher and Calabrese, 2007;Wilhelm et al., 2009;Grashow et al., 2010) may form part of the repertoire of compensatory mechanisms that endow networks with remarkable robustness characteristics. Homeostatic plasticity has been reported in single isolated neurons (Swensen and Bean, 2005;Ransdell et al., 2012), in artificial networks formed in culture (Desai et al., 1999;Ibata et al., 2008), or in thin sections of the CNS (Karmarkar and Buonomano, 2006;Lambo and Turrigiano, 2013).
Cortical slice preparations have demonstrated homeostatic plasticity of intrinsic conductances coordinated with chemical synaptic plasticity (Karmarkar and Buonomano, 2006;Lambo and Turrigiano, 2013), and Mauthner cells in fish have been shown to exhibit coordinated activitydependent changes in both the chemical and electrical components of its mixed synapses (Yang et al., 1990;Pereda and Faber, 1996). While membrane conductances and properties of electrical coupling are known to interact in critically important ways to promote synchrony (Curti et al., 2012;Gutierrez and Marder, 2013), to our knowledge the present study is the first demonstration of coordinated homeostatic plasticity of intrinsic and electrical synaptic conductances in a comprehensive network-level example of homeostatic compensation.
We hypothesized that plasticity distributed throughout the network might provide robustness of network output. Using a combination of pharmacology, electrophysiology, and modelling approaches in intact neural networks, we first uncovered a striking vulnerability of neural networks to underlying cellular variability. We then demonstrated that an intact network has mechanisms that allow for robustness in the output via synergistic compensatory changes across the individual cell (ionic conductance) as well as trans-cellular (electrical coupling) properties. The interplay of intrinsic and synaptic parameters, and distributed plasticity in determining network output, is difficult, if not impossible, to unravel experimentally. This makes investigations via biophysical models an attractive alternative to narrow possibilities, and complement experiments, as we have demonstrated.

Multi-component compensation can synergistically restore network output
Previous modelling studies found that K + currents can increase or help restore synchrony between electrically coupled neurons (Pfeuty et al., 2003), so we first hypothesized that a compensatory increase in A-Type K + membrane conductance could be a mechanism underlying both restored excitability and resynchronization. Over the course of 30-60 min, increased I A was associated with decreasing cellular excitability [see also (Golowasch et al., 1999)] and improvement of coordinated motor neuron firing. However, intrinsic compensation alone was insufficient to fully restore synchrony across LCs. A concomitant increase in electrotonic coupling ensured virtually complete resynchronization. Our modelling results suggest that although a sufficient increase in electrical coupling alone could restore full synchrony (140% increase), it could not simultaneously restore the original level of excitability. Only a 70% increase was necessary when accompanied by a concomitant increase in G A . Therefore, we conclude that multi-component mechanisms are not only necessary for full compensation, but also that their synergistic action is potentially more efficient than either mechanism operating in isolation.

Comparison of our results in the context of behavioural/motor plasticity
While our results occurred in a compensatory context, the underlying mechanisms bear striking similarity to motor output plasticity induced by operant conditioning in Aplysia. In a series of elegant experiments, it has been shown that chaotic exploratory and consummatory radula biting movements of Aplysia during food searching behaviour can be stably modified by operant conditioning, leading to prolonged bouts of radula movements with increased frequency and more stereotyped rhythmic organization (Nargeot et al., 2007;Nargeot and Simmers, 2011). Chaotic biting patterns result from inherently variable and unsynchronized bursting of CPG neurons that are each randomly capable of triggering bites (Nargeot et al., 2009). Following operant conditioning, induction of regular and synchronized bursting of pattern-initiating cells can be attributed to changes in both intrinsic excitability and electrical coupling strength. Specifically, changes in intrinsic excitability attributed to changes in leak conductance underlie the increase in frequency of motor output, while increases in coupling strength allow for the synchronization and regularization of bursting (Nargeot et al., 2009;Sieling et al., 2014). The full shift in behavioural and circuit output is therefore the additive influence of both intrinsic and electrical synaptic conductances. The striking similarity in these underlying mechanisms suggests that these kinds of circuit-level mechanisms may be a conserved strategy for stabilization of synchrony within network output, be it compensatory or in the context of behavioural plasticity.

Physiological regulation of coupling conductance
The speed (within 30 min) and magnitude (up to a doubling of effective coupling) of physiological changes seen in electrical coupling was remarkable. Although electrical coupling has long been known to promote synchrony in many systems, including the CG (Tazaki, 1972;Bennett and Zukin, 2004), the physiological interaction of electrical coupling with intrinsic conductances to affect a compensatory output has not been examined. Previous work in the crustacean STG has demonstrated how the synchronized activity of pacemaker cells is dependent on an interaction of intrinsic conductances and electrical coupling (Szücs et al., 2000(Szücs et al., , 2001Soto-Treviño et al., 2005), and that distinct circuits can be brought into synchrony via manipulations of electrical and chemical synapses (Elson et al., 1998;Szücs et al., 2000Szücs et al., , 2009. But none of these studies have addressed the interaction of membrane conductance and electrical coupling in a compensatory context. Similarly, plasticity of electrical synapses has drawn considerable attention after being discovered in the mammalian central nervous system, including the thalamic reticular nucleus (Landisman and Connors, 2005;Haas et al., 2011), inferior olive (Lefler et al., 2014;Mathy et al., 2014), and retina (Kothmann et al., 2009;Vö lgyi et al., 2013). Studies in the thalamic reticular nucleus have suggested that potentially compensatory changes in coupling are important to maintain network stability as large changes in intrinsic excitability occur across development (Parker et al., 2009). These discoveries increased awareness of the complex functional roles and plasticity of coupling (Pereda et al., 2013;O'Brien, 2014;Haas, 2015), and also spurred research to identify molecular mechanisms that underlie plasticity and maintenance of these structures (Flores et al., 2012;Li et al., 2012;Turecek et al., 2014). Our study adds to this growing appreciation for plasticity of electrical synaptic connections in the context of homeostatic plasticity.
Variability, plasticity, and network output -the bigger picture Stable levels of intrinsic neuronal excitability and temporal coordination within networks are critical features across all nervous systems. Underlying both neuronal and network outputs are complex, and often highly variable intrinsic and synaptic properties of constituent neurons (Marder, 2011;Norris et al., 2011). Our data demonstrate that the intrinsic variability between cells of the same type can make networks vulnerable to loss of temporal coordination, in this case desynchronization of motor neuron output. Although LC activity was fully resynchronized within 1 hr, recovered LCs never recapitulated their original voltage waveforms. While the intrinsic conductances involved in our manipulation and compensation (G Kd , G BKKCa , G A ) have overlapping functions and characteristics (Ransdell et al., 2012), our findings demonstrate that individual conductances are not truly redundant. Degeneracy of ion channel properties leading to this type of relationship has been put forth as a mechanism underlying robustness and adaptability in neural networks (Tononi et al., 1999;Marder and Goaillard, 2006), but our study suggests physiological limits to neural network compensation and robustness. These limits may themselves be a contributing factor to the nature and progression of pathology in neurodegenerative diseases (Trasande and Ramirez, 2007;Beck and Yaari, 2008;Small, 2008).
We induce desynchronization and compensation in our studies through pharmacological block of a subset of K + conductances with TEA. However, the precise role these mechanisms play in fully intact biological networks is unclear. Intrinsic conductances can be differentially affected by ubiquitous natural mechanisms such as neuromodulation (Marder, 2011(Marder, , 2012 or temperature changes (Tang et al., 2012;Marder et al., 2015). Further, the effectiveness of electrical coupling can be affected by the modulation of intrinsic cellular conductances (Szabo et al., 2010). Maintaining reliable synchronization of the output under changing conditions is not trivial, and understanding the robustness and the constraints of homeostatic systems that cope with such perturbations remains an important area for future investigation (Marder et al., 2014).

Animals
Adult male Jonah crabs, Cancer borealis, were shipped overnight from The Fresh Lobster Company (Gloucester, MA). Crabs were maintained in artificial seawater at 12˚C until used. Crabs were anaesthetized by keeping them on ice for 30 min prior to dissection. The complete CG was dissected from the animal and pinned out in a Sylgard-lined petri dish in chilled physiological saline (440 mM NaCl, 26 mM MgCl 2 , 13 mM CaCl 2 , 11 mM KCl, and 10 mM HEPES, pH 7.4-7.5, 12˚C). Chemicals were obtained from Fisher Scientific unless otherwise noted.

Biological methods: Electrophysiology
The CG network is comprised of 9 cells: 4 Small Cell (SC) pacemaker interneurons which give simultaneous excitatory input to 5 Large Cell (LC) motor neurons. Superfusate of SCs can be separated from the anterior LC somata using petroleum jelly wells ( Figure 1A). Intact network activity was monitored with intracellular recordings from anterior LCs along with extracellular recording of the network output. For most experiments, the posterior end of the ganglion was maintained in normal physiological saline and protected from TEA superfusate with a barrier of petroleum jelly. The anterior end of the preparation was superfused at a rate of approximately 2 ml/min. A schematic of this experimental setup is shown in Figure 1A. All experiments were performed at 12˚C. Extracellular recordings using a Model 1700 Differential AC Amplifier (A-M Systems, Carlsborg, WA) were taken with stainless steel pin electrodes from a petroleum jelly well on the ganglionic trunk containing axons of all 9 neurons in the CG. LC spikes on the extracellular traces are easily distinguishable by their large amplitude. The LC somata were desheathed for sharp electrode recordings. Intracellular recordings were made using glass electrodes containing 3 M KCl (8-25 MW) and AxoClamp 900A and AxoClamp 2B amplifiers (Molecular Devices, Sunnyvale, CA). Two-electrode voltage clamp (TEVC) and two-electrode current clamp (TECC) protocols were created and run using Clampex 10.3 software (Molecular Devices).
Somata were isolated for dynamic clamp experiments by tightening a thread ligature past the anterior branch point on the nerve containing the LC soma. Isolated cells were simultaneously driven with the same current stimulus, described previously (Ransdell et al., 2013a). Briefly, a stimulus protocol was generated by recording the voltage waveform from a LC somata during intact network activity. This consisted of a 20 s recording from a LC3 soma which included four burst potentials with both pacemaker cell EPSPs and LC back-propagating APs present. In addition to the stimulus current, dynamic clamp artificial coupling current was applied with NetClamp software (developed in the Fishberg Department of Neuroscience of the Mount Sinai School of Medicine and available at http://gothamsci.com/NetClamp/) at a sampling rate of 50 kHz according to the equation: I gap = g* (V m 1 -V m 2) where g is a non-rectifying coupling conductance under experimental control, and the voltage difference between the cells (V m 1 -V m 2) determines the driving force.
I A was measured before and after compensation using voltage clamp protocols as described previously (Ransdell et al., 2012). Briefly, outward currents were measured from a holding potential of À30 mV and stepped from À50 mV to +5 mV in 5 mV increments in order to measure the high threshold K + current I HTK which is blocked by TEA. In LCs, I HTK is predominantly a mix of BKKCa and delayed rectifier currents (Ransdell et al., 2012). A-Type K + current (I A ) was measured by performing an identical voltage clamp steps from a holding potential of -80mV and subtracting I HTK . P/N leak subtraction was used for all TEVC. Coupling in the intact network was measured using TECC in both cells during the sucrose block. Negative current steps ranging from 1-6 nA were injected into one cell at a time while measuring voltage changes in both cells. Coupling coefficients were calculated as the ratio: (DV coupled cell / DV Injected Cell ).

Biological methods: Measurements of intrinsic compensation
To determine whether intrinsic compensation may be contributing to restoration of synchrony in the biological network, we used a current stimulus protocol simulating realistic network inputs to LCs in order to deliver the same biologically relevant stimulus at 3 time points (Ransdell et al., 2013a). Using this reversible sucrose block to suspend pacemaking activity, we were able to compare the similarity of responses of LC3 and LC5 to the same current injection at three time points. This allows us to test each cell in isolation with respect to its output waveform, but compensation occurs with full network activity after removal of the sucrose over the course of 1 hr. We measured individual LC responses to current injection in control saline, repeated after 5 min of TEA exposure (acute), and again after 1 hr of TEA exposure.
Our modelling results suggested that an increase in I A may act in a manner to accomplish both the decrease in excitability and the restoration of synchrony seen in our biological data. In order to determine whether compensatory changes in I A occur in the intact network, we performed two-electrode voltage clamp in the same LC before and after 1 hr of TEA exposure. In order to track changes in I A in individual LCs we again used reversible sucrose block to perform voltage clamp immediately after acute TEA application and again after 1 hr of TEA perfusion. After voltage clamping anterior LCs, the sucrose was washed out and replaced with physiological saline, allowing the network to resume normal activity. This process was repeated after 1 hr of compensation.

Biological methods: Data analysis
Intracellular burst waveforms were considered to begin with the first EPSP from pacemaker activity and ended upon return of the waveform to resting membrane potential. Recordings were analyzed using Clampfit 10.3 (Molecular Devices) and Spike 2 version 7 (CED, Cambridge, UK) software. Statistical analyses were performed using SigmaPlot 11.0. Correlation coefficients (R-values) were obtained by a Pearson correlation, and squared to calculate the coefficient of determination (R 2 ). Most data are 'before and after' effects within the same ganglion or cell, and therefore any two groups were compared with paired t-tests when the data were normally distributed, or Wilcoxon signed rank tests in the case of non-normality. The sample sizes for comparison of waveform synchrony were calculated with power analyses based on projected means and standard deviations from data reported in our previous study with very similar experimental manipulations of TEA exposure of LCs (Ransdell et al., 2013a), which yielded target sample size of N=6-10 to yield a power of 0.8 to 0.97. Sample sizes for changes in network output and changes in I A after TEA exposure were based on similar data in our previous work (Ransdell et al., 2012), and yielded target sample sizes of N=5 to achieve a power of 0.909. Power analyses were conducted based on the use of paired ttests to analyze the data. However, when data were not normally distributed we ended up using a Wilcoxon signed rank test, which was not utilized in our initial power analyses. All sample sizes used in our studies are reported in Figure Legends and/or in the Results section when significance values are reported.
To quantify synchrony of LC voltage waveforms, we performed a cross-correlation of the digitized voltage signal from LC pairs, as shown in Figure 2B. The first pacemaker spike was used to define the start of each LC burst, and bursts were considered to have terminated upon return to V Rest in the LCs. The coefficient of determination from this cross-correlation (R 2 ) was used to examine how accurately one burst waveform could predict the waveform in another LC (see Ransdell et al., 2013a). This cross-correlation was performed for every burst across the full time-course of the experiment ( Figure 2C, left). A decrease in R 2 therefore indicates a loss or reduction in waveform synchrony. R 2 values from 10 consecutive bursts were averaged for all data points presented as waveform synchrony data, save for the individual points found in Figure 2B and C.
Five different measures of excitability of LCs were calculated using both extracellular and intracellular recordings ( Figure 2D). Extracellular recordings from the ganglionic trunk were used to calculate the number of spikes per burst and spike frequency. Because our hypothesis predicts both increased LC excitability and desynchronization with TEA exposure, it should be noted that there is a potential confound in distinguishing these effects based on extracellular analysis alone. Axons of all LCs run through the ganglionic trunk, thus the increased spike count observed could result from an increase in the total number of action potentials, desynchronization of action potentials across LCs, or both. We therefore included three additional measures of excitability that helped clarify the effect. We measured the latency between the onset of SC pacemaker bursting to the first LC spike in each burst ( Figure 2D; SC-LC phase delay). Two other measures of excitability were calculated from intracellular voltage changes: burst amplitude and total burst depolarization. Burst amplitude was defined as the maximal voltage change from V Rest to the highest peak of the burst, and total depolarization is the area under the curve above V Rest , measured in mV*sec. For all measures of excitability, values from 10 consecutive bursts were averaged at each time-point in each preparation (N=8).

Model methods: Single cell models
A detailed model was created with eight voltage-dependent conductances (G A , G Kd, G NaP, G CaS, G CaT , G CAN, G SK(Ca) , G BK(Ca) ) and passive leak channels. These conductances are defined as: A-type potassium (G A ), delayed rectifier (G Kd ), persistent sodium (G NaP ), transient calcium (G CaT ), slow persistent calcium (G CaS ), calcium-dependent non-selective cation (G CAN ), two calcium-dependent potassium currents (G SKKCa and G BKKCa ), and leak (G Leak ). A single compartment model with biological dimensions for soma (Ransdell et al., 2010(Ransdell et al., , 2013b was created in NEURON and its capacitance (C m ) and leak conductance (G leak ) were tuned to match the observed biological membrane time constant (t) and input resistance (R in ). This resulted in a soma with a length of 284.87 mm and a diameter of 125 mm, with a capacitance of 2.719 mF/cm 2 . Channels were then added and their maximal conductances were tuned to match three biological properties observed, i.e., a) Total outward current b) Response to synaptic drive and c) Response to synaptic drive in the presence of TEA. These results were obtained from experiments performed on ligatured somata of C. borealis LCs (Ransdell et al., 2012(Ransdell et al., , 2013a(Ransdell et al., , 2013b. For the network studies another compartment termed Spike Initiation Zone (SIZ) was added to the model. This compartment was modelled as a cylinder with a capacitance of 1 mF/cm 2 , a length of 400 mm, and a diameter of 8 mm, and contained only sodium, potassium and passive leak channels. This compartment allowed us to obtain spiking activity in the model cells for spike synchrony analysis. The resulting model equations were as follows: C dV dt ¼ ÀI A À I Kd À I NaP À I CaS À I CaT À I CAN À I SKKCa À I BKKCa À I Leak ðSomaÞ C dV dt ¼ ÀI Na À I Kdr À I Leak ðSIZÞ The individual currents were modelled as I c ¼ g max;c m p h q V À E c ð Þ, where g max;c is its maximal conductance, m its activation variable (with exponent p), h its inactivation variable (with exponent q), and E c its reversal potential (a similar equation is used for the synaptic current but without m and h). The kinetic equation for each of the gating functions x (m or h) takes the form where x ¥ is the steady state gating voltage-and/or Ca 2þ -dependent gating variable and t x is the voltage-and/or Ca 2þdependent time constant. The equations for the active channels in the soma compartment were fit using biological recordings for these currents from the cardiac ganglion of Cancer borealis. These currents were fit as follows: Voltage clamp data obtained with Clampfit were imported into MATLAB (Mathworks, Natick, MA) and fit using the MATLAB curve-fitting toolbox. Current data were converted to conductance data by dividing by (V m -E Rev ), where E Rev was as follows: E Na = +55 mV, E K = À80 mV, E Ca = +45 mV, E Leak = À50 mV, and E CAN = À30 mV. The time axis was adjusted to start from 0 for the beginning of the clamp. The following parameterization was used: In this equation, A i ¼ G i;max Â m i was the maximal conductance of the current i multiplied by its voltage-dependent steady-state activation (m i ), h i was the steady-state inactivation value, and t m;i and t h;i were the time constants with which activation and inactivation reached steady-state, respectively. This fitting procedure assumed that ionic currents were completely deactivated (m = 0) and de-inactivated (h = 1) prior to the onset of the voltage clamp. This was fit to each trace in a voltage clamp experiment, giving the values of each of the four parameters for each test clamp voltage (V c ). These values were then fit for each current as functions of V c using the general forms as stated below. This procedure yielded equations for the currents recorded in voltage clamp that could be used in simulations according to the Hodgkin-Huxley mathematical formalism.
All the maximal conductances (G i;max ) were in S, time constants in ms and voltages in mV.

Model methods: Calcium dynamics
Intracellular calcium modulates the conductance of the calcium-activated potassium currents (BKKCa and SKKCa), calcium-activated nonselective cation current (CAN), and influences the magnitude of the inward calcium current in the LC (Tazaki and Cooke, 1990). A calcium pool was modelled in the LC with its concentration governed by the first-order dynamics (Prinz et al., 2003;Soto-Treviño et al., 2005) below: where F ¼ 0:256 M=nA is the constant specifying the amount of calcium influx that results per unit (nanoampere) inward calcium current; t Ca represents the calcium removal rate from the pool; and ½Ca 2þ rest ¼ 0:5 M. Voltage-clamp experiments of the calcium current (Ransdell et al., 2013b) showed the calcium buffering time constant to be around 690 ms (t Ca ).  (Tables 1,2), we wanted to search the conductance space for other possible conductance combinations that might exhibit appropriate LC output. The properties that had to be maintained were; a) Input Resistance (R in ) and Resting Membrane Potential, b) Pre-TEA and Post TEA response to current injection c) Response to Synaptic drive obtained from biological cell.
The rules used to select the potential parameters were as follows (based on biological recordings): Synaptic Drive response should have an R 2 value of at least 0.8 or higher when compared to biological Synaptic Drive response. The duration of the pre-TEA response to a 6 nA, 50 ms current injection should be less than 120 ms. Also the peak should be less than À22 mV. The duration of the post-TEA (G BKKCa and G Kd reduced by 90%) response should be between 255-667 ms and its peak should be greater than À15 mV. A 9-D max conductance parameter space (5-fold variation over each conductance except G Leak ) was searched randomly for sets that satisfied the constraints above. We searched 20,000 different combinations of parameter sets with these criteria, and most of those which passed did not have a proper termination of activity following current injection (i.e., did not return to V rest ). We concluded this was due to an inappropriate relationship between I CAN and I SKKCa . Subsequent trials revealed that a given ratio range (~1:0.83 respectively) of these two currents was necessary for proper termination of the activity. Larger ratios cause V rest to be higher due to the reversal potential (À30 mV) of CAN current. A higher fraction of I SKKCa (reversal potential À80mV) caused a large AHP after termination and reduced the duration of the post-TEA response. Using the updated selection criteria with a ratio I CAN to I SKKCa , we found 180 parameter sets that passed. Of these 180 potential model sets, we selected only the ones that had Synaptic Drive response R 2 value > 0.9 compared to the biological Synaptic Drive response. This resulted in 49 potential parameter sets.
Biological data showed that I A and I BKKCa had a negative correlation in their magnitudes in LCs (Ransdell et al., 2012). We added this to our criteria for screening potential parameter sets for the network studies. We converted biological I A -I BKKCa current data into factor data by dividing I A and I BKKCa by their respective factor average. G A and G BKKCa values of passed parameters were similarly divided by its average to get its factor data. The biological data were fitted using a linear polynomial from 95% to 70% confidence intervals, in steps of 10%. For network studies we used a 70% confidence interval, which left us with 14 potential parameter sets that represented LC model neurons for use in modeling studies.
Model methods: Development and validation of a population of conductance-based model networks for studying mechanisms restoring network synchrony Our results demonstrate that intact cardiac ganglia are able to compensate for the loss of highthreshold K + currents and restore both excitability and synchrony within one hour of TEA blockade. We next set out to explore the mechanisms by which excitability and synchrony could be restored in this network. To maximize our ability to interrogate multiple parameters that may be responsible for compensation in this system, we constructed a population of conductance-based biophysical models of the CG network. This allowed us to simulate the TEA conductance blockade and then manipulate individual conductances, both voltage-gated and synaptic, to examine their effects on network excitability and synchrony.
Our 14 parameter sets for LCs were used to create 50 random 5-cell networks of LCs, ensuring that the same model LC never appeared twice in the same network. The five cells within a network were then electrically coupled using conductance values tuned to reflect experimental observations of coupling coefficients. Small cell (SC) pacemaker drive was simulated as excitatory synapses via the NetStim function in NEURON. Parameters for the model of the synaptic drive onto LCs were tuned to get 6 to 9 spikes in the nominal LC model. It was observed biologically that frequency of SC firing increases within the slow wave oscillation cycle of LCs. Based on these recordings, the model SC burst initially fired at 18 Hz for first 440 ms and then increased to 25 Hz for 560 ms, with the burst terminating at 1000 ms. Our experimental TEA block was simulated in these networks by reducing G BKKCa and G Kd conductances by 90% in the 3 anterior LCs based on biological data in LCs (Ransdell et al., 2013a). We imposed a final set of selection criteria on the randomly generated model networks, rejecting networks that increased synchrony or decreased the total number of spikes after the simulated TEA block, as this was never observed in biological networks. This left 27 networks that reproduced the biological trends and these were used in subsequent analyses to explore potential conductance changes that could restore network synchrony.
Somatic burst potentials drive action potentials in LC axons, so divergent burst waveforms would be expected to cause desynchronized spiking. Our biological data qualitatively agreed with this, but a precise quantification of synchrony for all spikes within a burst is subject to many ambiguities. Our model networks easily provided precise spike times for each cell in the network, so we chose to examine actual spike synchrony in the model to complement the burst waveform analysis in the biological preparation. Our analysis considered synchrony for paired anterior LCs with a nominal coupling conductance of 0.0182 S using a 25 ms bin width for spike-times (Wang and Buzsáki, 1996). Spikes occurring in both cells during the same bin were considered synchronized, while spikes that did not bin together were tallied as desynchronized. Using the definition of synchrony listed in the next section, these randomly generated model networks exhibited 'control' synchrony scores ranging from 0.642 to 1.0 with a median value of 0.915 (matching data in biology), where 1.0 represents perfect spike synchrony.

Model methods: Data analysis
In the models, spike synchrony between two cells was calculated based on spike times (Wang and Buzsáki, 1996). Spike times were recorded from each LC's Spike Initiation Zone (SIZ). The simulation time was divided into 25 ms bins. After initializing all bins to zero, each cell spike was added to the corresponding bin. Synchrony (SY) between two cells A and B was calculated using the following equation: where l is the current bin and k is the maximum number of bins. Spikes occurring in both cells during the same bin were considered synchronized, while spikes that did not bin together were tallied as desynchronized.
To compare the measures of spike and waveform synchrony in model networks, model waveform synchrony measures (Figure 3) were performed as described above on filtered voltage traces (Gaussian lowpass filter, 15 Hz cutoff frequency; Clampfit 10.3) that remove the influence of the axonal spikes on the voltage waveforms. Statistical analyses were performed in SigmaPlot v11.0. The effects of changing G max on the number of spikes and synchrony among TEA-'treated' model neurons were tested with paired t-tests. Analyses changes in spike number and spike synchrony with changing coupling and synaptic strengths ( Figure 5) were analyzed with One-Way ANOVAs with post-hoc pairwise comparisons between a given percent change and the TEA case conducted via Holm-Sidak tests.