A thalamo-cortico-thalamic neural mass model to study alpha rhythms in Alzheimer's disease

Recent Progress in Brain Network Models for Medical Applications: A Review

Importance: Pathological perturbations of the brain often spread via connectome to fundamentally alter functional consequences. By integrating multimodal neuroimaging data with mathematical neural mass modeling, brain network models (BNMs) enable to quantitatively characterize aberrant network dynamics underlying multiple neurological and psychiatric disorders. We delved into the advancements of BNM-based medical applications, discussed the prevalent challenges within this field, and provided possible solutions and future directions. Highlights: This paper reviewed the theoretical foundations and current medical applications of computational BNMs. Composed of neural mass models, the BNM framework allows to investigate large-scale brain dynamics behind brain diseases by linking the simulated functional signals to the empirical neurophysiological data, and has shown promise in exploring neuropathological mechanisms, elucidating therapeutic effects, and predicting disease outcome. Despite that several limitations existed, one promising trend of this research field is to precisely guide clinical neuromodulation treatment based on individual BNM simulation. Conclusion: BNM carries the potential to help understand the mechanism underlying how neuropathology affects brain network dynamics, further contributing to decision-making in clinical diagnosis and treatment. Several constraints must be addressed and surmounted to pave the way for its utilization in the clinic.

Modeling the contribution of theta-gamma coupling to sequential memory, imagination, and dreaming

Gamma oscillations nested in a theta rhythm are observed in the hippocampus, where are assumed to play a role in sequential episodic memory, i.e., memorization and retrieval of events that unfold in time. In this work, we present an original neurocomputational model based on neural masses, which simulates the encoding of sequences of events in the hippocampus and subsequent retrieval by exploiting the theta-gamma code. The model is based on a three-layer structure in which individual Units oscillate with a gamma rhythm and code for individual features of an episode. The first layer (working memory in the prefrontal cortex) maintains a cue in memory until a new signal is presented. The second layer (CA3 cells) implements an auto-associative memory, exploiting excitatory and inhibitory plastic synapses to recover an entire episode from a single feature. Units in this layer are disinhibited by a theta rhythm from an external source (septum or Papez circuit). The third layer (CA1 cells) implements a hetero-associative net with the previous layer, able to recover a sequence of episodes from the first one. During an encoding phase, simulating high-acetylcholine levels, the network is trained with Hebbian (synchronizing) and anti-Hebbian (desynchronizing) rules. During retrieval (low-acetylcholine), the network can correctly recover sequences from an initial cue using gamma oscillations nested inside the theta rhythm. Moreover, in high noise, the network isolated from the environment simulates a mind-wandering condition, randomly replicating previous sequences. Interestingly, in a state simulating sleep, with increased noise and reduced synapses, the network can "dream" by creatively combining sequences, exploiting features shared by different episodes. Finally, an irrational behavior (erroneous superimposition of features in various episodes, like "delusion") occurs after pathological-like reduction in fast inhibitory synapses. The model can represent a straightforward and innovative tool to help mechanistically understand the theta-gamma code in different mental states.

Decrease in phase slip rates and phase cone structures during seizure evolution and epileptogenic activities derived from microgrid ECoG data

Sudden phase changes are related to cortical phase transitions, which likely change in frequency and spatial distribution as epileptogenic activity evolves. A 100 s long section of micro-ECoG data obtained before and during a seizure was selected and analyzed. In addition, nine other short-duration epileptic events were also examined. The data was collected at 420 Hz, imported into MATLAB, downsampled to 200 Hz, and filtered in the 1-50 Hz band. The Hilbert transform was applied to compute the analytic phase, which was then unwrapped, and detrended to look for sudden phase changes. The phase slip rate (counts/s) and its acceleration (counts/s2) were computed with a stepping window of 1-s duration and with a step size of 5 ms. The analysis was performed for theta (3-7 Hz), alpha (7-12 Hz), and beta (12-30 Hz) bands. The phase slip rate on all electrodes in the theta band decreased while it increased for the alpha and beta bands during the seizure period. Similar patterns were observed for isolated epileptogenic events. Spatiotemporal contour plots of the phase slip rates were also constructed using a montage layout of 8 × 8 electrode positions. These plots exhibited dynamic and oscillatory formation of phase cone-like structures which were higher in the theta band and lower in the alpha and beta bands during the seizure period and epileptogenic events. These results indicate that the formation of phase cones might be an excellent biomarker to study the evolution of a seizure and also the cortical dynamics of isolated epileptogenic events.

Distinct Spectral Profiles of Awake Resting EEG in Disorders of Consciousness: The Role of Frequency and Topography of Oscillations

The prolonged disorders of consciousness (PDOC) pose a challenge for an accurate clinical diagnosis, mainly due to patients' scarce or ambiguous behavioral responsiveness. Measurement of brain activity can support better diagnosis, independent of motor restrictions. Methods based on spectral analysis of resting-state EEG appear as a promising path, revealing specific changes within the internal brain dynamics in PDOC patients. In this study we used a robust method of resting-state EEG power spectrum parameter extraction to identify distinct spectral properties for different types of PDOC. Sixty patients and 37 healthy volunteers participated in this study. Patient group consisted of 22 unresponsive wakefulness patients, 25 minimally conscious patients and 13 patients emerging from the minimally conscious state. Ten minutes of resting EEG was acquired during wakefulness and transformed into individual power spectra. For each patient, using the spectral decomposition algorithm, we extracted maximum peak frequency within 1-14 Hz range in the centro-parietal region, and the antero-posterior (AP) gradient of the maximal frequency peak. All patients were behaviorally diagnosed using coma recovery scale-revised (CRS-R). The maximal peak frequency in the 1-14 Hz range successfully predicted both neurobehavioral capacity of patients as indicated by CRS-R total score and PDOC diagnosis. Additionally, in patients in whom only one peak within the 1-14 Hz range was observed, the AP gradient significantly contributed to the accuracy of prediction. We have identified three distinct spectral profiles of patients, likely representing separate neurophysiological modes of thalamocortical functioning. Etiology did not have significant influence on the obtained results.

The important role of glial transmitters released by astrocytes in Alzheimer's disease: A perspective from dynamical modeling

This paper aims to establish a coupling model of neuronal populations and astrocytes and, on this basis, explore the possible mechanism of electroencephalography (EEG) slowing in Alzheimer's disease (AD) from the viewpoint of dynamical modeling. First and foremost, excitatory and inhibitory time constants are shown to induce the early symptoms of AD. The corresponding dynamic nature is mainly due to changes in the amplitude and frequency of the oscillatory behavior. However, there are also a few cases that can be attributed to the change of the oscillation mode caused by the limit cycle bifurcation and birhythmicity. Then, an improved neural mass model influenced by astrocytes is proposed, considering the important effects of glutamate and adenosine triphosphate (ATP) released by astrocytes on the synaptic transmission process reported in experiments. The results show that a dysfunctional astrocyte disrupts the physiological state, causing three typical EEG slowing phenomena reported clinically: the decreased dominant frequency, the decreased rhythmic activity in the α band, and the increased rhythmic activity in the δ+θ band. In addition, astrocytes may control AD when the effect of ATP on synaptic connections is greater than that of glutamate. The control rate depends on the ratio of the effect of glutamate on excitatory and inhibitory synaptic connections. These modeling results can not only reproduce some experimental and clinical results, but, more importantly, may offer a prediction of some underlying phenomena, helping to inspire the disease mechanisms and therapeutic methods of targeting astrocytes.

Frontal EEG alterations induced by hippocampal amyloid pathology in rats

PURPOSE

In this study, it was aimed to determine the dose-dependent effects of hippocampal amyloid beta (Aβ) on frontal EEG activity and to elucidate the possible non-invasive biomarkers by recording spontaneous EEG in free-moving rats.

MATERIAL AND METHODS

Male albino Wistar rats aged 3 months were randomly divided into 4 groups (n ​= ​8 for each group), obtained by intrahippocampal injection of saline or different doses of Aβ1-42 i.e. 0.01 ​μg/μl, 0.1 ​μg/μl, and 1 ​μg/μl. After two weeks of recovery period, spontaneous EEG recordings were obtained from frontal regions and spectral power analyses were performed.

RESULTS

We detected a general slowdown in the brain activity two weeks after Aβ1-42 injection. We observed significant increases in frontal alpha power (p ​= ​0.0021) and significant decreases in frontal beta power (p ​= ​0.0003) between the Sh and Aβ1-42-injected groups. More specifically, the ratio of the frontal EEG beta and alpha power (rFBA) was significantly affected by the intrahippocampal injection of Aβ1-42 (p ​< ​0.0001). Also, we observed that rFBA was negatively and strongly correlated with hippocampal Aβ1-42 peptide levels (r ​= ​-0.781, p ​< ​0.0001).

CONCLUSION

Our findings indicate that spontaneous frontal EEG beta and alpha activity were significantly affected by the intrahippocampal injection of Aβ1-42. However, the results suggest that the power ratios of these bands are more sensitive to the hippocampal amyloid pathology. As such, it is proposed that the rFBA may be a more effective biomarker for diagnosing hippocampal pathology induced by Aβ1-42.

Evolution of brain network dynamics in early Parkinson's disease with mild cognitive impairment

How mild cognitive impairment (MCI) is instantiated in dynamically interacting and spatially distributed functional brain networks remains an unexplored mystery in early Parkinson's disease (PD). We applied a machine-learning technology based on personalized sliding-window algorithm to track continuously time-varying and overlapping subnetworks under the functional brain networks calculated form resting state electroencephalogram data within a sample of 33 early PD patients (13 early PD patients with MCI and 20 early PD patients without MCI). We decoded a set of subnetworks that captured surprisingly dynamically varying and integrated interactions among certain brain lobes. We observed that the master expressed subnetworks were particularly transient, and flexibly switching between high and low expression during integration into a dynamic brain network. This transience was particularly salient in a subnetwork predominantly linking temporal-parietal-occipital lobes, which decreases in both expression and flexibility in early PD patients with MCI and expresses their degree of cognitive impairment. Moreover, MCI induced a regularly interrupted, slow evolution of subnetworks in functional brain network dynamics in early PD at the individual level, and the dynamic expression characteristics of subnetworks also reflected the degree of cognitive impairment in patients with early PD. Collectively, these results provide novel and deeper insights regarding MCI-induced abnormal dynamical interaction and large-scale changes in functional brain network of early PD.

Turning Back the Clock: A Retrospective Single-Blind Study on Brain Age Change in Response to Nutraceuticals Supplementation vs. Lifestyle Modifications

BACKGROUND

There is a growing consensus that chronological age (CA) is not an accurate indicator of the aging process and that biological age (BA) instead is a better measure of an individual's risk of age-related outcomes and a more accurate predictor of mortality than actual CA. In this context, BA measures the "true" age, which is an integrated result of an individual's level of damage accumulation across all levels of biological organization, along with preserved resources. The BA is plastic and depends upon epigenetics. Brain state is an important factor contributing to health- and lifespan.

METHODS AND OBJECTIVE

Quantitative electroencephalography (qEEG)-derived brain BA (BBA) is a suitable and promising measure of brain aging. In the present study, we aimed to show that BBA can be decelerated or even reversed in humans (N = 89) by using customized programs of nutraceutical compounds or lifestyle changes (mean duration = 13 months).

RESULTS

We observed that BBA was younger than CA in both groups at the end of the intervention. Furthermore, the BBA of the participants in the nutraceuticals group was 2.83 years younger at the endpoint of the intervention compared with their BBA score at the beginning of the intervention, while the BBA of the participants in the lifestyle group was only 0.02 years younger at the end of the intervention. These results were accompanied by improvements in mental-physical health comorbidities in both groups. The pre-intervention BBA score and the sex of the participants were considered confounding factors and analyzed separately.

CONCLUSIONS

Overall, the obtained results support the feasibility of the goal of this study and also provide the first robust evidence that halting and reversal of brain aging are possible in humans within a reasonable (practical) timeframe of approximately one year.

A multi-scale model explains oscillatory slowing and neuronal hyperactivity in Alzheimer's disease

Alzheimer's disease is the most common cause of dementia and is linked to the spreading of pathological amyloid-β and tau proteins throughout the brain. Recent studies have highlighted stark differences in how amyloid-β and tau affect neurons at the cellular scale. On a larger scale, Alzheimer's patients are observed to undergo a period of early-stage neuronal hyperactivation followed by neurodegeneration and frequency slowing of neuronal oscillations. Herein, we model the spreading of both amyloid-β and tau across a human connectome and investigate how the neuronal dynamics are affected by disease progression. By including the effects of both amyloid-β and tau pathology, we find that our model explains AD-related frequency slowing, early-stage hyperactivation and late-stage hypoactivation. By testing different hypotheses, we show that hyperactivation and frequency slowing are not due to the topological interactions between different regions but are mostly the result of local neurotoxicity induced by amyloid-β and tau protein.

Analysis of the neural mechanism of spectra decrease in MCI by a thalamo-cortical coupled neural mass model

Objective.In order to deeply understand the neurophysiological mechanism of the spectra decrease in mild cognitive impairment (MCI), this paper studies a new neural mass model, which can simulate various intracerebral electrophysiological activities.Approach. In this study, a thalamo-cortical coupled neural mass model (TCC-NMM) is proposed. The influences of the coupling coefficients and other key parameters on the model spectra are simulated. Then, the unscented Kalman filter (UKF) algorithm is used to reversely identify the parameters in the TCC-NMM. Furthermore, the TCC-NMM and UKF are combined to analyze the spectra reduction mechanism of electroencephalogram (EEG) signals in MCI patients. The independent sample t-test is carried out to statistical analyze the differences of the identified parameters between MCI and normal controls. The Pearson correlation analysis is used to analyze the intrinsic relationship between parameters and the scores of the comprehensive competence assessment scale.Main results.The simulation results show that the decreased cortical synaptic connectivity constantsC₁can result in spectra decrease of the TCC-NMM outputs. The real EEG analysis results show that the identified values of parameterC₁are significant lower in the MCI group than in control group in frontal and occipital areas and the parametersC₁are positively correlated with the Montreal Cognitive Assessment (MoCA) scores in the two areas. This consistency suggests that the cortical synaptic connectivity loss from pyramidal cells to excitatory interneurons (eIN) may be one of the neural mechanisms of EEG spectra decrease in MCI.Significance. (a) In this study, a new mathematical model TCCNMM based on anatomy and neurophysiology is proposed. (b) All key parameters in TCC-NMM are studied in detail through the forward and reverse analysis and the influence of these parameters on the output spectra of the model is pointed out. (c) The possible neural mechanism of the decreased spectra in MCI patients is pointed out by the joint analysis of simulation in forward with TCC-NMM and analysis of the actual EEG signals in reverse with UKF identification algorithm. (d) We find that the identified parameter C1 of MCI patients is significantly lower than that of the control group, which is consistent with the simulation analysis of TCC-NMM. So, we suggest that the decreased MCI alpha power spectrum is likely related to the cortical synaptic connection loss from pyramidal cells to eIN.

Signal processing and computational modeling for interpretation of SEEG-recorded interictal epileptiform discharges in epileptogenic and non-epileptogenic zones

OBJECTIVE

In partial epilepsies, interictal epileptiform discharges (IEDs) are paroxysmal events observed in epileptogenic and non-epileptogenic zones. IEDs' generation and recurrence are subject to different hypotheses: they appear through glutamatergic and GABAergic processes; they may trigger seizures or prevent seizure propagation. This paper focuses on a specific class of IEDs, spike-waves (SWs), characterized by a short-duration spike followed by a longer duration wave, both of the same polarity. Signal analysis and neurophysiological mathematical models are used to interpret puzzling IED generation.

APPROACH

Interictal activity was recorded by intracranial stereo-electroencephalography (SEEG) electrodes in five different patients. SEEG experts identified the epileptic and non-epileptic zones in which IEDs were detected. After quantifying spatial and temporal features of the detected IEDs, the most significant features for classifying epileptic and non-epileptic zones were determined. A neurophysiologically-plausible mathematical model was then introduced to simulate the IEDs and understand the underlying differences observed in epileptic and non-epileptic zone IEDs.

MAIN RESULTS

Two classes of SWs were identified according to subtle differences in morphology and timing of the spike and wave component. Results showed that type-1 SWs were generated in epileptogenic regions also involved at seizure onset, while type-2 SWs were produced in the propagation or non-involved areas. The modeling study indicated that synaptic kinetics, cortical organization, and network interactions determined the morphology of the simulated SEEG signals. Modeling results suggested that the IED morphologies were linked to the degree of preserved inhibition.

SIGNIFICANCE

This work contributes to the understanding of different mechanisms generating IEDs in epileptic networks. The combination of signal analysis and computational models provides an efficient framework for exploring IEDs in partial epilepsies and classifying epileptogenic and non-epileptogenic zones.

Effect of acetylcholine deficiency on neural oscillation in a brainstem-thalamus-cortex neurocomputational model related with Alzheimer's disease

Previous works imply that involving brainstem in neuropathological studies of Alzheimer's disease (AD) is of clinically significant. This work constructs a comprehensive neural mass model for cholinergic neuropathogenesis that involves brainstem, thalamus and cortex, wherein how acetylcholine deficiency in AD affects neural oscillation of the model output is systematically explored from the perspective of neurocomputation. By decreasing synapse connectivity parameters in direct cholinergic pathway from brainstem to thalamus or in indirect glutamatergic synapse pathway from cortex to brainstem to mimic the pathological condition of reduced acetylcholine release in patients with AD, the property of neural oscillation in this model is numerically investigated by means of power spectrum in frequency domain and amplitude distribution in time domain. Simulated results demonstrate that decreasing synapse connectivity whether in the direct cholinergic pathway or in the indirect glutamatergic synapse pathway can alter the neural oscillation significantly in three aspects: it induces an obvious decrease of dominant frequency; it leads to a degraded rhythmic activity in the alpha frequency band as well as an enhanced rhythmic activity in the theta frequency band; it results in reduced oscillation amplitude of the model output. These results are agreement with the characteristic of electrophysiological EEG measurement recorded in AD, especially support the hypothesis that cholinergic deficiency is a promising pathophysiological origin of EEG slowing in AD. Our analysis indicates that targeting the cholinergic system may have potential prospects in early diagnosis and treatment of AD.

Normalized Theta but Increased Gamma Activity after Acetylcholinesterase Inhibitor Treatment in Alzheimer's Disease: Preliminary qEEG Study

Acetylcholinesterase inhibitors (AChE-I) are the core treatment of mild to severe Alzheimer's disease (AD). However, the efficacy of AChE-I treatment on electroencephalography (EEG) and cognition remains unclear. We aimed to investigate the EEG power and coherence changes, in addition to neuropsychological performance, following a one-year treatment. Nine de-novo AD patients and demographically-matched healthy controls (HC) were included. After baseline assessments, all AD participants started cholinergic therapy. We found that baseline and follow-up gamma power analyzes were similar between groups. Yet, within the AD group after AChE-I intake, individuals with AD displayed higher gamma power compared to their baselines (P < .039). Also, baseline gamma coherence analysis showed lower values in the AD than in HC (P < .048), while these differences disappeared with increased gamma values of AD patients at the follow-up. Within the AD group after AChE-I intake, individuals with AD displayed higher theta and alpha coherence compared to their baselines (all, P < .039). These increased results within the AD group may result from a subclinical epileptiform activity. Even though AChE-I is associated with lower mortality, our results showed a significant effect on EEG power yet can increase the subclinical epileptiform activity. It is essential to be conscious of the seizure risk that treatment may cause.

A novelty detection approach to effectively predict conversion from mild cognitive impairment to Alzheimer’s disease

Accurately recognising patients with progressive mild cognitive impairment (pMCI) who will develop Alzheimer’s disease (AD) in subsequent years is very important, as early identification of those patients will enable interventions to potentially reduce the number of those transitioning from MCI to AD. Most studies in this area have concentrated on high-dimensional neuroimaging data with supervised binary/multi-class classification algorithms. However, neuroimaging data is more costly to obtain than non-imaging, and healthcare datasets are normally imbalanced which may reduce classification performance and reliability. To address these challenges, we proposed a new strategy that employs unsupervised novelty detection (ND) techniques to predict pMCI from the AD neuroimaging initiative non-imaging data. ND algorithms, including the k-nearest neighbours (kNN), k-means, Gaussian mixture model (GMM), isolation forest (IF) and extreme learning machine (ELM), were employed and compared with supervised binary support vector machine (SVM) and random forest (RF). We introduced optimisation with nested cross-validation and focused on maximising the adjusted F measure to ensure maximum generalisation of the proposed system by minimising false negative rates. Our extensive experimental results show that ND algorithms (0.727 ± 0.029 kNN, 0.7179 ± 0.0523 GMM, 0.7276 ± 0.0281 ELM) obtained comparable performance to supervised binary SVM (0.7359 ± 0.0451) with 20% stable MCI misclassification tolerance and were significantly better than RF (0.4771 ± 0.0167). Moreover, we found that the non-invasive, readily obtainable, and cost-effective cognitive and functional assessment was the most efficient predictor for predicting the pMCI within 2 years with ND techniques. Importantly, we presented an accessible and cost-effective approach to pMCI prediction, which does not require labelled data.

Biophysical mechanism underlying compensatory preservation of neural synchrony over the adult lifespan

We propose that the preservation of functional integration, estimated from measures of neural synchrony, is a key objective of neurocompensatory mechanisms associated with healthy human ageing. To support this proposal, we demonstrate how phase-locking at the peak alpha frequency in Magnetoencephalography recordings remains invariant over the lifespan in a large cohort of human participants, aged 18-88 years. Using empirically derived connection topologies from diffusion tensor imaging data, we create an in-silico model of whole-brain alpha dynamics. We show that enhancing inter-areal coupling can cancel the effect of increased axonal transmission delays associated with age-related degeneration of white matter tracts, albeit at slower network frequencies. By deriving analytical solutions for simplified connection topologies, we further establish the theoretical principles underlying compensatory network re-organization. Our findings suggest that frequency slowing with age- frequently observed in the alpha band in diverse populations- may be viewed as an epiphenomenon of the underlying compensatory mechanism.

Analysis of MEG data from healthy participants and whole-brain network modeling suggests that the brain compensates for age-related disruptions in connectivity by slowing down the frequency of neural synchronization.

In silico Effects of Synaptic Connections in the Visual Thalamocortical Pathway

We have studied brain connectivity using a biologically inspired in silico model of the visual pathway consisting of the lateral geniculate nucleus (LGN) of the thalamus, and layers 4 and 6 of the primary visual cortex. The connectivity parameters in the model are informed by the existing anatomical parameters from mammals and rodents. In the base state, the LGN and layer 6 populations in the model oscillate with dominant alpha frequency, while the layer 4 oscillates in the theta band. By changing intra-cortical hyperparameters, specifically inhibition from layer 6 to layer 4, we demonstrate a transition to alpha mode for all the populations. Furthermore, by increasing the feedforward connectivities in the thalamo-cortico-thalamic loop, we could transition into the beta band for all the populations. On looking closely, we observed that the origin of this beta band is in the layer 6 (infragranular layers); lesioning the thalamic feedback from layer 6 removed the beta from the LGN and the layer 4. This agrees with existing physiological studies where it is shown that beta rhythm is generated in the infragranular layers. Lastly, we present a case study to demonstrate a neurological condition in the model. By changing connectivities in the network, we could simulate the condition of significant (P < 0.001) decrease in beta band power and a simultaneous increase in the theta band power, similar to that observed in Schizophrenia patients. Overall, we have shown that the connectivity changes in a simple visual thalamocortical in silico model can simulate state changes in the brain corresponding to both health and disease conditions.

Subcortical and Cerebellar Neural Correlates of Prodromal Alzheimer’s Disease with Prolonged Sleep Latency

Background: Despite the important associations among sleep, Alzheimer’s disease (AD), subcortical structures, and the cerebellum, structural and functional magnetic resonance imaging (MRI) with regard to these regions and sleep on patients in AD trajectory are scarce. Objective: This study aimed to evaluate the influence of prolonged sleep latency on the structural and functional alterations in the subcortical and cerebellar neural correlates in amyloid-β positive amnestic mild cognitive impairment patients (Aβ+aMCI). Methods: A total of 60 patients with aMCI who were identified as amyloid positive ([18F] flutemetamol+) were recruited in the study, 24 patients with normal sleep latency (aMCI-n) and 36 patients prolonged sleep latency (aMCI-p). Cortical thickness and volumes between the two groups were compared. Volumetric analyses were implemented on the brainstem, thalamus, and hippocampus. Subcortical and cerebellar resting state functional connectivity (FC) differences were measured between the both groups through seed-to-voxel analysis. Additionally, group x Aβ interactive effects on FC values were tested with a general linear model. Result: There was a significantly decreased brainstem volume in aMCI-p subjects. We observed a significant reduction of the locus coeruleus (LC) FC with frontal, temporal, insular cortices, hippocampus, and left thalamic FC with occipital cortex. Moreover, the LC FC with occipital cortex and left hippocampal FC with frontal cortex were increased in aMCI-p subjects. In addition, there was a statistically significant group by regional standardized uptake value ratio interactions discovered in cerebro-cerebellar networks. Conclusion: The aforementioned findings suggest that prolonged sleep latency may be a detrimental factor in compromising structural and functional correlates of subcortical structures and the cerebellum, which may accelerate AD pathophysiology.

Controlling Alzheimer's Disease Through the Deep Brain Stimulation to Thalamic Relay Cells

Experimental and clinical studies have shown that the technique of deep brain stimulation (DBS) plays a potential role in the regulation of Alzheimer’s disease (AD), yet it still desires for ongoing studies including clinical trials, theoretical approach and action mechanism. In this work, we develop a modified thalamo-cortico-thalamic (TCT) model associated with AD to explore the therapeutic effects of DBS on AD from the perspective of neurocomputation. First, the neuropathological state of AD resulting from synapse loss is mimicked by decreasing the synaptic connectivity strength from the Inter-Neurons (IN) neuron population to the Thalamic Relay Cells (TRC) neuron population. Under such AD condition, a specific deep brain stimulation voltage is then implanted into the neural nucleus of TRC in this TCT model. The symptom of AD is found significantly relieved by means of power spectrum analysis and nonlinear dynamical analysis. Furthermore, the therapeutic effects of DBS on AD are systematically examined in different parameter space of DBS. The results demonstrate that the controlling effect of DBS on AD can be efficient by appropriately tuning the key parameters of DBS including amplitude A, period P and duration D. This work highlights the critical role of thalamus stimulation for brain disease, and provides a theoretical basis for future experimental and clinical studies in treating AD.

Canard solutions in neural mass models: consequences on critical regimes

Mathematical models at multiple temporal and spatial scales can unveil the fundamental mechanisms of critical transitions in brain activities. Neural mass models (NMMs) consider the average temporal dynamics of interconnected neuronal subpopulations without explicitly representing the underlying cellular activity. The mesoscopic level offered by the neural mass formulation has been used to model electroencephalographic (EEG) recordings and to investigate various cerebral mechanisms, such as the generation of physiological and pathological brain activities. In this work, we consider a NMM widely accepted in the context of epilepsy, which includes four interacting neuronal subpopulations with different synaptic kinetics. Due to the resulting three-time-scale structure, the model yields complex oscillations of relaxation and bursting types. By applying the principles of geometric singular perturbation theory, we unveil the existence of the canard solutions and detail how they organize the complex oscillations and excitability properties of the model. In particular, we show that boundaries between pathological epileptic discharges and physiological background activity are determined by the canard solutions. Finally we report the existence of canard-mediated small-amplitude frequency-specific oscillations in simulated local field potentials for decreased inhibition conditions. Interestingly, such oscillations are actually observed in intracerebral EEG signals recorded in epileptic patients during pre-ictal periods, close to seizure onsets.

Decreased Alpha Peak Frequency Is Linked to Episodic Memory Impairment in Pathological Aging

The Free and Cued Selective Reminding Test (FCSRT) is a largely validated neuropsychological test for the identification of amnestic syndrome from the early stage of Alzheimer's disease (AD). Previous electrophysiological data suggested a slowing down of the alpha rhythm in the AD-continuum as well as a key role of this rhythmic brain activity for episodic memory processes. This study therefore investigates the link between alpha brain activity and alterations in episodic memory as assessed by the FCSRT. For that purpose, 37 patients with altered FCSRT performance underwent a comprehensive neuropsychological assessment, supplemented by 18F-fluorodeoxyglucose positron emission tomography/structural magnetic resonance imaging (18FDG-PET/MR), and 10 min of resting-state magnetoencephalography (MEG). The individual alpha peak frequency (APF) in MEG resting-state data was positively correlated with patients' encoding efficiency as well as with the efficacy of semantic cues in facilitating patients' retrieval of previous stored word. The APF also correlated positively with patients' hippocampal volume and their regional glucose consumption in the posterior cingulate cortex. Overall, this study demonstrates that alterations in the ability to learn and store new information for a relatively short-term period are related to a slowing down of alpha rhythmic activity, possibly due to altered interactions in the extended mnemonic system. As such, a decreased APF may be considered as an electrophysiological correlate of short-term episodic memory dysfunction accompanying pathological aging.

Computational neuronal correlation with enhanced synchronized activity in the basal ganglia and the slowing of thalamic theta and alpha rhythms in Parkinson's disease

The aim of this work is computationally to correlate the synchronized neuronal activity of basal ganglia and slowing in theta and alpha rhythms in electroencephalogram (EEG) signal in thalamic region in case of dopamine depletion and decrease of synaptic connections. The used network topology is a scale-free network with constant node degree. The dopamine-modulated type Izikhevich neuron model is used for modeling the striatal region, consisting of fast-spiking interneurons, D1 and D2 type dopamine expressing medium spiny neurons. On the other hand, the ordinary Izikhevich neuron model is used in the modeling of extrastriatal basal ganglia (BG) regions where globus pallidus (GP) subregion neurons have also dopamine-dependent parameters. The thalamic region of the network is mass modeled including inhibitory input from basal ganglia. Depending on the decrease of synaptic connections and dopamine level, the synchronization among basal ganglia neuron populations is investigated. The effect of synaptic delay on synchronization is also considered. It is observed that the decrease of dopamine neurotransmitter and decrease in the number of synaptic connections cause an increased synchronous activity in BG. Also, slowing in theta and alpha bands in thalamus EEG signals is observed. This shows the causal relation between synchronization and power shifting to lower frequency components in the case of neurodegenerative diseases such as Parkinson's disease (PD).

A novel density-based neural mass model for simulating neuronal network dynamics with conductance-based synapses and membrane current adaptation

Despite its success in understanding brain rhythms, the neural mass model, as a low-dimensional mean-field network model, is phenomenological in nature, so that it cannot replicate some of rich repertoire of responses seen in real neuronal tissues. Here, using a colored-synapse population density method, we derived a novel neural mass model, termed density-based neural mass model (dNMM), as the mean-field description of network dynamics of adaptive exponential integrate-and-fire (aEIF) neurons, in which two critical neuronal features, i.e., voltage-dependent conductance-based synaptic interactions and adaptation of firing rate responses, were included. Our results showed that the dNMM was capable of correctly estimating firing rate responses of a neuronal population of aEIF neurons receiving stationary or time-varying excitatory and inhibitory inputs. Finally, it was also able to quantitatively describe the effect of spike-frequency adaptation in the generation of asynchronous irregular activity of excitatory-inhibitory cortical networks. We conclude that in terms of its biological reality and calculation efficiency, the dNMM is a suitable candidate to build significantly large-scale network models involving multiple brain areas, where the neuronal population is the smallest dynamic unit.

NMMGenerator: an automatic neural mass model generator from population graphs

Neural mass models are among the most popular mathematical models of brain activity, since they enable the rapid simulation of large-scale networks involving different neural types at a spatial scale compatible with electrophysiological experiments (e.g. local field potentials). However, establishing neural mass model (NMM) equations associated with specific neuronal network architectures can be tedious and is an error-prone process, restricting their use to scientists who are familiar with mathematics. In order to overcome this challenge, we have developed a user-friendly software that enables a user to construct rapidly, under the form of a graph, a neuronal network with its populations and connectivity patterns. The resulting graph is then automatically translated into the corresponding set of differential equations, which can be solved and displayed within the same software environment. The software is proposed as open access, and should assist in offering the possibility for a wider audience of scientists to develop NMM corresponding to their specific neuroscience research questions.

Modeling Neurotransmission: Computational Tools to Investigate Neurological Disorders

The investigation of synaptic functions remains one of the most fascinating challenges in the field of neuroscience and a large number of experimental methods have been tuned to dissect the mechanisms taking part in the neurotransmission process. Furthermore, the understanding of the insights of neurological disorders originating from alterations in neurotransmission often requires the development of (i) animal models of pathologies, (ii) invasive tools and (iii) targeted pharmacological approaches. In the last decades, additional tools to explore neurological diseases have been provided to the scientific community. A wide range of computational models in fact have been developed to explore the alterations of the mechanisms involved in neurotransmission following the emergence of neurological pathologies. Here, we review some of the advancements in the development of computational methods employed to investigate neuronal circuits with a particular focus on the application to the most diffuse neurological disorders.

Dysregulation of excitatory neural firing replicates physiological and functional changes in aging visual cortex

The mammalian visual system has been the focus of countless experimental and theoretical studies designed to elucidate principles of neural computation and sensory coding. Most theoretical work has focused on networks intended to reflect developing or mature neural circuitry, in both health and disease. Few computational studies have attempted to model changes that occur in neural circuitry as an organism ages non-pathologically. In this work we contribute to closing this gap, studying how physiological changes correlated with advanced age impact the computational performance of a spiking network model of primary visual cortex (V1). Our results demonstrate that deterioration of homeostatic regulation of excitatory firing, coupled with long-term synaptic plasticity, is a sufficient mechanism to reproduce features of observed physiological and functional changes in neural activity data, specifically declines in inhibition and in selectivity to oriented stimuli. This suggests a potential causality between dysregulation of neuron firing and age-induced changes in brain physiology and functional performance. While this does not rule out deeper underlying causes or other mechanisms that could give rise to these changes, our approach opens new avenues for exploring these underlying mechanisms in greater depth and making predictions for future experiments.

Neural mass modeling of slow-fast dynamics of seizure initiation and abortion

Epilepsy is a dynamic and complex neurological disease affecting about 1% of the worldwide population, among which 30% of the patients are drug-resistant. Epilepsy is characterized by recurrent episodes of paroxysmal neural discharges (the so-called seizures), which manifest themselves through a large-amplitude rhythmic activity observed in depth-EEG recordings, in particular in local field potentials (LFPs). The signature characterizing the transition to seizures involves complex oscillatory patterns, which could serve as a marker to prevent seizure initiation by triggering appropriate therapeutic neurostimulation methods. To investigate such protocols, neurophysiological lumped-parameter models at the mesoscopic scale, namely neural mass models, are powerful tools that not only mimic the LFP signals but also give insights on the neural mechanisms related to different stages of seizures. Here, we analyze the multiple time-scale dynamics of a neural mass model and explain the underlying structure of the complex oscillations observed before seizure initiation. We investigate population-specific effects of the stimulation and the dependence of stimulation parameters on synaptic timescales. In particular, we show that intermediate stimulation frequencies (>20 Hz) can abort seizures if the timescale difference is pronounced. Those results have the potential in the design of therapeutic brain stimulation protocols based on the neurophysiological properties of tissue.

Thalamocortical inhibitory dynamics support conscious perception

Whether thalamocortical interactions play a decisive role in conscious perception remains an open question. We presented rapid red/green color flickering stimuli, which induced the mental perception of either an illusory orange color or non-fused red and green colors. Using magnetoencephalography, we observed 6-Hz thalamic activity associated with thalamocortical inhibitory coupling only during the conscious perception of the illusory orange color. This sustained thalamic disinhibition was temporally coupled with higher visual cortical activation during the conscious perception of the orange color, providing neurophysiological evidence of the role of thalamocortical synchronization in conscious awareness of mental representation. Bayesian model comparison consistently supported the thalamocortical model in conscious perception. Taken together, experimental and theoretical evidence established the thalamocortical inhibitory network as a gateway to conscious mental representations.

Transfer Entropy as a Measure of Brain Connectivity: A Critical Analysis With the Help of Neural Mass Models

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.

The Human Body as a Super Network: Digital Methods to Analyze the Propagation of Aging

Biological aging is a complex process involving multiple biological processes. These can be understood theoretically though considering them as individual networks-e.g., epigenetic networks, cell-cell networks (such as astroglial networks), and population genetics. Mathematical modeling allows the combination of such networks so that they may be studied in unison, to better understand how the so-called "seven pillars of aging" combine and to generate hypothesis for treating aging as a condition at relatively early biological ages. In this review, we consider how recent progression in mathematical modeling can be utilized to investigate aging, particularly in, but not exclusive to, the context of degenerative neuronal disease. We also consider how the latest techniques for generating biomarker models for disease prediction, such as longitudinal analysis and parenclitic analysis can be applied to as both biomarker platforms for aging, as well as to better understand the inescapable condition. This review is written by a highly diverse and multi-disciplinary team of scientists from across the globe and calls for greater collaboration between diverse fields of research.

Alpha rhythm slowing in a modified thalamo-cortico-thalamic model related with Alzheimer's disease

A decrease in alpha band power is defined as a hallmark of electroencephalogram (EEG) in Alzheimer’s disease (AD). This study devotes to understanding the neuronal correlates of alpha rhythm slowing associated with AD from the view of neurocomputation. Firstly, a modified computational model of thalamo-cortico-thalamic (TCT) circuitry is constructed by incorporating two important biologically plausible ingredients. One is the disinhibition property between different inhibitory interneurons in the cortical module. The other is the full relay function of thalamic relay nucleus (TCR) to the cortical module. Then, by decreasing synaptic connectivity parameters to mimic the neuropathological condition of synapse loss in AD, the correlation between neuronal synaptic behavior and abnormal alpha rhythm is simulated by means of power spectral analysis. The results indicate that these decreases of synaptic activity, i.e., not only the excitatory synaptic connections from TCR to fast inhibitory interneurons C_fte and from excitatory interneurons to pyramidal neurons C_pxe but also the inhibitory synaptic connections from fast inhibitory interneurons to slow inhibitory interneurons C_lfi and from inhibitory interneurons to TCR C_tii, can significantly diminish the peak power density over the alpha band of the thalamic output, which implies that there is a slowing of alpha band. Furthermore, the underlying mechanism behind the alpha rhythmic changes is analyzed using nonlinear dynamical technique. The results reveal that decreases of C_fte, C_pxe, C_lfi and C_tii can make the thalamic module transfer from a limit cycle mode to a point attractor mode, which may lead to the alpha rhythm slowing in the modified TCT model. We expect this work can be helpful in identifying early biomarkers of AD’s EEG and understanding potential pathogenesis of AD.

Experimental Study of Reinforcement Learning in Mobile Robots Through Spiking Architecture of Thalamo-Cortico-Thalamic Circuitry of Mammalian Brain

In this paper, the behavioral learning of robots through spiking neural networks is studied in which the architecture of the network is based on the thalamo-cortico-thalamic circuitry of the mammalian brain. According to a variety of neurons, the Izhikevich model of single neuron is used for the representation of neuronal behaviors. One thousand and ninety spiking neurons are considered in the network. The spiking model of the proposed architecture is derived and prepared for the learning problem of robots. The reinforcement learning algorithm is based on spike-timing-dependent plasticity and dopamine release as a reward. It results in strengthening the synaptic weights of the neurons that are involved in the robot’s proper performance. Sensory and motor neurons are placed in the thalamus and cortical module, respectively. The inputs of thalamo-cortico-thalamic circuitry are the signals related to distance of the target from robot, and the outputs are the velocities of actuators. The target attraction task is used as an example to validate the proposed method in which dopamine is released when the robot catches the target. Some simulation studies, as well as experimental implementation, are done on a mobile robot named Tabrizbot. Experimental studies illustrate that after successful learning, the meantime of catching target is decreased by about 36%. These prove that through the proposed method, thalamo-cortical structure could be trained successfully to learn to perform various robotic tasks.

COALIA: A Computational Model of Human EEG for Consciousness Research

Understanding the origin of the main physiological processes involved in consciousness is a major challenge of contemporary neuroscience, with crucial implications for the study of Disorders of Consciousness (DOC). The difficulties in achieving this task include the considerable quantity of experimental data in this field, along with the non-intuitive, nonlinear nature of neuronal dynamics. One possibility of integrating the main results from the experimental literature into a cohesive framework, while accounting for nonlinear brain dynamics, is the use of physiologically-inspired computational models. In this study, we present a physiologically-grounded computational model, attempting to account for the main micro-circuits identified in the human cortex, while including the specificities of each neuronal type. More specifically, the model accounts for thalamo-cortical (vertical) regulation of cortico-cortical (horizontal) connectivity, which is a central mechanism for brain information integration and processing. The distinct neuronal assemblies communicate through feedforward and feedback excitatory and inhibitory synaptic connections implemented in a template brain accounting for long-range connectome. The EEG generated by this physiologically-based simulated brain is validated through comparison with brain rhythms recorded in humans in two states of consciousness (wakefulness, sleep). Using the model, it is possible to reproduce the local disynaptic disinhibition of basket cells (fast GABAergic inhibition) and glutamatergic pyramidal neurons through long-range activation of vasoactive intestinal-peptide (VIP) interneurons that induced inhibition of somatostatin positive (SST) interneurons. The model (COALIA) predicts that the strength and dynamics of the thalamic output on the cortex control the local and long-range cortical processing of information. Furthermore, the model reproduces and explains clinical results regarding the complexity of transcranial magnetic stimulation TMS-evoked EEG responses in DOC patients and healthy volunteers, through a modulation of thalamo-cortical connectivity that governs the level of cortico-cortical communication. This new model provides a quantitative framework to accelerate the study of the physiological mechanisms involved in the emergence, maintenance and disruption (sleep, anesthesia, DOC) of consciousness.

Computational modeling and biomarker studies of pharmacological treatment of Alzheimer's disease (Review)

Alzheimer's disease (AD) is a complex and multifactorial disease. In order to understand the genetic influence in the progression of AD, and to identify novel pharmaceutical agents and their associated targets, the present study discusses computational modeling and biomarker evaluation approaches. Based on mechanistic signaling pathway approaches, various computational models, including biochemical and morphological models, are discussed to explore the strategies that may be used to target AD treatment. Different biomarkers are interpreted on the basis of morphological and functional features of amyloid β plaques and unstable microtubule‑associated tau protein, which is involved in neurodegeneration. Furthermore, imaging and cerebrospinal fluids are also considered to be key methods in the identification of novel markers for AD. In conclusion, the present study reviews various biochemical and morphological computational models and biomarkers to interpret novel targets and agonists for the treatment of AD. This review also highlights several therapeutic targets and their associated signaling pathways in AD, which may have potential to be used in the development of novel pharmacological agents for the treatment of patients with AD. Computational modeling approaches may aid the quest for the development of AD treatments with enhanced therapeutic efficacy and reduced toxicity.

Analytically determining frequency and amplitude of spontaneous alpha oscillation in Jansen's neural mass model using the describing function method

Spontaneous alpha oscillations are a ubiquitous phenomenon in the brain and play a key role in neural information processing and various cognitive functions. Jansen's neural mass model (NMM) was initially proposed to study the origin of alpha oscillations. Most of previous studies of the spontaneous alpha oscillations in the NMM were conducted using numerical methods. In this study, we aim to propose an analytical approach using the describing function method to elucidate the spontaneous alpha oscillation mechanism in the NMM. First, the sigmoid nonlinear function in the NMM is approximated by its describing function, allowing us to reformulate the NMM and derive its standard form composed of one nonlinear part and one linear part. Second, by conducting a theoretical analysis, we can assess whether or not the spontaneous alpha oscillation would occur in the NMM and, furthermore, accurately determine its amplitude and frequency. The results reveal analytically that the interaction between linearity and nonlinearity of the NMM plays a key role in generating the spontaneous alpha oscillations. Furthermore, strong nonlinearity and large linear strength are required to generate the spontaneous alpha oscillations.

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.

Disrupted Thalamus White Matter Anatomy and Posterior Default Mode Network Effective Connectivity in Amnestic Mild Cognitive Impairment

Alzheimer’s disease (AD) and its prodromal state amnestic mild cognitive impairment (aMCI) are characterized by widespread abnormalities in inter-areal white matter fiber pathways and parallel disruption of default mode network (DMN) resting state functional and effective connectivity. In healthy subjects, DMN and task positive network interaction are modulated by the thalamus suggesting that abnormal task-based DMN deactivation in aMCI may be a consequence of impaired thalamo-cortical white matter circuitry. Thus, this article uses a multimodal approach to assess white matter integrity between thalamus and DMN components and associated effective connectivity in healthy controls (HCs) relative to aMCI patients. Twenty-six HC and 20 older adults with aMCI underwent structural, functional and diffusion MRI scanning using the high angular resolution diffusion-weighted acquisition protocol. The DMN of each subject was identified using independent component analysis (ICA) and resting state effective connectivity was calculated between thalamus and DMN nodes. White matter integrity changes between thalamus and DMN were investigated with constrained spherical deconvolution (CSD) tractography. Significant structural deficits in thalamic white matter projection fibers to posterior DMN components posterior cingulate cortex (PCC) and lateral inferior parietal lobe (IPL) were identified together with significantly reduced effective connectivity from left thalamus to left IPL. Crucially, impaired thalamo-cortical white matter circuitry correlated with memory performance. Disrupted thalamo-cortical structure was accompanied by significant reductions in IPL and PCC cortico-cortical effective connectivity. No structural deficits were found between DMN nodes. Abnormal posterior DMN activity may be driven by changes in thalamic white matter connectivity; a view supported by the close anatomical and functional association of thalamic nuclei effected by AD pathology and the posterior DMN nodes. We conclude that dysfunctional posterior DMN activity in aMCI is consistent with disrupted cortico-thalamo-cortical processing and thalamic-based dissemination of hippocampal disease agents to cortical hubs.

Early Changes in Alpha Band Power and DMN BOLD Activity in Alzheimer's Disease: A Simultaneous Resting State EEG-fMRI Study

Simultaneous resting state functional magnetic resonance imaging (rsfMRI)-resting state electroencephalography (rsEEG) studies in healthy adults showed robust positive associations of signal power in the alpha band with BOLD signal in the thalamus, and more heterogeneous associations in cortical default mode network (DMN) regions. Negative associations were found in occipital regions. In Alzheimer's disease (AD), rsfMRI studies revealed a disruption of the DMN, while rsEEG studies consistently reported a reduced power within the alpha band. The present study is the first to employ simultaneous rsfMRI-rsEEG in an AD sample, investigating the association of alpha band power and BOLD signal, compared to healthy controls (HC). We hypothesized to find reduced positive associations in DMN regions and reduced negative associations in occipital regions in the AD group. Simultaneous resting state fMRI-EEG was recorded in 14 patients with mild AD and 14 HC, matched for age and gender. Power within the EEG alpha band (8-12 Hz, 8-10 Hz, and 10-12 Hz) was computed from occipital electrodes and served as regressor in voxel-wise linear regression analyses, to assess the association with the BOLD signal. Compared to HC, the AD group showed significantly decreased positive associations between BOLD signal and occipital alpha band power in clusters in the superior, middle and inferior frontal cortex, inferior temporal lobe and thalamus (p < 0.01, uncorr., cluster size ≥ 50 voxels). This group effect was more pronounced in the upper alpha sub-band, compared to the lower alpha sub-band. Notably, we observed a high inter-individual heterogeneity. Negative associations were only reduced in the lower alpha range in the hippocampus, putamen and cerebellum. The present study gives first insights into the relationship of resting-state EEG and fMRI characteristics in an AD sample. The results suggest that positive associations between alpha band power and BOLD signal in numerous regions, including DMN regions, are diminished in AD.

A Spiking Neural Network Model of the Lateral Geniculate Nucleus on the SpiNNaker Machine

We present a spiking neural network model of the thalamic Lateral Geniculate Nucleus (LGN) developed on SpiNNaker, which is a state-of-the-art digital neuromorphic hardware built with very-low-power ARM processors. The parallel, event-based data processing in SpiNNaker makes it viable for building massively parallel neuro-computational frameworks. The LGN model has 140 neurons representing a "basic building block" for larger modular architectures. The motivation of this work is to simulate biologically plausible LGN dynamics on SpiNNaker. Synaptic layout of the model is consistent with biology. The model response is validated with existing literature reporting entrainment in steady state visually evoked potentials (SSVEP)-brain oscillations corresponding to periodic visual stimuli recorded via electroencephalography (EEG). Periodic stimulus to the model is provided by: a synthetic spike-train with inter-spike-intervals in the range 10-50 Hz at a resolution of 1 Hz; and spike-train output from a state-of-the-art electronic retina subjected to a light emitting diode flashing at 10, 20, and 40 Hz, simulating real-world visual stimulus to the model. The resolution of simulation is 0.1 ms to ensure solution accuracy for the underlying differential equations defining Izhikevichs neuron model. Under this constraint, 1 s of model simulation time is executed in 10 s real time on SpiNNaker; this is because simulations on SpiNNaker work in real time for time-steps dt ⩾ 1 ms. The model output shows entrainment with both sets of input and contains harmonic components of the fundamental frequency. However, suppressing the feed-forward inhibition in the circuit produces subharmonics within the gamma band (>30 Hz) implying a reduced information transmission fidelity. These model predictions agree with recent lumped-parameter computational model-based predictions, using conventional computers. Scalability of the framework is demonstrated by a multi-node architecture consisting of three "nodes," where each node is the "basic building block" LGN model. This 420 neuron model is tested with synthetic periodic stimulus at 10 Hz to all the nodes. The model output is the average of the outputs from all nodes, and conforms to the above-mentioned predictions of each node. Power consumption for model simulation on SpiNNaker is ≪1 W.

A neural mass model of cross frequency coupling

Electrophysiological signals of cortical activity show a range of possible frequency and amplitude modulations, both within and across regions, collectively known as cross-frequency coupling. To investigate whether these modulations could be considered as manifestations of the same underlying mechanism, we developed a neural mass model. The model provides five out of the theoretically proposed six different coupling types. Within model components, slow and fast activity engage in phase-frequency coupling in conditions of low ambient noise level and with high noise level engage in phase-amplitude coupling. Between model components, these couplings can be coordinated via slow activity, giving rise to more complex modulations. The model, thus, provides a coherent account of cross-frequency coupling, both within and between components, with which regional and cross-regional frequency and amplitude modulations could be addressed.

Alpha band disruption in the AD-continuum starts in the Subjective Cognitive Decline stage: a MEG study

The consideration of Subjective Cognitive Decline (SCD) as a preclinical stage of AD remains still a matter of debate. Alpha band alterations represent one of the most significant changes in the electrophysiological profile of AD. In particular, AD patients exhibit reduced alpha relative power and frequency. We used alpha band activity measured with MEG to study whether SCD and MCI elders present these electrophysiological changes characteristic of AD, and to determine the evolution of the observed alterations across AD spectrum. The total sample consisted of 131 participants: 39 elders without SCD, 41 elders with SCD and 51 MCI patients. All of them underwent MEG and MRI scans and neuropsychological assessment. SCD and MCI patients exhibited a similar reduction in alpha band activity compared with the no SCD group. However, only MCI patients showed a slowing in their alpha peak frequency compared with both SCD and no SCD. These changes in alpha band were related to worse cognition. Our results suggest that AD-related alterations may start in the SCD stage, with a reduction in alpha relative power. It is later, in the MCI stage, where the slowing of the spectral profile takes place, giving rise to objective deficits in cognitive functioning.

Causal Role of Thalamic Interneurons in Brain State Transitions: A Study Using a Neural Mass Model Implementing Synaptic Kinetics

Experimental studies on the Lateral Geniculate Nucleus (LGN) of mammals and rodents show that the inhibitory interneurons (IN) receive around 47.1% of their afferents from the retinal spiking neurons, and constitute around 20–25% of the LGN cell population. However, there is a definite gap in knowledge about the role and impact of IN on thalamocortical dynamics in both experimental and model-based research. We use a neural mass computational model of the LGN with three neural populations viz. IN, thalamocortical relay (TCR), thalamic reticular nucleus (TRN), to study the causality of IN on LGN oscillations and state-transitions. The synaptic information transmission in the model is implemented with kinetic modeling, facilitating the linking of low-level cellular attributes with high-level population dynamics. The model is parameterized and tuned to simulate alpha (8–13 Hz) rhythm that is dominant in both Local Field Potential (LFP) of LGN and electroencephalogram (EEG) of visual cortex in an awake resting state with eyes closed. The results show that: First, the response of the TRN is suppressed in the presence of IN in the circuit; disconnecting the IN from the circuit effects a dramatic change in the model output, displaying high amplitude synchronous oscillations within the alpha band in both TCR and TRN. These observations conform to experimental reports implicating the IN as the primary inhibitory modulator of LGN dynamics in a cognitive state, and that reduced cognition is achieved by suppressing the TRN response. Second, the model validates steady state visually evoked potential response in humans corresponding to periodic input stimuli; however, when the IN is disconnected from the circuit, the output power spectra do not reflect the input frequency. This agrees with experimental reports underpinning the role of IN in efficient retino-geniculate information transmission. Third, a smooth transition from alpha to theta band is observed by progressive decrease of neurotransmitter concentrations in the synaptic clefts; however, the transition is abrupt with removal of the IN circuitry in the model. The results imply a role of IN toward maintaining homeostasis in the LGN by suppressing any instability that may arise due to anomalous synaptic attributes.

A neural mass model of phase-amplitude coupling

Brain activity shows phase-amplitude coupling between its slow and fast oscillatory components. We study phase-amplitude coupling as recorded at individual sites, using a modified version of the well-known Wendling neural mass model. To the population of fast inhibitory interneurons of this model, we added external modulatory input and dynamic self-feedback. These two modifications together are sufficient to let the inhibitory population serve as a limit-cycle oscillator, with frequency characteristics comparable to the beta and gamma bands. The frequency and power of these oscillations can be tuned through the time constant of the dynamic and modulatory input. Alpha band activity is generated, as is usual in such models, as a result of interactions of pyramidal neurons and a population of slow inhibitory interneurons. The slow inhibitory population activity directly influences the fast oscillations via the synaptic gain between slow and fast inhibitory populations. As a result, the amplitude envelope of the fast oscillation is coupled to the phase of the slow activity; this result is consistent with the notion that phase-amplitude coupling is effectuated by interactions between inhibitory interneurons.

Modeling the Generation of Phase-Amplitude Coupling in Cortical Circuits: From Detailed Networks to Neural Mass Models

Phase-amplitude coupling (PAC), the phenomenon where the amplitude of a high frequency oscillation is modulated by the phase of a lower frequency oscillation, is attracting an increasing interest in the neuroscience community due to its potential relevance for understanding healthy and pathological information processing in the brain. PAC is a diverse phenomenon, having been experimentally detected in at least ten combinations of rhythms: delta-theta, delta-alpha, delta-beta, delta-gamma, theta-alpha, theta-beta, theta-gamma, alpha-beta, alpha-gamma, and beta-gamma. However, a complete understanding of the biophysical mechanisms generating this diversity is lacking. Here we review computational models of PAC generation that range from detailed models of neuronal networks, where each cell is described by Hodgkin-Huxley-type equations, to neural mass models (NMMs) where only the average activities of neuronal populations are considered. We argue that NMMs are an appropriate mathematical framework (due to the small number of parameters and variables involved and the richness of the dynamics they can generate) to study the PAC phenomenon.

UKF-based closed loop iterative learning control of epileptiform wave in a neural mass model

A novel closed loop control framework is proposed to inhibit epileptiform wave in a neural mass model by external electric field, where the unscented Kalman filter method is used to reconstruct dynamics and estimate unmeasurable parameters of the model. Specifically speaking, the iterative learning control algorithm is introduced into the framework to optimize the control signal. In the proposed method, the control effect can be significantly improved based on the observation of the past attempts. Accordingly, the proposed method can effectively suppress the epileptiform wave as well as showing robustness to noises and uncertainties. Lastly, the simulation is carried out to illustrate the feasibility of the proposed method. Besides, this work shows potential value to design model-based feedback controllers for epilepsy treatment.