scSTAR reveals hidden heterogeneity with a real-virtual cell pair structure across conditions in single-cell RNA sequencing data

Cell-state transition can reveal additional information from single-cell ribonucleic acid (RNA)-sequencing data in time-resolved biological phenomena. However, most of the current methods are based on the time derivative of the gene expression state, which restricts them to the short-term evolution of cell states. Here, we present single-cell State Transition Across-samples of RNA-seq data (scSTAR), which overcomes this limitation by constructing a paired-cell projection between biological conditions with an arbitrary time span by maximizing the covariance between two feature spaces using partial least square and minimum squared error methods. In mouse ageing data, the response to stress in CD4+ memory T cell subtypes was found to be associated with ageing. A novel Treg subtype characterized by mTORC activation was identified to be associated with antitumour immune suppression, which was confirmed by immunofluorescence microscopy and survival analysis in 11 cancers from The Cancer Genome Atlas Program. On melanoma data, scSTAR improved immunotherapy-response prediction accuracy from 0.8 to 0.96.

Degeneracy and stability in neural circuits of dopamine and serotonin neuromodulators: A theoretical consideration

Degenerate neural circuits perform the same function despite being structurally different. However, it is unclear whether neural circuits with interacting neuromodulator sources can themselves degenerate while maintaining the same neuromodulatory function. Here, we address this by computationally modeling the neural circuits of neuromodulators serotonin and dopamine, local glutamatergic and GABAergic interneurons, and their possible interactions, under reward/punishment-based conditioning tasks. The neural modeling is constrained by relevant experimental studies of the VTA or DRN system using, e.g., electrophysiology, optogenetics, and voltammetry. We first show that a single parsimonious, sparsely connected neural circuit model can recapitulate several separate experimental findings that indicated diverse, heterogeneous, distributed, and mixed DRNVTA neuronal signaling in reward and punishment tasks. The inability of this model to recapitulate all observed neuronal signaling suggests potentially multiple circuits acting in parallel. Then using computational simulations and dynamical systems analysis, we demonstrate that several different stable circuit architectures can produce the same observed network activity profile, hence demonstrating degeneracy. Due to the extensive D2-mediated connections in the investigated circuits, we simulate the D2 receptor agonist by increasing the connection strengths emanating from the VTA DA neurons. We found that the simulated D2 agonist can distinguish among sub-groups of the degenerate neural circuits based on substantial deviations in specific neural populations' activities in reward and punishment conditions. This forms a testable model prediction using pharmacological means. Overall, this theoretical work suggests the plausibility of degeneracy within neuromodulator circuitry and has important implications for the stable and robust maintenance of neuromodulatory functions.

Multi-dimensional relationships among dementia, depression and prescribed drugs in England and Wales hospitals

Background

 Dementia is a group of symptoms that largely affects older people. The majority of patients face behavioural and psychological symptoms (BPSD) during the course of their illness. Alzheimer’s disease (AD) and vascular dementia (VaD) are two of the most prevalent types of dementia. Available medications provide symptomatic benefits and provide relief from BPSD and associated health issues. However, it is unclear how specific dementia, antidepressant, antipsychotic, antianxiety, and mood stabiliser drugs, used in the treatment of depression and dementia subtypes are prescribed in hospital admission, during hospital stay, and at the time of discharge. To address this, we apply multi-dimensional data analytical approaches to understand drug prescribing practices within hospitals in England and Wales.

Methods

We made use of the UK National Audit of Dementia (NAD) dataset and pre-processed the dataset. We evaluated the pairwise Pearson correlation of the dataset and selected key data features which are highly correlated with dementia subtypes. After that, we selected drug prescribing behaviours (e.g. specific medications at the time of admission, during the hospital stay, and upon discharge), drugs and disorders. Then to shed light on the relations across multiple features or dimensions, we carried out multiple regression analyses, considering the number of dementia, antidepressant, antipsychotic, antianxiety, mood stabiliser, and antiepileptic/anticonvulsant drug prescriptions as dependent variables, and the prescription of other drugs, number of patients with dementia subtypes (AD/VaD), and depression as independent variables.

Results

In terms of antidepressant drugs prescribed in hospital admission, during stay and discharge, the number of sertraline and venlafaxine prescriptions were associated with the number of VaD patients whilst the number of mirtazapine prescriptions was associated with frontotemporal dementia patients. During admission, the number of lamotrigine prescriptions was associated with frontotemporal dementia patients, and with the number of valproate and dosulepin prescriptions. During discharge, the number of mirtazapine prescriptions was associated with the number of donepezil prescriptions in conjunction with frontotemporal dementia patients. Finally, the number of prescriptions of donepezil/memantine at admission, during hospital stay and at discharge exhibited positive association with AD patients.

Conclusion

Our analyses reveal a complex, multifaceted set of interactions among prescribed drug types, dementia subtypes, and depression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12911-022-01892-9.

Distinguishing Lewy Body Dementia from Alzheimer's Disease using Machine Learning on Heterogeneous Data: A Feasibility Study

Dementia with Lewy Bodies (DLB) is the second most common form of dementia, but diagnostic markers for DLB can be expensive and inaccessible, and many cases of DLB are undiagnosed. This work applies machine learning techniques to determine the feasibility of distinguishing DLB from Alzheimer's Disease (AD) using heterogeneous data features. The Repeated Incremental Pruning to Produce Error Reduction (RIPPER) algorithm was first applied using a Leave-One-Out Cross-Validation protocol to a dataset comprising DLB and AD cases. Then, interpretable association rule-based diagnostic classifiers were obtained for distinguishing DLB from AD. The various diagnostic classifiers generated by this process had high accuracy over the whole dataset (mean accuracy of 94%). The mean accuracy in classifying their out-of-sample case was 80.5%. Every classifier generated consisted of very simple structure, each using 1-2 classification rules and 1-3 data features. As a group, the classifiers were heterogeneous and used several different data features. In particular, some of the classifiers used very simple and inexpensive diagnostic features, yet with high diagnostic accuracy. This work suggests that opportunities may exist for incorporating accessible diagnostic assessments while improving diagnostic rate for DLB. Clinical Relevance- Simple and interpretable high-performing machine learning algorithms identified a variety of readily available clinical assessments for differential diagnosis of dementia offering the opportunities to incorporate various simple and inexpensive screening tests for DLB and addressing the problem of DLB underdiagnosis.

Multiple Cost Optimisation for Alzheimer's Disease Diagnosis

Current machine learning techniques for dementia diagnosis often do not take into account real-world practical constraints, which may include, for example, the cost of diagnostic assessment time and financial budgets. In this work, we built on previous cost-sensitive feature selection approaches by generalising to multiple cost types, while taking into consideration that stakeholders attempting to optimise the dementia care pathway might face multiple non-fungible budget constraints. Our new optimisation algorithm involved the searching of cost-weighting hyperparameters while constrained by total budgets. We then provided a proof of concept using both assessment time cost and financial budget cost. We showed that budget constraints could control the feature selection process in an intuitive and practical manner, while adjusting the hyperparameter increased the range of solutions selected by feature selection. We further showed that our budget-constrained cost optimisation framework could be implemented in a user-friendly graphical user interface sandbox tool to encourage non-technical users and stakeholders to adopt and to further explore and audit the model - a humans-in-the-loop approach. Overall, we suggest that setting budget constraints initially and then fine tuning the cost-weighting hyperparameters can be an effective way to perform feature selection where multiple cost constraints exist, which will in turn lead to more realistic optimising and redesigning of dementia diagnostic assessments. Clinical Relevance-By optimising diagnostic accuracy against various costs (e.g. assessment administration time and financial budget) predictive yet practical dementia diagnostic assessments can be redesigned to suit clinical use.

Integrated space-frequency-time domain feature extraction for MEG-based Alzheimer's disease classification

Magnetoencephalography (MEG) has been combined with machine learning techniques, to recognize the Alzheimer's disease (AD), one of the most common forms of dementia. However, most of the previous studies are limited to binary classification and do not fully utilize the two available MEG modalities (extracted using magnetometer and gradiometer sensors). AD consists of several stages of progression, this study addresses this limitation by using both magnetometer and gradiometer data to discriminate between participants with AD, AD-related mild cognitive impairment (MCI), and healthy control (HC) participants in the form of a three-class classification problem. A series of wavelet-based biomarkers are developed and evaluated, which concurrently leverage the spatial, frequency and time domain characteristics of the signal. A bimodal recognition system based on an improved score-level fusion approach is proposed to reinforce interpretation of the brain activity captured by magnetometers and gradiometers. In this preliminary study, it was found that the markers derived from gradiometer tend to outperform the magnetometer-based markers. Interestingly, out of the total 10 regions of interest, left-frontal lobe demonstrates about 8% higher mean recognition rate than the second-best performing region (left temporal lobe) for AD/MCI/HC classification. Among the four types of markers proposed in this work, the spatial marker developed using wavelet coefficients provided the best recognition performance for the three-way classification. Overall, the proposed approach provides promising results for the potential of AD/MCI/HC three-way classification utilizing the bimodal MEG data.

Detection of Mild Cognitive Impairment with MEG Functional Connectivity Using Wavelet-Based Neuromarkers

Studies on developing effective neuromarkers based on magnetoencephalographic (MEG) signals have been drawing increasing attention in the neuroscience community. This study explores the idea of using source-based magnitude-squared spectral coherence as a spatial indicator for effective regions of interest (ROIs) localization, subsequently discriminating the participants with mild cognitive impairment (MCI) from a group of age-matched healthy control (HC) elderly participants. We found that the cortical regions could be divided into two distinctive groups based on their coherence indices. Compared to HC, some ROIs showed increased connectivity (hyper-connected ROIs) for MCI participants, whereas the remaining ROIs demonstrated reduced connectivity (hypo-connected ROIs). Based on these findings, a series of wavelet-based source-level neuromarkers for MCI detection are proposed and explored, with respect to the two distinctive ROI groups. It was found that the neuromarkers extracted from the hyper-connected ROIs performed significantly better for MCI detection than those from the hypo-connected ROIs. The neuromarkers were classified using support vector machine (SVM) and k-NN classifiers and evaluated through Monte Carlo cross-validation. An average recognition rate of 93.83% was obtained using source-reconstructed signals from the hyper-connected ROI group. To better conform to clinical practice settings, a leave-one-out cross-validation (LOOCV) approach was also employed to ensure that the data for testing was from a participant that the classifier has never seen. Using LOOCV, we found the best average classification accuracy was reduced to 83.80% using the same set of neuromarkers obtained from the ROI group with functional hyper-connections. This performance surpassed the results reported using wavelet-based features by approximately 15%. Overall, our work suggests that (1) certain ROIs are particularly effective for MCI detection, especially when multi-resolution wavelet biomarkers are employed for such diagnosis; (2) there exists a significant performance difference in system evaluation between research-based experimental design and clinically accepted evaluation standards.

Practical Strategies for Extreme Missing Data Imputation in Dementia Diagnosis

Accurate computational models for clinical decision support systems require clean and reliable data but, in clinical practice, data are often incomplete. Hence, missing data could arise not only from training datasets but also test datasets which could consist of a single undiagnosed case, an individual. This work addresses the problem of extreme missingness in both training and test data by evaluating multiple imputation and classification workflows based on both diagnostic classification accuracy and computational cost. Extreme missingness is defined as having ~50% of the total data missing in more than half the data features. In particular, we focus on dementia diagnosis due to long time delays, high variability, high attrition rates and lack of practical data imputation strategies in its diagnostic pathway. We identified and replicated the extreme missingness structure of data from a real-world memory clinic on a larger open dataset, with the original complete data acting as ground truth. Overall, we found that computational cost, but not accuracy, varies widely for various imputation and classification approaches. Particularly, we found that iterative imputation on the training dataset combined with a reduced-feature classification model provides the best approach, in terms of speed and accuracy. Taken together, this work has elucidated important factors to be considered when developing a predictive model for a dementia diagnostic support system.

A Computational Study of Astrocytic GABA Release at the Glutamatergic Synapse: EAAT-2 and GAT-3 Coupled Dynamics

Neurotransmitter dynamics within neuronal synapses can be controlled by astrocytes and reflect key contributors to neuronal activity. In particular, Glutamate (Glu) released by activated neurons is predominantly removed from the synaptic space by perisynaptic astrocytic transporters EAAT-2 (GLT-1). In previous work, we showed that the time course of Glu transport is affected by ionic concentration gradients either side of the astrocytic membrane and has the propensity for influencing postsynaptic neuronal excitability. Experimental findings co-localize GABA transporters GAT-3 with EAAT-2 on the perisynaptic astrocytic membrane. While these transporters are unlikely to facilitate the uptake of synaptic GABA, this paper presents simulation results which demonstrate the coupling of EAAT-2 and GAT-3, giving rise to the ionic-dependent reversed transport of GAT-3. The resulting efflux of GABA from the astrocyte to the synaptic space reflects an important astrocytic mechanism for modulation of hyperexcitability. Key results also illustrate an astrocytic-mediated modulation of synaptic neuronal excitation by released GABA at the glutamatergic synapse.

An Early Stage Researcher's Primer on Systems Medicine Terminology

Background: Systems Medicine is a novel approach to medicine, that is, an interdisciplinary field that considers the human body as a system, composed of multiple parts and of complex relationships at multiple levels, and further integrated into an environment. Exploring Systems Medicine implies understanding and combining concepts coming from diametral different fields, including medicine, biology, statistics, modeling and simulation, and data science. Such heterogeneity leads to semantic issues, which may slow down implementation and fruitful interaction between these highly diverse fields. Methods: In this review, we collect and explain more than100 terms related to Systems Medicine. These include both modeling and data science terms and basic systems medicine terms, along with some synthetic definitions, examples of applications, and lists of relevant references. Results: This glossary aims at being a first aid kit for the Systems Medicine researcher facing an unfamiliar term, where he/she can get a first understanding of them, and, more importantly, examples and references for digging into the topic.

Shaping a data-driven era in dementia care pathway through computational neurology approaches

BACKGROUND

Dementia is caused by a variety of neurodegenerative diseases and is associated with a decline in memory and other cognitive abilities, while inflicting an enormous socioeconomic burden. The complexity of dementia and its associated comorbidities presents immense challenges for dementia research and care, particularly in clinical decision-making.

MAIN BODY

Despite the lack of disease-modifying therapies, there is an increasing and urgent need to make timely and accurate clinical decisions in dementia diagnosis and prognosis to allow appropriate care and treatment. However, the dementia care pathway is currently suboptimal. We propose that through computational approaches, understanding of dementia aetiology could be improved, and dementia assessments could be more standardised, objective and efficient. In particular, we suggest that these will involve appropriate data infrastructure, the use of data-driven computational neurology approaches and the development of practical clinical decision support systems. We also discuss the technical, structural, economic, political and policy-making challenges that accompany such implementations.

CONCLUSION

The data-driven era for dementia research has arrived with the potential to transform the healthcare system, creating a more efficient, transparent and personalised service for dementia.

Multi-time-point data preparation robustly reveals MCI and dementia risk factors

Introduction

Conflicting results on dementia risk factors have been reported across studies. We hypothesize that variation in data preparation methods may partially contribute to this issue.

Methods

We propose a comprehensive data preparation approach comparing individuals with stable diagnosis over time to those who progress to mild cognitive impairment (MCI)/dementia. This was compared to the often-used "baseline" analysis. Multivariate logistic regression was used to evaluate both methods.

Results

The results obtained from sensitivity analyses were consistent with those from our multi-time-point data preparation approach, exhibiting its robustness. Compared to analysis using only baseline data, the number of significant risk factors identified in progression analyses was substantially lower. Additionally, we found that moderate depression increased healthy-to-MCI/dementia risk, while hypertension reduced MCI-to-dementia risk.

Discussion

Overall, multi-time-point-based data preparation approaches may pave the way for a better understanding of dementia risk factors, and address some of the reproducibility issues in the field.

Optimality and Limitations of Audio-Visual Integration for Cognitive Systems

Multimodal integration is an important process in perceptual decision-making. In humans, this process has often been shown to be statistically optimal, or near optimal: sensory information is combined in a fashion that minimizes the average error in perceptual representation of stimuli. However, sometimes there are costs that come with the optimization, manifesting as illusory percepts. We review audio-visual facilitations and illusions that are products of multisensory integration, and the computational models that account for these phenomena. In particular, the same optimal computational model can lead to illusory percepts, and we suggest that more studies should be needed to detect and mitigate these illusions, as artifacts in artificial cognitive systems. We provide cautionary considerations when designing artificial cognitive systems with the view of avoiding such artifacts. Finally, we suggest avenues of research toward solutions to potential pitfalls in system design. We conclude that detailed understanding of multisensory integration and the mechanisms behind audio-visual illusions can benefit the design of artificial cognitive systems.

Opportunities for multiscale computational modelling of serotonergic drug effects in Alzheimer's disease

Alzheimer's disease (AD) is an age-specific neurodegenerative disease that compromises cognitive functioning and impacts the quality of life of an individual. Pathologically, AD is characterised by abnormal accumulation of beta-amyloid (Aβ) and hyperphosphorylated tau protein. Despite research advances over the last few decades, there is currently still no cure for AD. Although, medications are available to control some behavioural symptoms and slow the disease's progression, most prescribed medications are based on cholinesterase inhibitors. Over the last decade, there has been increased attention towards novel drugs, targeting alternative neurotransmitter pathways, particularly those targeting serotonergic (5-HT) system. In this review, we focused on 5-HT receptor (5-HTR) mediated signalling and drugs that target these receptors. These pathways regulate key proteins and kinases such as GSK-3 that are associated with abnormal levels of Aβ and tau in AD. We then review computational studies related to 5-HT signalling pathways with the potential for providing deeper understanding of AD pathologies. In particular, we suggest that multiscale and multilevel modelling approaches could potentially provide new insights into AD mechanisms, and towards discovering novel 5-HTR based therapeutic targets.

•

Alzheimer's disease (AD) drug treatment is limited, and alternatives are needed.

•

Serotonin (5-HT) mediated signalling pathways may regulate Aβ and tau levels.

•

5-HT based drugs have the potential to provide as novel therapeutics for AD.

•

Complex 5-HT signalling mechanisms for AD and related drugs hinder understanding.

•

Multiscale models may offer insights into mechanisms and therapeutic targets.

Changes-of-mind in the absence of new post-decision evidence

Decisions are occasionally accompanied by changes-of-mind. While considered a hallmark of cognitive flexibility, the mechanisms underlying changes-of-mind remain elusive. Previous studies on perceptual decision making have focused on changes-of-mind that are primarily driven by the accumulation of additional noisy sensory evidence after the initial decision. In a motion discrimination task, we demonstrate that changes-of-mind can occur even in the absence of additional evidence after the initial decision. Unlike previous studies of changes-of-mind, the majority of changes-of-mind in our experiment occurred in trials with prolonged initial response times. This suggests a distinct mechanism underlying such changes. Using a neural circuit model of decision uncertainty and change-of-mind behaviour, we demonstrate that this phenomenon is associated with top-down signals mediated by an uncertainty-monitoring neural population. Such a mechanism is consistent with recent neurophysiological evidence showing a link between changes-of-mind and elevated top-down neural activity. Our model explains the long response times associated with changes-of-mind through high decision uncertainty levels in such trials, and accounts for the observed motor response trajectories. Overall, our work provides a computational framework that explains changes-of-mind in the absence of new post-decision evidence.

A practical computerized decision support system for predicting the severity of Alzheimer's disease of an individual

Computerized clinical decision support systems can help to provide objective, standardized, and timely dementia diagnosis. However, current computerized systems are mainly based on group analysis, discrete classification of disease stages, or expensive and not readily accessible biomarkers, while current clinical practice relies relatively heavily on cognitive and functional assessments (CFA). In this study, we developed a computational framework using a suite of machine learning tools for identifying key markers in predicting the severity of Alzheimer's disease (AD) from a large set of biological and clinical measures. Six machine learning approaches, namely Kernel Ridge Regression (KRR), Support Vector Regression, and k-Nearest Neighbor for regression and Support Vector Machine (SVM), Random Forest, and k-Nearest Neighbor for classification, were used for the development of predictive models. We demonstrated high predictive power of CFA. Predictive performance of models incorporating CFA was shown to consistently have higher accuracy than those based solely on biomarker modalities. We found that KRR and SVM were the best performing regression and classification methods respectively. The optimal SVM performance was observed for a set of four CFA test scores (FAQ, ADAS13, MoCA, MMSE) with multi-class classification accuracy of 83.0%, 95%CI = (72.1%, 93.8%) while the best performance of the KRR model was reported with combined CFA and MRI neuroimaging data, i.e., R ² = 0.874, 95%CI = (0.827, 0.922). Given the high predictive power of CFA and their widespread use in clinical practice, we then designed a data-driven and self-adaptive computerized clinical decision support system (CDSS) prototype for evaluating the severity of AD of an individual on a continuous spectrum. The system implemented an automated computational approach for data pre-processing, modelling, and validation and used exclusively the scores of selected cognitive measures as data entries. Taken together, we have developed an objective and practical CDSS to aid AD diagnosis.

Calcium Microdomain Formation at the Perisynaptic Cradle Due to NCX Reversal: A Computational Study

It has recently been proposed using a multi-compartmental mathematical model that negatively fixed charged membrane-associated sites constrain the flow of cations in perisynaptic astroglial processes. This restricted movement of ions between the perisynaptic cradle (PsC), principal astroglial processes and the astrocyte soma gives rise to potassium (K+) and sodium (Na+) microdomains at the PsC. The present paper extends the above model to demonstrate that the formation of an Na+ microdomain can reverse the Na+/Ca2+ exchanger (NCX) thus providing an additional source of calcium (Ca2+) at the PsC. Results presented clearly show that reversal of the Na+/Ca2+ exchanger is instigated by a glutamate transporter coupled increase in concentration of cytoplasmic [Na+]i at the PsC, which and instigates Ca2+ influx through the NCX. As the flow of Ca2+ along the astrocyte process and away from the PsC is also constrained by Ca2+ binding proteins, then a Ca2+ microdomain forms at the PsC. The paper also serves to demonstrate that the EAAT, NKA, and NCX represent the minimal requirement necessary and sufficient for the development of a Ca2+ microdomain and that these mechanisms directly link neuronal activity and glutamate release to the formation of localized Na+ and Ca2+ microdomains signals at the PsC. This local source of Ca2+ can provide a previously underexplored form of astroglial Ca2+ signaling.

Bridging Neural and Computational Viewpoints on Perceptual Decision-Making

Sequential sampling models have provided a dominant theoretical framework guiding computational and neurophysiological investigations of perceptual decision-making. While these models share the basic principle that decisions are formed by accumulating sensory evidence to a bound, they come in many forms that can make similar predictions of choice behaviour despite invoking fundamentally different mechanisms. The identification of neural signals that reflect some of the core computations underpinning decision formation offers new avenues for empirically testing and refining key model assumptions. Here, we highlight recent efforts to explore these avenues and, in so doing, consider the conceptual and methodological challenges that arise when seeking to infer decision computations from complex neural data.

A hybrid computational approach for efficient Alzheimer's disease classification based on heterogeneous data

There is currently a lack of an efficient, objective and systemic approach towards the classification of Alzheimer's disease (AD), due to its complex etiology and pathogenesis. As AD is inherently dynamic, it is also not clear how the relationships among AD indicators vary over time. To address these issues, we propose a hybrid computational approach for AD classification and evaluate it on the heterogeneous longitudinal AIBL dataset. Specifically, using clinical dementia rating as an index of AD severity, the most important indicators (mini-mental state examination, logical memory recall, grey matter and cerebrospinal volumes from MRI and active voxels from PiB-PET brain scans, ApoE, and age) can be automatically identified from parallel data mining algorithms. In this work, Bayesian network modelling across different time points is used to identify and visualize time-varying relationships among the significant features, and importantly, in an efficient way using only coarse-grained data. Crucially, our approach suggests key data features and their appropriate combinations that are relevant for AD severity classification with high accuracy. Overall, our study provides insights into AD developments and demonstrates the potential of our approach in supporting efficient AD diagnosis.

A multimodal interface to resolve the Midas-Touch problem in gaze controlled wheelchair

Human-computer interaction (HCI) research has been playing an essential role in the field of rehabilitation. The usability of the gaze controlled powered wheelchair is limited due to Midas-Touch problem. In this work, we propose a multimodal graphical user interface (GUI) to control a powered wheelchair that aims to help upper-limb mobility impaired people in daily living activities. The GUI was designed to include a portable and low-cost eye-tracker and a soft-switch wherein the wheelchair can be controlled in three different ways: 1) with a touchpad 2) with an eye-tracker only, and 3) eye-tracker with soft-switch. The interface includes nine different commands (eight directions and stop) and integrated within a powered wheelchair system. We evaluated the performance of the multimodal interface in terms of lap-completion time, the number of commands, and the information transfer rate (ITR) with eight healthy participants. The analysis of the results showed that the eye-tracker with soft-switch provides superior performance with an ITR of 37.77 bits/min among the three different conditions (p<;0.05). Thus, the proposed system provides an effective and economical solution to the Midas-Touch problem and extended usability for the large population of disabled users.

Potassium and sodium microdomains in thin astroglial processes: A computational model study

A biophysical model that captures molecular homeostatic control of ions at the perisynaptic cradle (PsC) is of fundamental importance for understanding the interplay between astroglial and neuronal compartments. In this paper, we develop a multi-compartmental mathematical model which proposes a novel mechanism whereby the flow of cations in thin processes is restricted due to negatively charged membrane lipids which result in the formation of deep potential wells near the dipole heads. These wells restrict the flow of cations to “hopping” between adjacent wells as they transverse the process, and this surface retention of cations will be shown to give rise to the formation of potassium (K⁺) and sodium (Na⁺) microdomains at the PsC. We further propose that a K⁺ microdomain formed at the PsC, provides the driving force for the return of K⁺ to the extracellular space for uptake by the neurone, thereby preventing K⁺ undershoot. A slow decay of Na⁺ was also observed in our simulation after a period of glutamate stimulation which is in strong agreement with experimental observations. The pathological implications of microdomain formation during neuronal excitation are also discussed.

During periods of neuronal activity, ionic homeostasis in the surrounding extracellular space (ECS) is disturbed. To provide a healthy environment for continued neuronal function, excess ions such as potassium must be buffered away from the ECS; a vital supportive role provided by astrocyte cells. It has long been thought that astrocytes not only removed ions from the ECS but also transport them to other areas of the brain where their concentrations are lower. However, while our computational model simulations agree that astrocytes do remove these ions from the ECS they also show that these ions are mainly stored locally at the PsC to be returned to the ECS, thus restoring ionic homeostasis. Furthermore, we detail in this paper that this happens because of a previously overlooked biophysical phenomenon that is only dominant in thin astrocyte processes. The flow of these cations within thin processes is primarily by surface conduction where they experience the attraction of fixed negative charge at the membrane inner surface. This negative charge constrains cation movement along the surface and so their flow rate is restricted. Consequently, ions such as potassium that are released during neuronal excitation enter the PsC and are stored locally due to the low conductance pathway between the PsC and the astrocyte soma. Our simulations also show that this local build-up of K⁺ is returned to the ECS after the neuronal activity dies off which could potentially explain why K⁺ undershoot has not been observed; this result agrees with experimental observations. Moreover, the same mechanism can also explain the transient behaviour of Na⁺ ions whereby in thin processes a slow decay time constant is experimentally observed. These findings have important implications for the role of astrocytes in regulating neuronal excitability under physiological and pathological conditions, and therefore highlight the significance of the work presented in this paper.

A computational study of astrocytic glutamate influence on post-synaptic neuronal excitability

The ability of astrocytes to rapidly clear synaptic glutamate and purposefully release the excitatory transmitter is critical in the functioning of synapses and neuronal circuits. Dysfunctions of these homeostatic functions have been implicated in the pathology of brain disorders such as mesial temporal lobe epilepsy. However, the reasons for these dysfunctions are not clear from experimental data and computational models have been developed to provide further understanding of the implications of glutamate clearance from the extracellular space, as a result of EAAT2 downregulation: although they only partially account for the glutamate clearance process. In this work, we develop an explicit model of the astrocytic glutamate transporters, providing a more complete description of the glutamate chemical potential across the astrocytic membrane and its contribution to glutamate transporter driving force based on thermodynamic principles and experimental data. Analysis of our model demonstrates that increased astrocytic glutamate content due to glutamine synthetase downregulation also results in increased postsynaptic quantal size due to gliotransmission. Moreover, the proposed model demonstrates that increased astrocytic glutamate could prolong the time course of glutamate in the synaptic cleft and enhances astrocyte-induced slow inward currents, causing a disruption to the clarity of synaptic signalling and the occurrence of intervals of higher frequency postsynaptic firing. Overall, our work distilled the necessity of a low astrocytic glutamate concentration for reliable synaptic transmission of information and the possible implications of enhanced glutamate levels as in epilepsy.

The role of astrocytes in the excitability and hyperexcitability of neurons is a subject which has gained a lot of attention, particularly in the pathology of neurological disorders including epilepsy. Although not completely understood, the control of glutamate homeostasis is believed to play a role in paroxysmal neuronal hyperexcitability known to precede seizure activity. We have developed a computational model which explores two of the astrocytic homeostatic mechanisms, namely glutamate clearance and gliotransmission, and connect them with a common controlling factor, astrocytic cytoplasmic glutamate concentration. In our model simulations we demonstrate both a slower clearance rate of synaptic glutamate and enhanced astrocytic glutamate release where cytoplasmic glutamate is elevated, both of which contribute to high frequency neuronal firing and conditions for seizure generation. We also describe a viable role for astrocytes as a “high pass” filter, where astrocytic activation in the form of intracellular calcium oscillations is possible for only a certain range of presynaptic neuronal firing rates, the lower bound of the range being reduced where astrocytic glutamate is elevated. In physiological terms this perhaps indicates not only neuronal but also astrocytic glutamate-mediated excitation in the neural-astrocytic network.

Prescriptive variability of drugs by general practitioners

Prescription drug spending is growing faster than any other sector of healthcare. However, very little is known about patterns of prescribing and cost of prescribing between general practices. In this study, we examined variation in prescription rates and prescription costs through time for 55 GP surgeries in Northern Ireland Western Health and Social Care Trust. Temporal changes in variability of prescribing rates and costs were assessed using the Mann–Kendall test. Outlier practices contributing to between practice variation in prescribing rates were identified with the interquartile range outlier detection method. The relationship between rates and cost of prescribing was explored with Spearman's statistics. The differences in variability and mean number of prescribing rates associated with the practice setting and socioeconomic deprivation were tested using t-test and F-test respectively. The largest between-practice difference in prescribing rates was observed for Apr-Jun 2015, with the number of prescriptions ranging from 3.34 to 8.36 per patient. We showed that practices with outlier prescribing rates greatly contributed to between-practice variability. The largest difference in prescribing costs was reported for Apr-Jun 2014, with the prescription cost per patient ranging from £26.4 to £64.5. In addition, the temporal changes in variability of prescribing rates and costs were shown to undergo an upward trend. We demonstrated that practice setting and socio-economic deprivation accounted for some of the between-practice variation in prescribing. Rural practices had higher between practice variability than urban practices at all time points. Practices situated in more deprived areas had higher prescribing rates but lower variability than those located in less deprived areas. Further analysis is recommended to assess if variation in prescribing can be explained by demographic characteristics of patient population and practice features. Identification of other factors contributing to prescribing variability can help us better address potential inappropriateness of prescribing.

An integrated modelling framework for neural circuits with multiple neuromodulators

Neuromodulators are endogenous neurochemicals that regulate biophysical and biochemical processes, which control brain function and behaviour, and are often the targets of neuropharmacological drugs. Neuromodulator effects are generally complex partly owing to the involvement of broad innervation, co-release of neuromodulators, complex intra- and extrasynaptic mechanism, existence of multiple receptor subtypes and high interconnectivity within the brain. In this work, we propose an efficient yet sufficiently realistic computational neural modelling framework to study some of these complex behaviours. Specifically, we propose a novel dynamical neural circuit model that integrates the effective neuromodulator-induced currents based on various experimental data (e.g. electrophysiology, neuropharmacology and voltammetry). The model can incorporate multiple interacting brain regions, including neuromodulator sources, simulate efficiently and easily extendable to large-scale brain models, e.g. for neuroimaging purposes. As an example, we model a network of mutually interacting neural populations in the lateral hypothalamus, dorsal raphe nucleus and locus coeruleus, which are major sources of neuromodulator orexin/hypocretin, serotonin and norepinephrine/noradrenaline, respectively, and which play significant roles in regulating many physiological functions. We demonstrate that such a model can provide predictions of systemic drug effects of the popular antidepressants (e.g. reuptake inhibitors), neuromodulator antagonists or their combinations. Finally, we developed user-friendly graphical user interface software for model simulation and visualization for both fundamental sciences and pharmacological studies.

Toward a multiscale modeling framework for understanding serotonergic function

Despite its importance in regulating emotion and mental wellbeing, the complex structure and function of the serotonergic system present formidable challenges toward understanding its mechanisms. In this paper, we review studies investigating the interactions between serotonergic and related brain systems and their behavior at multiple scales, with a focus on biologically-based computational modeling. We first discuss serotonergic intracellular signaling and neuronal excitability, followed by neuronal circuit and systems levels. At each level of organization, we will discuss the experimental work accompanied by related computational modeling work. We then suggest that a multiscale modeling approach that integrates the various levels of neurobiological organization could potentially transform the way we understand the complex functions associated with serotonin.

Multi-Kernel Learning with Dartel Improves Combined MRI-PET Classification of Alzheimer's Disease in AIBL Data: Group and Individual Analyses

Magnetic resonance imaging (MRI) and positron emission tomography (PET) are neuroimaging modalities typically used for evaluating brain changes in Alzheimer’s disease (AD). Due to their complementary nature, their combination can provide more accurate AD diagnosis or prognosis. In this work, we apply a multi-modal imaging machine-learning framework to enhance AD classification and prediction of diagnosis of subject-matched gray matter MRI and Pittsburgh compound B (PiB)-PET data related to 58 AD, 108 mild cognitive impairment (MCI) and 120 healthy elderly (HE) subjects from the Australian imaging, biomarkers and lifestyle (AIBL) dataset. Specifically, we combined a Dartel algorithm to enhance anatomical registration with multi-kernel learning (MKL) technique, yielding an average of >95% accuracy for three binary classification problems: AD-vs.-HE, MCI-vs.-HE and AD-vs.-MCI, a considerable improvement from individual modality approach. Consistent with t-contrasts, the MKL weight maps revealed known brain regions associated with AD, i.e., (para)hippocampus, posterior cingulate cortex and bilateral temporal gyrus. Importantly, MKL regression analysis provided excellent predictions of diagnosis of individuals by r² = 0.86. In addition, we found significant correlations between the MKL classification and delayed memory recall scores with r² = 0.62 (p < 0.01). Interestingly, outliers in the regression model for diagnosis were mainly converter samples with a higher likelihood of converting to the inclined diagnostic category. Overall, our work demonstrates the successful application of MKL with Dartel on combined neuromarkers from different neuroimaging modalities in the AIBL data. This lends further support in favor of machine learning approach in improving the diagnosis and risk prediction of AD.

Directed Functional Connectivity in Fronto-Centroparietal Circuit Correlates With Motor Adaptation in Gait Training

Lower-extremity robotic exoskeletons are used in gait rehabilitation to achieve functional motor recovery. To date, little is known about how gait training and post-training are characterized in brain signals and their causal connectivity. In this work, we used time-domain partial Granger causality (PGC) analysis to elucidate the directed functional connectivity of electroencephalogram (EEG) signals of healthy adults in robot-assisted gait training (RAGT). Our results confirm the presence of EEG rhythms and corticomuscular relationships during standing and walking using spectral and coherence analyses. The PGC analysis revealed enhanced connectivity close to sensorimotor areas ( C₃ and CP₃ ) during standing, whereas additional connectivities involve the centroparietal ( CP _z) and frontal ( F _z ) areas during walking with respect to standing. In addition, significant fronto-centroparietal causal effects were found during both training and post-training. Strong correlations were also found between kinematic errors and fronto-centroparietal connectivity during training and post-training. This study suggests fronto-centroparietal connectivity as a potential neuromarker for motor learning and adaptation in RAGT.

Integrated dopaminergic neuronal model with reduced intracellular processes and inhibitory autoreceptors

Dopamine (DA) is an important neurotransmitter for multiple brain functions, and dysfunctions of the dopaminergic system are implicated in neurological and neuropsychiatric disorders. Although the dopaminergic system has been studied at multiple levels, an integrated and efficient computational model that bridges from molecular to neuronal circuit level is still lacking. In this study, the authors aim to develop a realistic yet efficient computational model of a dopaminergic pre-synaptic terminal. They first systematically perturb the variables/substrates of an established computational model of DA synthesis, release and uptake, and based on their relative dynamical timescales and steady-state changes, approximate and reduce the model into two versions: one for simulating hourly timescale, and another for millisecond timescale. They show that the original and reduced models exhibit rather similar steady and perturbed states, whereas the reduced models are more computationally efficient and illuminate the underlying key mechanisms. They then incorporate the reduced fast model into a spiking neuronal model that can realistically simulate the spiking behaviour of dopaminergic neurons. In addition, they successfully include autoreceptor-mediated inhibitory current explicitly in the neuronal model. This integrated computational model provides the first step toward an efficient computational platform for realistic multiscale simulation of dopaminergic systems in in silico neuropharmacology.

On self-feedback connectivity in neural mass models applied to event-related potentials

Neural mass models (NMMs) applied to neuroimaging data often do not emphasise intrinsic self-feedback within a neural population. However, based on mean-field theory, any population of coupled neurons is intrinsically endowed with effective self-coupling. In this work, we examine the effectiveness of three cortical NMMs with different self-feedbacks using a dynamic causal modelling approach. Specifically, we compare the classic Jansen and Rit (1995) model (no self-feedback), a modified model by Moran et al. (2007) (only inhibitory self-feedback), and our proposed model with inhibitory and excitatory self-feedbacks. Using bifurcation analysis, we show that single-unit Jansen-Rit model is less robust in generating oscillatory behaviour than the other two models. Next, under Bayesian inversion, we simulate single-channel event-related potentials (ERPs) within a mismatch negativity auditory oddball paradigm. We found fully self-feedback model (FSM) to provide the best fit to single-channel data. By analysing the posterior covariances of model parameters, we show that self-feedback connections are less sensitive to the generated evoked responses than the other model parameters, and hence can be treated analogously to "higher-order" parameter corrections of the original Jansen-Rit model. This is further supported in the more realistic multi-area case where FSM can replicate data better than JRM and MoM in the majority of subjects by capturing the finer features of the ERP data more accurately. Our work informs how NMMs with full self-feedback connectivity are not only more consistent with the underlying neurophysiology, but can also account for more complex features in ERP data.

Reduced computational models of serotonin synthesis, release, and reuptake

Multiscale computational models can provide systemic evaluation and prediction of neuropharmacological drug effects. To date, little computational modeling work has been done to bridge from intracellular to neuronal circuit level. A complex model that describes the intracellular dynamics of the presynaptic terminal of a serotonergic neuron has been previously proposed. By systematically perturbing the model's components, we identify the slow and fast dynamical components of the model, and the reduced slow or fast mode of the model is computationally significantly more efficient with accuracy not deviating much from the original model. The reduced fast-mode model is particularly suitable for incorporating into neurobiologically realistic spiking neuronal models, and hence for large-scale realistic computational simulations. We also develop user-friendly software based on the reduced models to allow scientists to rapidly test and predict neuropharmacological drug effects at a systems level.

Neural circuit interactions between the dorsal raphe nucleus and the lateral hypothalamus: an experimental and computational study

Orexinergic/hypocretinergic (Ox) neurotransmission plays an important role in regulating sleep, as well as in anxiety and depression, for which the serotonergic (5-HT) system is also involved in. However, little is known regarding the direct and indirect interactions between 5-HT in the dorsal raphe nucleus (DRN) and Ox neurons in the lateral hypothalamus (LHA). In this study, we report the additional presence of 5-HT1BR, 5-HT2AR, 5-HT2CR and fast ligand-gated 5-HT3AR subtypes on the Ox neurons of transgenic Ox-enhanced green fluorescent protein (Ox-EGFP) and wild type C57Bl/6 mice using single and double immunofluorescence (IF) staining, respectively, and quantify the colocalization for each 5-HT receptor subtype. We further reveal the presence of 5-HT3AR and 5-HT1AR on GABAergic neurons in LHA. We also identify NMDAR1, OX1R and OX2R on Ox neurons, but none on adjacent GABAergic neurons. This suggests a one-way relationship between LHA's GABAergic and Ox neurons, wherein GABAergic neurons exerts an inhibitory effect on Ox neurons under partial DRN's 5-HT control. We also show that Ox axonal projections receive glutamatergic (PSD-95 immunopositive) and GABAergic (Gephyrin immunopositive) inputs in the DRN. We consider these and other available findings into our computational model to explore possible effects of neural circuit connection types and timescales on the DRN-LHA system's dynamics. We find that if the connections from 5-HT to LHA's GABAergic neurons are weakly excitatory or inhibitory, the network exhibits slow oscillations; not observed when the connection is strongly excitatory. Furthermore, if Ox directly excites 5-HT neurons at a fast timescale, phasic Ox activation can lead to an increase in 5-HT activity; no significant effect with slower timescale. Overall, our experimental and computational approaches provide insights towards a more complete understanding of the complex relationship between 5-HT in the DRN and Ox in the LHA.

Directed neural connectivity changes in robot-assisted gait training: a partial Granger causality analysis

Now-a-days robotic exoskeletons are often used to help in gait training of stroke patients. However, such robotic systems have so far yielded only mixed results in benefiting the clinical population. Therefore, there is a need to investigate how gait learning and de-learning get characterised in brain signals and thus determine neural substrate to focus attention on, possibly, through an appropriate brain-computer interface (BCI). To this end, this paper reports the analysis of EEG data acquired from six healthy individuals undergoing robot-assisted gait training of a new gait pattern. Time-domain partial Granger causality (PGC) method was applied to estimate directed neural connectivity among relevant brain regions. To validate the results, a power spectral density (PSD) analysis was also performed. Results showed a strong causal interaction between lateral motor cortical areas. A frontoparietal connection was found in all robot-assisted training sessions. Following training, a causal "top-down" cognitive control was evidenced, which may indicate plasticity in the connectivity in the respective brain regions.

A mathematical model to explore the interdependence between the serotonin and orexin/hypocretin systems

Among their multitude of physiological and behavioral effects, the neurochemicals serotonin (5-HT) and orexin (Ox) have been closely linked to major depressive disorders (MDD) and sleep alterations. The dorsal raphe nucleus (DRN) and the lateral hypothalamus area (LHA) are brain regions that are sources of 5-HT and Ox, and there is evidence that suggests a reciprocal interaction between them. This lends support to the hypothesis of a close relationship between MDD and sleep disorders. Based on various experimental data, and appropriate assumptions, we construct a mathematical model of the coupled DRN-LHA neural circuit. Our model relates the dynamics of four important variables that can be experimentally measured: (i) the firing rate of 5-HT-containing neurons in DRN, (ii) the firing rate of Ox-containing neurons in the LHA, (iii) 5-HT concentration level in LHA, and (iv) Ox concentration level in DRN. Simulations show that our model supports the co-existence of baseline activities and concentration levels as observed in various separate experiments. It also allows circuit-level exploration of various parameters not yet identified experimentally, e.g. the rise and decay of Ox concentration levels due to Ox neural activity, and the exact dependence of Ox neural activity on 5-HT level. Finally we have made some model predictions regarding the effects of the 5-HT antagonist on the circuit. Our model, which can be subjected to verification and refinement as new experimental data accumulates, provides unified quantitative relationships and predictions between two important connected brain regions strongly tied to MDD and sleep disorders.

Can post-error dynamics explain sequential reaction time patterns?

We investigate human error dynamics in sequential two-alternative choice tasks. When subjects repeatedly discriminate between two stimuli, their error rates and reaction times (RTs) systematically depend on prior sequences of stimuli. We analyze these sequential effects on RTs, separating error and correct responses, and identify a sequential RT tradeoff: a sequence of stimuli which yields a relatively fast RT on error trials will produce a relatively slow RT on correct trials and vice versa. We reanalyze previous data and acquire and analyze new data in a choice task with stimulus sequences generated by a first-order Markov process having unequal probabilities of repetitions and alternations. We then show that relationships among these stimulus sequences and the corresponding RTs for correct trials, error trials, and averaged over all trials are significantly influenced by the probability of alternations; these relationships have not been captured by previous models. Finally, we show that simple, sequential updates to the initial condition and thresholds of a pure drift diffusion model can account for the trends in RT for correct and error trials. Our results suggest that error-based parameter adjustments are critical to modeling sequential effects.

Network properties of a computational model of the dorsal raphe nucleus

Serotonin (5-HT) plays an important role in regulating mood, cognition and behaviour. The midbrain dorsal raphe nucleus (DRN) is one of the primary sources of 5-HT. Recent studies show that DRN neuronal activities can encode rewarding (e.g., appetitive) and unrewarding (e.g., aversive) behaviours. Experiments have also shown that DRN neurons can exhibit heterogeneous spiking behaviours. In this work, we build and study a basic spiking neuronal network model of the DRN constrained by neuronal properties observed in experiments. We use an efficient adaptive quadratic integrate-and-fire neuronal model to capture slow afterhyperpolarization current, occasional bursting behaviours in 5-HT neurons, and fast spiking activities in the non-5-HT inhibitory neurons. Provided that our noisy and heterogeneous spiking neuronal network model adopts a feedforward inhibitory network architecture, it is able to replicate the main features of DRN neuronal activities recorded in monkeys performing a reward-based memory-guided saccade task. The model exhibits theta band oscillation, especially among the non-5-HT inhibitory neurons during the rewarding outcome of a simulated trial, thus forming a model prediction. By varying the inhibitory synaptic strengths and the afferent inputs, we find that the network model can oscillate over a range of relatively low frequencies, allow co-existence of multiple stable frequencies, and spike synchrony can spread from within a local neural subgroup to global. Our model suggests plausible network architecture, provides interesting model predictions that can be experimentally tested, and offers a sufficiently realistic multi-scale model for 5-HT neuromodulation simulations.

Dimension Reduction and Dynamics of a Spiking Neural Network Model for Decision Making under Neuromodulation()

Previous models of neuromodulation in cortical circuits have used either physiologically based networks of spiking neurons or simplified gain adjustments in low-dimensional connectionist models. Here we reduce a high-dimensional spiking neuronal network model, first to a four-population mean-field model and then to a two-population model. This provides a realistic implementation of neuromodulation in low-dimensional decision-making models, speeds up simulations by three orders of magnitude, and allows bifurcation and phase-plane analyses of the reduced models that illuminate neuromodulatory mechanisms. As modulation of excitation-inhibition varies, the network can move from unaroused states, through optimal performance to impulsive states, and eventually lose inhibition-driven winner-take-all behavior: all are clear outcomes of the bifurcation structure. We illustrate the value of reduced models by a study of the speed-accuracy tradeoff in decision making. The ability of such models to recreate neuromodulatory dynamics of the spiking network will accelerate the pace of future experiments linking behavioral data to cellular neurophysiology.