Original Research ARTICLE
Feasibility of simultaneous whole-brain imaging on an integrated PET-MRI system using an enhanced 2-point Dixon attenuation correction method
- 1Lawson Health Research Institute, London, ON, Canada
- 2Medical Biophysics, Western University, London, ON, Canada
- 3Diagnostic Imaging, St. Joseph's Health Care, London, ON, Canada
- 4Fraunhofer Institute for Medical Image Computing MEVIS, Bremen, Germany
Purpose: To evaluate a potential approach for improved attenuation correction (AC) of PET in simultaneous PET and MRI brain imaging, a straightforward approach that adds bone information missing on Dixon AC was explored.
Methods: Bone information derived from individual T1-weighted MRI data using segmentation tools in SPM8, were added to the standard Dixon AC map. Percent relative difference between PET reconstructed with Dixon+bone and with Dixon AC maps were compared across brain regions of 13 oncology patients. The clinical potential of the improved Dixon AC was investigated by comparing relative perfusion (rCBF) measured with arterial spin labeling to relative glucose uptake (rPETdxbone) measured simultaneously with 18F-flurodexoyglucose in several regions across the brain.
Results: A gradual increase in PET signal from center to the edge of the brain was observed in PET reconstructed with Dixon+bone. A 5–20% reduction in regional PET signals were observed in data corrected with standard Dixon AC maps. These regional underestimations of PET were either reduced or removed when Dixon+bone AC was applied. The mean relative correlation coefficient between rCBF and rPETdxbone was r = 0.53 (p < 0.001). Marked regional variations in rCBF-to-rPET correlation were observed, with the highest associations in the caudate and cingulate and the lowest in limbic structures. All findings were well matched to observations from previous studies conducted with PET data reconstructed with computed tomography derived AC maps.
Conclusion: Adding bone information derived from T1-weighted MRI to Dixon AC maps can improve underestimation of PET activity in hybrid PET-MRI neuroimaging.
The development of positron emission tomography (PET) and magnetic resonance imaging (MRI) hybrid systems that provide the combined advantages of high molecular sensitivity of PET and high spatial resolution of MRI, among other benefits, has led to ~400% increase in PET-MRI imaging publications in the last 4 years (Pubmed search “PET/MRI, PET-MRI, PET/MR, or PET-MR”). The clinical potential of this hybrid approach in neuroimaging extends beyond image fusion considering the ability of MRI to collect data associated with not only anatomy but also function, diffusion and metabolite concentrations (Pichler et al., 2010). In dementia, for instance, accuracy of early and differential diagnosis is improved with combined information from volumetric MRI and PET-specific tracers (Morbelli et al., 2010; Dukart et al., 2011), where MRI is not only used for anatomical localization of PET tracer but also to correct misclassification of gray matter voxels (Gutierrez et al., 2012) and to aid in retrospective motion correction of PET data (Catana et al., 2011). In neuro-oncology, metabolically active tumors can be more clearly delineated and better characterized with PET-MRI (Catana et al., 2012). In neuroscience, PET-MRI provides the ability to explore neurophysiologic processes such as neurovascular coupling or neuromodulation within one integrated dynamic model (Catana et al., 2012). These neurological applications are achieved with scanners that offer simultaneous hybrid imaging where diagnostic accuracy is maximized as image registration errors and total image acquisition times are minimized, and physiological and metabolic states are identical.
Accurate correction of PET signal attenuation in PET-MRI hybrid systems using information derived from MRI is challenging. Inherent limitation(s) of various proposed methods can lead to PET quantification errors of up to 20% (Keereman et al., 2013; Andersen et al., 2014; Dickson et al., 2014; Hitz et al., 2014), which limits adoption of PET-MRI in neuroimaging. Commercial PET-MRI hybrid systems rely on MRI-derived attenuation correction (MRAC) using a 2-point Dixon method (Dixon, 1984) based on segmentation of MR signals from water and fat into air, lung, fat and soft tissue (Coombs et al., 1997) but ignores bone, rendering this method less ideal for brain imaging. Various methods have been explored to produce optimal MRAC maps for neuroimaging that are as accurate as clinically accepted AC maps produced by computed tomography (CT). Ultra-short echo (UTE) MRI (Catana et al., 2010) where bone is imaged and incorporated into AC maps, can suffer from long imaging time (~5 times longer than Dixon - too long for whole-body imaging but still feasible for brain imaging), low spatial resolution and inaccurate signal segmentation (Dickson et al., 2014). Methods that employ CT templates/atlases mapped to individual MRI (Beyer et al., 2008; Marshall et al., 2013) or to guide MRI segmentation (Hofmann et al., 2008; Kops et al., 2009; Poynton et al., 2014), require complex algorithms that have limited reproducibility and at best still retain some level of inter-modality misregistration (Pappas et al., 2005), especially where large signal intensity differences exist between MRI and CT (e.g., cortical bone).
In this study, we explored the feasibility of simultaneous measurements of brain function with PET and MRI using an enhanced Dixon-based MRAC for PET signal attenuation correction in an integrated whole-body PET-MRI system. We combined bone segmentation from high-resolution three-dimensional (3D) T1-weighted anatomical images, routinely acquired in brain imaging, with the standard Dixon MRAC method implemented by the manufacturer. Bone segmentation was achieved with segmentation tools available from SPM8 (http://www.fil.ion.ucl.ac.uk). Comparison between the Dixon and enhanced Dixon method was made across a number of brain regions to assess regional differences and performance. In addition, PET 18F-fluoro-deoxyglucose (18F-FDG) data corrected with the enhanced method were compared to cerebral blood flow (CBF), measured simultaneously with 3D pseudo-continuous arterial spin labeling (pCASL), to investigate the potential clinical application of this MRAC method, for instance, in future neurodegenerative studies where hybrid PET-MRI is highly recommended (Drzezga et al., 2014).
Materials and Methods
This study was approved by the University Research Ethics Board, and written informed consent was obtained from all subjects. PET-MRI images were acquired from oncology patients recruited following the completion of a clinical PET-CT exam as part of their healthcare management. A total of 13 oncology patients were recruited 7 males and 6 females (58 ± 8 years old). All subjects had no history of neurological disorder and no radiation therapy or chemotherapy prior to imaging and were referred for whole-body PET-CT exam for staging of oncology. To minimize errors in MRI image quality from potential image artifacts produced by metal implants such as MR signal loss, only patients free of metallic implants and fixtures including dental implants were included in this study.
PET-MRI Image Acquisition
Whole-brain PET and MRI data were acquired simultaneously on a Siemens Biograph® mMR system. The PET subsystem consists of 8 rings of 56 detector blocks, each housing 8 × 8 lutetium oxyorthosilicate (LSO) crystals (size = 4 × 4 × 20 mm3) coupled to a 3 × 3 array of avalanche photodiodes (APDs) producing a transaxial field of view (FOV) of 59.4 cm and an axial FOV of 25.8 cm. The 3T MRI subsystem is similar to a Siemens 3T Verio with a maximum gradient strength of 45 mT/m and slew rate of 200 T/m/s, but a reduced bore size of 60 cm. For detailed description of the Biograph® mMR specifications cf. Delso et al. (2011). All MR images were acquired using a PET-compatible 16-channel phased array head (12-channel) and neck (4-channel) radiofrequency (RF) coil.
PET-MRI data were acquired immediately following the completion of a 20-min whole-body PET-CT scan and following 60-min uptake of a single intravenous injection of 18F-FDG (5 MBq/kg) administered in a dimly lit room. All subjects fasted for a minimum of 6 h and had blood glucose levels between 4.8 and 5.8 mmol/L. Clinical whole-body PET-CT were acquired following standard oncology imaging protocol from the level of the eyes to the mid-thigh. A series of MR images were acquired during a 15 min PET list-mode acquisition with subjects in supine position. The MRI data included:
(a) Coronal T1-weighted dual-echo 3D VIBE-Dixon for MRI Dixon-based attenuation correction of PET data [repetition time (TR) = 3.60 ms, echo times (TE) = 1.23 and 2.46 ms, flip angle = 10°, FOV = 500 × 328 mm2, 128 slices, voxel size = 4.1 × 2.6 × 3.1 mm3, bandwidth = 965 Hz/Px, acceleration factor of 2, and total acquisition (TA) = 0.19 s] (Martinez-Möller et al., 2009).
(b) Sagittal T1-weighted 3D magnetization-prepared rapid gradient-echo (MPRAGE) sequence (TR/TE: 2000/2.98 ms, inversion time (TI) = 900 ms, flip angle = 9°, FOV = 256 × 256 mm2, 176 slices, isotropic voxel size = 1.0 mm3, bandwidth = 240 Hz/Px, an acceleration factor of 2 and TA = 4.48 min).
(c) Transverse 3D single shot gradient-and-spin-echo (GRASE) pCASL sequence (Günther et al., 2005) (TR/TE = 3500/22.76 ms, FOV = 240 × 240 mm2, 24 slices, voxel size = 3.8 × 3.8 × 5 mm3, bandwidth = 2298 Hz/Px, and acceleration factor of 2). The pCASL label consisted of 1.5 s train of RF pulses applied 9 cm below the center of the imaging volume with a mean gradient of 0.6 mT/m. Sixty-four label and control pairs were acquired after a post-label delay of 1.2 s for TA = 7 min. Two nonselective inversion pulses were applied for background suppression during the post-labeling delay. Proton density images (M0) were acquired with the pCASL sequences for CBF quantification using a TR of 5 s, with no labeling or background suppression pulses.
MR attenuation map (μ-map) Generation and PET Image Reconstruction
Two MRAC segmentation methods were used for attenuation correction of the PET data: the standard 2-point Dixon method (Martinez-Möller et al., 2009) and 2-point Dixon method plus bone segmentation (Dixon+bone), which is described below. The Dixon MRAC μmaps were generated online from segmentation of fat and water signals in the in- and out-of-phase Dixon images into fat, soft tissue and air/background, as implemented by the manufacturer. Linear attenuation coefficients of 0, 0.086, and 0.10 cm−1 were assigned to air, fat and soft tissue, respectively (Martinez-Möller et al., 2009). Attenuation correction factors for hardware including the RF coil and scanning bed were included in the Dixon μ-map. The Dixon+bone MRAC μ-maps were generated offline using SPM8 and in-house MATLAB (2012a, The MathWorks, Natick, MA) scripts. This method overlays a bone mask created from bone segmentation of the MPRAGE image to a Dixon μ-map and applies a linear attenuation coefficient (μa) of 0.143 cm−1 (Catana et al., 2010), the μa value of bone at 511 kev, to the bone mask.
The bone mask for each subject was created as follows. (1) The T1-weighted MRPAGE dataset was coregistered to the Dixon in-phase image in SPM8 using rigid-body transformation with a normalized mutual information cost function to place the MPRAGE in the same voxel space as the Dixon μ-map. (2) The coregistered MPRAGE image set was segmented into gray matter, white matter, cerebrospinal fluid, bone, soft tissue and air/background probability tissue maps in native space using the new segment function in SPM8 and ICBM Tissue Probabilistic Atlases (http://www.loni.usc.edu/ICBM/). The new segment toolbox employs the unified segmentation method (Ashburner and Friston, 2005) that combines affine registration to ICBM atlas, bias field correction and segmentation within one integrated model. (3) The bone probability map was down-sampled to the Dixon images matrix size, smoothed with a 4 mm Gaussian filter (Wagenknecht et al., 2011) and thresholded for probabilities above 80% to minimize inclusion of non-bone signals. (4) The resulting binary image was eroded with morphological filtering and connected component analysis using a 3 × 3 × 3 voxel size of ones before assigning μa of 0.143 cm−1. Erosion was performed to minimize inclusion of non-bone voxels and the size of the structure elements were determined after empirical evaluation of bone voxels within a region known to have true bone content.
Each subject's PET list-mode data were reconstructed using the Siemens e7 tools to one image volume 344 × 344 × 127 matrix with an iterative algorithm (ordered-subsets expectation maximization: 3 iterations, 21 subsets; 3D Gaussian filter with a full width half maximum (FWHM) = 2.0 mm and 2.5 zoom factor) and corrected for decay, dead time, scatter and attenuation. The reconstruction was performed twice using either the Dixon (PETdx) or the Dixon+bone (PETdxbone)MRAC μ-map.
To allow for group comparison between PETdx and PETdxbone, the PET images were spatially normalized to the Montreal Neurological Institute (MNI) PET template in SPM8 (http://www.fil.ion.ucl.ac.uk) using an affine transformation and non-linear warps, and smoothed with a 6-mm FWHM Gaussian filter. The gray matter (GM) probability map for each subject was spatially normalized to the T1 MNI template using the unified segment approach in SPM, also smoothed with a 6-mm Gaussian filter, transformed to a binary mask at a threshold of 80% and applied to the PETdx and PETdxbone images.
The difference between the Dixon+bone and the Dixon MRAC method across the brain was computed as the percent relative difference (% RD) in mean activity concentration (Bq/ml) in thirteen a priori regions of interest (ROIs). . ROI brain masks were created using the automated anatomical library (Tzourio-Mazoyer et al., 2002) incorporated within WFU PickAtlas toolbox version 3.0 (Maldjian et al., 2003). The ROIs were frontal lobe, cingulate gyrus, insula, parietal lobe, temporal lobe, hippocampus, amygdala, thalamus, caudate, globus pallidus, putamen, and cerebellum. The PETdx and PETdxbone images were scaled by their respective global gray matter mean value to allow for voxel-by-voxel comparison between glucose uptake and blood flow, and since absolute quantification of PET-FDG with tracer kinetic modeling was not feasible in this study and generally not feasible in clinical practice. Regional %RD in relative mean activity was also calculated in the scaled PET data (rPET). To investigate potential regional differences, the mean values extracted from individual ROIs were compared between the two PET data sets using paired t-tests with SPSS 20.0 statistical software (IBM Corp. Armonk, NY, USA). Statistical significant differences were set at threshold p < 0.05. Evaluation of the regional differences between PETdx/PETdxbone and a ground truth method such as PET-CT were not made (Supplementary Material illustrates the comparison in one subject).
pCASL preprocessing was performed with SPM8 and scripts written in MATLAB. The pCASL time series were motion corrected, pair-wise subtracted, and time-averaged. The mean signal was registered to the MPRAGE images using a rigid-body transformation and converted to CBF using a single-compartment flow model (Wang et al., 2003). The CBF images were spatially normalized to the MNI template and smoothed using a Gaussian filter with a FWHM of 6 mm. A GM mask was applied to the CBF images that were normalized by the global mean. To determine the temporal signal-to-noise ratio (tSNR) of the perfusion-weighted time series when acquired with a PET-compatible head coil, tSNR was calculated for each subject as defined by the mean pixel signal in the whole brain relative to the mean pixel standard deviation.
Whole brain voxel-by-voxel independent samples t-tests restricted within a GM mask were performed on the individual relative images in SPM8 to investigate differences between rCBF and rPET images. Statistical significant differences were identified for clusters greater than 50 voxels at p < 0.05 after correcting for multiple comparisons using the False Discovery Rate (FDR). Voxel-by-voxel Pearson correlation was conducted on relative images averaged across subjects within the thirteen ROIs to correlate rCBF to rPETdx and to rPETdxbone.
A full description of the study demographics is listed in Table 1, including age, gender and diagnosis. No neurological lesions or gross neuropathological abnormalities were observed on PET-MRI brain images. Images were reviewed by a Board Certified Radiologist and Nuclear Medicine physician. Mean and standard deviation of the time from 18F-FDG injection to PET-MR examinations were 103 ± 11.18 min.
Evaluation of MRAC Methods
An example of Dixon+bone μ-map created by adding bone to a Dixon μ-map from a representative subject is displayed in Figure 1. A gradual radial increase in relative difference between PET images reconstructed with Dixon+bone μ-map compared to Dixon μ-map were observed in all subjects. Figure 2A, illustrates the typical line profile (Figure 2B) through the center slice of the brain of one subject when the PETdx and PETdxbone images are compared. An increase of 5% in the center to nearly 20% in the areas around the cortex was seen. Group mean values and standard deviation of absolute and relative activity in PETdx and PETdxbone for whole brain and regions of interest are listed in Table 2. Paired student t-test showed that the absolute mean activity in PET data reconstructed with Dixon+bone μ-map was statistically higher than PET reconstruction with Dixon μ-map in all thirteen ROIs. Conversely, lower relative activity was found in rPETdxbone when compared to rPETdx in all ROIs except, occipital and cerebellum. These regional differences are depicted in Figure 3A as the percent relative difference in mean activity between PETdx and PETdxbone, and in Figure 3B as the percent relative difference in rPETdx and rPETdxbone for all thirteen ROIs.
Figure 1. Illustration of a Dixon+bone μ-map generated in a representative subject. Images (A) T1-weighted MPRAGE, (B) bone probability map, (C) Dixon μ-map, and (D) Dixon+bone μ-map are shown in axial, coronal and sagittal views of a single slice.
Figure 2. Line profile (A) across the center slice of a map (B) of relative difference (% RD) between PET images reconstructed with Dixon and Dixon+bone μ-maps. Images are presented from a representative subject.
Table 2. Group mean and standard deviation (std) of absolute and relative PETdx and PETdxbone activity concentration (kBq/ml) in regions of interest across the brain.
Figure 3. Regions of interest group means percent relative difference (% RD) in (A) mean activity concentration and (B) mean relative activity between PET signals reconstructed with Dixon and Dixon+bone μ-maps. Errors bars represent standard deviation on the means. Regions abbreviations: Fro, Frontal; Cng, Cingulate; Ins, Insula; Par, Parietal; Tem, Temporal; Hpp, Hippocampus; Amg, Amygdala; Tha, Thalamus; Cad, Caudate; Pal, Globus Pallidus; Put, Putamen; Occ, Occipital; Ceb, Cerebellum.
Correlations of Perfusion to Glucose Uptake
Images of relative CBF, PETdx, and PETdxbone averaged across all subjects are displayed in Figure 4. In general, these group-wise images are similar in appearance. On visual inspection, the group-averaged PET data reconstructed with Dixon+bone μ-map had an apparent higher intensity compared to PET data reconstructed with Dixon μ-map, particularly around the cortex. An additional difference between the data sets was the aliasing artifact observed in the CBF images, which are most visible in distal images. The mean group gray matter CBF was 43.82 ± 4.46 ml/100 g/min. The mean tSNR of whole brain pCASL signal in gray matter across all subjects was 2.83 ± 1.18.
Figure 4. Whole-brain group maps of relative perfusion and relative glucose uptake from; left to right: pCASL-CBF, 18F-FDG-PET reconstructed with Dixon+bone MRAC and 18F-FDG-PET reconstructed with Dixon MRAC. Spatial blurring seen in the bottom row images of the perfusion maps could be minimized with segmented multi-shot 3D GRASE acquisitions.
The average voxel-by-voxel correlation coefficient across whole-brain gray matter was r = 0.63 (p < 0.001) between rCBF and rPETdx, and r = 0.53 (p < 0.001) between rCBF and rPETdxbone. The correlation coefficient for each ROI, derived from voxel-wise comparison between rCBF and rPETdx and between rCBF and rPETdxbone, are presented in Table 3 along with corresponding total number of voxels. For the comparison between rCBF and rPETdx, moderate-to-high correlations were observed in all ROIs with the highest in the caudate and the lowest in the deep-lying structures of the limbic system such as the hippocampus and amygdala. Similar trends were observed between rCBF and rPETdxbone. Results of voxel-by-voxel independent samples t-test listed on Table 4, showed areas of significant increase or decrease in relative perfusion compared to relative 18F-FDG-PET signals reconstructed with Dixon and with Dixon+bone μ-maps for clusters that met the set statistical threshold. Maps of statistical differences between rCBF and either rPETdx or PETdxbone overlaid on axial slices of a single subject T1-weighted image are displayed in Figure 5. Regions with higher relative perfusion are marked in red while regions with higher relative glucose uptake are marked in blue.
Table 3. Results of voxel-by-voxel Pearson product moment correlation within thirteen regions of interest (p < 0.001).
Table 4. Results of independent samples t-test in gray matter comparing rCBF to rPETdx and comparing rCBF to rPETdxbone for clusters >50 voxels (p < 0.05, FDR). MNI coordinates (x–z) and corresponding anatomical location are included.
Figure 5. Whole brain voxel-by-voxel comparison between relative cerebral blood flow and relative 18F-FDG-PET activity reconstructed with Dixon (A) and Dixon+bone (B) μ-maps. Areas with greater relative perfusion are shown in red while areas with greater relative glucose uptake are shown in blue. Statistical differences are set at threshold for cluster >50 voxels, p < 0.05 corrected for multiple comparisons with FDR.
The feasibility of an integrated whole-body PET-MRI system in simultaneous acquisition of perfusion MRI and metabolic PET neuroimaging with an enhanced attenuation correction method was investigated in this study. Following efforts by other groups to minimize errors in PET signal attenuation corrected with MRI (Kops et al., 2009; Catana et al., 2010; Izquierdo-Garcia et al., 2014; Poynton et al., 2014), and given the lack of consensus on acceptable MRAC methods for bone attenuation (Bailey et al., 2013), we explored a practical workaround for MRAC prior to the start of a neurodegenerative study with PET-MRI.
Direct segmentation of high resolution T1-weighted MR data into various tissue classes for MR attenuation correction of PET have been proposed by Zaidi et al. (2003) and Wagenknecht et al. (2009). Segmentation of voxel intensities were performed with sophisticated algorithms that use fuzzy clustering or a priori knowledge of tissue location and shape to derive MRAC μ-maps with ~6% maximum absolute relative difference compared to ground truth (Zaidi et al., 2003; Wagenknecht et al., 2011). However, these methods are time consuming, require post-processing algorithms that are not readily available, and in Zaidi et al.'s (2003) method; require user intervention to improve bone segmentation. In addition, these proposed approaches are sensitive to MR intensity inhomogeneity errors which, when left unaccounted for can exacerbate inherent MR signal issues such as partial volume effects and motion artifacts, potentially inflating bias in MRAC μ-maps. In this study, bone segmentation of T1-weighted data was easily derived with SPM, an automated pipeline routinely used for brain segmentation in neuroimaging. For each voxel, the unified segmentation method in SPM determines intensity distributions of tissue types using a mixed model cluster analysis and spatial priors from a T1-weighted atlas derived from large number of subjects (Ashburner and Friston, 2005). Unlike aforementioned direct segmentation methods, the unified method accounts for non-uniformities in intensity distributions, improving segmentation accuracy (Tsang et al., 2008).
In general, a ~10% increase in whole brain mean 18F-FDG activity corrected with Dixon+bone compared to standard Dixon MRAC was observed. A gradual increase in mean 18F-FDG activity from ~5% in the center to ~20% in the cortex was observed in PET data corrected with Dixon+bone compared to the same PET data corrected with just Dixon alone, as illustrated in Figure 2 for one individual. This finding was consistent with a recent study comparing PET images corrected with the standard Dixon method to images corrected with Dixon plus bone information derived from individual CT (Andersen et al., 2014). In this study, Andersen et al. (2014) reported a 15% difference in mean activity from the center to the edge of the brain in PET signal reconstructed with Dixon compared to Dixon plus CT bone. Errors caused by not including bone information are more pronounced at the edge than the center of the brain as photons from the edge of the brain travel longer average path lengths through skull than photons from the center (Andersen et al., 2014).
Regional variability in mean activity between PETdx and PETdxbone was also observed in the ROI analysis. Cortical structures such as frontal and temporal lobe had significantly higher relative differences than central brain regions such as basal ganglia and limbic systems (Figure 3A). All thirteen brain regions investigated had statistical significant lower mean activity when reconstructed with the standard Dixon compared to the reconstruction with the Dixon plus bone (Table 2). Similar regional variability in PET-corrected with standard Dixon MRAC compared to CT attenuation maps have recently been observed by other groups (Dickson et al., 2014; Hitz et al., 2014). Scaling the PET images by their global mean removed the overall underestimation but resulted in regional over- and under-estimations (Figure 3B). These regional trends matched results from a previous study (Dickson et al., 2014) and indicate that irrespective of the chosen reference region, cerebellum (Dickson et al., 2014) or global mean as used in the current study, relative PET measurements reconstructed with Dixon still possess MRAC-related bias, albeit slightly reduced. Although signal normalization to a reference value can reduce systemic errors in PET reconstruction, it can inflate regional values particularly in cases where group differences actually exist between the reference values (Borghammer et al., 2008). This might explain the apparent increase in relative PET activity with standard Dixon compared to Dixon+bone (Table 2) or UTE (Dickson et al., 2014) in nearly all regions of the brain. Altogether, the results are in line with previous studies showing that the addition of bone information, either from MRI (UTE/MPRAGE) or CT, to the MRAC maps removes the radial bias (Andersen et al., 2014) and reduces differences in PET signals from ~20 to 10% when compared to the gold standard, CTAC (Kops et al., 2009; Catana et al., 2010; Marshall et al., 2013; Dickson et al., 2014).
The clinical significance of incorporating bone attenuation in PET reconstruction was highlighted by correlating relative cerebral glucose uptake acquired with 18F-FDG-PET to relative cerebral perfusion acquired with pCASL across regions in the brain. Recent evidence suggest that a pattern of regional distribution of perfusion to glucose metabolism exist in healthy brains (Newberg et al., 2005; Vaishnavi et al., 2010; Cha et al., 2013) and disruptions to regional blood flow-glucose metabolism coupling are potential markers of brain dysfunction (Vlassenko et al., 2010). In general, a good correlation between perfusion and glucose uptake was observed in this study over the entire brain. Regional variability in relative correlations was found, with the highest correlations in the occipital and caudate, and the lowest in the limbic structures regardless of the MRAC method used (Table 3). However, associations between perfusion and PET signals reconstructed with Dixon+bone closely resemble patterns of regional relative perfusion to glucose metabolism reported elsewhere, with 18F-FDG-PET and PET perfusion tracers (Bentourkia et al., 2000; Gur et al., 2009; Musiek et al., 2012), and with 18F-FDG-PET and ASL using standalone PET and MRI scanners (Newberg et al., 2005; Chen et al., 2011; Cha et al., 2013).
When compared to Cha et al. (2013), where similar regions were explored with a similar pCASL labeling scheme, higher regional correlations of rCBF to rPETdxbone were observed in the current study, except in the insula and putamen (Table 3). This apparent increase (~30%) in regional coupling of perfusion to glucose metabolism is likely a result of the improved registration accuracy provided by simultaneous imaging and the improved tSNR provided by using a 3D pCASL method (Günther et al., 2005; Vidorreta et al., 2013). The temporal SNR of 3D-GRASE pCASL was not compromised by the use of a PET-compatible head coil, and the mean tSNR reported here was in line with values reported elsewhere (Günther et al., 2005; Vidorreta et al., 2013).
Results from voxel-by-voxel t-tests between rCBF and rPET demonstrated regions of hyper- or hypo-perfusion not matched to glucose metabolism (Figure 5 and Table 4). These findings, specifically from PET data corrected with Dixon+bone MRAC, are in line with previous findings of increased resting perfusion to resting metabolism in areas of the brain associated with the default mode network (Gur et al., 2009; Cha et al., 2013), suggesting that these brain regions require increased blood flow due to the sustained state of arousal (Vaishnavi et al., 2010). Apparent bias in bone attenuation from standard Dixon MRAC was evident in the comparison between relative perfusion and glucose uptake as demonstrated in Figure 5 and in comparison to previous studies (Newberg et al., 2005; Chen et al., 2011; Cha et al., 2013). In particular, unexpected increases in rCBF were observed in the left paracentral, left inferior parietal and right superior temporal gyrus when compared to PETdx. The spatial pattern of these differences, notably their occurrence closer to the edge of brain, is suggestive of attenuation errors. Indeed, these regions of elevated rCBF were either removed or significantly reduced when compared to PETdxbone.
It is possible that segmentation errors such as voxel misclassifications in Dixon+bone μ-maps could account for some differences in correlations between the current study and observations reported elsewhere. Differences in μa for bone assigned to MR or CT derived μ-maps, and differences in PET scanner geometry and reconstruction parameters including number of iterations, subsets and smoothing kernels could also influence final observations. In addition, limitations inherent to pCASL and to PET-FDG imaging such as partial voluming, differences in cellular localization of blood water and FDG tracers, and optimal selection of post-labeling delay in ASL with minimal compromise to SNR could affect accuracy of regional comparisons (cf. Cha et al., 2013 for detailed description). The evaluation of Dixon+bone μ-maps would have benefited from comparison to ground truth using CT data, which was not feasible in this study. The PET-CT whole body survey acquired on all subjects prior to PET-MRI imaging, permitted coverage solely from the level of the eyes to the thighs as per clinical protocol, limiting whole brain comparisons. As such, absolute and regional performance of the enhanced Dixon μ-maps was not evaluated. Over estimation of bone information was observed in the inferior regions of the brain and can be seen in the Dixon+bone μ-maps (Figure 1 and Supplementary Figures) and in the comparisons of PETdxbone to pCASL CBF (Figure 5). Regions close to the base of the skull where bone anatomy is thin can yield low MR bone signals with MPRAGE (Wagenknecht et al., 2011) or UTE (Dickson et al., 2014) sequences, and these regions are prone to signal misclassification errors in soft tissue/bone/air interfaces in areas around the sinuses- an issue that affects most atlas-based MR-derived attenuation correction methods. Efforts were made to minimize bone misclassification by morphological filtering of non-bone voxels. A recent study demonstrated that SPM8-based method for generating MRAC is robust with a reported accuracy of under 4% when compared to CTAC (Izquierdo-Garcia et al., 2014). The results presented here were in good agreement with previous studies, suggesting that these issues had little or no effect on the study outcome. The regional correlation of pCASL CBF to 18F-FDG-PET demonstrated here for the first time with simultaneous imaging can be further improved with imaging of concurrent uptake of the pCASL and FDG tracer.
In general, the agreement between this study and previous studies using CT attenuation maps indicate that Dixon attenuation maps enhanced with bone information from high resolution T1-weighted images can provide a feasible method for correction of PET signal attenuation in PET-MRI neuroimaging. Errors in MRI-derived attenuation correction maps including tissue misclassification can be further minimized by methods that employ tissue classifiers from CT atlas to better guide segmentation of UTE or MPRAGE data (Poynton et al., 2014).
Conceived and designed the experiments: Udunna C. Anazodo, Keith S. St. Lawrence. Performed the experiments: Udunna C. Anazodo, John Butler. Analyzed the data: Udunna C. Anazodo, Jonathan D. Thiessen, Tracy Ssali, Jonathan Mandel, William Pavlosky, Keith S. St. Lawrence. Interpreted data: Udunna C. Anazodo, Jonathan Mandel, William Pavlosky, Keith S. St. Lawrence. Contributed reagents/materials/analysis tools: Matthias Günther, Frank S. Prato, R. Terry Thompson. Wrote the paper: Udunna C. Anazodo, Keith S. St. Lawrence, Frank S. Prato, R. Terry Thompson.
Conflict of Interest Statement
Dr. M. Gunther receives research and development funds from Siemens Healthcare. For the remaining authors none were declared. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank Dr. Gerry Moran of Siemens Canada for his technical advice and Dr. Elizabeth Finger for her contributions with the dementia data. This work was supported by a personnel award to K. St. Lawrence from the Heart and Stroke Foundation, Ontario Provincial Office, and funded in part by the Canadian Institutes of Health Research, Ontario Research Fund, the Canadian Foundation for Innovation, Lawson Health Research Institute and St. Joseph's Health Care, London.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fnins.2014.00434/abstract
Andersen, F. L., Ladefoged, C. N., Beyer, T., Keller, S. H., Hansen, A. E., Højgaard, L., et al. (2014). Combined PET/MR imaging in neurology: MR-based attenuation correction implies a strong spatial bias when ignoring bone. Neuroimage 84, 206–216. doi: 10.1016/j.neuroimage.2013.08.042
Bailey, D. L., Barthel, H., Beyer, T., Boellaard, R., Gückel, B., Hellwig, D., et al. (2013). Summary report of the first international workshop on PET/MR imaging, March 19-23, 2012, Tübingen, Germany. Mol. Imaging Biol. 15, 361–371. doi: 10.1007/s11307-013-0623-1
Bentourkia, M., Bol, A., Ivanoiu, A., Labar, D., Sibomana, M., Coppens, A., et al. (2000). Comparison of regional cerebral blood flow and glucose metabolism in the normal brain: effect of aging. J. Neurol. Sci. 181, 19–28. doi: 10.1016/S0022-510X(00)00396-8
Beyer, T., Weigert, M., Quick, H. H., Pietrzyk, U., Vogt, F., Palm, C., et al. (2008). MR-based attenuation correction for torso-PET/MR imaging: pitfalls in mapping MR to CT data. Eur. J. Nucl. Med. Mol. Imaging 35, 1142–1146. doi: 10.1007/s00259-008-0734-0
Borghammer, P., Jonsdottir, K. Y., Cumming, P., Ostergaard, K., Vang, K., Ashkanian, M., et al. (2008). Normalization in PET group comparison studies—the importance of a valid reference region. Neuroimage 40, 529–540. doi: 10.1016/j.neuroimage.2007.12.057
Catana, C., Benner, T., van der Kouwe, A., Byars, L., Hamm, M., Chonde, D. B., et al. (2011). MRI-assisted PET motion correction for neurologic studies in an integrated MR-PET scanner. J. Nucl. Med. 52, 154–161. doi: 10.2967/jnumed.110.079343
Catana, C., van der Kouwe, A., Benner, T., Michel, C. J., Hamm, M., Fenchel, M., et al. (2010). Toward implementing an MRI-based PET attenuation-correction method for neurologic studies on the MR-PET brain prototype. J. Nucl. Med. 51, 1431–1438. doi: 10.2967/jnumed.109.069112
Cha, Y.-H. K., Jog, M. A., Kim, Y.-C., Chakrapani, S., Kraman, S. M., and Wang, D. J. (2013). Regional correlation between resting state FDG PET and pCASL perfusion MRI. J. Cereb. Blood Flow Metab. 33, 1909–1914. doi: 10.1038/jcbfm.2013.147
Chen, Y., Wolk, D. A., Reddin, J. S., Korczykowski, M., Martinez, P. M., Musiek, E. S., et al. (2011). Voxel-level comparison of arterial spin-labeled perfusion MRI and FDG-PET in Alzheimer disease. Neurology 77, 1977–1985. doi: 10.1212/WNL.0b013e31823a0ef7
Delso, G., Fürst, S., Jakoby, B., Ladebeck, R., Ganter, C., Nekolla, S. G., et al. (2011). Performance measurements of the Siemens mMR integrated whole-body PET/MR scanner. J. Nucl. Med. 52, 1914–1922. doi: 10.2967/jnumed.111.092726
Dickson, J. C., O'Meara, C., and Barnes, A. (2014). A comparison of CT- and MR-based attenuation correction in neurological PET. Eur. J. Nucl. Med. Mol. Imaging. 41, 1176–1189. doi: 10.1007/s00259-013-2652-z
Drzezga, A., Barthel, H., Minoshima, S., and Sabri, O. (2014). Potential clinical applications of PET/MR imaging in neurodegenerative diseases. J. Nucl. Med. 55, 47S–55S. doi: 10.2967/jnumed.113.129254
Dukart, J., Mueller, K., Horstmann, A., Barthel, H., Möller, H. E., Villringer, A., et al. (2011). Combined evaluation of FDG-PET and MRI improves detection and differentiation of dementia. PLoS ONE 6:e18111. doi: 10.1371/journal.pone.0018111
Gur, R. C., Ragland, J. D., Reivich, M., Greenberg, J. H., Alavi, A., and Gur, R. E. (2009). Regional differences in the coupling between resting cerebral blood flow and metabolism may indicate action preparedness as a default state. Cereb. Cortex 19, 375–382. doi: 10.1093/cercor/bhn087
Gutierrez, D., Montandon, M.-L., Assal, F., Allaoua, M., Ratib, O., Lövblad, K.-O., et al. (2012). Anatomically guided voxel-based partial volume effect correction in brain PET: impact of MRI segmentation. Comput. Med. Imaging Graph. 36, 610–619. doi: 10.1016/j.compmedimag.2012.09.001
Hitz, S., Habekost, C., Fürst, S., Delso, G., Förster, S., Ziegler, S., et al. (2014). Systematic comparison of the performance of integrated whole-body PET/MR imaging to conventional PET/CT for 18F-FDG brain imaging in patients examined for suspected dementia. J. Nucl. Med. 55, 923–931. doi: 10.2967/jnumed.113.126813
Hofmann, M., Steinke, F., Scheel, V., Charpiat, G., Farquhar, J., Aschoff, P., et al. (2008). MRI-based attenuation correction for PET/MRI: a novel approach combining pattern recognition and atlas registration. J. Nucl. Med. 49, 1875–1883. doi: 10.2967/jnumed.107.049353
Izquierdo-Garcia, D., Hansen, A. E., Förster, S., Benoit, D., Schachoff, S., Fürst, S., et al. (2014). An SPM8-based approach for attenuation correction combining segmentation and nonrigid template formation: application to simultaneous PET/MR brain imaging. J. Nucl. Med. 55, 1825–1830. doi: 10.2967/jnumed.113.136341
Kops, E. R., Wagenknecht, G., Scheins, J., Tellmann, L., and Herzog, H. (2009). “Attenuation correction in MR-PET scanners with segmented T1-weighted MR images,” in 2009 IEEE Nuclear Science Symposium Conference Record (Orlando, FL), 2530–2533.
Maldjian, J. A., Laurienti, P. J., Kraft, R. A., and Burdette, J. H. (2003). An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets. Neuroimage 19, 1233–1239. doi: 10.1016/S1053-8119(03)00169-1
Marshall, H. R., Patrick, J., Laidley, D., Prato, F. S., Butler, J., Théberge, J., et al. (2013). Description and assessment of a registration-based approach to include bones for attenuation correction of whole-body PET/MRI. Med. Phys. 40, 082509. doi: 10.1118/1.4816301
Martinez-Möller, A., Souvatzoglou, M., Delso, G., Bundschuh, R. A., Chefd'hotel, C., Ziegler, S. I., et al. (2009). Tissue classification as a potential approach for attenuation correction in whole-body PET/MRI: evaluation with PET/CT data. J. Nucl. Med. 50, 520–526. doi: 10.2967/jnumed.108.054726
Morbelli, S., Piccardo, A., Villavecchia, G., Dessi, B., Brugnolo, A., Piccini, A., et al. (2010). Mapping brain morphological and functional conversion patterns in amnestic MCI: a voxel-based MRI and FDG-PET study. Eur. J. Nucl. Med. Mol. Imaging 37, 36–45. doi: 10.1007/s00259-009-1218-6
Musiek, E. S., Chen, Y., Korczykowski, M., Saboury, B., Martinez, P. M., Reddin, J. S., et al. (2012). Direct comparison of fluorodeoxyglucose positron emission tomography and arterial spin labeling magnetic resonance imaging in Alzheimer's disease. Alzheimers. Dement. 8, 51–59. doi: 10.1016/j.jalz.2011.06.003
Newberg, A. B., Wang, J., Rao, H., Swanson, R. L., Wintering, N., Karp, J. S., et al. (2005). Concurrent CBF and CMRGlc changes during human brain activation by combined fMRI-PET scanning. Neuroimage 28, 500–506. doi: 10.1016/j.neuroimage.2005.06.040
Pappas, I. P. I., Styner, M., Malik, P., Remonda, L., and Caversaccio, M. (2005). Automatic method to assess local CT-MR imaging registration accuracy on images of the head. AJNR Am. J. Neuroradiol. 26, 137–44.
Pichler, B. J., Kolb, A., Nägele, T., and Schlemmer, H.-P. (2010). PET/MRI: paving the way for the next generation of clinical multimodality imaging applications. J. Nucl. Med. 51, 333–336. doi: 10.2967/jnumed.109.061853
Poynton, C. B., Chen, K. T., Chonde, D. B., Izquierdo-Garcia, D., Gollub, R. L., Gerstner, E. R., et al. (2014). Probabilistic atlas-based segmentation of combined T1-weighted and DUTE MRI for calculation of head attenuation maps in integrated PET/MRI scanners. Am. J. Nucl. Med. Mol. Imaging 4, 160–171.
Tsang, O., Gholipour, A., Kehtarnavaz, N., Gopinath, K., Briggs, R., and Panahi, I. (2008). Comparison of tissue segmentation algorithms in neuroimage analysis software tools. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2008, 3924–3928. doi: 10.1109/IEMBS.2008.4650068
Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15, 273–289. doi: 10.1006/nimg.2001.0978
Vaishnavi, S. N., Vlassenko, A. G., Rundle, M. M., Snyder, A. Z., Mintun, M. A., and Raichle, M. E. (2010). Regional aerobic glycolysis in the human brain. Proc. Natl. Acad. Sci. U.S.A. 107, 17757–17762. doi: 10.1073/pnas.1010459107
Vidorreta, M., Wang, Z., Rodríguez, I., Pastor, M. A., Detre, J. A., and Fernández-Seara, M. A. (2013). Comparison of 2D and 3D single-shot ASL perfusion fMRI sequences. Neuroimage 66, 662–671. doi: 10.1016/j.neuroimage.2012.10.087
Vlassenko, A. G., Vaishnavi, S. N., Couture, L., Sacco, D., Shannon, B. J., Mach, R. H., et al. (2010). Spatial correlation between brain aerobic glycolysis and amyloid-β (Aβ) deposition. Proc. Natl. Acad. Sci. U.S.A. 107, 17763–17767. doi: 10.1073/pnas.1010461107
Wagenknecht, G., Kops, E. R., Mantlik, F., Fried, E., Pilz, T., Hautzel, H., et al. (2011). “Attenuation correction in MR-BrainPET with segmented Tl-weighted MR images of the patient' s head—a comparative study with CT,” in 2011 IEEE Nuclear Science Symposium Conference Record (Valencia: IEEE), 2261–2266.
Wagenknecht, G., Kops, E. R., Tellmann, L., and Herzog, H. (2009). “Knowledge-based segmentation of attenuation- relevant regions of the head in T1-weighted MR images for attenuation correction in MR/PET systems,” in 2009 IEEE Nuclear Science Symposium Conference Record (Orlando, FL: IEEE), 3338–3343.
Wang, J., Alsop, D. C., Song, H. K., Maldjian, J. A., Tang, K., Salvucci, A. E., et al. (2003). Arterial transit time imaging with flow encoding arterial spin tagging (FEAST). Magn. Reson. Med. 50, 599–607. doi: 10.1002/mrm.10559
Keywords: PET-MRI, attenuation correction, arterial spin labeling, ASL, 18F-fluorodexoyglucose, FDG, cerebral blood flow, glucose uptake
Citation: Anazodo UC, Thiessen JD, Ssali T, Mandel J, Günther M, Butler J, Pavlosky W, Prato FS, Thompson RT and St. Lawrence KS (2015) Feasibility of simultaneous whole-brain imaging on an integrated PET-MRI system using an enhanced 2-point Dixon attenuation correction method. Front. Neurosci. 8:434. doi: 10.3389/fnins.2014.00434
Received: 28 August 2014; Accepted: 10 December 2014;
Published online: 05 January 2015.
Edited by:Jeanette Mumford, University of Texas at Austin, USA
Reviewed by:Valentina Garibotto, Geneva University and Geneva University Hospital, Switzerland
Holger Schmidt, University Hospital Tuebingen, Germany
Flemming Littrup Andersen, University of Copenhagen, Denmark
Copyright © 2015 Anazodo, Thiessen, Ssali, Mandel, Günther, Butler, Pavlosky, Prato, Thompson and St. Lawrence. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Udunna C. Anazodo, Lawson Health Research Institute, 268 Grosvenor St., London, ON N6A4V2, Canada e-mail: firstname.lastname@example.org