The benefits of skull stripping in the normalization of clinical fMRI data

CT-guided spatial normalization of nuclear hybrid imaging adapted to enlarged ventricles: Impact on striatal uptake quantification

INTRODUCTION

Spatial normalization is a prerequisite step for the quantitative analysis of SPECT or PET brain images using volume-of-interest (VOI) template or voxel-based analysis. MRI-guided spatial normalization is the gold standard, but the wide use of PET/CT or SPECT/CT in routine clinical practice makes CT-guided spatial normalization a necessary alternative. Ventricular enlargement is observed with aging, and it hampers the spatial normalization of the lateral ventricles and striatal regions, limiting their analysis. The aim of the present study was to propose a robust spatial normalization method based on CT scans that takes into account features of the aging brain to reduce bias in the CT-guided striatal analysis of SPECT images.

METHODS

We propose an enhanced CT-guided spatial normalization pipeline based on SPM12. Performance of the proposed pipeline was assessed on visually normal [123I]-FP-CIT SPECT/CT images. SPM12 default CT-guided spatial normalization was used as reference method. The metrics assessed were the overlap between the spatially normalized lateral ventricles and caudate/putamen VOIs, and the computation of caudate and putamen specific binding ratios (SBR).

RESULTS

In total 231 subjects (mean age ± SD = 61.9 ± 15.5 years) were included in the statistical analysis. The mean overlap between the spatially normalized lateral ventricles of subjects and the caudate VOI and the mean SBR of caudate were respectively 38.40 % (± SD = 19.48 %) of the VOI and 1.77 (± 0.79) when performing SPM12 default spatial normalization. The mean overlap decreased to 9.13 % (± SD = 1.41 %, P < 0.001) of the VOI and the SBR of caudate increased to 2.38 (± 0.51, P < 0.0001) when performing the proposed pipeline. Spatially normalized lateral ventricles did not overlap with putamen VOI using either method. The mean putamen SBR value derived from the proposed spatial normalization (2.75 ± 0.54) was not significantly different from that derived from the default SPM12 spatial normalization (2.83 ± 0.52, P > 0.05).

CONCLUSION

The automatic CT-guided spatial normalization used herein led to a less biased spatial normalization of SPECT images, hence an improved semi-quantitative analysis. The proposed pipeline could be implemented in clinical routine to perform a more robust SBR computation using hybrid imaging.

Anatomy-aware and acquisition-agnostic joint registration with SynthMorph

Affine image registration is a cornerstone of medical-image analysis. While classical algorithms can achieve excellent accuracy, they solve a time-consuming optimization for every image pair. Deep-learning (DL) methods learn a function that maps an image pair to an output transform. Evaluating the function is fast, but capturing large transforms can be challenging, and networks tend to struggle if a test-image characteristic shifts from the training domain, such as the resolution. Most affine methods are agnostic to the anatomy the user wishes to align, meaning the registration will be inaccurate if algorithms consider all structures in the image. We address these shortcomings with SynthMorph, a fast, symmetric, diffeomorphic, and easy-to-use DL tool for joint affine-deformable registration of any brain image without preprocessing. First, we leverage a strategy that trains networks with widely varying images synthesized from label maps, yielding robust performance across acquisition specifics unseen at training. Second, we optimize the spatial overlap of select anatomical labels. This enables networks to distinguish anatomy of interest from irrelevant structures, removing the need for preprocessing that excludes content which would impinge on anatomy-specific registration. Third, we combine the affine model with a deformable hypernetwork that lets users choose the optimal deformation-field regularity for their specific data, at registration time, in a fraction of the time required by classical methods. This framework is applicable to learning anatomy-aware, acquisition-agnostic registration of any anatomy with any architecture, as long as label maps are available for training. We analyze how competing architectures learn affine transforms and compare state-of-the-art registration tools across an extremely diverse set of neuroimaging data, aiming to truly capture the behavior of methods in the real world. SynthMorph demonstrates high accuracy and is available at https://w3id.org/synthmorph, as a single complete end-to-end solution for registration of brain magnetic resonance imaging (MRI) data.

Recognizing Pediatric Tuberous Sclerosis Complex Based on Multi-Contrast MRI and Deep Weighted Fusion Network

Multi-contrast magnetic resonance imaging (MRI) is wildly applied to identify tuberous sclerosis complex (TSC) children in a clinic. In this work, a deep convolutional neural network with multi-contrast MRI is proposed to diagnose pediatric TSC. Firstly, by combining T2W and FLAIR images, a new synthesis modality named FLAIR3 was created to enhance the contrast between TSC lesions and normal brain tissues. After that, a deep weighted fusion network (DWF-net) using a late fusion strategy is proposed to diagnose TSC children. In experiments, a total of 680 children were enrolled, including 331 healthy children and 349 TSC children. The experimental results indicate that FLAIR3 successfully enhances the visibility of TSC lesions and improves the classification performance. Additionally, the proposed DWF-net delivers a superior classification performance compared to previous methods, achieving an AUC of 0.998 and an accuracy of 0.985. The proposed method has the potential to be a reliable computer-aided diagnostic tool for assisting radiologists in diagnosing TSC children.

Automatic brain extraction and brain tissues segmentation on multi-contrast animal MRI

For many neuroscience applications, brain extraction in MRI images is the first pre-processing step of a quantification pipeline. Once the brain is extracted, further post-processing calculations become faster, more specific and easier to implement and interpret. It is the case, for example, of functional MRI brain studies, or relaxation time mappings and brain tissues classifications to characterise brain pathologies. Existing brain extraction tools are mostly adapted to work on the human anatomy, this gives poor results when applied to animal brain images. We have developed an atlas-based Veterinary Images Brain Extraction (VIBE) algorithm that encompasses a pre-processing step to adapt the atlas to the patient's image, and a subsequent registration step. We show that the brain extraction is achieved with excellent results in terms of Dice and Jaccard metrics. The algorithm is automatic, with no need to adapt the parameters in a broad range of situations: we successfully tested multiple MRI contrasts (T1-weighted, T2-weighted, T2-weighted FLAIR), all the acquisition planes (sagittal, dorsal, transverse), different animal species (dogs and cats) and canine cranial conformations (brachycephalic, mesocephalic, dolichocephalic). VIBE can be successfully extended to other animal species, provided that an atlas for that specific species exists. We show also how brain extraction, as a preliminary step, can help to segment brain tissues with a K-Means clustering algorithm.

A functional MRI pre-processing and quality control protocol based on statistical parametric mapping (SPM) and MATLAB

Functional MRI (fMRI) has become a popular technique to study brain functions and their alterations in psychiatric and neurological conditions. The sample sizes for fMRI studies have been increasing steadily, and growing studies are sourced from open-access brain imaging repositories. Quality control becomes critical to ensure successful data processing and valid statistical results. Here, we outline a simple protocol for fMRI data pre-processing and quality control based on statistical parametric mapping (SPM) and MATLAB. The focus of this protocol is not only to identify and remove data with artifacts and anomalies, but also to ensure the processing has been performed properly. We apply this protocol to the data from fMRI Open quality control (QC) Project, and illustrate how each quality control step can help to identify potential issues. We also show that simple steps such as skull stripping can improve coregistration between the functional and anatomical images.

SynthStrip: skull-stripping for any brain image

The removal of non-brain signal from magnetic resonance imaging (MRI) data, known as skull-stripping, is an integral component of many neuroimage analysis streams. Despite their abundance, popular classical skull-stripping methods are usually tailored to images with specific acquisition properties, namely near-isotropic resolution and T1-weighted (T1w) MRI contrast, which are prevalent in research settings. As a result, existing tools tend to adapt poorly to other image types, such as stacks of thick slices acquired with fast spin-echo (FSE) MRI that are common in the clinic. While learning-based approaches for brain extraction have gained traction in recent years, these methods face a similar burden, as they are only effective for image types seen during the training procedure. To achieve robust skull-stripping across a landscape of imaging protocols, we introduce SynthStrip, a rapid, learning-based brain-extraction tool. By leveraging anatomical segmentations to generate an entirely synthetic training dataset with anatomies, intensity distributions, and artifacts that far exceed the realistic range of medical images, SynthStrip learns to successfully generalize to a variety of real acquired brain images, removing the need for training data with target contrasts. We demonstrate the efficacy of SynthStrip for a diverse set of image acquisitions and resolutions across subject populations, ranging from newborn to adult. We show substantial improvements in accuracy over popular skull-stripping baselines - all with a single trained model. Our method and labeled evaluation data are available at https://w3id.org/synthstrip.

A fuzzy membership based comparison of the grey matter (GM) in cognitively normal (CN), mild cognitive impairment (MCI), and Alzheimer’s disease (AD) using brain images

Grey matter (GM) in human brain contains most of the important cells covering the regions involved in neurophysiological operations such as memory, emotions, decision making, etc. Alzheimer’s disease (AD) is a neurological disease that kills the brain cells in regions which are mostly involved in the neurophysiological operations. Mild Cognitive Impairment (MCI) is a stage between Cognitively Normal (CN) and AD, where a significant cognitive declination can be observed. The destruction of brain cells causes a reduction in the size of GM. Evaluation of changes in GM, may help in studying the overall brain transformations and accurate classification of different stages of AD. In this work, firstly skull of brain images is stripped for 5 different slices, then segmentation of GM is performed. Finally, the average number of pixels in grey region and the average atrophy in grey pixels per year is calculated and compared amongst CN, MCI, and AD patients of various ages and genders. It is observed that, for some subjects (in some particular ages) from different dementia stages, pattern of GM changes is almost identical. To solve this issue, we have used the concept of fuzzy membership functions to classify the dementia stages more accurately. It is observed from the comparison that average difference in the number of pixels between CN and MCI= 10.01%, CN and AD= 19.63%, MCI and AD= 10.72%. It can be also observed from the comparison that, the average atrophy in grey matter per year in CN= 1.92%, MCI= 3.13%, and AD= 4.33%.

Non-parametric combination of multimodal MRI for lesion detection in focal epilepsy

•

Multivariate voxel-based analysis useful for lesion detection in focal epilepsy.

•

Non-parametric combination algorithm used to combine data from various MR sequences.

•

Successful lesion detection demonstrated in MRI-positive and MRI-negative patients.

•

Multimodal analysis detected abnormalities from diverse epileptogenic pathologies.

•

Sensitivity of multivariate analysis notably higher than univariate analyses.

One third of patients with medically refractory focal epilepsy have normal-appearing MRI scans. This poses a problem as identification of the epileptogenic region is required for surgical treatment. This study performs a multimodal voxel-based analysis (VBA) to identify brain abnormalities in MRI-negative focal epilepsy. Data was collected from 69 focal epilepsy patients (42 with discrete lesions on MRI scans, 27 with no visible findings on scans), and 62 healthy controls. MR images comprised T1-weighted, fluid-attenuated inversion recovery (FLAIR), fractional anisotropy (FA) and mean diffusivity (MD) from diffusion tensor imaging, and neurite density index (NDI) from neurite orientation dispersion and density imaging. These multimodal images were coregistered to T1-weighted scans, normalized to a standard space, and smoothed with 8 mm FWHM. Initial analysis performed voxel-wise one-tailed t-tests separately on grey matter concentration (GMC), FLAIR, FA, MD, and NDI, comparing patients with epilepsy to controls. A multimodal non-parametric combination (NPC) analysis was also performed simultaneously on FLAIR, FA, MD, and NDI. Resulting p-maps were family-wise error rate corrected, threshold-free cluster enhanced, and thresholded at p < 0.05. Sensitivity was established through visual comparison of results to manually drawn lesion masks or seizure onset zone (SOZ) from stereoelectroencephalography. A leave-one-out cross-validation with the same analysis protocols was performed on controls to determine specificity. NDI was the best performing individual modality, detecting focal abnormalities in 38% of patients with normal MRI and conclusive SOZ. GMC demonstrated the lowest sensitivity at 19%. NPC provided superior performance to univariate analyses with 50% sensitivity. Specificity in controls ranged between 96 and 100% for all analyses. This study demonstrated the utility of a multimodal VBA utilizing NPC for detecting epileptogenic lesions in MRI-negative focal epilepsy. Future work will apply this approach to datasets from other centres and will experiment with different combinations of MR sequences.

Anxiety and the Neurobiology of Temporally Uncertain Threat Anticipation

When extreme, anxiety-a state of distress and arousal prototypically evoked by uncertain danger-can be debilitating. Uncertain anticipation is a shared feature of situations that elicit signs and symptoms of anxiety across psychiatric disorders, species, and assays. Despite the profound significance of anxiety for human health and wellbeing, the neurobiology of uncertain-threat anticipation remains unsettled. Leveraging a paradigm adapted from animal research and optimized for fMRI signal decomposition, we examined the neural circuits engaged during the anticipation of temporally uncertain and certain threat in 99 men and women. Results revealed that the neural systems recruited by uncertain and certain threat anticipation are anatomically colocalized in frontocortical regions, extended amygdala, and periaqueductal gray. Comparison of the threat conditions demonstrated that this circuitry can be fractionated, with frontocortical regions showing relatively stronger engagement during the anticipation of uncertain threat, and the extended amygdala showing the reverse pattern. Although there is widespread agreement that the bed nucleus of the stria terminalis and dorsal amygdala-the two major subdivisions of the extended amygdala-play a critical role in orchestrating adaptive responses to potential danger, their precise contributions to human anxiety have remained contentious. Follow-up analyses demonstrated that these regions show statistically indistinguishable responses to temporally uncertain and certain threat anticipation. These observations provide a framework for conceptualizing anxiety and fear, for understanding the functional neuroanatomy of threat anticipation in humans, and for accelerating the development of more effective intervention strategies for pathological anxiety.SIGNIFICANCE STATEMENT Anxiety-an emotion prototypically associated with the anticipation of uncertain harm-has profound significance for public health, yet the underlying neurobiology remains unclear. Leveraging a novel neuroimaging paradigm in a relatively large sample, we identify a core circuit responsive to both uncertain and certain threat anticipation, and show that this circuitry can be fractionated into subdivisions with a bias for one kind of threat or the other. The extended amygdala occupies center stage in neuropsychiatric models of anxiety, but its functional architecture has remained contentious. Here we demonstrate that its major subdivisions show statistically indistinguishable responses to temporally uncertain and certain threat. Collectively, these observations indicate the need to revise how we think about the neurobiology of anxiety and fear.

Brain compensation during response inhibition in premanifest Huntington's disease

Premanifest Huntington's disease (pre-HD) individuals typically show increased task-related functional magnetic resonance imaging (fMRI), suggested to reflect compensatory strategies. Despite the evidence, no study has attempted to understand the compensatory process in light of 'formal' models of compensation. We used a quantitative model of compensation - the Compensation-Related Utilization of Neural Circuits Hypothesis (CRUNCH), to characterise compensation in pre-HD using fMRI. Pre-HD individuals (n = 15) and controls (n = 15) performed a modified stop-signal task that incremented in four levels of stop difficulty. Our results did not support the critical assumption of the CRUNCH model - controls did not show increased fMRI activity with increased level of stop difficulty; however, controls showed decreased fMRI activity with increased stop difficulty in right inferior frontal gyrus and right caudate nucleus. Relative to controls, pre-HD individuals showed increased fMRI activity in right inferior frontal gyrus and in right caudate nucleus at higher levels of stop difficulty, which is the opposite effect to that predicted by the model. Our findings suggest a compensatory process of the response inhibition network in pre-HD; however, the pattern of fMRI activity was not in the manner expected by CRUNCH.

Conventional and Deep Learning Methods for Skull Stripping in Brain MRI

Skull stripping in brain magnetic resonance volume has recently been attracting attention due to an increased demand to develop an efficient, accurate, and general algorithm for diverse datasets of the brain. Accurate skull stripping is a critical step for neuroimaging diagnostic systems because neither the inclusion of non-brain tissues nor removal of brain parts can be corrected in subsequent steps, which results in unfixed error through subsequent analysis. The objective of this review article is to give a comprehensive overview of skull stripping approaches, including recent deep learning-based approaches. In this paper, the current methods of skull stripping have been divided into two distinct groups—conventional or classical approaches, and convolutional neural networks or deep learning approaches. The potentials of several methods are emphasized because they can be applied to standard clinical imaging protocols. Finally, current trends and future developments are addressed giving special attention to recent deep learning algorithms.

Quantitative brain relaxation atlases for personalized detection and characterization of brain pathology

PURPOSE

To exploit the improved comparability and hardware independency of quantitative MRI, databases of MR physical parameters in healthy tissue are required, to which tissue properties of patients can be compared. In this work, normative values for longitudinal and transverse relaxation times in the brain were established and tested in single-subject comparisons for detection of abnormal relaxation times.

METHODS

Relaxometry maps of the brain were acquired from 52 healthy volunteers. After spatially normalizing the volumes into a common space, T₁ and T₂ inter-subject variability within the healthy cohort was modeled voxel-wise. A method for a single-subject comparison against the atlases was developed by computing z-scores with respect to the established healthy norms. The comparison was applied to two multiple sclerosis and one clinically isolated syndrome cases for a proof of concept.

RESULTS

The established atlases exhibit a low variation in white matter structures (median RMSE of models equal to 32 ms for T₁ and 4 ms for T₂ ), indicating that relaxation times are in a narrow range for normal tissues. The proposed method for single-subject comparison detected relaxation time deviations from healthy norms in the example patient data sets. Relaxation times were found to be increased in brain lesions (mean z-scores >5). Moreover, subtle and confluent differences (z-scores ~2-4) were observed in clinically plausible regions (between lesions, corpus callosum).

CONCLUSIONS

Brain T₁ and T₂ quantitative norms were derived voxel-wise with low variability in healthy tissue. Example patient deviation maps demonstrated good sensitivity of the atlases for detecting relaxation time alterations.

A Spatial Registration Toolbox for Structural MR Imaging of the Aging Brain

During aging the brain undergoes a series of structural changes, in size, shape as well as tissue composition. In particular, cortical atrophy and ventricular enlargement are often present in the brain of elderly individuals. This poses serious challenges in the spatial registration of structural MR images. In this study, we addressed this open issue by proposing an enhanced framework for MR registration and segmentation. Our solution was compared with other approaches based on the tools available in SPM12, a widely used software package. Performance of the different methods was assessed on 229 T1-weighted images collected in healthy individuals, with age ranging between 55 and 90 years old. Our method showed a consistent improvement as compared to other solutions, especially for subjects with enlarged lateral ventricles. It also provided a superior inter-subject alignment in cortical regions, with the most marked improvement in the frontal lobe. We conclude that our method is a valid alternative to standard approaches based on SPM12, and is particularly suitable for the processing of structural MR images of brains with cortical atrophy and ventricular enlargement. The method is integrated in our software toolbox MRTool, which is freely available to the scientific community.

fMRI data processing in MRTOOL: to what extent does anatomical registration affect the reliability of functional results?

Spatial registration is an essential step in the analysis of fMRI data because it enables between-subject analyses of brain activity, measured either during task performance or in the resting state. In this study, we investigated how anatomical registration with MRTOOL affects the reliability of task-related fMRI activity. We used as a benchmark the results from two other spatial registration methods implemented in SPM12: the Unified Segmentation algorithm and the DARTEL toolbox. Structural alignment accuracy and the impact on functional activation maps were assessed with high-resolution T1-weighted images and a set of task-related functional volumes acquired in 10 healthy volunteers. Our findings confirmed that anatomical registration is a crucial step in fMRI data processing, contributing significantly to the total inter-subject variance of the activation maps. MRTOOL and DARTEL provided greater registration accuracy than Unified Segmentation. Although DARTEL had superior gray matter and white matter tissue alignment than MRTOOL, there were no significant differences between DARTEL and MRTOOL in test-retest reliability. Likewise, we found only limited differences in BOLD activation morphology between MRTOOL and DARTEL. The test-retest reliability of task-related responses was comparable between MRTOOL and DARTEL, and both proved superior to Unified Segmentation. We conclude that MRTOOL, which is suitable for single-subject processing of structural and functional MR images, is a valid alternative to other SPM12-based approaches that are intended for group analysis. MRTOOL now includes a normalization module for fMRI data and is freely available to the scientific community.

Functional Connectivity within the Primate Extended Amygdala Is Heritable and Associated with Early-Life Anxious Temperament

Children with an extremely inhibited, anxious temperament (AT) are at increased risk for anxiety disorders and depression. Using a rhesus monkey model of early-life AT, we previously demonstrated that metabolism in the central extended amygdala (EAc), including the central nucleus of the amygdala (Ce) and bed nucleus of the stria terminalis (BST), is associated with trait-like variation in AT. Here, we use fMRI to examine relationships between Ce-BST functional connectivity and AT in a large multigenerational family pedigree of rhesus monkeys (n = 170 females and 208 males). Results demonstrate that Ce-BST functional connectivity is heritable, accounts for a significant but modest portion of the variance in AT, and is coheritable with AT. Interestingly, Ce-BST functional connectivity and AT-related BST metabolism were not correlated and accounted for non-overlapping variance in AT. Exploratory analyses suggest that Ce-BST functional connectivity is associated with metabolism in the hypothalamus and periaqueductal gray. Together, these results suggest the importance of coordinated function within the EAc for determining individual differences in AT and metabolism in brain regions associated with its behavioral and neuroendocrine components.SIGNIFICANCE STATEMENT Anxiety disorders directly impact the lives of nearly one in five people, accounting for substantial worldwide suffering and disability. Here, we use a nonhuman primate model of anxious temperament (AT) to understand the neurobiology underlying the early-life risk to develop anxiety disorders. Leveraging the same kinds of neuroimaging measures routinely used in human studies, we demonstrate that coordinated activation between the central nucleus of the amygdala and the bed nucleus of the stria terminalis is correlated with, and coinherited with, early-life AT. Understanding how these central extended amygdala regions work together to produce extreme anxiety provides a neural target for early-life interventions with the promise of preventing lifelong disability in at-risk children.

Intrinsic functional connectivity of the central extended amygdala

The central extended amygdala (EAc)-including the bed nucleus of the stria terminalis (BST) and central nucleus of the amygdala (Ce)-plays a critical role in triggering fear and anxiety and is implicated in the development of a range of debilitating neuropsychiatric disorders. Although it is widely believed that these disorders reflect the coordinated activity of distributed neural circuits, the functional architecture of the EAc network and the degree to which the BST and the Ce show distinct patterns of functional connectivity is unclear. Here, we used a novel combination of imaging approaches to trace the connectivity of the BST and the Ce in 130 healthy, racially diverse, community-dwelling adults. Multiband imaging, high-precision registration techniques, and spatially unsmoothed data maximized anatomical specificity. Using newly developed seed regions, whole-brain regression analyses revealed robust functional connectivity between the BST and Ce via the sublenticular extended amygdala, the ribbon of subcortical gray matter encompassing the ventral amygdalofugal pathway. Both regions displayed coupling with the ventromedial prefrontal cortex (vmPFC), midcingulate cortex (MCC), insula, and anterior hippocampus. The BST showed stronger connectivity with the thalamus, striatum, periaqueductal gray, and several prefrontal territories. The only regions showing stronger functional connectivity with the Ce were neighboring regions of the dorsal amygdala, amygdalohippocampal area, and anterior hippocampus. These observations provide a baseline against which to compare a range of special populations, inform our understanding of the role of the EAc in normal and pathological fear and anxiety, and showcase image registration techniques that are likely to be useful for researchers working with "deidentified" neuroimaging data.

Comparison of arterial spin labeling registration strategies in the multi-center GENetic frontotemporal dementia initiative (GENFI)

PURPOSE

To compare registration strategies to align arterial spin labeling (ASL) with 3D T1-weighted (T1w) images, with the goal of reducing the between-subject variability of cerebral blood flow (CBF) images.

MATERIALS AND METHODS

Multi-center 3T ASL data were collected at eight sites with four different sequences in the multi-center GENetic Frontotemporal dementia Initiative (GENFI) study. In a total of 48 healthy controls, we compared the following image registration options: (I) which images to use for registration (perfusion-weighted images [PWI] to the segmented gray matter (GM) probability map (pGM) (CBF-pGM) or M0 to T1w (M0-T1w); (II) which transformation to use (rigid-body or non-rigid); and (III) whether to mask or not (no masking, M0-based FMRIB software library Brain Extraction Tool [BET] masking). In addition to visual comparison, we quantified image similarity using the Pearson correlation coefficient (CC), and used the Mann-Whitney U rank sum test.

RESULTS

CBF-pGM outperformed M0-T1w (CC improvement 47.2% ± 22.0%; P < 0.001), and the non-rigid transformation outperformed rigid-body (20.6% ± 5.3%; P < 0.001). Masking only improved the M0-T1w rigid-body registration (14.5% ± 15.5%; P = 0.007).

CONCLUSION

The choice of image registration strategy impacts ASL group analyses. The non-rigid transformation is promising but requires validation. CBF-pGM rigid-body registration without masking can be used as a default strategy. In patients with expansive perfusion deficits, M0-T1w may outperform CBF-pGM in sequences with high effective spatial resolution. BET-masking only improves M0-T1w registration when the M0 image has sufficient contrast.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:131-140.

Robust presurgical functional MRI at 7 T using response consistency

Functional MRI is valuable in presurgical planning due to its non-invasive nature, repeatability, and broad availability. Using ultra-high field MRI increases the specificity and sensitivity, increasing the localization reliability and reducing scan time. Ideally, fMRI analysis for this application should identify unreliable runs and work even if the patient deviates from the prescribed task timing or if there are changes to the hemodynamic response due to pathology. In this study, a model-free analysis method-UNBIASED-based on the consistency of fMRI responses over runs was applied, to ultra-high field fMRI localizations of the hand area. Ten patients with brain tumors and epilepsy underwent 7 Tesla fMRI with multiple runs of a hand motor task in a block design. FMRI data were analyzed with the proposed approach (UNBIASED) and the conventional General Linear Model (GLM) approach. UNBIASED correctly identified and excluded fMRI runs that contained little or no activation. Generally, less motion artifact contamination was present in UNBIASED than in GLM results. Some cortical regions were identified as activated in UNBIASED but not GLM results. These were confirmed to show reproducible delayed or transient activation, which was time-locked to the task. UNBIASED is a robust approach to generating activation maps without the need for assumptions about response timing or shape. In presurgical planning, UNBIASED can complement model-based methods to aid surgeons in making prudent choices about optimal surgical access and resection margins for each patient, even if the hemodynamic response is modified by pathology. Hum Brain Mapp 38:3163-3174, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

CerebroMatic: A Versatile Toolbox for Spline-Based MRI Template Creation

Brain image spatial normalization and tissue segmentation rely on prior tissue probability maps. Appropriately selecting these tissue maps becomes particularly important when investigating "unusual" populations, such as young children or elderly subjects. When creating such priors, the disadvantage of applying more deformation must be weighed against the benefit of achieving a crisper image. We have previously suggested that statistically modeling demographic variables, instead of simply averaging images, is advantageous. Both aspects (more vs. less deformation and modeling vs. averaging) were explored here. We used imaging data from 1914 subjects, aged 13 months to 75 years, and employed multivariate adaptive regression splines to model the effects of age, field strength, gender, and data quality. Within the spm/cat12 framework, we compared an affine-only with a low- and a high-dimensional warping approach. As expected, more deformation on the individual level results in lower group dissimilarity. Consequently, effects of age in particular are less apparent in the resulting tissue maps when using a more extensive deformation scheme. Using statistically-described parameters, high-quality tissue probability maps could be generated for the whole age range; they are consistently closer to a gold standard than conventionally-generated priors based on 25, 50, or 100 subjects. Distinct effects of field strength, gender, and data quality were seen. We conclude that an extensive matching for generating tissue priors may model much of the variability inherent in the dataset which is then not contained in the resulting priors. Further, the statistical description of relevant parameters (using regression splines) allows for the generation of high-quality tissue probability maps while controlling for known confounds. The resulting CerebroMatic toolbox is available for download at http://irc.cchmc.org/software/cerebromatic.php.

A Hitchhiker's Guide to Functional Magnetic Resonance Imaging

Functional Magnetic Resonance Imaging (fMRI) studies have become increasingly popular both with clinicians and researchers as they are capable of providing unique insights into brain functions. However, multiple technical considerations (ranging from specifics of paradigm design to imaging artifacts, complex protocol definition, and multitude of processing and methods of analysis, as well as intrinsic methodological limitations) must be considered and addressed in order to optimize fMRI analysis and to arrive at the most accurate and grounded interpretation of the data. In practice, the researcher/clinician must choose, from many available options, the most suitable software tool for each stage of the fMRI analysis pipeline. Herein we provide a straightforward guide designed to address, for each of the major stages, the techniques, and tools involved in the process. We have developed this guide both to help those new to the technique to overcome the most critical difficulties in its use, as well as to serve as a resource for the neuroimaging community.

Lawrence K. Mapping Long-Term Functional Changes in Cerebral Blood Flow by Arterial Spin Labeling

Although arterial spin labeling (ASL) is appealing for mapping long-term changes in functional activity, inter-sessional variations in basal blood flow, arterial transit times (ATTs), and alignment errors, can result in significant false activation when comparing images from separate sessions. By taking steps to reduce these sources of noise, this study assessed the ability of ASL to detect functional CBF changes between sessions. ASL data were collected in three sessions to image ATT, resting CBF and CBF changes associated with motor activation (7 participants). Activation maps were generated using rest and task images acquired in the same session and from sessions separated by up to a month. Good agreement was found when comparing between-session activation maps to within-session activation maps with only a 16% decrease in precision (within-session: 90 ± 7%) and a 13% decrease in the Dice similarity (within-session: 0.75 ± 0.07) coefficient after a month. In addition, voxel-wise reproducibility (within-session: 4.7 ± 4.5%) and reliability (within-session: 0.89 ± 0.20) of resting grey-matter CBF decreased by less than 18% for the between-session analysis relative to within-session values. ATT variability between sessions (5.0 ± 2.7%) was roughly half the between-subject variability, indicating that its effects on longitudinal CBF were minimal. These results demonstrate that conducting voxel-wise analysis on CBF images acquired on different days is feasible with only modest loss in precision, highlighting the potential of ASL for longitudinal studies.

Self-similarity and recursion as default modes in human cognition

Humans generate recursive hierarchies in a variety of domains, including linguistic, social and visuo-spatial modalities. The ability to represent recursive structures has been hypothesized to increase the efficiency of hierarchical processing. Theoretical work together with recent empirical findings suggests that the ability to represent the self-similar structure of hierarchical recursive stimuli may be supported by internal neural representations that compress raw external information and increase efficiency. In order to explicitly test whether the representation of recursive hierarchies depends on internalized rules we compared the processing of visual hierarchies represented either as recursive or non-recursive, using task-free resting-state fMRI data. We aimed to evaluate the relationship between task-evoked functional networks induced by cognitive representations with the corresponding resting-state architecture. We observed increased connectivity within Default Mode Network (DMN) related brain areas during the representation of recursion, while non-recursive representations yielded increased connectivity within the Fronto-Parietal Control-Network. Our results suggest that human hierarchical information processing using recursion is supported by the DMN. In particular, the representation of recursion seems to constitute an internally-biased mode of information-processing that is mediated by both the core and dorsal-medial subsystems of the DMN. Compressed internal rule representations mediated by the DMN may help humans to represent and process hierarchical structures in complex environments by considerably reducing information processing load.

Improving the clinical potential of ultra-high field fMRI using a model-free analysis method based on response consistency

Objective

To develop an analysis method that is sensitive to non-model-conform responses often encountered in ultra-high field presurgical planning fMRI. Using the consistency of time courses over a number of experiment repetitions, it should exclude low quality runs and generate activation maps that reflect the reliability of responses.

Materials and methods

7 T fMRI data were acquired from six healthy volunteers: three performing purely motor tasks and three a visuomotor task. These were analysed with the proposed approach (UNBIASED) and the GLM.

Results

UNBIASED results were generally less affected by false positive results than the GLM. Runs that were identified as being of low quality were confirmed to contain little or no activation. In two cases, regions were identified as activated in UNBIASED but not GLM results. Signal changes in these areas were time-locked to the task, but were delayed or transient.

Conclusion

UNBIASED is shown to be a reliable means of identifying consistent task-related signal changes regardless of response timing. In presurgical planning, UNBIASED could be used to rapidly generate reliable maps of the consistency with which eloquent brain regions are activated without recourse to task timing and despite modified hemodynamics.

Electronic supplementary material

The online version of this article (doi:10.1007/s10334-016-0533-8) contains supplementary material, which is available to authorized users.

Feasibility of Computed Tomography-Guided Methods for Spatial Normalization of Dopamine Transporter Positron Emission Tomography Image

Background

Spatial normalization is a prerequisite step for analyzing positron emission tomography (PET) images both by using volume-of-interest (VOI) template and voxel-based analysis. Magnetic resonance (MR) or ligand-specific PET templates are currently used for spatial normalization of PET images. We used computed tomography (CT) images acquired with PET/CT scanner for the spatial normalization for [¹⁸F]-N-3-fluoropropyl-2-betacarboxymethoxy-3-beta-(4-iodophenyl) nortropane (FP-CIT) PET images and compared target-to-cerebellar standardized uptake value ratio (SUVR) values with those obtained from MR- or PET-guided spatial normalization method in healthy controls and patients with Parkinson’s disease (PD).

Methods

We included 71 healthy controls and 56 patients with PD who underwent [¹⁸F]-FP-CIT PET scans with a PET/CT scanner and T1-weighted MR scans. Spatial normalization of MR images was done with a conventional spatial normalization tool (cvMR) and with DARTEL toolbox (dtMR) in statistical parametric mapping software. The CT images were modified in two ways, skull-stripping (ssCT) and intensity transformation (itCT). We normalized PET images with cvMR-, dtMR-, ssCT-, itCT-, and PET-guided methods by using specific templates for each modality and measured striatal SUVR with a VOI template. The SUVR values measured with FreeSurfer-generated VOIs (FSVOI) overlaid on original PET images were also used as a gold standard for comparison.

Results

The SUVR values derived from all four structure-guided spatial normalization methods were highly correlated with those measured with FSVOI (P < 0.0001). Putaminal SUVR values were highly effective for discriminating PD patients from controls. However, the PET-guided method excessively overestimated striatal SUVR values in the PD patients by more than 30% in caudate and putamen, and thereby spoiled the linearity between the striatal SUVR values in all subjects and showed lower disease discrimination ability. Two CT-guided methods showed comparable capability with the MR-guided methods in separating PD patients from controls and showed better correlation between putaminal SUVR values and the parkinsonian motor severity than the PET-guided method.

Conclusion

CT-guided spatial normalization methods provided reliable striatal SUVR values comparable to those obtained with MR-guided methods. CT-guided methods can be useful for analyzing dopamine transporter PET images when MR images are unavailable.

Stochastic process for white matter injury detection in preterm neonates

Preterm births are rising in Canada and worldwide. As clinicians strive to identify preterm neonates at greatest risk of significant developmental or motor problems, accurate predictive tools are required. Infants at highest risk will be able to receive early developmental interventions, and will also enable clinicians to implement and evaluate new methods to improve outcomes. While severe white matter injury (WMI) is associated with adverse developmental outcome, more subtle injuries are difficult to identify and the association with later impairments remains unknown. Thus, our goal was to develop an automated method for detection and visualization of brain abnormalities in MR images acquired in very preterm born neonates. We have developed a technique to detect WMI in T1-weighted images acquired in 177 very preterm born infants (24–32 weeks gestation). Our approach uses a stochastic process that estimates the likelihood of intensity variations in nearby pixels; with small variations being more likely than large variations. We first detect the boundaries between normal and injured regions of the white matter. Following this we use a measure of pixel similarity to identify WMI regions. Our algorithm is able to detect WMI in all of the images in the ground truth dataset with some false positives in situations where the white matter region is not segmented accurately.

A computed tomography-based spatial normalization for the analysis of [18F] fluorodeoxyglucose positron emission tomography of the brain

Objective

We developed a new computed tomography (CT)-based spatial normalization method and CT template to demonstrate its usefulness in spatial normalization of positron emission tomography (PET) images with [¹⁸F] fluorodeoxyglucose (FDG) PET studies in healthy controls.

Materials and Methods

Seventy healthy controls underwent brain CT scan (120 KeV, 180 mAs, and 3 mm of thickness) and [¹⁸F] FDG PET scans using a PET/CT scanner. T1-weighted magnetic resonance (MR) images were acquired for all subjects. By averaging skull-stripped and spatially-normalized MR and CT images, we created skull-stripped MR and CT templates for spatial normalization. The skull-stripped MR and CT images were spatially normalized to each structural template. PET images were spatially normalized by applying spatial transformation parameters to normalize skull-stripped MR and CT images. A conventional perfusion PET template was used for PET-based spatial normalization. Regional standardized uptake values (SUV) measured by overlaying the template volume of interest (VOI) were compared to those measured with FreeSurfer-generated VOI (FSVOI).

Results

All three spatial normalization methods underestimated regional SUV values by 0.3-20% compared to those measured with FSVOI. The CT-based method showed slightly greater underestimation bias. Regional SUV values derived from all three spatial normalization methods were correlated significantly (p < 0.0001) with those measured with FSVOI.

Conclusion

CT-based spatial normalization may be an alternative method for structure-based spatial normalization of [¹⁸F] FDG PET when MR imaging is unavailable. Therefore, it is useful for PET/CT studies with various radiotracers whose uptake is expected to be limited to specific brain regions or highly variable within study population.