Multimodal MEMPRAGE, FLAIR, and [Formula: see text] Segmentation to Resolve Dura and Vessels from Cortical Gray Matter

Brains in space: impact of microgravity and cosmic radiation on the CNS during space exploration

Solar system exploration is a grand endeavour of humankind. Space agencies have been planning crewed missions to the Moon and Mars for several decades. However, several environmental stress factors in space, such as microgravity and cosmic radiation, confer health risks for human explorers. This Review examines the effects of microgravity and exposure to cosmic radiation on the CNS. Microgravity presents challenges for the brain, necessitating the development of adaptive movement and orientation strategies to cope with alterations in sensory information. Exposure to microgravity also affects cognitive function to a certain extent. Recent MRI results show that microgravity affects brain structure and function. Post-flight recovery from these changes is gradual, with some lasting up to a year. Regarding cosmic radiation, animal experiments suggest that the brain could be much more sensitive to this stressor than may be expected from experiences on Earth. This may be due to the presence of energetic heavy ions in space that have an impact on cognitive function, even at low doses. However, all data about space radiation risk stem from rodent experiments, and extrapolation of these data to humans carries a high degree of uncertainty. Here, after presenting an overview of current knowledge in the above areas, we provide a concise description of possible counter-measures to protect the brain against microgravity and cosmic radiation during future space missions.

Brain Morphometrics Correlations With Age Among 350 Participants Imaged With Both 3T and 7T MRI: 7T Improves Statistical Power and Reduces Required Sample Size

Magnetic resonance imaging (MRI) at 7T has a superior signal-to-noise ratio to 3T but also presents higher signal inhomogeneities and geometric distortions. A key knowledge gap is to robustly investigate the sensitivity and accuracy of 3T and 7T MRI in assessing brain morphometrics. This study aims to (a) aggregate a large number of paired 3T and 7T scans to evaluate their differences in quantitative brain morphological assessment using a widely available brain segmentation tool, FreeSurfer, as well as to (b) examine the impact of normalization methods for subject variability and smaller sample sizes on data analysis. A total of 401 healthy participants aged 29-68 were imaged at both 3T and 7T. Structural T1-weighted magnetization-prepared rapid gradient-echo (MPRAGE) images were processed and segmented using FreeSurfer. To account for head size variability, the brain volumes underwent intracranial volume (ICV) correction using the Residual (regression model) and Proportional (simple division to ICV) methods. The resulting volumes and thicknesses were correlated with age using Pearson's correlation and false discovery rate correction. The correlations were also calculated in increasing sample size from three to the whole sample to estimate the sample size required to detect aging-related brain variation. Three hundred and fifty subjects (208 females) passed the image quality control, with 51 subjects excluded due to excessive motion artifacts on 3T, 7T, or both. 7T MRI showed an overall stronger correlation between morphometrics and age and a larger number of significantly correlated brain volumes and cortical thicknesses. While the ICV is consistent between both field strengths, the Residual normalization method shows markedly higher correlation with age for 3T when compared with the Proportional normalization method. The 7T results are consistent regardless of the normalization method used. In a large cohort of healthy participants with paired 3T and 7T scans, we compared the statistical performance in assessing age-related brain morphological changes. Our study reaffirmed the inverse correlation between brain volumes and cortical thicknesses and age and highlighted varying correlations in different brain regions and normalization methods at 3T and 7T. 7T imaging significantly improves statistical power and thus reduces the required sample size.

Versatile MRI acquisition and processing protocol for population-based neuroimaging

Neuroimaging has an essential role in studies of brain health and of cerebrovascular and neurodegenerative diseases, requiring the availability of versatile magnetic resonance imaging (MRI) acquisition and processing protocols. We designed and developed a multipurpose high-resolution MRI protocol for large-scale and long-term population neuroimaging studies that includes structural, diffusion-weighted and functional MRI modalities. This modular protocol takes almost 1 h of scan time and is, apart from a concluding abdominal scan, entirely dedicated to the brain. The protocol links the acquisition of an extensive set of MRI contrasts directly to the corresponding fully automated data processing pipelines and to the required quality assurance of the MRI data and of the image-derived phenotypes. Since its successful implementation in the population-based Rhineland Study (ongoing, currently more than 11,000 participants, target participant number of 20,000), the proposed MRI protocol has proved suitable for epidemiological and clinical cross-sectional and longitudinal studies, including multisite studies. The approach requires expertise in magnetic resonance image acquisition, in computer science for the data management and the execution of processing pipelines, and in brain anatomy for the quality assessment of the MRI data. The protocol takes ~1 h of MRI acquisition and ~20 h of data processing to complete for a single dataset, but parallelization over multiple datasets using high-performance computing resources reduces the processing time. By making the protocol, MRI sequences and pipelines available, we aim to contribute to better comparability, interoperability and reusability of large-scale neuroimaging data.

Brain morphometrics correlations with age among 352 participants imaged with both 3T and 7T MRI: 7T improves statistical power and reduces required sample size

Introduction

Magnetic resonance imaging (MRI) at 7 Telsa (7T) has superior signal-to-noise ratio to 3 Telsa (3T) but also presents higher signal inhomogeneities and geometric distortions. A key knowledge gap is to robustly investigate the sensitivity and accuracy of 3T and 7T MRI in assessing brain morphometrics. This study aims to (a) aggregate a large number of paired 3T and 7T scans to evaluate their differences in quantitative brain morphological assessment using a widely available brain segmentation tool, FreeSurfer, as well as to (b) examine the impact of normalization methods for subject variability and smaller sample sizes on data analysis.

Methods

A total of 452 healthy participants aged 29 to 68 were imaged at both 3T and 7T. Structural T1-weighted magnetization-prepared rapid gradient-echo (MPRAGE) images were processed and segmented using FreeSurfer. To account for head size variability, the brain volumes underwent intracranial volume (ICV) correction using the Residual (regression model) and Proportional (simple division to ICV) methods. The resulting volumes and thicknesses were correlated with age using Pearson correlation and false discovery rate correction. The correlations were also calculated in increasing sample size from 3 to the whole sample to estimate the sample size required to detect aging-related brain variation.

Results

352 subjects (210 females) passed the image quality control with 100 subjects excluded due to excessive motion artifacts on 3T, 7T, or both. 7T MRI showed an overall stronger correlation between morphometrics and age and a larger number of significantly correlated brain volumes and cortical thicknesses. While the ICV is consistent between both field strengths, the Residual normalization method shows markedly higher correlation with age for 3T when compared with the Proportional normalization method. The 7T results are consistent regardless of the normalization method used.

Conclusion

In a large cohort of healthy participants with paired 3T and 7T scans, we compared the statistical performance in assessing age-related brain morphological changes. Our study reaffirmed the inverse correlation between brain volumes and cortical thicknesses and age and highlighted varying correlations in different brain regions and normalization methods at 3T and 7T. 7T imaging significantly improves statistical power and thus reduces required sample size.

Different MRI structural processing methods do not impact functional connectivity computation

Resting-state functional magnetic resonance imaging (rs-fMRI) has become an increasingly popular technique. This technique can assess several features of brain connectivity, such as inter-regional temporal correlation (functional connectivity), from which graph measures of network organization can be derived. However, these measures are prone to a certain degree of variability depending on the analytical steps during preprocessing. Many studies have investigated the effect of different preprocessing steps on functional connectivity measures; however, no study investigated whether different structural reconstructions lead to different functional connectivity metrics. Here, we evaluated the impact of different structural segmentation strategies on functional connectivity outcomes. To this aim, we compared different metrics computed after two different registration strategies. The first strategy used structural information from the 3D T1-weighted image (unimodal), while the second strategy implemented a multimodal approach, where an additional registration step used the information from the T2-weighted image. The impact of these different approaches was evaluated on a sample of 58 healthy adults. As expected, different approaches led to significant differences in structural measures (i.e., cortical thickness, volume, and gyrification index), with the maximum impact on the insula cortex. However, these differences were only slightly translated to functional metrics. We reported no differences in graph measures and seed-based functional connectivity maps, but slight differences in the insula when we compared the mean functional strength for each parcel. Overall, these results suggested that functional metrics are only slightly different when using a unimodal compared to a multimodal approach, while the structural output can be significantly affected.

Big Field of View MRI T1w and FLAIR Template - NMRI225

Image templates are a common tool for neuroscience research. Often, they are used for spatial normalization of magnetic resonance imaging (MRI) data, which is a necessary procedure for analyzing brain morphology and function via voxel-based analysis. This allows the researcher to reduce individual shape differences across images and make inferences across multiple subjects. Many templates have a small field-of-view typically focussed on the brain, limiting the use for applications requiring detailed information about other extra-cranial structures in the head and neck area. However, there are several applications where such information is important, for example source reconstruction of electroencephalography (EEG) and/or magnetoencephalography (MEG). We have constructed a new template based on 225 T1w and FLAIR images with a big field-of-view that can serve both as target for across subject spatial normalization as well as a basis to build high-resolution head models. This template is based on and iteratively re-registered to the MNI152 space to provide maximal compatibility with the most commonly used brain MRI template.

The Unresolved Methodological Challenge of Detecting Neuroplastic Changes in Astronauts

After completing a spaceflight, astronauts display a salient upward shift in the position of the brain within the skull, accompanied by a redistribution of cerebrospinal fluid. Magnetic resonance imaging studies have also reported local changes in brain volume following a spaceflight, which have been cautiously interpreted as a neuroplastic response to spaceflight. Here, we provide evidence that the grey matter volume changes seen in astronauts following spaceflight are contaminated by preprocessing errors exacerbated by the upwards shift of the brain within the skull. While it is expected that an astronaut's brain undergoes some neuroplastic adaptations during spaceflight, our findings suggest that the brain volume changes detected using standard processing pipelines for neuroimaging analyses could be contaminated by errors in identifying different tissue types (i.e., tissue segmentation). These errors may undermine the interpretation of such analyses as direct evidence of neuroplastic adaptation, and novel or alternate preprocessing or experimental paradigms are needed in order to resolve this important issue in space health research.

Analysis of synthetic magnetic resonance images by multi-channel segmentation increases accuracy of volumetry in the putamen and decreases mis-segmentation in the dural sinuses

BACKGROUND

Voxel-based morphometry (VBM) using magnetic resonance imaging (MR) has been used to estimate cortical atrophy associated with various diseases. However, there are mis-segmentations of segmented gray matter image in VBM.

PURPOSE

To study a twofold evaluation of single- and multi-channel segmentation using synthetic MR images: (1) mis-segmentation of segmented gray matter images in transverse and cavernous sinuses; and (2) accuracy and repeatability of segmented gray matter images.

MATERIAL AND METHODS

A total of 13 healthy individuals were scanned with 3D quantification using an interleaved Look-Locker acquisition sequence with a T2 preparation pulse (3D-QALAS) sequence on a 1.5-T scanner. Three of the 13 healthy participants were scanned five consecutive times for evaluation of repeatability. We used SyMRI software to create images with three contrasts: T1-weighted (T1W), T2-weighted (T2W), and proton density-weighted (PDW) images. Manual regions of interest (ROI) on T1W imaging were individually set as the gold standard in the transverse sinus, cavernous sinus, and putamen. Single-channel (T1W) and multi-channel (T1W + T2W, T1W + PDW, and T1W + T2W + PDW imaging) segmentations were performed with statistical parametric mapping 12 software.

RESULTS

We found that mis-segmentations in both the transverse and cavernous sinuses were large in single-channel segmentation compared with multi-channel segmentations. Furthermore, the accuracy of segmented gray matter images in the putamen was high in both multi-channel T1W + PDW and T1W + T2W + PDW segmentations compared with other segmentations. Finally, the highest repeatability of left putamen volumetry was found with multi-channel segmentation T1WI + PDWI.

CONCLUSION

Multi-channel segmentation with T1WI + PDWI provides good results for VBM compared with single-channel and other multi-channel segmentations.

Increased Resolution of Structural MRI at 3T Improves Estimation of Regional Cortical Degeneration in Individual Dementia Patients Using Surface Thickness Maps

BACKGROUND

Objective measurement of regional cortical atrophy in individual patients would be a highly desirable adjunct for diagnosis of degenerative dementias.

OBJECTIVE

We hypothesized that increasing the resolution of magnetic resonance scans would improve the sensitivity of cortical atrophy detection for individual patients.

METHODS

46 participants including 8 semantic-variant primary progressive aphasia (svPPA), seven posterior cortical atrophy (PCA), and 31 cognitively unimpaired participants underwent clinical assessment and 3.0T brain scans. SvPPA and PCA were chosen because there is overwhelming prior knowledge of the expected atrophy pattern. Two sets of T1-weighted images with 0.8 mm3 (HighRes) and conventional 1.0 mm3 (ConvRes) resolution were acquired. The cortical ribbon was segmented using FreeSurfer software to obtain surface-based thickness maps. Inter-sequence performance was assessed in terms of cortical thickness and sub-cortical volume reproducibility, signal-to-noise and contrast-to-noise ratios. For clinical cases, diagnostic effect size (Cohen's d) and lesion distribution (z-score and t-value maps) were compared between HighRes and ConvRes scans.

RESULTS

The HighRes scans produced higher image quality scores at 90 seconds extra scan time. The effect size of cortical thickness differences between patients and cognitively unimpaired participants was 15-20% larger for HighRes scans. HighRes scans showed more robust patterns of atrophy in expected regions in each and every individual patient.

CONCLUSIONS

HighRes T1-weighted scans showed superior precision for identifying the severity of cortical atrophy in individual patients, offering a proof-of-concept for clinical translation. Studying svPPA and PCA, two syndromes with well-defined focal atrophy patterns, offers a method to clinically validate and contrast automated algorithms.

MIL normalization -- prerequisites for accurate MRI radiomics analysis

The quality of magnetic resonance (MR) images obtained with different instruments and imaging parameters varies greatly. A large number of heterogeneous images are collected, and they suffer from acquisition variation. Such imaging quality differences will have a great impact on the radiomics analysis. The main differences in MR images include modality mismatch (M), intensity distribution variance (I), and layer-spacing differences (L), which are referred to as MIL differences in this paper for convenience. An MIL normalization system is proposed to reconstruct uneven MR images into high-quality data with complete modality, a uniform intensity distribution and consistent layer spacing. Three radiomics tasks, including tumor segmentation, pathological grading and genetic diagnosis of glioma, were used to verify the effect of MIL normalization on radiomics analysis. Three retrospective glioma datasets were analyzed in this study: BraTs (285 cases), TCGA (112 cases) and HuaShan (403 cases). They were used to test the effect of MIL on three different radiomics tasks, including tumor segmentation, pathological grading and genetic diagnosis. MIL normalization included three components: multimodal synthesis based on an encoder-decoder network, intensity normalization based on CycleGAN, and layer-spacing unification based on Statistical Parametric Mapping (SPM). The Dice similarity coefficient, areas under the curve (AUC) and six other indicators were calculated and compared after different normalization steps. The MIL normalization system can improved the Dice coefficient of segmentation by 9% (P < .001), the AUC of pathological grading by 32% (P < .001), and IDH1 status prediction by 25% (P < .001) when compared to non-normalization. The proposed MIL normalization system provides high-quality standardized data, which is a prerequisite for accurate radiomics analysis.

Automatic multispectral MRI segmentation of human hippocampal subfields: an evaluation of multicentric test-retest reproducibility

Accurate and reproducible automated segmentation of human hippocampal subfields is of interest to study their roles in cognitive functions and disease processes. Multispectral structural MRI methods have been proposed to improve automated hippocampal subfield segmentation accuracy, but the reproducibility in a multicentric setting is, to date, not well characterized. Here, we assessed test-retest reproducibility of FreeSurfer 6.0 hippocampal subfield segmentations using multispectral MRI analysis pipelines (22 healthy subjects scanned twice, a week apart, at four 3T MRI sites). The harmonized MRI protocol included two 3D-T1, a 3D-FLAIR, and a high-resolution 2D-T2. After within-session T1 averaging, subfield volumes were segmented using three pipelines with different multispectral data: two longitudinal ("long_T1s" and "long_T1s_FLAIR") and one cross-sectional ("long_T1s_FLAIR_crossT2"). Volume reproducibility was quantified in magnitude (reproducibility error-RE) and space (DICE coefficient). RE was lower in all hippocampal subfields, except for hippocampal fissure, using the longitudinal pipelines compared to long_T1s_FLAIR_crossT2 (average RE reduction of 0.4-3.6%). Similarly, the longitudinal pipelines showed a higher spatial reproducibility (1.1-7.8% of DICE improvement) in all hippocampal structures compared to long_T1s_FLAIR_crossT2. Moreover, long_T1s_FLAIR provided a small but significant RE improvement in comparison to long_T1s (p = 0.015), whereas no significant DICE differences were found. In addition, structures with volumes larger than 200 mm3 had better RE (1-2%) and DICE (0.7-0.95) than smaller structures. In summary, our study suggests that the most reproducible hippocampal subfield FreeSurfer segmentations are derived from a longitudinal pipeline using 3D-T1s and 3D-FLAIR. Adapting a longitudinal pipeline to include high-resolution 2D-T2 may lead to further improvements.

Distinct Interplay Between Atrophy and Hypometabolism in Alzheimer's Versus Semantic Dementia

Multimodal neuroimaging analyses offer additional information beyond that provided by each neuroimaging modality. Thus, direct comparisons and correlations between neuroimaging modalities allow revealing disease-specific topographic relationships. Here, we compared the topographic discrepancies between atrophy and hypometabolism in two neurodegenerative diseases characterized by distinct pathological processes, namely Alzheimer's disease (AD) versus semantic dementia (SD), to unravel their specific influence on local and global brain structure-function relationships. We found that intermodality topographic discrepancies clearly distinguished the two patient groups: AD showed marked discrepancies between both alterations, with greater hypometabolism than atrophy in large posterior associative neocortical regions, while SD showed more topographic consistency between atrophy and hypometabolism across brain regions. These findings likely reflect the multiple pathologies versus the relatively unitary pathological process underlying AD versus SD respectively. Our results evidence that multimodal neuroimaging-derived indexes can provide clinically relevant information to discriminate the two diseases, and potentially reveal distinct neuropathological processes.

Amygdalar nuclei and hippocampal subfields on MRI: Test-retest reliability of automated volumetry across different MRI sites and vendors

BACKGROUND

The amygdala and the hippocampus are two limbic structures that play a critical role in cognition and behavior, however their manual segmentation and that of their smaller nuclei/subfields in multicenter datasets is time consuming and difficult due to the low contrast of standard MRI. Here, we assessed the reliability of the automated segmentation of amygdalar nuclei and hippocampal subfields across sites and vendors using FreeSurfer in two independent cohorts of older and younger healthy adults.

METHODS

Sixty-five healthy older (cohort 1) and 68 younger subjects (cohort 2), from the PharmaCog and CoRR consortia, underwent repeated 3D-T1 MRI (interval 1-90 days). Segmentation was performed using FreeSurfer v6.0. Reliability was assessed using volume reproducibility error (ε) and spatial overlapping coefficient (DICE) between test and retest session.

RESULTS

Significant MRI site and vendor effects (p ​< ​.05) were found in a few subfields/nuclei for the ε, while extensive effects were found for the DICE score of most subfields/nuclei. Reliability was strongly influenced by volume, as ε correlated negatively and DICE correlated positively with volume size of structures (absolute value of Spearman's r correlations >0.43, p ​< ​1.39E-36). In particular, volumes larger than 200 ​mm³ (for amygdalar nuclei) and 300 ​mm³ (for hippocampal subfields, except for molecular layer) had the best test-retest reproducibility (ε ​< ​5% and DICE ​> ​0.80).

CONCLUSION

Our results support the use of volumetric measures of larger amygdalar nuclei and hippocampal subfields in multisite MRI studies. These measures could be useful for disease tracking and assessment of efficacy in drug trials.

Identifying errors in Freesurfer automated skull stripping and the incremental utility of manual intervention

Quality assurance (QA) is vital for ensuring the integrity of processed neuroimaging data for use in clinical neurosciences research. Manual QA (visual inspection) of processed brains for cortical surface reconstruction errors is resource-intensive, particularly with large datasets. Several semi-automated QA tools use quantitative detection of subjects for editing based on outlier brain regions. There were two project goals: (1) evaluate the assumption that statistical outliers are related to errors of cortical extension, and (2) examine whether error identification and correction significantly impacts estimation of cortical parameters and established brain-behavior relationships. T1 MPRAGE images (N = 530) of healthy adults were obtained from the NKI-Rockland Sample and reconstructed using Freesurfer 5.3. Visual inspection of T1 images was conducted for: (1) participants (n = 110) with outlier values (z scores ±3 SD) for subcortical and cortical segmentation volumes (outlier group), and (2) a random sample of remaining participants (n = 110) with segmentation values that did not meet the outlier criterion (non-outlier group). The outlier group had 21% more participants with visual inspection-identified errors than participants in the non-outlier group, with a medium effect size (Φ = 0.22). Nevertheless, a considerable portion of images with errors of cortical extension were found in the non-outlier group (41%). Although nine brain regions significantly changed size from pre- to post-editing (with effect sizes ranging from 0.26 to 0.59), editing did not substantially change the correlations of neurocognitive tasks and brain volumes (ps > 0.05). Statistically-based QA, although less resource intensive, is not accurate enough to supplant visual inspection. We discuss practical implications of our findings to guide resource allocation decisions for image processing.

A Method to Estimate Brain Volume from Head CT Images and Application to Detect Brain Atrophy in Alzheimer Disease

BACKGROUND AND PURPOSE

Total brain volume and total intracranial volume are important measures for assessing whole-brain atrophy in Alzheimer disease, dementia, and other neurodegenerative diseases. Unlike MR imaging, which has a number of well-validated fully-automated methods, only a handful of methods segment CT images. Available methods either use enhanced CT, do not estimate both volumes, or require formal validation. Reliable computation of total brain volume and total intracranial volume from CT is needed because head CTs are more widely used than head MRIs in the clinical setting. We present an automated head CT segmentation method (CTseg) to estimate total brain volume and total intracranial volume.

MATERIALS AND METHODS

CTseg adapts a widely used brain MR imaging segmentation method from the Statistical Parametric Mapping toolbox using a CT-based template for initial registration. CTseg was tested and validated using head CT images from a clinical archive.

RESULTS

CTseg showed excellent agreement with 20 manually segmented head CTs. The intraclass correlation was 0.97 (P < .001) for total intracranial volume and 0.94 (P < .001) for total brain volume. When CTseg was applied to a cross-sectional Alzheimer disease dataset (58 with Alzheimer disease patients and 58 matched controls), CTseg detected a loss in percentage total brain volume (as a percentage of total intracranial volume) with age (P < .001) as well as a group difference between patients with Alzheimer disease and controls (P < .01). We observed similar results when total brain volume was modeled with total intracranial volume as a confounding variable.

CONCLUSIONS

In current clinical practice, brain atrophy is assessed by inaccurate and subjective "eyeballing" of CT images. Manual segmentation of head CT images is prohibitively arduous and time-consuming. CTseg can potentially help clinicians to automatically measure total brain volume and detect and track atrophy in neurodegenerative diseases. In addition, CTseg can be applied to large clinical archives for a variety of research studies.

Test-retest variability of brain morphometry analysis: an investigation of sequence and coil effects

Background

Precise and reliable brain morphometry analysis is critical for clinical and research purposes. The magnetization-prepared rapid gradient echo (MPRAGE), multi-echo MPRAGE (MEMPRAGE) and magnetization-prepared 2 rapid acquisition gradient echo (MP2RAGE) sequences have all been used to acquire brain structural images, but it is unclear which of these sequences is the most suitable for brain morphometry and whether the number of coil channels (20 or 32) affects scan precision. This study aimed to assess the impact of T1-weighted image acquisition variables (sequence and head coil) on the repeatability of resultant automated volumetric measurements.

Methods

Twenty-four healthy volunteers underwent back-to-back scanning protocols with three sequences and two different coils (i.e., six scanning conditions in total) presented in a randomized order in a single session. MorphoBox prototype and FreeSurfer were used for brain segmentation. Brain structures were divided into cortical and subcortical regions for more precise analysis. The acquired volume and thickness values were used to calculate test-retest variability (TRV) values. TRV values from the six different combinations were compared for total brain structures, total cortical structures, total subcortical structures, and every single structure.

Results

The median TRV value for all brain regions was 1.23% with MorphoBox and 3.14% with FreeSurfer. When using FreeSurfer results to compare the six combinations, for total brain structures volume and total cortical structures volume and thickness, the MEMPRAGE-32 channel combination showed significantly lower TRV values than the others (P<0.01). Similar results were observed with MorphoBox. For total subcortical structures, the MP2RAGE-32 channel combination showed the lowest TRV values with both MorphoBox (lower about 0.01% to 0.17%) and FreeSurfer analyses (lower about 0.02% to 0.37%).

Conclusions

TRV values were generally low, indicating generally high reliability for every region. The MEMPRAGE sequence was the most reliable of the three sequences for total brain structures and cortical structures. However, MP2RAGE was the most reliable for subcortical structures. The 32-channel coil showed better repeatability results than the 20-channel coil.

MP2RAGE multispectral voxel-based morphometry in focal epilepsy

We assessed the applicability of MP2RAGE for voxel‐based morphometry. To this end, we analyzed its brain tissue segmentation characteristics in healthy subjects and the potential for detecting focal epileptogenic lesions (previously visible and nonvisible). Automated results and expert visual interpretations were compared with conventional VBM variants (i.e., T1 and T1 + FLAIR). Thirty‐one healthy controls and 21 patients with focal epilepsy were recruited. 3D T1‐, T2‐FLAIR, and MP2RAGE images (consisting of INV1, INV2, and MP2 maps) were acquired on a 3T MRI. The effects of brain tissue segmentation and lesion detection rates were analyzed among single‐ and multispectral VBM variants. MP2‐single‐contrast gave better delineation of deep, subcortical nuclei but was prone to misclassification of dura/vessels as gray matter, even more than conventional‐T1. The addition of multispectral combinations (INV1, INV2, or FLAIR) could markedly reduce such misclassifications. MP2 + INV1 yielded generally clearer gray matter segmentation allowing better differentiation of white matter and neighboring gyri. Different models detected known lesions with a sensitivity between 60 and 100%. In non lesional cases, MP2 + INV1 was found to be best with a concordant rate of 37.5%, specificity of 51.6% and concordant to discordant ratio of 0.60. In summary, we show that multispectral MP2RAGE VBM (e.g., MP2 + INV1, MP2 + INV2) can improve brain tissue segmentation and lesion detection in epilepsy.

Examining the identification of age-related atrophy between T1 and T1 + T2-FLAIR cortical thickness measurements

Cortical thickness is traditionally derived from T1-weighted MRI images. Recent studies have shown an improvement in segmentation with the combination of T1 + T2-FLAIR images. MRI data from 54 adults (mean: 71 years, 65–81 years, 48% females) that are part of an ongoing cohort study were analyzed to investigate whether T1 + T2-FLAIR cortical thickness measurements were superior to those derived from T1-weighted images in identifying age-related atrophy. T1-weighted and T2-FLAIR MRI images were processed through FreeSurfer v6.0. Data was extracted using the Desikan-Killiany (DKT) atlas. FreeSurfer’s GUI QDEC examined age-related atrophy. Nonparametric tests, effect sizes, and Pearson correlations examined differences between T1-only and T1 + T2-FLAIR cortical thickness data. These analyses demonstrated that T1 + T2-FLAIR processed images significantly improved the segmentation of gray matter (chi-square x², p < 0.05) and demonstrated significantly thicker cortical thickness means (p < 0.05) with medium to large effect sizes. Significant regions of age-related cortical atrophy were identified within the T1 + T2-FLAIR data (FDR corrected, p < 0.05). This is in contrast to the T1-only data where no regions survived FDR correction. In summary, T1 + T2-FLAIR data were associated with significant improvement in cortical segmentation and the identification of age-related cortical atrophy. Future studies should consider employing this imaging strategy to obtain cortical thickness measurements sensitive to age-related changes.

Brain tissue segmentation based on MP2RAGE multi-contrast images in 7 T MRI

We proposed a method for segmentation of brain tissues-gray matter, white matter, and cerebrospinal fluid-using multi-contrast images, including a T1 map and a uniform T1-weighted image, from a magnetization-prepared 2 rapid acquisition gradient echoes (MP2RAGE) sequence at 7 Tesla. The proposed method was evaluated with respect to the processing time and the similarity of the segmented masks of brain tissues with those obtained using FSL, FreeSurfer, and SPM12. The processing time of the proposed method (28 ± 0 s) was significantly shorter than those of FSL and SPM12 (444 ± 4 s and 159 ± 2 s for FSL and SPM12, respectively). In the similarity assessment, the tissue mask of the brain obtained by the proposed method showed higher consistency with those obtained using FSL than with those obtained using SPM12. The proposed method misclassified the subcortical structures and large vessels since it is based on the intensities of multi-contrast images obtained using MP2RAGE, which uses a similar segmentation approach as FSL but is not based on a template image or a parcellated brain atlas, which are used for FreeSurfer and SPM12, respectively. However, the proposed method showed good segmentation in the cerebellum and white matter in the medial part of the brain in comparison with the other methods. Thus, because the proposed method using different contrast images of MP2RAGE sequence showed the shortest processing time and similar segmentation ability as the other methods, it may be useful for both neuroimaging research and clinical diagnosis.

Systematic Assessment of Multispectral Voxel-Based Morphometry in Previously MRI-Negative Focal Epilepsy

BACKGROUND AND PURPOSE

Voxel-based morphometry is widely used for detecting gray matter abnormalities in epilepsy. However, its performance with changing parameters, smoothing and statistical threshold, is debatable. More important, the potential yield of combining multiple MR imaging contrasts (multispectral voxel-based morphometry) is still unclear. Our aim was to objectify smoothing and statistical cutoffs and systematically compare the performance of multispectral voxel-based morphometry with existing T1 voxel-based morphometry in patients with focal epilepsy and previously negative MRI.

MATERIALS AND METHODS

3D T1-, T2-, and T2-weighted FLAIR scans were acquired for 62 healthy volunteers and 13 patients with MR imaging negative for focal epilepsy on a Magnetom Skyra 3T scanner with an isotropic resolution of 0.9 mm³. We systematically optimized the main voxel-based morphometry parameters, smoothing level and statistical cutoff, with T1 voxel-based morphometry as a reference. As a next step, the performance of multispectral voxel-based morphometry models, T1+T2, T1+FLAIR, and T1+T2+FLAIR, was compared with that of T1 voxel-based morphometry using gray matter concentration and gray matter volume analysis.

RESULTS

We found the best performance of T1 at 12 mm and a T-threshold (statistical cutoff) of 3.7 for gray matter concentration analysis. When we incorporated these parameters, after expert visual interpretation of concordant and discordant findings, we identified T1+FLAIR as the best model with a concordant rate of 46.2% and a concordant rate/discordant rate of 1.20 compared with T1 with 30.8% and 0.67, respectively. Visual interpretation of voxel-based morphometry findings decreased concordant rates from 38.5%-46.2% to 15.4%-46.2% and discordant rates from 53.8%-84.6% to 30.8%-46.2% and increased specificity across models from 33.9%-40.3% to 46.8%-54.8%.

CONCLUSIONS

Multispectral voxel-based morphometry, especially T1+FLAIR, can yield superior results over single-channel T1 in focal epilepsy patients with a negative conventional MR imaging.

The morphology of the human cerebrovascular system

While several methodologies exist for quantifying gray and white matter properties in humans, relatively little is known regarding the spatial organization and the intersubject variability of cerebral vessels. To resolve this, we developed a fast, open-source processing algorithm using advanced vessel segmentation schemes and iterative nonlinear registration to isolate, extract, and quantify cerebral vessels in susceptibility weighting imaging (SWI) and time-of-flight angiography (TOF-MRA) datasets acquired in a large cohort (n = 42) of healthy individuals. From this, whole-brain venous and arterial probabilistic maps were generated along with the computation of regional densities and diameters within regions based on popular anatomical and functional atlases. The results show that cerebral vasculature is highly heterogeneous, displaying disproportionally large vessel densities in brain areas such as the anterior and posterior cingulate, cuneus, precuneus, parahippocampus, insula, and temporal gyri. On average, venous densities were slightly higher and less variable across subjects than arterial. Moreover, regional variations in both venous and arterial density were significantly correlated to cortical thickness (R = 0.42). This publicly available new atlas of the human cerebrovascular system provides a first step toward quantifying morphological changes in the diseased brain and serving as a potential regression tool in fMRI analysis.

A scalable method to improve gray matter segmentation at ultra high field MRI

High-resolution (functional) magnetic resonance imaging (MRI) at ultra high magnetic fields (7 Tesla and above) enables researchers to study how anatomical and functional properties change within the cortical ribbon, along surfaces and across cortical depths. These studies require an accurate delineation of the gray matter ribbon, which often suffers from inclusion of blood vessels, dura mater and other non-brain tissue. Residual segmentation errors are commonly corrected by browsing the data slice-by-slice and manually changing labels. This task becomes increasingly laborious and prone to error at higher resolutions since both work and error scale with the number of voxels. Here we show that many mislabeled, non-brain voxels can be corrected more efficiently and semi-automatically by representing three-dimensional anatomical images using two-dimensional histograms. We propose both a uni-modal (based on first spatial derivative) and multi-modal (based on compositional data analysis) approach to this representation and quantify the benefits in 7 Tesla MRI data of nine volunteers. We present an openly accessible Python implementation of these approaches and demonstrate that editing cortical segmentations using two-dimensional histogram representations as an additional post-processing step aids existing algorithms and yields improved gray matter borders. By making our data and corresponding expert (ground truth) segmentations openly available, we facilitate future efforts to develop and test segmentation algorithms on this challenging type of data.

The impact of B1+ correction on MP2RAGE cortical T1 and apparent cortical thickness at 7T

Determination of cortical thickness using MRI has often been criticized due to the presence of various error sources. Specifically, anatomical MRI relying on T₁ contrast may be unreliable due to spatially variable image contrast between gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF). Especially at ultra‐high field (≥ 7T) MRI, transmit and receive B₁‐related image inhomogeneities can hamper correct classification of tissue types. In the current paper, we demonstrate that residual B1+ (transmit) inhomogeneities in the T₁‐weighted and quantitative T₁ images using the MP2RAGE sequence at 7T lead to biases in cortical thickness measurements. As expected, post‐hoc correction for the spatially varying B1+ profile reduced the apparent T₁ values across the cortex in regions with low B1+, and slightly increased apparent T₁ in regions with high B1+. As a result, improved contrast‐to‐noise ratio both at the GM‐CSF and GM‐WM boundaries can be observed leading to more accurate surface reconstructions and cortical thickness estimates. Overall, the changes in cortical thickness ranged between a 5% decrease to a 70% increase after B1+ correction, reducing the variance of cortical thickness values across the brain dramatically and increasing the comparability with normative data. More specifically, the cortical thickness estimates increased in regions characterized by a strong decrease of apparent T₁ after B1+ correction in regions with low B1+ due to improved detection of the pial surface_. The current results suggest that cortical thickness can be more accurately determined using MP2RAGE data at 7T if B1+ inhomogeneities are accounted for.