Lack of reproducibility of resting-state functional MRI findings in migraine with aura

BACKGROUND

Several studies have applied resting-state functional MRI to examine whether functional brain connectivity is altered in migraine with aura patients. These studies had multiple limitations, including small sample sizes, and reported conflicting results. Here, we performed a large, cross-sectional brain imaging study to reproduce previous findings.

METHODS

We recruited women aged 30-60 years from the nationwide Danish Twin Registry. Resting-state functional MRI of women with migraine with aura, their co-twins, and unrelated migraine-free twins was performed at a single centre. We carried out an extensive series of brain connectivity data analyses. Patients were compared to migraine-free controls and to co-twins.

RESULTS

Comparisons were based on data from 160 patients, 30 co-twins, and 136 controls. Patients were similar to controls with regard to age, and several lifestyle characteristics. We replicated clear effects of age on resting-state networks. In contrast, we failed to detect any differences, and to replicate previously reported differences, in functional connectivity between migraine patients with aura and non-migraine controls or their co-twins in any of the analyses.

CONCLUSION

Given the large sample size and the unbiased population-based design of our study, we conclude that women with migraine with aura have normal resting-state brain connectivity outside of migraine attacks.

Anterior cingulate glutamate levels associate with functional activation and connectivity during sensory integration in schizophrenia: a multimodal 1H-MRS and fMRI study

BACKGROUND

Glutamatergic dysfunction has been implicated in sensory integration deficits in schizophrenia, yet how glutamatergic function contributes to behavioural impairments and neural activities of sensory integration remains unknown.

METHODS

Fifty schizophrenia patients and 43 healthy controls completed behavioural assessments for sensory integration and underwent magnetic resonance spectroscopy (MRS) for measuring the anterior cingulate cortex (ACC) glutamate levels. The correlation between glutamate levels and behavioural sensory integration deficits was examined in each group. A subsample of 20 pairs of patients and controls further completed an audiovisual sensory integration functional magnetic resonance imaging (fMRI) task. Blood Oxygenation Level Dependent (BOLD) activation and task-dependent functional connectivity (FC) were assessed based on fMRI data. Full factorial analyses were performed to examine the Group-by-Glutamate Level interaction effects on fMRI measurements (group differences in correlation between glutamate levels and fMRI measurements) and the correlation between glutamate levels and fMRI measurements within each group.

RESULTS

We found that schizophrenia patients exhibited impaired sensory integration which was positively correlated with ACC glutamate levels. Multimodal analyses showed significantly Group-by-Glutamate Level interaction effects on BOLD activation as well as task-dependent FC in a 'cortico-subcortical-cortical' network (including medial frontal gyrus, precuneus, ACC, middle cingulate gyrus, thalamus and caudate) with positive correlations in patients and negative in controls.

CONCLUSIONS

Our findings indicate that ACC glutamate influences neural activities in a large-scale network during sensory integration, but the effects have opposite directionality between schizophrenia patients and healthy people. This implicates the crucial role of glutamatergic system in sensory integration processing in schizophrenia.

Optimized flip angle schemes for the split acquisition of fast spin-echo signals (SPLICE) sequence and application to diffusion-weighted imaging

PURPOSE

The diffusion-weighted SPLICE (split acquisition of fast spin-echo signals) sequence employs split-echo rapid acquisition with relaxation enhancement (RARE) readout to provide images almost free of geometric distortions. However, due to the varying T 2 $$ {}_2 $$ -weighting during k-space traversal, SPLICE suffers from blurring. This work extends a method for controlling the spatial point spread function (PSF) while optimizing the signal-to-noise ratio (SNR) achieved by adjusting the flip angles in the refocusing pulse train of SPLICE.

METHODS

An algorithm based on extended phase graph (EPG) simulations optimizes the flip angles by maximizing SNR for a flexibly chosen predefined target PSF that describes the desired k-space density weighting and spatial resolution. An optimized flip angle scheme and a corresponding post-processing correction filter which together achieve the target PSF was tested by healthy subject brain imaging using a clinical 1.5 T scanner.

RESULTS

Brain images showed a clear and consistent improvement over those obtained with a standard constant flip angle scheme. SNR was increased and apparent diffusion coefficient estimates were more accurate. For a modified Hann k-space weighting example, considerable benefits resulted from acquisition weighting by flip angle control.

CONCLUSION

The presented flexible method for optimizing SPLICE flip angle schemes offers improved MR image quality of geometrically accurate diffusion-weighted images that makes the sequence a strong candidate for radiotherapy planning or stereotactic surgery.

Decomposition-based framework for tumor classification and prediction of treatment response from longitudinal MRI

Objective.In the field of radiation oncology, the benefit of MRI goes beyond that of providing high soft-tissue contrast images for staging and treatment planning. With the recent clinical introduction of hybrid MRI linear accelerators it has become feasible to map physiological parameters describing diffusion, perfusion, and relaxation during the entire course of radiotherapy, for example. However, advanced data analysis tools are required for extracting qualified prognostic and predictive imaging biomarkers from longitudinal MRI data. In this study, we propose a new prediction framework tailored to exploit temporal dynamics of tissue features from repeated measurements. We demonstrate the framework using a newly developed decomposition method for tumor characterization.Approach.Two previously published MRI datasets with multiple measurements during and after radiotherapy, were used for development and testing:T₂-weighted multi-echo images obtained for two mouse models of pancreatic cancer, and diffusion-weighted images for patients with brain metastases. Initially, the data was decomposed using the novel monotonous slope non-negative matrix factorization (msNMF) tailored for MR data. The following processing consisted of a tumor heterogeneity assessment using descriptive statistical measures, robust linear modelling to capture temporal changes of these, and finally logistic regression analysis for stratification of tumors and volumetric outcome.Main Results.The framework was able to classify the two pancreatic tumor types with an area under curve (AUC) of 0.999,P< 0.001 and predict the tumor volume change with a correlation coefficient of 0.513,P= 0.034. A classification of the human brain metastases into responders and non-responders resulted in an AUC of 0.74,P= 0.065.Significance.A general data processing framework for analyses of longitudinal MRI data has been developed and applications were demonstrated by classification of tumor type and prediction of radiotherapy response. Further, as part of the assessment, the merits of msNMF for tumor tissue decomposition were demonstrated.

Altered empathy-related resting-state functional connectivity in patients with bipolar disorder

Empathy is the ability to generate emotional responses (i.e., cognitive empathy) and to make cognitive inferences (i.e., affective empathy) to other people's emotions. Empirical evidence suggests that patients with bipolar disorder (BD) exhibit impairment in cognitive empathy, but findings on affective empathy are inconsistent. Few studies have examined the neural mechanisms of cognitive and affective empathy in patients with BD. In this study, we examined the empathy-related resting-state functional connectivity (rsFC) in BD patients. Thirty-seven patients with BD and 42 healthy controls completed the self-report Questionnaires of Cognitive and Affective Empathy (QCAE), the Yoni behavioural task, and resting-sate fMRI brain scans. Group comparison of empathic ability was conducted. The interactions between group and empathic ability on seed-based whole brain rsFC were examined. BD patients scored lower on the Online Simulation subscale of the QCAE and showed positive correlations between cognitive empathy and the rsFC of the dorsal Medial Prefrontal Cortex (dmPFC) with the lingual gyrus. The correlations between cognitive empathy and the rsFC of the temporal-parietal junction (TPJ) with the fusiform gyrus, the cerebellum and the parahippocampus were weaker in BD patients than that in healthy controls. These findings highlight the underlying neural mechanisms of empathy impairments in BD patients.

Combining electro- and magnetoencephalography data using directional archetypal analysis

Metastable microstates in electro- and magnetoencephalographic (EEG and MEG) measurements are usually determined using modified k-means accounting for polarity invariant states. However, hard state assignment approaches assume that the brain traverses microstates in a discrete rather than continuous fashion. We present multimodal, multisubject directional archetypal analysis as a scale and polarity invariant extension to archetypal analysis using a loss function based on the Watson distribution. With this method, EEG/MEG microstates are modeled using subject- and modality-specific archetypes that are representative, distinct topographic maps between which the brain continuously traverses. Archetypes are specified as convex combinations of unit norm input data based on a shared generator matrix, thus assuming that the timing of neural responses to stimuli is consistent across subjects and modalities. The input data is reconstructed as convex combinations of archetypes using a subject- and modality-specific continuous archetypal mixing matrix. We showcase the model on synthetic data and an openly available face perception event-related potential data set with concurrently recorded EEG and MEG. In synthetic and unimodal experiments, we compare our model to conventional Euclidean multisubject archetypal analysis. We also contrast our model to a directional clustering model with discrete state assignments to highlight the advantages of modeling state trajectories rather than hard assignments. We find that our approach successfully models scale and polarity invariant data, such as microstates, accounting for intersubject and intermodal variability. The model is readily extendable to other modalities ensuring component correspondence while elucidating spatiotemporal signal variability.

Tracking of rigid head motion during MRI using an EEG system

Purpose

To demonstrate a novel method for tracking of head movements during MRI using electroencephalography (EEG) hardware for recording signals induced by native imaging gradients.

Theory and Methods

Gradient switching during simultaneous EEG–fMRI induces distortions in EEG signals, which depend on subject head position and orientation. When EEG electrodes are interconnected with high‐impedance carbon wire loops, the induced voltages are linear combinations of the temporal gradient waveform derivatives. We introduce head tracking based on these signals (CapTrack) involving 3 steps: (1) phantom scanning is used to characterize the target sequence and a fast calibration sequence; (2) a linear relation between changes of induced signals and head pose is established using the calibration sequence; and (3) induced signals recorded during target sequence scanning are used for tracking and retrospective correction of head movement without prolonging the scan time of the target sequence. Performance of CapTrack is compared directly to interleaved navigators.

Results

Head‐pose tracking at 27.5 Hz during echo planar imaging (EPI) was demonstrated with close resemblance to rigid body alignment (mean absolute difference: [0.14 0.38 0.15]‐mm translation, [0.30 0.27 0.22]‐degree rotation). Retrospective correction of 3D gradient‐echo imaging shows an increase of average edge strength of 12%/−0.39% for instructed/uninstructed motion with CapTrack pose estimates, with a tracking interval of 1561 ms and high similarity to interleaved navigator estimates (mean absolute difference: [0.13 0.33 0.12] mm, [0.28 0.15 0.22] degrees).

Conclusion

Motion can be estimated from recordings of gradient switching with little or no sequence modification, optionally in real time at low computational burden and synchronized to image acquisition, using EEG equipment already found at many research institutions.

Uncovering Cortical Units of Processing From Multi-Layered Connectomes

Modern diffusion and functional magnetic resonance imaging (dMRI/fMRI) provide non-invasive high-resolution images from which multi-layered networks of whole-brain structural and functional connectivity can be derived. Unfortunately, the lack of observed correspondence between the connectivity profiles of the two modalities challenges the understanding of the relationship between the functional and structural connectome. Rather than focusing on correspondence at the level of connections we presently investigate correspondence in terms of modular organization according to shared canonical processing units. We use a stochastic block-model (SBM) as a data-driven approach for clustering high-resolution multi-layer whole-brain connectivity networks and use prediction to quantify the extent to which a given clustering accounts for the connectome within a modality. The employed SBM assumes a single underlying parcellation exists across modalities whilst permitting each modality to possess an independent connectivity structure between parcels thereby imposing concurrent functional and structural units but different structural and functional connectivity profiles. We contrast the joint processing units to their modality specific counterparts and find that even though data-driven structural and functional parcellations exhibit substantial differences, attributed to modality specific biases, the joint model is able to achieve a consensus representation that well accounts for both the functional and structural connectome providing improved representations of functional connectivity compared to using functional data alone. This implies that a representation persists in the consensus model that is shared by the individual modalities. We find additional support for this viewpoint when the anatomical correspondence between modalities is removed from the joint modeling. The resultant drop in predictive performance is in general substantial, confirming that the anatomical correspondence of processing units is indeed present between the two modalities. Our findings illustrate how multi-modal integration admits consensus representations well-characterizing each individual modality despite their biases and points to the importance of multi-layered connectomes as providing supplementary information regarding the brain's canonical processing units.

A checklist for assessing the methodological quality of concurrent tES-fMRI studies (ContES checklist): a consensus study and statement

Low-intensity transcranial electrical stimulation (tES), including alternating or direct current stimulation, applies weak electrical stimulation to modulate the activity of brain circuits. Integration of tES with concurrent functional MRI (fMRI) allows for the mapping of neural activity during neuromodulation, supporting causal studies of both brain function and tES effects. Methodological aspects of tES-fMRI studies underpin the results, and reporting them in appropriate detail is required for reproducibility and interpretability. Despite the growing number of published reports, there are no consensus-based checklists for disclosing methodological details of concurrent tES-fMRI studies. The objective of this work was to develop a consensus-based checklist of reporting standards for concurrent tES-fMRI studies to support methodological rigor, transparency and reproducibility (ContES checklist). A two-phase Delphi consensus process was conducted by a steering committee (SC) of 13 members and 49 expert panelists through the International Network of the tES-fMRI Consortium. The process began with a circulation of a preliminary checklist of essential items and additional recommendations, developed by the SC on the basis of a systematic review of 57 concurrent tES-fMRI studies. Contributors were then invited to suggest revisions or additions to the initial checklist. After the revision phase, contributors rated the importance of the 17 essential items and 42 additional recommendations in the final checklist. The state of methodological transparency within the 57 reviewed concurrent tES-fMRI studies was then assessed by using the checklist. Experts refined the checklist through the revision and rating phases, leading to a checklist with three categories of essential items and additional recommendations: (i) technological factors, (ii) safety and noise tests and (iii) methodological factors. The level of reporting of checklist items varied among the 57 concurrent tES-fMRI papers, ranging from 24% to 76%. On average, 53% of checklist items were reported in a given article. In conclusion, use of the ContES checklist is expected to enhance the methodological reporting quality of future concurrent tES-fMRI studies and increase methodological transparency and reproducibility.

Locus Coeruleus Shows a Spatial Pattern of Structural Disintegration in Parkinson's Disease

BACKGROUND

Parkinson's disease (PD) causes a loss of neuromelanin-positive, noradrenergic neurons in the locus coeruleus (LC), which has been implicated in nonmotor dysfunction.

OBJECTIVES

We used "neuromelanin sensitive" magnetic resonance imaging (MRI) to localize structural disintegration in the LC and its association with nonmotor dysfunction in PD.

METHODS

A total of 42 patients with PD and 24 age-matched healthy volunteers underwent magnetization transfer weighted (MTw) MRI of the LC. The contrast-to-noise ratio of the MTw signal (CNR_MTw ) was used as an index of structural LC integrity. We performed slicewise and voxelwise analyses to map spatial patterns of structural disintegration, complemented by principal component analysis (PCA). We also tested for correlations between regional CNR_MTw and severity of nonmotor symptoms.

RESULTS

Mean CNR_MTw of the right LC was reduced in patients relative to controls. Voxelwise and slicewise analyses showed that the attenuation of CNR_MTw was confined to the right mid-caudal LC and linked regional CNR_MTw to nonmotor symptoms. CNR_MTw attenuation in the left mid-caudal LC was associated with the orthostatic drop in systolic blood pressure, whereas CNR_MTw attenuation in the caudal most portion of right LC correlated with apathy ratings. PCA identified a bilateral component that was more weakly expressed in patients. This component was characterized by a gradient in CNR_MTw along the rostro-caudal and dorso-ventral axes of the nucleus. The individual expression score of this component reflected the overall severity of nonmotor symptoms.

CONCLUSION

A spatially heterogeneous disintegration of LC in PD may determine the individual expression of specific nonmotor symptoms such as orthostatic dysregulation or apathy. © 2022 International Parkinson and Movement Disorder Society.

Data-driven separation of MRI signal components for tissue characterization

PURPOSE

MRI can be utilized for quantitative characterization of tissue. To assess e.g. water fractions or diffusion coefficients for compartments in the brain, a decomposition of the signal is necessary. Imposing standard models carries the risk of estimating biased parameters if model assumptions are violated. This work introduces a data-driven multicomponent analysis, the monotonous slope non-negative matrix factorization (msNMF), tailored to extract data features expected in MR signals.

METHODS

The msNMF was implemented by extending the standard NMF with monotonicity constraints on the signal profiles and their first derivatives. The method was validated using simulated data, and subsequently applied to both ex vivo DWI data and in vivo relaxometry data. Reproducibility of the method was tested using the latter.

RESULTS

The msNMF recovered the multi-exponential signals in the simulated data and showed superiority to standard NMF (based on the explained variance, area under the ROC curve, and coefficient of variation). Diffusion components extracted from the DWI data reflected the cell density of the underlying tissue. The relaxometry analysis resulted in estimates of edema water fractions (EWF) highly correlated with published results, and demonstrated acceptable reproducibility.

CONCLUSION

The msNMF can robustly separate MR signals into components with relation to the underlying tissue composition, and may potentially be useful for e.g. tumor tissue characterization.

Mapping cortico-subcortical sensitivity to 4 Hz amplitude modulation depth in human auditory system with functional MRI

Temporal modulations in the envelope of acoustic waveforms at rates around 4 Hz constitute a strong acoustic cue in speech and other natural sounds. It is often assumed that the ascending auditory pathway is increasingly sensitive to slow amplitude modulation (AM), but sensitivity to AM is typically considered separately for individual stages of the auditory system. Here, we used blood oxygen level dependent (BOLD) fMRI in twenty human subjects (10 male) to measure sensitivity of regional neural activity in the auditory system to 4 Hz temporal modulations. Participants were exposed to AM noise stimuli varying parametrically in modulation depth to characterize modulation-depth effects on BOLD responses. A Bayesian hierarchical modeling approach was used to model potentially nonlinear relations between AM depth and group-level BOLD responses in auditory regions of interest (ROIs). Sound stimulation activated the auditory brainstem and cortex structures in single subjects. BOLD responses to noise exposure in core and belt auditory cortices scaled positively with modulation depth. This finding was corroborated by whole-brain cluster-level inference. Sensitivity to AM depth variations was particularly pronounced in the Heschl's gyrus but also found in higher-order auditory cortical regions. None of the sound-responsive subcortical auditory structures showed a BOLD response profile that reflected the parametric variation in AM depth. The results are compatible with the notion that early auditory cortical regions play a key role in processing low-rate modulation content of sounds in the human auditory system.

Ergodicity-breaking reveals time optimal decision making in humans

Ergodicity describes an equivalence between the expectation value and the time average of observables. Applied to human behaviour, ergodic theories of decision-making reveal how individuals should tolerate risk in different environments. To optimize wealth over time, agents should adapt their utility function according to the dynamical setting they face. Linear utility is optimal for additive dynamics, whereas logarithmic utility is optimal for multiplicative dynamics. Whether humans approximate time optimal behavior across different dynamics is unknown. Here we compare the effects of additive versus multiplicative gamble dynamics on risky choice. We show that utility functions are modulated by gamble dynamics in ways not explained by prevailing decision theories. Instead, as predicted by time optimality, risk aversion increases under multiplicative dynamics, distributing close to the values that maximize the time average growth of in-game wealth. We suggest that our findings motivate a need for explicitly grounding theories of decision-making on ergodic considerations.

The effect of effort-reward imbalance on brain structure and resting-state functional connectivity in individuals with high levels of schizotypal traits

INTRODUCTION

Effort-reward imbalance (ERI) is a typical psychosocial stress. Schizotypal traits are attenuated features of schizophrenia in the general population. According to the diathesis-stress model, schizotypal traits and psychosocial stress contribute to the onset of schizophrenia. However, few studies examined the effects of these factors on brain alterations. This study aimed to examine relationships between ERI, schizotypal traits and brain structures and functions.

METHODS

We recruited 37 (13 male, 24 female) participants with high levels of schizotypal traits and 36 (12 male, 24 female) participants with low levels of schizotypal traits by the Schizotypal Personality Questionnaire (SPQ). The Chinese school version of the effort-reward imbalance questionnaire (C-ERI-S) was used to measure ERI. We conducted the voxel-based morphometry (VBM) and whole brain resting-state functional connectivity (rsFC) analysis using reward or stress-related regions as seeds.

RESULTS

Participants with high levels of schizotypal traits were more likely to perceive ERI. The severity of ERI was correlated with grey matter volume (GMV) reduction of the left pallidum and altered rsFC among the prefrontal, striatum and cerebellum in participants with high levels of schizotypal traits.

CONCLUSION

ERI is associated with GMV reduction and altered rsFC in individuals with high levels of schizotypal traits.

Concurrent TMS-fMRI for causal network perturbation and proof of target engagement

The experimental manipulation of neural activity by neurostimulation techniques overcomes the inherent limitations of correlative recordings, enabling the researcher to investigate causal brain-behavior relationships. But only when stimulation and recordings are combined, the direct impact of the stimulation on neural activity can be evaluated. In humans, this can be achieved non-invasively through the concurrent combination of transcranial magnetic stimulation (TMS) with functional magnetic resonance imaging (fMRI). Concurrent TMS-fMRI allows the assessment of the neurovascular responses evoked by TMS with excellent spatial resolution and full-brain coverage. This enables the functional mapping of both local and remote network effects of TMS in cortical as well as deep subcortical structures, offering unique opportunities for basic research and clinical applications. The purpose of this review is to introduce the reader to this powerful tool. We will introduce the technical challenges and state-of-the art solutions and provide a comprehensive overview of the existing literature and the available experimental approaches. We will highlight the unique insights that can be gained from concurrent TMS-fMRI, including the state-dependent assessment of neural responsiveness and inter-regional effective connectivity, the demonstration of functional target engagement, and the systematic evaluation of stimulation parameters. We will also discuss how concurrent TMS-fMRI during a behavioral task can help to link behavioral TMS effects to changes in neural network activity and to identify peripheral co-stimulation confounds. Finally, we will review the use of concurrent TMS-fMRI for developing TMS treatments of psychiatric and neurological disorders and suggest future improvements for further advancing the application of concurrent TMS-fMRI.

Reward signalling in brainstem nuclei under fluctuating blood glucose

Phasic dopamine release from mid-brain dopaminergic neurons is thought to signal errors of reward prediction (RPE). If reward maximisation is to maintain homeostasis, then the value of primary rewards should be coupled to the homeostatic errors they remediate. This leads to the prediction that RPE signals should be configured as a function of homeostatic state and thus diminish with the attenuation of homeostatic error. To test this hypothesis, we collected a large volume of functional MRI data from five human volunteers on four separate days. After fasting for 12 hours, subjects consumed preloads that differed in glucose concentration. Participants then underwent a Pavlovian cue-conditioning paradigm in which the colour of a fixation-cross was stochastically associated with the delivery of water or glucose via a gustometer. This design afforded computation of RPE separately for better- and worse-than expected outcomes during ascending and descending trajectories of serum glucose fluctuations. In the parabrachial nuclei, regional activity coding positive RPEs scaled positively with serum glucose for both ascending and descending glucose levels. The ventral tegmental area and substantia nigra became more sensitive to negative RPEs when glucose levels were ascending. Together, the results suggest that RPE signals in key brainstem structures are modulated by homeostatic trajectories of naturally occurring glycaemic flux, revealing a tight interplay between homeostatic state and the neural encoding of primary reward in the human brain.

Accurate and robust whole-head segmentation from magnetic resonance images for individualized head modeling

Transcranial brain stimulation (TBS) has been established as a method for modulating and mapping the function of the human brain, and as a potential treatment tool in several brain disorders. Typically, the stimulation is applied using a one-size-fits-all approach with predetermined locations for the electrodes, in electric stimulation (TES), or the coil, in magnetic stimulation (TMS), which disregards anatomical variability between individuals. However, the induced electric field distribution in the head largely depends on anatomical features implying the need for individually tailored stimulation protocols for focal dosing. This requires detailed models of the individual head anatomy, combined with electric field simulations, to find an optimal stimulation protocol for a given cortical target. Considering the anatomical and functional complexity of different brain disorders and pathologies, it is crucial to account for the anatomical variability in order to translate TBS from a research tool into a viable option for treatment. In this article we present a new method, called CHARM, for automated segmentation of fifteen different head tissues from magnetic resonance (MR) scans. The new method compares favorably to two freely available software tools on a five-tissue segmentation task, while obtaining reasonable segmentation accuracy over all fifteen tissues. The method automatically adapts to variability in the input scans and can thus be directly applied to clinical or research scans acquired with different scanners, sequences or settings. We show that an increase in automated segmentation accuracy results in a lower relative error in electric field simulations when compared to anatomical head models constructed from reference segmentations. However, also the improved segmentations and, by implication, the electric field simulations are affected by systematic artifacts in the input MR scans. As long as the artifacts are unaccounted for, this can lead to local simulation differences up to 30% of the peak field strength on reference simulations. Finally, we exemplarily demonstrate the effect of including all fifteen tissue classes in the field simulations against the standard approach of using only five tissue classes and show that for specific stimulation configurations the local differences can reach 10% of the peak field strength.

Two Coarse Spatial Patterns of Altered Brain Microstructure Predict Post-traumatic Amnesia in the Subacute Stage of Severe Traumatic Brain Injury

Introduction: Diffuse traumatic axonal injury (TAI) is one of the key mechanisms leading to impaired consciousness after severe traumatic brain injury (TBI). In addition, preferential regional expression of TAI in the brain may also influence clinical outcome. Aim: We addressed the question whether the regional expression of microstructural changes as revealed by whole-brain diffusion tensor imaging (DTI) in the subacute stage after severe TBI may predict the duration of post-traumatic amnesia (PTA). Method: Fourteen patients underwent whole-brain DTI in the subacute stage after severe TBI. Mean fractional anisotropy (FA) and mean diffusivity (MD) were calculated for five bilateral brain regions: fronto-temporal, parieto-occipital, and midsagittal hemispheric white matter, as well as brainstem and basal ganglia. Region-specific calculation of mean FA and MD only considered voxels that showed no tissue damage, using an exclusive mask with all voxels that belonged to local brain lesions or microbleeds. Mean FA or MD of the five brain regions were entered in separate partial least squares (PLS) regression analyses to identify patterns of regional microstructural changes that account for inter-individual variations in PTA. Results: For FA, PLS analysis revealed two spatial patterns that significantly correlated with individual PTA. The lower the mean FA values in all five brain regions, the longer that PTA lasted. A pattern characterized by lower FA values in the deeper brain regions relative to the FA values in the hemispheric regions also correlated with longer PTA. Similar trends were found for MD, but opposite in sign. The spatial FA changes as revealed by PLS components predicted the duration of PTA. Individual PTA duration, as predicted by a leave-one-out cross-validation analysis, correlated with true PTA values (Spearman r = 0.68, p permutation = 0.008). Conclusion: Two coarse spatial patterns of microstructural damage, indexed as reduction in FA, were relevant to recovery of consciousness after TBI. One pattern expressed was consistent with diffuse microstructural damage across the entire brain. A second pattern was indicative of a preferential damage of deep midline brain structures.

Discrete finger sequences are widely represented in human striatum

Research in primates and rodents ascribes the striatum a critical role in integrating elementary movements into unitary action sequences through reinforcement-based learning. Yet it remains to be shown whether the human striatum represents action sequence-specific information. Young right-handed volunteers underwent functional magnetic resonance imaging while they performed four discrete finger sequences with their right hand, consisting of five button presses. Specific finger sequences could be discriminated based on the distributed activity patterns in left and right striatum, but not by average differences in single-voxel activity. Multiple bilateral clusters in putamen and caudate nucleus belonging to motor, associative, parietal and limbic territories contributed to classification sensitivity. The results show that individual finger movement sequences are widely represented in human striatum, supporting functional integration rather than segregation. The findings are compatible with the idea that the basal ganglia simultaneously integrate motor, associative and limbic aspects in the control of complex overlearned behaviour.

Linking brain activity during sequential gambling to impulse control in Parkinson's disease

Dopaminergic treatment may impair the ability to suppress impulsive behaviours in patients with Parkinson's disease, triggering impulse control disorders. It is unclear how dopaminergic medication affects the neural networks that contribute to withholding inappropriate actions. To address this question, we mapped task-related brain activity with whole-brain functional magnetic resonance imaging at 3 Tesla in 26 patients with Parkinson's disease. Patients performed a sequential gambling task while being ON and OFF their regular dopaminergic treatment. During a gambling round, patients repeatedly decided between the option to continue with gambling and accumulate more monetary reward under increasing risk or the option to bank the current balance and start a new round. 13 patients had an impulse control disorder (ICD + group). These patients did not differ in risk-taking attitude during sequential gambling from 13 patients without impulse control disorder (ICD - group), but they displayed differences in gambling-related activity in cortico-subcortical brain areas supporting inhibitory control. First, the ICD + group showed reduced "continue-to-gamble" activity in right inferior frontal gyrus and subthalamic nucleus. Second, the individual risk-attitude scaled positively with "continue-to-gamble" activity in right subthalamic nucleus and striatum in the ICD - group only. Third, ICD + patients differed in their functional neural responses to dopaminergic treatment from ICD - patients: dopaminergic therapy reduced functional connectivity between inferior frontal gyrus and subthalamic nucleus during "continue-to-gamble" decisions and attenuated striatal responses towards accumulating reward and risk. Together, the medication-independent (trait) and medication-related (state) differences in neural activity may set a permissive stage for the emergence of impulse control disorders during dopamine replacement therapy in Parkinson's disease.

Processing of Positive Visual Stimuli Before and After Symptoms Provocation in Posttraumatic Stress Disorder - A Functional Magnetic Resonance Imaging Study of Trauma-Affected Male Refugees

Background

Symptoms of anhedonia are often central to posttraumatic stress disorder (PTSD), but it is unclear how anhedonia is affected by processes induced by reliving past traumatic memories.

Methods

Sixty-nine male refugees (PTSD = 38) were interviewed and scanned with functional magnetic resonance imaging while viewing positive, neutral and Scrambled Pictures after being read personalized scripts evoking an emotionally neutral memory and a traumatic memory. We further measured postprovocation state symptoms, physiological measures and PTSD symptoms. We tested whether neural activity associated with positive picture viewing in participants with PTSD was differentially affected by symptom provocation compared to controls.

Results

For the pictures > scrambled contrast (Positive contrast), PTSD participants had significantly less activity than controls in fusiform gyrus, right inferior temporal gyrus and left middle occipital gyrus. The Positive contrast activity in fusiform gyrus scaled negatively with anhedonia symptoms in PTSD participants after controlling for total PTSD severity. Relative to the emotionally Neutral Script, the Trauma Script decreased positive picture viewing activity in posterior cingulate cortex, precuneus and left calcarine gyrus, but there was no difference between PTSD participants and controls.

Conclusions

We found reduced responsiveness of higher visual processing of emotionally positive pictures in PTSD. The significant correlation found between positive picture viewing activity and anhedonia suggests the reduced responsiveness to be due to the severity of anhedonia.

Generalizability of machine learning for classification of schizophrenia based on resting-state functional MRI data

Machine learning has increasingly been applied to classification of schizophrenia in neuroimaging research. However, direct replication studies and studies seeking to investigate generalizability are scarce. To address these issues, we assessed within-site and between-site generalizability of a machine learning classification framework which achieved excellent performance in a previous study using two independent resting-state functional magnetic resonance imaging data sets collected from different sites and scanners. We established within-site generalizability of the classification framework in the main data set using cross-validation. Then, we trained a model in the main data set and investigated between-site generalization in the validated data set using external validation. Finally, recognizing the poor between-site generalization performance, we updated the unsupervised algorithm to investigate if transfer learning using additional unlabeled data were able to improve between-site classification performance. Cross-validation showed that the published classification procedure achieved an accuracy of 0.73 using majority voting across all selected components. External validation found a classification accuracy of 0.55 (not significant) and 0.70 (significant) using the direct and transfer learning procedures, respectively. The failure of direct generalization from one site to another demonstrates the limitation of within-site cross-validation and points toward the need to incorporate efforts to facilitate application of machine learning across multiple data sets. The improvement in performance with transfer learning highlights the importance of taking into account the properties of data when constructing predictive models across samples and sites. Our findings suggest that machine learning classification result based on a single study should be interpreted cautiously.

Classification of social anhedonia using temporal and spatial network features from a social cognition fMRI task

Previous studies have suggested that the degree of social anhedonia reflects the vulnerability for developing schizophrenia. However, only few studies have investigated how functional network changes are related to social anhedonia. The aim of this fMRI study was to classify subjects according to their degree of social anhedonia using supervised machine learning. More specifically, we extracted both spatial and temporal network features during a social cognition task from 70 subjects, and used support vector machines for classification. Since impairment in social cognition is well established in schizophrenia-spectrum disorders, the subjects performed a comic strip task designed to specifically probe theory of mind (ToM) and empathy processing. Features representing both temporal (time series) and network dynamics were extracted using task activation maps, seed region analysis, independent component analysis (ICA), and a newly developed multi-subject archetypal analysis (MSAA), which here aimed to further bridge aspects of both seed region analysis and decomposition by incorporating a spotlight approach.We found significant classification of subjects with elevated levels of social anhedonia when using the times series extracted using MSAA, indicating that temporal dynamics carry important information for classification of social anhedonia. Interestingly, we found that the same time series yielded the highest classification performance in a task classification of the ToM condition. Finally, the spatial network corresponding to that time series included both prefrontal and temporal-parietal regions as well as insula activity, which previously have been related schizotypy and the development of schizophrenia.

Limited Colocalization of Microbleeds and Microstructural Changes after Severe Traumatic Brain Injury

Severe traumatic brain injury (TBI) produces shearing forces on long-range axons and brain vessels, causing axonal and vascular injury. To examine whether microbleeds and axonal injury colocalize after TBI, we performed whole-brain susceptibility-weighted imaging (SWI) and diffusion tensor imaging (DTI) in 14 patients during the subacute phase after severe TBI. SWI was used to determine the number and volumes of microbleeds in five brain regions: the frontotemporal lobe; parieto-occipital lobe; midsagittal region (cingular cortex, parasagittal white matter, and corpus callosum); deep nuclei (basal ganglia and thalamus); and brainstem. Averaged fractional anisotropy (FA) and mean diffusivity (MD) were measured to assess microstructural changes in the normal appearing white matter attributed to axonal injury in the same five regions. Regional expressions of microbleeds and microstructure were used in a partial least-squares model to predict the impairment of consciousness in the subacute stage after TBI as measured with the Coma Recovery Scale-Revised (CRS-R). Only in the midsagittal region, the expression of microbleeds was correlated with regional changes in microstructure as revealed by DTI. Microbleeds and microstructural DTI-based metrics of deep, but not superficial, brain regions were able to predict individual CRS-R. Our results suggest that microbleeds are not strictly related to axonal pathology in other than the midsagittal region. While each measure alone was predictive, the combination of both metrics scaled best with individual CRS-R. Structural alterations in deep brain structures are relevant in terms of determining the severity of impaired consciousness in the acute stage after TBI.

Whole-Brain Exploratory Analysis of Functional Task Response Following Erythropoietin Treatment in Mood Disorders: A Supervised Machine Learning Approach

A core symptom of mood disorders is cognitive impairment in attention, memory and executive functions. Erythropoietin (EPO) is a candidate treatment for cognitive impairment in unipolar and bipolar disorders (UD and BD) and modulates cognition-related neural activity across a fronto-temporo-parietal network. This report investigates predicting the pharmacological treatment from functional magnetic resonance imaging (fMRI) data using a supervised machine learning approach. A total of 84 patients with UD or BD were included in a randomized double-blind parallel-group study in which they received eight weekly infusions of either EPO (40 000 IU) or saline. Task fMRI data were collected before EPO/saline infusions started (baseline) and 6 weeks after last infusion (follow-up). During the scanning sessions, participants were given an n-back working memory and a picture encoding task. Linear classification models with different regularization techniques were used to predict treatment status from both cross-sectional data (at follow-up) and longitudinal data (difference between baseline and follow-up). For the n-back and picture encoding tasks, data were available and analyzed for 52 (EPO; n = 28, Saline; n = 24) and 59 patients (EPO; n = 31, Saline; n = 28), respectively. We found limited evidence that the classifiers used could predict treatment status at a reliable level of performance (≤60% accuracy) when tested using repeated cross-validation. There was no difference in using cross-sectional versus longitudinal data. Whole-brain multivariate decoding applied to pharmaco-fMRI in small to moderate samples seems to be suboptimal for exploring data driven neuronal treatment mechanisms.

Electric field simulations for transcranial brain stimulation using FEM: an efficient implementation and error analysis

OBJECTIVE

Transcranial magnetic stimulation (TMS) and transcranial electric stimulation (TES) modulate brain activity non-invasively by generating electric fields either by electromagnetic induction or by injecting currents via skin electrodes. Numerical simulations based on anatomically detailed head models of the TMS and TES electric fields can help us to understand and optimize the spatial stimulation pattern in the brain. However, most realistic simulations are still slow, and the role of anatomical fidelity on simulation accuracy has not been evaluated in detail so far.

APPROACH

We present and validate a new implementation of the finite element method (FEM) for TMS and TES that is based on modern algorithms and libraries. We also evaluate the convergence of the simulations and estimate errors stemming from numerical and modelling aspects.

MAIN RESULTS

Comparisons with analytical solutions for spherical phantoms validate our new FEM implementation, which is three to six times faster than previous implementations. The convergence results suggest that accurately capturing the tissue geometry in addition to choosing a sufficiently accurate numerical method is of fundamental importance for accurate simulations.

SIGNIFICANCE

The new implementation allows for a substantial increase in computational efficiency of FEM TMS and TES simulations. This is especially relevant for applications such as the systematic assessment of model uncertainty and the optimization of multi-electrode TES montages. The results of our systematic error analysis allow the user to select the best tradeoff between model resolution and simulation speed for a specific application. The new FEM code is openly available as a part of our open-source software SimNIBS 3.0.

A principled approach to conductivity uncertainty analysis in electric field calculations

Uncertainty surrounding ohmic tissue conductivity impedes accurate calculation of the electric fields generated by non-invasive brain stimulation. We present an efficient and generic technique for uncertainty and sensitivity analyses, which quantifies the reliability of field estimates and identifies the most influential parameters. For this purpose, we employ a non-intrusive generalized polynomial chaos expansion to compactly approximate the multidimensional dependency of the field on the conductivities. We demonstrate that the proposed pipeline yields detailed insight into the uncertainty of field estimates for transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS), identifies the most relevant tissue conductivities, and highlights characteristic differences between stimulation methods. Specifically, we test the influence of conductivity variations on (i) the magnitude of the electric field generated at each gray matter location, (ii) its normal component relative to the cortical sheet, (iii) its overall magnitude (indexed by the 98th percentile), and (iv) its overall spatial distribution. We show that TMS fields are generally less affected by conductivity variations than tDCS fields. For both TMS and tDCS, conductivity uncertainty causes much higher uncertainty in the magnitude as compared to the direction and overall spatial distribution of the electric field. Whereas the TMS fields were predominantly influenced by gray and white matter conductivity, the tDCS fields were additionally dependent on skull and scalp conductivities. Comprehensive uncertainty analyses of complex systems achieved by the proposed technique are not possible with classical methods, such as Monte Carlo sampling, without extreme computational effort. In addition, our method has the advantages of directly yielding interpretable and intuitive output metrics and of being easily adaptable to new problems.

Perspectives on Machine Learning for Classification of Schizotypy Using fMRI Data

Functional magnetic resonance imaging is capable of estimating functional activation and connectivity in the human brain, and lately there has been increased interest in the use of these functional modalities combined with machine learning for identification of psychiatric traits. While these methods bear great potential for early diagnosis and better understanding of disease processes, there are wide ranges of processing choices and pitfalls that may severely hamper interpretation and generalization performance unless carefully considered. In this perspective article, we aim to motivate the use of machine learning schizotypy research. To this end, we describe common data processing steps while commenting on best practices and procedures. First, we introduce the important role of schizotypy to motivate the importance of reliable classification, and summarize existing machine learning literature on schizotypy. Then, we describe procedures for extraction of features based on fMRI data, including statistical parametric mapping, parcellation, complex network analysis, and decomposition methods, as well as classification with a special focus on support vector classification and deep learning. We provide more detailed descriptions and software as supplementary material. Finally, we present current challenges in machine learning for classification of schizotypy and comment on future trends and perspectives.

Task-Modulated Cortical Representations of Natural Sound Source Categories

In everyday sound environments, we recognize sound sources and events by attending to relevant aspects of an acoustic input. Evidence about the cortical mechanisms involved in extracting relevant category information from natural sounds is, however, limited to speech. Here, we used functional MRI to measure cortical response patterns while human listeners categorized real-world sounds created by objects of different solid materials (glass, metal, wood) manipulated by different sound-producing actions (striking, rattling, dropping). In different sessions, subjects had to identify either material or action categories in the same sound stimuli. The sound-producing action and the material of the sound source could be decoded from multivoxel activity patterns in auditory cortex, including Heschl's gyrus and planum temporale. Importantly, decoding success depended on task relevance and category discriminability. Action categories were more accurately decoded in auditory cortex when subjects identified action information. Conversely, the material of the same sound sources was decoded with higher accuracy in the inferior frontal cortex during material identification. Representational similarity analyses indicated that both early and higher-order auditory cortex selectively enhanced spectrotemporal features relevant to the target category. Together, the results indicate a cortical selection mechanism that favors task-relevant information in the processing of nonvocal sound categories.

Working Memory Modulation of Frontoparietal Network Connectivity in First-Episode Schizophrenia

Working memory (WM) impairment is regarded as a core aspect of schizophrenia. However, the neural mechanisms behind this cognitive deficit remain unclear. The connectivity of a frontoparietal network is known to be important for subserving WM. Using functional magnetic resonance imaging, the current study investigated whether WM-dependent modulation of effective connectivity in this network is affected in a group of first-episode schizophrenia (FES) patients compared with similarly performing healthy participants during a verbal n-back task. Dynamic causal modeling (DCM) of the coupling between regions (left inferior frontal gyrus (IFG), left inferior parietal lobe (IPL), and primary visual area) identified in a psychophysiological interaction (PPI) analysis was performed to characterize effective connectivity during the n-back task. The PPI analysis revealed that the connectivity between the left IFG and left IPL was modulated by WM and that this modulation was reduced in FES patients. The subsequent DCM analysis confirmed this modulation by WM and found evidence that FES patients had reduced forward connectivity from IPL to IFG. These findings provide evidence for impaired WM modulation of frontoparietal effective connectivity in the early phase of schizophrenia, even with intact WM performance, suggesting a failure of context-sensitive coupling in the schizophrenic brain.

Predictive assessment of models for dynamic functional connectivity

In neuroimaging, it has become evident that models of dynamic functional connectivity (dFC), which characterize how intrinsic brain organization changes over time, can provide a more detailed representation of brain function than traditional static analyses. Many dFC models in the literature represent functional brain networks as a meta-stable process with a discrete number of states; however, there is a lack of consensus on how to perform model selection and learn the number of states, as well as a lack of understanding of how different modeling assumptions influence the estimated state dynamics. To address these issues, we consider a predictive likelihood approach to model assessment, where models are evaluated based on their predictive performance on held-out test data. Examining several prominent models of dFC (in their probabilistic formulations) we demonstrate our framework on synthetic data, and apply it on two real-world examples: a face recognition EEG experiment and resting-state fMRI. Our results evidence that both EEG and fMRI are better characterized using dynamic modeling approaches than by their static counterparts, but we also demonstrate that one must be cautious when interpreting dFC because parameter settings and modeling assumptions, such as window lengths and emission models, can have a large impact on the estimated states and consequently on the interpretation of the brain dynamics.

Risk for affective disorders is associated with greater prefrontal gray matter volumes: A prospective longitudinal study

Background

Major depression and bipolar disorders aggregates in families and are linked with a wide range of neurobiological abnormalities including cortical gray matter (GM) alterations. Prospective studies of individuals at familial risk may expose the neural mechanisms underlying risk transmission.

Methods

We used voxel based morphometry to investigate changes in regional GM brain volume, over a seven-year period, in 37 initially healthy individuals having a mono- or di-zygotic twin diagnosed with major depression or bipolar disorder (high-risk group; mean age 41.6 yrs.) as compared to 36 individuals with no history of affective disorders in the index twin and first-degree relatives (low-risk group; mean age 38.5 yrs.).

Results

Groups did not differ in regional GM volume changes over time. However, independent of time, high-risk twins had significantly greater GM volumes in bilateral dorsal anterior cingulate, inferior frontal gyrus and temporoparietal regions as compared to low-risk twins. Further, individuals who developed an affective disorder at follow-up (n = 12), had relatively the largest GM volumes, both at baseline and follow-up, in the right dorsal anterior cingulate cortex and right inferior frontal cortex compared to high- and low-risk twins who remained well at follow-up.

Conclusion

This pattern of apparently stable grater regional GM volume may constitute a neural marker of an increased risk for developing an affective disorder in individuals at familial risk.

Simultaneous representation of a spectrum of dynamically changing value estimates during decision making

Decisions are based on value expectations derived from experience. We show that dorsal anterior cingulate cortex and three other brain regions hold multiple representations of choice value based on different timescales of experience organized in terms of systematic gradients across the cortex. Some parts of each area represent value estimates based on recent reward experience while others represent value estimates based on experience over the longer term. The value estimates within these areas interact with one another according to their temporal scaling. Some aspects of the representations change dynamically as the environment changes. The spectrum of value estimates may act as a flexible selection mechanism for combining experience-derived value information with other aspects of value to allow flexible and adaptive decisions in changing environments.

Infinite von Mises–Fisher Mixture Modeling of Whole Brain fMRI Data

Cluster analysis of functional magnetic resonance imaging (fMRI) data is often performed using gaussian mixture models, but when the time series are standardized such that the data reside on a hypersphere, this modeling assumption is questionable. The consequences of ignoring the underlying spherical manifold are rarely analyzed, in part due to the computational challenges imposed by directional statistics. In this letter, we discuss a Bayesian von Mises–Fisher (vMF) mixture model for data on the unit hypersphere and present an efficient inference procedure based on collapsed Markov chain Monte Carlo sampling. Comparing the vMF and gaussian mixture models on synthetic data, we demonstrate that the vMF model has a slight advantage inferring the true underlying clustering when compared to gaussian-based models on data generated from both a mixture of vMFs and a mixture of gaussians subsequently normalized. Thus, when performing model selection, the two models are not in agreement. Analyzing multisubject whole brain resting-state fMRI data from healthy adult subjects, we find that the vMF mixture model is considerably more reliable than the gaussian mixture model when comparing solutions across models trained on different groups of subjects, and again we find that the two models disagree on the optimal number of components. The analysis indicates that the fMRI data support more than a thousand clusters, and we confirm this is not a result of overfitting by demonstrating better prediction on data from held-out subjects. Our results highlight the utility of using directional statistics to model standardized fMRI data and demonstrate that whole brain segmentation of fMRI data requires a very large number of functional units in order to adequately account for the discernible statistical patterns in the data.

Are Movement Artifacts in Magnetic Resonance Imaging a Real Problem?-A Narrative Review

Movement artifacts compromise image quality and may interfere with interpretation, especially in magnetic resonance imaging (MRI) applications with low signal-to-noise ratio such as functional MRI or diffusion tensor imaging, and when imaging small lesions. High image resolution has high sensitivity to motion artifacts and often prolongs scan time that again aggravates movement artifacts. During the scan fast imaging techniques and sequences, optimal receiver coils, careful patient positioning, and instruction may minimize movement artifacts. Physiological noise sources are motion from respiration, flow and pulse coupled to cardiac cycles, from the swallowing reflex and small spontaneous head movements. Par example, in resting-state functional MRI spontaneous neuronal activity adds 1–2% of signal change, even under optimal conditions signal contributions from physiological noise remain a considerable fraction hereof. Movement tracking during imaging may allow for prospective correction or postprocessing steps separating signal and noise.

Task relevance differentially shapes ventral visual stream sensitivity to visible and invisible faces

Top-down modulations of the visual cortex can be driven by task relevance. Yet, several accounts propose that the perceptual inferences underlying conscious recognition involve similar top-down modulations of sensory responses. Studying the pure impact of task relevance on sensory responses requires dissociating it from the top-down influences underlying conscious recognition. Here, using visual masking to abolish perceptual consciousness in humans, we report that functional magnetic resonance imaging (fMRI) responses to invisible faces in the fusiform gyrus are enhanced when they are task-relevant, but suppressed when they are task-irrelevant compared to other object categories. Under conscious perceptual conditions, task-related modulations were also present but drastically reduced, with visible faces always eliciting greater activity in the fusiform gyrus compared to other object categories. Thus, task relevance crucially shapes the sensitivity of fusiform regions to face stimuli, leading from enhancement to suppression of neural activity when the top-down influences accruing from conscious recognition are prevented.

Quantifying functional connectivity in multi-subject fMRI data using component models

Functional magnetic resonance imaging (fMRI) is increasingly used to characterize functional connectivity between brain regions. Given the vast number of between-voxel interactions in high-dimensional fMRI data, it is an ongoing challenge to detect stable and generalizable functional connectivity in the brain among groups of subjects. Component models can be used to define subspace representations of functional connectivity that are more interpretable. It is, however, unclear which component model provides the optimal representation of functional networks for multi-subject fMRI datasets. A flexible cross-validation approach that assesses the ability of the models to predict voxel-wise covariance in new data, using three different measures of generalization was proposed. This framework is used to compare a range of component models with varying degrees of flexibility in their representation of functional connectivity, evaluated on both simulated and experimental resting-state fMRI data. It was demonstrated that highly flexible subject-specific component subspaces, as well as very constrained average models, are poor predictors of whole-brain functional connectivity, whereas the best-generalizing models account for subject variability within a common spatial subspace. Within this set of models, spatial Independent Component Analysis (sICA) on concatenated data provides more interpretable brain patterns, whereas a consistent-covariance model that accounts for subject-specific network scaling (PARAFAC2) provides greater stability in functional connectivity relationships between components and their spatial representations. The proposed evaluation framework is a promising quantitative approach to evaluating component models, and reveals important differences between subspace models in terms of predictability, robustness, characterization of subject variability, and interpretability of the model parameters. Hum Brain Mapp 38:882-899, 2017. © 2016 Wiley Periodicals, Inc.

Chasing probabilities - Signaling negative and positive prediction errors across domains

Adaptive actions build on internal probabilistic models of possible outcomes that are tuned according to the errors of their predictions when experiencing an actual outcome. Prediction errors (PEs) inform choice behavior across a diversity of outcome domains and dimensions, yet neuroimaging studies have so far only investigated such signals in singular experimental contexts. It is thus unclear whether the neuroanatomical distribution of PE encoding reported previously pertains to computational features that are invariant with respect to outcome valence, sensory domain, or some combination of the two. We acquired functional MRI data while volunteers performed four probabilistic reversal learning tasks which differed in terms of outcome valence (reward-seeking versus punishment-avoidance) and domain (abstract symbols versus facial expressions) of outcomes. We found that ventral striatum and frontopolar cortex coded increasingly positive PEs, whereas dorsal anterior cingulate cortex (dACC) traced increasingly negative PEs, irrespectively of the outcome dimension. Individual reversal behavior was unaffected by context manipulations and was predicted by activity in dACC and right inferior frontal gyrus (IFG). The stronger the response to negative PEs in these areas, the lower was the tendency to reverse choice behavior in response to negative events, suggesting that these regions enforce a rule-based strategy across outcome dimensions. Outcome valence influenced PE-related activity in left amygdala, IFG, and dorsomedial prefrontal cortex, where activity selectively scaled with increasingly positive PEs in the reward-seeking but not punishment-avoidance context, irrespective of sensory domain. Left amygdala displayed an additional influence of sensory domain. In the context of avoiding punishment, amygdala activity increased with increasingly negative PEs, but only for facial stimuli, indicating an integration of outcome valence and sensory domain during probabilistic choices.

Resting-state connectivity predicts levodopa-induced dyskinesias in Parkinson's disease

Background

Levodopa‐induced dyskinesias are a common side effect of dopaminergic therapy in PD, but their neural correlates remain poorly understood.

Objectives

This study examines whether dyskinesias are associated with abnormal dopaminergic modulation of resting‐state cortico‐striatal connectivity.

Methods

Twelve PD patients with peak‐of‐dose dyskinesias and 12 patients without dyskinesias were withdrawn from dopaminergic medication. All patients received a single dose of fast‐acting soluble levodopa and then underwent resting‐state functional magnetic resonance imaging before any dyskinesias emerged. Levodopa‐induced modulation of cortico‐striatal resting‐state connectivity was assessed between the putamen and the following 3 cortical regions of interest: supplementary motor area, primary sensorimotor cortex, and right inferior frontal gyrus. These functional connectivity measures were entered into a linear support vector classifier to predict whether an individual patient would develop dyskinesias after levodopa intake. Linear regression analysis was applied to test which connectivity measures would predict dyskinesia severity.

Results

Dopaminergic modulation of resting‐state connectivity between the putamen and primary sensorimotor cortex in the most affected hemisphere predicted whether patients would develop dyskinesias with a specificity of 100% and a sensitivity of 91% (P < .0001). Modulation of resting‐state connectivity between the supplementary motor area and putamen predicted interindividual differences in dyskinesia severity (R² = 0.627, P = .004). Resting‐state connectivity between the right inferior frontal gyrus and putamen neither predicted dyskinesia status nor dyskinesia severity.

Conclusions

The results corroborate the notion that altered dopaminergic modulation of cortico‐striatal connectivity plays a key role in the pathophysiology of dyskinesias in PD. © 2016 International Parkinson and Movement Disorder Society

Measuring motion-induced B0 -fluctuations in the brain using field probes

PURPOSE

Fluctuations of the background magnetic field (B0 ) due to body and breathing motion can lead to significant artifacts in brain imaging at ultrahigh field. Corrections based on real-time sensing using external field probes show great potential. This study evaluates different aspects of field interpolation from these probes into the brain which is implicit in such methods. Measurements and simulations were performed to quantify how well B0 -fluctuations in the brain due to body and breathing motion are reflected in external field probe measurements.

METHODS

Field probe measurements were compared with scanner acquired B0 -maps from experiments with breathing and shoulder movements. A realistic simulation of B0 -fluctuations caused by breathing was performed, and used for testing different sets of field probe positions.

RESULTS

The B0 -fluctuations were well reflected in the field probe measurements in the shoulder experiments, while the breathing experiments showed only moderate correspondence. The simulations showed the importance of the probe positions, and that performing full 3(rd) order corrections based on 16 field probes is not recommended.

CONCLUSION

Methods for quantitative assessment of the field interpolation problem were developed and demonstrated. Field corrections based on external field measurements show great potential, although potential pitfalls were identified.

Abnormal dopaminergic modulation of striato-cortical networks underlies levodopa-induced dyskinesias in humans

Dopaminergic signalling in the striatum contributes to reinforcement of actions and motivational enhancement of motor vigour. Parkinson's disease leads to progressive dopaminergic denervation of the striatum, impairing the function of cortico-basal ganglia networks. While levodopa therapy alleviates basal ganglia dysfunction in Parkinson's disease, it often elicits involuntary movements, referred to as levodopa-induced peak-of-dose dyskinesias. Here, we used a novel pharmacodynamic neuroimaging approach to identify the changes in cortico-basal ganglia connectivity that herald the emergence of levodopa-induced dyskinesias. Twenty-six patients with Parkinson's disease (age range: 51-84 years; 11 females) received a single dose of levodopa and then performed a task in which they had to produce or suppress a movement in response to visual cues. Task-related activity was continuously mapped with functional magnetic resonance imaging. Dynamic causal modelling was applied to assess levodopa-induced modulation of effective connectivity between the pre-supplementary motor area, primary motor cortex and putamen when patients suppressed a motor response. Bayesian model selection revealed that patients who later developed levodopa-induced dyskinesias, but not patients without dyskinesias, showed a linear increase in connectivity between the putamen and primary motor cortex after levodopa intake during movement suppression. Individual dyskinesia severity was predicted by levodopa-induced modulation of striato-cortical feedback connections from putamen to the pre-supplementary motor area (Pcorrected = 0.020) and primary motor cortex (Pcorrected = 0.044), but not feed-forward connections from the cortex to the putamen. Our results identify for the first time, aberrant dopaminergic modulation of striatal-cortical connectivity as a neural signature of levodopa-induced dyskinesias in humans. We argue that excessive striato-cortical connectivity in response to levodopa produces an aberrant reinforcement signal producing an abnormal motor drive that ultimately triggers involuntary movements.

Non-parametric Bayesian graph models reveal community structure in resting state fMRI

Modeling of resting state functional magnetic resonance imaging (rs-fMRI) data using network models is of increasing interest. It is often desirable to group nodes into clusters to interpret the communication patterns between nodes. In this study we consider three different nonparametric Bayesian models for node clustering in complex networks. In particular, we test their ability to predict unseen data and their ability to reproduce clustering across datasets. The three generative models considered are the Infinite Relational Model (IRM), Bayesian Community Detection (BCD), and the Infinite Diagonal Model (IDM). The models define probabilities of generating links within and between clusters and the difference between the models lies in the restrictions they impose upon the between-cluster link probabilities. IRM is the most flexible model with no restrictions on the probabilities of links between clusters. BCD restricts the between-cluster link probabilities to be strictly lower than within-cluster link probabilities to conform to the community structure typically seen in social networks. IDM only models a single between-cluster link probability, which can be interpreted as a background noise probability. These probabilistic models are compared against three other approaches for node clustering, namely Infomap, Louvain modularity, and hierarchical clustering. Using 3 different datasets comprising healthy volunteers' rs-fMRI we found that the BCD model was in general the most predictive and reproducible model. This suggests that rs-fMRI data exhibits community structure and furthermore points to the significance of modeling heterogeneous between-cluster link probabilities.

The acute brain response to levodopa heralds dyskinesias in Parkinson disease

Objective

In Parkinson disease (PD), long‐term treatment with the dopamine precursor levodopa gradually induces involuntary “dyskinesia” movements. The neural mechanisms underlying the emergence of levodopa‐induced dyskinesias in vivo are still poorly understood. Here, we applied functional magnetic resonance imaging (fMRI) to map the emergence of peak‐of‐dose dyskinesias in patients with PD.

Methods

Thirteen PD patients with dyskinesias and 13 PD patients without dyskinesias received 200mg fast‐acting oral levodopa following prolonged withdrawal from their normal dopaminergic medication. Immediately before and after levodopa intake, we performed fMRI, while patients produced a mouse click with the right or left hand or no action (No‐Go) contingent on 3 arbitrary cues. The scan was continued for 45 minutes after levodopa intake or until dyskinesias emerged.

Results

During No‐Go trials, PD patients who would later develop dyskinesias showed an abnormal gradual increase of activity in the presupplementary motor area (preSMA) and the bilateral putamen. This hyperactivity emerged during the first 20 minutes after levodopa intake. At the individual level, the excessive No‐Go activity in the predyskinesia period predicted whether an individual patient would subsequently develop dyskinesias (p < 0.001) as well as severity of their day‐to‐day symptomatic dyskinesias (p < 0.001).

Interpretation

PD patients with dyskinesias display an immediate hypersensitivity of preSMA and putamen to levodopa, which heralds the failure of neural networks to suppress involuntary dyskinetic movements. Ann Neurol 2014;75:829–836