Functional connectivity in a rhythmic inhibitory circuit using Granger causality
© Kispersky et al; licensee BioMed Central Ltd. 2011
Received: 17 March 2011
Accepted: 25 May 2011
Published: 25 May 2011
Skip to main content
© Kispersky et al; licensee BioMed Central Ltd. 2011
Received: 17 March 2011
Accepted: 25 May 2011
Published: 25 May 2011
Understanding circuit function would be greatly facilitated by methods that allow the simultaneous estimation of the functional strengths of all of the synapses in the network during ongoing network activity. Towards that end, we used Granger causality analysis on electrical recordings from the pyloric network of the crab Cancer borealis, a small rhythmic circuit with known connectivity, and known neuronal intrinsic properties.
Granger causality analysis reported a causal relationship where there is no anatomical correlate because of the strong oscillatory behavior of the pyloric circuit. Additionally, we failed to find a direct relationship between synaptic strength and Granger causality in a set of pyloric circuit models.
We conclude that the lack of a relationship between synaptic strength and functional connectivity occurs because Granger causality essentially collapses the direct contribution of the synapse with the intrinsic properties of the postsynaptic neuron. We suggest that the richness of the dynamical properties of most biological neurons complicates the simple interpretation of the results of functional connectivity analyses using Granger causality.
The goal of much of neuroscience is to understand how circuit dynamics arise from the properties of the circuit neurons and their connectivity . Unfortunately, the number of circuits whose anatomical connectivity is well-established is relatively low. Additionally, even when available, the connectivity diagram only provides a static description of synaptic connectivity, while in reality synaptic strength varies as a function of neuromodulation [2–5] and time-dependent processes such as facilitation and depression [1, 6–8].
In most circuits, the anatomical connectivity is not completely known, and functional connectivity techniques are often used to give a directed, dynamically changing assessment of network interaction [9–14]. Functional connectivity is defined as the strength of the influence one network member has on another during ongoing circuit behavior [15, 16]. Such a metric is desirable because it provides a tool to describe circuits in terms of their functional interactions rather than only the presence or absence of synaptic connectivity. Additionally, functional connectivity metrics are desirable because unlike static, anatomical descriptions, they can be used to measure changes in how network elements influence one another when repeated over time and over different biological states. In the present work, we apply functional connectivity analyses to a biological circuit, the pyloric network of the crustacean stomatogastric ganglion (STG) in which the anatomical connectivity is known  to study the relationship between these two network descriptors.
One method of functional connectivity analysis that is widely used in neuroscience is Granger causality (GC), a technique from the field of economics [18, 19]. GC is based on the idea that, given two signals X and Y, if knowing the past of Y is useful for predicting the future of X then Y must have a causal influence on X. Because GC is capable of estimating directed interactions between multiple time series, it can be a useful tool for determining the functional strength and direction of connections between members of a neural network [10, 19, 20]. Granger causality has been applied to neural networks of diverse sizes and used to measure dynamically changing connectivity in functional magnetic resonance imaging (fMRI) learning paradigms , the interaction of closely related nuclei in a simulated hippocampus during environmental exploration , and the plastic changes in the connectivity of single neurons in a cortical cell culture [23, 24]. With few exceptions [25–27], however, functional network estimations made with GC have not been directly compared with underlying anatomical connectivity because the connection pattern is either not known or highly complex.
In the present work we apply GC analysis to electrophysiological recordings and model data from the STG to study the relationship between functional and anatomical connectivity more rigorously. The STG has been a model system for studies of neuromodulation and motor pattern generation for many years [28–31]. One major advantage of the pyloric circuit is its small size, which allowed investigators to map its connectivity [17, 32]. While the anatomical connectivity of the circuit is known, less is known about how each synaptic connection contributes to the functional output of the system, under many different modulatory conditions. This is complicated by the fact that every synapse in the pyloric circuit can be modulated by one or more amines and neuropeptides [3, 28, 33, 34] and that they show synaptic depression . Classically, synaptic strength between two neurons can be well-studied using electrophysiological methods when the network's activity is silenced [35, 36], but it is more difficult to measure synaptic strength during ongoing network activity. Consequently, the effects of a neuromodulator can be effectively studied on single synapses [34, 37], but it can be difficult to assess the functional importance of those changes when the network is active, and each neuron is influenced by many others in the circuit. For these reasons, we were interested to ask whether GC would be a useful tool with which to study functional connectivity during ongoing network activity.
Granger values numerically represent how much the prediction of future values of a particular trace is improved by the inclusion of the past values of an additional trace. Mathematically, GC values are computed by taking the logarithm of the ratio of the prediction error when not considering the additional trace to the prediction error when considering the additional trace. This ratio is always greater than one (error can only decrease when considering more information) and thus the logarithm of this ratio, Granger causality, is always greater than zero. A Granger value of zero would represent a prediction error ratio of 1 (for example, no error reduction with the inclusion of more information, e0 = 1) and a Granger value of 1 would represent a prediction error ratio of about 3 (for example, a 3 fold reduction in error, e1 = 2.71). Typical Granger values from biological data sets were around 0.2 (e0.2 = 1.2, prediction error reduced by about 20%) and for highly correlated traces (see below) around 1.0.
We prepared data sets as described in the Methods by computing a smooth rate function from spike-detected extra-cellular recordings and then applying the GC algorithm to each of these sets of recordings (Figure 2, color code as in Figure 1; in this paper, solid lines with diamond tips represent statistically significant Granger interactions). Because there was heterogeneity in the traces in terms of period length, number of spikes per burst and burst duration (Figure 2), we allowed the total data length analyzed to vary, but ensured that in all cases the total data length was at least 3 times the period. Further, we mandated that results were qualitatively insensitive to small changes in the data length. In this context, we defined parameter insensitivity to mean equivalent statistically significant Granger relationships when data length was changed by one period. Excepting total data length, all other parameters were kept constant between different data sets. Important among the constant parameters was the number of past time points that the GC algorithm would consider in constructing each auto-regressive model (the model order). If this time window was allowed to be too long, the auto-regressive model derived from a single trace would become highly effective at predicting its own future. This is because after an entire period has elapsed in an oscillator, a highly accurate auto-regressive model can be constructed. In this situation, additional external data sources can no longer aid in producing a more accurate model and thus are not considered to be Granger causal. For all biological data, we fixed the time lag parameter at approximately 40% of a pyloric period, or about 400 ms in keeping with the general principle of selecting model orders that adequately capture the timescale of the interactions of interest . This time window captures the slow oscillation and resulting graded neurotransmitter release seen between cells during a pyloric cycle [44, 45]. Importantly, graded neurotransmitter release (and thus causal interaction) due to the slow voltage oscillation is known to be functionally relevant in the pyloric circuit [45, 46].
Statistically significant functional connections were found with GC from the PD to the LP neurons and from the LP to the PY neurons, both connections that are known to have an anatomical analog in biological preparations (Figure 1). Notably, GC operated successfully on these data even though all neurons are coupled by inhibition. In all four data sets (Figure 2A-D), there is a Granger causal relationship from the PY to the PD neuron. While it is true that PY neuron activity always precedes PD neuron activity in the temporal progression of the pyloric rhythm, the PY neurons do not directly inhibit the PD neurons and thus a functional connection cannot be related to direct neurotransmitter release from PY to PD. Importantly, we had selected a model order of 400 ms, a window of time too short to allow for GC to report functional interactions from PY to PD related to the disynaptic connection via the LP neuron. The consistency of the functional connection between the PY and PD neuron between datasets led us to hypothesize that in data sets with strong oscillations, where temporally structured relationships between neurons may not necessarily be due to direct synaptic interactions, a fundamental confound may emerge in interpreting functional connectivity computed with Granger causality.
Given the functional connection between two neurons that was not attributable to a known synaptic connection, we set out to account for this unexpected finding by analyzing a variety of artificial data (with known properties), to gain intuition for what classes of data sets lend themselves well to GC and under what conditions the analysis could break down.
The correlated noise traces were then used as inputs to the GC algorithm to test how GC would change given different levels of correlation and anti-correlation between the two traces. In the case of zero correlation (that is, independent noise traces) the computed value of GC did not meet the requirements for significance and thus we concluded that there was no GC relationship as expected (Figure 3C, network diagram, GC = 0.048, P = 0.33, in this paper, lines are drawn as dashed to indicate the lack of statistical significance). When the correlation value was changed to either -20% or +20% significant GC values were obtained (Figure 3B, c = -0.2, GC = 0.108, P = 1.22 × 10-6; Figure 3D, c = 0.2, GC = 0.089, P = 1.63 × 10-4). Finally, when correlation values were raised to c = -80% and c = +80% highly significant GC values were obtained (Figure 3A, c = -0.8, GC = 1.056, P < 1 × 10-6; Figure 3E, c = 0.8, GC = 1.073, P < 1 × 10-6). Diamond-tipped lines are drawn with a thickness proportional to the strength of the GC relationship between S1 and S2. Thus, GC analysis detected directed relationships in correlated noise and the computed GC values were larger for traces that were more strongly correlated. Furthermore, GC analysis was equally effective at detecting functional relationships between negatively correlated traces (an 'inhibitory interaction') and positively correlated traces (an 'excitatory interaction').
When analyzing data of this nature an advantage of GC over cross-correlation analysis is that GC is less sensitive to variability especially at lower levels of correlation. While the peaks of the cross-correlation functions matched the predetermined values of correlation in the examples shown, in general the cross-correlation curves tended to have a large amount of variability. This is evident in the cross correlation curve of the two fully independent noise traces (Figure 3C, top right plot) which reports correlation values as high as -50% attributable to chance. If cross-correlation functions from repeated trials were averaged, the peak in the averaged function would converge upon the input value (data not shown). In contrast, GC analysis reports a highly significant relationship between the two traces correlated at only ± 20% (Figure 3B and 3D, network diagram) from a single trial. We hypothesize that this is due to the fact that GC considers a window of data over the entire data trace (thus implicitly averaging) while cross-correlation is a point-wise measure and thus sensitive to small fluctuations.
A disadvantage of GC is that it does not report the time lag between the two input waveforms, information that is reconstructed by sequentially taking the cross-correlation over the entire range of time lags. Time lag information could be extracted from GC analysis by looking at the coefficients of the auto-regressive model. In these examples, the largest coefficients (the most predictive past values of each trace) would fall at the artificially introduced time lag of 200 ms.
As expected, the effect of this coupling pattern was that spikes in each cell evoked hyperpolarizations of a few mVs in the respective postsynaptic neuron, as shown in spike-triggered averages (Figure 4B, traces have same grayscale as in A). This indicates that, as expected, spikes hyperpolarized the membrane voltage of target cells and thus could be considered to be causally influencing the spike times of their post-synaptic partner cells.
Spike trains were first turned into binary spike trains and then convolved with a half-Gaussian (see Methods) to produce a smooth rate function (as was done previously with the experimental data). In these data, functional connectivity computed with GC matched the synaptic structure of the network (Figure 4C). Although each cell was firing slowly and with variability, the GC relationship reported between each cell pair was highly significant (Figure 4D, Cell 1 to 2: GC = 0.012, P = 8.11 × 10-5; Cell 2 to 3: GC = 0.017, P = 3.98 × 10-7; Cell 3 to 1: GC = 0.011, P = 3.82 × 10-4).
Notably, because the input data are based upon an approximation of a rate function, they contain none of the subthreshold voltage deflections directly caused by inhibition. Thus, causal interactions are being calculated purely based on the timing of spikes. This is significant because it suggests that the analysis is insensitive to whether inhibition or excitation mediates interactions between neurons as long as an appropriate model order is chosen. Indeed, when the same cyclic network was simulated with excitatory synapses, GC analysis produced equivalent output (data not shown). When using excitatory coupling instead of inhibition, the same connections were found to have statistically significant Granger values (values were: Cell 1→2: 0.012, Cell 2→3: 0.017, Cell 3→1: 0.011, P < 0.01 for all three numbers).
We also wondered what functional relationships would exist if the model network was not connected in a cyclical fashion but instead connected analogously to the pyloric circuit (Figure 1). This result was of interest because we wanted to know whether the functional interactions we observed in biological data (Figure 2), that matched anatomical synaptic connections (PD/LP and LP/PY), were actually based on the underlying synapses or whether those functional interactions were also a result of fixed temporal relationships due to oscillations in the spike trains (as we hypothesized for PY to PD). To answer this, we reasoned that in a network coupled analogously to the pyloric circuit, but in which cells fired with variability (that is, without oscillations), there should be no functional interaction between the PY and the PD cell while the PD/LP and the LP/PY interactions ought to remain. Thus, we coupled the same model cells (Cells 1, 2 and 3) as in the pyloric circuit in which the PD/LP pair and the LP/PY pair are connected with reciprocal inhibition while PD inhibits PY without a return synapse from PY to PD. For all synapses we used synaptic maximal conductances of 0.1 nS (Figure 4E). Statistically significant GC relationships were found between the cell 1/2 pair (analogous to PD/LP) and the cell 2/3 pair (analogous to LP/PY) all of which match the synaptic architecture (Figure 4F, GC and P values given in the figure). Additionally, a functional connection from cell 1 to cell 3 was present but fell slightly above the threshold for significance (Figure 4F, cell 1 to 3, GC = 0.002, P = 0.06, dashed line indicates lack of significance). Importantly, the GC value for the PY to PD synapse was far from achieving statistical significance (GC = 0.001, P = 0.33). Thus, we conclude that due to the lack of oscillatory activity in this model, functional connections matched synaptic architecture even when model cells were coupled analogously to the pyloric circuit. This finding reinforces the hypothesis that functional connections unrelated to synaptic connections seen in the biological data (Figure 2) are related to the strongly oscillatory behavior of the AB/PD neurons.
To investigate the relationship between synaptic conductance magnitude and the strength of the functional relationships we applied GC analysis to a set of 84 pyloric model networks with highly variable intrinsic and synaptic properties. These 84 networks are a subset of a much larger pool of model networks  that had a wide range of known parameters in their synaptic conductances and membrane conductances. To further increase the variability of synaptic and intrinsic conductance values in the selected model networks we introduced a kinetic temperature dependency into each synapse and ion channel and then simulated each network at different temperatures (see Methods). Therefore, a large variety of synaptic conductances and channel time constants was represented by the selected model networks.
As part of the GC calculations, we computed GC values for the PY to PD synapses in all 84 networks. It is known that this synaptic connection does not exist in biological preparations (Figure 1) and thus, by design, gsyn for the model PY to PD synapse was zero in all model networks considered in this analysis. Nonetheless, many networks had GC values for the PY to PD interaction that were comparable to values computed for the other synapses (Figure 5F, PY to PD).
Previously, when considering biological data (Figure 1) we attributed functional connections unrelated to synaptic connectivity to the highly oscillatory nature of pyloric activity. We hypothesized that a functional relationship existed due to the fixed temporal relationship between the PY and PD neurons. To test whether this was the case in the set of model pyloric oscillators we subdivided the 84 networks into two groups, those that displayed fully pyloric behavior (triphasic, correct burst order, with appropriate phase relationships and burst durations ) and those that were pyloric-like (triphasic but with abnormally long burst durations, tonic firing instead of bursting, or phase relationships not seen in normal biological data). The mean GC values of the two populations were significantly different from one another (Figure 5, bottom right panel, P < 0.01, Kruskal-Wallis test, npyloric = 29, npyloric-like = 55). GC values for the PY to PD synapse were significantly higher when regular oscillations are present (Figure 5, black dots, pyloric networks only) as compared with more diverse firing behavior (Figure 5, gray dots, pyloric-like networks only) providing support for the conclusion that oscillations can lead to functional connectivity unrelated to synaptic connectivity.
While it is well appreciated that functional connectivity methods are not equivalent to direct anatomical and physiological measures of connectivity, it is clear that the two are related . The nature of the relationship is difficult to determine because the strength of a biological synaptic connection will be altered by the synapse's history of activation and neuromodulation, and the ability of a synaptic input to influence the firing of its target neuron will be altered by changes in the intrinsic postsynaptic excitability. Reliable functional connectivity metrics would be extremely useful because they can be used during ongoing circuit performance to assess how important each synapse is to the circuit's output under different conditions. Here we demonstrate that Granger causality, a measure of functional connectivity, can be used with synaptic inhibition, unlike previous claims to the contrary . Further, we point out a serious limitation in interpreting results of Granger methods when there are strongly oscillatory components in the data. The ability of an input to influence a target neuron's activity depends both on the input strength and the target neuron's excitability. These two factors are effectively collapsed in Granger calculations leading to ambiguity about the mechanism underlying a changed functional connection.
We used the pyloric network of the stomatogastric ganglion as the test circuit for these studies, precisely because the mechanisms underlying the generation of the pyloric rhythm are well understood . Rhythm generation in the pyloric circuit depends on a strongly oscillatory pacemaker kernel of three electrically coupled neurons (the AB and two PD neurons) , and the specific pattern of firing depends on the postinhibitory rebound properties of the LP and PY neurons [41, 42], as all of the synaptic connections among the pyloric circuit neurons are inhibitory (Figure 1).
Based on GC analysis, a functional relationship exists from the PY neuron to the PD neuron in biological extracellular recordings, an interaction which is not based on a synaptic connection in biological networks. We argue that this functional connection arises from the highly regular oscillations displayed by the circuit that causes the PY and PD neurons to have a fixed, repeating temporal relationship that is not related to their direct interaction. In other words, in each data set, the PD neurons, because of their intrinsic oscillation period, always fired with a predictable latency after the PY neurons, although the onset of PD activity was in no way causally related to the PY neuron activity. We provided support for this hypothesis by generating a model network that had pyloric-like connectivity but in which the neurons fired with variability instead of in a regularly oscillating pattern (Figure 4E, F). In this case, functional relationships matched the synaptic architecture closely. The interpretation of our results that oscillations introduce functional relationships unrelated to synaptic architecture is further strengthened by the higher Granger values computed for the PY to PD synapse in regularly oscillating pyloric networks when compared to networks that have pyloric-like activity, but are not reliably oscillatory (Figure 5).
Given the prevalence of oscillations in the brain [50–53] and the confound posed by these oscillations for interpreting GC analyses, care should be taken when applying functional connectivity metrics to neuronal data. Neuronal oscillations could complicate the interpretation of data gathered with electroencephalography (EEG)  or other methods [22, 54]. Presumably, the stronger and more regular the oscillation, the more serious this potential confound is likely to be. Certainly, in many brain areas oscillations are less regular and weaker than those found in the pyloric rhythm. Nonetheless, given the widespread use and appeal of GC methods, it is important to recognize that oscillations in the data could lead to functional connectivity unrelated to synaptic interactions.
Granger causality is one method in a family of functional network estimation methods each with its own advantages and drawbacks. GC, which is based on auto-regressive modeling [18, 19], can also be computed in the spectral domain where one GC value is computed for each frequency [55, 56]. This can be useful for separating the interactions of signals that contain multiple frequency components. Additionally, the concept of partial Granger causality provides methods for removing the effects of external inputs common to all network members or unmeasured network elements . Throughout this study, we compared our calculations with those returned by partial Granger causality and found no difference in the resulting values. We selected standard GC for this study because of the maturity of the existing tools  and the large set of previous studies that have employed the method. Other related methods include dynamic causal modeling that relies on Baysian inference to measure the causal interactions between signals  or Transfer Entropy that relies on mutual information to compute interactions between network members . We suspect that many measures of functional or effective connectivity will share the same potential interpretation issues as those outlined here. These issues in no way preclude the utility of these methods for understanding how different brain regions interact, as long as mechanistic conclusions are not made from them.
Inhibition is critically important in all brain regions and all functional brain circuits. We show here that GC works equivalently with inhibition and excitation. Importantly, if one uses test model neurons that are relatively silent unless driven by excitation and have no postinhibitory rebound, then a functional relationship may not appear because the postsynaptic neurons are not firing. However, if the postsynaptic neuron is either tonically active unless inhibited, or has postinhibitory rebound , then the change in the postsynaptic neuron's firing will form the basis for a functional interaction. Given that postinhibitory rebounds can be either relatively short or quite long , care should be taken to ensure that an appropriate model order is used to capture potential rebound firing.
In previous reports GC was found to be well-correlated with the amplitude of the synaptic potential . In contrast, we found little correlation between GC and synaptic strength in a population of neuronal models (Figure 5). This result may, at first, seem puzzling. Nonetheless, it is important to remember that the spike timing of a postsynaptic neuron will always be an interaction between the strength of the input synapse and the intrinsic properties of the follower neuron . The population of models studied in Figure 5 contains model pyloric networks with highly variable synaptic and intrinsic conductances . Therefore, it is likely that the variation of synaptic strength is compensated by variations in the intrinsic membrane properties of the follower neurons [36, 49, 62], as synaptic and intrinsic parameters can compensate for each other in a functional network . We conclude that because Granger causality represents the total effect of one cell on another, it effectively collapses the relative contribution of synaptic strength with any postsynaptic modulation due to intrinsic properties. Thus, relating Granger causality directly to synaptic strength is not possible as has been previously noted , even if the two numbers are related.
While we did not observe a relationship between pyloric synaptic strength and functional coupling strength (Figure 5) we did, however, observe such a relationship in the correlated noise data (Figure 3) and spiking model data (data not shown). This is consistent with the interpretation that the complex membrane properties of the follower neurons in the pyloric networks account for the lack of correspondence between synaptic conductance and GC values. Because many real biological neurons in the brain express complex membrane properties [64–68] that are not well-captured by simple integrate and fire or rate models (such as those used in many tests of the functional connectivity algorithms), this will potentially complicate the simple interpretation of GC for biological data. In particular, changes in neuromodulatory tone across the brain may alter both target neuron excitability and the strength of the inputs to the target. Thus, Granger causality may show that two brain regions are more 'effectively connected' when the strength of the input itself may be unchanged or even decreased.
We also found that functional interactions computed from spike timing information alone could mirror synaptic architecture. In our spiking model data the inhibitory interactions of the neurons caused hyperpolarizing deflections in post-synaptic membrane voltages (Figure 4) but this information was filtered out when spike trains were smoothed before GC analysis leaving only the relative timing of spikes. Functional relationships based on spike timing information matched synaptic architecture even when model synaptic connectivity was more complicated as in the pyloric-like model network (Figure 4E). In this network, some synaptic interactions are reciprocal and single cells receive multiple inputs that can lead to complicated interactions. Computing functional interactions that match network topology from spike timing information alone is non-trivial and speaks to the utility of the Granger method. Functional connection patterns based on spike times that closely match network topologies could also be useful for in vivo data from behaving animals that, similar to our recordings in the pyloric circuit, rely heavily on extracellular recordings of spike times. Spike times were recently shown to be sufficient to fully mathematically infer synaptic weights in a small network of integrate and fire neurons [26, 69].
We investigated the utility of Granger causality analysis for understanding functional interactions in a rhythmic, inhibitory circuit with a known connectivity diagram. We found that oscillations in the input data can yield functional connections that are not predicated on biological interactions. Additionally, we found that GC functions well in the context of inhibitory coupling and that when postsynaptic neurons have complex and realistic intrinsic membrane properties, the relationship between GC results and actual synaptic strength becomes more complicated.
where parameter values are a = 0.02, b = -0.1, c = -55, d = 6, I app = 4.65 where I app represents the applied current. Synapses are modeled as double exponential waveforms  with τ rise = 0.7 and τ fall = 5.8.
Simulation of STG networks followed established approaches . All model neurons were selected from an existing database of STG cells . Twelve total networks were selected such that the output voltage trace displayed either fully pyloric behavior or pyloric-like behavior. Fully pyloric behavior was defined by phase, duty cycle and relative burst times of the three cells that matched in vitro recordings . Pyloric-like networks had burst frequencies, spike rates or phase relationships that fell outside of the normal range of behaviors observed in biological preparations. Synaptic strengths between model cells could take on values of 0, 3, 10, 30, 100, or 200 nS. Each network was run over a range of temperatures (7, 11, 15, 19, 23, 27, or 31°C) to yield a total of 84 different networks. Each ionic conductance and each synapse was subject to kinetic modification based on the temperature of the simulation . Conductances were updated based on the equation where Q10 = 1.5, R2 is the updated conductance, R1 was the existing value of the conductance, T2 was the simulation temperature and T1 was the base temperature set to 11°C. This equation was solved for R2 to yield . All activation and inactivation rates were given Q10 = 2. Logarithms are natural logarithms (base e). After simulation, voltage traces were preprocessed and analyzed with Granger causality in the same manner as all other data sets.
Cancer borealis were purchased from Commercial Lobster (Boston, MA, USA) and maintained in artificial seawater tanks at 10 - 12°C. Before dissection, animals were cold-anesthetized by packing them in ice for 30·min. The stomach was removed from the crab and placed in chilled physiological saline NaCl, 440 mM; KCl, 13 mM; MgCl2, 26 mM; CaCl2, 13 mM; Trizma base, 11 mM; maleic acid, 5 mM; pH·7.45) while the stomatogastric nervous system (STNS) was dissected and pinned to a Sylgard-coated petri dish containing saline . The stomatogastric ganglion (STG) was desheathed to allow for intracellular recording of the cell bodies.
Extracellular recordings were obtained by making petroleum jelly wells around branching nerves of the STNS and placing a stainless steel electrode inside the well and one in the bath. The extracellular traces were amplified and recorded with an AM Systems Model 1700 AC Amplifier. During recordings, preparations were continuously superfused with chilled physiological saline (12°C) by means of a Peltier cooling system (Warner Instruments/Harvard Apparatus Hamden, CT) and the temperature was monitored using a thermoelectric probe in the bath.
All traces were recorded by a PC running pClamp software (Version 10.2). Spike times were extracted from the files using Spike 2 software (version 6.04) and a thresholding script. Data files were converted to a MATLAB compatible format for further processing.
GC analysis, as any data analysis method, makes several assumptions about the data upon which it operates. We ensured that for all data sets analyzed, both biological and model, these assumptions were not violated. First, GC requires that input data are covariance stationary meaning that the mean and variance cannot change over time. The pyloric circuit is a canonical example of covariance stationarity in biological data. Second, the auto-regressive models that are constructed and form the basis of GC analysis must represent the data well. For all data we verified that the root mean square error in the models did not exceed 5% and in many cases model fits were significantly better.
Important parameter choices were the total length of the input data and the model order (the total number of past observations used to construct the auto-regressive model). Total input data length was varied between data sets but constrained to always include at least three oscillation periods. Excepting total data length analyzed, all other parameters were kept constant between different data sets. Model order was selected manually to best reflect the time scale of functional interactions in the data. For simulations of spiking neurons, for example, model order was chosen to span 100 ms, to capture the kinetics of the synaptic waveform. For pyloric data, a model order of 400 ms was used to capture a significant portion of the slow oscillation. Statistics were computed using the Kruskal-Wallis test. Plots were produced in MATLAB and then finished in the Inkscape (http://inkscape.org) graphics program.
anterior bursting neuron
coefficient of variation
functional magnetic resonance imaging
lateral pyloric neuron
lateral ventricular nerve
standard error of the mean
stomatogastric nervous system
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.