On the Validity of Neural Mass Models.
Ontology highlight
ABSTRACT:
SUBMITTER: Deschle N
PROVIDER: S-EPMC7814001 | biostudies-literature | 2020
REPOSITORIES: biostudies-literature
Publications
On the Validity of Neural Mass Models.
Deschle Nicolás N Ignacio Gossn Juan J Tewarie Prejaas P Schelter Björn B Daffertshofer Andreas A
Frontiers in computational neuroscience 20210105
Modeling the dynamics of neural masses is a common approach in the study of neural populations. Various models have been proven useful to describe a plenitude of empirical observations including self-sustained local oscillations and patterns of distant synchronization. We discuss the extent to which mass models really resemble the mean dynamics of a neural population. In particular, we question the validity of neural mass models if the population under study comprises a mixture of excitatory and ...[more]
Similar Datasets
Project description:Neural Mass Models provide a compact description of the dynamical activity of cell populations in neocortical regions. Moreover, models of regional activity can be connected together into networks, and inferences made about the strength of connections, using M/EEG data and Bayesian inference. To date, however, Bayesian methods have been largely restricted to the Variational Laplace (VL) algorithm which assumes that the posterior distribution is Gaussian and finds model parameters that are only locally optimal. This paper explores the use of Annealed Importance Sampling (AIS) to address these restrictions. We implement AIS using proposals derived from Langevin Monte Carlo (LMC) which uses local gradient and curvature information for efficient exploration of parameter space. In terms of the estimation of Bayes factors, VL and AIS agree about which model is best but report different degrees of belief. Additionally, AIS finds better model parameters and we find evidence of non-Gaussianity in their posterior distribution.
| S-EPMC4778905 | biostudies-literature
Project description:Theta-nested gamma oscillations have been reported in many areas of the brain and are believed to represent a fundamental mechanism to transfer information across spatial and temporal scales. In a series of recent experiments in vitro it has been possible to replicate with an optogenetic theta frequency stimulation several features of cross-frequency coupling (CFC) among theta and gamma rhythms observed in behaving animals. In order to reproduce the main findings of these experiments we have considered a new class of neural mass models able to reproduce exactly the macroscopic dynamics of spiking neural networks. In this framework, we have examined two set-ups able to support collective gamma oscillations: namely, the pyramidal interneuronal network gamma (PING) and the interneuronal network gamma (ING). In both set-ups we observe the emergence of theta-nested gamma oscillations by driving the system with a sinusoidal theta-forcing in proximity of a Hopf bifurcation. These mixed rhythms always display phase amplitude coupling. However, two different types of nested oscillations can be identified: one characterized by a perfect phase locking between theta and gamma rhythms, corresponding to an overall periodic behavior; another one where the locking is imperfect and the dynamics is quasi-periodic or even chaotic. From our analysis it emerges that the locked states are more frequent in the ING set-up. In agreement with the experiments, we find theta-nested gamma oscillations for forcing frequencies in the range [1:10] Hz, whose amplitudes grow proportionally to the forcing intensity and which are clearly modulated by the theta phase. Furthermore, analogously to the experiments, the gamma power and the frequency of the gamma-power peak increase with the forcing amplitude. At variance with experimental findings, the gamma-power peak does not shift to higher frequencies by increasing the theta frequency. This effect can be obtained, in our model, only by incrementing, at the same time, also the stimulation power. An effect achieved by increasing the amplitude either of the noise or of the forcing term proportionally to the theta frequency. On the basis of our analysis both the PING and the ING mechanism give rise to theta-nested gamma oscillations with almost identical features.
| S-EPMC7270590 | biostudies-literature
Project description:New brain atlases with high spatial resolution and whole-brain coverage have rapidly advanced our knowledge of the brain's neural architecture, including the systematic variation of excitatory and inhibitory cell densities across the mammalian cortex. But understanding how the brain's microscale physiology shapes brain dynamics at the macroscale has remained a challenge. While physiologically based mathematical models of brain dynamics are well placed to bridge this explanatory gap, their complexity can form a barrier to providing clear mechanistic interpretation of the dynamics they generate. In this work, we develop a neural-mass model of the mouse cortex and show how bifurcation diagrams, which capture local dynamical responses to inputs and their variation across brain regions, can be used to understand the resulting whole-brain dynamics. We show that strong fits to resting-state functional magnetic resonance imaging (fMRI) data can be found in surprisingly simple dynamical regimes-including where all brain regions are confined to a stable fixed point-in which regions are able to respond strongly to variations in their inputs, consistent with direct structural connections providing a strong constraint on functional connectivity in the anesthetized mouse. We also use bifurcation diagrams to show how perturbations to local excitatory and inhibitory coupling strengths across the cortex, constrained by cell-density data, provide spatially dependent constraints on resulting cortical activity, and support a greater diversity of coincident dynamical regimes. Our work illustrates methods for visualizing and interpreting model performance in terms of underlying dynamical mechanisms, an approach that is crucial for building explanatory and physiologically grounded models of the dynamical principles that underpin large-scale brain activity.
| S-EPMC9081874 | biostudies-literature
Project description:Neural mass models (NMMs) are important for helping us interpret observations of brain dynamics. They provide a means to understand data in terms of mechanisms such as synaptic interactions between excitatory and inhibitory neuronal populations. To interpret data using NMMs we need to quantitatively compare the output of NMMs with data, and thereby find parameter values for which the model can produce the observed dynamics. Mapping dynamics to NMM parameter values in this way has the potential to improve our understanding of the brain in health and disease. Though abstract, NMMs still comprise of many parameters that are difficult to constrain a priori. This makes it challenging to explore the dynamics of NMMs and elucidate regions of parameter space in which their dynamics best approximate data. Existing approaches to overcome this challenge use a combination of linearising models, constraining the values they can take and exploring restricted subspaces by fixing the values of many parameters a priori. As such, we have little knowledge of the extent to which different regions of parameter space of NMMs can yield dynamics that approximate data, how nonlinearities in models can affect parameter mapping or how best to quantify similarities between model output and data. These issues need to be addressed in order to fully understand the potential and limitations of NMMs, and to aid the development of new models of brain dynamics in the future. To begin to overcome these issues, we present a global nonlinear approach to recovering parameters of NMMs from data. We use global optimisation to explore all parameters of nonlinear NMMs simultaneously, in a minimally constrained way. We do this using multi-objective optimisation (multi-objective evolutionary algorithm, MOEA) so that multiple data features can be quantified. In particular, we use the weighted horizontal visibility graph (wHVG), which is a flexible framework for quantifying different aspects of time series, by converting them into networks. We study EEG alpha activity recorded during the eyes closed resting state from 20 healthy individuals and demonstrate that the MOEA performs favourably compared to single objective approaches. The addition of the wHVG objective allows us to better constrain the model output, which leads to the recovered parameter values being restricted to smaller regions of parameter space, thus improving the practical identifiability of the model. We then use the MOEA to study differences in the alpha rhythm observed in EEG recorded from 20 people with epilepsy. We find that a small number of parameters can explain this difference and that, counterintuitively, the mean excitatory synaptic gain parameter is reduced in people with epilepsy compared to control. In addition, we propose that the MOEA could be used to mine for the presence of pathological rhythms, and demonstrate the application of this to epileptiform spike-wave discharges.
| S-EPMC10075456 | biostudies-literature
Project description:Neural mass models (NMMs) are increasingly used to uncover the large-scale mechanisms of brain rhythms in health and disease. The dynamics of these models is dependent upon the choice of parameters, and therefore it is crucial to be able to understand how dynamics change when parameters are varied. Despite being considered low dimensional in comparison to micro-scale, neuronal network models, with regards to understanding the relationship between parameters and dynamics, NMMs are still prohibitively high dimensional for classical approaches such as numerical continuation. Therefore, we need alternative methods to characterise dynamics of NMMs in high dimensional parameter spaces. Here, we introduce a statistical framework that enables the efficient exploration of the relationship between model parameters and selected features of the simulated, emergent model dynamics of NMMs. We combine the classical machine learning approaches of trees and random forests to enable studying the effect that varying multiple parameters has on the dynamics of a model. The method proceeds by using simulations to transform the mathematical model into a database. This database is then used to partition parameter space with respect to dynamic features of interest, using random forests. This allows us to rapidly explore dynamics in high dimensional parameter space, capture the approximate location of qualitative transitions in dynamics and assess the relative importance of all parameters in the model in all dimensions simultaneously. We apply this method to a commonly used NMM in the context of transitions to seizure dynamics. We find that the inhibitory sub-system is most crucial for the generation of seizure dynamics, confirm and expand previous findings regarding the ratio of excitation and inhibition, and demonstrate that previously overlooked parameters can have a significant impact on model dynamics. We advocate the use of this method in future to constrain high dimensional parameter spaces enabling more efficient, person-specific, model calibration.
| S-EPMC5851637 | biostudies-literature
Project description:Brain rhythms emerge from the mean-field activity of networks of neurons. There have been many efforts to build mathematical and computational embodiments in the form of discrete cell-group activities-termed neural masses-to understand in particular the origins of evoked potentials, intrinsic patterns of activities such as theta, regulation of sleep, Parkinson's disease related dynamics, and mimic seizure dynamics. As originally utilized, standard neural masses convert input through a sigmoidal function to a firing rate, and firing rate through a synaptic alpha function to other masses. Here we define a process to build mechanistic neural masses (mNMs) as mean-field models of microscopic membrane-type (Hodgkin Huxley type) models of different neuron types that duplicate the stability, firing rate, and associated bifurcations as function of relevant slow variables - such as extracellular potassium - and synaptic current; and whose output is both firing rate and impact on the slow variables - such as transmembrane potassium flux. Small networks composed of just excitatory and inhibitory mNMs demonstrate expected dynamical states including firing, runaway excitation and depolarization block, and these transitions change in biologically observed ways with changes in extracellular potassium and excitatory-inhibitory balance.
| S-EPMC9980379 | biostudies-literature
Project description:The extraction of electrophysiological features that reliably forecast the occurrence of seizures is one of the most challenging goals in epilepsy research. Among possible approaches to tackle this problem is the use of active probing paradigms in which responses to stimuli are used to detect underlying system changes leading up to seizures. This work evaluates the theoretical and mechanistic underpinnings of this strategy using two coupled populations of the well-studied Wendling neural mass model. Different model settings are evaluated, shifting parameters (excitability, slow inhibition, or inter-population coupling gains) from normal towards ictal states while probing stimuli are applied every 2 seconds to the input of either one or both populations. The correlation between the extracted features and the ictogenic parameter shifting indicates if the impending transition to the ictal state may be identified in advance. Results show that not only can the response to the probing stimuli forecast seizures but this is true regardless of the altered ictogenic parameter. That is, similar feature changes are highlighted by probing stimuli responses in advance of the seizure including: increased response variance and lag-1 autocorrelation, decreased skewness, and increased mutual information between the outputs of both model subsets. These changes were mostly restricted to the stimulated population, showing a local effect of this perturbational approach. The transition latencies from normal activity to sustained discharges of spikes were not affected, suggesting that stimuli had no pro-ictal effects. However, stimuli were found to elicit interictal-like spikes just before the transition to the ictal state. Furthermore, the observed feature changes highlighted by probing the neuronal populations may reflect the phenomenon of critical slowing down, where increased recovery times from perturbations may signal the loss of a systems' resilience and are common hallmarks of an impending critical transition. These results provide more evidence that active probing approaches highlight information about underlying system changes involved in ictogenesis and may be able to play a role in assisting seizure forecasting methods which can be incorporated into early-warning systems that ultimately enable closing the loop for targeted seizure-controlling interventions.
| S-EPMC7861539 | biostudies-literature
Project description:Recently, in studies examining fibroblasts obtained from the tissues of one set of monozygotic twins (i.e. fetuses derived from the same egg) discordant for trisomy 21 (Down syndrome; DS), Letourneau et al. reported the presence of a defined pattern of dysregulation within specific genomic domains they referred to as Gene Expression Dysregulated Domains (GEDDs). GEDDs were described as alternating segments of increased or decreased gene expression affecting all chromosomes. Strikingly, GEDDs in fibroblasts were largely conserved in induced pluripotent cells (iPSCs) generated from the twin’s fibroblasts as well as in fibroblasts from the Ts65Dn mouse model of DS. Our recent analysis failed to find GEDDs. We reexamined human iPSCs RNAseq data from Letourneau et al., and data from this same research group published earlier examining iPSCs from the same monozygotic twins. In addition, an independent analysis of RNAseq data from Ts65Dn fibroblasts also failed to confirm presence of GEDDs. Our analysis questions the validity of GEDDs in DS.
2015-02-21 | GSE64840 | GEO
Project description:Objective: Assessing brain connectivity from electrophysiological signals is of great relevance in neuroscience, but results are still debated and depend crucially on how connectivity is defined and on mathematical instruments utilized. Aim of this work is to assess the capacity of bivariate Transfer Entropy (TE) to evaluate connectivity, using data generated from simple neural mass models of connected Regions of Interest (ROIs). Approach: Signals simulating mean field potentials were generated assuming two, three or four ROIs, connected via excitatory or by-synaptic inhibitory links. We investigated whether the presence of a statistically significant connection can be detected and if connection strength can be quantified. Main Results: Results suggest that TE can reliably estimate the strength of connectivity if neural populations work in their linear regions, and if the epoch lengths are longer than 10 s. In case of multivariate networks, some spurious connections can emerge (i.e., a statistically significant TE even in the absence of a true connection); however, quite a good correlation between TE and synaptic strength is still preserved. Moreover, TE appears more robust for distal regions (longer delays) compared with proximal regions (smaller delays): an approximate a priori knowledge on this delay can improve the procedure. Finally, non-linear phenomena affect the assessment of connectivity, since they may significantly reduce TE estimation: information transmission between two ROIs may be weak, due to non-linear phenomena, even if a strong causal connection is present. Significance: Changes in functional connectivity during different tasks or brain conditions, might not always reflect a true change in the connecting network, but rather a change in information transmission. A limitation of the work is the use of bivariate TE. In perspective, the use of multivariate TE can improve estimation and reduce some of the problems encountered in the present study.
| S-EPMC7292208 | biostudies-literature
Project description:BackgroundDeleterious mutations of the BRCA1 and BRCA2 genes confer susceptibility to breast and ovarian cancer. At least 7 models for estimating the probabilities of having a mutation are used widely in clinical and scientific activities; however, the merits and limitations of these models are not fully understood.ObjectiveTo systematically quantify the accuracy of the following publicly available models to predict mutation carrier status: BRCAPRO, family history assessment tool, Finnish, Myriad, National Cancer Institute, University of Pennsylvania, and Yale University.DesignCross-sectional validation study, using model predictions and BRCA1 or BRCA2 mutation status of patients different from those used to develop the models.SettingMulticenter study across Cancer Genetics Network participating centers.Patients3 population-based samples of participants in research studies and 8 samples from genetic counseling clinics.MeasurementsDiscrimination between individuals testing positive for a mutation in BRCA1 or BRCA2 from those testing negative, as measured by the c-statistic, and sensitivity and specificity of model predictions.ResultsThe 7 models differ in their predictions. The better-performing models have a c-statistic around 80%. BRCAPRO has the largest c-statistic overall and in all but 2 patient subgroups, although the margin over other models is narrow in many strata. Outside of high-risk populations, all models have high false-negative and false-positive rates across a range of probability thresholds used to refer for mutation testing.LimitationThree recently published models were not included.ConclusionsAll models identify women who probably carry a deleterious mutation of BRCA1 or BRCA2 with adequate discrimination to support individualized genetic counseling, although discrimination varies across models and populations.
| S-EPMC2423214 | biostudies-literature