Measurement of the Mapping between Intracranial EEG and fMRI Recordings in the Human Brain

There are considerable gaps in our understanding of the relationship between human brain activity measured at different temporal and spatial scales. Here, electrocorticography (ECoG) measures were used to predict functional MRI changes in the sensorimotor cortex in two brain states: at rest and during motor performance. The specificity of this relationship to spatial co-localisation of the two signals was also investigated. We acquired simultaneous ECoG-fMRI in the sensorimotor cortex of three patients with epilepsy. During motor activity, high gamma power was the only frequency band where the electrophysiological response was co-localised with fMRI measures across all subjects. The best model of fMRI changes across states was its principal components, a parsimonious description of the entire ECoG spectrogram. This model performed much better than any others that were based either on the classical frequency bands or on summary measures of cross-spectral changes. The region-specific fMRI signal is reflected in spatially and spectrally distributed EEG activity.

Classification and characterisation of brain network changes in chronic back pain: A multicenter study

Background. Chronic pain is a common, often disabling condition thought to involve a combination of peripheral and central neurobiological factors. However, the extent and nature of changes in the brain is poorly understood.

Methods. We investigated brain network architecture using resting-state fMRI data in chronic back pain patients in the UK and Japan (41 patients, 56 controls), as well as open data from USA. We applied machine learning and deep learning (conditional variational autoencoder architecture) methods to explore classification of patients/controls based on network connectivity. We then studied the network topology of the data, and developed a multislice modularity method to look for consensus evidence of modular reorganisation in chronic back pain.

Results. Machine learning and deep learning allowed reliable classification of patients in a third, independent open data set with an accuracy of 63%, with 68% in cross validation of all data. We identified robust evidence of network hub disruption in chronic pain, most consistently with respect to clustering coefficient and betweenness centrality. We found a consensus pattern of modular reorganisation involving extensive, bilateral regions of sensorimotor cortex, and characterised primarily by negative reorganisation - a tendency for sensorimotor cortex nodes to be less inclined to form pairwise modular links with other brain nodes. Furthermore, these regions were found to display increased connectivity with the pregenual anterior cingulate cortex, a region known to be involved in endogenous pain control. In contrast, intraparietal sulcus displayed a propensity towards positive modular reorganisation, suggesting that it might have a role in forming modules associated with the chronic pain state.

Conclusion. The results provide evidence of consistent and characteristic brain network changes in chronic pain, characterised primarily by extensive reorganisation of the network architecture of the sensorimotor cortex.

Associations between polygenic risk scores for four psychiatric illnesses and brain structure using multivariate pattern recognition

Psychiatric illnesses are complex and polygenic. They are associated with widespread alterations in the brain, which are partly influenced by genetic factors. There have been some attempts to relate polygenic risk scores (PRS) - a measure of the overall genetic risk an individual carries for a disorder - to brain structure using univariate methods. However, PRS are likely associated with distributed and covarying effects across the brain. We therefore used multivariate machine learning in this proof-of-principle study to investigate associations between brain structure and PRS for four psychiatric disorders; attention deficit-hyperactivity disorder (ADHD), autism, bipolar disorder and schizophrenia. The sample included 213 individuals comprising patients with depression (69), bipolar disorder (33), and healthy controls (111). The five psychiatric PRSs were calculated based on summary data from the Psychiatric Genomics Consortium. T1-weighted magnetic resonance images were obtained and voxel-based morphometry was implemented in SPM12. Multivariate relevance vector regression was implemented in the Pattern Recognition for Neuroimaging Toolbox (PRoNTo). Across the whole sample, a multivariate pattern of grey matter significantly predicted the PRS for autism (r = 0.20, pFDR = 0.03; MSE = 4.20 × 10-5, pFDR = 0.02). For the schizophrenia PRS, the MSE was significant (MSE = 1.30 × 10-5, pFDR = 0.02) although the correlation was not (r = 0.15, pFDR = 0.06). These results lend support to the hypothesis that polygenic liability for autism and schizophrenia is associated with widespread changes in grey matter concentrations. These associations were seen in individuals not affected by these disorders, indicating that this is not driven by the expression of the disease, but by the genetic risk captured by the PRSs.

Classification and characterisation of brain network changes in chronic back pain: A multicenter study

Background. Chronic pain is a common, often disabling condition thought to involve a combination of peripheral and central neurobiological factors. However, the extent and nature of changes in the brain is poorly understood. Methods. We investigated brain network architecture using resting-state fMRI data in chronic back pain patients in the UK and Japan (41 patients, 56 controls), as well as open data from USA. We applied machine learning and deep learning (conditional variational autoencoder architecture) methods to explore classification of patients/controls based on network connectivity. We then studied the network topology of the data, and developed a multislice modularity method to look for consensus evidence of modular reorganisation in chronic back pain. Results. Machine learning and deep learning allowed reliable classification of patients in a third, independent open data set with an accuracy of 63%, with 68% in cross validation of all data. We identified robust evidence of network hub disruption in chronic pain, most consistently with respect to clustering coefficient and betweenness centrality. We found a consensus pattern of modular reorganisation involving extensive, bilateral regions of sensorimotor cortex, and characterised primarily by negative reorganisation - a tendency for sensorimotor cortex nodes to be less inclined to form pairwise modular links with other brain nodes. In contrast, intraparietal sulcus displayed a propensity towards positive modular reorganisation, suggesting that it might have a role in forming modules associated with the chronic pain state. Conclusion. The results provide evidence of consistent and characteristic brain network changes in chronic pain, characterised primarily by extensive reorganisation of the network architecture of the sensorimotor cortex.

Embedding Anatomical or Functional Knowledge in Whole-Brain Multiple Kernel Learning Models

Pattern recognition models have been increasingly applied to neuroimaging data over the last two decades. These applications have ranged from cognitive neuroscience to clinical problems. A common limitation of these approaches is that they do not incorporate previous knowledge about the brain structure and function into the models. Previous knowledge can be embedded into pattern recognition models by imposing a grouping structure based on anatomically or functionally defined brain regions. In this work, we present a novel approach that uses group sparsity to model the whole brain multivariate pattern as a combination of regional patterns. More specifically, we use a sparse version of Multiple Kernel Learning (MKL) to simultaneously learn the contribution of each brain region, previously defined by an atlas, to the decision function. Our application of MKL provides two beneficial features: (1) it can lead to improved overall generalisation performance when the grouping structure imposed by the atlas is consistent with the data; (2) it can identify a subset of relevant brain regions for the predictive model. In order to investigate the effect of the grouping in the proposed MKL approach we compared the results of three different atlases using three different datasets. The method has been implemented in the new version of the open-source Pattern Recognition for Neuroimaging Toolbox (PRoNTo).

Auditory prediction errors as individual biomarkers of schizophrenia

Schizophrenia is a complex psychiatric disorder, typically diagnosed through symptomatic evidence collected through patient interview. We aim to develop an objective biologically-based computational tool which aids diagnosis and relies on accessible imaging technologies such as electroencephalography (EEG). To achieve this, we used machine learning techniques and a combination of paradigms designed to elicit prediction errors or Mismatch Negativity (MMN) responses. MMN, an EEG component elicited by unpredictable changes in sequences of auditory stimuli, has previously been shown to be reduced in people with schizophrenia and this is arguably one of the most reproducible neurophysiological markers of schizophrenia.

EEG data were acquired from 21 patients with schizophrenia and 22 healthy controls whilst they listened to three auditory oddball paradigms comprising sequences of tones which deviated in 10% of trials from regularly occurring standard tones. Deviant tones shared the same properties as standard tones, except for one physical aspect: 1) duration - the deviant stimulus was twice the duration of the standard; 2) monaural gap - deviants had a silent interval omitted from the standard, or 3) inter-aural timing difference, which caused the deviant location to be perceived as 90° away from the standards.

We used multivariate pattern analysis, a machine learning technique implemented in the Pattern Recognition for Neuroimaging Toolbox (PRoNTo) to classify images generated through statistical parametric mapping (SPM) of spatiotemporal EEG data, i.e. event-related potentials measured on the two-dimensional surface of the scalp over time. Using support vector machine (SVM) and Gaussian processes classifiers (GPC), we were able classify individual patients and controls with balanced accuracies of up to 80.48% (p-values = 0.0326, FDR corrected) and an ROC analysis yielding an AUC of 0.87. Crucially, a GP regression revealed that MMN predicted global assessment of functioning (GAF) scores (correlation = 0.73, R² = 0.53, p = 0.0006).

•

The diagnostic utility of multiple auditory oddball stimulus paradigms is assessed.

•

Greatest classification accuracy was achieved using a monaural gap stimulus paradigm.

•

The full post-stimulus epoch contains relevant discriminatory components.

Performance indicators and indices of sludge management in urban wastewater treatment plants

Sludge (or biosolids) management is highly complex and has a significant cost associated with the biosolids disposal, as well as with the energy and flocculant consumption in the sludge processing units. The sludge management performance indicators (PIs) and indices (PXs) are thus core measures of the performance assessment system developed for urban wastewater treatment plants (WWTPs). The key PIs proposed cover the sludge unit production and dry solids concentration (DS), disposal/beneficial use, quality compliance for agricultural use and costs, whereas the complementary PIs assess the plant reliability and the chemical reagents' use. A key PI was also developed for assessing the phosphorus reclamation, namely through the beneficial use of the biosolids and the reclaimed water in agriculture. The results of a field study with 17 Portuguese urban WWTPs in a 5-year period were used to derive the PI reference values which are neither inherent to the PI formulation nor literature-based. Clusters by sludge type (primary, activated, trickling filter and mixed sludge) and by digestion and dewatering processes were analysed and the reference values for sludge production and dry solids were proposed for two clusters: activated sludge or biofilter WWTPs with primary sedimentation, sludge anaerobic digestion and centrifuge dewatering; activated sludge WWTPs without primary sedimentation and anaerobic digestion and with centrifuge dewatering. The key PXs are computed for the DS after each processing unit and the complementary PXs for the energy consumption and the operating conditions DS-determining. The PX reference values are treatment specific and literature based. The PI and PX system was applied to a WWTP and the results demonstrate that it diagnosis the situation and indicates opportunities and measures for improving the WWTP performance in sludge management.

A comprehensive approach for diagnosing opportunities for improving the performance of a WWTP

High quality services of wastewater treatment require a continuous assessment and improvement of the technical, environmental and economic performance. This paper demonstrates a comprehensive approach for benchmarking wastewater treatment plants (WWTPs), using performance indicators (PIs) and indices (PXs), in a 'plan-do-check-act' cycle routine driven by objectives. The performance objectives herein illustrated were to diagnose the effectiveness and energy performance of an oxidation ditch WWTP. The PI and PX results demonstrated an effective and reliable oxidation ditch (good-excellent performance), and a non-reliable UV disinfection (unsatisfactory-excellent performance) related with influent transmittance and total suspended solids. The energy performance increased with the treated wastewater volume and was unsatisfactory below 50% of plant capacity utilization. The oxidation ditch aeration performed unsatisfactorily and represented 38% of the plant energy consumption. The results allowed diagnosing opportunities for improving the energy and economic performance considering the influent flows, temperature and concentrations, and for levering the WWTP performance to acceptable-good effectiveness, reliability and energy efficiency. Regarding the plant reliability for fecal coliforms, improvement of UV lamp maintenance and optimization of the UV dose applied and microscreen recommissioning were suggested.

A NEURAL BIOMARKER FOR CHRONIC PAIN BASED ON DECODED BRAIN NETWORKS

The lack of a biomarker for chronic pain remains an important impediment to clinical and translational pain research. The problem stems from the multiple parallel but subtle abnormalties thought to represent the chronic pain state, yielding the emerging view of chronic pain as a ‘network disorder’. This suggests analysis approaches that aim to identify distributed patterns of data (multivariate, machine learning methods) might offer the best opportunity to discover biomarkers. Here, we performed a multi-center functional brain imaging study to record state functional brain networks resting in 41 patients with chronic back pain and 33 healthy control subjects. We calculated with functional covariance matrix from 160 regions of interest, and used Sparse Multinomial Logistic Regression to classify subjects as patient or control using a leave-one-out cross validation. Diagnostic accuracy was 91.9%, with sensitivity and specificity 90.2% and 93.9% respectively. We then used graph theoretic measures to characterise the pattern of network differences between the groups, and showed that the chronic pain state was associated with disrupted network ‘assortativity’. These data provide evidence to support an accurate functional biomarker of chronic pain, and open the door to the development of translatable biomarkers using similar methodologies in animals.

Connectivity changes underlying neurofeedback training of visual cortex activity

Neurofeedback based on real-time functional magnetic resonance imaging (fMRI) is a new approach that allows training of voluntary control over regionally specific brain activity. However, the neural basis of successful neurofeedback learning remains poorly understood. Here, we assessed changes in effective brain connectivity associated with neurofeedback training of visual cortex activity. Using dynamic causal modeling (DCM), we found that training participants to increase visual cortex activity was associated with increased effective connectivity between the visual cortex and the superior parietal lobe. Specifically, participants who learned to control activity in their visual cortex showed increased top-down control of the superior parietal lobe over the visual cortex, and at the same time reduced bottom-up processing. These results are consistent with efficient employment of top-down visual attention and imagery, which were the cognitive strategies used by participants to increase their visual cortex activity.

Connectivity-based neurofeedback: dynamic causal modeling for real-time fMRI

Neurofeedback based on real-time fMRI is an emerging technique that can be used to train voluntary control of brain activity. Such brain training has been shown to lead to behavioral effects that are specific to the functional role of the targeted brain area. However, real-time fMRI-based neurofeedback so far was limited to mainly training localized brain activity within a region of interest. Here, we overcome this limitation by presenting near real-time dynamic causal modeling in order to provide feedback information based on connectivity between brain areas rather than activity within a single brain area. Using a visual–spatial attention paradigm, we show that participants can voluntarily control a feedback signal that is based on the Bayesian model comparison between two predefined model alternatives, i.e. the connectivity between left visual cortex and left parietal cortex vs. the connectivity between right visual cortex and right parietal cortex. Our new approach thus allows for training voluntary control over specific functional brain networks. Because most mental functions and most neurological disorders are associated with network activity rather than with activity in a single brain region, this novel approach is an important methodological innovation in order to more directly target functionally relevant brain networks.

•

We adapt DCM for use in neurofeedback experiments.

•

Participants can control a DCM-based neurofeedback signal.

•

Real-time DCM allows for voluntary control over brain connectivity.

EEG-fMRI integration: a critical review of biophysical modeling and data analysis approaches

The diverse nature of cerebral activity, as measured using neuroimaging techniques, has been recognised long ago. It seems obvious that using single modality recordings can be limited when it comes to capturing its complex nature. Thus, it has been argued that moving to a multimodal approach will allow neuroscientists to better understand the dynamics and structure of this activity. This means that integrating information from different techniques, such as electroencephalography (EEG) and the blood oxygenated level dependent (BOLD) signal recorded with functional magnetic resonance imaging (fMRI), represents an important methodological challenge. In this work, we review the work that has been done thus far to derive EEG/fMRI integration approaches. This leads us to inspect the conditions under which such an integration approach could work or fail, and to disclose the types of scientific questions one could (and could not) hope to answer with it.

Time scales of representation in the human brain: weighing past information to predict future events

The estimates that humans make of statistical dependencies in the environment and therefore their representation of uncertainty crucially depend on the integration of data over time. As such, the extent to which past events are used to represent uncertainty has been postulated to vary over the cortex. For example, primary visual cortex responds to rapid perturbations in the environment, while frontal cortices involved in executive control encode the longer term contexts within which these perturbations occur. Here we tested whether primary and executive regions can be distinguished by the number of past observations they represent. This was based on a decay-dependent model that weights past observations from a Markov process and Bayesian Model Selection to test the prediction that neuronal responses are characterized by different decay half-lives depending on location in the brain. We show distributions of brain responses for short and long term decay functions in primary and secondary visual and frontal cortices, respectively. We found that visual and parietal responses are released from the burden of the past, enabling an agile response to fluctuations in events as they unfold. In contrast, frontal regions are more concerned with average trends over longer time scales within which local variations are embedded. Specifically, we provide evidence for a temporal gradient for representing context within the prefrontal cortex and possibly beyond to include primary sensory and association areas.

Bayesian model selection maps for group studies

This technical note describes the construction of posterior probability maps (PPMs) for Bayesian model selection (BMS) at the group level. This technique allows neuroimagers to make inferences about regionally specific effects using imaging data from a group of subjects. These effects are characterised using Bayesian model comparisons that are analogous to the F-tests used in statistical parametric mapping, with the advantage that the models to be compared do not need to be nested. Additionally, an arbitrary number of models can be compared together. This note describes the integration of the Bayesian mapping approach with a random effects analysis model for BMS using group data. We illustrate the method using fMRI data from a group of subjects performing a target detection task.

Estimating the transfer function from neuronal activity to BOLD using simultaneous EEG-fMRI

Previous studies using combined electrical and hemodynamic measurements of brain activity, such as EEG and (BOLD) fMRI, have yielded discrepant results regarding the relationship between neuronal activity and the associated BOLD response. In particular, some studies suggest that this link, or transfer function, depends on the frequency content of neuronal activity, while others suggest that total neuronal power accounts for the changes in BOLD. Here we explored this dependency by comparing different frequency-dependent and -independent transfer functions, using simultaneous EEG-fMRI. Our results suggest that changes in BOLD are indeed associated with changes in the spectral profile of neuronal activity and that these changes do not arise from one specific spectral band. Instead they result from the dynamics of the various frequency components together, in particular, from the relative power between high and low frequencies. Understanding the nature of the link between neuronal activity and BOLD plays a crucial role in improving the interpretability of BOLD images as well as on the design of more robust and realistic models for the integration of EEG and fMRI.

Storage proteins from Lathyrus sativus seeds

The proteins from Lathyrus sativus Linn. (chickling vetch or grass pea) seeds were investigated. Protein constitutes approximately 20% of the seed dry weight, >60% of which is composed by globulins and 30% by albumins. A single, 24 kDa polypeptide comprises more than half of the protein present in the albumin fraction. The globulins may be fractionated into three main components, which were named alpha-lathyrin (the major globulin), beta-lathyrin, and gamma-lathyrin. alpha-Lathyrin, with a sedimentation coefficient of approximately 18S, is composed of three main types of unglycosylated subunits (50-66 kDa), each of which produce, upon reduction, a heavy and a light polypeptide chain, by analogy with 11S. beta-Lathyrin, with a sedimentation coefficient of 13S, is composed by a relatively large number of subunits (8-66 kDa). Two major polypeptides are glycosylated and exhibit structural similarity with beta-conglutin from Lupinus albus. One of these possesses an internal disulfide bond. gamma-Lathyrin, with a sedimentation coefficient of approximately 5S, contains two interacting, unglycosylated polypeptides, with no disulfide bonds: the major 24 kDa albumin and the heavier (20 kDa) polypeptide chain of La. sativus lectin.

Diabetic retinopathy in a population of 1,302 insulin dependent diabetics (IDDM) diagnosed before 30 years of age

The authors analyse the data resulting from the first ophthalmological observation of 1,302 insulin dependent diabetics whose age at diagnosis is less than 30 years and who have been observed regularly by the Portuguese Diabetic Association. The prevalence of retinopathy is 41, 6%; 34.3% is non-proliferative and 7.3% is proliferative. Retinopathy is more frequent in males (P < 0.001). The prevalence of retinopathy increases with the duration of diabetes and it is equal to or greater than 80% in people who have had diabetes for 10 years or more. 'Poor' glucose control, the coexistence of other late complications and arterial hypertension increase the risk of retinopathy (P < 0.001).