Abstract
Neurons and networks undergo a process of homeostatic plasticity that stabilizes output by integrating activity levels with network and cellular properties to counter longer-term perturbations. Here we describe a rapid compensatory interaction among a pair of potassium currents, IA and IKCa, that stabilizes both intrinsic excitability and network function in the cardiac ganglion of the crab, Cancer borealis. We determined that mRNA levels in single identified neurons for the channels which encode IA and IKCa are positively correlated, yet the ionic currents themselves are negatively correlated, across a population of motor neurons. We then determined that these currents are functionally coupled; decreasing levels of either current within a neuron causes a rapid increase in the other. This functional interdependence results in homeostatic stabilization of both the individual neuronal and the network output. Furthermore, these compensatory increases are mechanistically independent, suggesting robustness in the maintenance of neural network output that is critical for survival. Together, we generate a complete model for homeostatic plasticity from mRNA to network output where rapid post-translational compensatory mechanisms acting on a reservoir of channels proteins regulated at the level of gene expression provide homeostatic stabilization of both cellular and network activity.
Introduction
The balance of plasticity and stability in generating appropriate output is a matter of fundamental importance in the nervous system across all functional levels. These processes occur even at the most fundamental level, the excitability of individual neurons, and yet little is known about mechanisms governing these processes (Marder, 2011; Turrigiano, 2011). Early work identified such processes of “homeostatic plasticity” of intrinsic excitability (LeMasson et al., 1993; Turrigiano et al., 1994, 1995; Golowasch et al., 1999), but subsequent focus more intensely shifted to determining how stabilization of synapses is accomplished through synaptic scaling (Turrigiano, 2012). Recently, a resurgence of interest in plasticity of intrinsic excitability has accompanied work on synaptic scaling (Debanne and Poo, 2010; Misonou, 2010; Turrigiano, 2011).
Fewer studies on homeostatic plasticity have considered functional compensation in the context of endogenous network activity. The most dramatic example may be complete recovery of motor network output following loss of central inputs as a result of changes in conductances in the crustacean stomatogastric ganglion (STG) (Thoby-Brisson and Simmers, 1998, 2002). Knock-outs of K+ channels in mice have been shown to have modest effects on phenotype and cellular output as a result of compensation by other K+ channels (Guo et al., 2005; Nerbonne et al., 2008). Additionally, overexpression of A-type K+ channels in STG neurons results in little change in neuronal output as a result of compensatory increases in H-current (MacLean et al., 2003, 2005). However, these examples feature mechanisms that act over longer time scales of days to weeks. While initial reports of plasticity in intrinsic excitability were found over shorter time scales (Desai et al., 1999; Golowasch et al., 1999), surprisingly little is known of the role these mechanisms may play in short-term ongoing activity of biologically intact networks, specifically where an expectation for rapid conservation of output could be argued, such as in central pattern generators (CPGs).
Compensation in CPG circuits may be inferred from the fact that normal populations of unmanipulated motor neurons of two invertebrate CPGs, the cardiac and stomatogastric ganglia, show correlations in expression levels of mRNAs for ion channels (Schulz et al., 2007; Tobin et al., 2009) and membrane conductances (Khorkova and Golowasch, 2007; Temporal et al., 2012). One relationship detected in previous work (Tobin et al., 2009) is a positive correlation between BKKCa and shaker mRNA levels in neurons of the cardiac ganglion. These channels encode calcium-activated and A-type K+ currents, respectively. However, it is unclear how a neuron would use two similar hyperpolarizing conductances additively to generate or maintain its output. In this study we focused on elucidating the functional relationship between these two conductances in motor neurons of the crab cardiac ganglion. We discovered a striking discrepancy in the relationship of these channels across functional levels: mRNAs for these channels were positively correlated, while their conductances were negatively correlated. Therefore, we propose a comprehensive hypothesis for plasticity of excitability from mRNA to network output whereby rapid compensation provides stabilization of cellular and network activity.
Materials and Methods
Preparations.
Cancer borealis crabs of either sex were purchased and shipped overnight from The Fresh Lobster Company (Gloucester, MA). Crabs were kept between 24 h and 2 weeks in artificial sea water at 12°C before use. Crabs were anesthetized in ice for 15 min before the dissection. The dissection took place in chilled physiological saline comprised of 440 mm NaCl, 26 mm MgCl2 13 mm CaCl2, 11.2 mm Trizma base, 11 mm KCl, and 5 mm maleic acid (pH = 7.4). When we wanted to isolate individual large cells, individual strands of bulking nylon were used to ligate the nerve on both sides of a large cell soma. To impale large cells, each cell was individually desheathed using a tungsten needle (Fine Science Tools).
Quantitative single-cell RT-PCR.
Quantitative RT-PCR was performed as previously described (Schulz et al., 2006a; Tobin et al., 2009). Primers specific for real-time PCR detection of shal, BKKCa, shab, and shaker using Sybr Green were developed and designed using Primer3 software and are the same as previously reported (Schulz et al., 2006a; Tobin et al., 2009). Briefly, total RNA was isolated using RNeasy micro column-based RNA extraction kit (Qiagen), reverse transcribed using SuperScript III reverse transcriptase (Invitrogen), and used as a template in real-time RT-PCR with Sybr Green (SABiosciences) in a RotorGene 3000 real-time PCR machine (Corbett Research). Previous studies have determined that in LC motor neurons, correlations can be equally well detected among channel mRNA levels with and without normalization of real-time results to 18S rRNA (Tobin et al., 2009). Values reported here are total copy numbers from a single neuron, and are not normalized with respect to 18S levels.
Pharmacology.
Pharmacological blockers were dissolved in physiological saline and perfused onto the cardiac ganglion using a Rabbit peristaltic pump (Rainin Instruments) at a rate of 1.5 ml/min or added to the preparation from a stock solution via pipette. The following pharmacological agents were used: tetraethylammonium dissolved in saline at 25 mm, 4-aminopyridine dissolved in saline at 1 mm, cadmium chloride dissolved in saline at 250 μm (Acros Organics), tetrodotoxin dissolved in saline at 1 μm (Alomone Laboratories), BAPTA-AM dissolved DMSO and applied at 30 μm in saline, ryanodine dissolved DMSO and applied at 100 μm in saline, staurosporine dissolved in DMSO and applied at 5 μm in saline, okadaic acid dissolved in DMSO and applied at 500 nm in saline (Ascent Scientific), cyclosporine A dissolved in DMSO and applied at 2 μm in saline (Tocris Biosciences). All DMSO applications resulted in a final concentration of DMSO that was <1% (range: 0.000025% to 0.5%).
Pharmacological agents used to investigate intracellular mechanisms involved in the compensatory response (BAPTA-AM, ryanodine, okadaic acid, cyclosporine A, staurosporine) were applied to the cardiac ganglion 1 h prior (2 h prior with ryanodine) to the application of the blocker which caused the compensation (TEA or 4AP). A cell or preparation was exposed only to one channel blocker type (TEA or 4AP) for a given experiment.
Current measurements.
All experiments were performed in physiological saline cooled to 12°C. To measure current magnitudes and activation properties, two-electrode voltage-clamp (TEVC) experiments were performed by impaling a large cell with two glass electrodes filled with 3 m KCl (8–17 MΩ resistance) and an Axoclamp 2A amplifier (Molecular Devices). All recordings were made from anterior large cell somata; action potential conductances were blocked (unless noted otherwise) by tightening thread ligatures on both sides of the large cell soma, preserving space clamp. TEVC protocols were created, driven and recorded with clampex 9.2 software (Molecular Devices). Current recordings were analyzed with Clampfit 9.2 software (Molecular Devices). Current and voltage traces were sometimes filtered with a lowpass boxcar filter using 7 smoothing points. Most voltage clamps were modified from those used previously in STG preparations (Golowasch and Marder, 1992; Khorkova and Golowasch, 2007; Temporal et al., 2012). High threshold potassium current (IHTK) magnitude was measured using a leak subtracted TEVC protocol with a holding potential of −40 mV and 16 voltage steps from −55 mV to +20 mV (5 mV intervals). A-type potassium current (IA) magnitude was measured by subtracting the IHTK current traces from a TEVC protocol that is identical except for a holding potential of −80 mV. Calcium-activated potassium current (IKCa) was isolated by subtracting postcadmium (250 μm CdCl2, 1 h) IHTK current traces from precadmium IHTK current traces (isolating the cadmium-sensitive outward current). Delayed rectifier potassium current (IKd) was isolated using the IHTK TEVC protocol after cadmium exposure (250 μm CdCl2, 1 h). All current magnitude measurements were taken at 0 mV on an I-V plot made from the current traces.
Large cell excitability and network output.
Cardiac network output was monitored with a single intracellular recording (using same equipment as TEVC protocols) taken from one of the three anterior large cells and an extracellular differential recording made with a model 1700 A-M Systems AC amplifier and two stainless steel wires; one placed inside and one outside a vaseline well located around the central nerve of the CG (see Fig. 1). Pharmacological blockers (TEA or 4AP) were perfused on the entirety of the CG with the exception of the four small cell and two posterior large cell somata. These cells were isolated from the perfusion by a vaseline well placed around these cells and the posterior branch point containing regular physiological saline (see Fig. 1). Using this experimental configuration we monitored the effect of the LC compensatory response (isolated to the best of our ability from small cells) in the context of the functioning network. Network activity was recorded in 10 min intervals before and during blocker perfusion. These recordings were analyzed using Spike2 v6.00 software (Cambridge Electronic Design).
LC intrinsic excitability was examined under similar conditions except 10−6 m TTX saline was placed in the vaseline well around the four small cell and two posterior LC somata. This eliminated spontaneous network activity and small cell excitatory input into the anterior LCs. Excitability in the anterior LCs was then monitored using two-electrode current-clamp (TECC) protocols run before and every 5 min after TEA or 4AP perfusion. TECC protocol was a six step depolarizing current injection (from 1 to 6 nA) lasting six seconds per step and six seconds between steps.
Statistics.
All statistical tests were performed with SigmaPlot v11.0 (Systat, Aspire Software International). All data were confirmed to be of normal distribution as required by statistical analyses used. Relationships between channel mRNAs and ionic currents were analyzed using Pearson's correlation test, and coefficients of determination were calculated from the resulting correlation coefficients. In the case of Figure 2, C and D, a potential outlier was identified that could be anchoring a false positive for the Pearson's test (see arrows). Analyses on these datasets were performed both with and without the data point in question and both results reported. Bonferroni corrections were used for multiple comparisons in the correlation analyses, and the p value adjusted to 0.017 for statistical significance (three comparisons each for mRNA and ionic currents). Changes in current magnitude before and after pharmacological block were analyzed in one of two ways. Raw currents were analyzed before and after blockade via paired t tests, and these are reported in Figure 3. In Figure 5, significant changes in a current relative to baseline were expressed as a percent change from zero, and analyzed via one-sample t test with the hypothesized population mean set to 0. Overall changes in burst duration, spikes per burst, and spike frequency within the burst reported in Figure 4 were analyzed with repeated-measures ANOVA.
Results
The crustacean cardiac ganglion as a model for central pattern generator network activity
The rhythmic pumping of the heart in decapod crustaceans such as the crab, Cancer borealis (our model organism), is neurogenic in nature, and under the control of a simple central pattern generator network called the cardiac ganglion (Alexandrowicz, 1932) (Fig. 1). The ganglion consists of only nine neurons: four “small cell” interneurons (SCs) that generate the pacemaker activity and five “large cell” motor neurons (LCs) that innervate the heart musculature (Hartline, 1967; Tazaki and Cooke, 1983c). The SCs of the cardiac ganglion are endogenous oscillators, i.e., they undergo spontaneous and rhythmic generation of a depolarizing wave of membrane potential that leads to a bursting phenotype of multiple spikes per burst (Tazaki and Cooke, 1983b; Cooke, 2002). The LCs of the CG produce bursts of action potentials as a result of synaptic pacemaker input from the SCs, ultimately leading to muscle contraction (Hartline, 1967; Tazaki and Cooke, 1979, 1983a,c; Berlind, 1989). The behavioral output of the ganglion represents a direct correlation of the influence of LCs on heart muscle (Sakurai and Wilkens, 2003; García-Crescioni et al., 2010), and thus a direct measure of heart activity. Yet the entire network can be dissected out intact, and maintained in physiological saline for extended periods of recording while continuing to produce its endogenous rhythmic output. The motor neurons are all individually identifiable, and due to their distributed nature within the ganglion (Fig. 1), we can perform pharmacological manipulations on one or multiple motor neurons, either in isolation or in the intact, functional network (Tazaki and Cooke, 1983a; Cooke, 2002). The underlying burst potentials of the LCs represent functional output at the motor neuron level, and our preliminary modeling studies show how the simplicity of this model system can be used to study functional implications of the relationships between mRNAs and ionic conductances on cellular output (Ball et al., 2010; Franklin et al., 2010).
Relationship between IA and IKCa in a population of large cell motor neurons
Potassium currents were measured using protocols developed for the STG cells of the same species (Golowasch and Marder, 1992). The total outward current of LC motor neurons consists primarily of three K+ currents (Fig. 2B) (Golowasch and Marder, 1992): A-type transient K+ current (IA), calcium-activated K+ current (IKCa), and delayed rectifier K+ current (IKd). IKCa and IKd can be found in one combined current trace elicited from holding potentials at −40 mV or higher, and is termed here as the high-threshold K+ current (IHTK), while the A-type current can be measured by subtracting IHTK from the total outward current elicited from a holding potential of −80 mV.
We first examined the relationships among mRNA levels in single identified LC motor neurons for three channel genes that correspond to these three K+ currents (Atkinson et al., 1991; Tsunoda and Salkoff, 1995; Kim et al., 1998): BKKCa (IKCa), shaker (IA), and shab (IKd). We detected a significant correlation between BKKCa and shaker mRNA levels across a population of 20 LC motor neurons (Fig. 2A, left), but no correlations among any other channel mRNAs (Fig. 2A, middle, right). A fourth channel mRNA encoding an A-type K+ current, shal, also was significantly positively correlated to BKKCa mRNA levels (p < 0.005; R2 = 0.58), as well as to shaker (p < 0.002; R2 = 0.42) but not to shab, suggesting an overall relationship between IKCa and IA, but not with IKd, in these cells.
We next examined the relationships among the ionic currents encoded by these channel genes across a population of LC motor neurons. IA, IKCa, and IKd were all measured in each cell across a population of LC motor neurons. A similar pattern of correlated current levels was seen as with the mRNA with one striking distinction: only IA and IKCa showed a significant correlation (Fig. 2A), but the correlation was strongly negative as opposed to positive as seen in the mRNA measurements. No significant correlations were found between IA and IKd, or IKCa and IKd (Fig. 2A). Because of the striking difference between the mRNA correlation (positive), and the current correlation (negative), we were most interested in pursuing the relationship between IA and IKCa, but had concerns about the effects of pharmacological blockers used to measure IKCa (Fig. 3F); we use Cd2+ to block calcium currents that evoke IKCa to measure this current via subtraction from IHTK, as we have not identified any blockers specific to IKCa in our preparation despite many attempts. Therefore, we decided to use the peak IHTK as an indicator of IKCa abundance in these experiments.
IHTK (HighThreshold K+) is known from previous work (Golowasch and Marder, 1992; Haedo and Golowasch, 2006; Khorkova and Golowasch, 2007) to, in large part, consist of IKCa, particularly the transient portion, as IKd shows no transient peak that could account for the peak in IHTK (Fig. 2B). We also determined that measurements of the peak IHTK were likely sufficient to reveal the relationship between IKCa and IA, as the only transient peak seen in the currents that underlie IHTK belongs to IKCa (Fig. 2B, right). Peak IHTK also shows the same negative relationship with IA as does IKCa (Fig. 2C, left), and peak IKCa itself very strongly correlates with peak IHTK (Fig. 2C, middle), but not IKd (Fig. 2C, right). Accordingly, for the remainder of the study, we used IHTK peak current from baseline as a measure of IKCa abundance, to allow for measurement of all K+ currents under the least manipulative pharmacological blocker regime, i.e., without the need to block voltage-gated calcium currents.
IKCa and IA are rapidly upregulated within the same LC motor neuron
Because the focus of the study was on the negative relationship between IA and IKCa, we decided to repeat the initial measurements in a new population of LC motor neurons to confirm the original finding before examining the functional impact of this relationship. In a distinct population of 20 LC motor neurons, we once again found a strong negative relationship between IA and IKCa (as revealed by IHTK measurements; Fig. 3A), confirming our original finding. This strongly suggested a causal relationship between IA and IKCa. Specifically, we hypothesized that if these two currents are functionally interrelated, then a decreased level of IA should result in an increased level of IKCa and vice versa. We tested this hypothesis by blocking one current in this pair and determining whether there was an effect on the magnitude of the other. We measured baseline levels of IA in a given cell, then used tetraethylammonium (TEA) to block the HTK-current for 60 min, and then measured IA again post-TEA. 60 min of exposure to TEA significantly increased peak (p < 0.001; paired t test) A-current (Fig. 3B,C). The converse experiment was performed using 4-aminopyridine (4AP) for 60 min to block IA. The peak (p < 0.003) of IHTK was significantly increased following 60 min of 4-AP blockade (Fig. 3B,C). Additionally, these results appear to be due to an overall increase in conductance, since no change was seen in the voltages of activation of IA and IHTK concomitant with the changes in total current as a result of the block experiments (Fig. 3D). No consistent or significant effects of TEA and 4-AP on neuronal input resistance were observed in these cells (Δ RIN TEA: −1.87 ± 2.5 MΩ, Δ RIN 4-AP: −0.52 ± 1.7 MΩ).
Our initial data strongly implicate IKCa and IA as the key pair of currents involved in this response. Since IHTK is a mixed current comprised predominantly of IKCa and IKd, and both of these currents are blocked by TEA, we determined whether the changes seen in IHTK as a result of 4AP blockade were attributable to just one or to both of these currents. Most convincingly, while 4AP block significantly increases IHTK (Fig. 3C), 1 h of exposure to 4AP had no significant effect on IKd magnitude itself (Fig. 3E, left). Unfortunately, we have not identified a blocker specific to IKCa in our preparation to directly test the effects of blocking this current on changes in IA. However, we can indirectly block IKCa by using CdCl2 to block voltage-gated Ca2+ channels that trigger IKCa. Blocking with Cd2+ completely eliminates IKCa in these cells, and causes a significant increase in IA after 1 h of exposure that is very similar to that seen with TEA (Fig. 3E, right). Together, these data strongly suggest that the relationship between IA and IHTK identified in the experiments (Fig. 3) are the result of a causal relationship between the levels of IA and IKCa in a given motor neuron.
IA and IKCa act in a compensatory fashion to stabilize cellular excitability and network output
The negative correlation we measured between IA and IKCa in the normal population, together with the rapid change in currents seen within a cell after TEA and 4AP block, led us to hypothesize that levels of IA and IKCa may be acting in a compensatory fashion to stabilize both cellular and network output. To determine whether such a functional compensation exists, we followed over time the activity of both isolated LCs, as well as LCs in an intact network, with either 4AP or TEA blockade, and then determined whether the output of the cells and of the network changed over the course of the blockade and the potential compensation.
As seen in Figure 4, A and B, TEA blockade caused a substantial increase in the excitability of the LC motor neuron in the initial 10–20 min. However, over the time course of the measured increase in IA seen in previous experiments (Fig. 2C) there was a compensatory change in the output of the cell and of the network toward the baseline level of activity. In isolated LCs (Fig. 4A), both TEA and 4AP application shifted the cells from a less excitable state to a state characterized by large sustained burst potentials. However, over the time course of compensation, the excitability of isolated cells returned to control levels, losing the large burst potentials and showing firing patterns similar to those before the application of the blocker (Fig. 4A).
In the intact network, TEA often caused a complex multiphase bursting output in treated LCs (Fig. 4B; 10 min recording). Given this somewhat complex pattern of bursting, we determined that measurements of burst duration across a population of these cells represented the clearest means to quantify changes in output over time. Over the course of 60–90 min, TEA blockade resulted in a significant change in burst duration (p < 0.01; ANOVA; Fig. 4C, left), initially increasing to nearly twice as long as control before re-establishing a stable burst duration indistinguishable from that at the control level (Fig. 4B).
In the companion experiment with 4AP blockade in the intact network, we see a similar effect of compensation as in the LC motor neuron output (seen in Fig. 4B,C). 4AP block does not cause the same change in burst duration as with TEA. Rather, 4AP blockade initially causes a significant increase in spike frequency and the number of spikes fired per burst (p < 0.01; ANOVA). Over the time course of the blockade, this spike frequency and number of spikes per burst returns to control levels (Fig. 4D). These data indicate that the changes seen in the level of IA and IKCa as a result of TEA and 4AP block, respectively, are a functional compensatory response that stabilizes the output of the cells and of the network as a whole. However, the roles of IA and IKCa are not simply functionally redundant; blocking either current causes an initial change in the output of the cell that is not the same. These results implicate distinct roles for these two currents in addition to their ability to partially compensate for one another. However, we do not rule out the possibility that there could be other currents or mechanisms involved in the compensatory effect.
Compensatory changes in IA and IKCa follow distinct regulatory pathways
We then set out to determine the underlying mechanisms implementing the compensation in these currents. We first investigated whether the change in IKCa elicited by 4AP blockade was dependent on intracellular calcium. To check this, we co-applied the calcium chelator BAPTA with the blocker. When BAPTA was co-applied with 4AP, there was no change in IHTK, in contrast to the characteristically significant increase in IHTK with 4AP blockade (Fig. 5A). We also confirmed that this effect was not simply due to an effect of BAPTA on the calcium-dependence of IKCa; BAPTA alone did not result in a significant reduction of IHTK (data not shown; see also Turrigiano et al., 1994). More specifically, the calcium dependence of the increase in IHTK was attributable to the release from intracellular calcium stores since co-application of 4AP and the intracellular calcium release blocker ryanodine also prevented the compensatory increase in IHTK normally seen with 4AP blockade (Fig. 5A). These results suggest that the compensatory increase in IKCa is dependent on intracellular calcium signaling mechanisms that depend on the release of intracellular calcium stores.
Because relatively rapid changes in current magnitude can often be attributed to changes in phosphorylation states of ion channels and associated regulatory proteins, we examined the effects of phosphatase and kinase inhibitors on the ability of the blockers to elicit the compensatory changes in these currents. Okadaic acid is a broad spectrum inhibitor of serine/threonine protein phosphatases (Cohen et al., 1990). While application of okadaic acid alone did not affect baseline levels of IHTK (data not shown), co-application of okadaic acid with 4AP also abolished the compensatory increase in IHTK (Fig. 5A), suggesting the action of a serine/threonine protein phosphatase in this compensatory effect. Because calcineurin is a well known Ca2+-dependent protein serine/threonine phosphatase (Klee et al., 1998), we used cyclosporine A to inhibit calcineurin activity in our system. Cyclosporine alone did not affect baseline levels of IHTK (data not shown), but combined application of 4AP with cyclosporine eliminated the compensatory increase in IHTK (Fig. 5A), indicating that the increase in IHTK is dependent, at least in part, on the activity of calcineurin.
In complete contrast, the compensatory influence of TEA on IA appears to be implemented via mechanistically distinct pathways compared with those we found for the 4AP block on IHTK. The increase of IA in response to TEA blockade is neither calcium-dependent (Fig. 5B) nor dependent on the activity of serine/threonine phosphatases (Fig. 5B). Parallel experiments showed that BAPTA, ryanodine, okadaic acid and cyclosporine, when co-applied with TEA, failed to prevent a significant increase in IA (Fig. 5B). The only exception was that the application of ryanodine in conjunction with TEA appeared to reduce the magnitude of increase in IA, although this effect was not statistically significant with respect to TEA alone (p = 0.103; t test), or TEA plus BAPTA (p = 0.08).
We also examined the effects of inhibiting kinase activity on the compensatory effects of 4AP and TEA blockade. Staurosporine is a potent, cell-permeable protein kinase C inhibitor which also partially inhibits other kinases such as PKA, PKG, and CaMKII (Rüegg and Burgess, 1989). Co-application of staurosporine with 4AP appears to cause a significant decrease in IHTK (Fig. 5A). However, this is a more complex result than that for the other pharmacological blockers. Unlike any of the other pharmacological blockers used in this study, staurosporine was the only one to cause a change in baseline levels of IHTK when applied alone (Fig. 5C). Therefore, combined application of 4AP and staurosporine results in a relative decrease in IHTK compared to that with staurosporine alone, but the net effect of 4AP + staurosporine is a restoration of IHTK to baseline levels (Fig. 5D).
Once again, the effects of staurosporine on IA in the context of TEA blockade contrast with those for 4AP and IHTK. There is no effect of staurosporine on baseline levels of IA (Fig. 5C), and while staurosporine co-applied with TEA results in an apparent decrease in the magnitude of the effect on IA, there is still a significant increase in IA with staurosporine present (Fig. 5B), which is not statistically significant from the effect of TEA alone (p = 0.08; t test).
Discussion
Rapid compensation between IA and IKCa preserves both cellular and network outputs
We have identified a naturally occurring coregulatory relationship between potassium currents (IA and IKCa) in an intact CPG network that results in homeostatic compensation of neuronal excitability as well as network function. Furthermore, we have determined that these compensatory changes in K+ current magnitudes are independently regulated by distinct mechanisms. Compensatory increases in IKCa are calcium-dependent and due, at least in part, to the activity of calcineurin-based phosphatase activity. Conversely, compensatory increases in IA are independent of all of the regulatory pathways implicated in the IKCa response. These effects are also fairly rapid, acting over the course of 60–90 min to stabilize the activity of a critical CPG network responsible for cardiac muscle contraction and heart beat generation in the animal. While such homeostatic responses and their role in stabilizing synaptic function have been well studied (Bergquist et al., 2010; Turrigiano, 2011), we know much less about the mechanisms underlying homeostatic plasticity of intrinsic excitability and the role this form of plasticity plays in the stabilization of neuronal and motor network output (Turrigiano, 2011).
Even neurons and networks with extremely robust output display highly variable underlying physiological parameters responsible for neuronal output, particularly membrane conductances (Schulz et al., 2006a, 2007; Khorkova and Golowasch, 2007; Goaillard et al., 2009; Temporal et al., 2012). Our work reveals that embedded within this variability are coregulatory relationships that act to stabilize the excitability of the cell, in part by balancing the sum total of major transient outward currents, IA and IKCa. A similar relationship between IA and IKCa was reported in the STG. Artificial depolarization of inferior cardiac neurons results in increased IA and decreased IHTK, which is abolished by Cd2+ blockade of calcium channels (Golowasch et al., 1999). These experiments did not reveal a naturally existing correlation between these currents (Golowasch et al., 1999), suggesting that excitability in the STG may be a more complex interaction among multiple conductances, or perhaps between conductances and their constitutive neuromodulation (Harris-Warrick, 2011). Although it is not known what effect this has on network output in the STG, our complementary results suggest that such homeostatic mechanisms may be common among motor neurons in different CPGs.
Distinct intracellular pathways mediate compensation
The mechanisms involved in the regulation of the compensatory response mediated by IKCa after 4AP block in the CG network are consistent with work supporting changes in K+ current density in a homeostatic fashion (Desai et al., 1999; Schulz et al., 2006b; Debanne and Poo, 2010; Misonou, 2010). In particular, the mechanisms we see for the compensatory increase in IKCa are similar to those found in regulation of excitability in cultured hippocampal neurons (Misonou et al., 2004, 2005; Misonou and Trimmer, 2004). Increased excitability in these cultured neurons results in calcineurin-dependent dephosphorylation of Kv2.1 channels, leading to a functional potentiation of these channels via a shift in activation voltage, and to restoration of excitability (Misonou et al., 2004). While our data strongly suggest that IKCa is primarily responsible for the observed change in IHTK, we cannot conclusively rule out a role for IKd (similar to Kv2.1). However, cyclosporine has been shown to alter IKCa channels directly (Hay, 1998), and calcineurin is known to have widespread effects on neuronal plasticity and excitability beyond its effects on Kv2.1 activity (Groth et al., 2003). The effects seen in our experiments likely can be attributed to an overall increase in the maximal conductance of IKCa, as we see no changes in either activation curves or any consistent or significant changes in input resistances of these cells (O'Leary et al., 2010).
The mechanisms underlying compensatory increases in IA are less clear, but distinct from the pathways involved in the upregulation of IKCa. Unlike for IKCa, the increase in IA is neither calcium-dependent, nor influenced by the activity of calcineurin or any other phosphatases that would be affected by treatment with okadaic acid (Cohen et al., 1990). Indeed, blocking neither phosphatase nor kinase activity by two broad spectrum blockers was able to prevent an increase in IA as a result of TEA blockade, possibly precluding the contribution of phosphorylation state to this half of the compensation story. However, we were able to dampen the effect by blocking release of intracellular calcium with ryanodine, as well as by inhibiting kinase activity, suggesting a more complex path of regulation in this response. Given the potential nonspecific effects of pharmacological treatment, we cannot rule out at least some role for a kinase in this pathway. For example, Kv1.2 potassium channels are known to be affected by cAMP/protein kinase A pathways which enhance their conductance, in part by altering trafficking of these channels (Connors et al., 2008). Our staurosporine treatment may result in only partial inhibition of PKA (Rüegg and Burgess, 1989), which could explain the intermediate effects seen on the A-type current compensation. Regardless, the fact that IA and IKCa are regulated by distinct pathways demonstrates an inherent robustness to this homeostatic response, which may be characteristic of networks responsible for critical functions, such as CPGs.
Reservoirs of channel proteins might implement rapid compensation
Our data also begin to shed light on a striking disparity between levels of regulation in this system. Namely, we see an entirely different relationship at the level of mRNA for the channel genes BKKCA and shaker (positive correlation), than we do for the currents IKCa and IA that these channels most directly encode (negative correlation). Indeed, to date every correlation we have found between channel mRNAs at the single cell level has been a positive correlation (Schulz et al., 2007; Tobin et al., 2009; Temporal et al., 2012). These data include measurements of at least nine different channel genes in seven different motor neuron types in the crustacean stomatogastric and cardiac ganglia (S. Temporal and D. J. Schulz, unpublished data), for a total of no fewer than 60 correlations in channel expression detected among channel transcript levels, and every one has thus far been positive. In addition, our data impact the present, although admittedly limited, understanding of the relationship between channel mRNA and ionic current in a given neuron. Namely, that mRNA for a given K+ channel correlates positively with its corresponding ionic conductance (Schulz et al., 2006). Although both of these measurements were not made in the same cells in this study, logic dictates that there cannot simply be a positive relationship between mRNA and conductance for both K+ currents in this study. Clearly, and not particularly surprisingly, our data demonstrate that more complexity lies within the relationship of mRNA and ionic conductance than a one-for-one tracking between these disparate levels of function in all cell types. However, our data continue to support the concept that correlation among mRNA levels (Schulz et al., 2007; Tobin et al., 2009) underlies a functional relationship at the conductance level.
So then why do we see only positively correlated mRNA levels for channel genes in these motor networks? The linkage to cellular and network output in our study enables us to propose a new hypothesis in this regard. We hypothesize that rapid compensation of intrinsic excitability, such as that reported in this study, is implemented via an already existing pool of protein. That is, protein pools for the two opposing channels must already be available in relatively equivalent numbers, even if they are not functionally equally represented. This requires that mechanisms must occur at the level of gene regulation to ensure a sufficient “reservoir” of protein capable of compensating for the loss of a given channel. Thus, regardless of whether channels act in concert to set neuronal output (Ball et al., 2010; Franklin et al., 2010), or act in a compensatory fashion, the overall gene regulation between the two must be relatively balanced. Such a reservoir perspective is reminiscent of mechanisms seen for synaptic vesicle protein dynamics (Fernández-Alfonso et al., 2006), or extra-synaptic AMPA receptors that serve as a reservoir during synaptic plasticity (Hayashi et al., 2000; Zhu et al., 2000). Interestingly, insertion of AMPA receptor during synaptic plasticity is also under the control of phosphorylation-dependent processes (Lin et al., 2009). Therefore, this hypothesis provides a compelling framework for future investigation of the relationship between mRNA levels and ionic conductance, as well as the interplay between gene expression and post-transcriptional mechanisms involved in homeostatic plasticity of intrinsic excitability over different time scales.
Putting it all together
Our findings represent one of the first comprehensive demonstrations of rapid homeostatic plasticity of intrinsic excitability that results in a stabilization of output in a mature, intact network of a CPG. These rapid compensatory increases are mechanistically independent, suggesting robustness in the maintenance of neural network output that is critical for survival. Furthermore, this study reveals a distinct mechanism for compensation that leads us to at least one working model of homeostatic plasticity in this system (Fig. 6). We hypothesize that one pathway to functional compensation in this system relies on intracellular calcium concentration as a measure of cell excitability (Kennedy, 1989; Ross, 1989; LeMasson et al., 1993). A block of A-type K+ channels leads to an increase in the excitability of the cell, causing release of intracellular calcium stores. This calcium influx alters calcineurin activity, resulting in dephosphorylation of targets, perhaps the KCa-channels themselves, that increases IKCa. This ultimately restores the outward current balance that was lost during decreased IA and thus restores excitability. The corresponding mechanism responsible for compensatory increases in IA (Fig. 6) is presently unknown, but is likely to involve distinct mechanisms for monitoring excitability (Dirnagl et al., 2003), as well as other candidate mechanisms that alter A-type conductances (Connors et al., 2008). Finally, we hypothesize that positive coregulation of mRNA numbers for channel genes may ultimately provide a reservoir of channel protein (Fig. 6) for implementing compensation over rapid time scales during which de novo synthesis may be insufficient to implement full compensation.
Footnotes
This work was supported by a Craig H. Neilsen Foundation Grant 83026 (to D.J.S.), the Missouri Spinal Cord Injuries Program (to D.J.S.), and Department of Defense–Congressionally Directed Medical Research Programs Exploration-Hypothesis Development Award SC090555 (to D.J.S.).
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. David J. Schulz, University of Missouri, Department of Biological Sciences, Columbia, MO 65211. schulzd{at}missouri.edu