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.
•We introduce a new automated method for whole-head segmentation.•The method segments 15 different head tissues covering also the neck.•The segmentation accuracy and robustness compare favorably to existing tools.•Choice of scan parameters can cause segmentation and simulation errors up to 30%.•Including extra tissues into the simulation affects the electric field locally.
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.
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.
•Fields from different head modeling pipelines can differ up to 49%.•However, all modeling pipelines perform similarly in predicting the measurements.•Bayesian analysis robustly estimates the fit between simulations and recordings.•Previous standard regression analyses of the fit suffered from a systematic bias.•New results do not support a systematic “misestimation” of the fields by the models.
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.
Abstract For decades, fMRI data have been used to search for biomarkers for patients with schizophrenia. Still, firm conclusions are yet to be made, which is often attributed to the high internal ...heterogeneity of the disorder. A promising way to disentangle the heterogeneity is to search for subgroups of patients with more homogeneous biological profiles. We applied an unsupervised multiple co-clustering (MCC) method to identify subtypes using functional connectivity data from a multisite resting-state data set. We merged data from two publicly available databases and split the data into a discovery data set (143 patients and 143 healthy controls (HC)) and an external test data set (63 patients and 63 HC) from independent sites. On the discovery data, we investigated the stability of the clustering toward data splits and initializations. Subsequently we searched for cluster solutions, also called “views,” with a significant diagnosis association and evaluated these based on their subject and feature cluster separability, and correlation to clinical manifestations as measured with the positive and negative syndrome scale (PANSS). Finally, we validated our findings by testing the diagnosis association on the external test data. A major finding of our study was that the stability of the clustering was highly dependent on variations in the data set, and even across initializations, we found only a moderate subject clustering stability. Nevertheless, we still discovered one view with a significant diagnosis association. This view reproducibly showed an overrepresentation of schizophrenia patients in three subject clusters, and one feature cluster showed a continuous trend, ranging from positive to negative connectivity values, when sorted according to the proportions of patients with schizophrenia. When investigating all patients, none of the feature clusters in the view were associated with severity of positive, negative, and generalized symptoms, indicating that the cluster solutions reflect other disease related mechanisms.
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 T2$$ {}_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.
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.
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.
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.