The microcircuits of striatum in silico

Modeling the impact of neuromorphological alterations in Down syndrome on fast neural oscillations

Cognitive disorders, including Down syndrome (DS), present significant morphological alterations in neuron architectural complexity. However, the relationship between neuromorphological alterations and impaired brain function is not fully understood. To address this gap, we propose a novel computational model that accounts for the observed cell deformations in DS. The model consists of a cross-sectional layer of the mouse motor cortex, composed of 3000 neurons. The network connectivity is obtained by accounting explicitly for two single-neuron morphological parameters: the mean dendritic tree radius and the spine density in excitatory pyramidal cells. We obtained these values by fitting reconstructed neuron data corresponding to three mouse models: wild-type (WT), transgenic (TgDyrk1A), and trisomic (Ts65Dn). Our findings reveal a dynamic interplay between pyramidal and fast-spiking interneurons leading to the emergence of gamma activity (∼40 Hz). In the DS models this gamma activity is diminished, corroborating experimental observations and validating our computational methodology. We further explore the impact of disrupted excitation-inhibition balance by mimicking the reduction recurrent inhibition present in DS. In this case, gamma power exhibits variable responses as a function of the external input to the network. Finally, we perform a numerical exploration of the morphological parameter space, unveiling the direct influence of each structural parameter on gamma frequency and power. Our research demonstrates a clear link between changes in morphology and the disruption of gamma oscillations in DS. This work underscores the potential of computational modeling to elucidate the relationship between neuron architecture and brain function, and ultimately improve our understanding of cognitive disorders.

Serotonergic and dopaminergic neurons in the dorsal raphe are differentially altered in a mouse model for parkinsonism

Parkinson's disease (PD) is characterized by motor impairments caused by degeneration of dopamine neurons in the substantia nigra pars compacta. In addition to these symptoms, PD patients often suffer from non-motor comorbidities including sleep and psychiatric disturbances, which are thought to depend on concomitant alterations of serotonergic and noradrenergic transmission. A primary locus of serotonergic neurons is the dorsal raphe nucleus (DRN), providing brain-wide serotonergic input. Here, we identified electrophysiological and morphological parameters to classify serotonergic and dopaminergic neurons in the murine DRN under control conditions and in a PD model, following striatal injection of the catecholamine toxin, 6-hydroxydopamine (6-OHDA). Electrical and morphological properties of both neuronal populations were altered by 6-OHDA. In serotonergic neurons, most changes were reversed when 6-OHDA was injected in combination with desipramine, a noradrenaline (NA) reuptake inhibitor, protecting the noradrenergic terminals. Our results show that the depletion of both NA and dopamine in the 6-OHDA mouse model causes changes in the DRN neural circuitry.

Pathological cell assembly dynamics in a striatal MSN network model

Under normal conditions the principal cells of the striatum, medium spiny neurons (MSNs), show structured cell assembly activity patterns which alternate sequentially over exceedingly long timescales of many minutes. It is important to understand this activity since it is characteristically disrupted in multiple pathologies, such as Parkinson's disease and dyskinesia, and thought to be caused by alterations in the MSN to MSN lateral inhibitory connections and in the strength and distribution of cortical excitation to MSNs. To understand how these long timescales arise we extended a previous network model of MSN cells to include synapses with short-term plasticity, with parameters taken from a recent detailed striatal connectome study. We first confirmed the presence of sequentially switching cell clusters using the non-linear dimensionality reduction technique, Uniform Manifold Approximation and Projection (UMAP). We found that the network could generate non-stationary activity patterns varying extremely slowly on the order of minutes under biologically realistic conditions. Next we used Simulation Based Inference (SBI) to train a deep net to map features of the MSN network generated cell assembly activity to MSN network parameters. We used the trained SBI model to estimate MSN network parameters from ex-vivo brain slice calcium imaging data. We found that best fit network parameters were very close to their physiologically observed values. On the other hand network parameters estimated from Parkinsonian, decorticated and dyskinetic ex-vivo slice preparations were different. Our work may provide a pipeline for diagnosis of basal ganglia pathology from spiking data as well as for the design pharmacological treatments.

Data-driven multiscale computational models of cortical and subcortical regions

Data-driven computational models of neurons, synapses, microcircuits, and mesocircuits have become essential tools in modern brain research. The goal of these multiscale models is to integrate and synthesize information from different levels of brain organization, from cellular properties, dendritic excitability, and synaptic dynamics to microcircuits, mesocircuits, and ultimately behavior. This article surveys recent advances in the genesis of data-driven computational models of mammalian neural networks in cortical and subcortical areas. I discuss the challenges and opportunities in developing data-driven multiscale models, including the need for interdisciplinary collaborations, the importance of model validation and comparison, and the potential impact on basic and translational neuroscience research. Finally, I highlight future directions and emerging technologies that will enable more comprehensive and predictive data-driven models of brain function and dysfunction.

Synaptic and dendritic architecture of different types of hippocampal somatostatin interneurons

GABAergic inhibitory neurons fundamentally shape the activity and plasticity of cortical circuits. A major subset of these neurons contains somatostatin (SOM); these cells play crucial roles in neuroplasticity, learning, and memory in many brain areas including the hippocampus, and are implicated in several neuropsychiatric diseases and neurodegenerative disorders. Two main types of SOM-containing cells in area CA1 of the hippocampus are oriens-lacunosum-moleculare (OLM) cells and hippocampo-septal (HS) cells. These cell types show many similarities in their soma-dendritic architecture, but they have different axonal targets, display different activity patterns in vivo, and are thought to have distinct network functions. However, a complete understanding of the functional roles of these interneurons requires a precise description of their intrinsic computational properties and their synaptic interactions. In the current study we generated, analyzed, and make available several key data sets that enable a quantitative comparison of various anatomical and physiological properties of OLM and HS cells in mouse. The data set includes detailed scanning electron microscopy (SEM)-based 3D reconstructions of OLM and HS cells along with their excitatory and inhibitory synaptic inputs. Combining this core data set with other anatomical data, patch-clamp electrophysiology, and compartmental modeling, we examined the precise morphological structure, inputs, outputs, and basic physiological properties of these cells. Our results highlight key differences between OLM and HS cells, particularly regarding the density and distribution of their synaptic inputs and mitochondria. For example, we estimated that an OLM cell receives about 8,400, whereas an HS cell about 15,600 synaptic inputs, about 16% of which are GABAergic. Our data and models provide insight into the possible basis of the different functionality of OLM and HS cell types and supply essential information for more detailed functional models of these neurons and the hippocampal network.

Grid cells: the missing link in understanding Parkinson's disease?

The mechanisms underlying Parkinson's disease (PD) are complex and not fully understood, and the box-and-arrow model among other current models present significant challenges. This paper explores the potential role of the allocentric brain and especially its grid cells in several PD motor symptoms, including bradykinesia, kinesia paradoxa, freezing of gait, the bottleneck phenomenon, and their dependency on cueing. It is argued that central hubs, like the locus coeruleus and the pedunculopontine nucleus, often narrowly interpreted in the context of PD, play an equally important role in governing the allocentric brain as the basal ganglia. Consequently, the motor and secondary motor (e.g., spatially related) symptoms of PD linked with dopamine depletion may be more closely tied to erroneous computation by grid cells than to the basal ganglia alone. Because grid cells and their associated central hubs introduce both spatial and temporal information to the brain influencing velocity perception they may cause bradykinesia or hyperkinesia as well. In summary, PD motor symptoms may primarily be an allocentric disturbance resulting from virtual faulty computation by grid cells revealed by dopamine depletion in PD.

Cell-Type Specific Connectivity of Whisker-Related Sensory and Motor Cortical Input to Dorsal Striatum

The anterior dorsolateral striatum (DLS) is heavily innervated by convergent excitatory projections from the primary motor (M1) and sensory cortex (S1) and considered an important site of sensorimotor integration. M1 and S1 corticostriatal synapses have functional differences in their connection strength with striatal spiny projection neurons (SPNs) and fast-spiking interneurons (FSIs) in the DLS and, as a result, exert distinct influences on sensory-guided behaviors. In the present study, we tested whether M1 and S1 inputs exhibit differences in the subcellular anatomical distribution of striatal neurons. We injected adeno-associated viral vectors encoding spaghetti monster fluorescent proteins (sm.FPs) into M1 and S1 in male and female mice and used confocal microscopy to generate 3D reconstructions of corticostriatal inputs to single identified SPNs and FSIs obtained through ex vivo patch clamp electrophysiology. We found that M1 and S1 dually innervate SPNs and FSIs; however, there is a consistent bias towards the M1 input in SPNs that is not found in FSIs. In addition, M1 and S1 inputs were distributed similarly across the proximal, medial, and distal regions of SPN and FSI dendrites. Notably, closely localized M1 and S1 clusters of inputs were more prevalent in SPNs than FSIs, suggesting that cortical inputs are integrated through cell-type specific mechanisms. Our results suggest that the stronger functional connectivity from M1 to SPNs compared to S1, as previously observed, is due to a higher quantity of synaptic inputs. Our results have implications for how sensorimotor integration is performed in the striatum through cell-specific differences in corticostriatal connections.

The Basal Ganglia Downstream Control of Action - An Evolutionarily Conserved Strategy

The motor areas of the cortex and the basal ganglia both contribute to determining which motor actions will be recruited at any moment in time, and their functions are intertwined. Here, we review the basal ganglia mechanisms underlying the selection of behavior of the downstream control of motor centers in the midbrain and brainstem and show that the basic organization of the forebrain motor system is evolutionarily conserved throughout vertebrate phylogeny. The output level of the basal ganglia (e.g. substantia nigra pars reticulata) has GABAergic neurons that are spontaneously active at rest and inhibit a number of specific motor centers, each of which can be relieved from inhibition if the inhibitory output neurons themselves become inhibited. The motor areas of the cortex act partially via the dorsolateral striatum (putamen), which has specific modules for the forelimb, hindlimb, trunk, etc. Each module operates in turn through the two types of striatal projection neurons that control the output modules of the basal ganglia and thereby the downstream motor centers. The mechanisms for lateral inhibition in the striatum are reviewed as well as other striatal mechanisms contributing to action selection. The motor cortex also exerts a direct excitatory action on specific motor centers. An overview is given of the basal ganglia control exerted on the different midbrain/brainstem motor centers, and the efference copy information fed back via the thalamus to the striatum and cortex, which is of importance for the planning of future movements.

Distributed dopaminergic signaling in the basal ganglia and its relationship to motor disability in Parkinson's disease

The degeneration of mesencephalic dopaminergic neurons that innervate the basal ganglia is responsible for the cardinal motor symptoms of Parkinson's disease (PD). It has been thought that loss of dopaminergic signaling in one basal ganglia region - the striatum - was solely responsible for the network pathophysiology causing PD motor symptoms. While our understanding of dopamine (DA)'s role in modulating striatal circuitry has deepened in recent years, it also has become clear that it acts in other regions of the basal ganglia to influence movement. Underscoring this point, examination of a new progressive mouse model of PD shows that striatal dopamine DA depletion alone is not sufficient to induce parkinsonism and that restoration of extra-striatal DA signaling attenuates parkinsonian motor deficits once they appear. This review summarizes recent advances in the effort to understand basal ganglia circuitry, its modulation by DA, and how its dysfunction drives PD motor symptoms.

Spiny projection neurons exhibit transcriptional signatures within subregions of the dorsal striatum

The dorsal striatum is organized into functional territories defined by corticostriatal inputs onto both direct and indirect spiny projection neurons (SPNs), the major cell types within the striatum. In addition to circuit connectivity, striatal domains are likely defined by the spatially determined transcriptomes of SPNs themselves. To identify cell-type-specific spatiomolecular signatures of direct and indirect SPNs within dorsomedial, dorsolateral, and ventrolateral dorsal striatum, we used RNA profiling in situ hybridization with probes to >98% of protein coding genes. We demonstrate that the molecular identity of SPNs is mediated by hundreds of differentially expressed genes across territories of the striatum, revealing extraordinary heterogeneity in the expression of genes that mediate synaptic function in both direct and indirect SPNs. This deep insight into the complex spatiomolecular organization of the striatum provides a foundation for understanding both normal striatal function and for dissecting region-specific dysfunction in disorders of the striatum.

A universal workflow for creation, validation, and generalization of detailed neuronal models

Detailed single-neuron modeling is widely used to study neuronal functions. While cellular and functional diversity across the mammalian cortex is vast, most of the available computational tools focus on a limited set of specific features characteristic of a single neuron. Here, we present a generalized automated workflow for the creation of robust electrical models and illustrate its performance by building cell models for the rat somatosensory cortex. Each model is based on a 3D morphological reconstruction and a set of ionic mechanisms. We use an evolutionary algorithm to optimize neuronal parameters to match the electrophysiological features extracted from experimental data. Then we validate the optimized models against additional stimuli and assess their generalizability on a population of similar morphologies. Compared to the state-of-the-art canonical models, our models show 5-fold improved generalizability. This versatile approach can be used to build robust models of any neuronal type.

The roles of surround inhibition for the intrinsic function of the striatum, analyzed in silico

The basal ganglia are important for action initiation, selection, and motor learning. The input level, the striatum, receives input preferentially from the cortex and thalamus and is to 95% composed of striatal projection neurons (SPNs) with sparse GABAergic collaterals targeting distal dendrites of neighboring SPNs, in a distance-dependent manner. The remaining 5% are GABAergic and cholinergic interneurons. Our aim here is to investigate the role of surround inhibition for the intrinsic function of the striatum. Large-scale striatal networks of 20 to 40 thousand neurons were simulated with detailed multicompartmental models of different cell types, corresponding to the size of a module of the dorsolateral striatum, like the forelimb area (mouse). The effect of surround inhibition on dendritic computation and network activity was investigated, while groups of SPNs were activated. The SPN-induced surround inhibition in distal dendrites shunted effectively the corticostriatal EPSPs. The size of dendritic plateau-like potentials within the specific dendritic segment was both reduced and enhanced by inhibition, due to the hyperpolarized membrane potential of SPNs and the reversal-potential of GABA. On a population level, the competition between two subpopulations of SPNs was found to depend on the distance between the two units, the size of each unit, the activity level in each subgroup and the dopaminergic modulation of the dSPNs and iSPNs. The SPNs provided the dominating source of inhibition within the striatum, while the fast-spiking interneuron mainly had an initial effect due to short-term synaptic plasticity as shown in with ablation of the synaptic interaction.

Waxholm Space atlas of the rat brain: a 3D atlas supporting data analysis and integration

Volumetric brain atlases are increasingly used to integrate and analyze diverse experimental neuroscience data acquired from animal models, but until recently a publicly available digital atlas with complete coverage of the rat brain has been missing. Here we present an update of the Waxholm Space rat brain atlas, a comprehensive open-access volumetric atlas resource. This brain atlas features annotations of 222 structures, of which 112 are new and 57 revised compared to previous versions. It provides a detailed map of the cerebral cortex, hippocampal region, striatopallidal areas, midbrain dopaminergic system, thalamic cell groups, the auditory system and main fiber tracts. We document the criteria underlying the annotations and demonstrate how the atlas with related tools and workflows can be used to support interpretation, integration, analysis and dissemination of experimental rat brain data.

A GPU-based computational framework that bridges neuron simulation and artificial intelligence

Biophysically detailed multi-compartment models are powerful tools to explore computational principles of the brain and also serve as a theoretical framework to generate algorithms for artificial intelligence (AI) systems. However, the expensive computational cost severely limits the applications in both the neuroscience and AI fields. The major bottleneck during simulating detailed compartment models is the ability of a simulator to solve large systems of linear equations. Here, we present a novel Dendritic Hierarchical Scheduling (DHS) method to markedly accelerate such a process. We theoretically prove that the DHS implementation is computationally optimal and accurate. This GPU-based method performs with 2-3 orders of magnitude higher speed than that of the classic serial Hines method in the conventional CPU platform. We build a DeepDendrite framework, which integrates the DHS method and the GPU computing engine of the NEURON simulator and demonstrate applications of DeepDendrite in neuroscience tasks. We investigate how spatial patterns of spine inputs affect neuronal excitability in a detailed human pyramidal neuron model with 25,000 spines. Furthermore, we provide a brief discussion on the potential of DeepDendrite for AI, specifically highlighting its ability to enable the efficient training of biophysically detailed models in typical image classification tasks.

Rethinking the network determinants of motor disability in Parkinson's disease

For roughly the last 30 years, the notion that striatal dopamine (DA) depletion was the critical determinant of network pathophysiology underlying the motor symptoms of Parkinson's disease (PD) has dominated the field. While the basal ganglia circuit model underpinning this hypothesis has been of great heuristic value, the hypothesis itself has never been directly tested. Moreover, studies in the last couple of decades have made it clear that the network model underlying this hypothesis fails to incorporate key features of the basal ganglia, including the fact that DA acts throughout the basal ganglia, not just in the striatum. Underscoring this point, recent work using a progressive mouse model of PD has shown that striatal DA depletion alone is not sufficient to induce parkinsonism and that restoration of extra-striatal DA signaling attenuates parkinsonian motor deficits once they appear. Given the broad array of discoveries in the field, it is time for a new model of the network determinants of motor disability in PD.

Compartment-specific dendritic information processing in striatal cholinergic interneurons is reconfigured by peptide neuromodulation

Cholinergic interneurons are central hubs of the striatal neuronal network, controlling information processing in a behavioral-state-dependent manner. It remains unknown, however, how such state transitions influence the integrative properties of these neurons. To address this, we made simultaneous somato-dendritic recordings from identified rodent cholinergic interneurons, revealing that action potentials are initiated at dendritic sites because of a dendritic axonal origin. Functionally, this anatomical arrangement ensured that the action potential initiation threshold was lowest at axon-bearing dendritic sites, a privilege efficacy powerfully accentuated at the hyperpolarized membrane potentials achieved in cholinergic interneurons following salient behavioral stimuli. Experimental analysis revealed the voltage-dependent attenuation of the efficacy of non-axon-bearing dendritic excitatory input was mediated by the recruitment of dendritic potassium channels, a regulatory mechanism that, in turn, was controlled by the pharmacological activation of neurokinin receptors. Together, these results indicate that the neuropeptide microenvironment dynamically controls state- and compartment-dependent dendritic information processing in striatal cholinergic interneurons.

Striatal Subregion-selective Dysregulated Dopamine Receptor-mediated Intracellular Signaling in a Model of DOPA-responsive Dystonia

Although the mechanisms underlying dystonia are largely unknown, dystonia is often associated with abnormal dopamine neurotransmission. DOPA-responsive dystonia (DRD) is a prototype disorder for understanding dopamine dysfunction in dystonia because it is caused by mutations in genes necessary for the synthesis of dopamine and alleviated by the indirect-acting dopamine agonist l-DOPA. Although adaptations in striatal dopamine receptor-mediated intracellular signaling have been studied extensively in models of Parkinson's disease, another movement disorders associated with dopamine deficiency, little is known about dopaminergic adaptations in dystonia. To identify the dopamine receptor-mediated intracellular signaling associated with dystonia, we used immunohistochemistry to quantify striatal protein kinase A activity and extracellular signal-related kinase (ERK) phosphorylation after dopaminergic challenges in a knockin mouse model of DRD. l-DOPA treatment induced the phosphorylation of both protein kinase A substrates and ERK largely in D1 dopamine receptor-expressing striatal neurons. As expected, this response was blocked by pretreatment with the D1 dopamine receptor antagonist SCH23390. The D2 dopamine receptor antagonist raclopride also significantly reduced the phosphorylation of ERK; this contrasts with models of parkinsonism in which l-DOPA-induced ERK phosphorylation is not mediated by D2 dopamine receptors. Further, the dysregulated signaling was dependent on striatal subdomains whereby ERK phosphorylation was largely confined to dorsomedial (associative) striatum while the dorsolateral (sensorimotor) striatum was unresponsive. This complex interaction between striatal functional domains and dysregulated dopamine-receptor mediated responses has not been observed in other models of dopamine deficiency, such as parkinsonism, suggesting that regional variation in dopamine-mediated neurotransmission may be a hallmark of dystonia.

In vivo optogenetic inhibition of striatal parvalbumin-reactive interneurons induced genotype-specific changes in neuronal activity without dystonic signs in male DYT1 knock-in mice

The pathophysiology of early-onset torsion dystonia (TOR1A/DYT1) remains unclear. Like 70% of human mutation carriers, rodent models with ΔGAG mutation such as DYT1 knock-in (KI) mice do not show overt dystonia but have subtle sensorimotor deficits and pattern of abnormal synaptic plasticity within the striatal microcircuits. There is evidence that dysfunction of striatal parvalbumin-reactive (Parv+) fast-spiking interneurons (FSIs) can be involved in dystonic signs. To elucidate the relevance of these GABAergic interneurons in the pathophysiology of DYT1 dystonia, we used in vivo optogenetics to specifically inhibit Parv+ and to detect changes in motor behavior and neuronal activity. Optogenetic fibers were bilaterally implanted into the dorsal striatum of male DYT1 KI mice and wild-type (WT) littermates expressing halorhodopsin (eNpHR3.0) in Parv+ interneurons. While stimulations with yellow light pulses for up to 60 min at different pulse durations and interval lengths did not induce abnormal movements, such as dystonic signs, immunohistochemical examinations revealed genotype-dependent differences. In contrast to WT mice, stimulated DYT1 KI showed decreased striatal neuronal activity, that is, less c-Fos reactive neurons, and increased activation of cholinergic interneurons after optogenetic inhibition of Parv+ interneurons. These findings suggest an involvement of Parv+ interneurons in an impaired striatal network in DYT1 KI mice, but at least short-term inhibition of these GABAergic interneurons is not sufficient to trigger a dystonic phenotype, similar to previously shown optogenetic activation of cholinergic interneurons.

Heuristic Tree-Partition-Based Parallel Method for Biophysically Detailed Neuron Simulation

Biophysically detailed neuron simulation is a powerful tool to explore the mechanisms behind biological experiments and bridge the gap between various scales in neuroscience research. However, the extremely high computational complexity of detailed neuron simulation restricts the modeling and exploration of detailed network models. The bottleneck is solving the system of linear equations. To accelerate detailed simulation, we propose a heuristic tree-partition-based parallel method (HTP) to parallelize the computation of the Hines algorithm, the kernel for solving linear equations, and leverage the strong parallel capability of the graphic processing unit (GPU) to achieve further speedup. We formulate the problem of how to get a fine parallel process as a tree-partition problem. Next, we present a heuristic partition algorithm to obtain an effective partition to efficiently parallelize the equation-solving process in detailed simulation. With further optimization on GPU, our HTP method achieves 2.2 to 8.5 folds speedup compared to the state-of-the-art GPU method and 36 to 660 folds speedup compared to the typical Hines algorithm.

Cell-Type Specific Connectivity of Whisker-Related Sensory and Motor Cortical Input to Dorsal Striatum

The anterior dorsolateral striatum (DLS) is heavily innervated by convergent excitatory projections from the primary motor (M1) and sensory cortex (S1) and is considered an important site of sensorimotor integration. M1 and S1 corticostriatal synapses have functional differences in the strength of their connections with striatal spiny projection neurons (SPNs) and fast-spiking interneurons (FSIs) in the DLS, and as a result exert an opposing influence on sensory-guided behaviors. In the present study, we tested whether M1 and S1 inputs exhibit differences in the subcellular anatomical distribution onto striatal neurons. We injected adeno-associated viral vectors encoding spaghetti monster fluorescent proteins (sm.FPs) into M1 and S1, and used confocal microscopy to generate 3D reconstructions of corticostriatal inputs to single identified SPNs and FSIs obtained through ex-vivo patch-clamp electrophysiology. We found that SPNs are less innervated by S1 compared to M1, but FSIs receive a similar number of inputs from both M1 and S1. In addition, M1 and S1 inputs were distributed similarly across the proximal, medial, and distal regions of SPNs and FSIs. Notably, clusters of inputs were prevalent in SPNs but not FSIs. Our results suggest that SPNs have stronger functional connectivity to M1 compared to S1 due to a higher density of synaptic inputs. The clustering of M1 and S1 inputs onto SPNs but not FSIs suggest that cortical inputs are integrated through cell-type specific mechanisms and more generally have implications for how sensorimotor integration is performed in the striatum.

Significance Statement

The dorsolateral striatum (DLS) is a key brain area involved in sensorimotor integration due to its dense innervation by the primary motor (M1) and sensory cortex (S1). However, the quantity and anatomical distribution of these inputs to the striatal cell population has not been well characterized. In this study we demonstrate that corticostriatal projections from M1 and S1 differentially innervate spiny projection neurons (SPNs) and fast-spiking interneurons (FSIs) in the DLS. S1 inputs innervate SPNs less than M1 and are likely to form synaptic clusters in SPNs but not in FSIs. These findings suggest that sensorimotor integration is partly achieved by differences in the synaptic organization of corticostriatal inputs to local striatal microcircuits.

Full-scale scaffold model of the human hippocampus CA1 area

The increasing availability of quantitative data on the human brain is opening new avenues to study neural function and dysfunction, thus bringing us closer and closer to the implementation of digital twin applications for personalized medicine. Here we provide a resource to the neuroscience community: a computational method to generate full-scale scaffold model of human brain regions starting from microscopy images. We have benchmarked the method to reconstruct the CA1 region of a right human hippocampus, which accounts for about half of the entire right hippocampal formation. Together with 3D soma positioning we provide a connectivity matrix generated using a morpho-anatomical connection strategy based on axonal and dendritic probability density functions accounting for morphological properties of hippocampal neurons. The data and algorithms are supplied in a ready-to-use format, suited to implement computational models at different scales and detail.

A striatal circuit balances learned fear in the presence and absence of sensory cues

During fear learning, defensive behaviors like freezing need to be finely balanced in the presence or absence of threat-predicting cues (conditioned stimulus, CS). Nevertheless, the circuits underlying such balancing are largely unknown. Here, we investigate the role of the ventral tail striatum (vTS) in auditory-cued fear learning of male mice. In vivo Ca²⁺ imaging showed that sizable sub-populations of direct (D1R+) and indirect pathway neurons (Adora+) in the vTS responded to footshocks, and to the initiation of movements after freezing; moreover, a sub-population of D1R+ neurons increased its responsiveness to an auditory CS during fear learning. In-vivo optogenetic silencing shows that footshock-driven activity of D1R+ neurons contributes to fear memory formation, whereas Adora+ neurons modulate freezing in the absence of a learned CS. Circuit tracing identified the posterior insular cortex (pInsCx) as an important cortical input to the vTS, and recording of optogenetically evoked EPSCs revealed long-term plasticity with opposite outcomes at the pInsCx synapses onto D1R+ - and Adora+ neurons. Thus, direct- and indirect pathways neurons of the vTS show differential signs of plasticity after fear learning, and balance defensive behaviors in the presence and absence of learned sensory cues.

Solving brain circuit function and dysfunction with computational modeling and optogenetic fMRI

Can we construct a model of brain function that enables an understanding of whole-brain circuit mechanisms underlying neurological disease and use it to predict the outcome of therapeutic interventions? How are pathologies in neurological disease, some of which are observed to have spatial spreading mechanisms, associated with circuits and brain function? In this review, we discuss approaches that have been used to date and future directions that can be explored to answer these questions. By combining optogenetic functional magnetic resonance imaging (fMRI) with computational modeling, cell type-specific, large-scale brain circuit function and dysfunction are beginning to be quantitatively parameterized. We envision that these developments will pave the path for future therapeutics developments based on a systems engineering approach aimed at directly restoring brain function.

Perineuronal Nets in the Dorsomedial Striatum Contribute to Behavioral Dysfunction in Mouse Models of Excessive Repetitive Behavior

Background

Excessive repetitive behavior is a debilitating symptom of several neuropsychiatric disorders. Parvalbumin-positive inhibitory interneurons in the dorsal striatum have been linked to repetitive behavior, and a sizable portion of these cells are surrounded by perineuronal nets (PNNs), specialized extracellular matrix structures. Although PNNs have been associated with plasticity and neuropsychiatric disease, no previous studies have investigated their involvement in excessive repetitive behavior.

Methods

We used histochemistry and confocal imaging to investigate PNNs surrounding parvalbumin-positive cells in the dorsal striatum of 4 mouse models of excessive repetitive behavior (BTBR, Cntnap2, Shank3, prenatal valproate treatment). We then investigated one of these models, the BTBR mouse, in detail, with DiI labeling, in vivo and in vitro recordings, and behavioral analyses. We next degraded PNNs in the dorsomedial striatum (DMS) using the enzyme chondroitinase ABC and assessed dendritic spine density, electrophysiology, and repetitive behavior.

Results

We found a greater percentage of parvalbumin-positive interneurons with PNNs in the DMS of all 4 mouse models of excessive repetitive behavior compared with control mice. In BTBR mice, we found fewer dendritic spines on medium spiny neurons (targets of parvalbumin-positive interneurons) and differences in neuronal oscillations as well as inhibitory postsynaptic potentials compared with control mice. Reduction of DMS PNNs in BTBR mice altered dendritic spine density and inhibitory responses and normalized repetitive behavior.

Conclusions

These findings suggest that cellular abnormalities in the DMS are associated with maladaptive repetitive behaviors and that manipulating PNNs can restore normal levels of repetitive behavior while altering DMS dendritic spines and inhibitory signaling.

A realistic morpho-anatomical connection strategy for modelling full-scale point-neuron microcircuits

The modeling of extended microcircuits is emerging as an effective tool to simulate the neurophysiological correlates of brain activity and to investigate brain dysfunctions. However, for specific networks, a realistic modeling approach based on the combination of available physiological, morphological and anatomical data is still an open issue. One of the main problems in the generation of realistic networks lies in the strategy adopted to build network connectivity. Here we propose a method to implement a neuronal network at single cell resolution by using the geometrical probability volumes associated with pre- and postsynaptic neurites. This allows us to build a network with plausible connectivity properties without the explicit use of computationally intensive touch detection algorithms using full 3D neuron reconstructions. The method has been benchmarked for the mouse hippocampus CA1 area, and the results show that this approach is able to generate full-scale brain networks at single cell resolution that are in good agreement with experimental findings. This geometric reconstruction of axonal and dendritic occupancy, by effectively reflecting morphological and anatomical constraints, could be integrated into structured simulators generating entire circuits of different brain areas facilitating the simulation of different brain regions with realistic models.

The quest for multiscale brain modeling

Addressing the multiscale organization of the brain, which is fundamental to the dynamic repertoire of the organ, remains challenging. In principle, it should be possible to model neurons and synapses in detail and then connect them into large neuronal assemblies to explain the relationship between microscopic phenomena, large-scale brain functions, and behavior. It is more difficult to infer neuronal functions from ensemble measurements such as those currently obtained with brain activity recordings. In this article we consider theories and strategies for combining bottom-up models, generated from principles of neuronal biophysics, with top-down models based on ensemble representations of network activity and on functional principles. These integrative approaches are hoped to provide effective multiscale simulations in virtual brains and neurorobots, and pave the way to future applications in medicine and information technologies.

A Spiking Neural Network Builder for Systematic Data-to-Model Workflow

In building biological neural network models, it is crucial to efficiently convert diverse anatomical and physiological data into parameters of neurons and synapses and to systematically estimate unknown parameters in reference to experimental observations. Web-based tools for systematic model building can improve the transparency and reproducibility of computational models and can facilitate collaborative model building, validation, and evolution. Here, we present a framework to support collaborative data-driven development of spiking neural network (SNN) models based on the Entity-Relationship (ER) data description commonly used in large-scale business software development. We organize all data attributes, including species, brain regions, neuron types, projections, neuron models, and references as tables and relations within a database management system (DBMS) and provide GUI interfaces for data registration and visualization. This allows a robust “business-oriented” data representation that supports collaborative model building and traceability of source information for every detail of a model. We tested this data-to-model framework in cortical and striatal network models by successfully combining data from papers with existing neuron and synapse models and by generating NEST simulation codes for various network sizes. Our framework also helps to check data integrity and consistency and data comparisons across species. The framework enables the modeling of any region of the brain and is being deployed to support the integration of anatomical and physiological datasets from the brain/MINDS project for systematic SNN modeling of the marmoset brain.

Modernizing the NEURON Simulator for Sustainability, Portability, and Performance

The need for reproducible, credible, multiscale biological modeling has led to the development of standardized simulation platforms, such as the widely-used NEURON environment for computational neuroscience. Developing and maintaining NEURON over several decades has required attention to the competing needs of backwards compatibility, evolving computer architectures, the addition of new scales and physical processes, accessibility to new users, and efficiency and flexibility for specialists. In order to meet these challenges, we have now substantially modernized NEURON, providing continuous integration, an improved build system and release workflow, and better documentation. With the help of a new source-to-source compiler of the NMODL domain-specific language we have enhanced NEURON's ability to run efficiently, via the CoreNEURON simulation engine, on a variety of hardware platforms, including GPUs. Through the implementation of an optimized in-memory transfer mechanism this performance optimized backend is made easily accessible to users, providing training and model-development paths from laptop to workstation to supercomputer and cloud platform. Similarly, we have been able to accelerate NEURON's reaction-diffusion simulation performance through the use of just-in-time compilation. We show that these efforts have led to a growing developer base, a simpler and more robust software distribution, a wider range of supported computer architectures, a better integration of NEURON with other scientific workflows, and substantially improved performance for the simulation of biophysical and biochemical models.

Comprehensive Spatial Profile of the Orphan G Protein Coupled Receptor GPRC5B Expression in Mouse Brain

Orphan G Protein Coupled Receptors (GPCRs) are GPCRs whose endogenous ligands are unknown or still debated. Due to the lack of pharmacological modulators, the physiological function of orphan GPCRs is understudied. However, relevant physiological roles associated with orphan GPCRs have been revealed by analysis of animal models and genome wide association studies illuminating an untapped potential for drug discovery. G Protein Coupled Receptor class C Group 5 Member B (GPRC5B) is among the most expressed GPCRs in the central nervous system. Thus, the expression profiling of GPRC5B is an essential step toward understanding GPRC5B function in health and disease. In this study, we generated new GPRC5B polyclonal antibodies and investigated the expression levels of GPRC5B across different organs and brain regions. We identified high levels of GPRC5B glycosylation both in transfected cells and in mouse brain. Moreover, in situ hybridization imaging analysis indicated that Gprc5b was expressed at the highest level in olfactory bulb, hippocampus, cerebellum, and pons. To dissect expression within various neuronal populations, we conducted a comprehensive spatial profiling of Gprc5b across excitatory and inhibitory neuronal types in medial prefrontal cortex, motor cortex, hippocampal regions, hypothalamus, and cerebellum. Overall, we discovered that GABAergic neurons displayed higher Gprc5b expression levels than glutamatergic neurons in most of the analyzed regions with the important exception of the hippocampal dentate gyrus. Overall, the expression analysis of GPRC5B in mouse brain will guide functional studies ultimately positioning GPRC5B in pathophysiological mechanisms and drug discovery.

Non-monotonic effects of GABAergic synaptic inputs on neuronal firing

GABA is generally known as the principal inhibitory neurotransmitter in the nervous system, usually acting by hyperpolarizing membrane potential. However, GABAergic currents sometimes exhibit non-inhibitory effects, depending on the brain region, developmental stage or pathological condition. Here, we investigate the diverse effects of GABA on the firing rate of several single neuron models, using both analytical calculations and numerical simulations. We find that GABAergic synaptic conductance and output firing rate exhibit three qualitatively different regimes as a function of GABA reversal potential, EGABA: monotonically decreasing for sufficiently low EGABA (inhibitory), monotonically increasing for EGABA above firing threshold (excitatory); and a non-monotonic region for intermediate values of EGABA. In the non-monotonic regime, small GABA conductances have an excitatory effect while large GABA conductances show an inhibitory effect. We provide a phase diagram of different GABAergic effects as a function of GABA reversal potential and glutamate conductance. We find that noisy inputs increase the range of EGABA for which the non-monotonic effect can be observed. We also construct a micro-circuit model of striatum to explain observed effects of GABAergic fast spiking interneurons on spiny projection neurons, including non-monotonicity, as well as the heterogeneity of the effects. Our work provides a mechanistic explanation of paradoxical effects of GABAergic synaptic inputs, with implications for understanding the effects of GABA in neural computation and development.

Data-Driven Model of Postsynaptic Currents Mediated by NMDA or AMPA Receptors in Striatal Neurons

The majority of excitatory synapses in the brain uses glutamate as neurotransmitter, and the synaptic transmission is primarily mediated by AMPA and NMDA receptors in postsynaptic neurons. Here, we present data-driven models of the postsynaptic currents of these receptors in excitatory synapses in mouse striatum. It is common to fit two decay time constants to the decay phases of the current profiles but then compute a single weighted mean time constant to describe them. We have shown that this approach does not lead to an improvement in the fitting, and, hence, we present a new model based on the use of both the fast and slow time constants and a numerical calculation of the peak time using Newton's method. Our framework allows for a more accurate description of the current profiles without needing extra data and without overburdening the comptuational costs. The user-friendliness of the method, here implemented in Python, makes it easily applicable to other data sets.

Computational Concepts for Reconstructing and Simulating Brain Tissue

It has previously been shown that it is possible to derive a new class of biophysically detailed brain tissue models when one computationally analyzes and exploits the interdependencies or the multi-modal and multi-scale organization of the brain. These reconstructions, sometimes referred to as digital twins, enable a spectrum of scientific investigations. Building such models has become possible because of increase in quantitative data but also advances in computational capabilities, algorithmic and methodological innovations. This chapter presents the computational science concepts that provide the foundation to the data-driven approach to reconstructing and simulating brain tissue as developed by the EPFL Blue Brain Project, which was originally applied to neocortical microcircuitry and extended to other brain regions. Accordingly, the chapter covers aspects such as a knowledge graph-based data organization and the importance of the concept of a dataset release. We illustrate algorithmic advances in finding suitable parameters for electrical models of neurons or how spatial constraints can be exploited for predicting synaptic connections. Furthermore, we explain how in silico experimentation with such models necessitates specific addressing schemes or requires strategies for an efficient simulation. The entire data-driven approach relies on the systematic validation of the model. We conclude by discussing complementary strategies that not only enable judging the fidelity of the model but also form the basis for its systematic refinements.

Neuron types in the primate striatum: stereological analysis of projection neurons and interneurons in control and parkinsonian monkeys

AIMS

The striatum is mainly composed of projection neurons. It also contains interneurons, which modulate and control striatal output. The aim of the present study was to assess the percentages of projection neurons and interneuron populations in the striatum of control monkeys and of parkinsonian monkeys.

METHODS

Unbiased stereology was used to estimate the volume density of every neuron population in the caudate, putamen and ventral striatum of control monkeys and of monkeys treated with MPTP, which results in striatal dopamine depletion. The various neuron population phenotypes were identified by immunohistochemistry. All analyses were performed within the same subjects using similar processing and analysis parameters, thus allowing for reliable data comparisons.

RESULTS

In control monkeys, the projection neurons, which express the Dopamine-and-cAMP-Regulated-Phosphoprotein, 32-KDa (DARPP-32), were the most abundant: ~86% of the total neurons counted. The interneurons accounted for the remaining 14%. Among the interneurons, those expressing Calretinin were the most abundant (Cr+: ~57%; ~8% of the total striatal neurons counted), followed those expressing Parvalbumin (Pv+: ~18 %; 2.6%), Dinucleotide Phosphate-Diaphorase (NADPH+: ~13 %; 1.8%), Choline Acetyltransferase (ChAT+: ~11%; 1.5%) and Tyrosine Hydroxylase (TH+: ~0.5%; 0.1%). No significant changes in volume densities occurred in any population following dopamine depletion, except for the TH+ interneurons, which increased in parkinsonian non-symptomatic monkeys and even more in symptomatic monkeys.

CONCLUSIONS

These data are relevant for translational studies targeting specific neuron populations of the striatum. The fact that dopaminergic denervation does not cause neuron loss in any population has potential pathophysiological implications.

Biophysical Modeling of Dopaminergic Denervation Landscapes in the Striatum Reveals New Therapeutic Strategy

Parkinson’s disease (PD) results from a loss of dopaminergic neurons. What triggers the break-down of neuronal signaling, and how this might be compensated, is not understood. The age of onset, progression and symptoms vary between patients, and our understanding of the clinical variability remains incomplete. In this study, we investigate this, by characterizing the dopaminergic landscape in healthy and denervated striatum, using biophysical modeling. Based on currently proposed mechanisms, we model three distinct denervation patterns, and show how this affect the dopaminergic network. Depending on the denervation pattern, we show how local and global differences arise in the activity of striatal neurons. Finally, we use the mathematical formalism to suggest a cellular strategy for maintaining normal dopamine (DA) signaling following neuronal denervation. This strategy is characterized by dual enhancement of both the release and uptake capacity of DA in the remaining neurons. Overall, our results derive a new conceptual framework for the impaired dopaminergic signaling related to PD and offers testable predictions for future research directions.

Transient Response of Basal Ganglia Network in Healthy and Low-Dopamine State

The basal ganglia (BG) are crucial for a variety of motor and cognitive functions. Changes induced by persistent low-dopamine (e.g., in Parkinson’s disease; PD) result in aberrant changes in steady-state population activity (β band oscillations) and the transient response of the BG. Typically, a brief cortical stimulation results in a triphasic response in the substantia nigra pars reticulata (SNr; an output of the BG). The properties of the triphasic responses are shaped by dopamine levels. While mechanisms underlying aberrant steady state activity are well studied, it is still unclear which BG interactions are crucial for the aberrant transient responses in the BG. Moreover, it is also unclear whether mechanisms underlying the aberrant changes in steady-state activity and transient response are the same. Here, we used numerical simulations of a network model of BG to identify the key factors that determine the shape of the transient responses. We show that an aberrant transient response of the SNr in the low-dopamine state involves changes in the direct pathway and the recurrent interactions within the globus pallidus externa (GPe) and between GPe and subthalamic nucleus (STN). However, the connections from D2-type spiny projection neurons (D2-SPN) to GPe are most crucial in shaping the transient response and by restoring them to their healthy level, we could restore the shape of transient response even in low-dopamine state. Finally, we show that the changes in BG that result in aberrant transient response are also sufficient to generate pathologic oscillatory activity in the steady state.

Three-Dimensional Spatial Analyses of Cholinergic Neuronal Distributions Across The Mouse Septum, Nucleus Basalis, Globus Pallidus, Nucleus Accumbens, and Caudate-Putamen

Neuronal networks are regulated by three-dimensional spatial and structural properties. Despite robust evidence of functional implications in the modulation of cognition, little is known about the three-dimensional internal organization of cholinergic networks in the forebrain. Cholinergic networks in the forebrain primarily occur in subcortical nuclei, specifically the septum, nucleus basalis, globus pallidus, nucleus accumbens, and the caudate-putamen. Therefore, the present investigation analyzed the three-dimensional spatial organization of 14,000 cholinergic neurons that expressed choline acetyltransferase (ChAT) in these subcortical nuclei of the mouse forebrain. Point process theory and graph signal processing techniques identified three topological principles of organization. First, cholinergic interneuronal distance is not uniform across brain regions. Specifically, in the septum, globus pallidus, nucleus accumbens, and the caudate-putamen, the cholinergic neurons were clustered compared with a uniform random distribution. In contrast, in the nucleus basalis, the cholinergic neurons had a spatial distribution of greater regularity than a uniform random distribution. Second, a quarter of the caudate-putamen is composed of axonal bundles, yet the spatial distribution of cholinergic neurons remained clustered when axonal bundles were accounted for. However, comparison with an inhomogeneous Poisson distribution showed that the nucleus basalis and caudate-putamen findings could be explained by density gradients in those structures. Third, the number of cholinergic neurons varies as a function of the volume of a specific brain region but cell body volume is constant across regions. The results of the present investigation provide topographic descriptions of cholinergic somata distribution and axonal conduits, and demonstrate spatial differences in cognitive control networks. The study provides a comprehensive digital database of the total population of ChAT-positive neurons in the reported structures, with the x,y,z coordinates of each neuron at micrometer resolution. This information is important for future digital cellular atlases and computational models of the forebrain cholinergic system enabling models based on actual spatial geometry.

Autophagy in Spinocerebellar ataxia type 2, a dysregulated pathway, and a target for therapy

Spinocerebellar ataxia type 2 (SCA2) is an incurable and genetic neurodegenerative disorder. The disease is characterized by progressive degeneration of several brain regions, resulting in severe motor and non-motor clinical manifestations. The mutation causing SCA2 disease is an abnormal expansion of CAG trinucleotide repeats in the ATXN2 gene, leading to a toxic expanded polyglutamine segment in the translated ataxin-2 protein. While the genetic cause is well established, the exact mechanisms behind neuronal death induced by mutant ataxin-2 are not yet completely understood. Thus, the goal of this study is to investigate the role of autophagy in SCA2 pathogenesis and investigate its suitability as a target for therapeutic intervention. For that, we developed and characterized a new striatal lentiviral mouse model that resembled several neuropathological hallmarks observed in SCA2 disease, including formation of aggregates, neuronal marker loss, cell death and neuroinflammation. In this new model, we analyzed autophagic markers, which were also analyzed in a SCA2 cellular model and in human post-mortem brain samples. Our results showed altered levels of SQSTM1 and LC3B in cells and tissues expressing mutant ataxin-2. Moreover, an abnormal accumulation of these markers was detected in SCA2 patients' striatum and cerebellum. Importantly, the molecular activation of autophagy, using the compound cordycepin, mitigated the phenotypic alterations observed in disease models. Overall, our study suggests an important role for autophagy in the context of SCA2 pathology, proposing that targeting this pathway could be a potential target to treat SCA2 patients.

Personalized Medicine to Improve Treatment of Dopa-Responsive Dystonia-A Focus on Tyrosine Hydroxylase Deficiency

Dopa-responsive dystonia (DRD) is a rare movement disorder associated with defective dopamine synthesis. This impairment may be due to the fact of a deficiency in GTP cyclohydrolase I (GTPCHI, GCH1 gene), sepiapterin reductase (SR), tyrosine hydroxylase (TH), or 6-pyruvoyl tetrahydrobiopterin synthase (PTPS) enzyme functions. Mutations in GCH1 are most frequent, whereas fewer cases have been reported for individual SR-, PTP synthase-, and TH deficiencies. Although termed DRD, a subset of patients responds poorly to L-DOPA. As this is regularly observed in severe cases of TH deficiency (THD), there is an urgent demand for more adequate or personalized treatment options. TH is a key enzyme that catalyzes the rate-limiting step in catecholamine biosynthesis, and THD patients often present with complex and variable phenotypes, which results in frequent misdiagnosis and lack of appropriate treatment. In this expert opinion review, we focus on THD pathophysiology and ongoing efforts to develop novel therapeutics for this rare disorder. We also describe how different modeling approaches can be used to improve genotype to phenotype predictions and to develop in silico testing of treatment strategies. We further discuss the current status of mathematical modeling of catecholamine synthesis and how such models can be used together with biochemical data to improve treatment of DRD patients.

Dopaminergic and Cholinergic Modulation of Large Scale Networks in silico Using Snudda

Neuromodulation is present throughout the nervous system and serves a critical role for circuit function and dynamics. The computational investigations of neuromodulation in large scale networks require supportive software platforms. Snudda is a software for the creation and simulation of large scale networks of detailed microcircuits consisting of multicompartmental neuron models. We have developed an extension to Snudda to incorporate neuromodulation in large scale simulations. The extended Snudda framework implements neuromodulation at the level of single cells incorporated into large-scale microcircuits. We also developed Neuromodcell, a software for optimizing neuromodulation in detailed multicompartmental neuron models. The software adds parameters within the models modulating the conductances of ion channels and ionotropic receptors. Bath application of neuromodulators is simulated and models which reproduce the experimentally measured effects are selected. In Snudda, we developed an extension to accommodate large scale simulations of neuromodulation. The simulator has two modes of simulation – denoted replay and adaptive. In the replay mode, transient levels of neuromodulators can be defined as a time-varying function which modulates the receptors and ion channels within the network in a cell-type specific manner. In the adaptive mode, spiking neuromodulatory neurons are connected via integrative modulating mechanisms to ion channels and receptors. Both modes of simulating neuromodulation allow for simultaneous modulation by several neuromodulators that can interact dynamically with each other. Here, we used the Neuromodcell software to simulate dopaminergic and muscarinic modulation of neurons from the striatum. We also demonstrate how to simulate different neuromodulatory states with dopamine and acetylcholine using Snudda. All software is freely available on Github, including tutorials on Neuromodcell and Snudda-neuromodulation.

Making decisions in the dark basement of the brain: A look back at the GPR model of action selection and the basal ganglia

How does your brain decide what you will do next? Over the past few decades compelling evidence has emerged that the basal ganglia, a collection of nuclei in the fore- and mid-brain of all vertebrates, are vital to action selection. Gurney, Prescott, and Redgrave published an influential computational account of this idea in Biological Cybernetics in 2001. Here we take a look back at this pair of papers, outlining the "GPR" model contained therein, the context of that model's development, and the influence it has had over the past twenty years. Tracing its lineage into models and theories still emerging now, we are encouraged that the GPR model is that rare thing, a computational model of a brain circuit whose advances were directly built on by others.

Predicting Synaptic Connectivity for Large-Scale Microcircuit Simulations Using Snudda

Simulation of large-scale networks of neurons is an important approach to understanding and interpreting experimental data from healthy and diseased brains. Owing to the rapid development of simulation software and the accumulation of quantitative data of different neuronal types, it is possible to predict both computational and dynamical properties of local microcircuits in a ‘bottom-up’ manner. Simulated data from these models can be compared with experiments and ‘top-down’ modelling approaches, successively bridging the scales. Here we describe an open source pipeline, using the software Snudda, for predicting microcircuit connectivity and for setting up simulations using the NEURON simulation environment in a reproducible way. We also illustrate how to further ‘curate’ data on single neuron morphologies acquired from public databases. This model building pipeline was used to set up a first version of a full-scale cellular level model of mouse dorsal striatum. Model components from that work are here used to illustrate the different steps that are needed when modelling subcortical nuclei, such as the basal ganglia.

Comparative study of neuron-specific promoters in mouse brain transduced by intravenously administered AAV-PHP.eB

Adeno-associated virus (AAV)- PHP.B and AAV-PHP.eB (PHP.eB), a capsid variant of AAV serotype 9, efficiently penetrates the mouse blood-brain barrier and predominantly infects neurons. Thus, the PHP.B / PHP.eB capsid and a neuron-specific promoter is a reasonable combination for effective neuronal transduction. However, the transduction characteristics of intravenously administered PHP.B / PHP.eB carrying different neuron-specific promoters have not been studied systematically. In this study, using an intravenous infusion of PHP.eB in mice, we performed a comparative study of the ubiquitous CBh and three neuron-specific promoters, the Ca²⁺/calmodulin-dependent kinase subunit α (CaMKII) promoter, neuron-specific enolase (NSE) promoter, and synapsin I with a minimal CMV sequence (SynI-minCMV) promoter. Expression levels of a transgene by three neuron-specific promoters were comparable to or higher than those of the CBh promoter. Among the promoters examined, the NSE promoter showed the highest transgene expression. All neuron-specific promoters were activated specifically in the neurons. PHP.eB carrying the CaMKII promoter, which is generally believed to exert its function exclusively in the excitatory neurons, transduced both the excitatory and inhibitory neurons without bias, whereas PHP.eB with the NSE and SynI-minCMV promoters transduced neurons with significant bias toward inhibitory neurons. These results are useful in neuron-targeted broad transgene expression through systemic infusion of blood-brain-barrier-penetrating AAV vectors carrying the neuron-specific promoter.

Er81 Transcription Factor Fine-Tunes Striatal Cholinergic Interneuron Activity and Drives Habit Formation

The molecular mechanisms tuning cholinergic interneuron (CIN) activity, although crucial for striatal function and behavior, remain largely unexplored. Previous studies report that the Etv1/Er81 transcription factor is vital for regulating neuronal maturation and activity. While Er81 is known to be expressed in the striatum during development, its specific role in defining CIN properties and the resulting consequences on striatal function is unknown. We report here that Er81 is expressed in CINs and its specific ablation leads to prominent changes in their molecular, morphologic, and electrophysiological features. In particular, the lack of Er81 amplifies intrinsic delayed-rectifier and hyperpolarization-activated currents, which subsequently alters the tonic and phasic activity of CINs. We further reveal that Er81 expression is required for normal CIN pause and time-locked responses to sensorimotor inputs in awake mice. Overall, this study uncovers a new cell type-specific control of CIN function in the striatum which drives habit formation in adult male mice.SIGNIFICANCE STATEMENT Although previous studies have shown that cholinergic interneurons drive striatal activity and habit formation, the underlying molecular mechanisms controlling their function are unknown. Here we reveal that key cholinergic interneuron physiological properties are controlled by Er81, a transcription factor regulating neuronal activity and development in a cell-specific manner. Moreover, our findings uncover a link between the Er81-dependent molecular control of cholinergic interneuron function and habit formation in mice. These insights will contribute to the future enhancement of our understanding of disorders that involve behavioral inflexibility, such as autism and addiction.

On the structural connectivity of large-scale models of brain networks at cellular level

The brain’s structural connectivity plays a fundamental role in determining how neuron networks generate, process, and transfer information within and between brain regions. The underlying mechanisms are extremely difficult to study experimentally and, in many cases, large-scale model networks are of great help. However, the implementation of these models relies on experimental findings that are often sparse and limited. Their predicting power ultimately depends on how closely a model’s connectivity represents the real system. Here we argue that the data-driven probabilistic rules, widely used to build neuronal network models, may not be appropriate to represent the dynamics of the corresponding biological system. To solve this problem, we propose to use a new mathematical framework able to use sparse and limited experimental data to quantitatively reproduce the structural connectivity of biological brain networks at cellular level.

HippoUnit: A software tool for the automated testing and systematic comparison of detailed models of hippocampal neurons based on electrophysiological data

Anatomically and biophysically detailed data-driven neuronal models have become widely used tools for understanding and predicting the behavior and function of neurons. Due to the increasing availability of experimental data from anatomical and electrophysiological measurements as well as the growing number of computational and software tools that enable accurate neuronal modeling, there are now a large number of different models of many cell types available in the literature. These models were usually built to capture a few important or interesting properties of the given neuron type, and it is often unknown how they would behave outside their original context. In addition, there is currently no simple way of quantitatively comparing different models regarding how closely they match specific experimental observations. This limits the evaluation, re-use and further development of the existing models. Further, the development of new models could also be significantly facilitated by the ability to rapidly test the behavior of model candidates against the relevant collection of experimental data. We address these problems for the representative case of the CA1 pyramidal cell of the rat hippocampus by developing an open-source Python test suite, which makes it possible to automatically and systematically test multiple properties of models by making quantitative comparisons between the models and electrophysiological data. The tests cover various aspects of somatic behavior, and signal propagation and integration in apical dendrites. To demonstrate the utility of our approach, we applied our tests to compare the behavior of several different rat hippocampal CA1 pyramidal cell models from the ModelDB database against electrophysiological data available in the literature, and evaluated how well these models match experimental observations in different domains. We also show how we employed the test suite to aid the development of models within the European Human Brain Project (HBP), and describe the integration of the tests into the validation framework developed in the HBP, with the aim of facilitating more reproducible and transparent model building in the neuroscience community.

Toward Closed-Loop Electrical Stimulation of Neuronal Systems: A Review

Biological neuronal cells communicate using neurochemistry and electrical signals. The same phenomena also allow us to probe and manipulate neuronal systems and communicate with them. Neuronal system malfunctions cause a multitude of symptoms and functional deficiencies that can be assessed and sometimes alleviated by electrical stimulation. Our working hypothesis is that real-time closed-loop full-duplex measurement and stimulation paradigms can provide more in-depth insight into neuronal networks and enhance our capability to control diseases of the nervous system. In this study, we review extracellular electrical stimulation methods used in in vivo, in vitro, and in silico neuroscience research and in the clinic (excluding methods mainly aimed at neuronal growth and other similar effects) and highlight the potential of closed-loop measurement and stimulation systems. A multitude of electrical stimulation and measurement-based methods are widely used in research and the clinic. Closed-loop methods have been proposed, and some are used in the clinic. However, closed-loop systems utilizing more complex measurement analysis and adaptive stimulation systems, such as artificial intelligence systems connected to biological neuronal systems, do not yet exist. Our review promotes the research and development of intelligent paradigms aimed at meaningful communications between neuronal and information and communications technology systems, "dialogical paradigms," which have the potential to take neuroscience and clinical methods to a new level.

The Tail of the Striatum: From Anatomy to Connectivity and Function

The dorsal striatum, the largest subcortical structure of the basal ganglia, is critical in controlling motor, procedural, and reinforcement-based behaviors. Although in mammals the striatum extends widely along the rostro-caudal axis, current knowledge and derived theories about its anatomo-functional organization largely rely on results obtained from studies of its rostral sectors, leading to potentially oversimplified working models of the striatum as a whole. Recent findings indicate that the extreme caudal part of the striatum, also referred to as the tail of striatum (TS), represents an additional functional domain. Here, we provide an overview of past and recent studies revealing that the TS displays a heterogeneous cell-type-specific organization, and a unique input-output connectivity, which poises the TS as an integrator of sensory processing.

Modeling nucleus accumbens : A Computational Model from Single Cell to Circuit Level

Nucleus accumbens is part of the neural structures required for reward based learning and cognitive processing of motivation. Understanding its cellular dynamics and its role in basal ganglia circuits is important not only in diagnosing behavioral disorders and psychiatric problems as addiction and depression but also for developing therapeutic treatments for them. Building a computational model would expand our comprehension of nucleus accumbens. In this work, we are focusing on establishing a model of nucleus accumbens which has not been considered as much as dorsal striatum in computational neuroscience. We will begin by modeling the behavior of single cells and then build a holistic model of nucleus accumbens considering the effect of synaptic currents. We will verify the validity of the model by showing the consistency of simulation results with the empirical data. Furthermore, the simulation results reveal the joint effect of cortical stimulation and dopaminergic modulation on the activity of medium spiny neurons. This effect differentiates with the type of dopamine receptors.