Correlations in ion channel mRNA in rhythmically active neurons

Do calcium channel blockers applied to cardiomyocytes cause increased channel expression resulting in reduced efficacy?

In the initial hours following the application of the calcium channel blocker (CCB) nifedipine to microtissues consisting of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), we observe notable variations in the drug's efficacy. Here, we investigate the possibility that these temporal changes in CCB effects are associated with adaptations in the expression of calcium ion channels in cardiomyocyte membranes. To explore this, we employ a recently developed mathematical model that delineates the regulation of calcium ion channel expression by intracellular calcium concentrations. According to the model, a decline in intracellular calcium levels below a certain target level triggers an upregulation of calcium ion channels. Such an upregulation, if instigated by a CCB, would then counteract the drug's inhibitory effect on calcium currents. We assess this hypothesis using time-dependent measurements of hiPSC-CMs dynamics and by refining an existing mathematical model of myocyte action potentials incorporating the dynamic nature of the number of calcium ion channels. The revised model forecasts that the CCB-induced reduction in intracellular calcium concentrations leads to a subsequent increase in calcium ion channel expression, thereby attenuating the drug's overall efficacy. The data and fit models suggest that dynamic changes in cardiac cells in the presence of CCBs may be explainable by induced changes in protein expression, and that this may lead to challenges in understanding calcium based drug effects on the heart unless timings of applications are carefully considered.

Pairwise biosynthesis of ion channels stabilizes excitability and mitigates arrhythmias

Coordinated expression of ion channels is crucial for cardiac rhythms, neural signaling, and cell cycle progression. Perturbation of this balance results in many disorders including cardiac arrhythmias. Prior work revealed association of mRNAs encoding cardiac NaV1.5 (SCN5A) and hERG1 (KCNH2), but the functional significance of this association was not established. Here, we provide a more comprehensive picture of KCNH2, SCN5A, CACNA1C, and KCNQ1 transcripts collectively copurifying with nascent hERG1, NaV1.5, CaV1.2, or KCNQ1 channel proteins. Single-molecule fluorescence in situ hybridization (smFISH) combined with immunofluorescence reveals that the channel proteins are synthesized predominantly as heterotypic pairs from discrete molecules of mRNA, not as larger cotranslational complexes. Puromycin disrupted colocalization of mRNA with its encoded protein, as expected, but remarkably also pairwise mRNA association, suggesting that transcript association relies on intact translational machinery or the presence of the nascent protein. Targeted depletion of KCHN2 by specific shRNA resulted in concomitant reduction of all associated mRNAs, with a corresponding reduction in the encoded channel currents. This co-knockdown effect, originally described for KCNH2 and SCN5A, thus appears to be a general phenomenon among transcripts encoding functionally related proteins. In multielectrode array recordings, proarrhythmic behavior arose when IKr was reduced by the selective blocker dofetilide at IC50 concentrations, but not when equivalent reductions were mediated by shRNA, suggesting that co-knockdown mitigates proarrhythmic behavior expected from the selective reduction of a single channel species. We propose that coordinated, cotranslational association of functionally related ion channel mRNAs confers electrical stability by co-regulating complementary ion channels in macromolecular complexes.

A biophysical perspective on the resilience of neuronal excitability across timescales

Neuronal membrane excitability must be resilient to perturbations that can take place over timescales from milliseconds to months (or even years in long-lived animals). A great deal of attention has been paid to classes of homeostatic mechanisms that contribute to long-term maintenance of neuronal excitability through processes that alter a key structural parameter: the number of ion channel proteins present at the neuronal membrane. However, less attention has been paid to the self-regulating 'automatic' mechanisms that contribute to neuronal resilience by virtue of the kinetic properties of ion channels themselves. Here, we propose that these two sets of mechanisms are complementary instantiations of feedback control, together enabling resilience on a wide range of temporal scales. We further point to several methodological and conceptual challenges entailed in studying these processes - both of which involve enmeshed feedback control loops - and consider the consequences of these mechanisms of resilience.

Comodulation reduces interindividual variability of circuit output

Ionic current levels of identified neurons vary substantially across individual animals. Yet, under similar conditions, neural circuit output can be remarkably similar, as evidenced in many motor systems. All neural circuits are influenced by multiple neuromodulators which provide flexibility to their output. These neuromodulators often overlap in their actions by modulating the same channel type or synapse, yet have neuron-specific actions resulting from distinct receptor expression. Because of this different receptor expression pattern, in the presence of multiple convergent neuromodulators, a common downstream target would be activated more uniformly in circuit neurons across individuals. We therefore propose that a baseline tonic (non-saturating) level of comodulation by convergent neuromodulators can reduce interindividual variability of circuit output. We tested this hypothesis in the pyloric circuit of the crab, Cancer borealis. Multiple excitatory neuropeptides converge to activate the same voltage-gated current in this circuit, but different subsets of pyloric neurons have receptors for each peptide. We quantified the interindividual variability of the unmodulated pyloric circuit output by measuring the activity phases, cycle frequency and intraburst spike number and frequency. We then examined the variability in the presence of different combinations and concentrations of three neuropeptides. We found that at mid-level concentration (30 nM) but not at near-threshold (1 nM) or saturating (1 μM) concentrations, comodulation by multiple neuropeptides reduced the circuit output variability. Notably, the interindividual variability of response properties of an isolated neuron was not reduced by comodulation, suggesting that the reduction of output variability may emerge as a network effect.

Homeostatic regulation of neuronal function: importance of degeneracy and pleiotropy

Neurons maintain their average firing rate and other properties within narrow bounds despite changing conditions. This homeostatic regulation is achieved using negative feedback to adjust ion channel expression levels. To understand how homeostatic regulation of excitability normally works and how it goes awry, one must consider the various ion channels involved as well as the other regulated properties impacted by adjusting those channels when regulating excitability. This raises issues of degeneracy and pleiotropy. Degeneracy refers to disparate solutions conveying equivalent function (e.g., different channel combinations yielding equivalent excitability). This many-to-one mapping contrasts the one-to-many mapping described by pleiotropy (e.g., one channel affecting multiple properties). Degeneracy facilitates homeostatic regulation by enabling a disturbance to be offset by compensatory changes in any one of several different channels or combinations thereof. Pleiotropy complicates homeostatic regulation because compensatory changes intended to regulate one property may inadvertently disrupt other properties. Co-regulating multiple properties by adjusting pleiotropic channels requires greater degeneracy than regulating one property in isolation and, by extension, can fail for additional reasons such as solutions for each property being incompatible with one another. Problems also arise if a perturbation is too strong and/or negative feedback is too weak, or because the set point is disturbed. Delineating feedback loops and their interactions provides valuable insight into how homeostatic regulation might fail. Insofar as different failure modes require distinct interventions to restore homeostasis, deeper understanding of homeostatic regulation and its pathological disruption may reveal more effective treatments for chronic neurological disorders like neuropathic pain and epilepsy.

Electrophysiological heterogeneity in large populations of rabbit ventricular cardiomyocytes

AIMS

Cardiac electrophysiological heterogeneity includes: (i) regional differences in action potential (AP) waveform, (ii) AP waveform differences in cells isolated from a single region, (iii) variability of the contribution of individual ion currents in cells with similar AP durations (APDs). The aim of this study is to assess intra-regional AP waveform differences, to quantify the contribution of specific ion channels to the APD via drug responses and to generate a population of mathematical models to investigate the mechanisms underlying heterogeneity in rabbit ventricular cells.

METHODS AND RESULTS

APD in ∼50 isolated cells from subregions of the LV free wall of rabbit hearts were measured using a voltage-sensitive dye. When stimulated at 2 Hz, average APD90 value in cells from the basal epicardial region was 254 ± 25 ms (mean ± standard deviation) in 17 hearts with a mean interquartile range (IQR) of 53 ± 17 ms. Endo-epicardial and apical-basal APD90 differences accounted for ∼10% of the IQR value. Highly variable changes in APD occurred after IK(r) or ICa(L) block that included a sub-population of cells (HR) with an exaggerated (hyper) response to IK(r) inhibition. A set of 4471 AP models matching the experimental APD90 distribution was generated from a larger population of models created by random variation of the maximum conductances (Gmax) of 8 key ion channels/exchangers/pumps. This set reproduced the pattern of cell-specific responses to ICa(L) and IK(r) block, including the HR sub-population. The models exhibited a wide range of Gmax values with constrained relationships linking ICa(L) with IK(r), ICl, INCX, and INaK.

CONCLUSION

Modelling the measured range of inter-cell APDs required a larger range of key Gmax values indicating that ventricular tissue has considerable inter-cell variation in channel/pump/exchanger activity. AP morphology is retained by relationships linking specific ionic conductances. These interrelationships are necessary for stable repolarization despite large inter-cell variation of individual conductances and this explains the variable sensitivity to ion channel block.

Neuromodulation reduces interindividual variability of neuronal output

In similar states, neural circuits produce similar outputs across individuals despite substantial interindividual variability in neuronal ionic conductances and synapses. Circuit states are largely shaped by neuromodulators that tune ionic conductances. It is therefore possible that, in addition to producing flexible circuit output, neuromodulators also contribute to output similarity despite varying ion channel expression. We studied whether neuromodulation at saturating concentrations can increase the output similarity of a single identified neuron across individual animals. Using the LP neuron of the crab stomatogastric ganglion (STG), we compared the variability of f-I curves and rebound properties in the presence of neuropeptides. The two neuropeptides we used converge to activate the same target current, which increases neuronal excitability. Output variability was lower in the presence of the neuropeptides, regardless of whether the neuropeptides significantly changed the mean of the corresponding parameter or not. However, the addition of the second neuropeptide did not add further to the reduction of variability. With a family of computational LP-like models, we explored how increased excitability and target variability contribute to output similarity and found two mechanisms: Saturation of the responses and a differential increase in baseline activity. Saturation alone can reduce the interindividual variability only if the population shares a similar ceiling for the responses. In contrast, reduction of variability due to the increase in baseline activity is independent of ceiling effects.Significance StatementThe activity of single neurons and neural circuits can be very similar across individuals even though the ionic currents underlying activity are variable. The mechanisms that compensate for the underlying variability and promote output similarity are poorly understood but may involve neuromodulation. Using an identified neuron, we show that neuropeptide modulation of excitability can reduce interindividual variability of response properties at a single-neuron level in two ways. First, the neuropeptide increases baseline excitability in a differential manner, resulting in similar response thresholds. Second, the neuropeptide increases excitability towards a shared saturation level, promoting similar maximal firing rates across individuals. Such tuning of neuronal excitability could be an important mechanism compensating for interindividual variability of ion channel expression.

Mapping circuit dynamics during function and dysfunction

Neural circuits can generate many spike patterns, but only some are functional. The study of how circuits generate and maintain functional dynamics is hindered by a poverty of description of circuit dynamics across functional and dysfunctional states. For example, although the regular oscillation of a central pattern generator is well characterized by its frequency and the phase relationships between its neurons, these metrics are ineffective descriptors of the irregular and aperiodic dynamics that circuits can generate under perturbation or in disease states. By recording the circuit dynamics of the well-studied pyloric circuit in Cancer borealis, we used statistical features of spike times from neurons in the circuit to visualize the spike patterns generated by this circuit under a variety of conditions. This approach captures both the variability of functional rhythms and the diversity of atypical dynamics in a single map. Clusters in the map identify qualitatively different spike patterns hinting at different dynamic states in the circuit. State probability and the statistics of the transitions between states varied with environmental perturbations, removal of descending neuromodulatory inputs, and the addition of exogenous neuromodulators. This analysis reveals strong mechanistically interpretable links between complex changes in the collective behavior of a neural circuit and specific experimental manipulations, and can constrain hypotheses of how circuits generate functional dynamics despite variability in circuit architecture and environmental perturbations.

Minimal requirements for a neuron to co-regulate many properties and the implications for ion channel correlations and robustness

Neurons regulate their excitability by adjusting their ion channel levels. Degeneracy – achieving equivalent outcomes (excitability) using different solutions (channel combinations) – facilitates this regulation by enabling a disruptive change in one channel to be offset by compensatory changes in other channels. But neurons must coregulate many properties. Pleiotropy – the impact of one channel on more than one property – complicates regulation because a compensatory ion channel change that restores one property to its target value often disrupts other properties. How then does a neuron simultaneously regulate multiple properties? Here, we demonstrate that of the many channel combinations producing the target value for one property (the single-output solution set), few combinations produce the target value for other properties. Combinations producing the target value for two or more properties (the multioutput solution set) correspond to the intersection between single-output solution sets. Properties can be effectively coregulated only if the number of adjustable channels (n_in) exceeds the number of regulated properties (n_out). Ion channel correlations emerge during homeostatic regulation when the dimensionality of solution space (n_in − n_out) is low. Even if each property can be regulated to its target value when considered in isolation, regulation as a whole fails if single-output solution sets do not intersect. Our results also highlight that ion channels must be coadjusted with different ratios to regulate different properties, which suggests that each error signal drives modulatory changes independently, despite those changes ultimately affecting the same ion channels.

Reciprocally inhibitory circuits operating with distinct mechanisms are differently robust to perturbation and modulation

Reciprocal inhibition is a building block in many sensory and motor circuits. We studied the features that underly robustness in reciprocally inhibitory two neuron circuits. We used the dynamic clamp to create reciprocally inhibitory circuits from pharmacologically isolated neurons of the crab stomatogastric ganglion by injecting artificial graded synaptic (I_Syn) and hyperpolarization-activated inward (I_H) currents. There is a continuum of mechanisms in circuits that generate antiphase oscillations, with ‘release’ and ‘escape’ mechanisms at the extremes, and mixed mode oscillations between these extremes. In release, the active neuron primarily controls the off/on transitions. In escape, the inhibited neuron controls the transitions. We characterized the robustness of escape and release circuits to alterations in circuit parameters, temperature, and neuromodulation. We found that escape circuits rely on tight correlations between synaptic and H conductances to generate bursting but are resilient to temperature increase. Release circuits are robust to variations in synaptic and H conductances but fragile to temperature increase. The modulatory current (I_MI) restores oscillations in release circuits but has little effect in escape circuits. Perturbations can alter the balance of escape and release mechanisms and can create mixed mode oscillations. We conclude that the same perturbation can have dramatically different effects depending on the circuits’ mechanism of operation that may not be observable from basal circuit activity.

Degeneracy in the nervous system: from neuronal excitability to neural coding

Degeneracy is ubiquitous across biological systems where structurally different elements can yield a similar outcome. Degeneracy is of particular interest in neuroscience too. On the one hand, degeneracy confers robustness to the nervous system and facilitates evolvability: Different elements provide a backup plan for the system in response to any perturbation or disturbance. On the other, a difficulty in the treatment of some neurological disorders such as chronic pain is explained in light of different elements all of which contribute to the pathological behavior of the system. Under these circumstances, targeting a specific element is ineffective because other elements can compensate for this modulation. Understanding degeneracy in the physiological context explains its beneficial role in the robustness of neural circuits. Likewise, understanding degeneracy in the pathological context opens new avenues of discovery to find more effective therapies.

Ion Channel Degeneracy, Variability, and Covariation in Neuron and Circuit Resilience

The large number of ion channels found in all nervous systems poses fundamental questions concerning how the characteristic intrinsic properties of single neurons are determined by the specific subsets of channels they express. All neurons display many different ion channels with overlapping voltage- and time-dependent properties. We speculate that these overlapping properties promote resilience in neuronal function. Individual neurons of the same cell type show variability in ion channel conductance densities even though they can generate reliable and similar behavior. This complicates a simple assignment of function to any conductance and is associated with variable responses of neurons of the same cell type to perturbations, deletions, and pharmacological manipulation. Ion channel genes often show strong positively correlated expression, which may result from the molecular and developmental rules that determine which ion channels are expressed in a given cell type.

Modeling of lamprey reticulospinal neurons: multiple distinct parameter sets yield realistic simulations

For the lamprey and other vertebrates, reticulospinal (RS) neurons project descending axons to the spinal cord and activate motor networks to initiate locomotion and other behaviors. In the present study, a biophysically detailed computer model of lamprey RS neurons was constructed consisting of three compartments: dendritic, somatic, and axon initial segment (AIS). All compartments included passive channels. In addition, the soma and AIS had fast potassium and sodium channels. The soma included three additional voltage-gated ion channels (slow sodium and high- and low-voltage-activated calcium) and calcium-activated potassium channels. An initial manually adjusted default parameter set, which was based, in part, on modified parameters from models of lamprey spinal neurons, generated simulations of single action potentials and repetitive firing that scored favorably (0.658; maximum = 0.964) compared with experimentally derived properties of lamprey RS neurons. Subsequently, a dual-annealing search paradigm identified 4,302 viable parameter sets at local maxima within parameter space that yielded higher scores than the default parameter set, including many with much higher scores of approximately 0.85-0.87 (i.e., ~30% improvement). In addition, 5- and 2-conductance grid searches identified a relatively large number of viable parameters sets for which significant correlations were present between maximum conductances for pairs of ion channels. The present results indicated that multiple model parameter sets ("solutions") generated action potentials and repetitive firing that mimicked many of the properties of lamprey RS neurons. To our knowledge, this is the first study to systematically explore parameter space for a biophysically detailed model of lamprey RS neurons.NEW & NOTEWORTHY A computer model of lamprey reticulospinal neurons with a default parameter set produced simulations of action potentials and repetitive firing that scored favorably compared with the properties of these neurons. A dual-annealing search algorithm explored ~50 million parameter sets and identified 4,302 distinct viable parameter sets within parameter space that yielded higher/much higher scores than the default parameter set. In addition, 5- and 2-conductance grid searches identified significant correlations between maximum conductances for pairs of ion channels.

Circuit Stability to Perturbations Reveals Hidden Variability in the Balance of Intrinsic and Synaptic Conductances

Neurons and circuits each with a distinct balance of intrinsic and synaptic conductances can generate similar behavior but sometimes respond very differently to perturbation. Examining a large family of circuit models with non-identical neurons and synapses underlying rhythmic behavior, we analyzed the circuits' response to modifications in single and multiple intrinsic conductances in the individual neurons. To summarize these changes over the entire range of perturbed parameters, we quantified circuit output by defining a global stability measure. Using this measure, we identified specific subsets of conductances that when perturbed generate similar behavior in diverse individuals of the population. Our unbiased clustering analysis enabled us to quantify circuit stability when simultaneously perturbing multiple conductances as a nonlinear combination of single conductance perturbations. This revealed surprising conductance combinations that can predict the response to specific perturbations, even when the remaining intrinsic and synaptic conductances are unknown. Therefore, our approach can expose hidden variability in the balance of intrinsic and synaptic conductances of the same neurons across different versions of the same circuit solely from the circuit response to perturbations. Developed for a specific family of model circuits, our quantitative approach to characterizing high-dimensional degenerate systems provides a conceptual and analytic framework to guide future theoretical and experimental studies on degeneracy and robustness.

SIGNIFICANCE STATEMENT Neural circuits can generate nearly identical behavior despite neuronal and synaptic parameters varying several-fold between individual instantiations. Yet, when these parameters are perturbed through channel deletions and mutations or environmental disturbances, seemingly identical circuits can respond very differently. What distinguishes inconsequential perturbations that barely alter circuit behavior from disruptive perturbations that drastically disturb circuit output remains unclear. Focusing on a family of rhythmic circuits, we propose a computational approach to reveal hidden variability in the intrinsic and synaptic conductances in seemingly identical circuits based solely on circuit output to different perturbations. We uncover specific conductance combinations that work similarly to maintain stability and predict the effect of changing multiple conductances simultaneously, which often results from neuromodulation or injury.

Neuromodulation of central pattern generators and its role in the functional recovery of central pattern generator activity

Neuromodulators play an important role in how the nervous system organizes activity that results in behavior. Disruption of the normal patterns of neuromodulatory release or production is known to be related to the onset of severe pathologies such as Parkinson's disease, Rett syndrome, Alzheimer's disease, and affective disorders. Some of these pathologies involve neuronal structures that are called central pattern generators (CPGs), which are involved in the production of rhythmic activities throughout the nervous system. Here I discuss the interplay between CPGs and neuromodulatory activity, with particular emphasis on the potential role of neuromodulators in the recovery of disrupted neuronal activity. I refer to invertebrate and vertebrate model systems and some of the lessons we have learned from research on these systems and propose a few avenues for future research. I make one suggestion that may guide future research in the field: neuromodulators restrict the parameter landscape in which CPG components operate, and the removal of neuromodulators may enable a perturbed CPG in finding a new set of parameter values that can allow it to regain normal function.

Dual separable feedback systems govern firing rate homeostasis

Firing rate homeostasis (FRH) stabilizes neural activity. A pervasive and intuitive theory argues that a single variable, calcium, is detected and stabilized through regulatory feedback. A prediction is that ion channel gene mutations with equivalent effects on neuronal excitability should invoke the same homeostatic response. In agreement, we demonstrate robust FRH following either elimination of Kv4/Shal protein or elimination of the Kv4/Shal conductance. However, the underlying homeostatic signaling mechanisms are distinct. Eliminating Shal protein invokes Krüppel-dependent rebalancing of ion channel gene expression including enhanced slo, Shab, and Shaker. By contrast, expression of these genes remains unchanged in animals harboring a CRISPR-engineered, Shal pore-blocking mutation where compensation is achieved by enhanced IK_DR. These different homeostatic processes have distinct effects on homeostatic synaptic plasticity and animal behavior. We propose that FRH includes mechanisms of proteostatic feedback that act in parallel with activity-driven feedback, with implications for the pathophysiology of human channelopathies.

Combined mechanisms of neural firing rate homeostasis

Spikes in the membrane potential of neurons comprise the currency of information processing in the brain. The ability of neurons to convert any information present across their multiple inputs into a significant modification to the pattern of their emitted spikes depends on the rate at which they emit spikes. If the mean rate is near the neuron's maximum, or if the rate is near zero, then changes in the inputs have minimal impact on the neuron's firing rate. Therefore, a neuron needs to control its mean rate. Protocols that either dramatically increase or decrease a neuron's firing rate lead to multiple compensatory changes that return the neuron's mean rate toward its prior value. In this primer, first as a summary of our previous work (Cannon and Miller in J Neurophysiol 116(5):2004-2022, 2016; Cannon and Miller in J Math Neurosci 7(1):1, 2017), we describe the advantages and disadvantages of having more than one such control mechanism responding to the neuron's firing rate. We suggest how problems of two, coexisting, potentially competing mechanisms can be overcome. Key requirements are: (1) the control be of a distribution of values, which the controlled variable achieves over a fast timescale compared to the timescale of the control system; (2) at least one of the control mechanisms be nonlinear; and (3) the two control systems are satisfied by a stable distribution or range of values that can be achieved by the variable. We show examples of functional control systems, including the previously studied integral feedback controller and new simulations of a "bang-bang" controller, that allow for compensation when inputs to the system change. Finally, we present new results describing how the underlying signal processing pathways would produce mechanisms of dual control, as opposed to a single mechanism with two outputs, and compare the responses of these systems to changes of input statistics.

Ionic current correlations are ubiquitous across phyla

Ionic currents, whether measured as conductance amplitude or as ion channel transcript numbers, can vary many-fold within a population of identified neurons. In invertebrate neuronal types multiple currents can be seen to vary while at the same time their magnitudes are correlated. These conductance amplitude correlations are thought to reflect a tight homeostasis of cellular excitability that enhances the robustness and stability of neuronal activity over long stretches of time. Although such ionic conductance correlations are well documented in invertebrates, they have not been reported in vertebrates. Here we demonstrate with two examples, identified mouse hippocampal granule cells (GCs) and cholinergic basal forebrain neurons, that the correlation of ionic conductance amplitudes between different ionic currents also exists in vertebrates, and we argue that it is a ubiquitous phenomenon expressed by many species across phyla. We further demonstrate that in dentate gyrus GCs these conductance correlations are likely regulated in a circadian manner. This is reminiscent of the known conductance regulation by neuromodulators in crustaceans. However, in GCs we observe a more nuanced regulation, where for some conductance pairs the correlations are completely eliminated while for others the correlation is quantitatively modified but not obliterated.

Visualization of currents in neural models with similar behavior and different conductance densities

Conductance-based models of neural activity produce large amounts of data that can be hard to visualize and interpret. We introduce visualization methods to display the dynamics of the ionic currents and to display the models’ response to perturbations. To visualize the currents’ dynamics, we compute the percent contribution of each current and display them over time using stacked-area plots. The waveform of the membrane potential and the contribution of each current change as the models are perturbed. To represent these changes over a range of the perturbation control parameter, we compute and display the distributions of these waveforms. We illustrate these procedures in six examples of bursting model neurons with similar activity but that differ as much as threefold in their conductance densities. These visualization methods provide heuristic insight into why individual neurons or networks with similar behavior can respond widely differently to perturbations.

The nervous system contains networks of neurons that generate electrical signals to communicate with each other and the rest of the body. Such electrical signals are due to the flow of ions into or out of the neurons via proteins known as ion channels. Neurons have many different kinds of ion channels that only allow specific ions to pass. Therefore, for a neuron to produce an electrical signal, the activities of several different ion channels need to be coordinated so that they all open and close at certain times.

Researchers have previously used data collected from various experiments to develop detailed models of electrical signals in neurons. These models incorporate information about how and when the ion channels may open and close, and can produce numerical simulations of the different ionic currents. However, it can be difficult to display the currents and observe how they change when several different ion channels are involved.

Alonso and Marder used simple mathematical concepts to develop new methods to display ionic currents in computational models of neurons. These tools use color to capture changes in ionic currents and provide insights into how the opening and closing of ion channels shape electrical signals.

The methods developed by Alonso and Marder could be adapted to display the behavior of biochemical reactions or other topics in biology and may, therefore, be useful to analyze data generated by computational models of many different types of cells. Additionally, these methods may potentially be used as educational tools to illustrate the coordinated opening and closing of ion channels in neurons and other fundamental principles of neuroscience that are otherwise hard to demonstrate.

Two central pattern generators from the crab, Cancer borealis, respond robustly and differentially to extreme extracellular pH

The activity of neuronal circuits depends on the properties of the constituent neurons and their underlying synaptic and intrinsic currents. We describe the effects of extreme changes in extracellular pH – from pH 5.5 to 10.4 – on two central pattern generating networks, the stomatogastric and cardiac ganglia of the crab, Cancer borealis. Given that the physiological properties of ion channels are known to be sensitive to pH within the range tested, it is surprising that these rhythms generally remained robust from pH 6.1 to pH 8.8. The pH sensitivity of these rhythms was highly variable between animals and, unexpectedly, between ganglia. Animal-to-animal variability was likely a consequence of similar network performance arising from variable sets of underlying conductances. Together, these results illustrate the potential difficulty in generalizing the effects of environmental perturbation across circuits, even within the same animal.

Effects of Chronic Spinal Cord Injury on Relationships among Ion Channel and Receptor mRNAs in Mouse Lumbar Spinal Cord

Spinal cord injury (SCI) causes widespread changes in gene expression of the spinal cord, even in the undamaged spinal cord below the level of the lesion. Less is known about changes in the correlated expression of genes after SCI. We investigated gene co-expression networks among voltage-gated ion channel and neurotransmitter receptor mRNA levels using quantitative RT-PCR in longitudinal slices of the mouse lumbar spinal cord in control and chronic SCI animals. These longitudinal slices were made from the ventral surface of the cord, thus forming slices relatively enriched in motor neurons or interneurons. We performed absolute quantitation of mRNA copy number for 50 ion channel or receptor transcripts from each sample, and used multiple correlation analyses to detect patterns in correlated mRNA levels across all pairs of genes. The majority of channels and receptors changed in expression as a result of chronic SCI, but did so differently across slice levels. Furthermore, motor neuron-enriched slices experienced an overall loss of correlated channel and receptor expression, while interneuron slices showed a dramatic increase in the number of positively correlated transcripts. These correlation profiles suggest that spinal cord injury induces distinct changes across cell types in the organization of gene co-expression networks for ion channels and transmitter receptors.

Neurotransmitter identity and electrophysiological phenotype are genetically coupled in midbrain dopaminergic neurons

Most neuronal types have a well-identified electrical phenotype. It is now admitted that a same phenotype can be produced using multiple biophysical solutions defined by ion channel expression levels. This argues that systems-level approaches are necessary to understand electrical phenotype genesis and stability. Midbrain dopaminergic (DA) neurons, although quite heterogeneous, exhibit a characteristic electrical phenotype. However, the quantitative genetic principles underlying this conserved phenotype remain unknown. Here we investigated the quantitative relationships between ion channels’ gene expression levels in midbrain DA neurons using single-cell microfluidic qPCR. Using multivariate mutual information analysis to decipher high-dimensional statistical dependences, we unravel co-varying gene modules that link neurotransmitter identity and electrical phenotype. We also identify new segregating gene modules underlying the diversity of this neuronal population. We propose that the newly identified genetic coupling between neurotransmitter identity and ion channels may play a homeostatic role in maintaining the electrophysiological phenotype of midbrain DA neurons.

A Heart for Diversity: Simulating Variability in Cardiac Arrhythmia Research

In cardiac electrophysiology, there exist many sources of inter- and intra-personal variability. These include variability in conditions and environment, and genotypic and molecular diversity, including differences in expression and behavior of ion channels and transporters, which lead to phenotypic diversity (e.g., variable integrated responses at the cell, tissue, and organ levels). These variabilities play an important role in progression of heart disease and arrhythmia syndromes and outcomes of therapeutic interventions. Yet, the traditional in silico framework for investigating cardiac arrhythmias is built upon a parameter/property-averaging approach that typically overlooks the physiological diversity. Inspired by work done in genetics and neuroscience, new modeling frameworks of cardiac electrophysiology have been recently developed that take advantage of modern computational capabilities and approaches, and account for the variance in the biological data they are intended to illuminate. In this review, we outline the recent advances in statistical and computational techniques that take into account physiological variability, and move beyond the traditional cardiac model-building scheme that involves averaging over samples from many individuals in the construction of a highly tuned composite model. We discuss how these advanced methods have harnessed the power of big (simulated) data to study the mechanisms of cardiac arrhythmias, with a special emphasis on atrial fibrillation, and improve the assessment of proarrhythmic risk and drug response. The challenges of using in silico approaches with variability are also addressed and future directions are proposed.

Innexin expression in electrically coupled motor circuits

The many roles of innexins, the molecules that form gap junctions in invertebrates, have been explored in numerous species. Here, we present a summary of innexin expression and function in two small, central pattern generating circuits found in crustaceans: the stomatogastric ganglion and the cardiac ganglion. The two ganglia express multiple innexin genes, exhibit varying combinations of symmetrical and rectifying gap junctions, as well as gap junctions within and across different cell types. Past studies have revealed correlations in ion channel and innexin expression in coupled neurons, as well as intriguing functional relationships between ion channel conductances and electrical coupling. Together, these studies suggest a putative role for innexins in correlating activity between coupled neurons at the levels of gene expression and physiological activity during development and in the adult animal.

Homeostatic plasticity of excitability in crustacean central pattern generator networks

Plasticity of excitability can come in two general forms: changes in excitability that alter neuronal output (e.g. long-term potentiation of intrinsic excitability) or excitability changes that stabilize neuronal output (homeostatic plasticity). Here we discuss the latter form of plasticity in the context of the crustacean stomatogastric nervous system, and a second central pattern generator circuit, the cardiac ganglion. We discuss this plasticity at three levels: rapid homeostatic changes in membrane conductance, longer-term effects of neuromodulation on excitability, and the impacts of activity-dependent feedback on steady-state channel mRNA levels. We then conclude with thoughts on the implications of plasticity of excitability for variability of conductance levels across populations of motor neurons.

When complex neuronal structures may not matter

Much work has explored animal-to-animal variability and compensation in ion channel expression. Yet, little is known regarding the physiological consequences of morphological variability. We quantify animal-to-animal variability in cable lengths (CV = 0.4) and branching patterns in the Gastric Mill (GM) neuron, an identified neuron type with highly-conserved physiological properties in the crustacean stomatogastric ganglion (STG) of Cancer borealis. We examined passive GM electrotonic structure by measuring the amplitudes and apparent reversal potentials (E_revs) of inhibitory responses evoked with focal glutamate photo-uncaging in the presence of TTX. Apparent E_revs were relatively invariant across sites (mean CV ± SD = 0.04 ± 0.01; 7–20 sites in each of 10 neurons), which ranged between 100–800 µm from the somatic recording site. Thus, GM neurons are remarkably electrotonically compact (estimated λ > 1.5 mm). Electrotonically compact structures, in consort with graded transmission, provide an elegant solution to observed morphological variability in the STG.

DOI:http://dx.doi.org/10.7554/eLife.23508.001

Deep sequencing of transcriptomes from the nervous systems of two decapod crustaceans to characterize genes important for neural circuit function and modulation

BACKGROUND

Crustaceans have been studied extensively as model systems for nervous system function from single neuron properties to behavior. However, lack of molecular sequence information and tools have slowed the adoption of these physiological systems as molecular model systems. In this study, we sequenced and performed de novo assembly for the nervous system transcriptomes of two decapod crustaceans: the Jonah crab (Cancer borealis) and the American lobster (Homarus americanus).

RESULTS

Forty-two thousand, seven hundred sixty-six and sixty thousand, two hundred seventy-three contigs were assembled from C. borealis and H. americanus respectively, representing 9,489 and 11,061 unique coding sequences. From these transcripts, genes associated with neural function were identified and manually curated to produce a characterization of multiple gene families important for nervous system function. This included genes for 34 distinct ion channel types, 17 biogenic amine and 5 GABA receptors, 28 major transmitter receptor subtypes including glutamate and acetylcholine receptors, and 6 gap junction proteins - the Innexins.

CONCLUSION

With this resource, crustacean model systems are better poised for incorporation of modern genomic and molecular biology technologies to further enhance the interrogation of fundamentals of nervous system function.

Consequences of degeneracy in network function

Often distinct elements serve similar functions within a network. However, it is unclear whether this network degeneracy is beneficial, or merely a reflection of tighter regulation of overall network performance relative to individual neuronal properties. We review circumstances where data strongly suggest that degeneracy is beneficial in that it makes network function more robust. Importantly, network degeneracy is likely to have functional consequences that are not widely appreciated. This is likely to be true when network activity is configured by modulators with persistent actions, and the history of network activity potentially impacts subsequent functioning. Data suggest that degeneracy in this context may be important for the creation of latent memories, and for state-dependent task switching.

Analysis of Family Structures Reveals Robustness or Sensitivity of Bursting Activity to Parameter Variations in a Half-Center Oscillator (HCO) Model

The underlying mechanisms that support robustness in neuronal networks are as yet unknown. However, recent studies provide evidence that neuronal networks are robust to natural variations, modulation, and environmental perturbations of parameters, such as maximal conductances of intrinsic membrane and synaptic currents. Here we sought a method for assessing robustness, which might easily be applied to large brute-force databases of model instances. Starting with groups of instances with appropriate activity (e.g., tonic spiking), our method classifies instances into much smaller subgroups, called families, in which all members vary only by the one parameter that defines the family. By analyzing the structures of families, we developed measures of robustness for activity type. Then, we applied these measures to our previously developed model database, HCO-db, of a two-neuron half-center oscillator (HCO), a neuronal microcircuit from the leech heartbeat central pattern generator where the appropriate activity type is alternating bursting. In HCO-db, the maximal conductances of five intrinsic and two synaptic currents were varied over eight values (leak reversal potential also varied, five values). We focused on how variations of particular conductance parameters maintain normal alternating bursting activity while still allowing for functional modulation of period and spike frequency. We explored the trade-off between robustness of activity type and desirable change in activity characteristics when intrinsic conductances are altered and identified the hyperpolarization-activated (h) current as an ideal target for modulation. We also identified ensembles of model instances that closely approximate physiological activity and can be used in future modeling studies.

There and back again: Iterating between population-based modeling and experiments reveals surprising regulation of calcium transients in rat cardiac myocytes

While many ion channels and transporters involved in cardiac cellular physiology have been identified and described, the relative importance of each in determining emergent cellular behaviors remains unclear. Here we address this issue with a novel approach that combines population-based mathematical modeling with experimental tests to systematically quantify the relative contributions of different ion channels and transporters to the amplitude of the cellular Ca(2+) transient. Sensitivity analysis of a mathematical model of the rat ventricular cardiomyocyte quantified the response of cell behaviors to changes in the level of each ion channel and transporter, and experimental tests of these predictions were performed to validate or invalidate the predictions. The model analysis found that partial inhibition of the transient outward current in rat ventricular epicardial myocytes was predicted to have a greater impact on Ca(2+) transient amplitude than either: (1) inhibition of the same current in endocardial myocytes, or (2) comparable inhibition of the sarco/endoplasmic reticulum Ca(2+) ATPase (SERCA). Experimental tests confirmed the model predictions qualitatively but showed some quantitative disagreement. This guided us to recalibrate the model by adjusting the relative importance of several Ca(2+) fluxes, thereby improving the consistency with experimental data and producing a more predictive model. Analysis of human cardiomyocyte models suggests that the relative importance of outward currents to Ca(2+) transporters is generalizable to human atrial cardiomyocytes, but not ventricular cardiomyocytes. Overall, our novel approach of combining population-based mathematical modeling with experimental tests has yielded new insight into the relative importance of different determinants of cell behavior.

Molecular variability elicits a tunable switch with discrete neuromodulatory response phenotypes

Recent single cell studies show extensive molecular variability underlying cellular responses. We evaluated the impact of molecular variability in the expression of cell signaling components and ion channels on electrophysiological excitability and neuromodulation. We employed a computational approach that integrated neuropeptide receptor-mediated signaling with electrophysiology. We simulated a population of neurons in which expression levels of a neuropeptide receptor and multiple ion channels were simultaneously varied within a physiological range. We analyzed the effects of variation on the electrophysiological response to a neuropeptide stimulus. Our results revealed distinct response patterns associated with low versus high receptor levels. Neurons with low receptor levels showed increased excitability and neurons with high receptor levels showed reduced excitability. These response patterns were separated by a narrow receptor level range forming a separatrix. The position of this separatrix was dependent on the expression levels of multiple ion channels. To assess the relative contributions of receptor and ion channel levels to the response profiles, we categorized the responses into six phenotypes based on response kinetics and magnitude. We applied several multivariate statistical approaches and found that receptor and channel expression levels influence the neuromodulation response phenotype through a complex though systematic mapping. Our analyses extended our understanding of how cellular responses to neuromodulation vary as a function of molecular expression. Our study showed that receptor expression and biophysical state interact with distinct relative contributions to neuronal excitability.

De novo transcriptome assembly for the lobster Homarus americanus and characterization of differential gene expression across nervous system tissues

BACKGROUND

The American lobster, Homarus americanus, is an important species as an economically valuable fishery, a key member in marine ecosystems, and a well-studied model for central pattern generation, the neural networks that control rhythmic motor patterns. Despite multi-faceted scientific interest in this species, currently our genetic resources for the lobster are limited. In this study, we de novo assemble a transcriptome for Homarus americanus using central nervous system (CNS), muscle, and hybrid neurosecretory tissues and compare gene expression across these tissue types. In particular, we focus our analysis on genes relevant to central pattern generation and the identity of the neurons in a neural network, which is defined by combinations of genes distinguishing the neuronal behavior and phenotype, including ion channels, neurotransmitters, neuromodulators, receptors, transcription factors, and other gene products.

RESULTS

Using samples from the central nervous system (brain, abdominal ganglia), abdominal muscle, and heart (cardiac ganglia, pericardial organs, muscle), we used RNA-Seq to characterize gene expression patterns across tissues types. We also compared control tissues with those challenged with the neuropeptide proctolin in vivo. Our transcriptome generated 34,813 transcripts with known protein annotations. Of these, 5,000-10,000 of annotated transcripts were significantly differentially expressed (DE) across tissue types. We found 421 transcripts for ion channels and identified receptors and/or proteins for over 20 different neurotransmitters and neuromodulators. Results indicated tissue-specific expression of select neuromodulator (allostatin, myomodulin, octopamine, nitric oxide) and neurotransmitter (glutamate, acetylcholine) pathways. We also identify differential expression of ion channel families, including kainite family glutamate receptors, inward-rectifying K(+) (IRK) channels, and transient receptor potential (TRP) A family channels, across central pattern generating tissues.

CONCLUSIONS

Our transcriptome-wide profiles of the rhythmic pattern generating abdominal and cardiac nervous systems in Homarus americanus reveal candidates for neuronal features that drive the production of motor output in these systems.

Monoaminergic tone supports conductance correlations and stabilizes activity features in pattern generating neurons of the lobster, Panulirus interruptus

Experimental and computational studies demonstrate that different sets of intrinsic and synaptic conductances can give rise to equivalent activity patterns. This is because the balance of conductances, not their absolute values, defines a given activity feature. Activity-dependent feedback mechanisms maintain neuronal conductance correlations and their corresponding activity features. This study demonstrates that tonic nM concentrations of monoamines enable slow, activity-dependent processes that can maintain a correlation between the transient potassium current (I_A) and the hyperpolarization activated current (I_h) over the long-term (i.e., regulatory change persists for hours after removal of modulator). Tonic 5 nM DA acted through an RNA interference silencing complex (RISC)- and RNA polymerase II-dependent mechanism to maintain a long-term positive correlation between I_A and I_h in the lateral pyloric neuron (LP) but not in the pyloric dilator neuron (PD). In contrast, tonic 5 nM 5HT maintained a RISC-dependent positive correlation between I_A and I_h in PD but not LP over the long-term. Tonic 5 nM OCT maintained a long-term negative correlation between I_A and I_h in PD but not LP; however, it was only revealed when RISC was inhibited. This study also demonstrated that monoaminergic tone can also preserve activity features over the long-term: the timing of LP activity, LP duty cycle and LP spike number per burst were maintained by tonic 5 nM DA. The data suggest that low-level monoaminergic tone acts through multiple slow processes to permit cell-specific, activity-dependent regulation of ionic conductances to maintain conductance correlations and their corresponding activity features over the long-term.

Variability in State-Dependent Plasticity of Intrinsic Properties during Cell-Autonomous Self-Regulation of Calcium Homeostasis in Hippocampal Model Neurons

How do neurons reconcile the maintenance of calcium homeostasis with perpetual switches in patterns of afferent activity? Here, we assessed state-dependent evolution of calcium homeostasis in a population of hippocampal pyramidal neuron models, through an adaptation of a recent study on stomatogastric ganglion neurons. Calcium homeostasis was set to emerge through cell-autonomous updates to 12 ionic conductances, responding to different types of synaptically driven afferent activity. We first assessed the impact of theta-frequency inputs on the evolution of ionic conductances toward maintenance of calcium homeostasis. Although calcium homeostasis emerged efficaciously across all models in the population, disparate changes in ionic conductances that mediated this emergence resulted in variable plasticity to several intrinsic properties, also manifesting as significant differences in firing responses across models. Assessing the sensitivity of this form of plasticity, we noted that intrinsic neuronal properties and the firing response were sensitive to the target calcium concentration and to the strength and frequency of afferent activity. Next, we studied the evolution of calcium homeostasis when afferent activity was switched, in different temporal sequences, between two behaviorally distinct types of activity: theta-frequency inputs and sharp-wave ripples riding on largely silent periods. We found that the conductance values, intrinsic properties, and firing response of neurons exhibited differential robustness to an intervening switch in the type of afferent activity. These results unveil critical dissociations between different forms of homeostasis, and call for a systematic evaluation of the impact of state-dependent switches in afferent activity on neuronal intrinsic properties during neural coding and homeostasis.

Differential effects of conductances on the phase resetting curve of a bursting neuronal oscillator

The intrinsically oscillating neurons in the crustacean pyloric circuit have membrane conductances that influence their spontaneous activity patterns and responses to synaptic activity. The relationship between the magnitudes of these membrane conductances and the response of the oscillating neurons to synaptic input has not yet been fully or systematically explored. We examined this relationship using the phase resetting curve (PRC), which summarizes the change in the cycle period of a neuronal oscillator as a function of the input's timing within the oscillation. We first utilized a large database of single-compartment model neurons to determine the effect of individual membrane conductances on PRC shape; we found that the effects vary across conductance space, but on average, the hyperpolarization-activated and leak conductances advance the PRC. We next investigated how membrane conductances affect PRCs of the isolated pacemaker kernel in the pyloric circuit of Cancer borealis by: (1) tabulating PRCs while using dynamic clamp to artificially add varying levels of specific conductances, and (2) tabulating PRCs before and after blocking the endogenous hyperpolarization-activated current. We additionally used a previously described four-compartment model to determine how the location of the hyperpolarization-activated conductance influences that current's effect on the PRC. We report that while dynamic-clamp-injected leak current has much smaller effects on the PRC than suggested by the single-compartment model, an increase in the hyperpolarization-activated conductance both advances and reduces the noisiness of the PRC in the pacemaker kernel of the pyloric circuit in both modeling and experimental studies.

Neuromodulation of circuits with variable parameters: single neurons and small circuits reveal principles of state-dependent and robust neuromodulation

Neuromodulation underlies many behavioral states and has been extensively studied in small circuits. This has allowed the systematic exploration of how neuromodulatory substances and the neurons that release them can influence circuit function. The physiological state of a network and its level of activity can have profound effects on how the modulators act, a phenomenon known as state dependence. We provide insights from experiments and computational work that show how state dependence can arise and the consequences it can have for cellular and circuit function. These observations pose a general unsolved question that is relevant to all nervous systems: How is robust modulation achieved in spite of animal-to-animal variability and degenerate, nonlinear mechanisms for the production of neuronal and network activity?

Tonic 5nM DA stabilizes neuronal output by enabling bidirectional activity-dependent regulation of the hyperpolarization activated current via PKA and calcineurin

Volume transmission results in phasic and tonic modulatory signals. The actions of tonic dopamine (DA) at type 1 DA receptors (D1Rs) are largely undefined. Here we show that tonic 5nM DA acts at D1Rs to stabilize neuronal output over minutes by enabling activity-dependent regulation of the hyperpolarization activated current (I _h). In the presence but not absence of 5nM DA, I _h maximal conductance (G _max) was adjusted according to changes in slow wave activity in order to maintain spike timing. Our study on the lateral pyloric neuron (LP), which undergoes rhythmic oscillations in membrane potential with depolarized plateaus, demonstrated that incremental, bi-directional changes in plateau duration produced corresponding alterations in LP I _hG _max when preparations were superfused with saline containing 5nM DA. However, when preparations were superfused with saline alone there was no linear correlation between LP I _hG_max and duty cycle. Thus, tonic nM DA modulated the capacity for activity to modulate LP I _h G _max; this exemplifies metamodulation (modulation of modulation). Pretreatment with the Ca2+-chelator, BAPTA, or the specific PKA inhibitor, PKI, prevented all changes in LP I _h in 5nM DA. Calcineurin inhibitors blocked activity-dependent changes enabled by DA and revealed a PKA-mediated, activity-independent enhancement of LP I _hG _max. These data suggested that tonic 5nM DA produced two simultaneous, PKA-dependent effects: a direct increase in LP I _h G _max and a priming event that permitted calcineurin regulation of LP I _h. The latter produced graded reductions in LP I _hG _max with increasing duty cycles. We also demonstrated that this metamodulation preserved the timing of LP’s first spike when network output was perturbed with bath-applied 4AP. In sum, 5nM DA permits slow wave activity to provide feedback that maintains spike timing, suggesting that one function of low-level, tonic modulation is to stabilize specific features of a dynamic output.

Robust circuit rhythms in small circuits arise from variable circuit components and mechanisms

Small central pattern generating circuits found in invertebrates have significant advantages for the study of the circuit mechanisms that generate brain rhythms. Experimental and computational studies of small oscillatory circuits reveal that similar rhythms can arise from disparate mechanisms. Animal-to-animal variation in the properties of single neurons and synapses may underly robust circuit performance, and can be revealed by perturbations. Neuromodulation can produce altered circuit performance but also ensure reliable circuit function.

Electrical coupling and innexin expression in the stomatogastric ganglion of the crab Cancer borealis

Gap junctions are intercellular channels that allow for the movement of small molecules and ions between the cytoplasm of adjacent cells and form electrical synapses between neurons. In invertebrates, the gap junction proteins are coded for by the innexin family of genes. The stomatogastric ganglion (STG) in the crab Cancer borealis contains a small number of identified and electrically coupled neurons. We identified Innexin 1 (Inx1), Innexin 2 (Inx2), Innexin 3 (Inx3), Innexin 4 (Inx4), Innexin 5 (Inx5), and Innexin 6 (Inx6) members of the C. borealis innexin family. We also identified six members of the innexin family from the lobster Homarus americanus transcriptome. These innexins show significant sequence similarity to other arthropod innexins. Using in situ hybridization and reverse transcriptase-quantitative PCR (RT-qPCR), we determined that all the cells in the crab STG express multiple innexin genes. Electrophysiological recordings of coupling coefficients between identified pairs of pyloric dilator (PD) cells and PD-lateral posterior gastric (LPG) neurons show that the PD-PD electrical synapse is nonrectifying while the PD-LPG synapse is apparently strongly rectifying.

Inter-subject variability in human atrial action potential in sinus rhythm versus chronic atrial fibrillation

AIMS

Human atrial electrophysiology exhibits high inter-subject variability in both sinus rhythm (SR) and chronic atrial fibrillation (cAF) patients. Variability is however rarely investigated in experimental and theoretical electrophysiological studies, thus hampering the understanding of its underlying causes but also its implications in explaining differences in the response to disease and treatment. In our study, we aim at investigating the ability of populations of human atrial cell models to capture the inter-subject variability in action potential (AP) recorded in 363 patients both under SR and cAF conditions.

METHODS AND RESULTS

Human AP recordings in atrial trabeculae (n = 469) from SR and cAF patients were used to calibrate populations of computational SR and cAF atrial AP models. Three populations of over 2000 sampled models were generated, based on three different human atrial AP models. Experimental calibration selected populations of AP models yielding AP with morphology and duration in range with experimental recordings. Populations using the three original models can mimic variability in experimental AP in both SR and cAF, with median conductance values in SR for most ionic currents deviating less than 30% from their original peak values. All cAF populations show similar variations in G(K1), G(Kur) and G(to), consistent with AF-related remodeling as reported in experiments. In all SR and cAF model populations, inter-subject variability in I(K1) and I(NaK) underlies variability in APD90, variability in I(Kur), I(CaL) and I(NaK) modulates variability in APD50 and combined variability in Ito and I(Kur) determines variability in APD20. The large variability in human atrial AP triangulation is mostly determined by I(K1) and either I(NaK) or I(NaCa) depending on the model.

CONCLUSION

Experimentally-calibrated human atrial AP models populations mimic AP variability in SR and cAF patient recordings, and identify potential ionic determinants of inter-subject variability in human atrial AP duration and morphology in SR versus cAF.

Activity-dependent feedback regulates correlated ion channel mRNA levels in single identified motor neurons

Neurons generate cell-specific outputs via interactions of conductances carried by ion channel proteins that are homeostatically regulated to maintain key quantitative relationships among subsets of conductances. Given the challenges of both normal channel protein turnover and short-term plasticity, how is the balance of membrane conductances maintained over long-term timescales to ensure stable electrophysiological phenotype? One possible mechanism is to dynamically regulate production of channel protein via feedback that constrains relationships at the channel mRNA level. Recent modeling work has postulated that such mRNA relationships could emerge as a result of activity-dependent homeostatic tuning rules that ensure an appropriate ratio of mRNA for key ion channels is maintained to preserve robust cellular output. Yet, this has never been demonstrated in biological neurons. In this study, we quantified multiple ion channel mRNAs from single identified motor neurons of the stomatogastric ganglion to determine whether correlations among channel mRNAs are actively maintained, and, if so, by what form of feedback. In these neurons, we identified correlations among mRNAs for voltage-gated calcium and potassium channels. By performing experiments that decoupled activity, synaptic connectivity, and neuromodulatory state, we determined that correlated channel mRNAs are maintained by an activity-dependent process. This is the first study to demonstrate that distinct relationships across channel mRNAs are dynamically maintained in an activity-dependent manner. This feedback from cellular activity to coordinated transcriptome-level interactions represents a novel aspect of regulation of neuronal output with implications for long-term stability of neuron function.

Quantitative expression profiling in mouse spinal cord reveals changing relationships among channel and receptor mRNA levels across postnatal maturation

Neural networks ultimately arrive at functional output via interaction of the excitability of individual neurons and their synaptic interactions. We investigated the relationships between voltage-gated ion channel and neurotransmitter receptor mRNA levels in mouse spinal cord at four different postnatal time points (P5, P11, P17, and adult) and three different adult cord levels (cervical, thoracic, and lumbosacral) using quantitative reverse-transcription polymerase chain reaction (qRT-PCR). Our analysis and data visualization are novel in that we chose a focal group of voltage-gated channel subunits and transmitter receptor subunits, performed absolute quantitation of mRNA copy number for each gene from a sample, and used multiple correlation analyses and correlation matrices to detect patterns in correlated mRNA levels across all genes of interest. These correlation profiles suggest that postnatal maturation of the spinal cord includes changes among channel and receptor subunits that proceed from widespread co-regulation to more refined and distinct functional relationships.

Channeling the Central Dogma

How do neurons and networks achieve their characteristic electrical activity, regulate this activity homeostatically, and yet show population variability in expression? In this issue of Neuron, O'Leary et al. (2014) address some of these thorny questions in this theoretical analysis that starts with the Central Dogma.

Cell types, network homeostasis, and pathological compensation from a biologically plausible ion channel expression model

How do neurons develop, control, and maintain their electrical signaling properties in spite of ongoing protein turnover and perturbations to activity? From generic assumptions about the molecular biology underlying channel expression, we derive a simple model and show how it encodes an "activity set point" in single neurons. The model generates diverse self-regulating cell types and relates correlations in conductance expression observed in vivo to underlying channel expression rates. Synaptic as well as intrinsic conductances can be regulated to make a self-assembling central pattern generator network; thus, network-level homeostasis can emerge from cell-autonomous regulation rules. Finally, we demonstrate that the outcome of homeostatic regulation depends on the complement of ion channels expressed in cells: in some cases, loss of specific ion channels can be compensated; in others, the homeostatic mechanism itself causes pathological loss of function.

Ionic Current Variability and Functional Stability in the Nervous System

Identified neurons in different animals express ionic currents at highly variable levels (population variability). If neuronal identity is associated with stereotypical function, as is the case in genetically identical neurons or in unambiguously identified individual neurons, this variability poses a conundrum: How is activity the same if the components that generate it-ionic current levels-are different? In some cases, ionic current variability across similar neurons generates an output gradient. However, many neurons produce very similar output activity, despite substantial variability in ionic conductances. It appears that, in many such cells, conductance levels of one ionic current vary in proportion to the conductance levels of another current. As a result, in a population of neurons, these conductances appear to be correlated. Here, I review theoretical and experimental work that suggests that neuronal ionic current correlation can reduce the global ionic current variability and can contribute to functional stability.

Many parameter sets in a multicompartment model oscillator are robust to temperature perturbations

Neurons in cold-blooded animals remarkably maintain their function over a wide range of temperatures, even though the rates of many cellular processes increase twofold, threefold, or many-fold for each 10°C increase in temperature. Moreover, the kinetics of ion channels, maximal conductances, and Ca(2+) buffering each have independent temperature sensitivities, suggesting that the balance of biological parameters can be disturbed by even modest temperature changes. In stomatogastric ganglia of the crab Cancer borealis, the duty cycle of the bursting pacemaker kernel is highly robust between 7 and 23°C (Rinberg et al., 2013). We examined how this might be achieved in a detailed conductance-based model in which exponential temperature sensitivities were given by Q10 parameters. We assessed the temperature robustness of this model across 125,000 random sets of Q10 parameters. To examine how robustness might be achieved across a variable population of animals, we repeated this analysis across six sets of maximal conductance parameters that produced similar activity at 11°C. Many permissible combinations of maximal conductance and Q10 parameters were found over broad regions of parameter space and relatively few correlations among Q10s were observed across successful parameter sets. A significant portion of Q10 sets worked for at least 3 of the 6 maximal conductance sets (∼11.1%). Nonetheless, no Q10 set produced robust function across all six maximal conductance sets, suggesting that maximal conductance parameters critically contribute to temperature robustness. Overall, these results provide insight into principles of temperature robustness in neuronal oscillators.

Orchestrated increase of dopamine and PARK mRNAs but not miR-133b in dopamine neurons in Parkinson's disease

Progressive loss of substantia nigra dopamine neurons (SN DA) is a hallmark of aging and of Parkinson's disease (PD). Mutations in PARK genes cause familial PD forms. Increased expression of alpha-synuclein (PARK4) is a disease-triggering event in familial PD and also observed in SN DA neurons in sporadic PD but related transcriptional changes are unknown. With optimized single-cell quantitative real-time polymerase chain reaction analysis, we compared messenger RNA and microRNA levels in SN DA neurons from sporadic PD patients and controls. Non-optimally matched donor ages and RNA integrities are common problems when analyzing human samples. We dissected the influence of distinct ages and RNA integrities of our samples by applying a specifically-optimized, linear-mixed-effects model to quantitative real-time polymerase chain reaction-data. We identified that elevated alpha-synuclein messenger RNA levels in SN DA neurons of human PD brains were positively correlated with corresponding elevated levels of mRNAs for functional compensation of progressive SN DA loss and for enhanced proteasomal (PARK5/UCHL1) and lysosomal (PARK9/ATPase13A2) function, possibly counteracting alpha-synuclein toxicity. In contrast, microRNA miR-133b levels, previously implicated in transcriptional dysregulation in PD, were not altered in SN DA neurons in PD.

Correlated conductance parameters in leech heart motor neurons contribute to motor pattern formation

Neurons can have widely differing intrinsic membrane properties, in particular the density of specific conductances, but how these contribute to characteristic neuronal activity or pattern formation is not well understood. To explore the relationship between conductances, and in particular how they influence the activity of motor neurons in the well characterized leech heartbeat system, we developed a new multi-compartmental Hodgkin-Huxley style leech heart motor neuron model. To do so, we evolved a population of model instances, which differed in the density of specific conductances, capable of achieving specific output activity targets given an associated input pattern. We then examined the sensitivity of measures of output activity to conductances and how the model instances responded to hyperpolarizing current injections. We found that the strengths of many conductances, including those with differing dynamics, had strong partial correlations and that these relationships appeared to be linked by their influence on heart motor neuron activity. Conductances that had positive correlations opposed one another and had the opposite effects on activity metrics when perturbed whereas conductances that had negative correlations could compensate for one another and had similar effects on activity metrics.