Comparative validation of automated presurgical tractography based on constrained spherical deconvolution and diffusion tensor imaging with direct electrical stimulation.

OBJECTIVES: Accurate presurgical brain mapping enables preoperative risk assessment and intraoperative guidance. This cross-sectional study investigated whether constrained spherical deconvolution (CSD) methods were more accurate than diffusion tensor imaging (DTI)-based methods for presurgical white matter mapping using intraoperative direct electrical stimulation (DES) as the ground truth. METHODS: Five different tractography methods were compared (three DTI-based and two CSD-based) in 22 preoperative neurosurgical patients undergoing surgery with DES mapping. The corticospinal tract (CST, N = 20) and arcuate fasciculus (AF, N = 7) bundles were reconstructed, then minimum distances between tractograms and DES coordinates were compared between tractography methods. Receiver-operating characteristic (ROC) curves were used for both bundles. For the CST, binary agreement, linear modeling, and posthoc testing were used to compare tractography methods while correcting for relative lesion and bundle volumes. RESULTS: Distance measures between 154 positive (functional response, pDES) and negative (no response, nDES) coordinates, and 134 tractograms resulted in 860 data points. Higher agreement was found between pDES coordinates and CSD-based compared to DTI-based tractograms. ROC curves showed overall higher sensitivity at shorter distance cutoffs for CSD (8.5 mm) compared to DTI (14.5 mm). CSD-based CST tractograms showed significantly higher agreement with pDES, which was confirmed by linear modeling and posthoc tests (PFWE < .05). CONCLUSIONS: CSD-based CST tractograms were more accurate than DTI-based ones when validated using DES-based assessment of motor and sensory function. This demonstrates the potential benefits of structural mapping using CSD in clinical practice.

Isolating the sources of pipeline-variability in group-level task-fMRI results.

Task-fMRI researchers have great flexibility as to how they analyze their data, with multiple methodological options to choose from at each stage of the analysis workflow. While the development of tools and techniques has broadened our horizons for comprehending the complexities of the human brain, a growing body of research has highlighted the pitfalls of such methodological plurality. In a recent study, we found that the choice of software package used to run the analysis pipeline can have a considerable impact on the final group-level results of a task-fMRI investigation (Bowring et al., 2019, BMN). Here we revisit our work, seeking to identify the stages of the pipeline where the greatest variation between analysis software is induced. We carry out further analyses on the three datasets evaluated in BMN, employing a common processing strategy across parts of the analysis workflow and then utilizing procedures from three software packages (AFNI, FSL, and SPM) across the remaining steps of the pipeline. We use quantitative methods to compare the statistical maps and isolate the main stages of the workflow where the three packages diverge. Across all datasets, we find that variation between the packages' results is largely attributable to a handful of individual analysis stages, and that these sources of variability were heterogeneous across the datasets (e.g., choice of first-level signal model had the most impact for the balloon analog risk task dataset, while first-level noise model and group-level model were more influential for the false belief and antisaccade task datasets, respectively). We also observe areas of the analysis workflow where changing the software package causes minimal differences in the final results, finding that the group-level results were largely unaffected by which software package was used to model the low-frequency fMRI drifts.

Evaluation of movement and brain activity.

Clinical neurophysiology studies can contribute important information about the physiology of human movement and the pathophysiology and diagnosis of different movement disorders. Some techniques can be accomplished in a routine clinical neurophysiology laboratory and others require some special equipment. This review, initiating a series of articles on this topic, focuses on the methods and techniques. The methods reviewed include EMG, EEG, MEG, evoked potentials, coherence, accelerometry, posturography (balance), gait, and sleep studies. Functional MRI (fMRI) is also reviewed as a physiological method that can be used independently or together with other methods. A few applications to patients with movement disorders are discussed as examples, but the detailed applications will be the subject of other articles.

A simple permutation-based test of intermodal correspondence.

Many key findings in neuroimaging studies involve similarities between brain maps, but statistical methods used to measure these findings have varied. Current state-of-the-art methods involve comparing observed group-level brain maps (after averaging intensities at each image location across multiple subjects) against spatial null models of these group-level maps. However, these methods typically make strong and potentially unrealistic statistical assumptions, such as covariance stationarity. To address these issues, in this article we propose using subject-level data and a classical permutation testing framework to test and assess similarities between brain maps. Our method is comparable to traditional permutation tests in that it involves randomly permuting subjects to generate a null distribution of intermodal correspondence statistics, which we compare to an observed statistic to estimate a p-value. We apply and compare our method in simulated and real neuroimaging data from the Philadelphia Neurodevelopmental Cohort. We show that our method performs well for detecting relationships between modalities known to be strongly related (cortical thickness and sulcal depth), and it is conservative when an association would not be expected (cortical thickness and activation on the n-back working memory task). Notably, our method is the most flexible and reliable for localizing intermodal relationships within subregions of the brain and allows for generalizable statistical inference.

Optimal design of on-scalp electromagnetic sensor arrays for brain source localisation.

Optically pumped magnetometers (OPMs) are quickly widening the scopes of noninvasive neurophysiological imaging. The possibility of placing these magnetic field sensors on the scalp allows not only to acquire signals from people in movement, but also to reduce the distance between the sensors and the brain, with a consequent gain in the signal-to-noise ratio. These advantages make the technique particularly attractive to characterise sources of brain activity in demanding populations, such as children and patients with epilepsy. However, the technology is currently in an early stage, presenting new design challenges around the optimal sensor arrangement and their complementarity with other techniques as electroencephalography (EEG). In this article, we present an optimal array design strategy focussed on minimising the brain source localisation error. The methodology is based on the Cramér-Rao bound, which provides lower error bounds on the estimation of source parameters regardless of the algorithm used. We utilise this framework to compare whole head OPM arrays with commercially available electro/magnetoencephalography (E/MEG) systems for localising brain signal generators. In addition, we study the complementarity between EEG and OPM-based MEG, and design optimal whole head systems based on OPMs only and a combination of OPMs and EEG electrodes for characterising deep and superficial sources alike. Finally, we show the usefulness of the approach to find the nearly optimal sensor positions minimising the estimation error bound in a given cortical region when a limited number of OPMs are available. This is of special interest for maximising the performance of small scale systems to ad hoc neurophysiological experiments, a common situation arising in most OPM labs.

Degree of improving TMS focality through a geometrically stable solution of an inverse TMS problem.

The Transcranial Magnetic Stimulation (TMS) inverse problem (TMS-IP) investigated in this study aims to focus the TMS induced electric field close to a specified target point defined on the gray matter interface in the M1HAND area while otherwise minimizing it. The goal of the study is to numerically evaluate the degree of improvement of the TMS-IP solutions relative to the well-known sulcus-aligned mapping (a projection approach with the 90∘ local sulcal angle). In total, 1536 individual TMS-IP solutions have been analyzed for multiple target points and multiple subjects using the boundary element fast multipole method (BEM-FMM) as the forward solver. Our results show that the optimal TMS inverse-problem solutions improve the focality - reduce the size of the field "hot spot" and its deviation from the target - by approximately 21-33% on average for all considered subjects, all observation points, two distinct coil types, two segmentation types, two intracortical observation surfaces under study, and three tested values of the field threshold. The inverse-problem solutions with the maximized focality simultaneously improve the TMS mapping resolution (differentiation between neighbor targets separated by approximately 10 mm) although this improvement is quite modest. Coil position/orientation and conductivity uncertainties have been included into consideration as the corresponding de-focalization factors. The present results will change when the levels of uncertainties change. Our results also indicate that the accuracy of the head segmentation critically influences the expected TMS-IP performance.

An optimized registration workflow and standard geometric space for small animal brain imaging.

The reliability of scientific results critically depends on reproducible and transparent data processing. Cross-subject and cross-study comparability of imaging data in general, and magnetic resonance imaging (MRI) data in particular, is contingent on the quality of registration to a standard reference space. In small animal MRI this is not adequately provided by currently used processing workflows, which utilize high-level scripts optimized for human data, and adapt animal data to fit the scripts, rather than vice-versa. In this fully reproducible article we showcase a generic workflow optimized for the mouse brain, alongside a standard reference space suited to harmonize data between analysis and operation. We introduce four separate metrics for automated quality control (QC), and a visualization method to aid operator inspection. Benchmarking this workflow against common legacy practices reveals that it performs more consistently, better preserves variance across subjects while minimizing variance across sessions, and improves both volume and smoothness conservation RMSE approximately 2-fold. We propose this open source workflow and the QC metrics as a new standard for small animal MRI registration, ensuring workflow robustness, data comparability, and region assignment validity, all of which are indispensable prerequisites for the comparability of scientific results across experiments and centers.

Comparison of Whole-Head Functional Near-Infrared Spectroscopy With Functional Magnetic Resonance Imaging and Potential Application in Pediatric Neurology.

BACKGROUND: Changes in cerebral blood flow in response to neuronal activation can be measured by time-dependent fluctuations in hemoglobin species within the brain; this is the basis of functional magnetic resonance imaging (fMRI) and functional near-infrared spectroscopy (fNIRS). There is a clinical need for portable neural imaging systems, such as fNIRS, to accommodate patients who are unable to tolerate an MR environment. OBJECTIVE: Our objective was to compare task-related full-head fNIRS and fMRI signals across cortical regions. METHODS: Eighteen healthy adults completed a same-day fNIRS-fMRI study, in which they performed right- and left-hand finger tapping tasks and a semantic-decision tones-decision task. First- and second-level general linear models were applied to both datasets. RESULTS: The finger tapping task showed that significant fNIRS channel activity over the contralateral primary motor cortex corresponded to surface fMRI activity. Similarly, significant fNIRS channel activity over the bilateral temporal lobe corresponded to the same primary auditory regions as surface fMRI during the semantic-decision tones-decision task. Additional channels were significant for this task that did not correspond to surface fMRI activity. CONCLUSION: Although both imaging modalities showed left-lateralized activation for language processing, the current fNIRS analysis did not show concordant or expected localization at the level necessary for clinical use in individual pediatric epileptic patients. Future work is needed to show whether fNIRS and fMRI are comparable at the source level so that fNIRS can be used in a clinical setting on individual patients. If comparable, such an imaging approach could be applied to children with neurological disorders.

Visualization of iron-rich subcortical structures in non-human primates in vivo by quantitative susceptibility mapping at 3T MRI.

Magnetic resonance imaging (MRI) is now an essential tool in the field of neuroscience involving non-human primates (NHP). Structural MRI scanning using T1-weighted (T1w) or T2-weighted (T2w) images provides anatomical information, particularly for experiments involving deep structures such as the basal ganglia and cerebellum. However, for certain subcortical structures, T1w and T2w image contrasts are insufficient for their detection of important anatomical details. To better visualize such structures in the macaque brain, we applied a relatively new method called quantitative susceptibility mapping (QSM), which enhances tissue contrast based on the local tissue magnetic susceptibility. The QSM significantly improved the visualization of important structures, including the ventral pallidum (VP), globus pallidus external and internal segments (GPe and GPi), substantia nigra (SN), subthalamic nucleus (STN) in the basal ganglia and the dentate nucleus (DN) in the cerebellum. We quantified this the contrast enhancement by systematically comparing of contrast-to-noise ratios (CNRs) of QSM images relative to the corresponding T1w and T2w images. In addition, QSM values of some structures were correlated to the age of the macaque subjects. These results identify the QSM method as a straightforward and useful tool for clearly visualizing details of subcortical structures that are invisible with more traditional scanning sequences.

Common functional localizers to enhance NHP & cross-species neuroscience imaging research.

Functional localizers are invaluable as they can help define regions of interest, provide cross-study comparisons, and most importantly, allow for the aggregation and meta-analyses of data across studies and laboratories. To achieve these goals within the non-human primate (NHP) imaging community, there is a pressing need for the use of standardized and validated localizers that can be readily implemented across different groups. The goal of this paper is to provide an overview of the value of localizer protocols to imaging research and we describe a number of commonly used or novel localizers within NHPs, and keys to implement them across studies. As has been shown with the aggregation of resting-state imaging data in the original PRIME-DE submissions, we believe that the field is ready to apply the same initiative for task-based functional localizers in NHP imaging. By coming together to collect large datasets across research group, implementing the same functional localizers, and sharing the localizers and data via PRIME-DE, it is now possible to fully test their robustness, selectivity and specificity. To do this, we reviewed a number of common localizers and we created a repository of well-established localizer that are easily accessible and implemented through the PRIME-RE platform.

Perturbing fMRI brain dynamics using transcranial direct current stimulation.

The dynamic nature of resting-state functional magnetic resonance imaging (fMRI) brain activity and connectivity has drawn great interest in the past decade. Specific temporal properties of fMRI brain dynamics, including metrics such as occurrence rate and transitions, have been associated with cognition and behaviors, indicating the existence of mechanism distruption in neuropsychiatric disorders. The development of new methods to manipulate fMRI brain dynamics will advance our understanding of these pathophysiological mechanisms from native observation to experimental mechanistic manipulation. In the present study, we applied repeated transcranial direct current stimulation (tDCS) to the right dorsolateral prefrontal cortex (rDLPFC) and the left orbitofrontal cortex (lOFC), during multiple simultaneous tDCS-fMRI sessions from 81 healthy participants to assess the modulatory effects of stimulating target brain regions on fMRI brain dynamics. Using the rDLPFC and the lOFC as seeds, respectively, we first identified two reoccurring co-activation patterns (CAPs) and calculated their temporal properties (e.g., occurrence rate and transitions) before administering tDCS. The spatial maps of CAPs were associated with different cognitive and disease domains using meta-analytical decoding analysis. We then investigated how active tDCS compared to sham tDCS in the modulation of the occurrence rates of these different CAPs and perturbations of transitions between CAPs. We found that by enhancing neuronal excitability of the rDLPFC and the lOFC, the occurrence rate of one CAP was significantly decreased while that of another CAP was significantly increased during the first 6 min of stimulation. Furthermore, these tDCS-associated changes persisted over subsequent testing sessions (both during and before/after tDCS) across three consecutive days. Active tDCS could perturb transitions between CAPs and a non-CAP state (when the rDLPFC and the lOFC were not activated), but not the transitions within CAPs. These results demonstrate the feasibility of modulating fMRI brain dynamics, and open new possibilities for discovering stimulation targets and dynamic connectivity patterns that can ensure the propagation of tDCS-induced neuronal excitability, which may facilitate the development of new treatments for disorders with altered dynamics.

Flexible annotation atlas of the mouse brain: combining and dividing brain structures of the Allen Brain Atlas while maintaining anatomical hierarchy.

A brain atlas is necessary for analyzing structure and function in neuroimaging research. Although various annotation volumes (AVs) for the mouse brain have been proposed, it is common in magnetic resonance imaging (MRI) of the mouse brain that regions-of-interest (ROIs) for brain structures (nodes) are created arbitrarily according to each researcher's necessity, leading to inconsistent ROIs among studies. One reason for such a situation is the fact that earlier AVs were fixed, i.e. combination and division of nodes were not implemented. This report presents a pipeline for constructing a flexible annotation atlas (FAA) of the mouse brain by leveraging public resources of the Allen Institute for Brain Science on brain structure, gene expression, and axonal projection. A mere two-step procedure with user-specified, text-based information and Python codes constructs FAA with nodes which can be combined or divided objectively while maintaining anatomical hierarchy of brain structures. Four FAAs with total node count of 4, 101, 866, and 1381 were demonstrated. Unique characteristics of FAA realized analysis of resting-state functional connectivity (FC) across the anatomical hierarchy and among cortical layers, which were thin but large brain structures. FAA can improve the consistency of whole brain ROI definition among laboratories by fulfilling various requests from researchers with its flexibility and reproducibility.

Mapping of multiple muscles with transcranial magnetic stimulation: absolute and relative test-retest reliability.

The spatial accuracy of transcranial magnetic stimulation (TMS) may be as small as a few millimeters. Despite such great potential, navigated TMS (nTMS) mapping is still underused for the assessment of motor plasticity, particularly in clinical settings. Here, we investigate the within-limb somatotopy gradient as well as absolute and relative reliability of three hand muscle cortical representations (MCRs) using a comprehensive grid-based sulcus-informed nTMS motor mapping. We enrolled 22 young healthy male volunteers. Two nTMS mapping sessions were separated by 5-10 days. Motor evoked potentials were obtained from abductor pollicis brevis (APB), abductor digiti minimi, and extensor digitorum communis. In addition to individual MRI-based analysis, we studied normalized MNI MCRs. For the reliability assessment, we calculated intraclass correlation and the smallest detectable change. Our results revealed a somatotopy gradient reflected by APB MCR having the most lateral location. Reliability analysis showed that the commonly used metrics of MCRs, such as areas, volumes, centers of gravity (COGs), and hotspots had a high relative and low absolute reliability for all three muscles. For within-limb TMS somatotopy, the most common metrics such as the shifts between MCR COGs and hotspots had poor relative reliability. However, overlaps between different muscle MCRs were highly reliable. We, thus, provide novel evidence that inter-muscle MCR interaction can be reliably traced using MCR overlaps while shifts between the COGs and hotspots of different MCRs are not suitable for this purpose. Our results have implications for the interpretation of nTMS motor mapping results in healthy subjects and patients with neurological conditions.

Standardizing human brain parcellations.

Using brain atlases to localize regions of interest is a requirement for making neuroscientifically valid statistical inferences. These atlases, represented in volumetric or surface coordinate spaces, can describe brain topology from a variety of perspectives. Although many human brain atlases have circulated the field over the past fifty years, limited effort has been devoted to their standardization. Standardization can facilitate consistency and transparency with respect to orientation, resolution, labeling scheme, file storage format, and coordinate space designation. Our group has worked to consolidate an extensive selection of popular human brain atlases into a single, curated, open-source library, where they are stored following a standardized protocol with accompanying metadata, which can serve as the basis for future atlases. The repository containing the atlases, the specification, as well as relevant transformation functions is available in the neuroparc OSF registered repository or https://github.com/neurodata/neuroparc .

Tactile acuity of fingertips and hand representation size in human Area 3b and Area 1 of the primary somatosensory cortex.

Intracortical mapping in monkeys revealed a full body map in all four cytoarchitectonic subdivisions of the contralateral primary somatosensory cortex (S1), as well as positive associations between spatio-tactile acuity performance of the fingers and their representation field size especially within cytoarchitectonic Area 3b and Area 1. Previous non-invasive investigations on these associations in humans assumed a monotonous decrease of representation field size from index finger to little finger although the field sizes are known to change in response to training or in disease. Recent developments improved noninvasive functional mapping of S1 by a) adding a cognitive task during repetitive stimulation to decrease habituation to the stimuli, b) smaller voxel size of fMRI-sequences, c) surface-based analysis accounting for cortical curvature, and d) increase of spatial specificity for fMRI data analysis by avoidance of smoothing, partial volume effects, and pial vein signals. We here applied repetitive pneumatic stimulation of digit 1 (D1; thumb) and digit 5 (D5; little finger) on both hands to investigate finger/hand representation maps in the complete S1, but also in cytoarchitectonic Areas 1, 2, 3a, and 3b separately, in 21 healthy volunteers using 3T fMRI. The distances between activation maxima of D1 and D5 were evaluated by two independent raters, blinded for performance parameters. The fingertip representations showed a somatotopy and were localized in the transition region between the crown and the anterior wall of the post central gyrus agreeing with Area 1 and 3b. Participants were comprehensively tested for tactile performance using von Freyhair filaments to determine cutaneous sensory thresholds (CST) as well as grating orientation thresholds (GOT) and two-point resolution (TPD) for spatio-tactile acuity testing. Motor performance was evaluated with pinch grip performance (Roeder test). We found bilateral associations of D1-D5 distance for GOT thresholds and partially also for TPD in Area 3b and in Area 1, but not if using the complete S1 mask. In conclusion, we here demonstrate that 3T fMRI is capable to map associations between spatio-tactile acuity and the fingertip representation in Area 3b and Area 1 in healthy participants.

Variability and effect sizes of intracranial current source density estimations during pain: Systematic review, experimental findings, and future perspectives.

Pain arises from the integration of sensory and cognitive processes in the brain, resulting in specific patterns of neural oscillations that can be characterized by measuring electrical brain activity. Current source density (CSD) estimation from low-resolution brain electromagnetic tomography (LORETA) and its standardized (sLORETA) and exact (eLORETA) variants, is a common approach to identify the spatiotemporal dynamics of the brain sources in physiological and pathological pain-related conditions. However, there is no consensus on the magnitude and variability of clinically or experimentally relevant effects for CSD estimations. Here, we systematically examined reports of sample size calculations and effect size estimations in all studies that included the keywords pain, and LORETA, sLORETA, or eLORETA in Scopus and PubMed. We also assessed the reliability of LORETA CSD estimations during non-painful and painful conditions to estimate hypothetical sample sizes for future experiments using CSD estimations. We found that none of the studies included in the systematic review reported sample size calculations, and less than 20% reported measures of central tendency and dispersion, which are necessary to estimate effect sizes. Based on these data and our experimental results, we determined that sample sizes commonly used in pain studies using CSD estimations are suitable to detect medium and large effect sizes in crossover designs and only large effects in parallel designs. These results provide a comprehensive summary of the effect sizes observed using LORETA in pain research, and this information can be used by clinicians and researchers to improve settings and designs of future pain studies.

Localization accuracy of a common beamformer for the comparison of two conditions.

The linearly constrained minimum variance beamformer is frequently used to reconstruct sources underpinning neuromagnetic recordings. When reconstructions must be compared across conditions, it is considered good practice to use a single, "common" beamformer estimated from all the data at once. This is to ensure that differences between conditions are not ascribable to differences in beamformer weights. Here, we investigate the localization accuracy of such a common beamformer. Based on theoretical derivations, we first show that the common beamformer leads to localization errors in source reconstruction. We then turn to simulations in which we attempt to reconstruct a (genuine) source in a first condition, while considering a second condition in which there is an (interfering) source elsewhere in the brain. We estimate maps of mislocalization and assess statistically the difference between "standard" and "common" beamformers. We complement our findings with an application to experimental MEG data. The results show that the common beamformer may yield significant mislocalization. Specifically, the common beamformer may force the genuine source to be reconstructed closer to the interfering source than it really is. As the same applies to the reconstruction of the interfering source, both sources are pulled closer together than they are. This observation was further illustrated in experimental data. Thus, although the common beamformer allows for the comparison of conditions, in some circumstances it introduces localization inaccuracies. We recommend alternative approaches to the general problem of comparing conditions.

Brain topography beyond parcellations: Local gradients of functional maps.

Functional neuroimaging provides the unique opportunity to characterize brain regions based on their response to tasks or ongoing activity. As such, it holds the premise to capture brain spatial organization. Yet, the conceptual framework to describe this organization has remained elusive: on the one hand, parcellations build implicitly on a piecewise constant organization, i.e. flat regions separated by sharp boundaries; on the other hand, the recently popularized concept of functional gradient hints instead at a smooth structure. Noting that both views converge to a topographic scheme that pieces together local variations of functional features, we perform a quantitative assessment of local gradient-based models. Using as a driving case the prediction of functional Magnetic Resonance Imaging (fMRI) data -concretely, the prediction of task-fMRI from rest-fMRI maps across subjects- we develop a parcel-wise linear regression model based on a dictionary of reference topographies. Our method uses multiple random parcellations -as opposed to a single fixed parcellation- and aggregates estimates across these parcellations to predict functional features in left-out subjects. Our experiments demonstrate the existence of an optimal cardinality of the parcellation to capture local gradients of functional maps.

QSMART: Quantitative susceptibility mapping artifact reduction technique.

PURPOSE: Quantitative susceptibility mapping (QSM) is a novel MR technique that allows mapping of tissue susceptibility values from MR phase images. QSM is an ill-conditioned inverse problem, and although several methods have been proposed in the field, in the presence of a wide range of susceptibility sources, streaking artifacts appear around high susceptibility regions and contaminate the whole QSM map. QSMART is a post-processing pipeline that uses two-stage parallel inversion to reduce the streaking artifacts and remove banding artifact at the cortical surface and around the vasculature. METHOD: Tissue and vein susceptibility values were separately estimated by generating a mask of vasculature driven from the magnitude data using a Frangi filter. Spatially dependent filtering was used for the background field removal step and the two susceptibility estimates were combined in the final QSM map. QSMART was compared to RESHARP/iLSQR and V-SHARP/iLSQR inversion in a numerical phantom, 7T in vivo single and multiple-orientation scans, 9.4T ex vivo mouse data, and 4.7T in vivo rat brain with induced focal ischemia. RESULTS: Spatially dependent filtering showed better suppression of phase artifacts near cortex compared to RESHARP and V-SHARP, while preserving voxels located within regions of interest without brain edge erosion. QSMART showed successful reduction of streaking artifacts as well as improved contrast between different brain tissues compared to the QSM maps obtained by RESHARP/iLSQR and V-SHARP/iLSQR. CONCLUSION: QSMART can reduce QSM artifacts to enable more robust estimation of susceptibility values in vivo and ex vivo.