Data-driven multiscale model of macaque auditory thalamocortical circuits reproduces in vivo dynamics

A connectome manipulation framework for the systematic and reproducible study of structure-function relationships through simulations

Synaptic connectivity at the neuronal level is characterized by highly nonrandom features. Hypotheses about their role can be developed by correlating structural metrics to functional features. But, to prove causation, manipulations of connectivity would have to be studied. However, the fine-grained scale at which nonrandom trends are expressed makes this approach challenging to pursue experimentally. Simulations of neuronal networks provide an alternative route to study arbitrarily complex manipulations in morphologically and biophysically detailed models. Here, we present Connectome-Manipulator, a Python framework for rapid connectome manipulations of large-scale network models in Scalable Open Network Architecture TemplAte (SONATA) format. In addition to creating or manipulating the connectome of a model, it provides tools to fit parameters of stochastic connectivity models against existing connectomes. This enables rapid replacement of any existing connectome with equivalent connectomes at different levels of complexity, or transplantation of connectivity features from one connectome to another, for systematic study. We employed the framework in the detailed model of the rat somatosensory cortex in two exemplary use cases: transplanting interneuron connectivity trends from electron microscopy data and creating simplified connectomes of excitatory connectivity. We ran a series of network simulations and found diverse shifts in the activity of individual neuron populations causally linked to these manipulations.

Spike frequency adaptation in primate lateral prefrontal cortex neurons results from interplay between intrinsic properties and circuit dynamics

Cortical neurons in brain slices display intrinsic spike frequency adaptation (I-SFA) to constant current inputs, while extracellular recordings show extrinsic SFA (E-SFA) during sustained visual stimulation. Inferring how I-SFA contributes to E-SFA during behavior is challenging due to the isolated nature of slice recordings. To address this, we recorded macaque lateral prefrontal cortex (LPFC) neurons in vivo during a visually guided saccade task and in vitro in brain slices. Broad-spiking (BS) putative pyramidal cells and narrow-spiking (NS) putative inhibitory interneurons exhibit both E-SFA and I-SFA. Developing a data-driven hybrid circuit model comprising NS model neurons receiving BS input reveals that NS model neurons exhibit longer SFA than observed in vivo; however, adding feedforward inhibition corrects this in a manner dependent on I-SFA. Identification of this circuit motif shaping E-SFA in LPFC highlights the roles of both intrinsic and network mechanisms in neural activity underlying behavior.

Evaluation and comparison of methods for neuronal parameter optimization using the Neuroptimus software framework

Finding optimal parameters for detailed neuronal models is a ubiquitous challenge in neuroscientific research. In recent years, manual model tuning has been gradually replaced by automated parameter search using a variety of different tools and methods. However, using most of these software tools and choosing the most appropriate algorithm for a given optimization task require substantial technical expertise, which prevents the majority of researchers from using these methods effectively. To address these issues, we developed a generic platform (called Neuroptimus) that allows users to set up neural parameter optimization tasks via a graphical interface, and to solve these tasks using a wide selection of state-of-the-art parameter search methods implemented by five different Python packages. Neuroptimus also offers several features to support more advanced usage, including the ability to run most algorithms in parallel, which allows it to take advantage of high-performance computing architectures. We used the common interface provided by Neuroptimus to conduct a detailed comparison of more than twenty different algorithms (and implementations) on six distinct benchmarks that represent typical scenarios in neuronal parameter search. We quantified the performance of the algorithms in terms of the best solutions found and in terms of convergence speed. We identified several algorithms, including covariance matrix adaptation evolution strategy and particle swarm optimization, that consistently, without any fine-tuning, found good solutions in all of our use cases. By contrast, some other algorithms including all local search methods provided good solutions only for the simplest use cases, and failed completely on more complex problems. We also demonstrate the versatility of Neuroptimus by applying it to an additional use case that involves tuning the parameters of a subcellular model of biochemical pathways. Finally, we created an online database that allows uploading, querying and analyzing the results of optimization runs performed by Neuroptimus, which enables all researchers to update and extend the current benchmarking study. The tools and analysis we provide should aid members of the neuroscience community to apply parameter search methods more effectively in their research.

A Model-Driven Meta-Analysis Supports the Emerging Consensus View that Inhibitory Neurons Dominate BOLD-fMRI Responses

Functional magnetic resonance imaging (fMRI) is a pivotal tool for mapping neuronal activity in the brain. Traditionally, the observed hemodynamic changes are assumed to reflect the activity of the most common neuronal type: excitatory neurons. In contrast, recent experiments, using optogenetic techniques, suggest that the fMRI-signal instead reflects the activity of inhibitory interneurons. However, these data paint a complex picture, with numerous regulatory interactions, and where the different experiments display many qualitative differences. It is therefore not trivial how to quantify the relative contributions of the different cell types and to combine all observations into a unified theory. To address this, we present a new model-driven meta-analysis, which provides a unified and quantitative explanation for all data. This model-driven analysis allows for quantification of the relative contribution of different cell types: the contribution to the BOLD-signal from the excitatory cells is <20 % and 50-80 % comes from the interneurons. Our analysis also provides a mechanistic explanation for the observed experiment-to-experiment differences, e.g. a biphasic vascular response dependent on different stimulation intensities and an emerging secondary post-stimulation peak during longer stimulations. In summary, our study provides a new, emerging consensus-view supporting the larger role of interneurons in fMRI.

Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons

Understanding the brain requires studying its multiscale interactions from molecules to networks. The increasing availability of large-scale datasets detailing brain circuit composition, connectivity, and activity is transforming neuroscience. However, integrating and interpreting this data remains challenging. Concurrently, advances in supercomputing and sophisticated modeling tools now enable the development of highly detailed, large-scale biophysical circuit models. These mechanistic multiscale models offer a method to systematically integrate experimental data, facilitating investigations into brain structure, function, and disease. This review, based on a Society for Neuroscience 2024 MiniSymposium, aims to disseminate recent advances in large-scale mechanistic modeling to the broader community. It highlights (1) examples of current models for various brain regions developed through experimental data integration; (2) their predictive capabilities regarding cellular and circuit mechanisms underlying experimental recordings (e.g., membrane voltage, spikes, local-field potential, electroencephalography/magnetoencephalography) and brain function; and (3) their use in simulating biomarkers for brain diseases like epilepsy, depression, schizophrenia, and Parkinson's, aiding in understanding their biophysical underpinnings and developing novel treatments. The review showcases state-of-the-art models covering hippocampus, somatosensory, visual, motor, auditory cortical, and thalamic circuits across species. These models predict neural activity at multiple scales and provide insights into the biophysical mechanisms underlying sensation, motor behavior, brain signals, neural coding, disease, pharmacological interventions, and neural stimulation. Collaboration with experimental neuroscientists and clinicians is essential for the development and validation of these models, particularly as datasets grow. Hence, this review aims to foster interest in detailed brain circuit models, leading to cross-disciplinary collaborations that accelerate brain research.

Is there a ubiquitous spectrolaminar motif of local field potential power across primate neocortex?

Mendoza-Halliday, Major et al., 2024 ("The Paper")1 advocates a local field potential (LFP)-based approach to functional identification of cortical layers during "laminar" (simultaneous recordings from all cortical layers) multielectrode recordings in nonhuman primates (NHPs). The Paper describes a "ubiquitous spectrolaminar motif" in the primate neocortex: 1) 75-150 Hz power peaks in the supragranular layers, 2) 10-19 Hz power peaks in the infragranular layers and 3) the crossing point of their laminar power gradients identifies Layer 4 (L4). Identification of L4 is critical in general, but especially for The Paper as the "motif" discovery is couched within a framework whose central hypothesis is that gamma activity originates in the supragranular layers and reflects feedforward activity, while alpha-beta activity originates in the infragranular layers and reflects feedback activity. In an impressive scientific effort, The Paper analyzed laminar data from 14 cortical areas in 2 prior macaque studies and compared them to marmoset, mouse, and human data to further bolster the canonical nature of the motif. Identification of such canonical principles of brain operation is clearly a topic of broad scientific interest. Similarly, a reliable online method for L4 identification would be of broad scientific value for the rapidly increasing use of laminar recordings using numerous evolving technologies. Despite The Paper's strengths, and its potential for scientific impact, a series of concerns that are fundamental to the analysis and interpretation of laminar activity profile data in general, and local field potential (LFP) signals in particular, led us to question its conclusions. We thus evaluated the generality of The Paper's methods and findings using new sets of data comprised of stimulus-evoked laminar response profiles from primary and higher-order auditory cortices (A1 and belt cortex), and primary visual cortex (V1). The rationale for using these areas as a test bed for new methods is that their laminar anatomy and physiology have already been extensively characterized by prior studies, and there is general agreement across laboratories on key matters like L4 identification. Our analyses indicate that The Paper's findings do not generalize well to any of these cortical areas. In particular, we find The Paper's methods for L4 identification to be unreliable. Moreover, both methodological and statistical concerns, outlined below and in the supplement, question the stated prevalence of the motif in The Paper's published dataset. After summarizing our findings and related broader concerns, we briefly critique the evidence from biophysical modeling studies cited to support The Paper's conclusions. While our findings are at odds with the proposition of a ubiquitous spectrolaminar motif in the primate neocortex, The Paper already has, and will continue to spark debate and further experimentation. Hopefully this countervailing presentation will lead to robust collegial efforts to define optimal strategies for applying laminar recording methods in future studies.

Fibroblast growth factor homologous factors: canonical and non-canonical mechanisms of action

Since their discovery nearly 30 years ago, fibroblast growth factor homologous factors (FHFs) are now known to control the functionality of excitable tissues through a range of mechanisms. Nervous and cardiac system dysfunctions are caused by loss- or gain-of-function mutations in FHF genes. The best understood 'canonical' targets for FHF action are voltage-gated sodium channels, and recent studies have expanded the repertoire of ways that FHFs modulate sodium channel gating. Additional 'non-canonical' functions of FHFs in excitable and non-excitable cells, including cancer cells, have been reported over the past dozen years. This review summarizes and evaluates reported canonical and non-canonical FHF functions.

Mechanism of an Intrinsic Oscillation in Rat Geniculate Interneurons

Depolarizing current injections produced a rhythmic bursting of action potentials - a bursting oscillation - in a set of local interneurons in the lateral geniculate nucleus (LGN) of rats. The current dynamics underlying this firing pattern have not been determined, though this cell type constitutes an important cellular component of thalamocortical circuitry, and contributes to both pathologic and non-pathologic brain states. We thus investigated the source of the bursting oscillation using pharmacological manipulations in LGN slices in vitro and in silico. 1. Selective blockade of calcium channel subtypes revealed that high-threshold calcium currentsI L andI P contributed strongly to the oscillation. 2. Increased extracellular K+ concentration (decreased K+currents) eliminated the oscillation. 3. Selective blockade of K+ channel subtypes demonstrated that the calcium-sensitive potassium current (I A H P ) was of primary importance. A morphologically simplified, multicompartment model of the thalamic interneuron characterized the oscillation as follows: 1. The low-threshold calcium currentI T provided the strong initial burst characteristic of the oscillation. 2. Alternating fluxes through high-threshold calcium channels andI A H P then provided the continuing oscillation's burst and interburst periods respectively. This interplay betweenI L andI A H P contrasts with the current dynamics underlying oscillations in thalamocortical and reticularis neurons, which primarily involveI T andI H , orI T andI A H P respectively. These findings thus point to a novel electrophysiological mechanism for generating intrinsic oscillations in a major thalamic cell type. Because local interneurons can sculpt the behavior of thalamocortical circuits, these results suggest new targets for the manipulation of ascending thalamocortical network activity.

In-silico testing of new pharmacology for restoring inhibition and human cortical function in depression

Reduced inhibition by somatostatin-expressing interneurons is associated with depression. Administration of positive allosteric modulators of α5 subunit-containing GABAA receptor (α5-PAM) that selectively target this lost inhibition exhibit antidepressant and pro-cognitive effects in rodent models of chronic stress. However, the functional effects of α5-PAM on the human brain in vivo are unknown, and currently cannot be assessed experimentally. We modeled the effects of α5-PAM on tonic inhibition as measured in human neurons, and tested in silico α5-PAM effects on detailed models of human cortical microcircuits in health and depression. We found that α5-PAM effectively recovered impaired cortical processing as quantified by stimulus detection metrics, and also recovered the power spectral density profile of the microcircuit EEG signals. We performed an α5-PAM dose-response and identified simulated EEG biomarker candidates. Our results serve to de-risk and facilitate α5-PAM translation and provide biomarkers in non-invasive brain signals for monitoring target engagement and drug efficacy.

Detecting Spontaneous Neural Oscillation Events in Primate Auditory Cortex

Electrophysiological oscillations in the brain have been shown to occur as multicycle events, with onset and offset dependent on behavioral and cognitive state. To provide a baseline for state-related and task-related events, we quantified oscillation features in resting-state recordings. We developed an open-source wavelet-based tool to detect and characterize such oscillation events (OEvents) and exemplify the use of this tool in both simulations and two invasively-recorded electrophysiology datasets: one from human, and one from nonhuman primate (NHP) auditory system. After removing incidentally occurring event-related potentials (ERPs), we used OEvents to quantify oscillation features. We identified ∼2 million oscillation events, classified within traditional frequency bands: δ, θ, α, β, low γ, γ, and high γ. Oscillation events of 1-44 cycles could be identified in at least one frequency band 90% of the time in human and NHP recordings. Individual oscillation events were characterized by nonconstant frequency and amplitude. This result necessarily contrasts with prior studies which assumed frequency constancy, but is consistent with evidence from event-associated oscillations. We measured oscillation event duration, frequency span, and waveform shape. Oscillations tended to exhibit multiple cycles per event, verifiable by comparing filtered to unfiltered waveforms. In addition to the clear intraevent rhythmicity, there was also evidence of interevent rhythmicity within bands, demonstrated by finding that coefficient of variation of interval distributions and Fano factor (FF) measures differed significantly from a Poisson distribution assumption. Overall, our study provides an easy-to-use tool to study oscillation events at the single-trial level or in ongoing recordings, and demonstrates that rhythmic, multicycle oscillation events dominate auditory cortical dynamics.