Emergence of resting state networks in the preterm human brain.

The functions of the resting state networks (RSNs) revealed by functional MRI remain unclear, but it has seemed possible that networks emerge in parallel with the development of related cognitive functions. We tested the alternative hypothesis: that the full repertoire of resting state dynamics emerges during the period of rapid neural growth before the normal time of birth at term (around 40 wk of gestation). We used a series of independent analytical techniques to map in detail the development of different networks in 70 infants born between 29 and 43 wk of postmenstrual age (PMA). We characterized and charted the development of RSNs from recognizable but often fragmentary elements at 30 wk of PMA to full facsimiles of adult patterns at term. Visual, auditory, somatosensory, motor, default mode, frontoparietal, and executive control networks developed at different rates; however, by term, complete networks were present, several of which were integrated with thalamic activity. These results place the emergence of RSNs largely during the period of rapid neural growth in the third trimester of gestation, suggesting that they are formed before the acquisition of cognitive competencies in later childhood.

A method for rapid in vivo measurement of blood T1.

We present a technique to measure the longitudinal relaxation time constant of venous blood (T(1b) ) in vivo in a few seconds. The MRI sequence consists of a thick-slab adiabatic inversion, followed by a series of slice-selective excitations and single-shot echo planar imaging readouts. The time intervals between excitations were chosen so that blood in macroscopic vessels is fully refreshed between excitations, making the blood signal follow an unperturbed inversion recovery curve. Static tissue, which experiences the inversion and all excitation pulses, quickly reaches a steady state at a low signal as a result of partial saturation. This allows blood-filled voxels to be discriminated from those containing static tissue, and to be fitted voxel-by-voxel to a simple inversion recovery model. The sequence was tested on a flow phantom with the proposed method, yielding T(1) values consistent to within 3% of those obtained using a conventional inversion recovery sequence with a spin-echo readout. The method was applied to seven adult volunteers and 18 neonates. The blood T(1) of the neonates (1799 ± 206 ms; range, 1393-2035 ms) was found to be more variable than that of adults (1717 ± 39 ms; range, 1662-1779 ms). A linear correlation between the inverse of T(1b) and the haematocrit was established in 12 neonates (R(2) = 0.90).

Subject-specific water-selective imaging using parallel transmission.

Spectral-spatial excitation pulses are an efficient means of achieving water- or fat-only imaging and can be used in conjunction with a variety of pulse sequences. However, the approach lacks reliability since its performance is dependent on the homogeneity of the static magnetic field. Sensitivity to static magnetic field variation can be reduced by designing pulses with wider frequency stop bands, but these require longer pulse durations. In the proposed method, spectral-spatial pulses are optimized on a subject-dependent basis to take into account measured subject-specific static magnetic field variation. Extra control of the radiofrequency (RF) field from multichannel transmission is used to achieve this without increasing the length of the pulses. The method characterizes RF pulses using relatively few parameters and has been applied to abdominal imaging at 3 T with an eight-channel system. In a comparison of standard and subject-specific pulses on five healthy volunteers, the latter improved fat suppression in all subjects, with a reduction in RF power of 13% +/- 6%. A forward model suggests that the mean flip angle in fat was reduced from 0.72 degrees +/- 0.55 degrees to 0.12 degrees +/- 0.04 degrees for a 20 degrees excitation; uniformity of water excitation also improved, with the standard deviation divided by mean reduced from 0.26 +/- 0.05 to 0.16 +/- 0.05.

PROPELLER for motion-robust imaging of in vivo mouse abdomen at 9.4 T.

In vivo high-field MRI in the abdomen of small animals is technically challenging because of the small voxel sizes, short T(2) and physiological motion. In standard Cartesian sampling, respiratory and gastrointestinal motion can lead to ghosting artefacts. Although respiratory triggering and navigator echoes can either avoid or compensate for motion, they can lead to variable TRs, require invasive intubation and ventilation, or extend TEs. A self-navigated fast spin echo (FSE)-based periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) acquisition was implemented at 9.4 T to enable high-resolution in vivo MRI of mouse abdomen without the use of additional navigators or triggering. T(2)-weighted FSE-PROPELLER data were compared with single-shot FSE and multi-shot FSE data with and without triggering. Single-shot methods, although rapid and robust to motion, demonstrated strong blurring. Multi-shot FSE data showed better resolution, but suffered from marked blurring in the phase-encoding direction and motion in between shots, leading to ghosting artefacts. When respiratory triggering was used, motion artefacts were largely avoided. However, TRs and acquisition times were lengthened by up to approximately 20%. The PROPELLER data showed a 25% and 61% improvement in signal-to-noise ratio and contrast-to-noise ratio, respectively, compared with multi-shot FSE data, together with a 35% reduction in artefact power. A qualitative comparison between acquisition methods using diffusion-weighted imaging was performed. The results were similar, with the exception that respiratory triggering was unable to exclude major motion artefacts as a result of the sensitisation to motion by the diffusion gradients. The PROPELLER data were of consistently higher quality. Considerations specific to the use of PROPELLER at high field are discussed, including the selection of practical blade widths and the effects on contrast, resolution and artefacts.

Manganese enhancement in non-CNS organs.

Manganese-enhanced magnetic resonance imaging (MEMRI) is a novel imaging technique capable of monitoring calcium influx, in vivo. Manganese (Mn2+) ions, similar to calcium ions (Ca2+), are taken up by activated cells where their paramagnetic properties afford signal enhancement in T(1)-weighted MRI methodologies. In this study we have assessed Mn2+ distribution in mice using magnetization-prepared rapid gradient echo (MP-RAGE) based MRI, by measuring changes in T(1)-effective relaxation times (T(1)-eff), effective R(1)-relaxation rates (R(1)-eff) and signal intensity (SI) profiles over time. The manganese concentration in the tissue was also determined using inductively coupled plasma atomic emission spectrometry (ICP-AES). Our results show a strong positive correlation between infused dose of MnCl2 and the tissue manganese concentration. Furthermore, we demonstrate a linear relationship between R(1)-eff and tissue manganese concentration and tissue-specific Mn2+ distribution in murine tissues following dose-dependent Mn2+ administration. This data provides an optimized MnCl2 dose regimen for an MP-RAGE based sequence protocol for specific target organs and presents a potential 3D MRI technique for in vivo imaging of Ca2+ entry during Ca2+-dependent processes in a wide range of tissues.

Combining RF encoding with parallel imaging: a simulation study.

OBJECT: The aim of this work was to investigate combining spatial encoding by radio frequency (RF) excitation with conventional parallel imaging (PI) methods to determine whether this could improve overall imaging performance. MATERIALS AND METHODS: A simulation framework was developed to predict imaging performance for regular, central and random under-sampled parallel imaging methods augmented by RF spatial signal modulation. Optimisation methods were used to find the RF modulation patterns that produce optimal image reconstruction using the condition number of the PI encoding matrix as a quality metric. The diverse patterns of raw data sampling produced were compared using a measure of data uniformity across k-space. RESULTS: Regular under-sampling of k-space provided the best reconstruction quality. When other under-sampling schemes were employed then RF modulation could be used to improve reconstruction, with the optimum achieved by redistributing the signal in k-space to return to regular sub-sampling. For all tested under-sampling patterns, no further improvements in image quality were attained. CONCLUSION: Using the simulation framework and metrics described the interaction of different spatial encoding approaches could be investigated. Regular sub-sampling provided optimal reconstruction, independent of whether the spatial encoding was achieved by gradients only or a combination of gradient and RF.

Optimal linear combinations of array elements for B1 mapping.

Accuracy of B(1) mapping for array coils can be improved by mapping the fields produced by driving linear combinations of the array elements, chosen to produce a more uniform distribution of B(1) amplitude. Quality of the resulting single element B(1) maps is influenced by the transformation used both via the uniformity of the resulting linear combination fields and by the degree to which these linear combinations differ from one another. In this work we investigate the effect of using different transformations on the quality of B(1) maps by simulating the B(1) mapping process for two different techniques, using real data from a 3T 8-channel body transmit system. Different transformations are generated using a single complex parameter. It is demonstrated that the optimal transformation within this framework is different for different imaging targets (pelvis and brain of healthy volunteers, and water and oil phantoms). For the same target (pelvis) the optimum condition, however, is similar for a number of subjects, suggesting that optimal configurations to be used for calibrating coils in specific anatomical contexts can be determined in advance. Potential gains may be translated into significant reductions in scan time for equivalent signal-to-noise ratio coil maps.

An efficient automated z-shim based method to correct through-slice signal loss in EPI at 3T.

OBJECT: To develop an efficient, automated method to correct through-slice signal loss in gradient-echo EPI at 3T. MATERIAL AND METHODS: The optimal choice of two z-shim values for signal recovery was determined from simulations and experiments. The specific required z-shim values are determined using a rapid calibration method that combines information about the slice profile with a sparse set of measurements. The proposed correction method was implemented for a language fMRI study which suffers from signal loss near the auditory canals, and tested on 12 volunteers. RESULTS: Using a square root sum of squares combination of two z-shim values full signal restoration (to within 2% of the correct value) was achieved in 96% of all correctable brain pixels for 3 mm slices, and partial correction in pixels outside this range. In all subjects, language processing activation was recovered in the inferior and lateral areas of the left temporal lobe which was not detectable with conventional fMRI. CONCLUSION: The careful choice of two z-shim values by the proposed method achieves through-slice signal loss correction for the majority of pixels in the brain for 3 mm slices at 3T.

Parallel magnetic resonance imaging.

Parallel imaging has been the single biggest innovation in magnetic resonance imaging in the last decade. The use of multiple receiver coils to augment the time consuming Fourier encoding has reduced acquisition times significantly. This increase in speed comes at a time when other approaches to acquisition time reduction were reaching engineering and human limits. A brief summary of spatial encoding in MRI is followed by an introduction to the problem parallel imaging is designed to solve. There are a large number of parallel reconstruction algorithms; this article reviews a cross-section, SENSE, SMASH, g-SMASH and GRAPPA, selected to demonstrate the different approaches. Theoretical (the g-factor) and practical (coil design) limits to acquisition speed are reviewed. The practical implementation of parallel imaging is also discussed, in particular coil calibration. How to recognize potential failure modes and their associated artefacts are shown. Well-established applications including angiography, cardiac imaging and applications using echo planar imaging are reviewed and we discuss what makes a good application for parallel imaging. Finally, active research areas where parallel imaging is being used to improve data quality by repairing artefacted images are also reviewed.

Imaging molecular and cellular events in transplantation.

Imaging methods such as nuclear medicine (including positron emission tomography), magnetic resonance imaging, ultrasound, and optical imaging can be used to provide information about the expression of genes, and the location of molecules and cells in intact animals or patients. In the setting of transplantation, this will allow monitoring of inflammatory responses, as well as the state of the graft. In this review, the advantages and disadvantages of different approaches to imaging will be discussed, as well as their potential application to transplantation.

Current and future applications of magnetic resonance imaging and spectroscopy of the brain in hepatic encephalopathy.

Hepatic encephalopathy (HE) is a common neuro-psychiatric abnormality, which complicates the course of patients with liver disease and results from hepatocellular failure and/or portosystemic shunting. The manifestations of HE are widely variable and involve a spectrum from mild subclinical disturbance to deep coma. Research interest has focused on the role of circulating gut-derived toxins, particularly ammonia, the development of brain swelling and changes in cerebral neurotransmitter systems that lead to global CNS depression and disordered function. Until recently the direct investigation of cerebral function has been difficult in man. However, new magnetic resonance imaging (MRI) techniques provide a non-invasive means of assessment of changes in brain volume (coregistered MRI) and impaired brain function (fMRI), while proton magnetic resonance spectroscopy (1H MRS) detects changes in brain biochemistry, including direct measurement of cerebral osmolytes, such as myoinositol, glutamate and glutamine which govern processes intrinsic to cellular homeostasis, including the accumulation of intracellular water. The concentrations of these intracellular osmolytes alter with hyperammonaemia. MRS-detected metabolite abnormalities correlate with the severity of neuropsychiatric impairment and since MR spectra return towards normal after treatment, the technique may be of use in objective patient monitoring and in assessing the effectiveness of various treatment regimens.

Artifact reduction using parallel imaging methods.

Multiple receiver coils produce images with different but complementary views of a patient. This can be used to shorten scans times but there often remain image artifacts caused by patient motion or physiological processes such as flowing blood. This paper reviews how the extra information from the multiple coils can be used to reduce image artifacts. In one method, affected portions of data can be identified and discarded but enough information is still available to reconstruct an improved image. In other methods, the motion itself is determined and the corrupted data is then corrected, leading to an image with reduced artifacts. Results are presented from images corrupted by motion or by flowing blood.

Thalamo-cortical connectivity in children born preterm mapped using probabilistic magnetic resonance tractography.

Our aim was to investigate the feasibility of studying white matter tracts and connections between the thalamus and the cortex in 2-year-old infants who were born preterm by probabilistic magnetic resonance (MR) tractography. Using this approach, we were able to visualize and quantify connectivity distributions in a number of white matter tracts, including the corticospinal tracts, optic radiations, fibers of the genu and splenium of the corpus callosum, superior longitudinal fasciculus and inferior fronto-occipital fasciculus, and to map the distribution within thalamus of fibers connecting to specific cortical regions. In eleven infants with no MR evidence of focal cerebral lesions and appropriate neurodevelopment as shown by general quotient (GQ) scores above 100, we mapped cortical connections to the thalamus that appeared similar to those reported in adults. However, in a proof-of-principle experiment, we examined one further child with marked white matter abnormalities and found that the volume and pattern of thalamo-cortical connections were severely disrupted. This technique promises to be a useful tool for assessing connectivity in the developing brain and in infants with lesions.

Beyond the g-factor limit in sensitivity encoding using joint histogram entropy.

The maximum practical speed-up that can be achieved using parallel imaging methods is widely accepted to be limited by g-factor noise. An approximate expression for the g-factor noise as a function of the principal eigenvector of the inverse sensitivity matrix is derived. This formulation allows g-factor enhanced noise to be reduced by a constrained optimization procedure with joint image histogram entropy between a reference image and a SENSE image as an image quality metric. The reference image does not need to have identical resolution or contrast. The reference image may also be used for coil calibration. The limits of the method are explored using simulated and real array coil data with high g-factor using a variety of contrast and resolution combinations. The method preserves image structure, contrast, and lesions even when these were not observable in the reference data. In all cases g-factor was dramatically reduced.

x-f Choice: reconstruction of undersampled dynamic MRI by data-driven alias rejection applied to contrast-enhanced angiography.

A technique for reconstructing dynamic undersampled MRI data, termed "x-f choice," was developed and applied to dynamic contrast-enhanced MR angiography (DCE-MRA). Regular undersampling in k-t space (a hybrid of k-space and time) creates aliasing in the conjugate x-f space that must be resolved. When regions in the object containing fast dynamic change are sparse, as in DCE-MRA, signal overlap caused by aliasing is often much less than the undersample factor would imply. x-f Choice reconstruction identifies overlapping signals using a model of the full non-aliased x-f space that is automatically generated from the undersampled data, and applies parallel imaging (PI) to separate them. No extra reference scans are required to generate either the model or the coil sensitivity maps. At each location in the reconstructed images, g-factor noise amplification is compared with predicted reconstruction errors to obtain an optimized solution. Acceleration factors greater than the number of receiver coils are possible, but are limited by the sparseness of the dynamic content and the signal-to-noise ratio (SNR) (in DCE-MRA the latter is dominant). Temporal fidelity was validated for up to a factor 10 speed-up using retrospectively undersampled data from a six-coil array. The method was tested on volunteers using fivefold prospective undersampling.

Nonlinear phase correction of navigated multi-coil diffusion images.

Cardiac pulsatility causes a nonrigid motion of the brain. In multi-shot diffusion imaging this leads to spatially varying phase changes that must be corrected. A conjugate gradient based reconstruction is presented that includes phase changes measured using two-dimensional navigator echoes, coil sensitivity information, navigator-determined weightings, and data from multiple coils and averages.A multi-shot echo planar sequence was used to image brain regions where pulsatile motion is not uniform. Reduced susceptibility artifacts were observed compared to a clinical single-shot sequence. In a higher slice, fiber directions derived from single-shot data show distortions from anatomical scans by as much as 7 mm compared to less than 2 mm for our multi-shot reconstructions. The reduced distortions imply that phase encoding can be applied in the shorter left-right direction, enabling time savings through the use of a rectangular field of view. Higher resolution diffusion imaging in the spine permits visualization of a nerve root.

Detection of vascular expression of E-selectin in vivo with MR imaging.

PURPOSE: To develop a contrast agent for targeting E-selectin expressed on activated vascular endothelium and to evaluate detection of the agent with magnetic resonance (MR) imaging in an in vivo mouse model of inflammation. MATERIALS AND METHODS: All animal experiments were approved according to animal welfare and local ethics committee regulations. An anti-murine E-selectin F(ab')2 monoclonal antibody, MES-1, was conjugated with ultrasmall superparamagnetic iron oxide (USPIO) nanoparticles. Flow cytometry, Perl Prussian blue staining for iron, and MR imaging were performed by using Chinese hamster ovary (CHO) cells expressing mouse E-selectin to detect binding of the conjugate in vitro, and a mouse model of contact hypersensitivity to oxazolone in the ear was used to investigate the in vivo characteristics of the MES-1-USPIO. Serial imaging was performed by using a 9.4-T MR imaging system with a custom receive-only coil. Tissue slices were stained to define distribution of E-selectin expression and localization of the MES-1-USPIO conjugate. RESULTS: MES-1-USPIO was shown to bind to CHO cells expressing mouse E-selectin in vitro. After injection of MES-1-USPIO in vivo, distinct changes in R2 relaxation rate (1/T2) characteristics were detected in inflamed ears when they were compared with control ears. Histologic analysis confirmed the vascular endothelial distribution of MES-1-USPIO. CONCLUSION: E-selectin expression in vivo can be selectively and directly imaged noninvasively with MR. This has the potential to be useful in the study of inflammatory disease.

Axial and radial diffusivity in preterm infants who have diffuse white matter changes on magnetic resonance imaging at term-equivalent age.

OBJECTIVE: Diffuse excessive high signal intensity (DEHSI) is observed in the majority of preterm infants at term-equivalent age on conventional MRI, and diffusion-weighted imaging has shown that apparent diffusion coefficient values are elevated in the white matter (WM) in DEHSI. Our aim was to obtain diffusion tensor imaging on preterm infants at term-equivalent age and term control infants to test the hypothesis that radial diffusivity was significantly different in the WM in preterm infants with DEHSI compared with both preterm infants with normal-appearing WM on conventional MRI and term control infants. METHODS: Diffusion tensor imaging was obtained on 38 preterm infants at term-equivalent age and 8 term control infants. Values for axial (lambda1) and radial [(lambda2 + lambda3)/2] diffusivity were calculated in regions of interest positioned in the central WM at the level of the centrum semiovale, frontal WM, posterior periventricular WM, occipital WM, anterior and posterior portions of the posterior limb of the internal capsule, and the genu and splenium of the corpus callosum. RESULTS: Radial diffusivity was elevated significantly in the posterior portion of the posterior limb of the internal capsule and the splenium of the corpus callosum, and both axial and radial diffusivity were elevated significantly in the WM at the level of the centrum semiovale, the frontal WM, the periventricular WM, and the occipital WM in preterm infants with DEHSI compared with preterm infants with normal-appearing WM and term control infants. There was no significant difference between term control infants and preterm infants with normal-appearing WM in any region studied. CONCLUSIONS: These findings suggest that DEHSI represents an oligodendrocyte and/or axonal abnormality that is widespread throughout the cerebral WM.

Padé methods for reconstruction and feature extraction in magnetic resonance imaging.

Methods utilizing Padé approximants are investigated for implementation with magnetic resonance imaging data and are presented both for direct image reconstruction and for feature extraction. Padé approximants are a numerical tool that can be used to accelerate the convergence of a slowly converging sequence by estimating the fully converged sequence values from early data points. Padé approximants can be calculated directly from k-space data by solving a set of linear matrix equations to produce signal values for any desired location in the image domain. This gives an estimate of the fully converged signal intensity at each pixel location in the image, raising the possibility of reconstructing a better estimate of the object from a reduced data set. These methods have been tested on phantom and human data both for image reconstruction and for feature extraction. In image reconstruction, considerable convergence acceleration can be achieved, with steep intensity boundaries reproduced in keeping with higher resolution reconstructions and oscillatory truncation artifact characteristic of Fourier reconstruction removed. The convergence acceleration is variable and there is the possibility of fine detail suppression when insufficient data are included. The use of Padé methods as a tool for feature extraction has shown good agreement with extraction from high-resolution reference data. In this approach the edge information comes intrinsically from Padé reconstruction.

Generalized SMASH imaging.

A generalized parallel imaging method has been developed that uses coil profiles to generate missing k-space lines. The proposed method is an extension of SMASH, which uses linear combinations of coil sensitivity profiles to synthesize spatial harmonics. In the generalized SMASH approach described here, coil sensitivity profiles are represented directly in the Fourier domain to provide a general description of the spatial properties of the coils. This removes restrictions imposed by conventional SMASH, so that the choice and position of the receiver coils can be made on the basis of sensitivity to the volume of interest rather than suitability for constructing spatial harmonics. Generalized SMASH also intrinsically allows the freedom to accommodate acquisitions with uniform or nonuniform k-space sampling. The proposed method places SMASH on an equal footing with other parallel imaging techniques (SENSE and SPACE-RIP), while combining strengths from each. The method was tested on phantom and human data and provides a robust method of data recovery.