High-resolution MR imaging of the human brainstem in vivo at 7 Tesla

The human brainstem, which comprises a multitude of axonal nerve fibers and nuclei, plays an important functional role in the human brain. Depicting its anatomy non-invasively with high spatial resolution may thus in turn help to better relate normal and pathological anatomical variations to medical conditions as well as neurological and peripheral functions. We explored the potential of high-resolution magnetic resonance imaging (MRI) at 7 T for depicting the intricate anatomy of the human brainstem in vivo by acquiring and generating images with multiple contrasts: T₂-weighted images, quantitative maps of longitudinal relaxation rate (R₁ maps) and effective transverse relaxation rate ($R_{\mathrm{2}}^{\mathrm{*}}$ maps), magnetic susceptibility maps, and direction-encoded track-density images. Images and quantitative maps were compared with histological stains and anatomical atlases to identify nerve nuclei and nerve fibers. Among the investigated contrasts, susceptibility maps displayed the largest number of brainstem structures. Contrary to R₁ maps and T₂-weighted images, which showed rather homogeneous contrast, $R_{\mathrm{2}}^{\mathrm{*}}$ maps, magnetic susceptibility maps, and track-density images clearly displayed a multitude of smaller and larger fiber bundles. Several brainstem nuclei were identifiable in sections covering the pons and medulla oblongata, including the spinal trigeminal nucleus and the reticulotegmental nucleus on magnetic susceptibility maps as well as the inferior olive on R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps. The substantia nigra and red nuclei were visible in all contrasts. In conclusion, high-resolution, multi-contrast MR imaging at 7 T is a versatile tool to non-invasively assess the individual anatomy and tissue composition of the human brainstem.

Introduction

The brainstem is the primary relay center for afferent and efferent connections between the cerebral cortex, the cerebellum, and the spinal cord (Sutin and Carpenter, 1976). It plays an important role in the regulation of vital functions (e.g., respiratory function, cardio-vascular function, nausea) and coordinates motor control signals sent from the brain to the body. Anatomically, the human brainstem consists of three connected parts, medulla oblongata, pons, and midbrain, and is composed of a multitude of axonal nerve fibers as well as cranial and non-cranial nerve nuclei. Due to the spatial concentration of important neural structures in this relatively small brain region, pathological variations or lesions in the brainstem can introduce severe and multiple neurological effects (Donaldson et al., 2006). The brainstem is also known to be involved in a number of degenerative diseases, such as Alzheimer’s disease, Parkinson’s disease, and Huntington’s disease (Urban and Caplan, 2011). Hence, detailed non-invasive depiction of its morphology and anatomy may help to better relate normal and pathological anatomical variations to medical conditions as well as neurological and peripheral functions.

Magnetic resonance imaging (MRI) at ultrahigh magnetic field strength (B₀ ≥ 7 T) is a powerful means to assess non-invasively normal and abnormal brain tissue with high spatial resolution (Li et al., 2006; Duyn et al., 2007; Deistung et al., 2008; Trampel et al., 2011; Moser et al., 2012; Turner, 2012, 2013; Eichner et al., 2013; Marques and Gruetter, 2013). MRI provides a variety of qualitative and quantitative tissue contrasts, mainly reflecting nuclear relaxation (T₁, T₂, and $T_{\mathrm{2}}^{\mathrm{*}}$, or equivalently R₁, R₂, and $R_{\mathrm{2}}^{\mathrm{*}}$), diffusion, and magnetic susceptibility. In contrast-sensitized MR images (e.g., T₁-weighted, T₂-weighted, or $T_{\mathrm{2}}^{\mathrm{*}}$-weighted images), however, the measured signal is an intricate function of proton density, longitudinal (T₁) and transverse relaxation (T₂) and depends on the applied MRI sequence parameters. Thus, it may be difficult to interpret the underlying tissue composition and structure based on the intensity of these images. Quantitative MRI techniques, on the other hand, provide information that is intrinsically more tissue specific and consequently less dependent on the chosen scan parameters. Hence, quantification of different physical tissue properties by using quantitative MRI techniques also increases intra- and inter-individual comparability and enables objective measurements of disease related brain changes (Tofts, 2003).

The most common approach of quantitative MRI is relaxometry, i.e., mapping of the longitudinal (R₁) and/or transverse (R₂, $R_{\mathrm{2}}^{\prime }$, $R_{\mathrm{2}}^{\mathrm{*}}$) relaxation rates. In the brain these relaxation rates are influenced predominantly by both myelin (Koenig, 1991; Bender and Klose, 2010; Langkammer et al., 2012a; Lee et al., 2012) and tissue iron (Ogg and Steen, 1998; Langkammer et al., 2010). An additional, recently introduced novel quantitative MR method that is unique in its sensitivity to both tissue constituents is quantitative susceptibility mapping (QSM) (Schweser et al., 2011; Reichenbach, 2012). QSM is a post-processing method that uses the phase signal of complex-valued gradient (recalled) echo (GRE) data to produce maps of tissue magnetic susceptibility. Quantification of brain tissue iron content (Wharton et al., 2010; Schweser et al., 2011; Bilgic et al., 2012; Langkammer et al., 2012b; Zheng et al., 2013) and blood oxygenation in vivo (Haacke et al., 2010) as well as assessment of brain myelination (Liu et al., 2011; Li et al., 2012a,b) has recently been reported by applying QSM. Excellent anatomical delineation of cortical and deep gray matter structures and substructures has been demonstrated with QSM (Schweser et al., 2011, 2012; Deistung et al., 2013a), particularly in concert with relaxation rate mapping (Sigalovsky et al., 2006; Fukunaga et al., 2010; Geyer et al., 2011; Lebel et al., 2012; Deistung et al., 2013a), thus rendering this combined approach highly promising for improved delineation of the anatomical structure of the brainstem.

Nerve fibers and their complex network are usually imaged and analyzed by exploiting the translational molecular diffusion of protons. So far, most studies have applied diffusion tensor imaging (DTI) (Basser and Pierpaoli, 1996) for anatomical imaging of the brainstem (Stieltjes et al., 2001; Nagae-Poetscher et al., 2004; Salamon et al., 2005). DTI, however, is commonly limited by the relatively coarse spatial resolution available, which makes detailed analysis of smaller fiber tracts difficult, and often fails in regions of crossing fibers (Jones et al., 2013). One recently suggested approach to overcome this limitation is track-density imaging (TDI), which combines high angular resolution diffusion imaging (HARDI) with fiber tracking information to delineate structures beyond the image resolution of the acquired diffusion-weighted data (Calamante et al., 2010). TDI has been already applied successfully to delineate fiber pathways in the cerebrum and cerebellum, and the substructure of the thalamus (Calamante et al., 2010, 2013).

The present study explores the potential of these high-resolution MRI approaches to resolve non-invasively the intricate anatomy of the human brainstem at 7 T in vivo. To this end, T₂-weighted images, quantitative R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and magnetic susceptibility maps as well as track-density images were compared with respect to their ability of depicting nerve nuclei and nerve fibers. Furthermore, R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and magnetic susceptibility values are presented for selected anatomical regions. It is shown that the combination of multiple image contrasts provides a distinctly improved portrayal of the morphology of the brainstem.

Materials and Methods

Data Acquisition

The study was approved by the internal Institutional Review Board of the Max Planck Institute in Leipzig and written informed consent was obtained from all participating subjects. A total of six healthy volunteers (four male and two female; 27.7 ± 3 years) were examined on a 7 T human whole body MRI system (Siemens Healthcare, Erlangen, Germany) using a 24-channel head array coil (NOVA Medical Inc., Wilmington, MA, USA). Two dielectric pads containing deuterated water (99%, Sigma Aldrich GmbH, Germany) enriched with calcium titanate (Alfa Aesar GmbH and Co KG, Karlsruhe, Germany) were placed on the left and right side of the subjects’ head to increase both field strength and homogeneity of the transmit radio frequency field, $B_1^{+}$, in the area of the brainstem (Teeuwisse et al., 2012).

Coronal T₂-weighted images were acquired using a two-dimensional gradient-echo and spin-echo (GRASE) sequence (Feinberg and Oshio, 1991; Oshio and Feinberg, 1991) with echo time (TE) = 35 ms, repetition time (TR) = 10,000 ms, bandwidth (BW) = 343 Hz/px, acquisition matrix = 384 × 384 and in-plane resolution of 0.53 mm × 0.53 mm. Thirty-five images were collected with a slice thickness of 0.6 mm and a gap between two adjacent slices of 0.3 mm. For each TR interval there were three 180°-refocusing pulses, and between each 180°-pulse there were three gradient recalled echoes, resulting in an effective echo train length of nine and an acquisition time (TA) of 7:22 min:s. The GRASE sequence was applied to account for the high specific absorption rate (SAR) that substantially restricts T₂-weighted imaging with turbo-spin-echo sequences at ultra-high magnetic fields.

High-resolution R₁-mapping was performed based on the data acquired with the MP2RAGE sequence (Marques et al., 2010), a magnetization-prepared rapid gradient-echo (MP-RAGE) sequence with two different inversion times (TI). The acquisition parameters of the coronal MP2RAGE scan included TI₁ = 900 ms, TI₂ = 2750 ms, TE = 3.7 ms, TR = 5000 ms, BW = 240 Hz/px, acquisition matrix = 320 × 260 × 240, voxel size = 0.6 mm × 0.6 mm × 0.6 mm, and TA = 17:40 min:s. Magnetization preparation was achieved with a tailored radiofrequency pulse to take into account the heterogeneity of the radiofrequency transmit field (Hurley et al., 2010).

Coronal, multi-echo, three-dimensional gradient-echo imaging was carried out to compute quantitative $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps. To this end, three echoes with monopolar readout were recorded (TE₁ = 11 ms, TE₂ = 21 ms, and TE₃ = 31 ms), TR = 43 ms, flip angle (FA) = 12.5°, BW = 149 Hz/px, acquisition matrix = 448 × 364 × 104, and voxel size = 0.43 mm × 0.43 mm × 0.43 mm. Data were collected with 75% and 87.5% partial Fourier along phase and slice encoding direction, respectively, resulting in an acquisition time of 17:48 min:s. Choosing a readout BW of 149 Hz/px yielded single-echo images with a rather high signal-to-noise ratio (SNR), however, at the expense of a rather long inter-echo distance. The signal decay was sampled with three echoes only to keep acquisition time below 20 min.

Finally, diffusion-weighted imaging (DWI) data were acquired in sagittal orientation with two-dimensional, single-shot, spin-echo echo-planar imaging (EPI), 60 diffusion encoding directions each with a b-value of 1000 s/mm² and seven volumes with a b-value of 50 s/mm², TE = 64 ms, TR = 10,000 ms, BW = 1050 Hz/px, 86 contiguous slices with an acquisition matrix of 170 × 170, voxel size = 1.2 mm × 1.2 mm × 1.2 mm, and TA = 14:20 min:s. To reduce geometric distortions partial parallel under sampling [GRAPPA (Griswold et al., 2002)] with an acceleration factor of 3 and 45 reference lines was applied in the phase-encoding direction.

In addition, a field-map was acquired with a 2D double-echo gradient-echo sequence to correct residual geometric distortions of the diffusion-weighted images caused by susceptibility differences between air-bone and air-tissue interfaces. The field-map was recorded with TE₁ = 6 ms, TE₂ = 7 ms, TR = 1000 ms, FA = 50°, 51 contiguous slices, matrix = 102 × 102, voxel size = 2 mm × 2 mm × 2 mm, and TA = 3:22 min:s.

Data Processing

Relaxometry

Maps of the longitudinal relaxation times were calculated directly at the MRI scanner by applying the MP2RAGE reconstruction framework provided by the manufacturer. The relaxation times were subsequently converted to relaxation rates according to the relation R₁ = 1/T₁ to facilitate better depiction of anatomical structures.

Maps of the effective transverse relaxation rate, $R_{\mathrm{2}}^{\mathrm{*}}$, were obtained using the power method (McGibney and Smith, 1993; Miller and Joseph, 1993), i.e., the squared magnitude signal decay, S(r→, TE)², of the 3D multi-echo GRE scan was used for the regression (McGibney and Smith, 1993; Miller and Joseph, 1993):

$S{\left(\overrightarrow{r},TE\right)}^2=S_0{\left(\overrightarrow{r}\right)}^2\cdot \exp \left(-2\cdot TE\cdot R_2^{*}\left(\overrightarrow{r}\right)\right),$ (1)

where $S_0\left(\overrightarrow{r}\right)$ is the signal intensity at TE = 0 and $\overrightarrow{r}$ is the position vector. This approach reduces contamination of Rician noise to the fit and provides more accurate relaxation rates than fitting a mono-exponential model to the magnitude signal decay $S(\overrightarrow{r},\text{\hspace{0.17em}TE})$ (van der Weerd et al., 2000).

Quantitative susceptibility mapping

Single-channel GRE magnitude images were combined using the sum-of-squares method (Roemer et al., 1990), whereas single-channel GRE phase images were combined by taking into account the channel-dependent phase offset, which was estimated from the single-channel images at different echo times by voxel-wise linear fitting of the phase evolution (Robinson et al., 2011). Phase aliasing in the combined GRE phase data was resolved using a 3D Laplacian-based phase unwrapping algorithm (Schofield and Zhu, 2003) and phase images of different echo times were then combined in a CNR-optimized manner according to (Wu et al., 2012). Background phase contributions were eliminated with sophisticated harmonic artifact reduction for phase data (SHARP) (Schweser et al., 2011) (regularization parameter: 0.05) combining 10 different spherical kernels with varying radii ranging from 0.43 to 4.3 mm (Li et al., 2011). Susceptibility mapping was performed using homogeneity enabled incremental dipole inversion (HEIDI), which incorporates a priori information extracted from the complex GRE signal to address the ill-posed nature of QSM (Schweser et al., 2012). Since the calculated susceptibility values represent relative rather than absolute values (Schweser et al., 2011), susceptibility differences were specified with respect to a homogenous region of normal appearing white matter (NAWM, semioval center).

Track-density imaging

Geometric distortion of diffusion-weighted images was corrected using FUGUE (FSL toolbox, FMRIB, Oxford, England) based on the additionally acquired field-map. The unwarped DWI data were then processed using the MRtrix software package (Brain Research Institute, Melbourne, VIC, Australia) to produce track-density images as proposed in Calamante et al. (2010). Processing steps performed with MRtrix included the calculation of tensor data as well as fractional anisotropy maps, estimation of the response function, spherical deconvolution, probabilistic tractography and track-density mapping. Only voxels with fractional anisotropy, derived from the tensor model, that exceeded a threshold of 0.7 were taken into account for estimating the coefficients of the response function based on DWI data (Tournier et al., 2004). These coefficients and a maximum harmonic order of eight were used for constrained spherical deconvolution of DWI data (Tournier et al., 2007). This approach allowed modeling of multiple fiber populations within an imaging voxel, thereby overcoming the well-known limitation of the diffusion tensor model in regions with crossing fibers. Probabilistic streamline tracking was carried out using the second order integration over fiber orientation distributions (iFOD2) algorithm (Tournier et al., 2010) by seeding randomly throughout the predefined mask covering the whole brainstem and using the following relevant parameters: step size = 0.1 mm, number of fibers = 800,000, minimum curvature radius = 0.4 mm and a maximum number of 2000 trials for each point. In order to generate the track-density images, a virtual regularly spaced grid was superposed with a resolution of 0.43 mm × 0.43 mm × 0.43 mm on the tractography results and the total number of fiber tracks present in each grid element was calculated. Hereby, the super-resolution property is achieved by utilizing the additional information provided by modeling the fiber tracking results (Calamante et al., 2010). The track-density images were calculated on a virtual isotropic resolution of 0.43 mm to facilitate comparison with the $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps. Finally, the track-density images were direction-encoded by assigning each grid element an RGB color that represents the local fiber orientation as given by averaging the colors of all the fiber tracking (streamline) segments contained within each grid element.

Data Analysis

All MRI contrasts were converted into an identical space for visual inspection. To this end, the R₁ maps, track-density images and T₂-weighted images were linearly transformed into the space of the susceptibility maps based on orientation information encoded in the DICOM images. The susceptibility maps were chosen as target to maintain the high spatial resolution of the susceptibility and $R_{\mathrm{2}}^{\mathrm{*}}$ maps. Manual reformatting of quantitative susceptibility maps in multi-planar orientations was performed with Freeview of the Freesurfer software library¹ and the resulting rigid transformation matrix was applied to the previously registered MR image contrasts and to the $R_{\mathrm{2}}^{\mathrm{*}}$ maps. The MR images were averaged across three adjacent slices for signal-to-noise improvement, resulting in an effective slice thickness of 1.3 mm.

For qualitative comparison horizontal histological sections obtained from a 36-year-old male stained for myelin and for cells were consulted. These histological images were courtesy of the Brain Biodiversity Bank of the Michigan State University^{2 , 3} with support from the US National Science Foundation. An experienced neuroanatomist (UB; experience more than 15 years) assessed the MR images and identified anatomical structures and substructures if they coincided with a histoarchitectonic atlas (Bergman et al., 1989; Naidich et al., 2009; Paxinos et al., 2012) and if they were identifiable with respect to the surrounding tissue. For this qualitative analysis the window and level settings of the MR images were adjusted freely.

For quantitative characterization the contrast-to-noise ratio (CNR) of selected brain structures [red nucleus, substantia nigra, central tegmental tract, superior colliculus, inferior colliculus, reticulotegmental nucleus, middle cerebellar peduncle (MCP), superior cerebellar peduncle (SCP), and transverse pontine fibers] was determined across the different MR image contrasts (R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, susceptibility, and T₂-weighted contrast) according to:

$\mathrm{CN}{\mathrm{R}}_{\mathrm{cy}}=\left|\frac{S_{\mathrm{cy}}-S_{\mathrm{cys}}}{\mathrm{\sigma }}\right|\cdot {\mathrm{\eta }}_{\mathrm{c}}.$ (2)

CNR_cy denotes the CNR for MR image contrast c and anatomic region y. Mean image intensities measured in a volume-of-interest (VOI) of the investigated structure y and its surrounding tissue are denoted by S_cy and S_cys, respectively. The standard deviation of the signal intensities measured in a VOI of normal appearing white matter (semioval center) was used as an estimate of noise, σ. To account for variations in spatial resolution of the different investigated contrasts, the CNR was normalized to a voxel volume of 1 mm × 1 mm × 1 mm as suggested in (Nölte et al., 2012). To this end, the normalization factor η_c is introduced in Eq. 2, which is defined as the ratio between 1 mm³ and the acquired voxel volume of image contrast c. CNR was not determined on diffusion-weighted data because of their incompletely compensated geometric distortions.

To calculate the CNR and to measure both the relaxation rates and magnetic susceptibility, VOIs were identified based on magnetic susceptibility, R₁, and $R_{\mathrm{2}}^{\mathrm{*}}$ maps in the coordinate space of the susceptibility maps for each subject in both hemispheres. The VOIs were then transferred to the space of the R₁ maps and the T₂-weighted data, respectively, by employing orientation information encoded in the DICOM images and nearest-neighbor interpolation. Mean values and standard deviation of CNR and of the quantitative tissue parameters (R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and magnetic susceptibility) were then computed across hemispheres and subjects.

Results

All qualitative findings presented in this study were consistent across all subjects unless otherwise specified.

Figure 1 presents representative MR images and corresponding histological sections of the midbrain. The red nuclei and the substantia nigra were most strikingly discernible on the $R_{\mathrm{2}}^{\mathrm{*}}$ (Figure 1B) and susceptibility maps (Figure 1C) as well as on the direction-encoded track-density (Figure 1D) and T₂-weighted images (Figure 1E), but were barely visible on the R₁ map (Figure 1A). This qualitative finding was also supported by the CNR measurements presented in Figure 2. On the track-density image (Figure 1D) both nuclei showed lower anisotropy as reflected in the reduced color saturation (arrows b and d). Different fiber bundles traversing the crus cerebri [corticobulbar fibers (arrow c₁), corticospinal fibers (arrow c₂), and corticopontine fibers (arrow c₃)] could be identified on the track-density images. In the midbrain, interestingly, the central tegmental tract (arrow f) and the medial lemniscus (arrow e, five of six subjects) were only identifiable reliably on susceptibility maps, whereas the mammillary body (Figure 1; arrow a) was discernible on all MR image contrasts except the track-density images. The superior (not shown) and inferior colliculi (Figure 1; arrow g) were distinguishable across all subjects only on the R₁ maps. Both of these structures could be assessed on the other MRI contrasts in at least four of six subjects (superior colliculus on susceptibility maps in six of six subjects). The superior colliculus and inferior colliculus presented highest CNR on susceptibility maps and $R_{\mathrm{2}}^{\mathrm{*}}$ maps, respectively (Figure 2). Anatomic regions with low magnetic susceptibility (Figure 1C; arrows c, e, f), such as the hypointense rim around the red nuclei coincided with regions of high myelin content (Figure 1G; arrows c_1–3, e, f) and had no unique imaging correlate on any of the other MR images. The histological stains (Figures 1G,H) showed increased myelin and decreased cell density in the red nucleus compared to the substantia nigra.

FIGURE 1

Figure 1. Images of the midbrain. R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps as well as direction-encoded track-density and T₂-weighted images of the same brain region are presented in axial orientation from (A–E), respectively. Red, blue, and green in the track-density image represent anisotropy along medial-lateral, superior-inferior, and anterior-posterior directions, respectively. Note that the track-density image is slightly distorted compared to the other image contrasts (A–C,E). The sectional plane of the axial images is indicated by the dashed line overlaid on the sagittal R₁ map shown in (F). Axial histological sections from a different subject stained for myelin and cells are illustrated in (G) and (H), respectively. The arrows in the images indicate: (a) mammillary body, (b) substantia nigra, (c) crus cerebri, (c₁) corticobulbar fibers, (c₂) corticospinal fibers, (c₃) corticopontine fibers, (d) red nucleus, (e) medial lemniscus, (f) central tegmental tract, (g) inferior colliculus, and (h) medial longitudinal fasciculus. [The histological stains (G,H) were adapted with permission from http://www.brains.rad.msu.edu and http://brainmuseum.org, supported by the US National Science Foundation.]

FIGURE 2

Figure 2. Contrast-to-noise ratio (CNR) of selected anatomical regions. Mean values of CNR across six subjects and both hemispheres were extracted from R₁ (blue bars), $R_{\mathrm{2}}^{\mathrm{*}}$ (red bars), and susceptibility maps (green bars) as well as T₂-weighted images (black bars). Brain regions include: RN - red nucleus vs. surrounding tissue, SN - substantia nigra vs. tissue between red nucleus and substantia nigra, CTT - central tegmental tract vs. surrounding tissue (without red nucleus), sup col - superior colliculus vs. adjacent tissue, inf col - inferior colliculus vs. adjacent tissue, RTN - reticulotegmental nucleus vs. surrounding white matter tissue, MCP - middle cerebellar peduncle vs. adjacent cerebrospinal fluid, SCP - superior cerebellar peduncle vs. reticular formation, TPF - transverse pontine fibers vs. corticospinal tract. The error bars indicate the 95%-confidence interval.

MR images of the rostral and the middle part of the pons are depicted in Figures 3 and 4, respectively. The $R_{\mathrm{2}}^{\mathrm{*}}$ and magnetic susceptibility maps clearly displayed fibers traversing the pons in medial-lateral direction and fibers of the middle cerebellar peduncle (MCP) (arrow a in Figure 3) as also indicated by the myelin stains (Figures 3G and 4G). These fiber structures were barely visible on the R₁ maps [Figures 2 (TPF), 3A, and 4A]. The direction of the transverse pontine fibers could not be reliably separated from the superior-inferior running corticospinal tract on the track-density images (Figures 3D and 4D). In contrast, the susceptibility maps revealed even smaller fiber tracts [e.g., central tegmental tract (arrow f, two of six subjects), medial longitudinal fasciculus (arrow e)] and nerve nuclei [reticulotegmental nucleus (arrow c)] that were indiscernible on R₁ and $R_{\mathrm{2}}^{\mathrm{*}}$ maps as well as on the T₂-weighted images. In the direction-encoded track-density images larger [e.g., MCP (arrow a), superior cerebellar peduncle (arrow d), corticospinal tract (arrow h), pontocerebellar fibers (white solid outline in Figure 3D)] and smaller [e.g., central tegmental tract (arrow f)] fiber tracts were clearly identifiable. The MCP was also discernible on the R₁, $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps, whereas the superior cerebellar peduncle displayed on the R₁ and the susceptibility maps on all subjects and in two of six subjects on the $R_{\mathrm{2}}^{\mathrm{*}}$ maps. Both of these structures (MCP, SCP) were also barely visible on T₂-weighted images (MCP in three subjects, SCP in one subject; Figure 2). The susceptibility maps showed the highest CNR for nearly all investigated structures within the pons (Figure 2). Only the CNR value of MCP with respect to adjacent cerebrospinal fluid were outperformed by the $R_{\mathrm{2}}^{\mathrm{*}}$ maps. Although pontine veins were only seen on both $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps (Figures 4B,C; arrow n), their identification compared to adjacent tissue structures (e.g., transverse pontine fibers) was substantially improved on susceptibility maps. The solid white outline in Figures 4A–D encloses the caudal part of the pontine reticular nucleus, the facial nucleus, salivatory nucleus and nucleus abducens. Although R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps displayed heterogeneous signal intensities in this area, discrimination of these nuclei was not possible. The interruption of the fiber tracts between the MCP seen on the track-density image (dashed outline in Figure 4D) is most likely caused by data acquisition and TDI post-processing issues (incompletely corrected geometric distortions in the unwarped diffusion-weighted images that impeded reliable TDI post-processing) because no such interruption is observed on the susceptibility map (Figure 4C) and on the histological stains (Figures 4G,H). The T₂-weighted images (Figures 3E and 4E) displayed rather homogeneously with low contrast (Figures 2–4), which prevented reliable identification of any substructures across all subjects.

FIGURE 3

Figure 3. Images of the rostral part of the pons. R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps as well as direction-encoded track-density (same color-encoding as in Figure 1; slightly distorted compared to the other axial image contrasts) and T₂-weighted images of the same region are presented in axial orientation from (A–E), respectively. The sectional plane of the axial images is indicated by the dashed line overlaid on the sagittal R₁ map shown in (F). Axial histological sections from a different subject stained for myelin and cells are illustrated in (G) and (H), respectively. The arrows in the images indicate: (a) middle cerebellar peduncle, (b) medial lemniscus, (c) reticulotegmental nucleus, (d) superior cerebellar peduncle, (e) medial longitudinal fasciculus, (f) central tegmental tract, and (h) corticospinal tract. The solid white outline in (D) indicates the pontocerebellar fibers. [The histological stains (G,H) were adapted with permission from http://www.brains.rad.msu.edu and http://brainmuseum.org, supported by the US National Science Foundation.]

FIGURE 4

Figure 4. Images of the middle part of the pons. R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps as well as track-density (same color-encoding as in Figure 1; slightly distorted compared to the other axial image contrasts) and T₂-weighted images of the same region are presented in axial orientation from (A–E), respectively. The sectional plane of the axial images is indicated by the dashed line overlaid on the sagittal R₁ map shown in (F). Axial histological sections from a different subject stained for myelin and cells are illustrated in (G) and (H), respectively. The arrows in the images indicate: (a) middle cerebellar peduncle, (d) superior cerebellar peduncle, (e) medial longitudinal fasciculus, (h) corticospinal fibers, (i) pontocerebellar tract, (j) facial nucleus, (k) salivatory nucleus, (l) abducens nucleus, (m) solitary nucleus, (n) vein, (o) fiber bundle that traverses the pons in medial to lateral direction, and (p) fourth ventricle (cerebrospinal fluid). The solid white outline in (A–D) summarizes the caudal part of the pontine reticular nucleus, the facial nucleus, salivatory nucleus, and abducens nucleus. The dashed white outline in (D) indicates a region between the middle cerebellar peduncle that could not be resolved with TDI. [The histological stains (G,H) were adapted with permission from http://www.brains.rad.msu.edu and http://brainmuseum.org, supported by the US National Science Foundation.]

Figures 5 and 6 display the rostral and middle part of the medulla oblongata, respectively. The pyramid (arrow a) is located anterior to the inferior olive and was clearly visible on all contrasts (four of six subjects on $R_{\mathrm{2}}^{\mathrm{*}}$ maps) except the T₂-weighted images. It could be well discriminated from the inferior olive (arrow b) which was also clearly discernible on the R₁ (five subjects) and susceptibility map (four subjects), and less obvious on $R_{\mathrm{2}}^{\mathrm{*}}$ maps (three subjects). The inferior cerebellar peduncle was clearly visible on track-density images and was distinguishable on R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps in four subjects at least. Only the susceptibility maps enabled identification of the spinal trigeminal nucleus (arrow c in Figure 5; four subjects). Although the histological stains (Figures 5G,H) suggested the presence of several cranial nuclei, further nuclei were not visible on the MR images.

FIGURE 5

Figure 5. Images of the rostral part of the medulla oblongata. R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps as well as direction-encoded track-density (same color-encoding as in Figure 1; slightly distorted compared to the other axial image contrasts) and T₂-weighted images of the same region are presented in axial orientation from (A–E), respectively. The sectional plane of the axial images is indicated by the dashed line overlaid on the sagittal R₁ map shown in (F). Axial histological sections from a different subject stained for myelin and cells are illustrated in (G) and (H), respectively. The arrows indicate: (a) pyramid, (b) inferior olive, (c) spinal trigeminal nucleus, (d) inferior cerebellar peduncle, (e) medial longitudinal fasciculus, (f) medial lemniscus, (h) reticular formation, (i) accessory cuneate nucleus, (j) solitary nucleus, (k) dorsal motor nucleus of the vagus nerve, and (l) hypoglossal nucleus. [The histological stains (G,H) were adapted with permission from http://www.brains.rad.msu.edu and http://brainmuseum.org, supported by the US National Science Foundation.]

FIGURE 6

Figure 6. Images of the middle part of the medulla oblongata. R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility maps as well as track-density (same color-encoding as in Figure 1; slightly distorted compared to the other axial image contrasts) and T₂-weighted images of the same region are presented in axial orientation from (A–E), respectively. The sectional plane of the axial images is indicated by the dashed line overlaid on the sagittal R₁ map shown in (F). The arrows in the images indicate: (a) pyramid, (b) inferior olive, and (d) inferior cerebellar peduncle.

Finally, R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and susceptibility differences are listed in Table 1 for selected anatomical regions. The differences of the susceptibility and the $R_{\mathrm{2}}^{\mathrm{*}}$ values between fibers running nearly parallel (transverse pontine fibers) and nearly perpendicular (corticospinal tract) to the static magnetic field were (−49 ± 11) ppb and (8.9 ± 6) s⁻¹, respectively.

TABLE 1

Table 1. Mean values of longitudinal relaxation rate (R₁), effective transverse relaxation rate ($R_{\mathrm{2}}^{\mathrm{*}}$), and volume magnetic susceptibility difference (Δχ) with respect to normal appearing white matter over six subjects and both brain hemispheres for selected anatomical regions.

Discussion

We have investigated the anatomy of the human brainstem in vivo by applying MR imaging at 7 T with high spatial, isotropic resolution, and using multiple image contrasts. QSM and TDI have been applied jointly for the first time to identify anatomical substructures of the brainstem. Both MR methods directly reflect subtle variations in tissue composition that were found to be consistent with histology as demonstrated with corresponding histological stains selected from a data base. In contrast, longitudinal relaxation rate mapping and T₂-weighted imaging both performed inferiorly in their ability to delineate anatomical details (Figure 2).

$R_{\mathrm{2}}^{\mathrm{*}}$ relaxation comprises both intrinsic microscopic transverse relaxation (R₂) and relaxation due to heterogeneous magnetic susceptibility distributions on a mesoscopic or macroscopic scale $\left(R_{\mathrm{2}}^{\prime }\right)$. In contrast to the $R_{\mathrm{2}}^{\mathrm{*}}$ relaxation rate, which increases with the concentration of both iron (Langkammer et al., 2010) and myelin (Lee et al., 2012), magnetic susceptibility shows a different dependence on the concentration of these substances, i.e., higher magnetic susceptibility values with increasing iron concentration (Langkammer et al., 2012b), and lower magnetic susceptibility values with increasing myelin density (Liu et al., 2011; Schweser et al., 2011; Li et al., 2012a,b). Hence, image contrast on $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps appears complementary. The longitudinal relaxation rate, R₁, in turn is heavily influenced by myelin content, but shows little effect with respect to tissue iron concentration (Bridge and Clare, 2006).

Both $R_{\mathrm{2}}^{\mathrm{*}}$ maps and susceptibility maps provided detailed depiction of the non-cranial nuclei and nerve fibers in the midbrain in accordance with recent studies (Schweser et al., 2011; Deistung et al., 2013a). In particular, the red nuclei and the substantia nigra were clearly discernible on susceptibility and $R_{\mathrm{2}}^{\mathrm{*}}$ maps because both structures contain substantially higher iron concentrations compared to the surrounding brainstem tissue (Morris et al., 1992). Interestingly, the highly myelinated rim around the red nuclei, which was clearly seen on the myelin stain and on the susceptibility map was not visible on the R₁ map, potentially due to the lower spatial resolution of the MP2RAGE acquisition compared to the gradient-echo acquisition.

Brainstem nuclei were barely seen on the R₁, $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps sectioning the pons and medulla oblongata; however, the inferior olive (Figures 5 and 6, arrow b) could be identified and both the spinal trigeminal nucleus (Figure 5C, arrow c) and the reticulotegmental nucleus (Figure 3C, arrow c) were discernible on the susceptibility maps. The inferior olive is the sole source of the climbing fibers to the Purkinje cells in the cerebellar cortex and is involved in learning and timing of motor behavior (De Zeeuw et al., 1998). Lesions in the inferior olive can thus introduce restrictions in performing specialized motor tasks (Martin et al., 1996). The spinal trigeminal nucleus is the largest of all cranial nerve nuclei and one of the three paired sensory nuclei associated with the trigeminal nerve. It receives input from pain, temperature, and some tactile afferents in the trigeminal nerve and is involved in migraine headache (Kaube et al., 2002). The reticulotegmental nucleus of the pons is an important component of the oculomotor circuit that regulates horizontal saccades and smooth pursuit movements of the eyes (Keller and Crandall, 1981; Büttner-Ennever and Horn, 1997). Damage to the reticulotegmental nucleus has been observed in spinocerebellar ataxia (Rüb et al., 2004) and it most probably contributes to hypometric horizontal saccades and slowing of smooth pursuits that characteristically develop in patients suffering from Alzheimer’s disease (Rüb et al., 2001). Therefore, non-invasive assessment of these nuclei may help to improve diagnosis and understanding of these diseases.

These nuclei (inferior olive, spinal trigeminal nucleus, reticulotegmental nucleus) were visible on the susceptibility maps because of their different iron concentration and/or myelin content with respect to surrounding tissue. A histological study by Morris et al. (1992), for instance, demonstrated iron accumulation in the nuclei of the dorsal parts of the brainstem, the pyramids, the spinal trigeminal nucleus, and the inferior olive (although with substantially lower iron concentration than in the red nuclei and substantia nigra) and reported only little or no intrinsic iron reactivity for many other nuclei and tracts of the brainstem (e.g., solitary nucleus, hypoglossal nucleus, raphe pallidus nucleus, cuneate fasciculus, and gracile fasciculus). Brainstem nuclei have a similar, relatively low myelin density as seen from the histological myelin stains. Only the inferior olive (Figure 5G, arrow b) and the spinal trigeminal nucleus (Figure 5G, arrow c) and reticulotegmental nucleus [see (Paxinos et al., 2012) and Figure 3C, arrow c] are embedded in tissue with substantially increased myelin density.

In addition to myelin’s bulk diamagnetic nature, a dependence of the magnetic susceptibility on the axons’ alignment with respect to the external magnetic field has been recently reported (Liu, 2010; Li et al., 2012a; Schweser et al., 2012; Wharton and Bowtell, 2012). Similar direction dependent observations have been made for $R_{\mathrm{2}}^{\mathrm{*}}$ (Bender and Klose, 2010; Denk et al., 2011; Lee et al., 2011; Sati et al., 2012). We measured differences of the susceptibility and the $R_{\mathrm{2}}^{\mathrm{*}}$ values between fibers running nearly parallel (transverse pontine fibers) and nearly perpendicular (corticospinal tract) to the static magnetic field of (−49 ± 11) ppb and (8.9 ± 6) s⁻¹, respectively. This susceptibility difference is nearly threefold as large as the values that have been derived for a hollow cylinder model of myelin (−18 to −16 ppb) (Li et al., 2012a; Wharton and Bowtell, 2012). In vivo, however, Li et al., 2012a also reported susceptibility differences between parallel and perpendicular fibers of up to −40 ppb in individual white matter regions, a value which is almost the same as reported here. The result for the orientation dependence of $R_{\mathrm{2}}^{\mathrm{*}}$ is in good agreement with the recently reported ex vivo value of −6.44 ± 0.15 s⁻¹ at 7 T (Lee et al., 2011). The orientation-dependent differences of magnetic susceptibility and $R_{\mathrm{2}}^{\mathrm{*}}$in our study may be caused by the presence of additional contributors other than myelin, such as pontine nuclei and neuropil within the investigated volumes-of-interest. Nevertheless, consistent with the theory of myelin’s orientation dependency, the contrast between fibers running in medial-lateral direction (oriented nearly perpendicular to the magnetic field; e.g., transverse pontine fibers, arrow o in Figure 4) and fibers running inferior-superiorly (oriented nearly parallel to the magnetic field; corticospinal tract) is enhanced on the $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps. More sophisticated approaches, such as susceptibility tensor imaging (STI) (Liu, 2010; Li et al., 2012a) or multipole susceptibility tensor imaging (MSTI) (Liu and Li, 2013) may further improve discrimination between fibers and nuclei of the pons and medulla, particularly in combination with high spatial resolution.

Direction-encoded TDI enabled discrimination of the directionality of nerve fibers within the brainstem and yielded complementary information to the relaxation and susceptibility maps. The depiction of larger fiber bundles was consistent with the results of other DTI studies conducted in vivo at lower magnetic fields (Nagae-Poetscher et al., 2004; Salamon et al., 2005). We were, however, unable to resolve unambiguously the transverse pontine fibers, which were clearly visible on the $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility maps (Figure 4). This may be due to the lower base voxel resolution of 1.2 mm of the DWI data, which was chosen as a compromise between spatial resolution, SNR, and acquisition time. To overcome this limitation we interpolated the DWI data using a super-resolution approach (TDI) to obtain an isotropic virtual resolution of 0.43 mm. It should be noted, however, that the TDI technique is not able to fully recover the whole information that would be present if the data would be acquired with such a high spatial resolution (Calamante et al., 2011). Due to the larger errors in tractography by performing fiber tracking on data with a lower resolution, the higher resolved track-density images appeared slightly blurred and lost some fine details; explaining the difficulties in identifying the transverse pontine fibers in the presented track-density images (Figures 3D and 4D). As might be expected, the level of blurring in the track-density images depends on the desired virtual resolution (Calamante et al., 2011). Higher spatial resolution is thus needed to delineate transverse pontine fibers unambiguously, which requires sophisticated modifications of the diffusion MRI sequence such as combining zoomed imaging and parallel imaging (Heidemann et al., 2012; Eichner et al., 2013). Karampinos et al. (2009) were recently able to resolve transverse pontine fibers in direction-encoded fractional anisotropy maps with a resolution of 0.8 mm × 0.8 mm × 3 mm at 3T by applying a dedicated acquisition scheme (self-navigated, multi-shot, variable-density, spiral-imaging with outer volume suppression).

One particular problem of brainstem imaging with MRI are the induced field inhomogeneities caused by nearby air-tissue or bone-tissue interfaces. In the vicinity of such interfaces significant image distortions of EPI-based diffusion-weighted images occur and additional spin dephasing arises in gradient-echo images. Despite unwarping the diffusion-weighted images based on gradient-echo field-map information, it was not possible to completely compensate geometric distortions, thus impeding accurate depiction of ventral fiber structures of the middle pons and medulla oblongata on track-density images. These induced field inhomogeneities manifest themselves as steep gradients in gradient-echo phase images and lead to signal reductions in gradient-echo magnitude images and decreased SNRs, affecting both susceptibility maps [due to increased phase noise (Gudbjartsson and Patz, 1995)] and $R_{\mathrm{2}}^{\mathrm{*}}$ maps. In future studies, the impact of spin dephasing due to air-tissue or bone-tissue interfaces on $R_{\mathrm{2}}^{\mathrm{*}}$ maps may be minimized by taking into account the macroscopic field inhomogeneity reflected in the phase images (de Leeuw and Bakker, 2012). R₁ maps and T₂-weighted images, on the other hand, are distinctly more immune against these field inhomogeneity effects.

The quantitative analysis revealed similar values for the longitudinal relaxation rates (R₁) for nuclei and myelinated fibers (Table 1), whereas susceptibility and $R_{\mathrm{2}}^{\mathrm{*}}$ values varied more substantially. This broader spread of $R_{\mathrm{2}}^{\mathrm{*}}$ and susceptibility values is reflected in the large CNR values (Figure 2) underlining the good discrimination of nerve fibers on susceptibility maps and the exquisite delineation of midbrain nuclei on susceptibility and $R_{\mathrm{2}}^{\mathrm{*}}$ maps. It should, however, be noted that the discrimination of brain tissue with respect to cerebral spinal fluid (CSF) on susceptibility maps is not as striking as on the R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and T₂-weighted images [Table 1; Figure 4 (arrow p)]. Due to the quantitative nature of R₁, $R_{\mathrm{2}}^{\mathrm{*}}$, and magnetic susceptibility maps, these contrasts may be combined by projections along support vectors observed from discriminant analysis or support vector machines to generate composite images with improved depiction of anatomical features while providing improved discrimination of CSF (Deistung et al., 2013b).

In conclusion, maps of magnetic susceptibility displayed the largest number of brainstem structures, including larger and smaller fiber pathways as well as several nerve nuclei. Usage of multiple image contrasts enables a detailed non-invasive view into tissue structure and composition. Hence, multi-contrast MR imaging of the brainstem at ultra-high magnetic fields that utilizes relaxation, diffusion, and magnetic susceptibility information is a versatile tool to assess anatomy individually in great detail.

Conflict of Interest Statement

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.

Acknowledgments

We would like to thank Andrew Webb (Gorter Center for High Field MRI, Department of Radiology, Leiden University Medical Center) for providing the dielectric pads as well as Penny Gowland and Andrew Peters (Sir Peter Mansfield Magnetic Resonance Centre, School of Physics and Astronomy, University of Nottingham) for providing the tailored RF pulse. This research was funded by a German Research Foundation research grant (DFG, RE 1123/9-2), by the Max-Planck society, and by seed grants awarded by the Interdisciplinary Center for Clinical Research (IZKF) in Jena (to Andreas Deistung), by the International Society of Magnetic Resonance in Medicine (to Ferdinand Schweser), and by the Friedrich Schiller University Jena (to Ferdinand Schweser).

Abbreviations

DTI, diffusion tensor imaging; GRE, gradient recalled echo; QSM, quantitative susceptibility mapping; ppb, parts per billion; $R_{\mathrm{2}}^{\mathrm{*}}$, effective transverse relaxation rate; SHARP, sophisticated harmonic artifact reduction for phase data; TDI, track-density imaging.

Footnotes

1. ^http://surfer.nmr.mgh.harvard.edu/
2. ^http://www.brains.rad.msu.edu
3. ^http://brainmuseum.org

References

Basser, P. J., and Pierpaoli, C. (1996). Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J. Magn. Reson. B. 111, 209–219. doi: 10.1006/jmrb.1996.0086

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bender, B., and Klose, U. (2010). The in vivo influence of white matter fiber orientation towards B(0) on T2* in the human brain. NMR Biomed. 23, 1071–1076. doi:10.1002/nbm.1534

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bergman, R. A., Afifi, A. K., and Heidger, P. M. (eds). (1989). “Atlas of microscopic anatomy: companion to histology and neuroanatomy,” in Central Nervous System, 2nd Edn (Philadelphia: W. B. Saunders Company), 418–495.

Bilgic, B., Pfefferbaum, A., Rohlfing, T., Sullivan, E. V., and Adalsteinsson, E. (2012). MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. Neuroimage 59, 2625–2635. doi:10.1016/j.neuroimage.2011.08.077

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bridge, H., and Clare, S. (2006). High-resolution MRI: in vivo histology? Philos. Trans. R. Soc. Lond. B Biol. Sci. 361, 137–146. doi:10.1098/rstb.2005.1777

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Büttner-Ennever, J. A., and Horn, A. K. (1997). Anatomical substrates of oculomotor control. Curr. Opin. Neurobiol. 7, 872–879. doi:10.1016/S0959-4388(97)80149-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Calamante, F., Oh, S.-H., Tournier, J.-D., Park, S.-Y., Son, Y.-D., Chung, J.-Y., et al. (2013). Super-resolution track-density imaging of thalamic substructures: comparison with high-resolution anatomical magnetic resonance imaging at 7.0T. Hum. Brain Mapp. 34, 2538–2548. doi:10.1002/hbm.22083

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Calamante, F., Tournier, J.-D., Heidemann, R. M., Anwander, A., Jackson, G. D., and Connelly, A. (2011). Track density imaging (TDI): validation of super resolution property. Neuroimage 56, 1259–1266. doi:10.1016/j.neuroimage.2011.02.059

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Calamante, F., Tournier, J.-D., Jackson, G. D., and Connelly, A. (2010). Track-density imaging (TDI): super-resolution white matter imaging using whole-brain track-density mapping. Neuroimage 53, 1233–1243. doi:10.1016/j.neuroimage.2010.07.024

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

de Leeuw, H., and Bakker, C. J. G. (2012). Correction of gradient echo images for first and second order macroscopic signal dephasing using phase derivative mapping. Neuroimage 60, 818–829. doi:10.1016/j.neuroimage.2011.11.083

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

De Zeeuw, C. I., Simpson, J. I., Hoogenraad, C. C., Galjart, N., Koekkoek, S. K., and Ruigrok, T. J. (1998). Microcircuitry and function of the inferior olive. Trends Neurosci. 21, 391–400. doi:10.1016/S0166-2236(98)01310-1

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Deistung, A., Rauscher, A., Sedlacik, J., Stadler, J., Witoszynskyj, S., and Reichenbach, J. R. (2008). Susceptibility weighted imaging at ultra high magnetic field strengths: theoretical considerations and experimental results. Magn. Reson. Med. 60, 1155–1168. doi:10.1002/mrm.21754

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Deistung, A., Schäfer, A., Schweser, F., Biedermann, U., Turner, R., and Reichenbach, J. R. (2013a). Toward in vivo histology: a comparison of quantitative susceptibility mapping (QSM) with magnitude-, phase-, and R2*-imaging at ultra-high magnetic field strength. Neuroimage 65, 299–314. doi:10.1016/j.neuroimage.2012.09.055

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Deistung, A., Schweser, F., and Reichenbach, J. R. (2013b). “Optimal enhancement of brain structures by combining different MR contrasts: demonstration of venous vessel enhancement in multi-echo gradient-echo MRI,” in Proc. 21st Annual Meeting of the Intl. Soc. Mag. Reson. Med (Salt Lake City, UT), 695.

Denk, C., Torres, E. H., MacKay, A., and Rauscher, A. (2011). The influence of white matter fibre orientation on MR signal phase and decay. NMR Biomed. 24, 246–252. doi:10.1002/nbm.1581

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Donaldson, S. S., Laningham, F., and Fisher, P. G. (2006). Advances toward an understanding of brainstem gliomas. J. Clin. Oncol. 24, 1266–1272. doi:10.1200/JCO.2005.04.6599

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Duyn, J. H., van Gelderen, P., Li, T. Q., de Zwart, J. A., Koretsky, A. P., and Fukunaga, M. (2007). From the cover: high-field MRI of brain cortical substructure based on signal phase. Proc. Natl. Acad. Sci. U.S.A. 104, 11796–11801. doi:10.1073/pnas.0610821104

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Eichner, C., Setsompop, K., Koopmans, P. J., Lützkendorf, R., Norris, D. G., Turner, R., et al. (2013). Slice accelerated diffusion-weighted imaging at ultra-high field strength. Magn. Reson. Med. doi:10.1002/mrm.24809. [Epub ahead of print]

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Feinberg, D. A., and Oshio, K. (1991). GRASE (gradient- and spin-echo) MR imaging: a new fast clinical imaging technique. Radiology 181, 597–602.

Pubmed Abstract | Pubmed Full Text

Fukunaga, M., Li, T.-Q., van Gelderen, P., de Zwart, J. A., Shmueli, K., Yao, B., et al. (2010). Layer-specific variation of iron content in cerebral cortex as a source of MRI contrast. Proc. Natl. Acad. Sci. U.S.A. 107, 3834–3839. doi:10.1073/pnas.0911177107

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Geyer, S., Weiss, M., Reimann, K., Lohmann, G., and Turner, R. (2011). Microstructural parcellation of the human cerebral cortex – from Brodmann’s post-mortem map to in vivo mapping with high-field magnetic resonance imaging. Front. Hum. Neurosci. 5:19. doi:10.3389/fnhum.2011.00019

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Griswold, M. A., Jakob, P. M., Heidemann, R. M., Nittka, M., Jellus, V., Wang, J., et al. (2002). Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn. Reson. Med. 47, 1202–1210. doi:10.1002/mrm.10171

CrossRef Full Text

Gudbjartsson, H., and Patz, S. (1995). The Rician distribution of noisy MRI data. Magn. Res. Med. 34, 910–914. doi:10.1002/mrm.1910340618

CrossRef Full Text

Haacke, E. M., Tang, J., Neelavalli, J., and Cheng, Y. C. N. (2010). Susceptibility mapping as a means to visualize veins and quantify oxygen saturation. J. Magn. Reson. Imaging 32, 663–676. doi:10.1002/jmri.22276

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Heidemann, R. M., Anwander, A., Feiweier, T., Knösche, T. R., and Turner, R. (2012). k-space and q-space: combining ultra-high spatial and angular resolution in diffusion imaging using ZOOPPA at 7 T. Neuroimage 60, 967–978. doi:10.1016/j.neuroimage.2011.12.081

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hurley, A. C., Al-Radaideh, A., Bai, L., Aickelin, U., Coxon, R., Glover, P., et al. (2010). Tailored RF pulse for magnetization inversion at ultrahigh field. Magn. Reson. Med. 63, 51–58. doi:10.1002/mrm.22167

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Jones, D. K., Knösche, T. R., and Turner, R. (2013). White matter integrity, fiber count, and other fallacies: the do’s and don’ts of diffusion MRI. Neuroimage 73, 239–254. doi:10.1016/j.neuroimage.2012.06.081

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Karampinos, D. C., Van, A. T., Olivero, W. C., Georgiadis, J. G., and Sutton, B. P. (2009). High-resolution diffusion tensor imaging of the human pons with a reduced field-of-view, multishot, variable-density, spiral acquisition at 3T. Magn. Reson. Med. 62, 1007–1016. doi:10.1002/mrm.22105

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kaube, H., Katsarava, Z., Przywara, S., Drepper, J., Ellrich, J., and Diener, H.-C. (2002). Acute migraine headache: possible sensitization of neurons in the spinal trigeminal nucleus? Neurology 58, 1234–1238. doi:10.1212/WNL.58.8.1234

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Keller, E. L., and Crandall, W. F. (1981). Neural activity in the nucleus reticularis tegmenti pontis in the monkey related to eye movements and visual stimulation. Ann. N. Y. Acad. Sci. 374, 249–261. doi:10.1111/j.1749-6632.1981.tb30874.x

CrossRef Full Text

Koenig, S. H. (1991). Cholesterol of myelin is the determinant of gray-white contrast in MRI of brain. Magn. Reson. Med. 20, 285–291. doi:10.1002/mrm.1910200210

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Langkammer, C., Krebs, N., Goessler, W., Scheurer, E., Ebner, F., Yen, K., et al. (2010). Quantitative MR imaging of brain iron: a postmortem validation study. Radiology 257, 455–462. doi:10.1148/radiol.10100495

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Langkammer, C., Krebs, N., Goessler, W., Scheurer, E., Yen, K., Fazekas, F., et al. (2012a). Susceptibility induced gray-white matter MRI contrast in the human brain. Neuroimage 59, 1413–1419. doi:10.1016/j.neuroimage.2011.08.045

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Langkammer, C., Schweser, F., Krebs, N., Deistung, A., Goessler, W., Scheurer, E., et al. (2012b). Quantitative susceptibility mapping (QSM) as a means to measure brain iron? A post mortem validation study. Neuroimage 62, 1593–1599. doi:10.1016/j.neuroimage.2012.05.049

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lebel, R. M., Eissa, A., Seres, P., Blevins, G., and Wilman, A. H. (2012). Quantitative high-field imaging of sub-cortical gray matter in multiple sclerosis. Mult. Scler. 18, 433–441. doi:10.1177/1352458511428464

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lee, J., Shmueli, K., Kang, B.-T., Yao, B., Fukunaga, M., van Gelderen, P., et al. (2012). The contribution of myelin to magnetic susceptibility-weighted contrasts in high-field MRI of the brain. Neuroimage 59, 3967–3975. doi:10.1016/j.neuroimage.2011.10.076

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lee, J., van Gelderen, P., Kuo, L.-W., Merkle, H., Silva, A. C., and Duyn, J. H. (2011). T2*-based fiber orientation mapping. Neuroimage 57, 225–234. doi:10.1016/j.neuroimage.2011.04.026

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Li, T. Q., van Gelderen, P., Merkle, H., Talagala, L., Koretsky, A. P., and Duyn, J. (2006). Extensive heterogeneity in white matter intensity in high-resolution T2*-weighted MRI of the human brain at 7.0 T. Neuroimage 32, 1032–1040. doi:10.1016/j.neuroimage.2006.05.053

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Li, W., Wu, B., Avram, A. V., and Liu, C. (2012a). Magnetic susceptibility anisotropy of human brain in vivo and its molecular underpinnings. Neuroimage 59, 2088–2097. doi:10.1016/j.neuroimage.2011.10.038

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Li, X., Vikram, D. S. I, Lim, A. L., Jones, C. K., Farrell, J. A. D., and van Zijl, P. C. M. (2012b). Mapping magnetic susceptibility anisotropies of white matter in vivo in the human brain at 7T. Neuroimage 62, 314–330. doi:10.1016/j.neuroimage.2012.04.042

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Li, W., Wu, B., and Liu, C. (2011). Quantitative susceptibility mapping of human brain reflects spatial variation in tissue composition. Neuroimage 55, 1645–1656. doi:10.1016/j.neuroimage.2010.11.088

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Liu, C. (2010). Susceptibility tensor imaging. Magn. Reson. Med. 63, 1471–1477. doi:10.1002/mrm.22482

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Liu, C., and Li, W. (2013). Imaging neural architecture of the brain based on its multipole magnetic response. Neuroimage 67, 193–202. doi:10.1016/j.neuroimage.2012.10.050

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Liu, C., Li, W., Johnson, G. A., and Wu, B. (2011). High-field (9.4 T) MRI of brain dysmyelination by quantitative mapping of magnetic susceptibility. Neuroimage 56, 930–938. doi:10.1016/j.neuroimage.2011.02.024

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Marques, J. P., and Gruetter, R. (2013). New developments and applications of the MP2RAGE sequence – focusing the contrast and high spatial resolution R1 mapping. PLoS ONE 8:e69294. doi:10.1371/journal.pone.0069294

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Marques, J. P., Kober, T., Krueger, G., van der Zwaag, W., de Moortele, P.-F. V., and Gruetter, R. (2010). MP2RAGE, a self bias-field corrected sequence for improved segmentation and T1-mapping at high field. Neuroimage 49, 1271–1281. doi:10.1016/j.neuroimage.2009.10.002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Martin, T. A., Keating, J. G., Goodkin, H. P., Bastian, A. J., and Thach, W. T. (1996). Throwing while looking through prisms. I. Focal olivocerebellar lesions impair adaptation. Brain 119(Pt 4), 1183–1198. doi:10.1093/brain/119.4.1183

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

McGibney, G., and Smith, M. R. (1993). An unbiased signal-to-noise ratio measure for magnetic resonance images. Med. Phys. 20, 1077–1078. doi:10.1118/1.597004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Miller, A. J., and Joseph, P. M. (1993). The use of power images to perform quantitative analysis on low SNR MR images. Magn. Reson. Imaging 11, 1051–1056. doi:10.1016/0730-725X(93)90225-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Morris, C. M., Candy, J. M., Oakley, A. E., Bloxham, C. A., and Edwardson, J. A. (1992). Histochemical distribution of non-haem iron in the human brain. Acta Anat. (Basel) 144, 235–257. doi:10.1159/000147312

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Moser, E., Stahlberg, F., Ladd, M. E., and Trattnig, S. (2012). 7-T MR–from research to clinical applications? NMR Biomed. 25, 695–716. doi:10.1002/nbm.1794

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Nagae-Poetscher, L. M., Jiang, H., Wakana, S., Golay, X., van Zijl, P. C. M., and Mori, S. (2004). High-resolution diffusion tensor imaging of the brain stem at 3 T. AJNR Am. J. Neuroradiol. 25, 1325–1330.

Pubmed Abstract | Pubmed Full Text

Naidich, T. P., Duvernoy, H. M., Delman, B., Sorensen, A., Kollias, S., and Haacke, E. (2009). Duvernoy’s Atlas of the Human Brain Stem and Cerebellum: High-Field MRI, Surface Anatomy, Internal Structure, Vascularization and 3D Sectional Anatomy. New York: Springer.

Nölte, I. S., Gerigk, L., Al-Zghloul, M., Groden, C., and Kerl, H. U. (2012). Visualization of the internal globus pallidus: sequence and orientation for deep brain stimulation using a standard installation protocol at 3.0 Tesla. Acta Neurochir. (Wien) 154, 481–494. doi:10.1007/s00701-011-1242-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ogg, R. J., and Steen, R. G. (1998). Age-related changes in brain T1 are correlated with iron concentration. Magn. Reson. Med. 40, 749–753. doi:10.1002/mrm.1910400516

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Oshio, K., and Feinberg, D. A. (1991). GRASE (Gradient- and spin-echo) imaging: a novel fast MRI technique. Magn. Reson. Med. 20, 344–349. doi:10.1002/mrm.1910200219

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Paxinos, G., Huang, X., Sengul, G., and Watson, C. (2012). The Human Nervous System, Chapter Organization of Brainstem Nuclei. Amsterdam: Elsevier Academic Press, 260–327.

Reichenbach, J. R. (2012). The future of susceptibility contrast for assessment of anatomy and function. Neuroimage 62, 1311–1315. doi:10.1016/j.neuroimage.2012.01.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Robinson, S., Grabner, G., Witoszynskyj, S., and Trattnig, S. (2011). Combining phase images from multi-channel RF coils using 3D phase offset maps derived from a dual-echo scan. Magn. Reson. Med. 65, 1638–1648. doi:10.1002/mrm.22753

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Roemer, P. B., Edelstein, W. A., Hayes, C. E., Souza, S. P., and Mueller, O. M. (1990). The NMR phased array. Magn. Reson. Med. 16, 192–225. doi:10.1002/mrm.1910160203

CrossRef Full Text

Rüb, U., Bürk, K., Schöls, L., Brunt, E. R., de Vos, R. A. I., Diaz, G. O., et al. (2004). Damage to the reticulotegmental nucleus of the pons in spinocerebellar ataxia type 1, 2, and 3. Neurology 63, 1258–1263. doi:10.1212/01.WNL.0000140498.24112.8C

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Rüb, U., Schultz, C., Del Tredici, K., and Braak, H. (2001). Early involvement of the tegmentopontine reticular nucleus during the evolution of Alzheimer’s disease-related cytoskeletal pathology. Brain Res. 908, 107–112. doi:10.1016/S0006-8993(01)02598-7

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Salamon, N., Sicotte, N., Alger, J., Shattuck, D., Perlman, S., Sinha, U., et al. (2005). Analysis of the brain-stem white-matter tracts with diffusion tensor imaging. Neuroradiology 47, 895–902. doi:10.1007/s00234-005-1439-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sati, P., Silva, A. C., van Gelderen, P., Gaitan, M. I., Wohler, J. E., Jacobson, S., et al. (2012). In vivo quantification of T2* anisotropy in white matter fibers in marmoset monkeys. Neuroimage 59, 979–985. doi:10.1016/j.neuroimage.2011.08.064

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schofield, M. A., and Zhu, Y. (2003). Fast phase unwrapping algorithm for interferometric applications. Opt. Lett. 28, 1194–1196. doi:10.1364/OL.28.001194

CrossRef Full Text

Schweser, F., Deistung, A., Lehr, B. W., and Reichenbach, J. R. (2011). Quantitative imaging of intrinsic magnetic tissue properties using MRI signal phase: an approach to in vivo brain iron metabolism? Neuroimage 54, 2789–2807. doi:10.1016/j.neuroimage.2010.10.070

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schweser, F., Sommer, K., Deistung, A., and Reichenbach, J. R. (2012). Quantitative susceptibility mapping for investigating subtle susceptibility variations in the human brain. Neuroimage 62, 2083–2100. doi:10.1016/j.neuroimage.2012.05.067

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sigalovsky, I. S., Fischl, B., and Melcher, J. R. (2006). Mapping an intrinsic MR property of gray matter in auditory cortex of living humans: a possible marker for primary cortex and hemispheric differences. Neuroimage 32, 1524–1537. doi:10.1016/j.neuroimage.2006.05.023

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Stieltjes, B., Kaufmann, W. E., van Zijl, P. C., Fredericksen, K., Pearlson, G. D., Solaiyappan, M., et al. (2001). Diffusion tensor imaging and axonal tracking in the human brainstem. Neuroimage 14, 723–735. doi:10.1006/nimg.2001.0861

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sutin, J., and Carpenter, M. B. (1976). Human Neuroanatomy. Baltimore: Lippincott Williams & Wilkins.

Teeuwisse, W. M., Brink, W. M., and Webb, A. G. (2012). Quantitative assessment of the effects of high-permittivity pads in 7 Tesla MRI of the brain. Magn. Reson. Med. 67, 1285–1293. doi:10.1002/mrm.23108

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tofts, P. (ed.). (2003). Quantitative MRI of the Brain: Measuring Changes Caused by Disease. Chichester, UK: John Wiley & Sons, Ltd. doi:10.1002/0470869526

CrossRef Full Text

Tournier, J.-D., Calamante, F., and Connelly, A. (2007). Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage 35, 1459–1472. doi:10.1016/j.neuroimage.2007.02.016

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tournier, J.-D., Calamante, F., and Connelly, A. (2010). “Improved probabilistic streamlines tractography by 2nd order integration over fibre orientation distributions,” in Proc. 18th Annual Meeting of the Intl. Soc. Mag. Reson. Med (Stockholm), 1670.

Tournier, J.-D., Calamante, F., Gadian, D. G., and Connelly, A. (2004). Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution. Neuroimage 23, 1176–1185. doi:10.1016/j.neuroimage.2004.07.037

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Trampel, R., Ott, D. V. M., and Turner, R. (2011). Do the congenitally blind have a stria of Gennari? First intracortical insights in vivo. Cereb. Cortex 21, 2075–2081. doi:10.1093/cercor/bhq282

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Turner, R. (2012). “Neuroscientific applications of high-field MRI in humans,” in High-Field MR Imaging, eds J. Hennig and O. Speck (Berlin: Springer), 137–149.

Turner, R. (2013). Microstructural Parcellation of the Human Cerebral Cortex, Chapter MRI Methods for In-Vivo Cortical Parcellation. Berlin: Springer, 197–220.

Urban, P., and Caplan, L. R. (eds) (2011). Brainstem Disorders. Berlin: Springer.

van der Weerd, L., Vergeldt, F. J., de Jager, P. A., and As, H. V. (2000). Evaluation of algorithms for analysis of NMR relaxation decay curves. Magn. Reson. Imaging 18, 1151–1158. doi:10.1016/S0730-725X(00)00200-9

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wharton, S., and Bowtell, R. (2012). Fiber orientation-dependent white matter contrast in gradient echo MRI. Proc. Natl. Acad. Sci. U.S.A. 109, 18559–18564. doi:10.1073/pnas.1211075109

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wharton, S., Schäfer, A., and Bowtell, R. (2010). Susceptibility mapping in the human brain using threshold-based k-space division. Magn. Reson. Med. 63, 1292–1304. doi:10.1002/mrm.22334

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wu, B., Li, W., Avram, A. V., Gho, S.-M., and Liu, C. (2012). Fast and tissue-optimized mapping of magnetic susceptibility and T2* with multi-echo and multi-shot spirals. Neuroimage 59, 297–305. doi:10.1016/j.neuroimage.2011.07.019

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zheng, W., Nichol, H., Liu, S., Cheng, Y.-C. N., and Haacke, E. M. (2013). Measuring iron in the brain using quantitative susceptibility mapping and X-ray fluorescence imaging. Neuroimage 78, 68–74. doi:10.1016/j.neuroimage.2013.04.022

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text