Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Application of Quantitative MRI for Brain Tissue Segmentation at 1.5 T and 3.0 T Field Strengths

  • Janne West ,

    janne.west@liu.se

    Affiliations Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden, Radiation Physics, Department of Medicine and Health, Linköping University, UHL County Council of Östergötland, Linköping, Sweden

  • Ida Blystad,

    Affiliations Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden, Radiology, Department of Medicine and Health, Linköping University, UHL County Council of Östergötland, Linköping, Sweden

  • Maria Engström,

    Affiliations Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden, Radiology, Department of Medicine and Health, Linköping University, UHL County Council of Östergötland, Linköping, Sweden

  • Jan B. M. Warntjes,

    Affiliations Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden, Division of Clinical Physiology, Department of Medicine and Health, Linköping University, Linköping, Sweden

  • Peter Lundberg

    Affiliations Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden, Radiation Physics, Department of Medicine and Health, Linköping University, UHL County Council of Östergötland, Linköping, Sweden, Radiology, Department of Medicine and Health, Linköping University, UHL County Council of Östergötland, Linköping, Sweden

Abstract

Background

Brain tissue segmentation of white matter (WM), grey matter (GM), and cerebrospinal fluid (CSF) are important in neuroradiological applications. Quantitative Mri (qMRI) allows segmentation based on physical tissue properties, and the dependencies on MR scanner settings are removed. Brain tissue groups into clusters in the three dimensional space formed by the qMRI parameters R1, R2 and PD, and partial volume voxels are intermediate in this space. The qMRI parameters, however, depend on the main magnetic field strength. Therefore, longitudinal studies can be seriously limited by system upgrades. The aim of this work was to apply one recently described brain tissue segmentation method, based on qMRI, at both 1.5 T and 3.0 T field strengths, and to investigate similarities and differences.

Methods

In vivo qMRI measurements were performed on 10 healthy subjects using both 1.5 T and 3.0 T MR scanners. The brain tissue segmentation method was applied for both 1.5 T and 3.0 T and volumes of WM, GM, CSF and brain parenchymal fraction (BPF) were calculated on both field strengths. Repeatability was calculated for each scanner and a General Linear Model was used to examine the effect of field strength. Voxel-wise t-tests were also performed to evaluate regional differences.

Results

Statistically significant differences were found between 1.5 T and 3.0 T for WM, GM, CSF and BPF (p<0.001). Analyses of main effects showed that WM was underestimated, while GM and CSF were overestimated on 1.5 T compared to 3.0 T. The mean differences between 1.5 T and 3.0 T were -66 mL WM, 40 mL GM, 29 mL CSF and -1.99% BPF. Voxel-wise t-tests revealed regional differences of WM and GM in deep brain structures, cerebellum and brain stem.

Conclusions

Most of the brain was identically classified at the two field strengths, although some regional differences were observed.

Introduction

Brain tissue segmentation and volume estimation of white matter (WM), grey matter (GM), and cerebrospinal fluid (CSF) are important in many neuroradiological applications [1,2,3]. Volume estimation of segmented brain tissue can be used to assess local and global brain atrophy that is present in diseases such as multiple sclerosis [1,4], Alzheimer’s disease [5] and other forms of dementia [6], and this measure can in turn be used to monitor neurodegenerative disease progression [5,7], as well as to monitor the effects of therapy and rehabilitation [8].

There are two general approaches to automated brain tissue segmentation; either (1) each voxel is assigned to one specific tissue class [9,10,11,12], or (2) each voxel is assigned volume fractions of several different tissue classes [13,14,15]. Allowing more than one tissue class in each voxel has advantages since partial volume is common in the brain, even at high-resolution imaging, and especially close to tissue borders, such as in the cortex where voxels often contain mixtures of GM and WM, or CSF. Allowing only one tissue class in each voxel is prone to partial-volume errors. Furthermore, calculating tissue fractions decreases the dependence on high-resolution imaging, something that requires a large signal-to-noise-ratio (SNR).

Most brain tissue segmentation methods are based on conventional contrast-weighted MR images [10,11,12,13,14,15]. Although these images provide high neuro-anatomical detail, segmentation is complicated by contrast inhomogeneities and the arbitrary grey-scale of the images. Furthermore a variety of experimental parameters influence conventional contrast-weighted MR images, such as the repetition time (TR), echo time (TE), and inversion delay time (TI). The complexity of the acquired images problematises physical interpretation of the signal intensity. Therefore, segmentation is often performed using contrast differences in the image, and not the absolute pixel values. Different normalisation strategies or ad hoc filters may be necessary to compensate for variations of signal intensity in the imaged volume (see 10 or [15]).

Conversely, quantitative MRI (qMRI) allows the measurements of physical tissue parameters such as the relaxation rates R1 (inverse of the relaxation time T1, 1/T1) and R2 (1/T2), and the proton density (PD) [16,17,18,19,20,21,22,23]. By using qMRI for tissue segmentation, the dependence on MR scanner hardware and settings is largely removed. This allows for segmentation based on physical tissue characteristics [24]. In recent years quantification methods that are feasible within clinically acceptable scanning times have been presented [18,22]. It has also been shown that human brain tissue can be differentiated using qMRI parameters [9,25] and that tissues group into clusters in the three-dimensional space formed by R1, R2 and PD (the R1-R2-PD space) [22,24]. In those voxels where two tissues are present, both tissue types contribute to the MR-signal, and quantification results in a weighted average of the tissues’ individual R1, R2 and PD values. The qMRI parameters describing the relaxation rates of tissue, however, depend on the main magnetic field strength. The R1 and R2 parameters generally decrease when the main magnetic field strength is increased; this effect is not linear and depends on tissue composition. In particular, signal alteration in tissues rich in iron increases due to the ferromagnetic properties of the iron, causing decreased R1 and R2 [26,27,28].

Because of this, longitudinal studies, as well as quantitative patient follow-ups, can be seriously limited by system changes from lower to higher field strengths. Many widely used segmentation tools, such as FreeSurfer [29] and Bayesian based methods [30], report variable results across field strengths in combination with conventional MRI [27,31], and combining data acquired at different field strengths has for this reason not been endorsed. Studies have reported regional differences across field strengths when using automatic segmentation methods [27,31,32], as well as overall volume differences of brain tissues and intracranial volume [33]. One study by Jovicich et al reported differences of up to 10% for total intracranial volume when using the FreeFurfer tool to compare 1.5 T and 3.0 T [31]. In another study Keihaninejad et al found differences of up to 12.5% using Bayesian methods and up to 14.7% using SPM, also comparing intracranial volume between 1.5 T and 3.0 T systems. Moreover, one study even reported systematic volume differences when manual segmentation was used [33].

The aim of this work was to apply one recently described brain tissue segmentation method based on qMRI [24] at both 1.5 T and 3.0 T field strengths, and to investigate the differences and potential value of combining field strengths in research studies and clinical practice.

Materials and Methods

Overview

In this work one brain tissue segmentation method, previously only reported for 1.5 T, was implemented for both 1.5 T and 3.0 T scanner systems [24]. Data acquisition was performed on 10 healthy subjects at both field strength, and results were compared. The segmentation method required prior knowledge about (1) ‘pure tissue clusters’ in the three dimensional space formed by using R1, R2 and PD as basis (the R1-R2-PD space) and (2) the characteristics of partial volume voxels in this space. Prior to performing the brain tissue segmentation data from the 10 healthy subjects were used to obtain this prior knowledge, which was then used to calibrate the segmentation method. Pure tissue clusters were determined by region of interest (ROI) measurements, and the characteristics of partial volume voxels were estimated with simulations. In the simulations the measured pure tissue cluster values were combined to create synthetic tissue mixtures. Finally; data from the 10 healthy subjects were segmented and results from the two field strengths were compared. The study design is illustrated in Figure 1.

thumbnail
Figure 1. Flow chart shows the different steps performed in this study.

First (1) data acquisition was performed on 10 healthy subjects at both 1.5 T and 3.0 T field strengths (2), prior knowledge required in the segmentation procedure were determined (3), the brain tissue segmentation procedure were performed, and finally (4) the results were compared between the field strengths.

https://doi.org/10.1371/journal.pone.0074795.g001

Data Acquisition

In vivo measurements were performed with a 1.5 T Philips Achieva scanner (software version 3.2.1), using an 8-channel SENSE head coil, and a 3.0 T Philips Ingenia scanner (software version 4.1.1), using an 16-channel SENSE head coil, on 10 healthy subjects (four males, six females, aged 24.4 ± 2.5, range 21 to 29 years). Health status was assessed using a questionnaire where brain disorders were ruled out. QMRI data were acquired using the multi-slice, multi-echo and multi-saturation delay sequence ‘QRAPMASTER’ (also known as qMAP) [22]. Whole-brain qMRI maps of the tissue parameters R1, R2, and PD were calculated using SyMRI Diagnostics (0.9.11) software (SyntheticMR AB, Sweden, 2012) from the acquired raw data. The qMRI sequence settings were; repetition time (TR) = 5000 ms, saturation delay times (TD) = 132/632/2257/4882 ms, echo times (TE) = 16/32/48/64/80/96 ms. Voxel size was 3x1.3x1.3 mm3 and a gap of 0.5 mm was applied in order to reach whole-brain coverage. The sequence was identical for 1.5 T and 3.0 T, and each scan took 6 minutes. For additional details on the QRAPMASTER qMRI sequence see 22. In addition, conventional T2-weighted images were acquired. For each of the two scanners, the subjects were scanned in two separate sessions, and were removed from the scanner room between these sessions. In each session, two acquisitions were performed. This resulted in a total of eight data sets for each subject (two scanners * two sessions * two acquisitions). The study was approved by the Regional Ethical Review Board in Linköping (Dnr: M88-07 T93-08) and written informed consent was obtained from all subjects prior to study entry.

Prior Knowledge - Pure Brain Tissue Clusters

In order to measure the pure brain tissue clusters in R1-R2-PD space, required as prior knowledge in the segmentation procedure, a neuroradiologist (IB) manually positioned standardised regions of interest (ROIs). The ROI size was 3x3 voxels (corresponding to 3.9x3.9 mm2) and the ROIs were placed on synthetic T2-weighted images with conventional T2-weighted images used as a visual reference. Synthetic T2-weighted images were generated from the quantitative R1-R2-PD data using the SyMRI Diagnostics (0.9.11) software (SyntheticMR AB, Sweden, 2012) as described elsewhere [22]. Synthetic T2-weighted images were generated from the same data as the qMRI parameter maps and therefore perfectly registered, removing the need for additional data registration. Both 1.5 T and 3.0 T image ROIs were placed in the same session in order to assure equivalent ROI placement on images from both field strengths. ROIs were placed on both datasets from the first session of each subject in each scanner, resulting in a total of four datasets for each individual subject (two for each of the field strengths). ROIs were placed in the structures of GM, WM and CSF. GM R1-R2-PD was measured para-sagittally in the frontal lobe, anterior to the genu of the corpus callosum where the cortex is thick (left and right); CSF was measured in the frontal horns (left and right); WM was measured peri-triagonally (left and right), frontally (left and right), parietally in the centrum semiovale (left and right), and in the genu and splenium of the corpus callosum. The average ± SD for R1, R2 and PD of the tissue types was calculated from all the voxels in the ROIs, for each tissue type. For WM, each tissue region was calculated separately, as was an average of all eight WM ROIs.

Prior Knowledge - Partial Volume Simulations

In order to estimate the characteristics of partial volume voxels in R1-R2-PD space, required as prior knowledge in the segmentation procedure, a numerical Bloch simulator that could simulate mixtures of different tissue types was used. Mathematical details are described in Appendix S1.

The Bloch simulator was used to estimate the observed R1, R2 and PD values for combinations of two tissues in the transitions WM↔GM, WM↔CSF and GM↔CSF, for both 1.5 T and 3.0 T field strengths separately. When two tissues were combined, the tissue content was altered in incremental steps of 0.1% to cover a range of 0 to 100% of each tissue type. The Bloch simulation with the QRAPMASTER sequence was executed three times for each of the two field strengths, with R1 and R2 values from the three previously measured pure brain tissue clusters, in order to acquire the magnetization states, Etissue, for WM, GM and CSF. For each specific partial volume mixture of two tissues, the magnetisation state matrices, Etissue1 and Etissue2, were combined and weighted with the PD and the volume fraction, α, of the tissues according to:

Ecombined=α*PDtissue1*Etissue1+(1α)*PDtissue2*Etissue2(1)

This combination of two magnetization states yields the sum of the two tissues assuming that each tissue generates signals independent of the other and that there is no interaction between the tissues in this approximation.

In order to retrieve R1, R2 and PD for the mixture the QRAPMASTER fitting in the software SyMRI Diagnostics was used, for additional details on the fitting process see 22.

The partial volume simulations were all included in one Monte Carlo simulation and repeated 10,000 times, with the SNR set to 20 (corresponding to the noise setting suggested in Ref. [24]).

As tissue fractions were altered in 0.1% steps in the three partial volume transitions, the Monte Carlo simulation resulted in 3,000 multivariate tissue clusters in the R1–R2–PD space, for both field strengths separately. These clusters were tested for normal distribution, using skewness and kurtosis tests. The skewness is a measure of the asymmetry of a distribution, and the kurtosis is a measure of the relative concentration of data values, in the centre versus in the tails, of a distribution. For normal distribution the skewness and kurtosis are zero, but any value below 2 suggests that the data are normally distributed [34]. For normal distribution data, covariance matrices were calculated. These simulated tissue mixtures were used to describe partial volume characteristics in the segmentation procedure.

Brain Tissue Segmentation

Once the pure tissue clusters and the relation between partial volume tissue composition and corresponding R1, R2, PD values were established, this prior knowledge was used to construct tissue ‘lookup grids’ in the R1-R2-PD space, for 1.5 T and 3.0 T. The lookup grids translated in vivo qMRI data into tissue volume fractions of WM, GM and CSF using the R1, R2 and PD values by comparison to the prior knowledge. The lookup grids were constructed such that in vivo voxels with R1, R2 and PD values within the 95% prediction ellipses of the measured pure tissue clusters for WM, GM and CSF, were mapped to 100% pure tissue. All the remaining simulated clusters, with different partial volume fractions of two tissues, were overlaid on the lookup grid. Since the change in voxel mixture in the simulations were performed using 0.1% incremental steps, these clusters appeared close together and overlapped. In vivo voxels with R1, R2 and PD values within the 95% prediction ellipses of the simulated partial-volume tissue clusters were mapped to the tissue fractions of the closest simulated cluster. Remaining in vivo voxels, outside the defined lookup grids in R1-R2-PD space, were considered undefined and termed ‘Non-WM/GM/CSF’ (or ‘NON’). Brain segmentation was implemented in C++ as a module to SyMRI Diagnostics.

Statistical Analysis

For all acquisitions at 1.5 T and 3.0 T, volumes of WM, GM, CSF and NON were calculated using the corresponding lookup grid. The ‘Brain Parenchymal Volume’ (BPV) was calculated as the sum of WM, GM and NON. The ‘Intracranial Volume’ (ICV) was calculated as the sum of BPV and CSF, and finally the ‘Brain Parenchymal Fraction’ (BPF) was calculated as the ratio of BPV to ICV. Volume measurements were tested for normal distribution using the Shapiro-Wilk’s test. For images acquired in the first acquisition of each session the within-subject standard deviation, sw, was estimated as the square root of the mean within-subject variance. The repeatability was calculated as 2.77 * sw, as suggested by Bland and Altman, see 35. This definition of ‘repeatability’ was based on the assumption that the difference between any two measurements of the same subject is expected to be less than this for 95% of pairs of observations. Finally, the within-subject coefficient of variation (CV) was calculated.

The General Linear Model (GLM) was used to examine the effect of field strength on the determined volume measures. Field strength, subject and session were included as fixed factors in the model. Analyses of main effects were performed where significance was found. Statistical analyses were performed using SPSS Statistics 19 (IBM, USA, 2010).

In order to investigate regional differences in tissue characteristics between the 1.5 T and 3.0 T scanners, voxel-wise t-tests of each tissue type were performed using SPM8 (Wellcome Department of Imaging Neuroscience, University College, London, UK). Before calculating the statistics, maps of WM, GM, and CSF were normalized to a standard stereotactic space in Montreal Neurological Institute (MNI) coordinates using the T2-weighted template available in SPM8. Synthetic T2-weighted images, which were smoothed with an 8 mm Gaussian kernel to reduce the individual anatomical details, were used as source images to calculate the transformation matrix. A 12-parameter (translation, rotation, shear, zoom) affine registration followed by nonlinear deformations, defined by a linear combination of three dimensional discrete cosine basis functions, were used to transform the synthetic T2-weighted images of each subject to the template. The resulting transformation matrices were then applied to the WM, GM, and CSF maps. The resulting maps were re-gridded to 2 x 2 x 2 mm3 voxel size to obtain an isotropic dataset. Voxel-wise differences between maps acquired at 1.5 T and 3.0 T were estimated by two-tailed paired t-tests of each tissue type. The resulting t-maps were initially thresholded at p=0.001, uncorrected for multiple comparisons. Based on those initially thresholded maps, ROI analyses were made of regions representing the different divisions of the brain: the frontal, temporal, parietal, occipital, and the limbic lobes; the latter including the cingulate cortex and the hippocampal formation; as well as the cerebellum (anterior and posterior lobes), pons, medulla, midbrain, i.e. the upper brain stem, and the sub-lobar regions comprising the thalamus and the basal ganglia. The results are reported as significant if cluster p<0.001, Family Wise Error (FWE) corrected, and cluster size was at least 20 voxels.

Results

Prior Knowledge

The presented brain tissue segmentation method required prior knowledge on pure brain tissue cluster positions, which were measured using ROIs. Results from these measurements are presented in Table 1 for both 1.5 T and 3.0 T. These values were used in the subsequent partial volume simulations and in order to create the lookup grids, which subsequently were used to segment the in vivo data.

1.5 T3.0 T
AnatomynR1 (s-1)R2 (s-1)PD (%)R1 (s-1)R2 (s-1)PD (%)
Frontal white matter21.67 ± 0.0913.55 ± 0.6364.86 ± 2.041.39 ± 0.0614.49 ± 0.5766.31 ± 1.46
Peritrigonal white matter21.64 ± 0.1212.66 ± 0.9062.03 ± 3.071.41 ± 0.0813.57 ± 0.8963.14 ± 2.90
Centrum Semiovale21.51 ± 0.0811.80 ± 0.5668.53 ± 1.751.31 ± 0.0412.77 ± 0.5668.27 ± 1.51
Genu11.78 ± 0.1414.23 ± 0.7359.88 ± 2.261.46 ± 0.0715.29 ± 0.5863.72 ± 1.46
Splenium11.58 ± 0.1212.46 ± 0.7965.33 ± 2.921.38 ± 0.0913.11 ± 0.8365.61 ± 2.37
Average white matter81.62 ± 0.1412.84 ± 1.0964.51 ± 3.771.38 ± 0.0813.76 ± 1.1065.60 ± 2.84
Frontal Cortex20.81 ± 0.0810.29 ± 0.5988.06 ± 3.760.77 ± 0.0810.69 ± 0.7686.88 ± 2.86
CSF20.24 ± 0.010.87 ± 0.55103.76 ± 3.540.32 ± 0.090.79 ± 0.5994.40 ± 6.59

Table 1. qMRI tissue parameters from manually placed regions of interest in 10 healthy subjects, on both 1.5 T and 3.0 T systems.

n indicates the number of ROI’s used per region in each subject.

CSV
Download CSV

The partial volume simulations resulted in 3,000 R1-R2–PD clusters for each field strength, one for each specific mixture of two tissue types in the transitions WM↔GM, WM↔CSF and GM↔CSF. Results from the simulations at both field strengths are shown in Figure 2. In this figure only the pure tissue clusters and 50% mixtures are presented. These clusters are shown as projections on the R1-R2, R1-PD and R2-PD planes, the indicated ellipses enclose 95% of the R1–R2–PD estimates from the Monte Carlo simulation. The observed clusters are not tilted in the R1–R2 plot, suggesting that the correlation between R1 and R2 is low. In the R1–PD and the R2–PD plots; however, the clusters are tilted, which indicates a higher correlation of R1 and R2 with PD. Skewness and kurtosis measurements indicated that the tissue clusters were normal distributed in the R1-R2-PD space.

thumbnail
Figure 2. Pure tissue clusters WM, GM, CSF, and the mixtures of 50% WM/GM, 50% WM/CSF and 50% GM/CSF shown in the R1-R2-PD space.

Projections on a) R1-R2, b) R1-PD and c) R2-PD planes. The contour lines show the 95% prediction ellipses for simulations with 20. Coloured ellipses are for 3.0 T and corresponding grey clusters are for 1.5 T. This figure shows the differences in tissue characteristics between 1.5 T and 3.0 T field strengths.

https://doi.org/10.1371/journal.pone.0074795.g002

Brain Tissue Segmentation

The lookup grids, which were subsequently used in the brain tissue segmentation, were created from the measured pure tissue cluster positions in R1-R2-PD space and the partial volume tissue simulations. When the 95% prediction ellipses of the 3,000 simulated partial volume clusters were joined together in the R1-R2-PD space, curved bands were formed. In Figure 3 These lookup grids are shown for both field strengths. The significance of the bands is that 95% of brain voxels containing WM, GM and CSF, acquired with an SNR of 20, were expected to fall within the indicated region. Figure 3 also shows results from partial volume simulations without added noise. These results are shown as curves between the mean values of the pure tissue clusters.

thumbnail
Figure 3. Lookup grids formed by joining together the 95% prediction ellipses of the 3,000 simulated partial-volume clusters (dashed lines).

The transition GM↔WM is shown in green, WM↔CSF in red and WM↔CSF in blue. The top panels a-c shows the results for 1.5 T and the bottom panels d-f shows results for 3.0 T. The solid lines represent the observations in the absence of noise where each diamond corresponds to a change of 20% tissue fractions and each dot correspond to a change of 10% tissue fractions. R1-R2-PD space is projected on a, d) R1-R2, b, e) R1-PD and c, f) R2-PD planes.

https://doi.org/10.1371/journal.pone.0074795.g003

Representative results from WM, GM and CSF segmentation of the in vivo data acquired at both field strengths are shown in Figure 4. By visual inspection, the fully automated segmentation follows the contours of each tissue type as visualized in the T2-weighted images in the left panel. Table 2 lists segmentation results from the 1.5 T scanner, including the repeatability expressed in mL of tissue, and table 3 lists the corresponding segmentation results from the 3.0 T system. Within-subject standard deviation, repeatability and the coefficient of variation are also reported. Tissue volumes were normally distributed as assessed by the Shapiro-Wilk test. Figure 5 shows Bland-Altman plots comparing the 1.5 T and 3.0 T measurements for each subject and tissue type separately.

thumbnail
Figure 4. Brain tissue segmentation, of one slice in one healthy subject, at both 1.5 T (top panel) and 3.0 T (lower panel).

From left to right: T2-weighted conventional image, white matter (blue), grey matter (green) and CSF (purple). The red lines are the brain intracranial mask calculated automatically by the SyMRI software.

https://doi.org/10.1371/journal.pone.0074795.g004

Subject (sex/age)WM (mL)GM (mL)CSF (mL)NON (mL)ICV (mL)BPV (mL)BPF (%)
F21383 ± 0679 ± 280 ± 09 ± 01151 ± 21071 ± 293.05 ± 0.07
M25568 ± 0840 ± 1177 ± 012 ± 01598 ± 11420 ± 188.90 ± 0.00
M23556 ± 3839 ± 0170 ± 015 ± 01580 ± 21410 ± 289.20 ± 0.00
M21521 ± 7824 ± 7130 ± 111 ± 01486 ± 11357 ± 191.25 ± 0.07
F29540 ± 3697 ± 4155 ± 113 ± 01405 ± 21250 ± 288.95 ± 0.07
F24505 ± 3774 ± 0104 ± 013 ± 11396 ± 21292 ± 292.55 ± 0.07
M23596 ± 1911 ± 3206 ± 122 ± 11735 ± 01530 ± 288.15 ± 0.07
F26452 ± 8723 ± 9131 ± 28 ± 01314 ± 11183 ± 190.10 ± 0.14
F27448 ± 3713 ± 4156 ± 010 ± 01328 ± 21171 ± 188.20 ± 0.00
F25564 ± 1768 ± 1166 ± 312 ± 11511 ± 11345 ± 189.00 ± 0.14
sw3.754.231.190.501.661.720.08
Repeatability10.3811.733.291.394.594.760.22
CV0.00570.00430.00590.03290.00110.00130.0007

Table 2. Segmentation results for the 10 healthy subjects at the 1.5 T system.

The columns list mean value and standard deviation from the first scan in the two sessions, for each subject. The last rows tabulate the within-subject standard deviation, sw, the repeatability and the coefficient of variation.
CSV
Download CSV
Subject (sex/age)WM (mL)GM (mL)CSF (mL)NON (mL)ICV (mL)BPV (mL)BPF (%)
F21438 ± 6644 ± 457 ± 19 ± 01149 ± 41092 ± 295.00 ± 0.14
M25633 ± 9807 ± 13137 ± 08 ± 01586 ± 51448 ± 491.35 ± 0.07
M23645 ± 25763 ± 6142 ± 618 ± 11568 ± 131426 ± 2090.95 ± 0.49
M21592 ± 4769 ± 6104 ± 114 ± 01480 ± 111375 ± 1093.00 ± 0.00
F29611 ± 5659 ± 8127 ± 912 ± 21409 ± 81282 ± 191.00 ± 0.57
F24564 ± 7746 ± 1481 ± 210 ± 21402 ± 31321 ± 694.20 ± 0.14
M23671 ± 1871 ± 20167 ± 1523 ± 31732 ± 21565 ± 1790.35 ± 0.92
F26517 ± 17680 ± 4107 ± 08 ± 01312 ± 131205 ± 1391.85 ± 0.07
F27509 ± 4684 ± 2127 ± 99 ± 11329 ± 41202 ± 690.45 ± 0.64
F25620 ± 3741 ± 8135 ± 411 ± 01508 ± 11373 ± 591.05 ± 0.21
Sw10.599.956.781.347.8310.260.44
Repeatability29.3327.5718.773.7121.7028.421.21
CV0.01410.01110.03710.07160.00450.00610.0036

Table 3. Segmentation results for the 10 healthy subjects at the 3.0 T system.

The columns list mean value and standard deviation from the first scan in the two sessions, for each subject. The last rows tabulate the within-subject standard deviation, sw, the repeatability and the coefficient of variation.
CSV
Download CSV
thumbnail
Figure 5. Bland-Altman plots comparing the volumes between 1.5 T and 3.0 T.

Volumes of a) WM, b) GM, c) CSF and d) ICV for all 10 subjects measured at 1.5 T and 3.0 T. The dashed lines indicate the limits of agreement and the mean of all measurements (bias).

https://doi.org/10.1371/journal.pone.0074795.g005

GLM revealed statistically significant differences in total volume, between 1.5 T and 3.0 T, for WM (p < 0.001), GM (p < 0.001), CSF (p < 0.001), BPV (p < 0.001), and BPF (p < 0.001) but not for NON and ICV. Analyses of main effects showed that WM was underestimated, while GM and CSF were overestimated on 1.5 T compared to 3.0 T. The absolute differences and confidence intervals are listed in Table 4.

Dependent VariableMean Difference (mL) (1.5 T-3.0 T)RSMEMean Difference (% of 1.5 T)(1.5 T-3.0 T)95% Confidence Interval of Difference (mL)
WM-6668-12.85%[-72, -61]
GM40435.15% [34,47]
CSF293019.66% [26,33]
NON020%[-1, 1]
ICV370.21%[-1, 7]
BPV-2627-2.00%[-31, -22]
BPF (%)-1.992.00-2.21%[-2.18, -1.79]

Table 4. Results from pair-wise comparisons of tissue volumes between the 1.5 T and the 3.0 T system.

CSV
Download CSV

Relating the repeatability for each tissue type to the whole ICV (repeatability/ICV) is a relative measurement of repeatability with respect to the complete brain volume. Results for the 1.5 T system were 0.72% for WM, 0.81% for GM, 0.23% for CSF, 0.10% for NON, 0.32% for ICV and 0.33% for BPV. The corresponding numbers for 3.0 T were 2.03% for WM, 1.91% for GM, 1.30% for CSF, 0.26% for NON, 1.50% for ICV and 1.96% for BPV. The repeatability for BPF at 1.5 T was 0.22% and for 3.0 T 1.21%.

Regional Analysis

The voxel-wise comparisons of tissue maps at 1.5 T and 3.0 T revealed regional differences in WM, GM and CSF volume estimations. In general, differences occurred in deep brain structures where tissue classified as mostly GM at 1.5 T was instead classified as mostly WM at 3.0 T. Statistically significant differences for CSF were only found in one small region close to the midline anterior cingulate gyrus in the limbic lobe (cluster size 26 voxels). Statistically significant regional differences in WM and GM estimation between 1.5 T and 3.0 T were found in bilateral and symmetrical regions in the following structures: the cerebellum, the anterior parts of the hippocampal formation, dorsal parts of the medulla and pons, midbrain structures comprising the thalamus, substantia nigra, red nucleus, and the dorsal midbrain tegmentum, as well as in several clusters in sub-lobar white and grey matter including the putamen. Figure 6 shows the results of the voxel-wise t-test where the WM estimate was significantly larger at 3.0 T compared to 1.5 T, the regions were GM estimates was significantly larger at 1.5 T was virtually identical.

thumbnail
Figure 6. Results from voxel-wise t-tests comparing brain segmentation from all 10 subjects between 1.5 T and 3.0 T.

Statistically significant differences are indicated where WM volume is estimated to be larger on the 3.0 T system.

https://doi.org/10.1371/journal.pone.0074795.g006

Discussion

In this study, one brain tissue segmentation and volume estimation method, previously reported for 1.5 T, was applied at both 1.5 T and 3.0 T Philips MR systems. It was shown that it was necessary to implement separate lookup grids (Figure 3) for the two field strengths since the tissue characteristics for WM, GM and CSF have varying R1-R2-PD values depending on main magnetic field strength (Table 1). Since the differences in tissue characteristics between the two field strengths were compensated for in the segmentation method, it was expected that tissue volumes would be similar at both field strengths, but variations were found. These may be caused by differences in regional tissue composition affecting the R1 and R2 values, differences in ROI placement in the field strengths or differences in the qMRI sequence implementation. These factors will be discussed in detail below.

First, regional tissue characteristics vary between 1.5 T and 3.0 T in a non-linear manner. This has been attributed to iron deposits in different anatomical regions [26,28], affecting mainly the measured R2, but also R1. Schenck et al. [28] showed that contrast between iron-rich and iron-depleted regions is substantially enhanced when the main magnetic field strength is increased (from 1.5 T to 4.0 T). This suggests that the differences in quantitative measurements between iron-rich and iron-depleted regions are also enhanced with field strength. Schenck et al. also found that the R2 dependence on field strength is greatest in regions with high iron concentration, and pointed out that the substantia nigra and the red nuclei are prominent indicators of this phenomenon. These structures were also found to differ using the present segmentation method. Schenck et al. showed a higher increase in T1 and T2 in the globus palladius, red nucleus, substantia nigra, putamen, thalamus and the splenium compared to the relatively iron-depleted regions of the midbrain. These discrepancies could possibly explain the misclassifications of the present algorithm in these particular regions. When determining tissue characteristics for the bulk WM, GM and CSF, as was done in the present study, local differences in R1 and R2 depending on anatomical structures were not included. This could lead to regional misclassifications when comparing field strengths. In this study the voxels which where differently classified in 1.5 T and 3.0 T closely resembled these reported anatomical regions. In a similar fashion, Pfefferbaum et al. [27] found regional volume differences between 1.5 T and 3.0 T, where some structures were measured as being up to 15% larger at 3.0 T, while other structures were instead down to 19% smaller, when using an SRI24 atlas-based approach. Pfefferbaum et al. thus concluded that regional volume differences were not due to simple (global) volumetric scaling. In another study Jovicich et al. [31] found, when using the FreeSurfer tool that both the magnitude and sign of regional volume differences depend on brain structure and field strength.

Second, since tissue characteristics vary between 1.5 T and 3.0 T, the appearance of conventional contrast-weighted MR images also varies. This may have led to differently interpreted boundaries between tissue types when ROIs were manually placed to obtain the prior knowledge of pure tissue clusters. The present method is to a degree sensitive to correctly measured pure tissue clusters and even small variations in ROI positioning on the two field strengths may have an effect on the total tissue volumes, calculated by the segmentation method. Identical ROI placement may be especially difficult to achieve in cortical GM, which is a relatively thin layer between WM and CSF. In one study, Han et al. [32] found that cortical thickness was detected to be 0.11 mm thicker in 3.0 T compared to 1.5 T in 15 subjects; something that suggests that GM ROIs may be easier to position at 3.0 T, compared to 1.5 T, and the pure tissue cluster for GM used in this study may differ between the two field strengths.

Third, even though the QRAPMASTER sequence is identically implemented at both field strengths, small differences may still exist due to such effects as a larger water-fat shift at 3.0 T, and slightly higher geometrical distortion due to the EPI readout. Moreover, B1-field inhomogeneities are larger on the 3.0 T MR scanner; thus residual inhomogeneities may have an effect.

The largest difference was observed in WM, while GM and CSF differed less between the two field strengths (Table 4). WM differences were observed up to 13%, GM up to 6% and CSF up to 20%. Since overestimating WM and underestimating GM have opposite effects on the BPV, the differences in BPV were notably smaller, only about 2%, although still statistically significant. The difference in ICV did not reach statistical significance between the field strengths. The BPF, which is the ratio of the BPV to the ICV, including CSF, showed a significantly higher value at 3.0 T. Even so, the results from this study indicated that the present segmentation method may reduce the differences in volume measurements between 1.5 T and 3.0 T system compared to FreeSurfer [31], SPM and Bayesian methods [33] for whole brain volumes.

In addition to the differences observed between the two scanners, we also investigated the repeatability in the two systems separately (Table 2 and Table 3). These results show that the repeatability is larger at 1.5 T for all measurements. This may be caused by the higher dynamic range in R1-R2 values on this system leading to greater separation of pure tissue clusters. Also, the MR head coil in the 3.0 T system was somewhat larger, allowing more movement by the subject, which could potentially lead to higher motion artefacts, affecting the repeatability.

Limitations of This Study

This study was conducted on a relatively low number of healthy subjects in the age range of 20 to 30 years. The included subjects showed no signs of the degeneration or disease in the brain that would be expected in an elderly group or patient group. Such degeneration includes changes due to tissue atrophy, in particular around the ventricles, in addition to small, unspecific age-induced WM changes. It would therefore be interesting to extend the validation to other groups in future research.

This study was conducted using two different MR scanners. While effectively comparing these two scanners of different field strengths, the site-specific effect of the scanners themselves could not be separated in the statistical analysis. There is a possibility that some of the observed differences were not exclusively due to the field strength, but due to the individual scanners in question. Differences in scanners may include inhomogeneities, different gradient systems and differences in scanner design, system release and coils.

Conclusions

The present study investigated the extension of a recently published brain tissue segmentation method, based on qMRI, to include both 1.5 T and 3.0 T field strengths. Validation was performed in 10 young and healthy subjects. The results showed that when compensating for global R1 and R2 differences, by implementation of separate lookup grids, most of the brain segmentation was identical at both field strengths and each scanner had high inner-scanner repeatability. However, regional differences in segmentation results were found in deep brain structures, the cerebellum, and the brain stem. These regional differences led to volume differences for each tissue type in the order of 40-70 mL (corresponding to about 5-20%) and differences in total brain volume of about 25 mL (2%). The total brain volume difference was less than that of the individual tissues, since differences in WM and GM had different sign. Therefore the total brain volume, where WM and GM were added together, was a more robust measure across field strengths. We hypothesize that some of the differences may have been caused by iron deposition, which previously has been shown to affect R1 and R2 values, in a non-linear manner. The present method is highly promising for application at both field strengths, but the local differences should preferably be accounted for before implementing the method in clinical practice or before combining field strengths in research studies. In future research it would be interesting to extend the present method using an atlas-based approach, with a locally adapting segmentation. Furthermore, as we have only examined relatively young adult healthy subjects, it would be of interest to also validate the segmentation method in selected groups with elderly subjects or in patients suffering from different neurological conditions.

Supporting Information

Appendix S1.

Mathematical Details of Bloch Simulation.

https://doi.org/10.1371/journal.pone.0074795.s001

(DOCX)

Acknowledgments

We acknowledge the very valuable contribution and excellent assistance of Dr. Anders Tisell, who helped in recruiting, assessing and scanning the subjects.

Author Contributions

Conceived and designed the experiments: JW IB ME. Performed the experiments: JW IB ME JBMW. Analyzed the data: JW ME. Contributed reagents/materials/analysis tools: JW JBMW. Wrote the manuscript: JW IB ME JBMW PL.

References

  1. 1. Zipp F (2009) A new window in multiple sclerosis pathology: non-conventional quantitative magnetic resonance imaging outcomes. J Neurol Sci 287 Suppl 1: S24-S29. doi:https://doi.org/10.1016/S0022-510X(09)71297-3. PubMed: 20106345.
  2. 2. Schmidt-Wilcke T (2008) Variations in brain volume and regional morphology associated with chronic pain. Curr Rheumatol Rep 10: 467-474. doi:https://doi.org/10.1007/s11926-008-0077-7. PubMed: 19007538.
  3. 3. Bozzali M, Cercignani M, Caltagirone C (2008) Brain volumetrics to investigate aging and the principal forms of degenerative cognitive decline: a brief review. Magn Reson Imaging 26: 1065-1070. doi:https://doi.org/10.1016/j.mri.2008.01.044. PubMed: 18436405.
  4. 4. Rovira A, León A (2008) MR in the diagnosis and monitoring of multiple sclerosis: an overview. Eur J Radiol 67: 409-414. doi:https://doi.org/10.1016/j.ejrad.2008.02.044. PubMed: 18434066.
  5. 5. Simmons A, Westman E, Muehlboeck S, Mecocci P, Vellas B et al. (2009) MRI measures of Alzheimer’s disease and the AddNeuroMed study. Ann N Y Acad Sci 1180: 47-55. doi:https://doi.org/10.1111/j.1749-6632.2009.05063.x. PubMed: 19906260.
  6. 6. Barnes J, Ourselin S, Fox NC (2009) Clinical application of measurement of hippocampal atrophy in degenerative dementias. Hippocampus 19: 510-516. doi:https://doi.org/10.1002/hipo.20617. PubMed: 19405145.
  7. 7. Jack CR Jr., Slomkowski M, Gracon S, Hoover TM, Felmlee JP et al. (2003) MRI as a biomarker of disease progression in a therapeutic trial of milameline for AD. Neurology 60: 253-260. doi:https://doi.org/10.1212/01.WNL.0000042480.86872.03. PubMed: 12552040.
  8. 8. Moncrieff J, Leo J (2010) A systematic review of the effects of antipsychotic drugs on brain volume. Psychol Med 40: 1409-1422. doi:https://doi.org/10.1017/S0033291709992297. PubMed: 20085668.
  9. 9. Alfano B, Brunetti A, Covelli EM, Quarantelli M, Panico MR et al. (1997) Unsupervised, automated segmentation of the normal brain using a multispectral relaxometric magnetic resonance approach. Magn Reson Med 37: 84-93. doi:https://doi.org/10.1002/mrm.1910370113. PubMed: 8978636.
  10. 10. Andersen AH, Zhang Z, Avison MJ, Gash DM (2002) Automated segmentation of multispectral brain MR images. J Neurosci Methods 122: 13-23. doi:https://doi.org/10.1016/S0165-0270(02)00273-X. PubMed: 12535761.
  11. 11. Jackson EF, Narayana PA, Falconer JC (1994) Reproducibility of nonparametric feature map segmentation for determination of normal human intracranial volumes with MR imaging data. J Magn Reson Imaging 4: 692-700. doi:https://doi.org/10.1002/jmri.1880040512. PubMed: 7981514.
  12. 12. Kikinis R, Shenton ME, Gerig G, Martin J, Anderson M et al. (1992) Routine quantitative analysis of brain and cerebrospinal fluid spaces with MR imaging. J Magn Reson Imaging 2: 619-629. doi:https://doi.org/10.1002/jmri.1880020603. PubMed: 1446105.
  13. 13. Bonar DC, Schaper KA, Anderson JR, Rottenberg DA, Strother SC (1993) Graphical analysis of MR feature space for measurement of CSF, gray-matter, and white-matter volumes. J Comput Assist Tomogr 17: 461-470. doi:https://doi.org/10.1097/00004728-199305000-00024. PubMed: 8491912.
  14. 14. Choi HS, Haynor DR, Kim Y (1991) Partial volume tissue classification of multichannel magnetic resonance images-a mixel model. IEEE Trans Med Imaging 10: 395-407. doi:https://doi.org/10.1109/42.97590. PubMed: 18222842.
  15. 15. Shattuck DW, Sandor-Leahy SR, Schaper KA, Rottenberg DA, Leahy RM (2001) Magnetic resonance image tissue classification using a partial volume model. Neuroimage 13: 856-876. doi:https://doi.org/10.1016/S1053-8119(01)92198-6. PubMed: 11304082.
  16. 16. Clare S, Jezzard P (2001) Rapid T(1) mapping using multislice echo planar imaging. Magn Reson Med 45: 630-634. doi:https://doi.org/10.1002/mrm.1085. PubMed: 11283991.
  17. 17. Deoni SC (2009) Transverse relaxation time (T2) mapping in the brain with off-resonance correction using phase-cycled steady-state free precession imaging. J Magn Reson Imaging 30: 411-417. doi:https://doi.org/10.1002/jmri.21849. PubMed: 19629970.
  18. 18. Deoni SC, Peters TM, Rutt BK (2005) High-resolution T1 and T2 mapping of the brain in a clinically acceptable time with DESPOT1 and DESPOT2. Magn Reson Med 53: 237-241. doi:https://doi.org/10.1002/mrm.20314. PubMed: 15690526.
  19. 19. Deoni SC, Rutt BK, Peters TM (2003) Rapid combined T1 and T2 mapping using gradient recalled acquisition in the steady state. Magn Reson Med 49: 515-526. doi:https://doi.org/10.1002/mrm.10407. PubMed: 12594755.
  20. 20. Neeb H, Zilles K, Shah NJ (2006) A new method for fast quantitative mapping of absolute water content in vivo. NeuroImage 31: 1156-1168. doi:https://doi.org/10.1016/j.neuroimage.2005.12.063. PubMed: 16650780.
  21. 21. Warntjes JB, Dahlqvist O, Lundberg P (2007) Novel method for rapid, simultaneous T1, T*2, and proton density quantification. Magn Reson Med 57: 528-537. doi:https://doi.org/10.1002/mrm.21165. PubMed: 17326183.
  22. 22. Warntjes JB, Leinhard OD, West J, Lundberg P (2008) Rapid magnetic resonance quantification on the brain: Optimization for clinical usage. Magn Reson Med 60: 320-329. doi:https://doi.org/10.1002/mrm.21635. PubMed: 18666127.
  23. 23. Zhu DC, Penn RD (2005) Full-brain T1 mapping through inversion recovery fast spin echo imaging with time-efficient slice ordering. Magn Reson Med 54: 725-731. doi:https://doi.org/10.1002/mrm.20602. PubMed: 16086307.
  24. 24. West J, Warntjes JB, Lundberg P (2012) Novel whole brain segmentation and volume estimation using quantitative MRI. Eur Radiol 22: 998-1007. doi:https://doi.org/10.1007/s00330-011-2336-7. PubMed: 22113264.
  25. 25. Bottomley PA, Hardy CJ, Argersinger RE, Allen-Moore G (1987) A review of 1H nuclear magnetic resonance relaxation in pathology: are T1 and T2 diagnostic? Med Phys 14: 1-37. doi:https://doi.org/10.1118/1.596111. PubMed: 3031439.
  26. 26. Hocq A, Brouette N, Saussez S, Luhmer M, Gillis P et al. (2009) Variable-field relaxometry of iron-containing human tissues: a preliminary study. Contrast Media. Mol Imaging 4: 157-164.
  27. 27. Pfefferbaum A, Rohlfing T, Rosenbloom MJ, Sullivan EV (2012) Combining atlas-based parcellation of regional brain data acquired across scanners at 1.5 T and 3.0 T field strengths. Neuroimage 60: 940-951. doi:https://doi.org/10.1016/j.neuroimage.2012.01.092. PubMed: 22297204.
  28. 28. Schenck JF (1995) Imaging of brain iron by magnetic resonance: T2 relaxation at different field strengths. J Neurol Sci 134 Suppl: 10-18. doi:https://doi.org/10.1016/0022-510X(95)00203-E. PubMed: 8847539.
  29. 29. Fischl B (2012) FreeSurfer. NeuroImage 62: 774-781. doi:https://doi.org/10.1016/j.neuroimage.2012.01.021. PubMed: 22248573.
  30. 30. Zhang Y, Brady M, Smith S (2001) Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging 20: 45-57. doi:https://doi.org/10.1109/42.906424. PubMed: 11293691.
  31. 31. Jovicich J, Czanner S, Han X, Salat D, van der Kouwe A et al. (2009) MRI-derived measurements of human subcortical, ventricular and intracranial brain volumes: Reliability effects of scan sessions, acquisition sequences, data analyses, scanner upgrade, scanner vendors and field strengths. Neuroimage 46: 177-192. doi:https://doi.org/10.1016/j.neuroimage.2009.02.010. PubMed: 19233293.
  32. 32. Han X, Jovicich J, Salat D, van der Kouwe A, Quinn B et al. (2006) Reliability of MRI-derived measurements of human cerebral cortical thickness: the effects of field strength, scanner upgrade and manufacturer. NeuroImage 32: 180-194. doi:https://doi.org/10.1016/j.neuroimage.2006.02.051. PubMed: 16651008.
  33. 33. Keihaninejad S, Heckemann RA, Fagiolo G, Symms MR, Hajnal JV et al. (2010) A robust method to estimate the intracranial volume across MRI field strengths (1.5T and 3T). Neuroimage 50: 1427-1437. doi:https://doi.org/10.1016/j.neuroimage.2010.01.064. PubMed: 20114082.
  34. 34. Groeneveld RA, Meeden G (1984) Measuring Skewness and Kurtosis. Stat 33: 9.
  35. 35. Bland JM, Altman DG (1996) Measurement error. BMJ 313: 744. doi:https://doi.org/10.1136/bmj.313.7059.744a. PubMed: 8819450.