Temperature compensation in a small rhythmic circuit
Abstract
Temperature affects the conductances and kinetics of the ionic channels that underlie neuronal activity. Each membrane conductance has a different characteristic temperature sensitivity, which raises the question of how neurons and neuronal circuits can operate robustly over wide temperature ranges. To address this, we employed computational models of the pyloric network of crabs and lobsters. We produced multiple different models that exhibit a triphasic pyloric rhythm over a range of temperatures and explored the dynamics of their currents and how they change with temperature. Temperature can produce smooth changes in the relative contributions of the currents to neural activity so that neurons and networks undergo graceful transitions in the mechanisms that give rise to their activity patterns. Moreover, responses of the models to deletions of a current can be different at high and low temperatures, indicating that even a well-defined genetic or pharmacological manipulation may produce qualitatively distinct effects depending on the temperature.
Introduction
Biological systems depend on many interacting nonlinear processes that together produce complex outputs. In the nervous system, neuronal activity requires the coordinated activation and inactivation of many inward and outward currents. Temperature influences all biological processes, to a greater or lesser degree. This poses an inherent difficulty for neuronal signaling: if the currents involved in neuronal and network dynamics are differentially temperature-dependent a system that is well-tuned to work at one temperature may not function at a different temperature (Caplan et al., 2014; O'Leary and Marder, 2016; Tang et al., 2010; Tang et al., 2012). Nonetheless, many ectothermic animals have neurons and circuits that function well over an extended temperature range (Robertson and Money, 2012). It then becomes important to understand how and to what extent this can occur.
The effects of temperature on channel function and neuronal activity have been studied extensively for many years (Taylor and Kerkut, 1958). Increasing temperature generally results in increases of channel maximal conductances and faster activation/inactivation rates (Frankenhaeuser and Moore, 1963). But, importantly ionic channels of different types are affected by temperature to different extents: each of these processes has a different (Schauf, 1973; Kukita, 1982; Ruff, 1999; Tang et al., 2010; Ranjan et al., 2019). The effect of temperature on neuronal intrinsic excitability such as voltage and current thresholds can be diverse (Sjodin and Mullins, 1958; Guttman, 1962; Guttman, 1966; Fitzhugh, 1966). For example, identified neurons in locust showed reversible changes in spike amplitude and duration of spikes as temperature varies from 18 to 35°C (Heitler et al., 1977). Studies in the jump neural circuit of grasshopers showed that neural excitability is differently affected by temperature across neuronal types (Abrams and Pearson, 1982).
There are also examples of neuronal and circuit processes that are relatively temperature-compensated. The frequency/current (f/I) curves in both molluscan and locust neurons are not substantially affected by small temperature changes (6–8°C) (Connor, 1975; Heitler et al., 1977). The f/I curves of auditory sensory neurons in grasshoppers remain largely unaffected by changes in temperature between 21°C and 29°C (Roemschied et al., 2014). Temperature compensation also takes place across the behavioral level. For example, a recent study in hunting archerfish showed that the duration of the two major phases of the C-start—a fast escape reflex which follows prey release—were temperature compensated (Krupczynski and Schuster, 2013). The question then arises of how these behaviors are preserved when the currents in the cells and their intrinsic excitability properties are differentially modified by temperature.
Here, we explore these issues using computational models of the pyloric network—a subnetwork within the stomatogastric ganglion (STG) of crustaceans (Marder and Bucher, 2007; Maynard, 1972). The pyloric rhythm is a triphasic motor pattern that consists of bursts of action potentials in a specific sequence. This behavior is robust and stable, and the cells and their connections are well-characterized. Robustness to temperature in the pyloric network has been explored experimentally by Tang et al., 2010; Tang et al., 2012; Soofi et al., 2014 and Haddad and Marder, 2018. These studies show that as temperature increases the frequency of the pyloric rhythm increases but the phases of the cycle at which each cell is active remain approximately constant. Additional work on the pacemaker kernel of the pyloric network (three cells connected by gap junctions that burst synchronously) showed that temperature increases the frequency of these bursts (), but their duty cycle (the burst duration in units of the period) stays approximately constant (Rinberg et al., 2013).
Temperature robustness was explored in computational models of the pacemaking kernel (Soto-Treviño et al., 2005) by Caplan et al., 2014 and O'Leary and Marder, 2016. They showed that it was possible to find multiple sets of values for different membrane conductances processes, so that the duty cycle of the cells remained constant as their bursting frequency increased. In this work, we build on the results in Caplan et al., 2014 and implemented temperature sensitivity in a model of the pyloric network (Prinz et al., 2004). We show that in these models, there are multiple sets of maximal conductances and temperature sensitivities that reproduce much of the experimental phenomenology previously reported (Tang et al., 2010; Tang et al., 2012; Soofi et al., 2014; Haddad and Marder, 2018). In addition, we explored how the dynamics of the currents are modified to sustain the correct activity at each temperature. We performed this study for 36 different models and found that in all cases, the contributions of the currents to the activity can be significantly different across temperatures. The currents are not simply scaled up but instead become reorganized: a current that is important for burst termination at 10°C may no longer play that role at 25°C. Because the contribution of a given current to a neuronal process can be replaced by another at different temperatures, deletion or blockade of a current can produce qualitatively different effects at high and low temperatures. These results provide a plausible hypothesis for why interactions between temperature and a second perturbation can be observed experimentally (Haddad and Marder, 2018; Ratliff et al., 2018).
Results
Duty cycle and phase maintenance in model pyloric networks
The triphasic pyloric rhythm is produced by the periodic sequential activation of the pyloric dilator (PD) neurons, which are electrically coupled to the anterior burster (AB) neuron forming a pacemaking kernel, the lateral pyloric (LP) neuron and five to eight pyloric (PY) neurons. We modified the model of the pyloric circuit in Prinz et al., 2004 to include temperature sensitivity. Figure 1A shows a schematic representation of the model pyloric network studied here. The model network consists of three cells, each modeled by a single compartment with eight currents as in previous studies (Golowasch and Marder, 1992; Buchholtz et al., 1992; Goldman et al., 2001). The synaptic connections are given by seven chemical synapses of two types as in Prinz et al., 2004. In this model the pacemaking kernel is aggregated into a single compartment (here we refer to this compartment as ). Following Liu et al., 1998, each neuron has a sodium current, ; transient and slow calcium currents, and ; a transient potassium current, ; a calcium-dependent potassium current, ; a delayed rectifier potassium current, ; a hyperpolarization-activated inward current, ; and a leak current . The traces in Figure 1A show a solution of this model for one set of maximal conductances . The traces exhibit a triphasic pyloric rhythm that consists of the sequential bursting of the , and cells. This pattern remains approximately periodic so we can measure the phases of the cycle at which each burst begins and terminates using the burst start of the cell as reference (indicated as in magenta in Figure 1A), as indicated by the color labels.
Temperature sensitivity was modeled for each conductance as in Caplan et al., 2014 and O'Leary and Marder, 2016 (see Materials and methods). We increased each conductance by an Arrhenius type factor (with and = 10°C) and also increased the rates of the channel kinetics in similar fashion. We assume that for the maximal conductances and for the timescales as these ranges are consistent with experimental measurements of these quantities (Tang et al., 2010). The values are a property of the channel proteins, so in a given model we use the same for a given current (and synapse) type in all cells. Different values of the maximal conductances can produce an innumerable variety of activities. The space of solutions of the model—the different patterns produced for each set of maximal conductances —is formidably complex. This is in part due to the nonlinear nature of neuronal dynamics but also because the number of conductances that need to be specified is large (31 in total, 8 intrinsic conductances ×3 cells +7 synaptic conductances). The picture is more complicated as we consider the effect of temperature because it affects both the conductances and the time scales of the different processes (14 in total: 12 intrinsic, 2 synaptic). Therefore in the models, changes in temperature correspond to a coordinated change or a path in a 45-dimensional parameter space. The complexity of this problem led us to confine this study to specific questions suggested by biological observations over a permissible temperature range.
In the crab, as temperature is increased pyloric activity remains triphasic—the cells fire in the same order and at the same relative phases—while the pyloric network frequency increases by a two- to three-fold factor over a 15°C range (Tang et al., 2010). For this to happen in the model, there should be many regions in the parameter space for which the activity is triphasic at different frequencies and paths in the parameter space (sequences of values of and ) that interpolate these regions. One of the main results in this work is that such paths do exist and that temperature compensation is possible in this model. Finding these temperature compensated networks is nontrivial: it requires the specification of 31 maximal conductances and 24 values (see Materials and methods). Because it is virtually impossible to find these sets of parameters at random we employed a landscape optimization scheme described previously (Alonso and Marder, 2019) and adapted it to this particular scenario (see Materials and methods). We first searched for sets of maximal conductances that displayed triphasic rhythms at control temperature (10°C). Inspired by physiological recordings of the pyloric rhythm of crabs and lobsters we targeted control activities with a network frequency of , duty cycle ≈ 20% for , and duty cycle ≈ 20% for and (Bucher et al., 2006; Hamood et al., 2015). We then selected 36 of these models and for each of them we searched for values so that the activity is preserved over a temperature range. We found that, regardless of the maximal conductances , it was always possible to find multiple sets of values that produced a temperature robust pyloric rhythm.
Figure 1B shows the activity of one example network at several temperatures, similar to the experimental results reported by Tang et al., 2010. Figure 1C shows the same analysis as in Tang et al., 2010 performed on this particular model. We simulated the model at several temperatures over the working range (10°C to 25°C) and computed the average duty cycle, frequency and phases of the bursts. While the burst frequency increases with , the duty cycle of the cells, and the phases of the onsets/offsets of activity, remain approximately constant. Figure 1—figure supplement 1A shows the same analysis performed on three different models. In all cases, the duty cycle and the phases stay approximately constant while the frequency increases. Figure 1—figure supplement 1B shows the values of the maximal conductances and the values for each of these models.
While the duty cycle remains approximately constant on average, the spiking patterns of the cells show discernible differences across temperatures. We subjected the models to temperature ramps from 10°C to 25°C over 60 min and recorded the spike times (Figure 1—figure supplement 2). As temperature is increased the spiking patterns change in complicated ways and the largest decreases monotonically, consistent with the overall increase in bursting frequency of the cells. Different values of temperature result in slightly different spiking patterns which in turn result in different values of the duty cycle. For some temperatures, the duty cycle is nearly identical in every burst resulting in a single point in the y-axis (Figure 1—figure supplement 2D, blue box in PD cell). There are temperatures for which the duty cycle takes two values preferentially (Figure 1—figure supplement 2D, pink box in cell) and there are temperatures for which the duty cycle differs noticeably from burst to burst. Although in all models, the average duty cycle and phases stay approximately constant, as temperature is changed the precise way in which the spiking patterns change show marked differences across models: different maximal conductances and different values produce different patterns. Figure 1—figure supplement 3 shows the temperature dependence of the duty cycle in the cell over the working range for all 36 models and provides a notion of how variable the bursting patterns of the cells can be across temperatures and individuals.
Changes in membrane potential over temperature
Experimental studies on other systems reported measurable effects of temperature on resting membrane potentials (Klee et al., 1974), spiking thresholds, and amplitudes of spikes (Heitler et al., 1977). In the models studied here, some temporal properties of the activity remain approximately constant over a temperature range but the precise shape of the waveforms can also change. To inspect how the membrane potential changes over the working temperature range (10°C - 25°C), we computed the distribution of of each cell for 101 values of temperature. Temperature changes features of the voltage waveforms, resulting in changes in these distributions. Figure 2 shows example traces and the result of this analysis for the cell in one model. The membrane potential changes with temperature in several ways: the spike amplitude decreases and the spike threshold depolarizes as temperature increases. At each temperature, we computed the distribution of , and used a gray scale to indicate the fraction of time that the cell spends at each voltage value (y-axis). The color lines match features of the voltage waveforms at each temperature, with features in the distribution on the right. The distributions permit visualizing changes in the waveform as the control parameter is changed (Alonso and Marder, 2019). The effect of temperature on membrane potential is not consistent across models and depends on the precise values of conductances and temperature sensitivities. Figure 2—figure supplement 1 shows the membrane potential distributions of each cell for three example models. In all cases, the waveforms show visible differences across temperatures and there is considerable variability in the amplitudes of the oscillations across models. In some models, the peak voltage of spikes, as indicated by the upper envelope of the distributions, decreases (left), increases (middle), or remains relatively constant (right). In addition, we found that cells can either become more depolarized or more hyperpolarized as temperature increases.
Dynamics of the currents at different temperatures
Neuronal activity is governed by currents that result from the precise activation and inactivation of ion channels with different kinetics. Experimental access to these quantities is limited because it is hard to measure currents individually without blocking all other currents, and thus changing the activity. Here, we employ models to explore how these currents are differentially altered by temperature, and yet are able to remain balanced in such a way that the network activity speeds up while preserving its behavior and phase relationships.
We first compared the dynamics of the currents in two model neurons with different maximal conductances and temperature sensitivities by direct inspection of the currents’ time series. Figure 3A shows the time series of each intrinsic current in at 10°C and 25°C for one model. The corresponding membrane potential activity is shown in the blue traces on top. The red horizontal dashed lines indicate the peak amplitude of the currents at control temperature (10°C) to visualize whether it increases or decreases at 25°C. This is indicated by the arrows in each row (in colors: ⇑ red increases, ⇓ green decreases). In this example, the current increases its amplitude by more than two fold over the 15°C range. Despite the fact that increases with temperature, the amplitude of the current decreases to about half its control value at 10°C. Similar observations can be made for other currents (Figure 3A). Note that all values in these models are greater than 1. Therefore, when the current decreases it is because the relative timing of the currents’ activation is changed by the neurons’ firing properties. The way currents are modified by temperature is different across models. Figure 3B shows the currents in the PD cell at two temperatures for another model. In this case, the current decreases by a small amount despite the fact that the maximal conductance and timescales of this current increase with temperature. The current which in Figure 3A grows by a factor of six, is moderately decreased by temperature in Figure 3B.
Because the total number of currents in the circuits is 31, it becomes cumbersome to compare them separately across temperatures (and models). For this reason, we computed and inspected their currentscapes (Alonso and Marder, 2019). The currentscapes use colors to show the percent contribution of each current to the total inward (or outward) current and are useful to display the dynamics of the currents in a compact fashion. Figure 4 shows the currentscapes of each cell in one model at 10°C and 25°C. The shares of each current of the total inward and outward currents are displayed in colors over time. The total inward and outward currents are represented by the filled black curves on a logarithmic scale at the top and bottom. The currentscapes are noticeably different at 10°C and 25°C indicating that the currents contribute differentially across temperatures. For example, in the cell the current contributes a visible share of the inward current both at the beginning and the end of the burst at 10°C, but at 25°C its contribution is mainly confined to the beginning of the burst. The current is evenly distributed during the burst at 10°C but at 25°C its contribution is visibly larger at the beginning of the burst. In all cells, but most notably in the and neurons in this example model, temperature substantially increases the contribution of the current toward the end of the bursts.
The currentscapes show that in all models and at all temperatures, there are periods of time when the activity is dominated by one or two currents, and periods of time when several currents contribute to it by similar amounts. During spike production, the activity is first driven by and later by . During these periods, the total inward and outward currents differ by orders of magnitude and are entirely dominated by one current type. Similarly, when the cells are most hyperpolarized the dominant currents are and the synaptic currents . There are periods of time—such as the beginnings and ends of bursts—when the total inward and outward currents are composed of comparable contributions of several current types. It is during these periods, when multiple currents act together to change the membrane potential, where we observe the most dramatic changes with temperature.
Figure 5 shows 20 ms after the burst in a neuron over temperatures between 10°C and 25°C. Almost all of the currents change over this temperature range, most notably there is almost a tradeoff between the contributions of the and currents, with decreasing smoothly and increasing smoothly. The specifics of how the contributions change is different across models, but in all cases there is a smooth transition between different relative fractions of currents. This allows the neuron to slide smoothly through changes in burst termination mechanisms.
Different compensation mechanisms fail in different ways. Figure 6 shows full currentscapes for two network models as they ‘crash’ at high temperature. Note that the complete loss of activity in the neurons results from different currents being active in the two examples. Further examples of high temperature ‘crashes’ are shown in Figure 6—figure supplement 1. These examples show that the networks are crashing by different mechanisms in each case.
High temperatures produce disordered circuits or ‘crashes’
For temperatures higher than a critical temperature, the biological network displays anomalous regimes characterized by the emergence of slow time scales and intermittency between what appear to be metastable states (Tang et al., 2012; Haddad and Marder, 2018). At temperatures higher than 25°C some models do not display irregular behavior but just become quiescent, or produce subthreshold oscillations. Other models display irregular and seemingly aperiodic regimes. As expected, irregular regimes typically take place at values of the temperature near a transition between qualitatively different patterns of activity. Figure 7 shows the activity of five models at temperatures higher than their working regimes. The left panels show 30 s of data and the right panels show 5 s. In all cases, there are timescales in the behavior that are much longer than the period of the rhythm at 10°C (≈ 1 s). The dynamics of the models in these regimes is daunting and their characterization is beyond the scope of this work. However, qualitative features of these states such as the appearance of slow timescales or absence/presence of activity in one or more cells, are captured by these models. This, together with the multistable regimes shown before plus a source of noise, may provide a reasonable model to account for the irregular regimes in the biological network (Tang et al., 2012; Haddad and Marder, 2018).
Hysteresis and multistability
It is natural to ask whether these models show hysteresis when the temperature is returned to control levels. Moreover, multistability is also often observed in conductance-based models of neuronal activity and biological neurons (Cymbalyuk and Shilnikov, 2005; Malashchenko et al., 2011; Marin et al., 2013; Lechner et al., 1996; Paydarfar et al., 2006), and the network studied here also exhibits such features. We explored hysteresis in the models by increasing temperature linearly from 10°C to 35°C in 30 min and then decreasing it back to 10°C at the same rate. Figure 8 shows the spiking patterns for two neurons during such temperature ramps (schematized on top). In Figure 8A the ISI distributions show bursting activity up to a critical temperature at which the model becomes quiescent. The quiescent state remains stable as temperature continues to increase and eventually recovers activity during the decreasing ramp, but at a lower temperature than when it stopped spiking during the increasing ramp. This implies that there is a range of temperatures where the model can be either quiescent or spiking, and thus is multistable. In most models (34 of 36), the spiking patterns over the working temperature range are similar across the up and down ramps but are not identical, suggesting there are multiple stable attractors that correspond to the triphasic rhythm but differ slightly in their spiking patterns (Figure 1—figure supplement 2E). This is consistent with the findings in Tang et al., 2010 where temperature was ramped slowly and no significant differences were found when ramping up and down in temperature.
We also found models with pathological stable attractors that coexist with the pyloric rhythm at working temperatures. Figure 8B shows a model in which hysteresis is salient. At high temperatures (times between 25 and 35 min), this model produces values as long as ≈ 1 s and remains spiking at all temperatures. The spiking patterns differ noticeably on the up and down ramps. The spiking pattern of the cell changes from regular bursting to a pathological state that remains stable until the model is ramped back down to ≈ 20°C. In our simulations, hysteresis mainly occurs at high temperatures indicating that states such as those in Figure 7 can be multistable. The extent such hysteresis can be revealed in the biological experiments depends on the details of the perturbation protocol, and its currently under study.
Response to perturbations at different temperatures
In all models, the dynamics of the currents are different at different temperatures and their contributions become reorganized in complex ways. The combined activity of all the variables in the model—the dynamical attractor—is therefore expected to have different stability properties at different temperatures. This means that an extreme perturbation can have qualitatively different effects at different temperatures. The responses of different models to extreme perturbations such as partially or completely removing a current can be diverse at any temperature (Alonso and Marder, 2019). While this is expected due to the complexity of these models, the fact that these responses can also be qualitatively different across temperatures indicates that changes in the current contributions are meaningful and affect the stability of the network. To shed light on this issue, we performed several simple perturbations.
Figure 9 shows one perturbation in two models. In Figure 9A, we explored the interaction between temperature and the effect of removing the -type K+ channel at two temperatures: the models were simulated for 20 s in control conditions and then was completely removed for the next 20 s. Figure 9A (top) shows that removing the current at 10°C has catastrophic consequences for the network activity which becomes irregular. When the same perturbation is performed at 25°C (bottom), the triphasic rhythm is almost normal except for the irregular bursting of . Figure 9B shows a similar pertubation in which we remove the current. This model displayed normal triphasic behavior at 10°C (top) and when the current was removed its activity slowed but the triphasic rhythm was maintained. The same perturbation at 25°C results in quiescence. In this model, the current is necessary for the activity at 25°C but not at 10°C.
To further characterize the differential response to perturbations at different temperatures, we performed a second assay in which we gradually decreased a current from control (100%) to complete deletion (0%) in five steps. We performed these simulations for all 36 models and all current types at 10°C and 25°C and compared their responses. In general, gradually removing a current can result in qualitatively different states depending on the temperature at which this perturbation is performed. The responses are so diverse that a coarse classification scheme is sufficient to highlight this observation. We classified the network activity based on the spiking patterns and waveforms of each cell, and also by their relative activities. The activities of each cell were classified as: regular bursting, irregular bursting, tonic spiking, irregular spiking, single spike bursting, quiescent and other, while the network activity was classified as triphasic or not triphasic. The classification scheme consists of a simple decision tree based on the statistics of spiking and is sufficient to tease apart the different regimes we observe at a coarse level (Materials and methods).
Figure 10A shows the result of classifying the responses of one model to gradually decreasing a conductance in five steps, at two different temperatures, for four conductances. We employed 5 × 4 response grids to summarize the effects of partially or completely removing a conductance at 10°C and 25°C. The rows in the grids correspond to different values of conductance removals, with the control condition (100%) on top, and the completely removed condition (0%) at the bottom. The first three columns correspond to each of the cells (PD-LP-PY) and the fourth column corresponds to the network (NET). The colors indicate the type of activity of each cell and the network color coded for condition as indicated by the labels in the figure. The top left panel in Figure 10A shows the responses of one model to removing the sodium current . At 10°C, when the cell bursts irregularly and the activity of the network is not triphasic, but if the same perturbation is performed at 25°C the activity remains triphasic. Removing the current in this model has little effect at 10°C as the activity remains triphasic. However, this current plays an important role at 25°C and the model becomes quiescent upon complete removal of . The same observation that the responses to perturbation are different at each temperature holds regardless of the perturbation type. In the cases of the and currents complete removal results in different non-triphasic activity at each temperature. The responses of the model are in general qualitatively different at different temperatures, and they are also diverse across models (not shown). The fact that the models ‘break-down’ in different ways is consistent with the observation that these network mechanisms are diverse. In all models, the responses to extreme perturbation are temperature dependent for most perturbation types. In Figure 10B we quantify how many models respond differently to complete removal of a current at 10°C and 25°C. This corresponds to comparing the bottom row of the response grids at each temperature, for each conductance type. In most cases, about half of the models respond to complete removal in different ways at each temperature. The exceptions are the current for which most models respond in different ways, and the current where complete removal results in equivalent states (quiescence in this case) at both temperatures.
Discussion
Enzymes and ion channels show altered behavior when temperature is changed, as temperature influences their conformation. Yet many animals live successfully over wide temperature ranges. For example, North Atlantic lobsters and crabs routinely experience more than 25°C temperature swings during a year, and may see 8 - 10°C changes within a short time period (Haefner, 1977; Stehlik et al., 1991). This raises a series of fascinating questions of how behavior and circuit performance can be maintained while all biological molecules are, to a greater or lesser degree, differentially temperature dependent (Robertson and Money, 2012).
It is important to distinguish between altered activity that occurs rapidly in response to acute, or immediate changes in temperature and those changes that occur during slow adaptation and acclimation, processes that can take days, weeks or months. In this study, we have modeled the fast changes that we assume would occur ‘automatically’ as proteins change conformation and properties rapidly, and as are usually measured in their ’s. This is to be contrasted with the slower, long-term effects that result in acclimation in which a variety of molecular processes are triggered that may in turn change the proteins present in the membrane. These mechanisms, and much adaptation to long-term changes in temperature depend on the activation of heat shock proteins, RNA editing and a whole array of molecular mechanisms (Garrett and Rosenthal, 2012; Rosenthal, 2015; Martin Anduaga et al., 2019). We will first discuss the issues that arise from the examination of robustness to acute temperature changes, as modeled here, and then offer some speculations for additional processes that may be important to enhance robustness in response to long-term temperature changes.
Unfortunately, there are relatively few data relevant to whether the of a channel in a given animal changes in response to acclimation or seasonal conditions. Therefore, although it is often assumed that a for a given protein is invariant over the lifetime of the animal this may not necessarily be the case. Another difficulty is that there are few, if any, studies of many ion channel s in the same animal in the same temperature range. Computational modeling will not answer what the ‘real’ conductance and combinations are.
The assumption that these quantities are stationary throughout the life of an animal (and/or the same across species) may not be warranted. For example, long-term acclimation over months resulting in RNA editing, could alter some s. Moreover, it is possible that minor genetic variations in channel proteins exist in the population so that different animals in the population could show modestly different s. Therefore, the point of this study is to ask generically how compensation occurs acutely assuming different temperature sensitivities for the processes, within the experimental range. For these reasons, in the absence of data, it is valuable to explore the space to have on overview of how compensation takes place acutely at any point during the lifetime of an animal.
One of the most important and non-intuitive results from this work is seen in Figure 5, which demonstrates that a neuron can smoothly slide through mechanisms that govern a physiological process as temperature is changed. In the case of Figure 5, the repolarization of the last spike in the burst shows a large contribution from at 10°C but becomes small at 25°C while contributes little at 10°C and becomes dominant at 25°C. This shows that while both of these currents can contribute to spike and burst repolarization, their shares in these processes change, and the neuron is robust as it smoothly slides through these mechanisms. This smooth transition can occur precisely because these three K+ currents depend on voltage, time, and temperature differently.
We speculate that one advantage conferred by the expression of many different ion channels that differ in their properties is to provide many sets of such ‘sliding stability mechanisms’. Individual neurons in all species have large numbers of ion channel genes that encode proteins that give rise to the membrane currents that are usually measured and characterized in voltage-clamp. The model neurons used are loosely based on voltage-clamp measurements from crab and lobster stomatogastric ganglion (STG) neurons (Golowasch et al., 1992; Turrigiano et al., 1995) and have eight currents. That said, we now know that there are at least twenty to thirty and probably more voltage-gated ion channel genes in individual STG neurons (Northcutt et al., 2016). In some cases, several of these genes may contribute to one of the currents that we model here. For example, the delayed rectifier current modeled here may have contributions from at least three different genes, and , modeled here, also probably has contributions from at least two different genes. Given that these genes encode similar but not identical protein subunits, it is not unreasonable to expect that these proteins may respond differently to temperature, and that the actual composition of the activated currents may reflect a change in these contributions as a function of temperature. This again, could provide still another ‘sliding stability mechanism’, and it would be interesting in the future to create models using data from currents that have been measured in isolation, such as in Ranjan et al., 2019.
The same principle of sliding circuit mechanisms is likely to hold when we approach network dynamics. The phase relationships of the pyloric rhythm are maintained over a considerable temperature range (Tang et al., 2010; Soofi et al., 2014). Likewise, the pyloric phase relationships are also maintained over a considerable frequency range (Hamood et al., 2015). In this latter case, it has been argued that the phase compensation depends on the conjoint action of a number of different cellular mechanisms, including synaptic depression (Manor et al., 1997; Hooper, 1997), and the activation and inactivation of and (Nadim et al., 1999; Nadim and Manor, 2000; Tang et al., 2010; Harris-Warrick et al., 1995). So, here the principle also holds: resilience and robust function may require smoothly moving between a variety of different cellular mechanisms.
Long-term mechanisms for temperature adaptation and acclimation may extend robustness by changing the functional properties of ion channel proteins. In biological systems, sustained exposure to temperature often activates heat shock proteins and has effects on many cellular and molecular functions (Sharp et al., 1999). Interestingly, changes in ion channel composition can result from RNA editing, or changes in splice variants (Johnson et al., 2011; Lin and Baines, 2015). It is easy to imagine that these kinds of molecular mechanisms can result in shifts in the temperatures at which neurons and circuits may function. Acclimation to high temperatures usually requires at least 3–4 weeks of sustained alteration of the environmental temperature, and this is enough time to replace a large subset of ion channel proteins in the membrane. This raises a series of fascinating questions for how homeostatic regulation of excitability may be preserved if the molecular mechanisms that are engaged in homeostasis are also influenced by temperature.
Hysteresis in response to perturbation can result from a variety of time-dependent mechanisms. Figure 8 shows examples of hysteresis that can occur in response to temperature ramps, even in models in which the effects of temperature on the properties of each ion channel are viewed as essentially instantaneous. This hysteresis can occur because at high temperatures (> 25°C) the network can produce several different patterns of activity or dynamical states (multistability). The activation and inactivation rates that characterize ion channel function are different for each of these states and therefore, as temperature is lowered, they lose stability in different ways. Of course, any long-term molecular mechanisms (not modeled here) that alter channel expression or splicing or phosphorylation, could produce long-lasting hysteresis.
In this work, we used a genetic algorithm (Alonso and Marder, 2019) to find temperature robust networks. It is important to reiterate that temperature robust neurons and networks are difficult to find doing random searches (Caplan et al., 2014). Therefore, despite the significant animal-to-animal differences in conductance densities seen across the population, randomly sampling conductance values is unlikely to result in successful solutions for individual animals. This argues that evolution has arrived at a series of biological mechanisms that give rise to successful temperature compensation. These are likely to include rules by which correlated values of conductances are produced, such as those found in homeostatic models (O'Leary and Marder, 2016; O'Leary et al., 2013).
Taken at face value, models suggest that conductance densities drift throughout the life of an animal as a result of homeostatic processes (LeMasson et al., 1993; Liu et al., 1998; Golowasch et al., 1999a; O'Leary et al., 2014). Conductance densities can also change in an activity-dependent manner as a result of perturbations (Turrigiano et al., 1994; Golowasch et al., 1999a; Golowasch et al., 1999b; Golowasch et al., 1999a; Santin and Schulz, 2019; Golowasch, 2019). The temperature sensitivities of ion channels are dictated by their molecular structures, and for this reason it is tempting to assume that the temperature sensitivity of a given channel is similar throughout the lifetime of the animal. Nonetheless, as all channels are subject to splice variants and other mechanisms of molecular regulation, it is possible that environmental temperature results in changes in that are consistent with acute robustness, even if there have been changes in conductance densities. This suggests that conductance densities and regulation may be coordinately controlled in biological systems.
One of the most significant results of this paper is that temperature robustness of network function can take place in multiple different ways. In these models, both the membrane potential and the spiking patterns are affected by temperature changes and these changes appear to be different in every model we inspected. Across the models, different values of the maximal conductances G and temperature sensitivities result in consistent differences in the duty cycle distributions. Measurements of the values of the maximal conductances in the STG show large variability across individual animals (Goaillard et al., 2009; Schulz et al., 2006; Schulz et al., 2007; Temporal et al., 2014; Northcutt et al., 2016) so our expectation is that the duty cycle distributions of the biological cells will also display intricate dependencies with temperature, and that these distributions will be different across individuals.
Temperature is not the only perturbation that crabs and lobsters experience. As with temperature, the responses of the STG to changes in pH are diverse across individuals (Haley et al., 2018), again consistent with the large amount of animal-to-animal variability in the expression of ion channels (Schulz et al., 2006; Temporal et al., 2014; Tran et al., 2019). It is therefore reasonable to assume that the responses to a global perturbation are diverse across individuals because different compositions of channel densities—which produce similar pyloric rhythms—are differentially resilient to any given perturbation. By the same token, as currents change as a function of temperature we expect that a second global perturbation of an individual may have qualitatively different effects at different temperatures, as is illustrated with the models in Figure 10. The interaction between temperature and a second perturbation is the subject of recent experimental studies (Haddad and Marder, 2018; Ratliff et al., 2018). These studies are consistent with the interpretation that different current configurations have different stability properties and that temperature changes these configurations.
One of the general take-home lessons from this work is that models can be resistant to a pharmacological perturbation at one temperature but the same model can be sensitive to the same perturbation at another temperature (Figure 10), and that different individuals may differ in their responses to the same perturbation. The potential implications of this for therapeutics are evident. But also evident, is that the search for therapeutic agents can be fraught if the assay itself is differentially sensitive to temperature or another perturbation. Moreover, the ‘sliding mechanisms’ provided by multiple ion channel genes and their molecular regulation suggests that over time therapeutic sensitivity to a given pharmacological intervention could be quite variable. These findings may contribute to our understanding of individual variable therapeutic efficacies within the population and across time within an individual.
Materials and methods
The model
Request a detailed protocolThe activity of the cells was modeled using single-compartment models similar to those described previously (Turrigiano et al., 1995; Liu et al., 1998; Goldman et al., 2001; Alonso and Marder, 2019). Each neuron has a sodium current, ; transient and slow calcium currents, and ; a transient potassium current, ; a calcium-dependent potassium current, ; a delayed rectifier potassium current, ; a hyperpolarization-activated inward current, ; and a leak current . The number of state variables per cell is 13. The units for voltage are , the conductances are expressed in and currents in . The pyloric network consists of three cells with the same ion channel types but different conductance densities. The interactions in the network consist of seven chemical synapses and are similar to Prinz et al., 2004. The synaptic current is given by , where is the synapse strength, is the membrane potential of the postsynaptic neuron and is the reversal potential of the synapse. The activation of a synapse is given by
with,
and
These equations are identical to Prinz et al., 2004 except for the inclusion of a bound for the timescale of activation (we want to avoid the case that as then and ). All other parameters (except ) are identical to Prinz et al., 2004. Following Prinz et al., 2004, we set and for glutamatergic synapses, and and for cholinergic synapses. We set and for both synapse types.
Temperature effects were included in this model as done previously by Caplan et al., 2014 and O'Leary and Marder, 2016. Temperature dependence was introduced in the time constants of the channel-gating variables and , the maximal conductances , the time constants of calcium buffering , and the maximal conductances and time constants of the synapses. This was done by replacing all conductances by and all time scales by where
Here is the temperature, is the reference temperature and is defined as the fold change per 10°C from the reference temperature ( indicates the process type). Finally, the calcium reversal potential depends on temperature through the Nernst equation.
The total dimension of the model is 46 = 3 × 13 + 7 with 13 state variables per cell, and 7 variables for the state of the synapses. In this work some parameters were considered fixed while others were allowed to take values in a range. The number of parameters we varied per cell is 9: the 8 maximal conductances plus the calcium time constant . The temperature sensitivities were assumed to be the same for all cells and amount to 24 additional parameters. We need to specify 8 values for the intrinsic maximal conductances (one per channel type), plus 11 values for the time scales of each intrinsic process (not all currents inactivate), plus the of the calcium time constant. The number of parameters for the synapses is 4:2 for the maximal conductances (of each synapse type), and 2 for the timescales of activation. The models were simulated using a Runge-Kutta order 4 (RK4) method with a time step of (Press et al., 1988). We used the same set of initial conditions for all simulations in this work , and .
Finding parameters
Request a detailed protocolEach model is specified by the maximal conductances and calcium time constants of each cell (9 × 3 parameters), and the temperature sensitivities of each process (24 parameters). We recently introduced a function that upon minimization, results in values of the maximal conductances for which the activity of a single compartment corresponds to periodic bursting with a target frequency () and duty cycle (). The function uses thresholds to obtain temporal information of the waveform, such as spike times, and then uses it to assign a score or error that measures how close the solutions are to a target activity. This function is described in detail in Alonso and Marder, 2019 and was used here to find temperature robust networks.
For any given set of conductances we simulated the model for 20 s and dropped the first 10 s to minimize the effects of transient activity. We then computed the average () burst frequency , the average duty cycle , the number of crossings with a slow wave threshold , the number of bursts , and the average lags between bursts and . To discard unstable solutions we checked if the standard deviation of the burst frequency and duty cycle distributions was small; a solution was discarded if or . If a solution is not discarded we can use these quantities to measure how close it is to a target behavior,
Here measures the mismatch of the bursting frequency of each cell with a target frequency and accounts for the duty cycle. measures the difference between the number of bursts and the number of crossings with the slow wave threshold (we ask that ). compares the lags between bursts (in units of the bursting period ) to a target lag . These measures are discussed in more detail in Alonso and Marder, 2019.
Let denote a set of maximal conductances (and ), we can then define an objective function
where the coefficients define how the different sources of penalties are weighted. In this work we used , , , . These weights were found by trial and error and kept fixed. The penalties were calculated using secs with msecs. The target control behavior was defined as having all cells bursting with 20% duty cycle for () and 25% for the LP and PY cells (). The lag between bursts was targeted to be and . The target burst frequency of all cells was set to (Bucher et al., 2006; Tang et al., 2010; Hamood et al., 2015).
Minimization of the objective function produces sets of maximal conductances for which the resulting circuit activity mimics the pyloric rhythm. The minimization was performed over a search space of allowed values listed here: for each cell we searched for ([]), , , , , , , , (msec). All synaptic conductances were searched in the range ([]). The choice of search space contains our knowledge that some conductances like and are significantly larger than or the synaptic conductances. This choice makes optimizations faster and produces diverse solutions. We discretized our search space by taking 103 equally spaced values for each parameter. We minimized the objective function (Equation 6) using a custom genetic algorithm (Holland, 1992) on a desktop computer with 32 cores. For each network, we ran the algorithm using 1000 individuals taken at random in the search space for 10,000 generations. Finding a successful population takes several hours.
Temperature robustness is achieved in our models by searching sets of values that produce the target pyloric rhythm at each temperature. As in Caplan et al., 2014 and O'Leary and Marder, 2016 we searched for values of between 1 and 2 for conductances and between 1 and 4 for activation/inactivation time scales. We built a new objective function by evaluating the previous objective function (Equation 6) over a set of control temperatures taken at 15°C, 20°C, 25°C and 35°C,
Here, is the score of the solutions of the set of conductances at temperature and is the total number of spikes in all three cells. This last term is introduced to enforce that when the temperature is close to 35°C the network stops working as in experiments. Evaluation of the objective functions (Equation 6) and (Equation 7) requires that the model is simulated for a number of seconds and this is the part of the procedure that requires most computing power. Longer simulations will provide better estimations for the burst frequency and duty cycle of the solutions, but it will linearly increase the time it takes to evaluate them. If the simulations are shorter, evaluations of the objective function are faster but its minimization may be more difficult due to transient behaviors, and its minima may not correspond to stable pyloric rhythms.
Because we only target the activity at three temperature values, it is not guaranteed that the models will behave properly for temperatures in between. We found that if the system meets the target activity at the control temperatures, in the vast majority of the cases, the models also displayed the correct activity for temperatures in between, but there were exceptions (2 in 40 cases). This means that even if a solution achieves low scores, we are still required to screen the temperature values in between to make sure it is indeed temperature compensated over the full range. Increasing the number of control temperatures makes the estimation procedure slower and we met a good balance with 4 control temperatures. All models were subject to the same temperature perturbations and temperature compensation was confirmed in all cases by recording spikes and phases of bursting (Figures 1 and 3), membrane potential (Figure 2), and currents’ dynamics (Figures 7, 8 and 9).
Membrane potential and current shares distributions
Request a detailed protocolThe membrane potential distributions were computed for 101 values of temperature between 10°C and 25°C. For each temperature, the models were simulated from identical initial conditions for 30 s with high numerical resolution (). The distributions were computed using the last 10 s of each simulation. The number of samples of the distributions at each temperature is 5 × 105.
Classification scheme
Request a detailed protocolWe classified the activity of each cell into several categories by inspecting their spiking patterns and features of their waveforms. For this we first compute the spike times and their ISI distributions. If the coefficient of variation , we declare the activity as spiking. We then measure the lag between the spike onset and the crossing times with a slow wave threshold at . If this we label the trace as tonic spiking and else, we label it as single-spike bursting. If then we ask if the cell is bursting. For this we group spikes into bursts using a temporal threshold of 200 ms and compute the distributions of duty cycle and instantaneous burst frequency . If , and we label the trace as regular bursting. If and either or we label the trace as irregular bursting. If we label the trace as irregular spiker. If none of these conditions is met we label the trace as other and finally, if the total number of ISI values is not greater than 1 we label the trace as quiescent. The network state was classified as triphasic or not triphasic. We declared the network (NET) activity as triphasic if the frequency of the cells was similar (within 5%) and they fired in the right order, and not triphasic elsewhere.
Parameters used in this study
Request a detailed protocolEach model was assigned a six-character name. Here, we list which model was used in each figure. Code to simulate the network and the parameters of each model are supplemental to this work.
Figure 1 4 × 6W66. Figure 1—figure supplement 1: from left to right, WWZ3CN, FGFKQS, 71G6LA. Figure 1—figure supplement 2: PLKTEM. Figure 1—figure supplement 3: all models. Figure 2: G7P2WE. Figure 2—figure supplement 1: from left to right, G7P2WE, 4 × 6W66, MAXTTP. Figure 3: 4 × 6W66, SJR46Y.Figure 4: WWZ3CN. Figure 5: FGFKQS.Figure 6: FGFKQS, 71G6LA. Figure 6—figure supplement 1: in lexicographical order, B7SDOL, J9U8SQ, 4R2719, R3VIEW, Q8FBN1, RMZ8KZ, EF84RN, SJR46Y. Figure 7: (top to bottom), J9U8SQ, G81PUL, G7EOZ8, SAXT8Y, KNMLBC. Figure 8: OSHNMD, H5AIZR. Figure 9: WWZ3CN, FGFKQS. Figure 10: 4 × 6W66.
Data availability
Code and parameters to reproduce the results in this work was uploaded to Dryad, http://doi.org/10.5061/dryad.m31b1h7.
-
Dryad Digital RepositoryData from: Temperature compensation in a small rhythmic circuit.https://doi.org/10.5061/dryad.m31b1h7
References
-
Effects of temperature on identified central neurons that control jumping in the grasshopperThe Journal of Neuroscience 2:1538–1553.https://doi.org/10.1523/JNEUROSCI.02-11-01538.1982
-
Mathematical model of an identified stomatogastric ganglion neuronJournal of Neurophysiology 67:332–340.https://doi.org/10.1152/jn.1992.67.2.332
-
Many parameter sets in a multicompartment model oscillator are robust to temperature perturbationsThe Journal of Neuroscience 34:4963–4975.https://doi.org/10.1523/JNEUROSCI.0280-14.2014
-
Neural repetitive firing: a comparative study of membrane properties of crustacean walking leg axonsJournal of Neurophysiology 38:922–932.https://doi.org/10.1152/jn.1975.38.4.922
-
Coexistence of tonic spiking oscillations in a leech neuron modelJournal of Computational Neuroscience 18:255–263.https://doi.org/10.1007/s10827-005-0354-7
-
Theoretical effect of temperature on threshold in the Hodgkin-Huxley nerve modelThe Journal of General Physiology 49:989–1005.https://doi.org/10.1085/jgp.49.5.989
-
Functional consequences of animal-to-animal variation in circuit parametersNature Neuroscience 12:1424–1430.https://doi.org/10.1038/nn.2404
-
Global structure, robustness, and modulation of neuronal modelsThe Journal of Neuroscience 21:5229–5238.https://doi.org/10.1523/JNEUROSCI.21-14-05229.2001
-
Contribution of individual ionic currents to activity of a model stomatogastric ganglion neuronJournal of Neurophysiology 67:341–349.https://doi.org/10.1152/jn.1992.67.2.341
-
Network stability from activity-dependent regulation of neuronal conductancesNeural Computation 11:1079–1096.https://doi.org/10.1162/089976699300016359
-
Neuronal homeostasis: voltage brings it all togetherCurrent Biology 29:R641–R644.https://doi.org/10.1016/j.cub.2019.05.029
-
Ionic currents of the lateral pyloric neuron of the stomatogastric ganglion of the crabJournal of Neurophysiology 67:318–331.https://doi.org/10.1152/jn.1992.67.2.318
-
Effect of temperature on the potential and current thresholds of axon membraneThe Journal of General Physiology 46:257–266.https://doi.org/10.1085/jgp.46.2.257
-
Temperature characteristics of excitation in space-clamped squid axonsThe Journal of General Physiology 49:1007–1018.https://doi.org/10.1085/jgp.49.5.1007
-
Aspects of the biology of the jonah crab, Cancer borealis Stimpson, 1859 in the mid-Atlantic BightJournal of Natural History 11:303–320.https://doi.org/10.1080/00222937700770221
-
The effects of temperature on the threshold of identified neurons in the locustJournal of Comparative Physiology ? A 117:163–182.https://doi.org/10.1007/BF00612785
-
Phase maintenance in the pyloric pattern of the lobster (Panulirus interruptus) stomatogastric ganglionJournal of Computational Neuroscience 4:191–205.https://doi.org/10.1023/a:1008822218061
-
Temperature effects on resting potential and spike parameters of cat motoneuronsExperimental Brain Research 19:478–492.https://doi.org/10.1007/BF00236112
-
Precision of archerfish C-starts is fully temperature compensatedJournal of Experimental Biology 216:3450–3460.https://doi.org/10.1242/jeb.088856
-
Properties of sodium and potassium channels of the squid giant axon far below 0 degrees CThe Journal of Membrane Biology 68:151–160.https://doi.org/10.1007/BF01872261
-
Bistability and its regulation by serotonin in the endogenously bursting neuron R15 in AplysiaJournal of Neurophysiology 75:957–962.https://doi.org/10.1152/jn.1996.75.2.957
-
A model neuron with activity-dependent conductances regulated by multiple calcium sensorsThe Journal of Neuroscience 18:2309–2320.https://doi.org/10.1523/JNEUROSCI.18-07-02309.1998
-
Temporal dynamics of graded synaptic transmission in the lobster stomatogastric ganglionThe Journal of Neuroscience 17:5610–5621.https://doi.org/10.1523/JNEUROSCI.17-14-05610.1997
-
Understanding circuit dynamics using the stomatogastric nervous system of lobsters and crabsAnnual Review of Physiology 69:291–316.https://doi.org/10.1146/annurev.physiol.69.031905.161516
-
High prevalence of multistability of rest states and bursting in a database of a model neuronPLOS Computational Biology 9:e1002930.https://doi.org/10.1371/journal.pcbi.1002930
-
Simpler networksAnnals of the New York Academy of Sciences 193:59–72.https://doi.org/10.1111/j.1749-6632.1972.tb27823.x
-
The role of short-term synaptic dynamics in motor controlCurrent Opinion in Neurobiology 10:683–690.https://doi.org/10.1016/S0959-4388(00)00159-8
-
Noisy inputs and the induction of on-off switching behavior in a neuronal pacemakerJournal of Neurophysiology 96:3338–3348.https://doi.org/10.1152/jn.00486.2006
-
Similar network activity from disparate circuit parametersNature Neuroscience 7:1345–1352.https://doi.org/10.1038/nn1352
-
A kinetic map of the homomeric Voltage-Gated potassium channel (Kv) FamilyFrontiers in Cellular Neuroscience 13:358.https://doi.org/10.3389/fncel.2019.00358
-
The effects of temperature on the stability of a neuronal oscillatorPLOS Computational Biology 9:e1002857.https://doi.org/10.1371/journal.pcbi.1002857
-
Temperature and neuronal circuit function: compensation, tuning and toleranceCurrent Opinion in Neurobiology 22:724–734.https://doi.org/10.1016/j.conb.2012.01.008
-
The emerging role of RNA editing in plasticityJournal of Experimental Biology 218:1812–1821.https://doi.org/10.1242/jeb.119065
-
Effects of temperature on slow and fast inactivation of rat skeletal muscle na + channelsAmerican Journal of Physiology-Cell Physiology 277:C937–C947.https://doi.org/10.1152/ajpcell.1999.277.5.C937
-
Temperature dependence of the ionic current kinetics of Myxicola giant axonsThe Journal of Physiology 235:197–205.https://doi.org/10.1113/jphysiol.1973.sp010384
-
Heat-shock protein protectionTrends in Neurosciences 22:97–99.https://doi.org/10.1016/S0166-2236(98)01392-7
-
Oscillatory behavior of the squid axon membrane potentialThe Journal of General Physiology 42:39–47.https://doi.org/10.1085/jgp.42.1.39
-
Phase maintenance in a rhythmic motor pattern during temperature changes in vivoJournal of Neurophysiology 111:2603–2613.https://doi.org/10.1152/jn.00906.2013
-
Computational model of electrically coupled, intrinsically distinct pacemaker neuronsJournal of Neurophysiology 94:590–604.https://doi.org/10.1152/jn.00013.2005
-
Distribution and abundance of four brachyuran crabs on the northwest Atlantic shelfFishery Bulletin 89:473–492.
-
Robustness of a rhythmic circuit to short- and long-term temperature changesJournal of Neuroscience 32:10075–10085.https://doi.org/10.1523/JNEUROSCI.1443-12.2012
-
Ionic current correlations are ubiquitous across phylaScientific Reports 9:1687.https://doi.org/10.1038/s41598-018-38405-6
-
Selective regulation of current densities underlies spontaneous changes in the activity of cultured neuronsThe Journal of Neuroscience 15:3640–3652.https://doi.org/10.1523/JNEUROSCI.15-05-03640.1995
Article and author information
Author details
Funding
National Institutes of Health (T32 NS07292)
- Leandro M Alonso
National Institutes of Health (MH046742)
- Eve Marder
National Institutes of Health (R35 NS097343)
- Eve Marder
Swartz Foundation (Swartz Foundation 2017)
- Leandro M Alonso
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
LMA acknowledges useful conversations with Hillary Rodgers.
Copyright
© 2020, Alonso and Marder
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,011
- views
-
- 331
- downloads
-
- 53
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Medicine
- Neuroscience
Background:
Alcohol use disorder (AUD) is a global health problem with limited therapeutic options. The biochemical mechanisms that lead to this disorder are not yet fully understood, and in this respect, metabolomics represents a promising approach to decipher metabolic events related to AUD. The plasma metabolome contains a plethora of bioactive molecules that reflects the functional changes in host metabolism but also the impact of the gut microbiome and nutritional habits.
Methods:
In this study, we investigated the impact of severe AUD (sAUD), and of a 3-week period of alcohol abstinence, on the blood metabolome (non-targeted LC-MS metabolomics analysis) in 96 sAUD patients hospitalized for alcohol withdrawal.
Results:
We found that the plasma levels of different lipids ((lyso)phosphatidylcholines, long-chain fatty acids), short-chain fatty acids (i.e. 3-hydroxyvaleric acid) and bile acids were altered in sAUD patients. In addition, several microbial metabolites, including indole-3-propionic acid, p-cresol sulfate, hippuric acid, pyrocatechol sulfate, and metabolites belonging to xanthine class (paraxanthine, theobromine and theophylline) were sensitive to alcohol exposure and alcohol withdrawal. 3-Hydroxyvaleric acid, caffeine metabolites (theobromine, paraxanthine, and theophylline) and microbial metabolites (hippuric acid and pyrocatechol sulfate) were correlated with anxiety, depression and alcohol craving. Metabolomics analysis in postmortem samples of frontal cortex and cerebrospinal fluid of those consuming a high level of alcohol revealed that those metabolites can be found also in brain tissue.
Conclusions:
Our data allow the identification of neuroactive metabolites, from interactions between food components and microbiota, which may represent new targets arising in the management of neuropsychiatric diseases such as sAUD.
Funding:
Gut2Behave project was initiated from ERA-NET NEURON network (Joint Transnational Call 2019) and was financed by Academy of Finland, French National Research Agency (ANR-19-NEUR-0003-03) and the Fonds de la Recherche Scientifique (FRS-FNRS; PINT-MULTI R.8013.19, Belgium). Metabolomics analysis of the TSDS samples was supported by grant from the Finnish Foundation for Alcohol Studies.
-
- Neuroscience
Spinal cord interneurons play critical roles shaping motor output, but their precise identity and connectivity remain unclear. Focusing on the V1 interneuron cardinal class we defined four major V1 subsets in the mouse according to neurogenesis, genetic lineage-tracing, synaptic output to motoneurons, and synaptic inputs from muscle afferents. Sequential neurogenesis delineates different V1 subsets: two early born (Renshaw and Pou6f2) and two late born (Foxp2 and Sp8). Early born Renshaw cells and late born Foxp2-V1 interneurons are tightly coupled to motoneurons, while early born Pou6f2-V1 and late born Sp8-V1 interneurons are not, indicating that timing of neurogenesis does not correlate with motoneuron targeting. V1 clades also differ in cell numbers and diversity. Lineage labeling shows that the Foxp2-V1 clade contains over half of all V1 interneurons, provides the largest inhibitory input to motoneuron cell bodies, and includes subgroups that differ in birthdate, location, and proprioceptive input. Notably, one Foxp2-V1 subgroup, defined by postnatal Otp expression, is positioned near the LMC and receives substantial input from proprioceptors, consistent with an involvement in reciprocal inhibitory pathways. Combined tracing of ankle flexor sensory afferents and interneurons monosynaptically connected to ankle extensors confirmed placement of Foxp2-V1 interneurons in reciprocal inhibitory pathways. Our results validate previously proposed V1 clades as unique functional subtypes that differ in circuit placement, with Foxp2-V1 cells forming the most heterogeneous subgroup. We discuss how V1 organizational diversity enables understanding of their roles in motor control, with implications for their diverse ontogenetic and phylogenetic origins.