Search

Abstract Segmentation of brain structures on magnetic resonance imaging (MRI) is a highly relevant neuroimaging topic, as it is a prerequisite for different analyses such as volumetry or shape analysis. Automated segmentation facilitates the study of brain structures in larger cohorts when compared with manual segmentation, which is time-consuming. However, the development of most automated methods relies on large and manually annotated datasets, which limits the generalizability of these methods. Recently, new techniques using synthetic images have emerged, reducing the need for manual annotation. Here we provide a dataset composed of label maps built from publicly available ultra-high resolution ex vivo MRI from 10 whole hemispheres, which can be used to develop segmentation methods using synthetic data. The label maps are obtained with a combination of manual labels for the hypothalamic regions and automated segmentations for the rest of the brain, and mirrored to simulate entire brains. We also provide the pre-processed ex vivo scans, as this dataset can support future projects to include other structures after these are manually segmented.

Non-human primates (NHPs) have become key for translational research in noninvasive brain stimulation (NIBS). However, in order to create comparable stimulation conditions for humans it is vital to study the accuracy of current modeling practices across species. Numerical models to simulate electric fields are an important tool for experimental planning in NHPs and translation to human studies. It is thus essential whether and to what extent the anatomical details of NHP models agree with current modeling practices when calculating NIBS electric fields. Here, we create highly accurate head models of two non-human primates (NHP) MR data. We evaluate how muscle tissue and head field of view (depending on MRI parameters) affect simulation results in transcranial electric and magnetic stimulation (TES and TMS). Our findings indicate that the inclusion of anisotropic muscle can affect TES electric field strength up to 22% while TMS is largely unaffected. Additionally, comparing a full head model to a cropped head model illustrates the impact of head field of view on electric fields for both TES and TMS. We find opposing effects between TES and TMS with an increase up to 24.8% for TES and a decrease up to 24.6% for TMS for the cropped head model compared to the full head model. Our results provide important insights into the level of anatomical detail needed for NHP head models and can inform future translational efforts for NIBS studies.

Generating realistic volume conductor models for forward calculations in electroencephalography (EEG) is not trivial and several factors contribute to the accuracy of such models, two of which are its anatomical accuracy and the accuracy with which electrode positions are known. Here, we investigate effects of anatomical accuracy by comparing forward solutions from SimNIBS, a tool which allows state-of-the-art anatomical modeling, with well-established pipelines in MNE-Python and FieldTrip. We also compare different ways of specifying electrode locations when digitized positions are not available such as transformation of measured positions from standard space and transformation of a manufacturer layout.Substantial effects of anatomical accuracy were seen throughout the entire brain both in terms of field topography and magnitude with SimNIBS generally being more accurate than the pipelines in MNE-Python and FieldTrip. Topographic and magnitude effects were particularly pronounced for MNE-Python which uses a three-layer boundary element method (BEM) model. We attribute these mainly to the coarse representation of the anatomy used in this model, in particular differences in skull and cerebrospinal fluid (CSF). Effects of electrode specification method were evident in occipital and posterior areas when using a transformed manufacturer layout whereas transforming measured positions from standard space generally resulted in smaller errors.We suggest modeling the anatomy of the volume conductor as accurately possible and we hope to facilitate this by making it easy to export simulations from SimNIBS to MNE-Python and FieldTrip for further analysis. Likewise, if digitized electrode positions are not available, a set of measured positions on a standard head template may be preferable to those specified by the manufacturer.

Transcranial focused Ultrasound Stimulation (TUS) at low intensities is emerging as a novel non-invasive brain stimulation method with higher spatial resolution than established transcranial stimulation methods and the ability to selectively stimulate also deep brain areas. Accurate control of the focus position and strength of the TUS acoustic waves is important to enable a beneficial use of the high spatial resolution and to ensure safety. As the human skull causes strong attenuation and distortion of the waves, simulations of the transmitted waves are needed to accurately determine the TUS dose distribution inside the cranial cavity. The simulations require information of the skull morphology and its acoustic properties. Ideally, they are informed by computed tomography (CT) images of the individual head. However, suited individual imaging data is often not readily available. For this reason, we here introduce and validate a head template that can be used to estimate the average effects of the skull on the TUS acoustic wave in the population.The template was created from CT images of the heads of 29 individuals of different ages (between 20–50 years), gender and ethnicity using an iterative non-linear co-registration procedure. For validation, we compared acoustic and thermal simulations based on the template to the average of the simulation results of all 29 individual datasets. Acoustic simulations were performed for a model of a focused transducer driven at 500 kHz, placed at 24 standardized positions by means of the EEG 10–10 system. Additional simulations at 250 kHz and 750 kHz at 16 of the positions were used for further confirmation. The amount of ultrasound-induced heating at 500 kHz was estimated for the same 16 transducer positions. Our results show that the template represents the median of the acoustic pressure and temperature maps from the individuals reasonably well in most cases. This underpins the usefulness of the template for the planning and optimization of TUS interventions in studies of healthy young adults. Our results further indicate that the amount of variability between the individual simulation results depends on the position. Specifically, the simulated ultrasound-induced heating inside the skull exhibited strong interindividual variability for three posterior positions close to the midline, caused by a high variability of the local skull shape and composition. This should be taken into account when interpreting simulation results based on the template.

Background and objective: Transcranial direct current stimulation (tDCS) has wide ranging applications in neuro-behavioural and physiological research, and in neurological rehabilitation. However, it is currently limited by substantial inter-subject variability in responses, which may be explained, at least in part, by anatomical differences that lead to variability in the electric field (E-field) induced in the cortex. Here, we tested whether the variability in the E-field in the stimulated cortex during anodal tDCS, estimated using computational simulations, explains the variability in tDCS induced changes in GABA, a neurophysiological marker of stimulation effect. Methods: Data from five previously conducted MRS studies were combined. The anode was placed over the left primary motor cortex (M1, 3 studies, N = 24) or right temporal cortex (2 studies, N = 32), with the cathode over the contralateral supraorbital ridge. Single voxel spectroscopy was performed in a 2x2x2cm voxel under the anode in all cases. MRS data were acquired before and either during or after 1 mA tDCS using either a sLASER sequence (7T) or a MEGA-PRESS sequence (3T). sLASER MRS data were analysed using LCModel, and MEGA-PRESS using FID-A and Gannet. E-fields were simulated in a finite element model of the head, based on individual structural MR images, using SimNIBS. Separate linear mixed effects models were run for each E-field variable (mean and 95th percentile; magnitude, and components normal and tangential to grey matter surface, within the MRS voxel). The model included effects of time (pre or post tDCS), E-field, grey matter volume in the MRS voxel, and a 3-way interaction between time, E-field and grey matter volume. Additionally, we ran a permutation analysis using PALM to determine whether E-field anywhere in the brain, not just in the MRS voxel, correlated with GABA change. Results: In M1, higher mean E-field magnitude was associated with greater anodal tDCS-induced decreases in GABA (t(24) = 3.24, p = 0.003). Further, the association between mean E-field magnitude and GABA change was moderated by the grey matter volume in the MRS voxel (t(24) = −3.55, p = 0.002). These relationships were consistent across all E-field variables except the mean of the normal component. No significant relationship was found between tDCS-induced GABA decrease and E-field in the temporal voxel. No significant clusters were found in the whole brain analysis. Conclusions: Our data suggest that the electric field induced by tDCS within the brain is variable, and is significantly related to anodal tDCS-induced decrease in GABA, a key neurophysiological marker of stimulation. These findings strongly support individualised dosing of tDCS, at least in M1. Further studies examining E-fields in relation to other outcome measures, including behaviour, will help determine the optimal E-fields required for any desired effects.

Abstract Methodological studies investigating transcranial direct current stimulation (tDCS) over the left dorsolateral prefrontal cortex (lDLPFC) in paediatric populations are limited. Therefore, we investigated in a paediatric population whether stimulation success of multichannel tDCS over the lDLPFC depends on concurrent task performance and individual head anatomy. In a randomised, sham-controlled, double-blind crossover study 22 healthy participants (10–17 years) received 2 mA multichannel anodal tDCS (atDCS) over the lDLPFC with and without a 2-back working memory (WM) task. After stimulation, the 2-back task and a Flanker task were performed. Resting state and task-related EEG were recorded. In 16 participants we calculated the individual electric field (E-field) distribution. Performance and neurophysiological activity in the 2-back task were not affected by atDCS. atDCS reduced reaction times in the Flanker task, independent of whether atDCS had been combined with the 2-back task. Flanker task related beta oscillation increased following stimulation without 2-back task performance. atDCS effects were not correlated with the E-field. We found no effect of multichannel atDCS over the lDLPFC on WM in children/adolescents but a transfer effect on interference control. While this effect on behaviour was independent of concurrent task performance, neurophysiological activity might be more sensitive to cognitive activation during stimulation. However, our results are limited by the small sample size, the lack of an active control group and variations in WM performance.

Background: Head and brain anatomy have been related to e-field strength induced by transcranial electrical stimulation (tES). Individualization based on anatomic factors require high-quality structural magnetic resonance images, which are not always available. Head circumference (HC) can serve as an alternative means, but its linkage to electric field strength has not yet been established. Methods: We simulated electric fields induced by tES based on individual T1w- and T2w-images of 47 healthy adults, for four conventional (“standard”) and four corresponding focal (”4x1”) electrode montages. Associations of electric field strength with individual HC were calculated using linear mixed models. Results: Larger HC was associated with lower electric field strength across montages. We provide mathematical equations to estimate individual electric field strength based on the HC. Conclusion: HC can be used as an alternative to estimate interindividual differences of the tES-induced electric field strength and to prospectively individualize stimulation dose, e.g., in the clinical context.

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.

Magnetic resonance current density imaging (MRCDI) of the human brain aims to reconstruct the current density distribution caused by transcranial electric stimulation from MR-based measurements of the current-induced magnetic fields. So far, the MRCDI data acquisition achieves only a low signal-to-noise ratio, does not provide a full volume coverage and lacks data from the scalp and skull regions. In addition, it is only sensitive to the component of the current-induced magnetic field parallel to the scanner field. The reconstruction problem thus involves coping with noisy and incomplete data, which makes it mathematically challenging. Most existing reconstruction methods have been validated using simulation studies and measurements in phantoms with simplified geometries. Only one reconstruction method, the projected current density algorithm, has been applied to human in-vivo data so far, however resulting in blurred current density estimates even when applied to noise-free simulated data.We analyze the underlying causes for the limited performance of the projected current density algorithm when applied to human brain data. In addition, we compare it with an approach that relies on the optimization of the conductivities of a small number of tissue compartments of anatomically detailed head models reconstructed from structural MR data. Both for simulated ground truth data and human in-vivo MRCDI data, our results indicate that the estimation of current densities benefits more from using a personalized volume conductor model than from applying the projected current density algorithm. In particular, we introduce a hierarchical statistical testing approach as a principled way to test and compare the quality of reconstructed current density images that accounts for the limited signal-to-noise ratio of the human in-vivo MRCDI data and the fact that the ground truth of the current density is unknown for measured data. Our results indicate that the statistical testing approach constitutes a valuable framework for the further development of accurate volume conductor models of the head. Our findings also highlight the importance of tailoring the reconstruction approaches to the quality and specific properties of the available data.

Here we present a method for the simultaneous segmentation of white matter lesions and normal-appearing neuroanatomical structures from multi-contrast brain MRI scans of multiple sclerosis patients. The method integrates a novel model for white matter lesions into a previously validated generative model for whole-brain segmentation. By using separate models for the shape of anatomical structures and their appearance in MRI, the algorithm can adapt to data acquired with different scanners and imaging protocols without retraining. We validate the method using four disparate datasets, showing robust performance in white matter lesion segmentation while simultaneously segmenting dozens of other brain structures. We further demonstrate that the contrast-adaptive method can also be safely applied to MRI scans of healthy controls, and replicate previously documented atrophy patterns in deep gray matter structures in MS. The algorithm is publicly available as part of the open-source neuroimaging package FreeSurfer.

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.

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, ppermutation = 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.

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.

Migraine with aura is a highly prevalent disorder involving transient neurological disturbances associated with migraine headache. While the pathophysiology is incompletely understood, findings from clinical and basic science studies indicate a potential key role of the thalamus in the mechanisms underlying migraine with and without aura. Two recent, clinic-based MRI studies investigated the volumes of individual thalamic nuclei in migraine patients with and without aura using two different data analysis methods. Both studies found differences of thalamic nuclei volumes between patients and healthy controls, but the results of the studies were not consistent.Here, we investigated whether migraine with aura is associated with changes in thalamic volume by analysing MRI data obtained from a large, cross-sectional population-based study which specifically included women with migraine with aura (N = 156), unrelated migraine-free matched controls (N = 126), and migraine aura-free co-twins (N = 29) identified from the Danish Twin Registry. We used two advanced, validated analysis methods to assess the volume of the thalamus and its nuclei; the MAGeT Brain Algorithm and a recently developed FreeSurfer-based method based on a probabilistic atlas of the thalamic nuclei combining ex vivo MRI and histology. These approaches were very similar to the methods used in each of the two previous studies.Between-group comparisons were corrected for potential effects of age, educational level, BMI, smoking, alcohol, and hypertension using a linear mixed model. Further, we used linear mixed models and visual inspection of data to assess relations between migraine aura frequency and thalamic nuclei volumes in patients. In addition, we performed paired t-tests to compare volumes of twin pairs (N = 29) discordant for migraine with aura. None of our analyses showed any between-group differences in volume of the thalamus or of individual thalamic nuclei. Our results indicate that the pathophysiology of migraine with aura does not involve alteration of thalamic volume.

During the past decade, it became clear that the electric field elicited by non-invasive brain stimulation (NIBS) techniques such as transcranial direct current stimulation (tDCS) and transcranial magnetic stimulation (TMS) are substantially influenced by variations in individual head and brain anatomy. In addition to structural variations in the healthy, several psychiatric disorders are characterized by anatomical alterations that are likely to further constrain the intracerebral effects of NIBS. Here, we present high-resolution realistic head models derived from structural magnetic resonance imaging data of 19 healthy adults and 19 patients diagnosed with major depressive disorder (MDD). By using a freely available software package for modelling the electric fields induced by different NIBS protocols, we show that our head models are well-suited for assessing inter-individual and between-group variability in the magnitude and focality of tDCS-induced electric fields for two protocols targeting the left dorsolateral prefrontal cortex.