Global nonlinear approach for mapping parameters of neural mass models.

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.

Using deep clustering to improve fMRI dynamic functional connectivity analysis.

Dynamic functional connectivity (dFC) analysis of resting-state fMRI data is commonly performed by calculating sliding-window correlations (SWC), followed by k-means clustering in order to assign each window to a given state. Studies using synthetic data have shown that k-means performance is highly dependent on sliding window parameters and signal-to-noise ratio. Additionally, sources of heterogeneity between subjects may affect the accuracy of group-level clustering, thus affecting measurements of dFC state temporal properties such as dwell time and fractional occupancy. This may result in spurious conclusions regarding differences between groups (e.g. when comparing a clinical population to healthy controls). Therefore, is it important to quantify the ability of k-means to estimate dFC state temporal properties when applied to cohorts of multiple subjects, and to explore ways in which clustering performance can be maximised. Here, we explore the use of dimensionality reduction methods prior to clustering in order to map high-dimensional data to a lower dimensional space, providing salient features to the subsequent clustering step. We assess the use of deep autoencoders for dimensionality reduction prior to applying k-means clustering to the encoded data. We compare this deep clustering method to dimensionality reduction using principle component analysis (PCA), uniform manifold approximation and projection (UMAP), as well as applying k-means to the original feature space using either L1 or L2 distance. We provide extensive quantitative evaluation of clustering performance using synthetic datasets, representing data from multiple heterogeneous subjects. In synthetic data we find that deep clustering gives the best performance, while other approaches are often insufficient to capture temporal properties of dFC states. We then demonstrate the application of each method to real-world data from human subjects and show that the choice of dimensionality reduction method has a significant effect on group-level measurements of state temporal properties.

A large-scale brain network mechanism for increased seizure propensity in Alzheimer's disease.

People with Alzheimer's disease (AD) are 6-10 times more likely to develop seizures than the healthy aging population. Leading hypotheses largely consider hyperexcitability of local cortical tissue as primarily responsible for increased seizure prevalence in AD. However, in the general population of people with epilepsy, large-scale brain network organization additionally plays a role in determining seizure likelihood and phenotype. Here, we propose that alterations to large-scale brain network organization seen in AD may contribute to increased seizure likelihood. To test this hypothesis, we combine computational modelling with electrophysiological data using an approach that has proved informative in clinical epilepsy cohorts without AD. EEG was recorded from 21 people with probable AD and 26 healthy controls. At the time of EEG acquisition, all participants were free from seizures. Whole brain functional connectivity derived from source-reconstructed EEG recordings was used to build subject-specific brain network models of seizure transitions. As cortical tissue excitability was increased in the simulations, AD simulations were more likely to transition into seizures than simulations from healthy controls, suggesting an increased group-level probability of developing seizures at a future time for AD participants. We subsequently used the model to assess seizure propensity of different regions across the cortex. We found the most important regions for seizure generation were those typically burdened by amyloid-beta at the early stages of AD, as previously reported by in-vivo and post-mortem staging of amyloid plaques. Analysis of these spatial distributions also give potential insight into mechanisms of increased susceptibility to generalized (as opposed to focal) seizures in AD vs controls. This research suggests avenues for future studies testing patients with seizures, e.g. co-morbid AD/epilepsy patients, and comparisons with PET and MRI scans to relate regional seizure propensity with AD pathologies.

Classifying dynamic transitions in high dimensional neural mass models: A random forest approach.

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.

Control of clustered action potential firing in a mathematical model of entorhinal cortex stellate cells.

The entorhinal cortex is a crucial component of our memory and spatial navigation systems and is one of the first areas to be affected in dementias featuring tau pathology, such as Alzheimer's disease and frontotemporal dementia. Electrophysiological recordings from principle cells of medial entorhinal cortex (layer II stellate cells, mEC-SCs) demonstrate a number of key identifying properties including subthreshold oscillations in the theta (4-12 Hz) range and clustered action potential firing. These single cell properties are correlated with network activity such as grid firing and coupling between theta and gamma rhythms, suggesting they are important for spatial memory. As such, experimental models of dementia have revealed disruption of organised dorsoventral gradients in clustered action potential firing. To better understand the mechanisms underpinning these different dynamics, we study a conductance based model of mEC-SCs. We demonstrate that the model, driven by extrinsic noise, can capture quantitative differences in clustered action potential firing patterns recorded from experimental models of tau pathology and healthy animals. The differential equation formulation of our model allows us to perform numerical bifurcation analyses in order to uncover the dynamic mechanisms underlying these patterns. We show that clustered dynamics can be understood as subcritical Hopf/homoclinic bursting in a fast-slow system where the slow sub-system is governed by activation of the persistent sodium current and inactivation of the slow A-type potassium current. In the full system, we demonstrate that clustered firing arises via flip bifurcations as conductance parameters are varied. Our model analyses confirm the experimentally suggested hypothesis that the breakdown of clustered dynamics in disease occurs via increases in AHP conductance.

An optimal strategy for epilepsy surgery: Disruption of the rich-club?

Surgery is a therapeutic option for people with epilepsy whose seizures are not controlled by anti-epilepsy drugs. In pre-surgical planning, an array of data modalities, often including intra-cranial EEG, is used in an attempt to map regions of the brain thought to be crucial for the generation of seizures. These regions are then resected with the hope that the individual is rendered seizure free as a consequence. However, post-operative seizure freedom is currently sub-optimal, suggesting that the pre-surgical assessment may be improved by taking advantage of a mechanistic understanding of seizure generation in large brain networks. Herein we use mathematical models to uncover the relative contribution of regions of the brain to seizure generation and consequently which brain regions should be considered for resection. A critical advantage of this modeling approach is that the effect of different surgical strategies can be predicted and quantitatively compared in advance of surgery. Herein we seek to understand seizure generation in networks with different topologies and study how the removal of different nodes in these networks reduces the occurrence of seizures. Since this a computationally demanding problem, a first step for this aim is to facilitate tractability of this approach for large networks. To do this, we demonstrate that predictions arising from a neural mass model are preserved in a lower dimensional, canonical model that is quicker to simulate. We then use this simpler model to study the emergence of seizures in artificial networks with different topologies, and calculate which nodes should be removed to render the network seizure free. We find that for scale-free and rich-club networks there exist specific nodes that are critical for seizure generation and should therefore be removed, whereas for small-world networks the strategy should instead focus on removing sufficient brain tissue. We demonstrate the validity of our approach by analysing intra-cranial EEG recordings from a database comprising 16 patients who have undergone epilepsy surgery, revealing rich-club structures within the obtained functional networks. We show that the postsurgical outcome for these patients was better when a greater proportion of the rich club was removed, in agreement with our theoretical predictions.

Electrical and Network Neuronal Properties Are Preferentially Disrupted in Dorsal, But Not Ventral, Medial Entorhinal Cortex in a Mouse Model of Tauopathy.

The entorhinal cortex (EC) is one of the first areas to be disrupted in neurodegenerative diseases such as Alzheimer's disease and frontotemporal dementia. The responsiveness of individual neurons to electrical and environmental stimuli varies along the dorsal-ventral axis of the medial EC (mEC) in a manner that suggests this topographical organization plays a key role in neural encoding of geometric space. We examined the cellular properties of layer II mEC stellate neurons (mEC-SCs) in rTg4510 mice, a rodent model of neurodegeneration. Dorsoventral gradients in certain intrinsic membrane properties, such as membrane capacitance and afterhyperpolarizations, were flattened in rTg4510 mEC-SCs, while other cellular gradients [e.g., input resistance (Ri), action potential properties] remained intact. Specifically, the intrinsic properties of rTg4510 mEC-SCs in dorsal aspects of the mEC were preferentially affected, such that action potential firing patterns in dorsal mEC-SCs were altered, while those in ventral mEC-SCs were unaffected. We also found that neuronal oscillations in the gamma frequency band (30-80 Hz) were preferentially disrupted in the dorsal mEC of rTg4510 slices, while those in ventral regions were comparatively preserved. These alterations corresponded to a flattened dorsoventral gradient in theta-gamma cross-frequency coupling of local field potentials recorded from the mEC of freely moving rTg4510 mice. These differences were not paralleled by changes to the dorsoventral gradient in parvalbumin staining or neurodegeneration. We propose that the selective disruption to dorsal mECs, and the resultant flattening of certain dorsoventral gradients, may contribute to disturbances in spatial information processing observed in this model of dementia. SIGNIFICANCE STATEMENT: The medial entorhinal cortex (mEC) plays a key role in spatial memory and is one of the first areas to express the pathological features of dementia. Neurons of the mEC are anatomically arranged to express functional dorsoventral gradients in a variety of neuronal properties, including grid cell firing field spacing, which is thought to encode geometric scale. We have investigated the effects of tau pathology on functional dorsoventral gradients in the mEC. Using electrophysiological approaches, we have shown that, in a transgenic mouse model of dementia, the functional properties of the dorsal mEC are preferentially disrupted, resulting in a flattening of some dorsoventral gradients. Our data suggest that neural signals arising in the mEC will have a reduced spatial content in dementia.

Temporally correlated fluctuations drive epileptiform dynamics.

Macroscopic models of brain networks typically incorporate assumptions regarding the characteristics of afferent noise, which is used to represent input from distal brain regions or ongoing fluctuations in non-modelled parts of the brain. Such inputs are often modelled by Gaussian white noise which has a flat power spectrum. In contrast, macroscopic fluctuations in the brain typically follow a 1/fb spectrum. It is therefore important to understand the effect on brain dynamics of deviations from the assumption of white noise. In particular, we wish to understand the role that noise might play in eliciting aberrant rhythms in the epileptic brain. To address this question we study the response of a neural mass model to driving by stochastic, temporally correlated input. We characterise the model in terms of whether it generates "healthy" or "epileptiform" dynamics and observe which of these dynamics predominate under different choices of temporal correlation and amplitude of an Ornstein-Uhlenbeck process. We find that certain temporal correlations are prone to eliciting epileptiform dynamics, and that these correlations produce noise with maximal power in the δ and θ bands. Crucially, these are rhythms that are found to be enhanced prior to seizures in humans and animal models of epilepsy. In order to understand why these rhythms can generate epileptiform dynamics, we analyse the response of the model to sinusoidal driving and explain how the bifurcation structure of the model gives rise to these findings. Our results provide insight into how ongoing fluctuations in brain dynamics can facilitate the onset and propagation of epileptiform rhythms in brain networks. Furthermore, we highlight the need to combine large-scale models with noise of a variety of different types in order to understand brain (dys-)function.

A computational biomarker of idiopathic generalized epilepsy from resting state EEG.

Epilepsy is one of the most common serious neurologic conditions. It is characterized by the tendency to have recurrent seizures, which arise against a backdrop of apparently normal brain activity. At present, clinical diagnosis relies on the following: (1) case history, which can be unreliable; (2) observation of transient abnormal activity during electroencephalography (EEG), which may not be present during clinical evaluation; and (3) if diagnostic uncertainty occurs, undertaking prolonged monitoring in an attempt to observe EEG abnormalities, which is costly. Herein, we describe the discovery and validation of an epilepsy biomarker based on computational analysis of a short segment of resting-state (interictal) EEG. Our method utilizes a computer model of dynamic networks, where the network is inferred from the extent of synchrony between EEG channels (functional networks) and the normalized power spectrum of the clinical data. We optimize model parameters using a leave-one-out classification on a dataset comprising 30 people with idiopathic generalized epilepsy (IGE) and 38 normal controls. Applying this scheme to all 68 subjects we find 100% specificity at 56.7% sensitivity, and 100% sensitivity at 65.8% specificity. We believe this biomarker could readily provide additional support to the diagnostic process.

Dynamic mechanisms of neocortical focal seizure onset.

Recent experimental and clinical studies have provided diverse insight into the mechanisms of human focal seizure initiation and propagation. Often these findings exist at different scales of observation, and are not reconciled into a common understanding. Here we develop a new, multiscale mathematical model of cortical electric activity with realistic mesoscopic connectivity. Relating the model dynamics to experimental and clinical findings leads us to propose three classes of dynamical mechanisms for the onset of focal seizures in a unified framework. These three classes are: (i) globally induced focal seizures; (ii) globally supported focal seizures; (iii) locally induced focal seizures. Using model simulations we illustrate these onset mechanisms and show how the three classes can be distinguished. Specifically, we find that although all focal seizures typically appear to arise from localised tissue, the mechanisms of onset could be due to either localised processes or processes on a larger spatial scale. We conclude that although focal seizures might have different patient-specific aetiologies and electrographic signatures, our model suggests that dynamically they can still be classified in a clinically useful way. Additionally, this novel classification according to the dynamical mechanisms is able to resolve some of the previously conflicting experimental and clinical findings.

Resting-state functional connectivity in children cooled for neonatal encephalopathy.

Therapeutic hypothermia improves outcomes following neonatal hypoxic-ischaemic encephalopathy, reducing cases of death and severe disability such as cerebral palsy compared with normothermia management. However, when cooled children reach early school-age, they have cognitive and motor impairments which are associated with underlying alterations to brain structure and white matter connectivity. It is unknown whether these differences in structural connectivity are associated with differences in functional connectivity between cooled children and healthy controls. Resting-state functional MRI has been used to characterize static and dynamic functional connectivity in children, both with typical development and those with neurodevelopmental disorders. Previous studies of resting-state brain networks in children with hypoxic-ischaemic encephalopathy have focussed on the neonatal period. In this study, we used resting-state fMRI to investigate static and dynamic functional connectivity in children aged 6-8 years who were cooled for neonatal hypoxic-ischaemic without cerebral palsy [n = 22, median age (interquartile range) 7.08 (6.85-7.52) years] and healthy controls matched for age, sex and socioeconomic status [n = 20, median age (interquartile range) 6.75 (6.48-7.25) years]. Using group independent component analysis, we identified 31 intrinsic functional connectivity networks consistent with those previously reported in children and adults. We found no case-control differences in the spatial maps of these intrinsic connectivity networks. We constructed subject-specific static functional connectivity networks by measuring pairwise Pearson correlations between component time courses and found no case-control differences in functional connectivity after false discovery rate correction. To study the time-varying organization of resting-state networks, we used sliding window correlations and deep clustering to investigate dynamic functional connectivity characteristics. We found k = 4 repetitively occurring functional connectivity states, which exhibited no case-control differences in dwell time, fractional occupancy or state functional connectivity matrices. In this small cohort, the spatiotemporal characteristics of resting-state brain networks in cooled children without severe disability were too subtle to be differentiated from healthy controls at early school-age, despite underlying differences in brain structure and white matter connectivity, possibly reflecting a level of recovery of healthy resting-state brain function. To our knowledge, this is the first study to investigate resting-state functional connectivity in children with hypoxic-ischaemic encephalopathy beyond the neonatal period and the first to investigate dynamic functional connectivity in any children with hypoxic-ischaemic encephalopathy.

Towards a large-scale model of patient-specific epileptic spike-wave discharges.

Clinical electroencephalographic (EEG) recordings of the transition into generalised epileptic seizures show a sudden onset of spike-wave dynamics from a low-amplitude irregular background. In addition, non-trivial and variable spatio-temporal dynamics are widely reported in combined EEG/fMRI studies on the scale of the whole cortex. It is unknown whether these characteristics can be accounted for in a large-scale mathematical model with fixed heterogeneous long-range connectivities. Here, we develop a modelling framework with which to investigate such EEG features. We show that a neural field model composed of a few coupled compartments can serve as a low-dimensional prototype for the transition between irregular background dynamics and spike-wave activity. This prototype then serves as a node in a large-scale network with long-range connectivities derived from human diffusion-tensor imaging data. We examine multivariate properties in 42 clinical EEG seizure recordings from 10 patients diagnosed with typical absence epilepsy and 50 simulated seizures from the large-scale model using 10 DTI connectivity sets from humans. The model can reproduce the clinical feature of stereotypy where seizures are more similar within a patient than between patients, essentially creating a patient-specific fingerprint. We propose the approach as a feasible technique for the investigation of patient-specific large-scale epileptic features in space and time.

Self-organised transients in a neural mass model of epileptogenic tissue dynamics.

Stimulation of human epileptic tissue can induce rhythmic, self-terminating responses on the EEG or ECoG. These responses play a potentially important role in localising tissue involved in the generation of seizure activity, yet the underlying mechanisms are unknown. However, in vitro evidence suggests that self-terminating oscillations in nervous tissue are underpinned by non-trivial spatio-temporal dynamics in an excitable medium. In this study, we investigate this hypothesis in spatial extensions to a neural mass model for epileptiform dynamics. We demonstrate that spatial extensions to this model in one and two dimensions display propagating travelling waves but also more complex transient dynamics in response to local perturbations. The neural mass formulation with local excitatory and inhibitory circuits, allows the direct incorporation of spatially distributed, functional heterogeneities into the model. We show that such heterogeneities can lead to prolonged reverberating responses to a single pulse perturbation, depending upon the location at which the stimulus is delivered. This leads to the hypothesis that prolonged rhythmic responses to local stimulation in epileptogenic tissue result from repeated self-excitation of regions of tissue with diminished inhibitory capabilities. Combined with previous models of the dynamics of focal seizures this macroscopic framework is a first step towards an explicit spatial formulation of the concept of the epileptogenic zone. Ultimately, an improved understanding of the pathophysiologic mechanisms of the epileptogenic zone will help to improve diagnostic and therapeutic measures for treating epilepsy.

Modelling the role of tissue heterogeneity in epileptic rhythms.

Epileptic seizure activity manifests as complex spatio-temporal dynamics on the clinically relevant macroscopic scale. These dynamics are known to arise from spatially heterogeneous tissue, but the relationship between specific spatial abnormalities and epileptic rhythm generation is not well understood. We formulate a simplified macroscopic modelling framework with which to study the role of spatial heterogeneity in the generation of epileptiform spatio-temporal rhythms. We characterize the overall model dynamics in terms of spontaneous activity and excitability and demonstrate normal and abnormal spreading of activity. We introduce a means to systematically investigate the topology of abnormal sub-networks and explore its impact on spontaneous and stimulus-evoked rhythmic dynamics. This computationally efficient framework complements results from detailed biophysical models, and allows the testing of specific hypotheses about epileptic dynamics on the macroscopic scale.

On seeing the trees and the forest: single-signal and multisignal analysis of periictal intracranial EEG.

PURPOSE: Epileptic seizures are associated with a dysregulation of electrical brain activity on many different spatial scales. To better understand the dynamics of epileptic seizures, that is, how the seizures initiate, propagate, and terminate, it is important to consider changes of electrical brain activity on different spatial scales. Herein we set out to analyze periictal electrical brain activity on comparatively small and large spatial scales by assessing changes in single intracranial electroencephalography (EEG) signals and of averaged interdependences of pairs of EEG signals. METHODS: Single and multiple EEG signals are analyzed by combining methods from symbolic dynamics and information theory. This computationally efficient approach is chosen because at its core it consists of analyzing the occurrence of patterns and bears analogy to classical visual EEG reading. Symbolization is achieved by first mapping the EEG signals into bit strings, that is, long sequences of zeros and ones, depending solely on whether their amplitudes increase or decrease. Bit strings reflect relational aspects between consecutive values of the original EEG signals, but not the values themselves. For each bit string the relative frequencies of the different constituent short bit patterns are then determined and used to compute two information theoretical measures: (1) redundancy (R) of single bit strings characterizes electrical brain activity on a comparatively small spatial scale represented by a single EEG signal and (2) averaged pair-wise mutual information with all other bit strings (M), which allows tracking of larger-scale EEG dynamics. KEY FINDINGS: We analyzed 20 periictal intracranial EEG recordings from five patients with pharmacoresistant temporal lobe epilepsy. At seizure onset, R first strongly increased and then decreased toward seizure termination, whereas M gradually increased throughout the seizure. Bit strings with maximal R were always derived from EEG signals recorded from the visually identified seizure-onset zone. When compared to the bit strings derived from other EEG signals, their M was relatively smaller. These findings are consistent with a strong but transient occurrence of information-poor, that is, redundant electrical brain activity on a smaller spatial scale, which is particularly pronounced in the seizure-onset zone. On a larger spatial scale, a progressively more collective state emerges, as revealed by increasing amounts of mutual information. SIGNIFICANCE: Information theoretical analysis of bit patterns derived from EEG signals helps to characterize periictal brain activity on different spatial scales in a quantitative and efficient way and may provide clinically relevant results.

A prototypic mathematical model of the human hair cycle.

The human hair cycle is a complex, dynamic organ-transformation process during which the hair follicle repetitively progresses from a growth phase (anagen) to a rapid apoptosis-driven involution (catagen) and finally a relative quiescent phase (telogen) before returning to anagen. At present no theory satisfactorily explains the origin of the hair cycle rhythm. Based on experimental evidence we propose a prototypic model that focuses on the dynamics of hair matrix keratinocytes. We argue that a plausible feedback-control structure between two key compartments (matrix keratinocytes and dermal papilla) leads to dynamic instabilities in the population dynamics resulting in rhythmic hair growth. The underlying oscillation consists of an autonomous switching between two quasi-steady states. Additional features of the model, namely bistability and excitability, lead to new hypotheses about the impact of interventions on hair growth. We show how in silico testing may facilitate testing of candidate hair growth modulatory agents in human HF organ culture or in clinical trials.

What Models and Tools can Contribute to a Better Understanding of Brain Activity?

Despite impressive scientific advances in understanding the structure and function of the human brain, big challenges remain. A deep understanding of healthy and aberrant brain activity at a wide range of temporal and spatial scales is needed. Here we discuss, from an interdisciplinary network perspective, the advancements in physical and mathematical modeling as well as in data analysis techniques that, in our opinion, have potential to further advance our understanding of brain structure and function.

Cerebello-Thalamo-Cortical Network Dynamics in the Harmaline Rodent Model of Essential Tremor.

Essential Tremor (ET) is a common movement disorder, characterised by a posture or movement-related tremor of the upper limbs. Abnormalities within cerebellar circuits are thought to underlie the pathogenesis of ET, resulting in aberrant synchronous oscillatory activity within the thalamo-cortical network leading to tremors. Harmaline produces pathological oscillations within the cerebellum, and a tremor that phenotypically resembles ET. However, the neural network dynamics in cerebellar-thalamo-cortical circuits in harmaline-induced tremor remains unclear, including the way circuit interactions may be influenced by behavioural state. Here, we examined the effect of harmaline on cerebello-thalamo-cortical oscillations during rest and movement. EEG recordings from the sensorimotor cortex and local field potentials (LFP) from thalamic and medial cerebellar nuclei were simultaneously recorded in awake behaving rats, alongside measures of tremor using EMG and accelerometery. Analyses compared neural oscillations before and after systemic administration of harmaline (10 mg/kg, I.P), and coherence across periods when rats were resting vs. moving. During movement, harmaline increased the 9-15 Hz behavioural tremor amplitude and increased thalamic LFP coherence with tremor. Medial cerebellar nuclei and cerebellar vermis LFP coherence with tremor however remained unchanged from rest. These findings suggest harmaline-induced cerebellar oscillations are independent of behavioural state and associated changes in tremor amplitude. By contrast, thalamic oscillations are dependent on behavioural state and related changes in tremor amplitude. This study provides new insights into the role of cerebello-thalamo-cortical network interactions in tremor, whereby neural oscillations in thalamocortical, but not cerebellar circuits can be influenced by movement and/or behavioural tremor amplitude in the harmaline model.

Differential Electrographic Signatures Generated by Mechanistically-Diverse Seizurogenic Compounds in the Larval Zebrafish Brain.

We assessed similarities and differences in the electrographic signatures of local field potentials (LFPs) evoked by different pharmacological agents in zebrafish larvae. We then compared and contrasted these characteristics with what is known from electrophysiological studies of seizures and epilepsy in mammals, including humans. Ultimately, our aim was to phenotype neurophysiological features of drug-induced seizures in larval zebrafish for expanding knowledge on the translational potential of this valuable alternative to mammalian models. LFPs were recorded from the midbrain of 4-d-old zebrafish larvae exposed to a pharmacologically diverse panel of seizurogenic compounds, and the outputs of these recordings were assessed using frequency domain analysis. This included analysis of changes occurring within various spectral frequency bands of relevance to mammalian CNS circuit pathophysiology. From these analyses, there were clear differences in the frequency spectra of drug-exposed LFPs, relative to controls, many of which shared notable similarities with the signatures exhibited by mammalian CNS circuits. These similarities included the presence of specific frequency components comparable to those observed in mammalian studies of seizures and epilepsy. Collectively, the data presented provide important information to support the value of larval zebrafish as an alternative model for the study of seizures and epilepsy. These data also provide further insight into the electrophysiological characteristics of seizures generated in nonmammalian species by the action of neuroactive drugs.