Neuroanatomical Predictors of Transcranial Direct Current Stimulation (tDCS)-Induced Modifications in Neurocognitive Task Performance in Typically Developing Individuals.

Transcranial direct current stimulation (tDCS) is a noninvasive neuromodulation technique gaining more attention in neurodevelopmental disorders (NDDs). Due to the phenotypic heterogeneity of NDDs, tDCS is unlikely to be equally effective in all individuals. The present study aimed to establish neuroanatomical markers in typically developing (TD) individuals that may be used for the prediction of individual responses to tDCS. Fifty-seven male and female children received 2 mA anodal and sham tDCS, targeting the left dorsolateral prefrontal cortex (DLPFCleft), right inferior frontal gyrus, and bilateral temporoparietal junction. Response to tDCS was assessed based on task performance differences between anodal and sham tDCS in different neurocognitive tasks (N-back, flanker, Mooney faces detection, attentional emotional recognition task). Measures of cortical thickness (CT) and surface area (SA) were derived from 3 Tesla structural MRI scans. Associations between neuroanatomy and task performance were assessed using general linear models (GLM). Machine learning (ML) algorithms were employed to predict responses to tDCS. Vertex-wise estimates of SA were more closely linked to differences in task performance than measures of CT. Across ML algorithms, highest accuracies were observed for the prediction of N-back task performance differences following stimulation of the DLPFCleft, where 65% of behavioral variance was explained by variability in SA. Lower accuracies were observed for all other tasks and stimulated regions. This suggests that it may be possible to predict individual responses to tDCS for some behavioral measures and target regions. In the future, these models might be extended to predict treatment outcome in individuals with NDDs.

Segmentation of supragranular and infragranular layers in ultra-high resolution 7T ex vivo MRI of the human cerebral cortex.

Accurate labeling of specific layers in the human cerebral cortex is crucial for advancing our understanding of neurodevelopmental and neurodegenerative disorders. Leveraging recent advancements in ultra-high resolution ex vivo MRI, we present a novel semi-supervised segmentation model capable of identifying supragranular and infragranular layers in ex vivo MRI with unprecedented precision. On a dataset consisting of 17 whole-hemisphere ex vivo scans at 120 µm, we propose a multi-resolution U-Nets framework (MUS) that integrates global and local structural information, achieving reliable segmentation maps of the entire hemisphere, with Dice scores over 0.8 for supra- and infragranular layers. This enables surface modeling, atlas construction, anomaly detection in disease states, and cross-modality validation, while also paving the way for finer layer segmentation. Our approach offers a powerful tool for comprehensive neuroanatomical investigations and holds promise for advancing our mechanistic understanding of progression of neurodegenerative diseases.

Anatomical details affect electric field predictions for non-invasive brain stimulation in non-human primates.

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.

Evaluating the influence of anatomical accuracy and electrode positions on EEG forward solutions.

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.

Investigations of the subarachnoid space as a potential link between aura and headache in migraine: A case-control MRI study.

BACKGROUND: The connection between migraine aura and headache is poorly understood. Some patients experience migraine aura without headache, and patients with migraine aura with headache commonly experience milder headaches with age. The distance between the cerebral cortex and the overlying dura mater has been hypothesized to influence development of headache following aura. We tested this hypothesis by comparing approximated distances between visual cortical areas and overlying dura mater between female patients with migraine aura without headache and female patients with migraine aura with headache. METHODS: Twelve cases with migraine aura without headache and 45 age-matched controls with migraine aura with headache underwent 3.0 T MRI. We calculated average distances between the occipital lobes, between the calcarine sulci, and between the skull and visual areas V1, V2 and V3a. We also measured volumes of corticospinal fluid between the occipital lobes, between the calcarine sulci, and overlying visual areas V2 and V3a. We investigated the relationship between headache status, distances and corticospinal fluid volumes using conditional logistic regression. RESULTS: Distances between the occipital lobes, calcarine sulci and between the skull and V1, V2 and V3a did not differ between patients with migraine aura with headache and patients with migraine aura without headache. We found no differences in corticospinal fluid volumes between groups. CONCLUSION: We found no indication for a connection between visual migraine aura and headache based on cortico-cortical, cortex-to-skull distances, or corticospinal fluid volumes overlying visual cortical areas. Longitudinal studies with imaging sequences optimized for measuring the cortico-dural distance and a larger sample of patients are needed to further investigate the hypothesis.

A head template for computational dose modelling for transcranial focused ultrasound stimulation.

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.

SynthSeg: Segmentation of brain MRI scans of any contrast and resolution without retraining.

Despite advances in data augmentation and transfer learning, convolutional neural networks (CNNs) difficultly generalise to unseen domains. When segmenting brain scans, CNNs are highly sensitive to changes in resolution and contrast: even within the same MRI modality, performance can decrease across datasets. Here we introduce SynthSeg, the first segmentation CNN robust against changes in contrast and resolution. SynthSeg is trained with synthetic data sampled from a generative model conditioned on segmentations. Crucially, we adopt a domain randomisation strategy where we fully randomise the contrast and resolution of the synthetic training data. Consequently, SynthSeg can segment real scans from a wide range of target domains without retraining or fine-tuning, which enables straightforward analysis of huge amounts of heterogeneous clinical data. Because SynthSeg only requires segmentations to be trained (no images), it can learn from labels obtained by automated methods on diverse populations (e.g., ageing and diseased), thus achieving robustness to a wide range of morphological variability. We demonstrate SynthSeg on 5,000 scans of six modalities (including CT) and ten resolutions, where it exhibits unparallelled generalisation compared with supervised CNNs, state-of-the-art domain adaptation, and Bayesian segmentation. Finally, we demonstrate the generalisability of SynthSeg by applying it to cardiac MRI and CT scans.

tDCS induced GABA change is associated with the simulated electric field in M1, an effect mediated by grey matter volume in the MRS voxel.

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.

Linking lesions in sensorimotor cortex to contralateral hand function in multiple sclerosis: a 7†T MRI study.

Cortical lesions constitute a key manifestation of multiple sclerosis and contribute to clinical disability and cognitive impairment. Yet it is unknown whether local cortical lesions and cortical lesion subtypes contribute to domain-specific impairments attributable to the function of the lesioned cortex. In this cross-sectional study, we assessed how cortical lesions in the primary sensorimotor hand area relate to corticomotor physiology and sensorimotor function of the contralateral hand. Fifty relapse-free patients with relapsing-remitting or secondary-progressive multiple sclerosis and 28 healthy age- and sex-matched participants underwent whole-brain 7 T MRI to map cortical lesions. Brain scans were also used to estimate normalized brain volume, pericentral cortical thickness, white matter lesion fraction of the corticospinal tract, infratentorial lesion volume and the cross-sectional area of the upper cervical spinal cord. We tested sensorimotor hand function and calculated a motor and sensory composite score for each hand. In 37 patients and 20 healthy controls, we measured maximal motor-evoked potential amplitude, resting motor threshold and corticomotor conduction time with transcranial magnetic stimulation and the N20 latency from somatosensory-evoked potentials. Patients showed at least one cortical lesion in the primary sensorimotor hand area in 47 of 100 hemispheres. The presence of a lesion was associated with worse contralateral sensory (P = 0.014) and motor (P = 0.009) composite scores. Transcranial magnetic stimulation of a lesion-positive primary sensorimotor hand area revealed a decreased maximal motor-evoked potential amplitude (P < 0.001) and delayed corticomotor conduction (P = 0.002) relative to a lesion-negative primary sensorimotor hand area. Stepwise mixed linear regressions showed that the presence of a primary sensorimotor hand area lesion, higher white-matter lesion fraction of the corticospinal tract, reduced spinal cord cross-sectional area and higher infratentorial lesion volume were associated with reduced contralateral motor hand function. Cortical lesions in the primary sensorimotor hand area, spinal cord cross-sectional area and normalized brain volume were also associated with smaller maximal motor-evoked potential amplitude and longer corticomotor conduction times. The effect of cortical lesions on sensory function was no longer significant when controlling for MRI-based covariates. Lastly, we found that intracortical and subpial lesions had the largest effect on reduced motor hand function, intracortical lesions on reduced motor-evoked potential amplitude and leucocortical lesions on delayed corticomotor conduction. Together, this comprehensive multilevel assessment of sensorimotor brain damage shows that the presence of a cortical lesion in the primary sensorimotor hand area is associated with impaired corticomotor function of the hand, after accounting for damage at the subcortical level. The results also provide preliminary evidence that cortical lesion types may affect the various facets of corticomotor function differentially.

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.

Multichannel anodal tDCS over the left dorsolateral prefrontal cortex in a paediatric population.

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.

On the reconstruction of magnetic resonance current density images of the human brain: Pitfalls and perspectives.

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.

Estimation of individually induced e-field strength during transcranial electric stimulation using the head circumference.

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.

A contrast-adaptive method for simultaneous whole-brain and lesion segmentation in multiple sclerosis.

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.

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.

Migraine with aura in women is not associated with structural thalamic abnormalities.

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.

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.

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.

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.

A modality-adaptive method for segmenting brain tumors and organs-at-risk in radiation therapy planning.

In this paper we present a method for simultaneously segmenting brain tumors and an extensive set of organs-at-risk for radiation therapy planning of glioblastomas. The method combines a contrast-adaptive generative model for whole-brain segmentation with a new spatial regularization model of tumor shape using convolutional restricted Boltzmann machines. We demonstrate experimentally that the method is able to adapt to image acquisitions that differ substantially from any available training data, ensuring its applicability across treatment sites; that its tumor segmentation accuracy is comparable to that of the current state of the art; and that it captures most organs-at-risk sufficiently well for radiation therapy planning purposes. The proposed method may be a valuable step towards automating the delineation of brain tumors and organs-at-risk in glioblastoma patients undergoing radiation therapy.