Loss of neuronal heterogeneity in epileptogenic human tissue impairs network resilience to sudden changes in synchrony

Rigor and reproducibility in human brain organoid research: Where we are and where we need to go

Human brain organoid models have emerged as a promising tool for studying human brain development and function. These models preserve human genetics and recapitulate some aspects of human brain development, while facilitating manipulation in an in vitro setting. Despite their potential to transform biology and medicine, concerns persist about their fidelity. To fully harness their potential, it is imperative to establish reliable analytic methods, ensuring rigor and reproducibility. Here, we review current analytical platforms used to characterize human forebrain cortical organoids, highlight challenges, and propose recommendations for future studies to achieve greater precision and uniformity across laboratories.

The Neurovascular Unit Dysfunction in the Molecular Mechanisms of Epileptogenesis and Targeted Therapy

Epilepsy is a multifaceted neurological syndrome characterized by recurrent, spontaneous, and synchronous seizures. The pathogenesis of epilepsy, known as epileptogenesis, involves intricate changes in neurons, neuroglia, and endothelium, leading to structural and functional disorders within neurovascular units and culminating in the development of spontaneous epilepsy. Although current research on epilepsy treatments primarily centers around anti-seizure drugs, it is imperative to seek effective interventions capable of disrupting epileptogenesis. To this end, a comprehensive exploration of the changes and the molecular mechanisms underlying epileptogenesis holds the promise of identifying vital biomarkers for accurate diagnosis and potential therapeutic targets. Emphasizing early diagnosis and timely intervention is paramount, as it stands to significantly improve patient prognosis and alleviate the socioeconomic burden. In this review, we highlight the changes and molecular mechanisms of the neurovascular unit in epileptogenesis and provide a theoretical basis for identifying biomarkers and drug targets.

Neural heterogeneity controls computations in spiking neural networks

The brain is composed of complex networks of interacting neurons that express considerable heterogeneity in their physiology and spiking characteristics. How does this neural heterogeneity influence macroscopic neural dynamics, and how might it contribute to neural computation? In this work, we use a mean-field model to investigate computation in heterogeneous neural networks, by studying how the heterogeneity of cell spiking thresholds affects three key computational functions of a neural population: the gating, encoding, and decoding of neural signals. Our results suggest that heterogeneity serves different computational functions in different cell types. In inhibitory interneurons, varying the degree of spike threshold heterogeneity allows them to gate the propagation of neural signals in a reciprocally coupled excitatory population. Whereas homogeneous interneurons impose synchronized dynamics that narrow the dynamic repertoire of the excitatory neurons, heterogeneous interneurons act as an inhibitory offset while preserving excitatory neuron function. Spike threshold heterogeneity also controls the entrainment properties of neural networks to periodic input, thus affecting the temporal gating of synaptic inputs. Among excitatory neurons, heterogeneity increases the dimensionality of neural dynamics, improving the network's capacity to perform decoding tasks. Conversely, homogeneous networks suffer in their capacity for function generation, but excel at encoding signals via multistable dynamic regimes. Drawing from these findings, we propose intra-cell-type heterogeneity as a mechanism for sculpting the computational properties of local circuits of excitatory and inhibitory spiking neurons, permitting the same canonical microcircuit to be tuned for diverse computational tasks.

Diversity-induced trivialization and resilience of neural dynamics

Heterogeneity is omnipresent across all living systems. Diversity enriches the dynamical repertoire of these systems but remains challenging to reconcile with their manifest robustness and dynamical persistence over time, a fundamental feature called resilience. To better understand the mechanism underlying resilience in neural circuits, we considered a nonlinear network model, extracting the relationship between excitability heterogeneity and resilience. To measure resilience, we quantified the number of stationary states of this network, and how they are affected by various control parameters. We analyzed both analytically and numerically gradient and non-gradient systems modeled as non-linear sparse neural networks evolving over long time scales. Our analysis shows that neuronal heterogeneity quenches the number of stationary states while decreasing the susceptibility to bifurcations: a phenomenon known as trivialization. Heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in network size and connection probability by quenching the system's dynamic volatility.

Induced synchronization by endogenous noise modulation in finite-size random neural networks: A stochastic mean-field study

Event-related synchronization and desynchronization (ERS/ERD) are well-known features found experimentally in brain signals during cognitive tasks. Their understanding promises to have much better insights into neural information processes in cognition. Under the hypothesis that neural information affects the endogenous neural noise level in populations, we propose to employ a stochastic mean-field model to explain ERS/ERD in the γ-frequency range. The work extends previous mean-field studies by deriving novel effects from finite network size. Moreover, numerical simulations of ERS/ERD and their analytical explanation by the mean-field model suggest several endogenous noise modulation schemes, which may modulate the system's synchronization.

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.

Nonmonotonic enhancement of diversity-induced resonance in systems of mobile oscillators

Diversity is omnipresent in natural and synthetic extended systems, the phenomenon of diversity-induced resonance (DIR), wherein a moderate degree of the diversity can provoke an optimal collective response, provides researchers a brand-new strategy to amplify and utilize the weak signal. As yet the relevant advances focus mostly on the ideal situations where the interactions among elements are uncorrelated with the physical proximity of agents. Such a consideration overlooks interactions mediated by the motion of agents in space. Here, we investigate the signal response of an ensemble of spatial mobile heterogeneous bistable oscillators with two canonical interacting modes: dynamic and preset. The oscillators are considered as mass points and perform random walks in a two-dimensional square plane. Under the dynamic scheme, the oscillators can only interact with other oscillators within a fixed vision radius. For the preset circumstance, the interaction among oscillators occurs only when all of them are in a predefined region at the same moment. We find that the DIR can be obtained in both situations. Additionally, the strength of resonance nonmonotonically rises with respect to the increase of moving speed, and the optimal resonance is acquired by an intermediate magnitude of speed. Finally, we propose reduced equations to guarantee the occurrence of such mobility-optimized DIR on the basis of the fast switching approximation theory and also examine the robustness of such phenomenon through the excitable FitzHugh-Nagumo model and a different spatial motion mechanism. Our results reveal for the first time that the DIR can be optimized by the spatial mobility and thus has promising potential application in the communication of mobile agents.

Distinctive biophysical features of human cell-types: insights from studies of neurosurgically resected brain tissue

Electrophysiological characterization of live human tissue from epilepsy patients has been performed for many decades. Although initially these studies sought to understand the biophysical and synaptic changes associated with human epilepsy, recently, it has become the mainstay for exploring the distinctive biophysical and synaptic features of human cell-types. Both epochs of these human cellular electrophysiological explorations have faced criticism. Early studies revealed that cortical pyramidal neurons obtained from individuals with epilepsy appeared to function "normally" in comparison to neurons from non-epilepsy controls or neurons from other species and thus there was little to gain from the study of human neurons from epilepsy patients. On the other hand, contemporary studies are often questioned for the "normalcy" of the recorded neurons since they are derived from epilepsy patients. In this review, we discuss our current understanding of the distinct biophysical features of human cortical neurons and glia obtained from tissue removed from patients with epilepsy and tumors. We then explore the concept of within cell-type diversity and its loss (i.e., "neural homogenization"). We introduce neural homogenization to help reconcile the epileptogenicity of seemingly "normal" human cortical cells and circuits. We propose that there should be continued efforts to study cortical tissue from epilepsy patients in the quest to understand what makes human cell-types "human".

KCNQ2/3 Gain-of-Function Variants and Cell Excitability: Differential Effects in CA1 versus L2/3 Pyramidal Neurons

Gain-of-function (GOF) pathogenic variants in the potassium channels KCNQ2 and KCNQ3 lead to hyperexcitability disorders such as epilepsy and autism spectrum disorders. However, the underlying cellular mechanisms of how these variants impair forebrain function are unclear. Here, we show that the R201C variant in KCNQ2 has opposite effects on the excitability of two types of mouse pyramidal neurons of either sex, causing hyperexcitability in layer 2/3 (L2/3) pyramidal neurons and hypoexcitability in CA1 pyramidal neurons. Similarly, the homologous R231C variant in KCNQ3 leads to hyperexcitability in L2/3 pyramidal neurons and hypoexcitability in CA1 pyramidal neurons. However, the effects of KCNQ3 gain-of-function on excitability are specific to superficial CA1 pyramidal neurons. These findings reveal a new level of complexity in the function of KCNQ2 and KCNQ3 channels in the forebrain and provide a framework for understanding the effects of gain-of-function variants and potassium channels in the brain.SIGNIFICANCE STATEMENT KCNQ2/3 gain-of-function (GOF) variants lead to severe forms of neurodevelopmental disorders, but the mechanisms by which these channels affect neuronal activity are poorly understood. In this study, using a series of transgenic mice we demonstrate that the same KCNQ2/3 GOF variants can lead to either hyperexcitability or hypoexcitability in different types of pyramidal neurons [CA1 vs layer (L)2/3]. Additionally, we show that expression of the recurrent KCNQ2 GOF variant R201C in forebrain pyramidal neurons could lead to seizures and SUDEP. Our data suggest that the effects of KCNQ2/3 GOF variants depend on specific cell types and brain regions, possibly accounting for the diverse range of phenotypes observed in individuals with KCNQ2/3 GOF variants.

Hyperpolarization-Activated Cation Channels Shape the Spiking Frequency Preference of Human Cortical Layer 5 Pyramidal Neurons

Discerning the contribution of specific ionic currents to complex neuronal dynamics is a difficult, but important, task. This challenge is exacerbated in the human setting, although the widely characterized uniqueness of the human brain compared with preclinical models necessitates the direct study of human neurons. Neuronal spiking frequency preference is of particular interest given its role in rhythm generation and signal transmission in cortical circuits. Here, we combine the frequency-dependent gain (FDG), a measure of spiking frequency preference, and novel in silico analyses to dissect the contributions of individual ionic currents to the suprathreshold features of human layer 5 (L5) neurons captured by the FDG. We confirm that a contemporary model of such a neuron, primarily constrained to capture subthreshold activity driven by the hyperpolarization-activated cyclic nucleotide gated (h-) current, replicates key features of the in vitro FDG both with and without h-current activity. With the model confirmed as a viable approximation of the biophysical features of interest, we applied new analysis techniques to quantify the activity of each modeled ionic current in the moments before spiking, revealing unique dynamics of the h-current. These findings motivated patch-clamp recordings in analogous rodent neurons to characterize their FDG, which confirmed that a biophysically detailed model of these neurons captures key interspecies differences in the FDG. These differences are correlated with distinct contributions of the h-current to neuronal activity. Together, this interdisciplinary and multispecies study provides new insights directly relating the dynamics of the h-current to suprathreshold spiking frequency preference in human L5 neurons.

The impact of antidepressants on human neurodevelopment: Brain organoids as experimental tools

The use of antidepressants during pregnancy benefits the mother's well-being, but the effects of such substances on neurodevelopment remain poorly understood. Moreover, the consequences of early exposure to antidepressants may not be immediately apparent at birth. In utero exposure to selective serotonin reuptake inhibitors (SSRIs) has been related to developmental abnormalities, including a reduced white matter volume. Several reports have observed an increased incidence of autism spectrum disorder (ASD) and attention deficit hyperactivity disorder (ADHD) after prenatal exposure to SSRIs such as sertraline, the most widely prescribed SSRI. The advent of human-induced pluripotent stem cell (hiPSC) methods and assays now offers appropriate tools to test the consequences of such compounds for neurodevelopment in vitro. In particular, hiPSCs can be used to generate cerebral organoids - self-organized structures that recapitulate the morphology and complex physiology of the developing human brain, overcoming the limitations found in 2D cell culture and experimental animal models for testing drug efficacy and side effects. For example, single-cell RNA sequencing (scRNA-seq) and electrophysiological measurements on organoids can be used to evaluate the impact of antidepressants on the transcriptome and neuronal activity signatures in developing neurons. While the analysis of large-scale transcriptomic data depends on dimensionality reduction methods, electrophysiological recordings rely on temporal data series to discriminate statistical characteristics of neuronal activity, allowing for the rigorous analysis of the effects of antidepressants and other molecules that affect the developing nervous system, especially when applied in combination with relevant human cellular models such as brain organoids.

Intrinsic neural diversity quenches the dynamic volatility of neural networks

Heterogeneity is the norm in biology. The brain is no different: Neuronal cell types are myriad, reflected through their cellular morphology, type, excitability, connectivity motifs, and ion channel distributions. While this biophysical diversity enriches neural systems' dynamical repertoire, it remains challenging to reconcile with the robustness and persistence of brain function over time (resilience). To better understand the relationship between excitability heterogeneity (variability in excitability within a population of neurons) and resilience, we analyzed both analytically and numerically a nonlinear sparse neural network with balanced excitatory and inhibitory connections evolving over long time scales. Homogeneous networks demonstrated increases in excitability, and strong firing rate correlations-signs of instability-in response to a slowly varying modulatory fluctuation. Excitability heterogeneity tuned network stability in a context-dependent way by restraining responses to modulatory challenges and limiting firing rate correlations, while enriching dynamics during states of low modulatory drive. Excitability heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in population size, connection probability, strength and variability of synaptic weights, by quenching the volatility (i.e., its susceptibility to critical transitions) of its dynamics. Together, these results highlight the fundamental role played by cell-to-cell heterogeneity in the robustness of brain function in the face of change.

Heterogeneity of voltage gated sodium current density between neurons decorrelates spiking and suppresses network synchronization in Scn1b null mouse models

Voltage gated sodium channels (VGSCs) are required for action potential initiation and propagation in mammalian neurons. As with other ion channel families, VGSC density varies between neurons. Importantly, sodium current (INa) density variability is reduced in pyramidal neurons of Scn1b null mice. Scn1b encodes the VGSC β1/ β1B subunits, which regulate channel expression, trafficking, and voltage dependent properties. Here, we investigate how variable INa density in cortical layer 6 and subicular pyramidal neurons affects spike patterning and network synchronization. Constitutive or inducible Scn1b deletion enhances spike timing correlations between pyramidal neurons in response to fluctuating stimuli and impairs spike-triggered average current pattern diversity while preserving spike reliability. Inhibiting INa with a low concentration of tetrodotoxin similarly alters patterning without impairing reliability, with modest effects on firing rate. Computational modeling shows that broad INa density ranges confer a similarly broad spectrum of spike patterning in response to fluctuating synaptic conductances. Network coupling of neurons with high INa density variability displaces the coupling requirements for synchronization and broadens the dynamic range of activity when varying synaptic strength and network topology. Our results show that INa heterogeneity between neurons potently regulates spike pattern diversity and network synchronization, expanding VGSC roles in the nervous system.

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.

From phenomenological to biophysical models of seizures

Epilepsy is a complex disease that requires various approaches for its study. In this short review, we discuss the contribution of theoretical and computational models. The review presents theoretical frameworks that underlie the understanding of certain seizure properties and their classification based on their dynamical properties at the onset and offset of seizures. Dynamical system tools are valuable resources in the study of seizures. By analyzing the complex, dynamic behavior of seizures, these tools can provide insights into seizure mechanisms and offer a framework for their classification. Additionally, computational models have high potential for clinical applications, as they can be used to develop more accurate diagnostic and personalized medicine tools. We discuss various modeling approaches that span different scales and levels, while also questioning the neurocentric view, and emphasize the importance of considering glial cells. Finally, we explore the epistemic value provided by this type of approach.

Mapping thalamic innervation to individual L2/3 pyramidal neurons and modeling their 'readout' of visual input

The thalamus is the main gateway for sensory information from the periphery to the mammalian cerebral cortex. A major conundrum has been the discrepancy between the thalamus's central role as the primary feedforward projection system into the neocortex and the sparseness of thalamocortical synapses. Here we use new methods, combining genetic tools and scalable tissue expansion microscopy for whole-cell synaptic mapping, revealing the number, density and size of thalamic versus cortical excitatory synapses onto individual layer 2/3 (L2/3) pyramidal cells (PCs) of the mouse primary visual cortex. We find that thalamic inputs are not only sparse, but remarkably heterogeneous in number and density across individual dendrites and neurons. Most surprising, despite their sparseness, thalamic synapses onto L2/3 PCs are smaller than their cortical counterparts. Incorporating these findings into fine-scale, anatomically faithful biophysical models of L2/3 PCs reveals how individual neurons with sparse and weak thalamocortical synapses, embedded in small heterogeneous neuronal ensembles, may reliably 'read out' visually driven thalamic input.

Macroscopic dynamics of neural networks with heterogeneous spiking thresholds

Mean-field theory links the physiological properties of individual neurons to the emergent dynamics of neural population activity. These models provide an essential tool for studying brain function at different scales; however, for their application to neural populations on large scale, they need to account for differences between distinct neuron types. The Izhikevich single neuron model can account for a broad range of different neuron types and spiking patterns, thus rendering it an optimal candidate for a mean-field theoretic treatment of brain dynamics in heterogeneous networks. Here we derive the mean-field equations for networks of all-to-all coupled Izhikevich neurons with heterogeneous spiking thresholds. Using methods from bifurcation theory, we examine the conditions under which the mean-field theory accurately predicts the dynamics of the Izhikevich neuron network. To this end, we focus on three important features of the Izhikevich model that are subject here to simplifying assumptions: (i) spike-frequency adaptation, (ii) the spike reset conditions, and (iii) the distribution of single-cell spike thresholds across neurons. Our results indicate that, while the mean-field model is not an exact model of the Izhikevich network dynamics, it faithfully captures its different dynamic regimes and phase transitions. We thus present a mean-field model that can represent different neuron types and spiking dynamics. The model comprises biophysical state variables and parameters, incorporates realistic spike resetting conditions, and accounts for heterogeneity in neural spiking thresholds. These features allow for a broad applicability of the model as well as for a direct comparison to experimental data.

Cell to network computational model of the epileptic human hippocampus suggests specific roles of network and channel dysfunctions in the ictal and interictal oscillations

The mechanisms underlying the generation of hippocampal epileptic seizures and interictal events and their interactions with the sleep-wake cycle are not yet fully understood. Indeed, medial temporal lobe epilepsy is associated with hippocampal abnormalities both at the neuronal (channelopathies, impaired potassium and chloride dynamics) and network level (neuronal and axonal loss, mossy fiber sprouting), with more frequent seizures during wakefulness compared with slow-wave sleep. In this article, starting from our previous computational modeling work of the hippocampal formation based on realistic topology and synaptic connectivity, we study the role of micro- and mesoscale pathological conditions of the epileptic hippocampus in the generation and maintenance of seizure-like theta and interictal oscillations. We show, through the simulations of hippocampal activity during slow-wave sleep and wakefulness that: (i) both mossy fiber sprouting and sclerosis account for seizure-like theta activity, (ii) but they have antagonist effects (seizure-like activity occurrence increases with sprouting but decreases with sclerosis), (iii) though impaired potassium and chloride dynamics have little influence on the generation of seizure-like activity, they do play a role on the generation of interictal patterns, and (iv) seizure-like activity and fast ripples are more likely to occur during wakefulness and interictal spikes during sleep.

Molecular Mechanisms of Epilepsy: The Role of the Chloride Transporter KCC2

Epilepsy is a neurological disease characterized by abnormal or synchronous brain activity causing seizures, which may produce convulsions, minor physical signs, or a combination of symptoms. These disorders affect approximately 65 million people worldwide, from all ages and genders. Seizures apart, epileptic patients present a high risk to develop neuropsychological comorbidities such as cognitive deficits, emotional disturbance, and psychiatric disorders, which severely impair quality of life. Currently, the treatment for epilepsy includes the administration of drugs or surgery, but about 30% of the patients treated with antiepileptic drugs develop time-dependent pharmacoresistence. Therefore, further investigation about epilepsy and its causes is needed to find new pharmacological targets and innovative therapeutic strategies. Pharmacoresistance is associated to changes in neuronal plasticity and alterations of GABA_A receptor-mediated neurotransmission. The downregulation of GABA inhibitory activity may arise from a positive shift in GABA_A receptor reversal potential, due to an alteration in chloride homeostasis. In this paper, we review the contribution of K⁺-Cl^--cotransporter (KCC2) to the alterations in the Cl^- gradient observed in epileptic condition, and how these alterations are coupled to the increase in the excitability.