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. Building on 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 $\mu $m, we propose a Multi-resolution U-Nets framework 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Fast and sequence-adaptive whole-brain segmentation using parametric Bayesian modeling.

Quantitative analysis of magnetic resonance imaging (MRI) scans of the brain requires accurate automated segmentation of anatomical structures. A desirable feature for such segmentation methods is to be robust against changes in acquisition platform and imaging protocol. In this paper we validate the performance of a segmentation algorithm designed to meet these requirements, building upon generative parametric models previously used in tissue classification. The method is tested on four different datasets acquired with different scanners, field strengths and pulse sequences, demonstrating comparable accuracy to state-of-the-art methods on T1-weighted scans while being one to two orders of magnitude faster. The proposed algorithm is also shown to be robust against small training datasets, and readily handles images with different MRI contrast as well as multi-contrast data.

High-resolution segmentations of the hypothalamus and its subregions for training of segmentation models.

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.

A next-generation, histological atlas of the human brain and its application to automated brain MRI segmentation.

Magnetic resonance imaging (MRI) is the standard tool to image the human brain in vivo. In this domain, digital brain atlases are essential for subject-specific segmentation of anatomical regions of interest (ROIs) and spatial comparison of neuroanatomy from different subjects in a common coordinate frame. High-resolution, digital atlases derived from histology (e.g., Allen atlas [7], BigBrain [13], Julich [15]), are currently the state of the art and provide exquisite 3D cytoarchitectural maps, but lack probabilistic labels throughout the whole brain. Here we present NextBrain, a next-generation probabilistic atlas of human brain anatomy built from serial 3D histology and corresponding highly granular delineations of five whole brain hemispheres. We developed AI techniques to align and reconstruct ~10,000 histological sections into coherent 3D volumes with joint geometric constraints (no overlap or gaps between sections), as well as to semi-automatically trace the boundaries of 333 distinct anatomical ROIs on all these sections. Comprehensive delineation on multiple cases enabled us to build the first probabilistic histological atlas of the whole human brain. Further, we created a companion Bayesian tool for automated segmentation of the 333 ROIs in any in vivo or ex vivo brain MRI scan using the NextBrain atlas. We showcase two applications of the atlas: automated segmentation of ultra-high-resolution ex vivo MRI and volumetric analysis of Alzheimer's disease and healthy brain ageing based on ~4,000 publicly available in vivo MRI scans. We publicly release: the raw and aligned data (including an online visualisation tool); the probabilistic atlas; the segmentation tool; and ground truth delineations for a 100 µm isotropic ex vivo hemisphere (that we use for quantitative evaluation of our segmentation method in this paper). By enabling researchers worldwide to analyse brain MRI scans at a superior level of granularity without manual effort or highly specific neuroanatomical knowledge, NextBrain holds promise to increase the specificity of MRI findings and ultimately accelerate our quest to understand the human brain in health and disease.

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.

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.

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.

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.