Original Research ARTICLE
Quantitative multi-parameter mapping of R1, PD*, MT, and R2* at 3T: a multi-center validation
- 1Wellcome Trust Centre for Neuroimaging, UCL Institute of Neurology, University College London, London, UK
- 2Department of Psychiatry, University of Cambridge, Cambridge, UK
- 3Behavioural and Clinical Neuroscience Institute, University of Cambridge, Cambridge, UK
- 4Cambridgeshire and Peterborough NHS Foundation Trust, Cambridge, UK
- 5Department of Clinical Neuroscience, Wolfson Brain Imaging Centre, University of Cambridge, Cambridge, UK
- 6MRC Cognition and Brain Sciences Unit, Cambridge, UK
- 7GlaxoSmithKline, Clinical Unit Cambridge, Addenbrooke's Hospital, Cambridge, UK
- 8Laboratoire de recherche en neuroimagerie, Département des neurosciences cliniques, CHUV, University of Lausanne, Lausanne, Switzerland
Multi-center studies using magnetic resonance imaging facilitate studying small effect sizes, global population variance and rare diseases. The reliability and sensitivity of these multi-center studies crucially depend on the comparability of the data generated at different sites and time points. The level of inter-site comparability is still controversial for conventional anatomical T1-weighted MRI data. Quantitative multi-parameter mapping (MPM) was designed to provide MR parameter measures that are comparable across sites and time points, i.e., 1 mm high-resolution maps of the longitudinal relaxation rate (R1 = 1/T1), effective proton density (PD*), magnetization transfer saturation (MT) and effective transverse relaxation rate (R2* = 1/T2*). MPM was validated at 3T for use in multi-center studies by scanning five volunteers at three different sites. We determined the inter-site bias, inter-site and intra-site coefficient of variation (CoV) for typical morphometric measures [i.e., gray matter (GM) probability maps used in voxel-based morphometry] and the four quantitative parameters. The inter-site bias and CoV were smaller than 3.1 and 8%, respectively, except for the inter-site CoV of R2* (<20%). The GM probability maps based on the MT parameter maps had a 14% higher inter-site reproducibility than maps based on conventional T1-weighted images. The low inter-site bias and variance in the parameters and derived GM probability maps confirm the high comparability of the quantitative maps across sites and time points. The reliability, short acquisition time, high resolution and the detailed insights into the brain microstructure provided by MPM makes it an efficient tool for multi-center imaging studies.
Multi-center studies using magnetic resonance imaging (MRI) facilitate the detection of small effects, detailed estimation of neuroanatomical population variance and investigation of rare diseases. For example, recent multi-center studies identified reliable markers for Alzheimer's disease (Kloppel et al., 2008a,b), distributed anatomical differences in autism spectrum disorder (Ecker, 2012) or the relationship of inter-individual differences in character traits to anatomical differences (Schilling et al., 2013), which would have been difficult or impossible in a single-site setting.
However, the reliability and sensitivity of multi-center studies crucially depend on the comparability of structural MRI data generated at the different sites (Tofts and Collins, 2012). The level of inter-site comparability is still controversial for conventional MRI data. Some studies demonstrated systematic inter-site differences in structural T1-weighted images, which biased morphometric analyses (Focke et al., 2011), whereas other studies argued that typical pathology-related differences can be detected reliably by, for example, adjusting for potential inter-site differences (Pardoe et al., 2008; Stonnington et al., 2008; Suckling et al., 2012).
Quantitative anatomical MRI (qMRI) aims to overcome the inter-site bias issue, since it is specifically designed to provide absolute measures and thus data that are comparable across sites and time points (Tofts, 2003). Although various anatomical qMRI methods were developed (Tofts, 2003), there are only few studies that validated them in multiple centers. Deoni et al. (2008) validated the use of quantitative mapping of the longitudinal and the transverse relaxation time (T1 and T2) in a multi-center study at 1.5T. They demonstrated a high comparability between sites and reproducibility within a single site in scan-rescan experiments (<10% deviation). A magnetization transfer ratio (MTR) imaging protocol was optimized and validated in a multi-center study on multiple sclerosis at 1.5T (Barker et al., 2005; Ropele et al., 2005; Tofts et al., 2006). Careful alignment of imaging protocols and post-processing achieved a high comparability between sites [<4% deviation, (Ropele et al., 2005)]. The transverse relaxation time (T2) was quantified in a multi-center study on Alzheimer's disease and compared across sites at 1.5T (Bauer et al., 2010). Significant inter-site bias of up to 20% was observed, obscuring pathological changes.
These multi-center studies were performed at rather low resolutions (1.2–8.7 mm3), a field strength of 1.5T and investigated only few quantitative parameters. Spatial coverage was usually small and often did not allow for whole-brain imaging.
Recently, a comprehensive quantitative multi-parameter mapping (MPM) approach was developed at 3T (Helms et al., 2008a, 2009; Weiskopf et al., 2011), which provides high resolution maps of the longitudinal relaxation rate (R1 = 1/T1), effective proton density (PD*), magnetization transfer saturation (MT) and effective transverse relaxation rate (R2* = 1/T2*). Whole-brain maps are acquired with 1 mm3 isotropic resolution in a clinically feasible time of approximately 24 mins. The multiple parameter maps and the high resolution allow for a detailed assessment of the white matter (WM) and gray matter (GM) tissue microstructure (Draganski et al., 2011; Dick et al., 2012; Sereno et al., 2012). MPM was used to study a wide range of different tissue changes in e.g., healthy aging (Draganski et al., 2011) or prosopagnosia (Garrido et al., 2009). The high sensitivity and specificity of the approach also improves segmentation of subcortical structures in morphometric studies (Helms et al., 2009). An advanced MPM approach using a higher isotropic resolution of 800 μm even allowed for mapping of cortical myelination and parcellation of brain areas (Dick et al., 2012; Sereno et al., 2012).
We validated the MPM approach at 3T for use in multi-center studies. The inter-site variance across three sites and intra-site variance were determined for typical morphometric measures (i.e., GM probability maps) and the four quantitative parameters. We also compared MPM to conventional T1-weighted (T1w) imaging using FLASH (fast low angle shot) imaging.
Participants and Centers
The same five healthy volunteers [2 males, age 24.2 ± 1.6 yrs (mean ± SD)] were scanned at these three sites within 12 weeks: (1) Wellcome Trust Centre for Neuroimaging, London; (2) Medical Research Council Cognition and Brain Sciences Unit and (3) Wolfson Brain Imaging Centre, Cambridge (in the following called WTCN, MRC CBSU and WBIC, respectively). The study received ethical approval by the Cambridge Psychology Research Ethics Committee (ref: 2012.17) and all scanning sites obtained local approvals and written informed consent was obtained before scanning. This study was one of several pilot studies conducted to demonstrate the feasibility of the imaging approach used by the NeuroScience in Psychiatry Network (NSPN), which addresses how psychiatric disorders are related to abnormal maturation of brain systems.
All scans were acquired on 3T whole body MRI systems (Magnetom TIM Trio, Siemens Healthcare, Erlangen, Germany; VB17 software version) operated with the standard 32-channel radio-frequency (RF) receive head coil and RF body coil for transmission. The MPM comprised three multi-echo 3D FLASH (fast low angle shot) scans, one RF transmit field map and one static magnetic (B0) field map scan (Weiskopf et al., 2011). The MPM acquisition and post-processing were developed and optimized in previous studies and are briefly described here for convenience (Helms et al., 2008a,b, 2009, 2011; Lutti et al., 2010, 2012; Weiskopf et al., 2011).
Three different multi-echo FLASH scans were acquired with predominant T1-, PD-, and MT-weighting by appropriate choice of the repetition time (TR) and the flip angle α: TR/α = 18.7 ms/20° for the T1w scan and 23.7 ms/6° for the PDw and the MTw scans. MT-weighting was achieved by applying an off-resonance Gaussian-shaped RF pulse (4 ms duration, 220° nominal flip angle, 2 kHz frequency offset from water resonance) prior to the excitation. Multiple gradient echoes were acquired with alternating readout polarity at six equidistant echo times (TE) between 2.2 and 14.7 ms for the T1w and MTw acquisitions and at 8 equidistant TE between 2.2 ms and 19.7 ms for the PDw acquisition. Other acquisition parameters were: 1 mm isotropic resolution, 176 sagittal partitions, field of view (FOV) = 256 × 240 mm, matrix = 256 × 240 × 176, parallel imaging using GRAPPA factor 2 in phase-encoding (PE) direction, 6/8 partial Fourier in partition direction, non-selective RF excitation, readout bandwidth BW = 425 Hz/pixel, RF spoiling phase increment = 50°, total acquisition time ~19 min.
The design of the protocol took into account the following criteria. The primary goal was to acquire all three FLASH whole brain images with 1 mm resolution within ca. 20 min. The short acquisition time was achieved by combining GRAPPA parallel imaging with Partial Fourier acquisition. The moderate 2× GRAPPA acceleration avoided deterioration of image quality due to a poor geometry factor (Pruessmann et al., 1999). The echo train length was limited to ca. 20 ms to trade off emerging R2* contrast against susceptibility-induced signal losses, while it allowed for determining R2* and averaging of images for high SNR (Helms and Dechent, 2009). It also allowed for a high readout bandwidth (425 Hz/pixel) to minimize off-resonance and chemical shift artifacts. To achieve the same TR for the MTw and PDw scans, only 6 echoes were acquired for the MTw scan, to accommodate the 4 ms long off-resonance RF pulse. The 2 kHz frequency offset of the MT saturation pulse was chosen so that direct saturation effects were reduced and stability was increased (Helms et al., 2008b). The flip angle of the off-resonance pulse was titrated to keep the specific absorption rate (SAR) below ca. 75% of the normal mode SAR limit. The T1w scan was acquired with a shorter TR and fewer echoes to increase the signal-to-noise ratio (SNR) per time unit. The flip angles of the acquisitions were optimized using a semi-empirical approach, in order to maximize SNR while limiting bias due to imperfect RF spoiling (Yarnykh, 2010; Helms et al., 2011).
Maps of the local RF transmit field were measured and estimated from a 3D EPI acquisition of spin and stimulated echoes (SE and STE) with different refocusing flip angles (Lutti et al., 2010, 2012). Imaging parameters were: 4 mm isotropic resolution, matrix = 64 × 48 × 48 and FOV = 256 mm × 192 mm × 192 mm along readout × PE × partition direction, parallel imaging using GRAPPA factor 2 × 2 in PE and partition direction, TESE/TESTE/TM (mixing time)/TR = 37.06/37.06/31.2/500 ms, acquisition time 3 min. The flip angles of the SE/STE refocusing pulses were decreased from 230°/115° to 130°/65° and in steps of 10°/5°.
To correct the 3D EPI RF transmit field maps for geometric distortion and off-resonance effects, a map of the static magnetic field (B0) was acquired with the following parameters (Lutti et al., 2010, 2012): 2D double-echo FLASH sequence with 64 axial slices, slice thickness = 2 mm, inter-slice gap = 1 mm, TR = 1020 ms, TE1/TE2 = 10/12.46 ms, α = 90°, matrix = 64 × 64, FOV = 192 × 192 mm, left-right PE direction, BW = 260 Hz/pixel, flow compensation, acquisition time ~2 min.
Estimation of Parameter Maps
All data analyses and processing were performed in Matlab (The MathWorks Inc., Natick, MA, USA) using SPM8 (www.fil.ion.ucl.ac.uk/spm) and custom-made Matlab tools. Arithmetic mean T1w, PDw and MTw images were calculated from the first 6 multi-echo acquisitions (with the shortest TE), in order to increase the SNR. The resulting three mean images were used to calculate the parameter maps of the MT saturation, the apparent longitudinal relaxation rate R1 and the signal amplitude using previously developed models describing the image intensity of FLASH scans (Helms et al., 2008a,b; Weiskopf et al., 2011). The signal amplitude maps are proportional to PD but not corrected for RF receive field inhomogeneities and R2* related effects (Helms et al., 2008a). The effective transverse relaxation rate R2* was estimated from the logarithm of the signal intensities (from the 8 PDw multi-echo images) at different echo times using a linear regression.
Quantitative R1 maps were determined from the apparent R1 maps by correcting for local RF transmit field inhomogeneities and imperfect RF spoiling using the approach described by (Preibisch and Deichmann, 2009), which was adapted to the FLASH acquisition parameters used here. RF transmit field maps were calculated from the 3D EPI acquisition and corrected for off-resonance effects as described in (Lutti et al., 2012).
Effective PD* maps were estimated from the signal amplitude maps by adjusting for global and local receive sensitivity differences. Since the local receive sensitivity can be described as a multiplicative factor between the signal amplitude and actual PD maps and receive sensitivity profiles are smooth, the UNICORT post-processing approach can be used for correction (Weiskopf et al., 2011). UNICORT is based on the new unified segmentation approach implemented in SPM8 that combines image registration, tissue classification, and multiplicative bias correction in a single generative model and optimizes its log-likelihood objective function (Ashburner and Friston, 2005). Since the global mean PD cannot be estimated accurately with this post-processing approach, we calibrated the mean WM PD value to 69 percent units [p.u.; (Tofts, 2003)]. Note that we did not correct for R2* dependent signal decay by extrapolating the signal to TE = 0 but used the averaged multi-echo FLASH data with an effective TE = 8.45 ms. Because these resulting PD estimates still partly depended on R2*, we called this parameter the effective proton density (PD*) in line with previous reports (Lin et al., 1997).
The semi-quantitative MT saturation parameter is relatively robust against differences in relaxation times and RF transmit and receive field inhomogeneities—unlike the conventional MT ratio, which is affected by R1 and RF transmit field variations (Helms et al., 2008b, 2010). Additionally, small residual higher order dependencies of the MT saturation on the local RF transmit field were corrected using a semi-empirical approach, resulting in a corrected MT saturation value used in the further processing: is the original MT value and RFlocal the relative local flip angle compared to the nominal flip angle (derivation and details of this correction will be reported elsewhere).
Inter-Site and Intra-Site Variation
The R1, PD*, MT, R2* maps, mean PDw and mean MTw from the three different sites were registered to the mean T1w image for each volunteer and interpolated using 4th order B-splines. The registered MT maps and mean T1w images were partitioned into cerebrospinal fluid (CSF), GM and WM using unified segmentation (Ashburner and Friston, 2005). T1w images are commonly used for brain segmentation but MT maps were shown to improve the segmentation of subcortical areas (Helms et al., 2009; Tardif et al., 2009). A brain mask was defined as all voxels whose sum of CSF, GM and WM probability exceeded 90%.
The inter-site coefficient of variation (CoV) for each parameter and weighted image was calculated by dividing the standard deviation (SD, normalized by the N−1 sample size to avoid bias) by the mean across the three sites (CoV = SD/Mean) for each volunteer independently. Mean and CoV maps were generated for visual assessment.
The mean inter-site CoV was determined for four different regions-of-interest (ROIs): head of the caudate nucleus, genu of the corpus callosum, GM and WM. The first two ROIs were manually defined within the head of the caudate nucleus and genu of the corpus callosum (bilateral). The MTw images were used for the delineation of these two ROIs, because they were not used as an outcome measure in this study, which minimized any potential bias or circularity in the following analyses. The other two ROIs described the whole-brain GM and WM and were defined as all voxels with a probability of GM and WM tissue over 99% as determined by the unified segmentation of the MT maps. The high probability threshold was used to reduce partial volume effects. To formally assess whether the four different parameter maps showed a lower inter-site CoV than the T1w image, a repeated measures ANOVA (SPSS Statistics 17.0; IBM Corp., NY) was performed with the within-subject factors ROI (4 levels) and image type (parameter map/T1w image) separately for each parameter map (significance threshold of p < 0.05). Post-hoc paired t-tests were conducted in case of a significant main effect of image type. We did not correct for multiple comparisons due to comparison with the four MPM parameters, since we regarded them as separate studies/experiments.
The systematic inter-site bias was assessed as the percent difference from the mean across all three sites and volunteers. The percent bias was averaged across the GM and WM ROIs for the R1, MT, R2* parameters and T1w image. The measure for the PD* map was averaged across the GM ROI only, since the WM PD* value was calibrated to 69 p.u. in the post-processing and thus the bias might have been underestimated.
The intra-site CoV was calculated based on the SD and mean across all voxels within an ROI for each volunteer and site. The intra-site variation was estimated from the SD of the signal between voxels within the ROI [similar to previous studies (Helms et al., 2009)], since scan-rescan data were not available to estimate the intra-site variance. The validity of this particular measure of intra-site CoV is limited to areas with negligible physiological variations. Therefore, intra-site CoV estimates were only extracted in ROIs with homogeneous signal, i.e., the head of the caudate nucleus and the genu of the corpus callosum. The CoVs were averaged across sites for each volunteer.
The inter-site CoVs of GM probability maps derived from MT maps and T1w images were determined by calculating the mean and SD across sites within a GM mask. The GM mask included all voxels with a GM probability over 90%. To reduce potential bias against one of the two image types, the GM mask was derived from the corresponding type of map, i.e., either the T1w image or MT map. To formally assess whether the GM probability maps derived from MT maps had a lower inter-site CoV than the ones derived from the T1w images, a one-tailed paired t-test was performed with a significance threshold of p < 0.05.
Visual inspection of the parameter maps indicated a high quality of parameter maps with no gross artifacts (Figure 1). In contrast to the parameter maps, the T1w images showed an inhomogeneous and high inter-site CoV.
Figure 1. Axial slice through R1, PD*, MT, and R2* parameter maps and T1w images for a single volunteer. (A) Mean and (B) inter-site coefficient of variation (CoV) across the three different sites.
Table 1. Group mean and standard deviation of parameter values and T1w image intensity in different ROIs.
Figure 2. R1, PD*, MT, and R2* parameter and T1w image intensity values averaged across the three different sites. Values were determined separately for gray matter (GM), head of the caudate nucleus (CN), white matter (WM) and genu of the corpus callosum (CC). The distribution across volunteers is depicted as a whisker plot: blue box = 25/75% percentile; red line, median, black whisker, most extreme data value excluding outliers; red cross, outlier (probability < 0.01 under assumption of normally distributed data).
The inter-site CoV for R1, PD* and MT parameters ranged from 2.7 to 7.9% (Table 2 and Figure 3). The inter-site CoV of R2* ranged from 11.4% in WM to 20.3% in GM. The bias between the three different sites for the R1, PD*, MT and R2* did not exceed 3.1% (Figure 4). The intra-site CoV for R1, PD*, MT and R2* determined in the caudate head and genu of the corpus callosum were similar to the inter-site CoV, ranging from 2.4 to 15.7% (Table 3 and Figure 5).
Table 2. Group mean and standard deviation of the inter-site coefficient of variation (CoV) in percent for the parameter values and T1w image intensity in different ROIs.
Figure 3. Inter-site coefficient of variation (CoV) of R1, PD*, MT, and R2* parameter maps and T1w images intensity values. Values were determined separately for gray matter (GM), head of the caudate nucleus (CN), white matter (WM), and genu of the corpus callosum (CC). For explanation of the whisker plot, see Figure 2.
Figure 4. Inter-site bias of R1, PD*, MT, and R2* parameter and T1w image intensity values in gray and white matter for the three different sites (except for PD*, which was assessed in gray matter only). To enhance visibility of any potential inter-site bias, the percent deviation from the mean across sites is shown and plotted for all volunteers. For description of the whisker plot, see Figure 2.
Table 3. Group mean and standard deviation of the intra-site coefficient of variation (CoV) in percent for the parameter values and T1w image intensity in different ROIs.
Figure 5. Intra-site coefficient of variation (CoV) of R1, PD*, MT, and R2* parameter maps and T1w images intensity values. Values were determined separately for caudate nucleus (CN) and genu of the corpus callosum (CC). For explanation of the whisker plot, see Figure 2.
The T1w imaging showed a high inter-site CoV of ca. 15% but one of the smallest intra-site CoVs with ca. 3% (Tables 2, 3; Figures 3, 5). The inter-site CoV of the T1w images was significantly higher than that of the MT (p = 0.025, F = 12.1), R1 (p = 0.009, F = 22.5) and PD* parameter maps (p = 0.003, F = 44.0). The R2* maps did not show a significantly different inter-site CoV compared to the T1w images (p > 0.75, F = 0.1). In line with the higher inter-site CoV, an inter-site bias of up to 5.9% was measured for the T1w images, which was higher than for the quantitative parameters.
The inter-site CoV of the GM probability maps (Figure 6) was 1.40% for the MT maps [with a 95% confidence interval of CI = (1.35, 1.46%)] and 1.64% for the T1w images [with CI = (1.56, 1.73%)]. Thus, the maps derived from MT maps had a ca. 14% higher inter-site reproducibility (p < 0.05).
Figure 6. Inter-site coefficient of variation (CoV) of gray matter (GM) probabilities for all five volunteers across the three different sites. CoV was measured within the cortex as determined by thresholding the GM probability map (p > 0.9). For explanation of the whisker plot, see Figure 3.
This multi-center validation study demonstrated a high inter-site reproducibility of the MPM approach, which was significantly higher compared to conventional T1w imaging. The inter-site CoV was smaller than 8% for R1, PD* and MT maps. R2* maps exhibited a higher inter-site CoV of up to 20% (similar to standard T1w images) due to the rather short maximal echo times used to determine R2*. The inter-site bias (i.e., systematic offset) between the three different sites did not exceed 3.1% for any of the parameter values.
The inter-site CoV and inter-site bias for standard T1w images were significantly higher with 15% and 5.9%, respectively—as expected for a non quantitative imaging technique. GM probability maps based on MT parameter maps had a 14% higher inter-site reproducibility than maps based on T1w images.
Comparison of Parameter Maps and T1w Images
The quantitative parameter maps showed a higher inter-site reproducibility than the conventional T1w imaging normally used for neuoranatomical studies, since T1w imaging suffers from signal and contrast bias due to inhomogeneities in the RF transmit and receive field. At 3T RF transmit field inhomogeneities frequently exceed 20–30% of the nominal field (Lutti et al., 2010), causing significant signal and contrast bias (Thomas et al., 2005). The 32-channel receive head coil used in this study exhibits a highly inhomogeneous sensitivity profile varying by 200–300% across the brain (Wiggins et al., 2006). Both types of bias depend on the position and orientation of the head within the transmit and receive RF coils. We aimed at positioning the head as consistently as possible at the different sites by following an identical protocol. However, even relatively small deviations of 1–2 cm and 3–5° may have lead to significant signal variation in T1w images. This variation was most prominent in the superficial cortical areas close to the receiver coils reflected in a highly increased inter-site CoV (ca. 10-fold increased compared to the center; Figure 1).
Since MPM is designed to be insensitive to inhomogeneities in the RF transmit field and receive fields, the parameter maps showed a much lower CoV. Also the MPM maps did not show wide-spread spatial variation of inter-site CoV unlike the T1w images (Figure 1) or previous R1 mapping approaches that did not employ RF transmit field mapping [see e.g., Figure 3 in (Deoni et al., 2008)]. Also the inter-site bias of the quantitative parameter maps was between 2- to 10-fold smaller than of the T1w imaging (Figure 4).
Gray Matter Probability Maps
Due to their signal and contrast bias T1w image intensities are usually not directly used in analyses but used indirectly for morphometry (Ashburner et al., 2003). For example, voxel-based morphometry (VBM) segments the image into different tissue classes including GM and WM (Ashburner and Friston, 2005). After appropriate spatial normalization the local GM volume is compared voxel-wise between groups and volunteers (Ashburner, 2007). Obviously, accurate and precise segmentation is crucial but relies on T1w image intensities. Most segmentation methods account for bias in the signal intensity (Ashburner and Friston, 2005) but have to rely on assumptions about the smoothness of the bias field and cannot model or compensate for contrast bias. Helms et al. (2009) demonstrated improved segmentation of subcortical areas with MT maps compared to T1w images, since they have a high contrast-to-noise ratio and minimal bias. Similarly, we observed a 14% lower inter-site CoV in the GM probability maps based on the MT parameter maps instead of T1w images (Figure 6). Under the assumption of (independent and identically distributed) Gaussian noise the improved inter-site reproducibility would therefore reduce minimal group sizes by ca. 30% in multi-center VBM studies. Thus, segmentation and VBM results are expected to be not only improved for single site but also multi-center studies by the use of MT parameter maps.
Variations in Parameter Maps and Their Causes
Some of the residual inter-site bias and CoV in the MPM were caused by the following mechanisms. The dual flip angle mapping approach used in MPM (Helms et al., 2008a) provides signal amplitude (proportional to PD) and R1 maps that need to be corrected for RF transmit and receive field inhomogeneities. We acquired highly accurate and precise RF transmit field maps with a total error of less than ca. 3% (Lutti et al., 2010, 2012) and corrected for imperfect RF spoiling, which leads to deviations from the Ernst signal equation underlying the R1 estimation (Preibisch and Deichmann, 2009; Yarnykh, 2010). Although small, this allows for errors of up to 6% in the R1 maps due to the quadratic dependence of the estimated R1 on the local flip angle.
The RF receive field effect on the PD map was minimized by image post-processing. Unified segmentation (Ashburner and Friston, 2005) was adapted to robustly determine and correct for the multiplicative receive coil sensitivity profile in the PD maps, similar to the previously developed UNICORT approach for correcting R1 maps (Weiskopf et al., 2011). However, if multiple and small receive coils are used, the spatial sensitivity profiles may become difficult to model due to the low spatial smoothness, potentially causing insufficient correction. We calibrated the mean WM PD value to 69 p.u. (Tofts, 2003), since the global offset in PD values cannot be accurately determined by the unified segmentation step. The calibration step may introduce a bias if pathologies or physiological changes affect the PD in the majority of WM.
The PD maps were estimated from averaged images acquired at different echo times, resulting in a mean echo time of 8.45 ms. Although the mean echo time was comparatively short, it introduced a certain R2* image intensity weighting. To point out this potential bias, we called the estimated parameter effective proton density (PD*) in line with previous studies (Lin et al., 1997). Since the UNICORT post-processing step accounted for the global offset by calibrating the WM PD* to 69 p.u., the mean overall reduction in signal and thus PD* was corrected, which would otherwise spuriously reduce the PD* by ca. 15% (assuming T2* = 50 ms and TE = 8.45 ms). However, the PD* estimates may be still locally biased in regions with high R2* values due to high iron concentrations, such as parts of the basal ganglia and certain brainstem nuclei. To overcome this problem, the signal at TE = 0 may be extrapolated from the multi-echo dataset (Neeb et al., 2006). However, we decided not to apply this type of correction, since the signal extrapolation is potentially unstable and can increase the noise level significantly (Neeb et al., 2006).
The R2* parameter maps yielded an approximately 100% higher CoV than the other parameter maps. Since the longest echo time acquired in the PDw FLASH multi-echo readout was 19.7 ms, the estimation of long T2* (=1/R2*) found in GM, WM or CSF was complicated. The precision of the R2* maps may be improved by increasing the maximal echo time, but this would also prolong the total acquisition time. Moreover, R2* was estimated from the logarithm of the signal intensities (in the PDw images) at different echo times using a linear regression. Since the SNR of the different echoes varies, heteroscadisticity may have impacted the fit. The assumption of a mono-exponential signal decay described by R2* may be violated in some brain areas [e.g., suffering from susceptibility artifacts (Neeb et al., 2006)], although the high spatial resolution of 1 mm reduced the effects of susceptibility artifacts on the signal decay due to a smaller within voxel spin phase coherence loss (Weiskopf et al., 2007).
The MT saturation parameter is a measure of the saturation due to the applied off-resonance RF pulse, which is highly correlated with macromolecular content and myelin density (Helms et al., 2008b; Draganski et al., 2011). The semi-quantitative MT saturation is largely insensitive to changes in the excitation flip angle, repetition time, RF field inhomogeneities or R1 of the tissue but depends on the power and frequency offset of the MT saturation pulse (Helms et al., 2008b, 2010). For example, a higher powered saturation pulse will lead to higher MT values. Thus, the MT saturation pulse has to be kept identical for different implementations if direct comparability is desired. Calibration of the MT effect or full MT quantification may ensure data comparability even when the saturation pulses differ (Sinclair et al., 2010; Volz et al., 2010).
RF transmitter instability and differences in the RF transmitter adjustment may have also contributed to the inter-site CoV. Non-linearities in the RF amplifier may have caused additional variation (Balezeau et al., 2011), since the FLASH images were acquired with significantly different RF transmit voltages to achieve the different desired excitation flip angles. These instabilities are known to cause variations in PD, R1 and MT maps (Stikov, 2010).
Limitations and Considerations
The MPM approach poses significant challenges to the MRI scanner hardware. The fast bipolar multi-echo readout demands high gradient performance. Multi-channel RF receive coils are required to achieve good parallel imaging capability with high SNR. The image reconstruction software and hardware must be capable of handling an approximately 20-times increased data rate compared to standard T1w imaging [e.g., 3D MDEFT, (Deichmann et al., 2004)].
This multi-center validation was based on data from five healthy volunteers and three different sites. We do not believe that the small sample size compared to conventional morphometric studies affected the primary results, since most of the quantitative comparisons were based on aggregate measures (e.g., means over ROIs) providing high statistical power due to their low noise.
We note that intra-site and inter-site CoVs were found largely comparable in this study. This observation supports that little variance is added by using different MR scanners. Thus, multi-centre studies can be conducted with high sensitivity. However, we note that the MPM performance may be degraded for non-compliant volunteers. For example, patients may have difficulties to minimize head motion or body motion, which can change the magnetic field in the head and affect data quality (Versluis et al., 2010). The MPM parameters are estimated from up to three acquired datasets and are sensitive to artifacts present in any of these runs. If relevant, these problems can be addressed by prospective motion correction (Maclaren et al., 2012) and phase navigator techniques (Versluis et al., 2010).
The three sites participating in this study used the same type of MRI scanner and RF coils. Thus, inter-site bias due RF coil and gradient coil differences is expected to be lower than for studies using different scanner types. Gradient non-linearities were not addressed in this study, since the MRI scanner type used exhibits a very high linearity (Mohammadi et al., 2012). They may be addressed by post-processing if necessary (Jovicich et al., 2006). To avoid bias, special attention needs to be paid to the exact implementation of the MR pulse sequences including RF pulses and RF spoiling (Yarnykh, 2010), which can differ between software versions and scanner manufacturers. Thus, we used custom-made optimized pulse sequences in this study.
The comparability of parameter maps or images is low if they are acquired at different field strengths (e.g., at 3T and the typical clinical field strength of 1.5T). R1 and R2* significantly depend on the field strength (Oros-Peusquens et al., 2008), whereas MT and PD* show a much smaller dependence. Thus, a direct comparison and pooling of data across field strengths are complicated but may not be impossible using appropriate relaxometry models (Rooney et al., 2007). Further studies are needed to develop and validate such an approach. While multi-center MPM studies at 1.5T should not pose particular problems, studies at higher fields, such as 7T, will require additional validation, since artifacts and bias are exacerbated (Lutti et al., 2012).
It is well-known that for large sample sizes even a small bias may become significant when classic frequentist statistics (testing against a null effect) are used. Quantitative MRI may help us with this fundamental issue, since the absolute quantitative values can be used to define minimal biological effect sizes, which are site independent and more generally applicable. Thus, it can inform statistical analyses, e.g., Bayesian approaches, which take into account prior information and minimal effect sizes. This helps researchers to avoid over-interpreting biologically implausible or clinically irrelevant effects including residual small inter-site bias.
We note that the reported inter-site bias should not be regarded a measure of absolute accuracy of the parameter maps. We did not aim to determine the accuracy with this study but only comparability across sites. Determining the accuracy would have required comparisons to gold standard measurements, which usually suffer from long acquisition times and low resolution (Weiskopf et al., 2011). The measured parameter values were relatively well in line with previous studies. For example, the R1 = 0.61/1.04 s−1 in GM/WM was similar to the R1 measured previously 0.63/1.19 s−1 (Wright et al., 2008). PD* estimates of 84.4/82.7 p.u in GM/caudate nucleus were in line with previous studies reporting 81.1/81.5 p.u. (Volz et al., 2012) and 82.2/84.8 p.u. (Neeb et al., 2008) for the same structures. Measurements of R2* in WM were also similar to previous studies [21.0 s−1 compared to 19.5 s−1 (Baudrexel et al., 2009) and 21.7 s−1 (Martin et al., 2008)]. However, R2* in the head of the caudate appeared reduced compared to previous studies [18.2 s−1 compared to 26.2 s−1 (Martin et al., 2008)]. However, a comparison with literature values should be interpreted with care, since estimates vary considerably between studies. For example, even recent studies differ by more than 15% in estimated R1 values (Oros-Peusquens et al., 2008; Wright et al., 2008). This is probably caused by different resolutions, varying ROI definitions and varying biases present in the studies (Tofts, 2003).
We expect some bias in the PD* and R2* maps in areas suffering from susceptibility artifacts, since the mono-exponential decay model used here was violated in these areas and we only partially corrected the PD* maps for R2* related effects (e.g., increased inter-site CoV in the frontal cortex, Figure 1). Thus, e.g., differences in shim and head positioning may have affected R2* values.
Signal variations within ROIs were used as proxies for intra-site CoV due to the absence of scan-rescan data. Although we chose small homogeneous ROIs (e.g., head of the caudate nucleus), it is likely that some variability in tissue architecture within the ROIs contributed to the intra-site CoV, spuriously increasing the estimates of the intra-site noise level. Furthermore, correlations between receive channels of the 32-channel RF receive head coil might have biased noise estimates (Hutton et al., 2012). This also explains why the intra-site CoV was paradoxically higher than the inter-site CoV for R2* in the corpus callosum (Figures 3, 5).
The multi-echo FLASH based T1w imaging, which was compared to MPM here, is less widely used since the advent of MPRAGE and MDEFT based T1w imaging (Mugler and Brookeman, 1990; Deichmann et al., 2000, 2004), since they offer a higher GM/WM contrast-to-noise ratio in the cortex (Tardif et al., 2009). However, neither MPRAGE nor MDEFT completely corrects for RF field bias. Thus, we expect a similarly increased inter-site CoV for these MR pulse sequences as well.
We have introduced and validated a MPM approach for multi-center studies at 3T. It provides high-resolution maps of R1, PD*, MT and R2* and thus detailed insights into the brain microstructure in a clinically feasible acquisition time. An optimized multi-echo FLASH acquisition with low artifact level and high signal-to-noise ratio combined with correction for RF transmit and receive field inhomogeneities results in accurate and precise quantitative measures. The resulting high comparability of MPM data across sites and time points facilitates multi-center studies and federating large datasets.
Conflict of Interest Statement
Edward T. Bullmore is employed half-time by GlaxoSmithKline (GSK) and half-time by the University of Cambridge; he holds stock in GSK. The other 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.
This work was supported by a Strategic Award from the Wellcome Trust to the Neuroscience in Psychiatry Network (NSPN; I. M. Goodyer, P. B. Jones, Edward T. Bullmore, P. Fonagy and R. J. Dolan; 095844/Z/11/Z). The WTCN is supported by the Wellcome Trust (079866/Z/06/Z). Scanning at the MRC CBSU was funded by the Medical Research Council (Dr. Tim Dalgleish). We thank the volunteers for their time and effort participating in this study and travelling from site to site.
Ashburner, J., Csernansky, J. G., Davatzikos, C., Fox, N. C., Frisoni, G. B., and Thompson, P. M. (2003). Computer-assisted imaging to assess brain structure in healthy and diseased brains. Lancet Neurol. 2, 79–88. doi: 10.1016/S1474-4422(03)00304-1
Barker, G. J., Schreiber, W. G., Gass, A., Ranjeva, J. P., Campi, A., van Waesberghe, J. H. T. M., et al. (2005). A standardised method for measuring magnetisation transfer ratio on MR imagers from different manufacturers–the EuroMT sequence. MAGMA 18, 76–80. doi: 10.1007/s10334-004-0095-z
Baudrexel, S., Volz, S., Preibisch, C., Klein, J. C., Steinmetz, H., Hilker, R., et al. (2009). Rapid single-scan T2*-mapping using exponential excitation pulses and image-based correction for linear background gradients. Magn. Reson. Med. 62, 263–268. doi: 10.1002/mrm.21971
Bauer, C. M., Jara, H., and Killiany, R. (2010). Whole brain quantitative T2 MRI across multiple scanners with dual echo FSE: applications to AD, MCI, and normal aging. Neuroimage 52, 508–514. doi: 10.1016/j.neuroimage.2010.04.255
Deichmann, R., Schwarzbauer, C., and Turner, R. (2004). Optimisation of the 3D MDEFT sequence for anatomical brain imaging: technical implications at 1.5 and 3 T. Neuroimage 21, 757–767. doi: 10.1016/j.neuroimage.2003.09.062
Deoni, S. C. L., Williams, S. C. R., Jezzard, P., Suckling, J., Murphy, D. G. M., and Jones, D. K. (2008). Standardized structural magnetic resonance imaging in multicentre studies using quantitative T1 and T2 imaging at 1.5 T. Neuroimage 40, 662–671. doi: 10.1016/j.neuroimage.2007.11.052
Dick, F., Tierney, A., Lutti, A., Josephs, O., Sereno, M. I., and Weiskopf, N. (2012). In vivo functional and myeloarchitectonic mapping of human primary auditory areas. J. Neurosci. 32, 16095–16105. doi: 10.1523/JNEUROSCI.1712-12.2012
Draganski, B., Ashburner, J., Hutton, C., Kherif, F., Frackowiak, R. S. J., Helms, G., et al. (2011). Regional specificity of MRI contrast parameter changes in normal ageing revealed by voxel-based quantification (VBQ). Neuroimage 55, 1423–1434. doi: 10.1016/j.neuroimage.2011.01.052
Ecker, C. (2012). Brain anatomy and its relationship to behavior in adults with autism spectrum disorder: a multicenter magnetic resonance imaging study. Arch. Gen. Psychiatry 69, 195. doi: 10.1001/archgenpsychiatry.2011.1251
Focke, N. K., Helms, G., Kaspar, S., Diederich, C., Tóth, V., Dechent, P., et al. (2011). Multi-site voxel-based morphometry—not quite there yet. Neuroimage 56, 1164–1170. doi: 10.1016/j.neuroimage.2011.02.029
Garrido, L., Furl, N., Draganski, B., Weiskopf, N., Stevens, J., Tan, G. C.-Y., et al. (2009). Voxel-based morphometry reveals reduced grey matter volume in the temporal cortex of developmental prosopagnosics. Brain 132, 3443–3455. doi: 10.1093/brain/awp271
Helms, G., Dathe, H., Kallenberg, K., and Dechent, P. (2008b). High-resolution maps of magnetization transfer with inherent correction for RF inhomogeneity and T1 relaxation obtained from 3D FLASH MRI. Magn. Reson. Med. 60, 1396–1407. doi: 10.1002/mrm.21732
Helms, G., Dathe, H., and Dechent, P. (2010). Modeling the influence of TR and excitation flip angle on the magnetization transfer ratio (MTR) in human brain obtained from 3D spoiled gradient echo MRI. Magn. Reson. Med. 64, 177–185. doi: 10.1002/mrm.22379
Helms, G., Dathe, H., Weiskopf, N., and Dechent, P. (2011). Identification of signal bias in the variable flip angle method by linear display of the algebraic ernst equation. Magn. Reson. Med. 66, 669–677. doi: 10.1002/mrm.22849
Helms, G., and Dechent, P. (2009). Increased SNR and reduced distortions by averaging multiple gradient echo signals in 3D FLASH imaging of the human brain at 3T. J. Magn. Reson. Imaging 29, 198–204. doi: 10.1002/jmri.21629
Helms, G., Draganski, B., Frackowiak, R., Ashburner, J., and Weiskopf, N. (2009). Improved segmentation of deep brain grey matter structures using magnetization transfer (MT) parameter maps. Neuroimage 47, 194–198. doi: 10.1016/j.neuroimage.2009.03.053t
Hutton, C., Balteau, E., Lutti, A., Josephs, O., and Weiskopf, N. (2012). Modelling temporal stability of EPI time series using magnitude images acquired with multi-channel receiver coils. PLoS ONE 7:e52075. doi: 10.1371/journal.pone.0052075
Jovicich, J., Czanner, S., Greve, D., Haley, E., van der Kouwe, A., Gollub, R., et al. (2006). Reliability in multi-site structural MRI studies: effects of gradient non-linearity correction on phantom and human data. Neuroimage 30, 436–443. doi: 10.1016/j.neuroimage.2005.09.046
Kloppel, S., Stonnington, C. M., Barnes, J., Chen, F., Chu, C., Good, C. D., et al. (2008a). Accuracy of dementia diagnosis: a direct comparison between radiologists and a computerized method. Brain 131, 2969–2974. doi: 10.3410/f.1141939.599046
Kloppel, S., Stonnington, C. M., Chu, C., Draganski, B., Scahill, R. I., Rohrer, J. D., et al. (2008b). Automatic classification of MR scans in Alzheimer's disease. Brain 131, 681–689. doi: 10.1093/brain/awm319
Lin, W., Paczynski, R. P., Venkatesan, R., He, Y. Y., Powers, W. J., Hsu, C. Y., et al. (1997). Quantitative regional brain water measurement with magnetic resonance imaging in a focal ischemia model. Magn. Reson. Med. 38, 303–310. doi: 10.1002/mrm.1910380221
Lutti, A., Hutton, C., Finsterbusch, J., Helms, G., and Weiskopf, N. (2010). Optimization and validation of methods for mapping of the radiofrequency transmit field at 3T. Magn. Reson. Med. 64, 229–238. doi: 10.1002/mrm.22421
Lutti, A., Stadler, J., Josephs, O., Windischberger, C., Speck, O., Bernarding, J., et al. (2012). Robust and fast whole brain mapping of the RF transmit field B1 at 7T. PLoS ONE 7:e32379. doi: 10.1371/journal.pone.0032379
Maclaren, J., Armstrong, B. S. R., Barrows, R. T., Danishad, K. A., Ernst, T., Foster, C. L., et al. (2012). Measurement and correction of microscopic head motion during magnetic resonance imaging of the brain. PLoS ONE 7:e48088. doi: 10.1371/journal.pone.0048088
Martin, W. R. W., Wieler, M., and Gee, M. (2008). Midbrain iron content in early Parkinson disease A potential biomarker of disease status. Neurology 70, 1411–1417. doi: 10.1212/01.wnl.0000286384.31050.b5
Mohammadi, S., Nagy, Z., Möller, H. E., Symms, M. R., Carmichael, D. W., Josephs, O., et al. (2012). The effect of local perturbation fields on human DTI: characterisation, measurement and correction. Neuroimage 60, 562–570. doi: 10.1016/j.neuroimage.2011.12.009
Oros-Peusquens, A. M., Laurila, M., and Shah, N. J. (2008). Magnetic field dependence of the distribution of NMR relaxation times in the living human brain. MAGMA 21, 131–147. doi: 10.1007/s10334-008-0107-5
Pardoe, H., Pell, G. S., Abbott, D. F., Berg, A. T., and Jackson, G. D. (2008). Multi-site voxel-based morphometry: methods and a feasibility demonstration with childhood absence epilepsy. Neuroimage 42, 611–616. doi: 10.1016/j.neuroimage.2008.05.007
Preibisch, C., and Deichmann, R. (2009). Influence of RF spoiling on the stability and accuracy of T1 mapping based on spoiled FLASH with varying flip angles. Magn. Reson. Med. 61, 125–135. doi: 10.1002/mrm.21776
Pruessmann, K. P., Weiger, M., Scheidegger, M. B., and Boesiger, P. (1999). SENSE: sensitivity encoding for fast MRI. Magn. Reson. Med. 42, 952–962. doi: 10.1002/(SICI)1522-2594(199911)42:5<952::AID-MRM16>3.0.CO;2-S
Rooney, W. D., Johnson, G., Li, X., Cohen, E. R., Kim, S.-G., Ugurbil, K., et al. (2007). Magnetic field and tissue dependencies of human brain longitudinal 1H2O relaxation in vivo. Magn. Reson. Med. 57, 308–318. doi: 10.1002/mrm.21122
Ropele, S., Filippi, M., Valsasina, P., Korteweg, T., Barkhof, F., Tofts, P. S., et al. (2005). Assessment and correction of B1-induced errors in magnetization transfer ratio measurements. Magn. Reson. Med. 53, 134–140. doi: 10.1002/mrm.20310
Schilling, C., Kühn, S., Paus, T., Romanowski, A., Banaschewski, T., and Barbot, A. (2013). Cortical thickness of superior frontal cortex predicts impulsiveness and perceptual reasoning in adolescence. Mol. Psychiatry 18, 624–630. doi: 10.1038/mp.2012.56
Sereno, M. I., Lutti, A., Weiskopf, N., and Dick, F. (2012). Mapping the human cortical surface by combining quantitative T1 with retinotopy. Cereb. Cortex. doi: 10.1093/cercor/bhs213. [Epub ahead of print].
Sinclair, C. D. J., Samson, R. S., Thomas, D. L., Weiskopf, N., Lutti, A., Thornton, J. S., et al. (2010). Quantitative magnetization transfer in in vivo healthy human skeletal muscle at 3 T. Magn. Reson. Med. 64, 1739–1748. doi: 10.1002/mrm.22562
Stonnington, C. M., Tan, G., Klöppel, S., Chu, C., Draganski, B., Jack, C. R. Jr., et al. (2008). Interpreting scan data acquired from multiple scanners: a study with Alzheimer's disease. Neuroimage 39, 1180–1185. doi: 10.1016/j.neuroimage.2007.09.066
Tardif, C. L., Collins, D. L., and Pike, G. B. (2009). Sensitivity of voxel-based morphometry analysis to choice of imaging protocol at 3 T. Neuroimage 44, 827–838. doi: 10.1016/j.neuroimage.2008.09.053
Thomas, D. L., De Vita, E., Deichmann, R., Turner, R., and Ordidge, R. J. (2005). 3D MDEFT imaging of the human brain at 4.7 T with reduced sensitivity to radiofrequency inhomogeneity. Magn Reson. Med. 53, 1452–1458. doi: 10.1002/mrm.20482
Tofts, P. S., Steens, S. C. A., Cercignani, M., Admiraal-Behloul, F., Hofman, P. A. M., van Osch, M. J. P., et al. (2006). Sources of variation in multi-centre brain MTR histogram studies: body-coil transmission eliminates inter-centre differences. MAGMA 19, 209–222. doi: 10.1007/s10334-006-0049-8
Versluis, M. J., Peeters, J. M., van Rooden, S., van der Grond, J., van Buchem, M. A., Webb, A. G., et al. (2010). Origin and reduction of motion and f0 artifacts in high resolution T2*-weighted magnetic resonance imaging: application in Alzheimer's disease patients. Neuroimage 51, 1082–1088. doi: 10.1016/j.neuroimage.2010.03.048
Volz, S., Nöth, U., Jurcoane, A., Ziemann, U., Hattingen, E., and Deichmann, R. (2012). Quantitative proton density mapping: correcting the receiver sensitivity bias via pseudo proton densities. Neuroimage 63, 540–552. doi: 10.1016/j.neuroimage.2012.06.076
Volz, S., Nöth, U., Rotarska-Jagiela, A., and Deichmann, R. (2010). A fast B1-mapping method for the correction and normalization of magnetization transfer ratio maps at 3 T. Neuroimage 49, 3015–3026. doi: 10.1016/j.neuroimage.2009.11.054
Weiskopf, N., Hutton, C., Josephs, O., Turner, R., and Deichmann, R. (2007). Optimized EPI for fMRI studies of the orbitofrontal cortex: compensation of susceptibility-induced gradients in the readout direction. MAGMA 20, 39–49. doi: 10.1007/s10334-006-0067-6
Weiskopf, N., Lutti, A., Helms, G., Novak, M., Ashburner, J., and Hutton, C. (2011). Unified segmentation based correction of R1 brain maps for RF transmit field inhomogeneities (UNICORT). Neuroimage 54, 2116–2124. doi: 10.1016/j.neuroimage.2010.10.023
Wiggins, G. C., Triantafyllou, C., Potthast, A., Reykowski, A., Nittka, M., and Wald, L. L. (2006). 32-channel 3 Tesla receive-only phased-array head coil with soccer-ball element geometry. Magn. Reson. Med. 56, 216–223. doi: 10.1002/mrm.20925
Wright, P. J., Mougin, O. E., Totman, J. J., Peters, A. M., Brookes, M. J., Coxon, R., et al. (2008). Water proton T1 measurements in brain tissue at 7, 3, and 1.5 T using IR-EPI, IR-TSE, and MPRAGE: results and optimization. MAGMA 21, 121–130. doi: 10.1007/s10334-008-0104-8
Keywords: multi-center, T1, PD, MT, T2*, 3T, MPM, qMRI
Citation: Weiskopf N, Suckling J, Williams G, Correia MM, Inkster B, Tait R, Ooi C, Bullmore ET and Lutti A (2013) Quantitative multi-parameter mapping of R1, PD*, MT, and R2* at 3T: a multi-center validation. Front. Neurosci. 7:95. doi: 10.3389/fnins.2013.00095
Received: 04 February 2013; Paper pending published: 25 March 2013;
Accepted: 18 May 2013; Published online: 10 June 2013.
Edited by:Ching-Po Lin, National Yang-Ming University, Taiwan
Reviewed by:Pierre Bellec, University of Montreal, Canada
Chao Yi-Ping, Chang Gung University, Taiwan
Copyright © 2013 Weiskopf, Suckling, Williams, Correia, Inkster, Tait, Ooi, Bullmore and Lutti. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Nikolaus Weiskopf, Wellcome Trust Centre for Neuroimaging, UCL Institute of Neurology, University College London, 12 Queen Square, London WC1N 3BG, UK e-mail: email@example.com