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.

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.

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.

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.

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.

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 (CNRMTw ) 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 CNRMTw and severity of nonmotor symptoms. RESULTS: Mean CNRMTw of the right LC was reduced in patients relative to controls. Voxelwise and slicewise analyses showed that the attenuation of CNRMTw was confined to the right mid-caudal LC and linked regional CNRMTw to nonmotor symptoms. CNRMTw attenuation in the left mid-caudal LC was associated with the orthostatic drop in systolic blood pressure, whereas CNRMTw 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 CNRMTw 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 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson Movement Disorder Society.

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.

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.

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.

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.

Value and limitations of intracranial recordings for validating electric field modeling for transcranial brain stimulation.

Comparing electric field simulations from individualized head models against in-vivo intra-cranial recordings is considered the gold standard for direct validation of computational field modeling for transcranial brain stimulation and brain mapping techniques such as electro- and magnetoencephalography. The measurements also help to improve simulation accuracy by pinning down the factors having the largest influence on the simulations. Here we compare field simulations from four different automated pipelines against intracranial voltage recordings in an existing dataset of 14 epilepsy patients. We show that modeling differences in the pipelines lead to notable differences in the simulated electric field distributions that are often large enough to change the conclusions regarding the dose distribution and strength in the brain. Specifically, differences in the automatic segmentations of the head anatomy from structural magnetic resonance images are a major factor contributing to the observed field differences. However, the differences in the simulated fields are not reflected in the comparison between the simulations and intra-cranial measurements. This apparent mismatch is partly explained by the noisiness of the intra-cranial measurements, which renders comparisons between the methods inconclusive. We further demonstrate that a standard regression analysis, which ignores uncertainties in the simulations, leads to a strong bias in the estimated linear relationship between simulated and measured fields. Ignoring this bias leads to the incorrect conclusion that the models systematically misestimate the field strength in the brain. We propose a new Bayesian regression analysis of the data that yields unbiased parameter estimates, along with their uncertainties, and gives further insights to the fit between simulations and measurements. Specifically, the unbiased results give only weak support for systematic misestimations of the fields by the models.

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.

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.

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.

Functional neuroimaging of recovery from motor conversion disorder: A case report.

A patient with motor conversion disorder presented with a functional paresis of the left hand. After exclusion of structural brain damage, she was repeatedly examined with whole-brain functional magnetic resonance imaging, while she performed visually paced finger-tapping tasks. The dorsal premotor cortex showed a bilateral deactivation in the acute-subacute phase. Recovery from unilateral hand paresis was associated with a gradual increase in task-based activation of the dorsal premotor cortex bilaterally. The right medial prefrontal cortex displayed the opposite pattern, showing initial task-based activation that gradually diminished with recovery. The inverse dynamics of premotor and medial prefrontal activity over time were found during unimanual finger-tapping with the affected and non-affected hand as well as during bimanual finger-tapping. These observations suggest that reduced premotor and increased medial prefrontal activity reflect an effector-independent cortical dysfunction in conversion paresis which gradually disappears in parallel with clinical remission of paresis. The results link the medial prefrontal and dorsal premotor areas to the generation of intentional actions. We hypothesise that an excessive 'veto' signal generated in medial prefrontal cortex along with decreased premotor activity might constitute the functional substrate of conversion disorder. This notion warrants further examination in a larger group of affected patients.

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.

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.

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.

Automatic skull segmentation from MR images for realistic volume conductor models of the head: Assessment of the state-of-the-art.

Anatomically realistic volume conductor models of the human head are important for accurate forward modeling of the electric field during transcranial brain stimulation (TBS), electro- (EEG) and magnetoencephalography (MEG). In particular, the skull compartment exerts a strong influence on the field distribution due to its low conductivity, suggesting the need to represent its geometry accurately. However, automatic skull reconstruction from structural magnetic resonance (MR) images is difficult, as compact bone has a very low signal in magnetic resonance imaging (MRI). Here, we evaluate three methods for skull segmentation, namely FSL BET2, the unified segmentation routine of SPM12 with extended spatial tissue priors, and the skullfinder tool of BrainSuite. To our knowledge, this study is the first to rigorously assess the accuracy of these state-of-the-art tools by comparison with CT-based skull segmentations on a group of ten subjects. We demonstrate several key factors that improve the segmentation quality, including the use of multi-contrast MRI data, the optimization of the MR sequences and the adaptation of the parameters of the segmentation methods. We conclude that FSL and SPM12 achieve better skull segmentations than BrainSuite. The former methods obtain reasonable results for the upper part of the skull when a combination of T1- and T2-weighted images is used as input. The SPM12-based results can be improved slightly further by means of simple morphological operations to fix local defects. In contrast to FSL BET2, the SPM12-based segmentation with extended spatial tissue priors and the BrainSuite-based segmentation provide coarse reconstructions of the vertebrae, enabling the construction of volume conductor models that include the neck. We exemplarily demonstrate that the extended models enable a more accurate estimation of the electric field distribution during transcranial direct current stimulation (tDCS) for montages that involve extraencephalic electrodes. The methods provided by FSL and SPM12 are integrated into pipelines for the automatic generation of realistic head models based on tetrahedral meshes, which are distributed as part of the open-source software package SimNIBS for field calculations for transcranial brain stimulation.

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.