Pareto optimality, economy-effectiveness trade-offs and ion channel degeneracy: improving population modelling for single neurons

The emergence of nonlinear evolutionary trade-offs and the maintenance of genetic polymorphisms

Evolutionary models of quantitative traits often assume trade-offs between beneficial and detrimental traits, requiring modellers to specify a function linking trait values. The choice of trade-off function can be consequential; functions that assume diminishing returns (accelerating costs) typically lead to single equilibrium genotypes, while decelerating costs often lead to genetic polymorphisms. Despite their importance, our current theory has little to say on which trade-off functions are the most biologically plausible. To address this gap, we explored how the genetic determination of quantitative traits can lead to different trade-off functions, using resistance to infectious diseases as an example trait. We developed a model where alleles at separate loci pleiotropically increase resistance while decreasing fecundity. We then used this model to generate genotype landscapes and investigate how epistasis effects the trade-off function. Regardless of the strength of epistasis, our model consistently led to accelerating costs. We then incorporated our genotype model into an eco-evolutionary model of disease resistance. Unlike other models with accelerating costs, our approach often led to genetic polymorphisms. Our results suggest that accelerating costs are a strong null model for evolutionary trade-offs and that the eco-evolutionary conditions required for polymorphism may be more nuanced than previously thought.

Ion-channel degeneracy and heterogeneities in the emergence of signature physiological characteristics of dentate gyrus granule cells

Complex systems are neither fully determined nor completely random. Biological complex systems, including single neurons, manifest intermediate regimes of randomness that recruit integration of specific combinations of functionally specialized subsystems. Such emergence of biological function provides the substrate for the expression of degeneracy, the ability of disparate combinations of subsystems to yield similar function. Here, we present evidence for the expression of degeneracy in morphologically realistic models of dentate gyrus granule cells (GCs) through functional integration of disparate ion-channel combinations. We performed a 45-parameter randomized search spanning 16 active and passive ion channels, each biophysically constrained by their gating kinetics and localization profiles, to search for valid GC models. Valid models were those that satisfied 17 sub- and suprathreshold cellular-scale electrophysiological measurements from rat GCs. A vast majority (>99%) of the 15,000 random models were not electrophysiologically valid, demonstrating that arbitrarily random ion-channel combinations would not yield GC functions. The 141 valid models (0.94% of 15,000) manifested heterogeneities in and cross-dependencies across local and propagating electrophysiological measurements, which matched with their respective biological counterparts. Importantly, these valid models were widespread throughout the parametric space and manifested weak cross-dependencies across different parameters. These observations together showed that GC physiology could neither be obtained by entirely random ion-channel combinations nor is there an entirely determined single parametric combination that satisfied all constraints. The complexity, the heterogeneities in measurement and parametric spaces, and degeneracy associated with GC physiology should be rigorously accounted for while assessing GCs and their robustness under physiological and pathological conditions.NEW & NOTEWORTHY A recent study from our laboratory had demonstrated pronounced heterogeneities in a set of 17 electrophysiological measurements obtained from a large population of rat hippocampal granule cells. Here, we demonstrate the manifestation of ion-channel degeneracy in a heterogeneous population of morphologically realistic conductance-based granule cell models that were validated against these measurements and their cross-dependencies. Our analyses show that single neurons are complex entities whose functions emerge through intricate interactions among several functionally specialized subsystems.

Tetrapod sperm length evolution in relation to body mass is shaped by multiple trade-offs

Sperm length is highly variable across species and many questions about its variation remain open. Although variation in body mass may affect sperm length evolution through its influence on multiple factors, the extent to which sperm length variation is linked to body mass remains elusive. Here, we use the Pareto multi-task evolution framework to investigate the relationship between sperm length and body mass across tetrapods. We find that tetrapods occupy a triangular Pareto front, indicating that trade-offs shape the evolution of sperm length in relation to body mass. By exploring the factors predicted to influence sperm length evolution, we find that sperm length evolution is mainly driven by sperm competition and clutch size, rather than by genome size. Moreover, the triangular Pareto front is maintained within endotherms, internal fertilizers, mammals and birds, suggesting similar evolutionary trade-offs within tetrapods. Finally, we demonstrate that the Pareto front is robust to phylogenetic dependencies and finite sampling bias. Our findings provide insights into the evolutionary mechanisms driving interspecific sperm length variation and highlight the importance of considering multiple trade-offs in optimizing reproductive traits.

Network motifs in cellular neurophysiology

Concepts from network science and graph theory, including the framework of network motifs, have been frequently applied in studying neuronal networks and other biological complex systems. Network-based approaches can also be used to study the functions of individual neurons, where cellular elements such as ion channels and membrane voltage are conceptualized as nodes within a network, and their interactions are denoted by edges. Network motifs in this context provide functional building blocks that help to illuminate the principles of cellular neurophysiology. In this review we build a case that network motifs operating within neurons provide tools for defining the functional architecture of single-neuron physiology and neuronal adaptations. We highlight the presence of such computational motifs in the cellular mechanisms underlying action potential generation, neuronal oscillations, dendritic integration, and neuronal plasticity. Future work applying the network motifs perspective may help to decipher the functional complexities of neurons and their adaptation during health and disease.

A Multimodal Fitting Approach to Construct Single-Neuron Models With Patch Clamp and High-Density Microelectrode Arrays

In computational neuroscience, multicompartment models are among the most biophysically realistic representations of single neurons. Constructing such models usually involves the use of the patch-clamp technique to record somatic voltage signals under different experimental conditions. The experimental data are then used to fit the many parameters of the model. While patching of the soma is currently the gold-standard approach to build multicompartment models, several studies have also evidenced a richness of dynamics in dendritic and axonal sections. Recording from the soma alone makes it hard to observe and correctly parameterize the activity of nonsomatic compartments. In order to provide a richer set of data as input to multicompartment models, we here investigate the combination of somatic patch-clamp recordings with recordings of high-density microelectrode arrays (HD-MEAs). HD-MEAs enable the observation of extracellular potentials and neural activity of neuronal compartments at subcellular resolution. In this work, we introduce a novel framework to combine patch-clamp and HD-MEA data to construct multicompartment models. We first validate our method on a ground-truth model with known parameters and show that the use of features extracted from extracellular signals, in addition to intracellular ones, yields models enabling better fits than using intracellular features alone. We also demonstrate our procedure using experimental data by constructing cell models from in vitro cell cultures. The proposed multimodal fitting procedure has the potential to augment the modeling efforts of the computational neuroscience community and provide the field with neuronal models that are more realistic and can be better validated.

The Emergence of Non-Linear Evolutionary Trade-offs and the Maintenance of Genetic Polymorphisms

Evolutionary models of quantitative traits often assume trade-offs between beneficial and detrimental traits, requiring modelers to specify a function linking costs to benefits. The choice of trade-off function is often consequential; functions that assume diminishing returns (accelerating costs) typically lead to single equilibrium genotypes, while decelerating costs often lead to evolutionary branching. Despite their importance, we still lack a strong theoretical foundation to base the choice of trade-off function. To address this gap, we explore how trade-off functions can emerge from the genetic architecture of a quantitative trait. We developed a multi-locus model of disease resistance, assuming each locus had random antagonistic pleiotropic effects on resistance and fecundity. We used this model to generate genotype landscapes and explored how additive versus epistatic genetic architectures influenced the shape of the trade-off function. Regardless of epistasis, our model consistently led to accelerating costs. We then used our genotype landscapes to build an evolutionary model of disease resistance. Unlike other models with accelerating costs, our approach often led to genetic polymorphisms at equilibrium. Our results suggest that accelerating costs are a strong null model for evolutionary trade-offs and that the eco-evolutionary conditions required for polymorphism may be more nuanced than previously believed.

Robust and consistent measures of pattern separation based on information theory and demonstrated in the dentate gyrus

Pattern separation is a valuable computational function performed by neuronal circuits, such as the dentate gyrus, where dissimilarity between inputs is increased, reducing noise and increasing the storage capacity of downstream networks. Pattern separation is studied from both in vivo experimental and computational perspectives and, a number of different measures (such as orthogonalisation, decorrelation, or spike train distance) have been applied to quantify the process of pattern separation. However, these are known to give conclusions that can differ qualitatively depending on the choice of measure and the parameters used to calculate it. We here demonstrate that arbitrarily increasing sparsity, a noticeable feature of dentate granule cell firing and one that is believed to be key to pattern separation, typically leads to improved classical measures for pattern separation even, inappropriately, up to the point where almost all information about the inputs is lost. Standard measures therefore both cannot differentiate between pattern separation and pattern destruction, and give results that may depend on arbitrary parameter choices. We propose that techniques from information theory, in particular mutual information, transfer entropy, and redundancy, should be applied to penalise the potential for lost information (often due to increased sparsity) that is neglected by existing measures. We compare five commonly-used measures of pattern separation with three novel techniques based on information theory, showing that the latter can be applied in a principled way and provide a robust and reliable measure for comparing the pattern separation performance of different neurons and networks. We demonstrate our new measures on detailed compartmental models of individual dentate granule cells and a dentate microcircuit, and show how structural changes associated with epilepsy affect pattern separation performance. We also demonstrate how our measures of pattern separation can predict pattern completion accuracy. Overall, our measures solve a widely acknowledged problem in assessing the pattern separation of neural circuits such as the dentate gyrus, as well as the cerebellum and mushroom body. Finally we provide a publicly available toolbox allowing for easy analysis of pattern separation in spike train ensembles.

Controlling morpho-electrophysiological variability of neurons with detailed biophysical models

Variability, which is known to be a universal feature among biological units such as neuronal cells, holds significant importance, as, for example, it enables a robust encoding of a high volume of information in neuronal circuits and prevents hypersynchronizations. While most computational studies on electrophysiological variability in neuronal circuits were done with single-compartment neuron models, we instead focus on the variability of detailed biophysical models of neuron multi-compartmental morphologies. We leverage a Markov chain Monte Carlo method to generate populations of electrical models reproducing the variability of experimental recordings while being compatible with a set of morphologies to faithfully represent specifi morpho-electrical type. We demonstrate our approach on layer 5 pyramidal cells and study the morpho-electrical variability and in particular, find that morphological variability alone is insufficient to reproduce electrical variability. Overall, this approach provides a strong statistical basis to create detailed models of neurons with controlled variability.

Heterogeneous off-target impact of ion-channel deletion on intrinsic properties of hippocampal model neurons that self-regulate calcium

How do neurons that implement cell-autonomous self-regulation of calcium react to knockout of individual ion-channel conductances? To address this question, we used a heterogeneous population of 78 conductance-based models of hippocampal pyramidal neurons that maintained cell-autonomous calcium homeostasis while receiving theta-frequency inputs. At calcium steady-state, we individually deleted each of the 11 active ion-channel conductances from each model. We measured the acute impact of deleting each conductance (one at a time) by comparing intrinsic electrophysiological properties before and immediately after channel deletion. The acute impact of deleting individual conductances on physiological properties (including calcium homeostasis) was heterogeneous, depending on the property, the specific model, and the deleted channel. The underlying many-to-many mapping between ion channels and properties pointed to ion-channel degeneracy. Next, we allowed the other conductances (barring the deleted conductance) to evolve towards achieving calcium homeostasis during theta-frequency activity. When calcium homeostasis was perturbed by ion-channel deletion, post-knockout plasticity in other conductances ensured resilience of calcium homeostasis to ion-channel deletion. These results demonstrate degeneracy in calcium homeostasis, as calcium homeostasis in knockout models was implemented in the absence of a channel that was earlier involved in the homeostatic process. Importantly, in reacquiring homeostasis, ion-channel conductances and physiological properties underwent heterogenous plasticity (dependent on the model, the property, and the deleted channel), even introducing changes in properties that were not directly connected to the deleted channel. Together, post-knockout plasticity geared towards maintaining homeostasis introduced heterogenous off-target effects on several channels and properties, suggesting that extreme caution be exercised in interpreting experimental outcomes involving channel knockouts.

Biological complexity facilitates tuning of the neuronal parameter space

The electrical and computational properties of neurons in our brains are determined by a rich repertoire of membrane-spanning ion channels and elaborate dendritic trees. However, the precise reason for this inherent complexity remains unknown, given that simpler models with fewer ion channels are also able to functionally reproduce the behaviour of some neurons. Here, we stochastically varied the ion channel densities of a biophysically detailed dentate gyrus granule cell model to produce a large population of putative granule cells, comparing those with all 15 original ion channels to their reduced but functional counterparts containing only 5 ion channels. Strikingly, valid parameter combinations in the full models were dramatically more frequent at -6% vs. -1% in the simpler model. The full models were also more stable in the face of perturbations to channel expression levels. Scaling up the numbers of ion channels artificially in the reduced models recovered these advantages confirming the key contribution of the actual number of ion channel types. We conclude that the diversity of ion channels gives a neuron greater flexibility and robustness to achieve a target excitability.

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.

Degeneracy in epilepsy: multiple routes to hyperexcitable brain circuits and their repair

Due to its complex and multifaceted nature, developing effective treatments for epilepsy is still a major challenge. To deal with this complexity we introduce the concept of degeneracy to the field of epilepsy research: the ability of disparate elements to cause an analogous function or malfunction. Here, we review examples of epilepsy-related degeneracy at multiple levels of brain organisation, ranging from the cellular to the network and systems level. Based on these insights, we outline new multiscale and population modelling approaches to disentangle the complex web of interactions underlying epilepsy and to design personalised multitarget therapies.