Simultaneous Acquisition of Diffusion Tensor and Dynamic Diffusion MRI.

BACKGROUND: Dynamic diffusion magnetic resonance imaging (ddMRI) metrics can assess transient microstructural alterations in tissue diffusivity but requires additional scan time hindering its clinical application. PURPOSE: To determine whether a diffusion gradient table can simultaneously acquire data to estimate dynamic and diffusion tensor imaging (DTI) metrics. STUDY TYPE: Prospective. SUBJECTS: Seven healthy subjects, 39 epilepsy patients (15 female, 31 male, age ± 15). FIELD STRENGTH/SEQUENCE: Two-dimensional diffusion MRI (b = 1000 s/mm2 ) at a field strength of 3 T. Sessions in healthy subjects-standard ddMRI (30 directions), standard DTI (15 and 30 directions), and nested cubes scans (15 and 30 directions). Sessions in epilepsy patients-two 30 direction (standard ddMRI, 10 nested cubes) or two 15 direction scans (standard DTI, 5 nested cubes). ASSESSMENT: Fifteen direction DTI was repeated twice for within-session test-retest measurements in healthy subjects. Bland-Altman analysis computed bias and limits of agreement for DTI metrics using test-retest scans and standard 15 direction vs. 5 nested cubes scans. Intraclass correlation (ICC) analysis compared tensor metrics between 15 direction DTI scans (standard vs. 5 nested cubes) and the coefficients of variation (CoV) of trace and apparent diffusion coefficient (ADC) between 30 direction ddMRI scans (standard vs. 10 nested cubes). STATISTICAL TESTS: Bland-Altman and ICC analysis using a P-value of 0.05 for statistical significance. RESULTS: Correlations of mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were strong and significant in gray (ICC > 0.95) and white matter (ICC > 0.95) between standard vs. nested cubes DTI acquisitions. Correlation of white matter fractional anisotropy was also strong (ICC > 0.95) and significant. ICCs of the CoV of dynamic ADC measured using repeated cubes and nested cubes acquisitions were modest (ICC >0.60), but significant in gray matter. CONCLUSION: A nested cubes diffusion gradient table produces tensor-based and dynamic diffusion measurements in a single acquisition. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 1.

Highly Efficient and Accurate Deep Learning-Based Classification of MRI Contrast on a CPU and GPU.

Classifying MR images based on their contrast mechanism can be useful in image segmentation where additional information from different contrast mechanisms can improve intensity-based segmentation and help separate the class distributions. In addition, automated processing of image type can be beneficial in archive management, image retrieval, and staff training. Different clinics and scanners have their own image labeling scheme, resulting in ambiguity when sorting images. Manual sorting of thousands of images would be a laborious task and prone to error. In this work, we used the power of transfer learning to modify pretrained residual convolution neural networks to classify MRI images based on their contrast mechanisms. Training and validation were performed on a total of 5169 images belonging to 10 different classes and from different MRI vendors and field strengths. Time for training and validation was 36 min. Testing was performed on a different data set with 2474 images. Percentage of correctly classified images (accuracy) was 99.76%. (A deeper version of the residual network was trained for 103 min and showed slightly lower accuracy of 99.68%.) In consideration of model deployment in the real world, performance on a single CPU computer was compared with GPU implementation. Highly accurate classification, training, and testing can be achieved without use of a GPU in a relatively short training time, through proper choice of a convolutional neural network and hyperparameters, making it feasible to improve accuracy by repeated training with cumulative training sets. Techniques to improve accuracy further are discussed and demonstrated. Derived heatmaps indicate areas of image used in decision making and correspond well with expert human perception. The methods used can be easily extended to other classification tasks with minimal changes.

Nonhuman primates exposed to Zika virus in utero are not protected against reinfection at 1 year postpartum.

There is limited information about the impact of Zika virus (ZIKV) exposure in utero on the anti-ZIKV immune responses of offspring. We infected six rhesus macaque dams with ZIKV early or late in pregnancy and studied four of their offspring over the course of a year postpartum. Despite evidence of ZIKV exposure in utero, we observed no structural brain abnormalities in the offspring. We detected infant-derived ZIKV-specific immunoglobulin A antibody responses and T cell memory responses during the first year postpartum in the two offspring born to dams infected with ZIKV early in pregnancy. Critically, although the infants had acquired some immunological memory of ZIKV, it was not sufficient to protect them against reinfection with ZIKV at 1 year postpartum. The four offspring reexposed to ZIKV at 1 year postpartum all survived but exhibited acute viremia and viral tropism to lymphoid tissues; three of four reexposed offspring exhibited spinal cord pathology. These data suggest that macaque infants born to dams infected with ZIKV during pregnancy remain susceptible to postnatal infection and consequent neuropathology.

Determining the optimal postlabeling delay for arterial spin labeling using subject-specific estimates of blood velocity in the carotid artery.

BACKGROUND: Arterial spin labeling with 3D acquisition requires determining a single postlabeling delay (PLD) value. PLD affects the signal-to-noise ratio (SNR) per unit time as well as quantitative cerebral blood flow (CBF) values due to its bearing on the presence of a vascular signal. PURPOSE: To search for an optimal PLD for pseudocontinuous arterial spin labeling (pCASL) using patient-specific carotid artery blood velocity measurements. STUDY TYPE: Prospective. SUBJECTS: A control group of 11 volunteers with no known pathology. Corroboration was through a separate group of six volunteers and a noncontrol group of five sickle cell disease (SCD) patients. FIELD STRENGTH/SEQUENCE: Pseudocontinuous arterial spin labeling with 3D nonsegmented echo planar imaging acquisition at 3T. ASSESSMENT: A perfusion-based measure was determined over a range of PLDs for each of 11 volunteers. A third-order polynomial was used to find the optimal PLD where the defined measure was maximum. This was plotted against the corresponding carotid artery velocity to determine a relationship between the perfusion measure and velocity. Corroboration was done using a group of six volunteers and a noncontrol group of five patients with SCD. PLD was determined from the carotid artery velocity and derived relationship and compared with optimal PLD obtained from measured perfusion over a range of PLD values. Error between the perfusion measure at predicted and measured optimal PLD was determined. STATISTICAL TESTS: Chi-squared goodness of fit; Pearson correlation; Bland-Altman. RESULTS: Carotid artery velocity was 63.8 ± 6.6 cm/s (53.1 ≤ v ≤ 72.3 cm/s) while optimal PLD was 1374 ± 226.5 msec (1102 ≤ PLD ≤ 1787 msec) across the 11 volunteers. PLD as a function of carotid velocity was determined to be PLD = -31.94. v + 3410 msec (Pearson correlation -0.93). In six volunteers, mean error between the perfusion measure at predicted and measured optimal PLD was 1.35%. Pearson correlation between the perfusion measure at the predicted PLD and the measure obtained experimentally was r = 0.96 (P < 0.001). Bland-Altman revealed a slight bias of 1.3%. For the test case of five SCD patients, the mean error was 1.3%. DATA CONCLUSION: Carotid artery velocity was used to determine optimal PLD for pCASL with 3D acquisition. The derived relationship was used to predict optimal PLD and the associated perfusion measure, which was found to be accurate when compared with its measured counterpart. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:951-960.

Reduced distortion artifact whole brain CBF mapping using blip-reversed non-segmented 3D echo planar imaging with pseudo-continuous arterial spin labeling.

PURPOSE: To implement and evaluate interleaved blip-up, blip-down, non-segmented 3D echo planar imaging (EPI) with pseudo-continuous arterial spin labeling (pCASL) and post-processing for reduced susceptibility artifact cerebral blood flow (CBF) maps. MATERIALS AND METHODS: 3D EPI non-segmented acquisition with a pCASL labeling sequence was modified to include alternating k-space coverage along phase encoding direction (referred to as "blip-reversed") for alternating dynamic acquisitions of control and label pairs. Eight volunteers were imaged on a 3T scanner. Images were corrected for distortion using spatial shifting transformation of the underlying field map. CBF maps were calculated and compared with maps obtained without blip reversal using matching gray matter (GM) images from a high resolution 3D scan. Additional benefit of using the correction for alternating blip-up and blip-down acquisitions was assessed by comparing to corrected blip-up only and corrected blip-down only CBF maps. Matched Student t-test of overlapping voxels for the eight volunteers was done to ascertain statistical improvement in distortion. RESULTS: Mean CBF value in GM for the eight volunteers from distortion corrected CBF maps was 50.8±9.9ml/min/100 gm tissue. Corrected CBF maps had 6.3% and 4.1% more voxels in GM when compared with uncorrected blip up (BU) and blip down (BD) images, respectively. Student t-test showed significant reduction in distortion when compared with blip-up images and blip-down images (p<0.001). When compared with corrected BU and corrected BD only CBF maps, BU and BD corrected maps had 2.3% and 1% more voxels (p=0.006 and 0.04, respectively). CONCLUSION: Pseudo-continuous arterial spin labeling with non-segmented 3D EPI acquisition using alternating blip-reversed k-space traversal and distortion correction provided significantly better matching GM CBF maps. In addition, employing alternating blip-reversed acquisitions during pCASL acquisition resulted in statistically significant improvement over corrected blip-up and blip-down CBF maps.

Long T2 suppression in native lung 3-D imaging using k-space reordered inversion recovery dual-echo ultrashort echo time MRI.

OBJECTIVE: Long T2 species can interfere with visualization of short T2 tissue imaging. For example, visualization of lung parenchyma can be hindered by breathing artifacts primarily from fat in the chest wall. The purpose of this work was to design and evaluate a scheme for long T2 species suppression in lung parenchyma imaging using 3-D inversion recovery double-echo ultrashort echo time imaging with a k-space reordering scheme for artifact suppression. MATERIALS AND METHODS: A hyperbolic secant (HS) pulse was evaluated for different tissues (T1/T2). Bloch simulations were performed with the inversion pulse followed by segmented UTE acquisition. Point spread function (PSF) was simulated for a standard interleaved acquisition order and a modulo 2 forward-reverse acquisition order. Phantom and in vivo images (eight volunteers) were acquired with both acquisition orders. Contrast to noise ratio (CNR) was evaluated in in vivo images prior to and after introduction of the long T2 suppression scheme. RESULTS: The PSF as well as phantom and in vivo images demonstrated reduction in artifacts arising from k-space modulation after using the reordering scheme. CNR measured between lung and fat and lung and muscle increased from -114 and -148.5 to +12.5 and 2.8 after use of the IR-DUTE sequence. Paired t test between the CNRs obtained from UTE and IR-DUTE showed significant positive change (p < 0.001 for lung-fat CNR and p = 0.03 for lung-muscle CNR). CONCLUSION: Full 3-D lung parenchyma imaging with improved positive contrast between lung and other long T2 tissue types can be achieved robustly in a clinically feasible time using IR-DUTE with image subtraction when segmented radial acquisition with k-space reordering is employed.

Three-dimensional T1 and T2* mapping of human lung parenchyma using interleaved saturation recovery with dual echo ultrashort echo time imaging (ITSR-DUTE).

PURPOSE: To develop and assess a new technique for three-dimensional (3D) full lung T1 and T2* mapping using a single free breathing scan during a clinically feasible time. MATERIALS AND METHODS: A 3D stack of dual-echo ultrashort echo time (UTE) radial acquisition interleaved with and without a WET (water suppression enhanced through T1 effects) saturation pulse was used to map T1 and T2* simultaneously in a single scan. Correction for modulation due to multiple views per segment was derived. Bloch simulations were performed to study saturation pulse excitation profile on lung tissue. Optimization of the saturation delay time (for T1 mapping) and echo time (for T2* mapping) was performed. Monte Carlo simulation was done to predict accuracy and precision of the sequence with signal-to-noise ratio of in vivo images used in the simulation. A phantom study was carried out using the 3D interleaved saturation recovery with dual echo ultrashort echo time imaging (ITSR-DUTE) sequence and reference standard inversion recovery spin echo sequence (IR-SE) to compare accuracy of the sequence. Nine healthy volunteers were imaged and mean (SD) of T1 and T2* in lung parenchyma at 3T were estimated through manually assisted segmentation. 3D lung coverage with a resolution of 2.5 × 2.5 × 6 mm3 was performed and nominal scan time was recorded for the scans. Repeatability was assessed in three of the volunteers. Regional differences in T1/T2* values were also assessed. RESULTS: The phantom study showed accuracy of T1 values to be within 2.3% of values obtained from IR-SE. Mean T1 value in lung parenchyma was 1002 ± 82 ms while T2* was 0.85 ± 0.1 ms. Scan time was ∼10 min for volunteer scans. Mean coefficient of variation (CV) across slices was 0.057 and 0.09, respectively. Regional variation along the gravitational direction and between right and left lung were not significant (P = 0.25 and P = 0.06, respectively) for T1. T2* showed significant variation (P = 0.03) along the gravitational direction. Repeatability for three volunteers was within 0.7% for T1 and 1.9% for T2*. CONCLUSION: 3D T1 and T2* maps of the entire lung can be obtained in a single scan of ∼10 min with a resolution of 2.5 × 2.5 × 6 mm3 . LEVEL OF EVIDENCE: 2 J. Magn. Reson. Imaging 2017;45:1097-1104.

Dose correction for post-contrast T1 mapping of the heart: the MESA study.

Post-contrast myocardial T1 (T1(myo,c)) values have been shown to be sensitive to myocardial fibrosis. Recent studies have shown differences in results obtained from T1(myo,c) and extracellular volume fraction (ECV) with respect to percentage fibrosis. By exploring the relationship between blood plasma volume and T1(myo,c), the underlying basis for the divergence can be explained. Furthermore, dose administration based on body mass index (BMI), age and gender can mitigate the divergence in results. Inter-subject comparison of T1(myo,c) required adjustment for dose (in mmol/kg), time and glomerular filtration rate. Further adjustment for effective dose based on lean muscle mass reflected by blood/plasma volume was performed. A test case of 605 subjects from the MESA study who had undergone pre- and post-contrast T1 mapping was studied. T1(myo,c) values were compared between subjects with and without metabolic syndrome (MetS), between smoking and non-smoking subjects, and subjects with and without impaired glucose tolerance, before and after dose adjustment based on plasma volume. Comparison with ECV (which is dose independent), pre-contrast myocardial T1 and blood normalized myocardial T1 values was also performed to validate the correction. There were significant differences in T1(myo,c) (post plasma volume correction) and ECV between current and former smokers (p value 0.017 and 0.01, respectively) but not T1(myo,c) prior to correction (p = 0.12). Prior to dose adjustment for plasma volume, p value was <0.001 for T1(myo,c) between MetS and non-MetS groups and was 0.13 between subjects with and without glucose intolerance; after adjustment for PV, p value was 0.63 and 0.99. Corresponding ECV p values were 0.44 and 0.99, respectively. Overall, ECV results showed the best agreement with PV corrected T1(myo,c) (mean absolute difference in p values = 0.073) and pre-contrast myocardial T1 in comparison with other measures (T1(myo,c( prior to correction, blood/plasma T1 value normalized myocardium). Weight-based contrast dosing administered in mmol/kg results in a bias in T1 values which can lead to erroneous conclusions. After adjustment for lean muscle mass based on plasma volume, results from T1(myo,c) were in line with ECV derived results. Furthermore, the use of a modified equivalent dose adjusted for BMI, age, sex and hematocrit can be adopted for quantitative imaging.

Evaluation of optimized breath-hold and free-breathing 3D ultrashort echo time contrast agent-free MRI of the human lung.

PURPOSE: To evaluate an optimized stack of radials ultrashort echo time (UTE) 3D magnetic resonance imaging (MRI) sequence for breath-hold and free-breathing imaging of the human lung. MATERIALS AND METHODS: A 3D stack of ultrashort echo time radials trajectory was optimized for coronal and axial lower-resolution breath-hold and higher-resolution free-breathing scans using Bloch simulations. The sequence was evaluated in 10 volunteers, without the use of contrast agents. Signal-to-noise ratio (SNR) mean and 95% confidence interval (CI) were determined from separate signal and noise images in a semiautomated fashion. The four scanning schemes were evaluated for significant differences in image quality using Student's t-test. Ten clinical patients were scanned with the sequence and findings were compared with concomitant computed tomography (CT) in nine patients. Breath-hold 3D spokes images were compared with 3D stack of radials in five volunteers. A Mann-Whitney U-test was performed to test significance in both cases. RESULTS: Breath-hold imaging of the entire lung in volunteers was performed with SNR (mean = 42.5 [CI]: 35.5-49.5; mean = 34.3 [CI]: 28.6-40) in lung parenchyma for coronal and axial scans, respectively, which can be used as a quick scout scan. Longer respiratory triggered free-breathing scan enabled high-resolution UTE scanning with mean SNR of 14.2 ([CI]: 12.9-15.5) and 9.2 ([CI]: 8.2-10.2) for coronal and axial scans, respectively. Axial free-breathing scans showed significantly higher image quality (P = 0.008) than the three other scanning schemes. The mean score for comparison with CT was 1.67 (score 0: n = 0; 1: n = 3; 2: n = 6). There was no significant difference between CT and MRI (P = 0.25). 3D stack of radials images were significantly better than 3D spokes images (P < 0.001). CONCLUSION: The optimized 3D stack of radials trajectory was shown to provide high-quality MR images of the lung parenchyma without the use of MRI contrast agents. The sequence may offer the possibility of breath-hold imaging and provides greater flexibility in trading off slice thickness and parallel imaging for scan time.

Optimized three-dimensional sodium imaging of the human heart on a clinical 3T scanner.

PURPOSE: Optimization of sequence and sequence parameters to allow three-dimensional (3D) sodium imaging of the entire human heart in vivo in a clinically reasonable time. THEORY AND METHODS: A stack of spirals pulse sequence was optimized for cardiac imaging by considering factors such as spoiling, nutation angles, repetition time, echo time, T1/T2 relaxation, off-resonance, data acquisition window, motion, and segmented k-space acquisition. Simulations based on Bloch equations as well as the exact trajectory used for data acquisition provided the basis for choice of parameter combinations for sodium imaging. Sodium phantom scanning was used to validate the choice of parameters and for corroboration with simulations. In vivo cardiac imaging in six volunteers was also performed with an optimized sequence. RESULTS: Phantom studies showed good correlation with simulation results. Images obtained from human volunteers showed that the heart can be imaged with a nominal resolution of 5 × 5 × 10 mm(3) and with a signal-to-noise ratio >15 (in the septum) in about 6-10 minutes. Long axis views of the reformatted human heart show true 3D imaging capability. CONCLUSION: Optimization of the sequence and its parameters allowed in vivo 3D sodium imaging of the entire human heart in a clinically reasonable time.

Reduced scan time three-dimensional FLAIR using modulated inversion and repetition time.

BACKGROUND: The purpose of this study is to design and evaluate a new reduced scan time three-dimensional (3D) FLuid Attenuated Inversion Recovery (FLAIR) sequence. METHODS: The 3D FLAIR sequence was modified so that the repetition time was modulated in a predetermined smooth manner (3D mFLAIR). Inversion times were adjusted accordingly to maintain cerebrospinal fluid (CSF) suppression. Simulations were performed to determine SNR for gray matter (GM), white matter (WM), and CSF. Fourteen volunteers were imaged using the modified and product sequence. SNR measurements were performed in GM, WM, and CSF. Mean value and the 95% confidence interval ([CI]) were assessed. Scan time for the 3D FLAIR and 3D mFLAIR sequences was measured. RESULTS: There was no statistically significant difference in the SNR measured in GM (P value = 0.5; mean SNR = 42.8 [CI]: 38.2-45.5 versus 42.2 [CI]: 38.3-46.1 for 3D FLAIR and 3D mFLAIR, respectively) and WM (P value = 0.25; mean SNR = 32.1 [CI]: 30.3-33.8 versus 32.9 [CI]: 31.1-34.7). Scan time reduction greater than 30% was achieved for the given parameter set with the 3D mFLAIR sequence. CONCLUSION: Scan time for 3D FLAIR can be effectively reduced by modulating repetition and inversion time in a predetermined manner while maintaining the SNR and CNR of a constant TR sequence.

Modified Look-Locker T1 evaluation using Bloch simulations: human and phantom validation.

Modified Look-Locker imaging is frequently used for T(1) mapping of the myocardium. However, the specific effect of various MRI parameters (e.g., encoding scheme, modifications of flip angle, heart rate, T(2), and inversion times) on the accuracy of T(1) measurement has not been studied through Bloch simulations. In this work, modified Look-Locker imaging was characterized through a numerical solution for Bloch equations. MRI sequence parameters that may affect T(1) accuracy were systematically varied in the simulation. For validation, phantoms were constructed with various T(2) and T(1) times and compared with Bloch equation simulations. Human volunteers were also evaluated with various pulse sequences parameters to assess the validity of the numerical simulations. There was close agreement between simulated T(1) times and T(1) times measured in phantoms and volunteers. Lower T(2) times (i.e., <30 ms) resulted in errors greater than 5% for T(1) determination. Increasing maximum inversion time value improved T(1) accuracy particularly for precontrast myocardial T(1). Balanced steady-state free precession k space centric encoding improved accuracy for short T(1) times (post gadolinium), but linear encoding provided improved accuracy for precontrast T(1) values. Lower flip angles are preferred if the signal-to-noise ratio is sufficiently high. Bloch simulations for modified Look-Locker imaging provide an accurate method to comprehensively quantify the effect of pulse sequence parameters on T(1) accuracy. As an alternative to otherwise lengthy phantom studies or human studies, such simulations may be useful to optimize the modified Look-Locker imaging sequence and compare differences in T(1)-derived measurements from different scanners or institutions.

Fast T1 mapping determined using incomplete inversion recovery look-locker 3D balanced SSFP acquisition and a simple two-parameter model fit.

PURPOSE: To evaluate a fast T1 mapping technique using incomplete inversion recovery 3D balanced steady-state free precession acquisition along with a two-parameter model fit. MATERIALS AND METHODS: Using Bloch simulations, we explored the two-parameter model fit for data acquired using such an acquisition scheme. The parameter space over which the fit holds good was determined through simulations. A linear correction was derived for the R1* (1/T1*) values so determined. Two phantoms and six volunteers were scanned using the described technique. Comparison scans using full recovery as well as gold standard inversion recovery spin echo were also performed. RESULTS: The two-parameter fit works exceedingly well over a large parameter space. T1 values in the phantoms showed an error of 4.9% and 39% before correction and 0.9% and 1.6% after correction. For the six volunteers, error in T1 value was 5.3% for white matter (WM) and 2.4% for gray matter (GM) after correction, while it was 11.2% and 18.2% before correction. CONCLUSION: The work presented here allows for T1 map determination with higher resolution and shorter acquisition time than previously possible. The technique is especially well suited for GM/WM T1 mapping.

Tracking accuracy of T2- and diffusion-weighted magnetic resonance imaging for infusate distribution by convection-enhanced delivery.

OBJECT: Because convection-enhanced delivery relies on bulk flow of fluid in the interstitial spaces, MR imaging techniques that detect extracellular fluid and fluid movement may be useful for tracking convective drug distribution. To determine the tracking accuracy of T2-weighted and diffusion-weighted MR imaging sequences, the authors followed convective distribution of radiolabeled compounds using these imaging sequences in nonhuman primates. METHODS: Three nonhuman primates underwent thalamic convective infusions (5 infusions) with (14)C-sucrose (MW 342 D) or (14)C-dextran (MW 70,000 D) during serial MR imaging (T2- and diffusion-weighted imaging). Imaging, histological, and autoradiographic findings were analyzed. RESULTS: Real-time T2- and diffusion-weighted imaging clearly demonstrated the region of infusion, and serial images revealed progressive filling of the bilateral thalami during infusion. Imaging analysis for T2- and diffusion-weighted sequences revealed that the tissue volume of distribution (Vd) increased linearly with volume of infusion (Vi; R(2) = 0.94, R(2) = 0.91). Magnetic resonance imaging analysis demonstrated that the mean ± SD Vd/Vi ratios for T2-weighted (3.6 ± 0.5) and diffusion-weighted (3.3 ± 0.4) imaging were similar (p = 0.5). While (14)C-sucrose and (14)C-dextran were homogeneously distributed over the infused region, autoradiographic analysis revealed that T2-weighted and diffusion-weighted imaging significantly underestimated the Vd of both (14)C-sucrose (mean differences 51.3% and 52.3%, respectively; p = 0.02) and (14)C-dextran (mean differences 49.3% and 59.6%; respectively, p = 0.001). CONCLUSIONS: Real-time T2- and diffusion-weighted MR imaging significantly underestimate tissue Vd during convection-enhanced delivery over a wide range of molecular sizes. Application of these imaging modalities may lead to inaccurate estimation of convective drug distribution.

Whole-brain cerebral blood flow mapping using 3D echo planar imaging and pulsed arterial tagging.

PURPOSE: To quantitate cerebral blood flow (CBF) in the entire brain using the 3D echo planar imaging (EPI) PULSAR (pulsed star labeling) technique. MATERIALS AND METHODS: The PULSAR technique was modified to 1) incorporate a nonselective inversion pulse to suppress background signal; 2) to use 3D EPI acquisition; and 3) to modulate flip angle in such a manner as to minimize the blurring resulting from T1 modulation along the slice encoding direction. Computation of CBF was performed using the general kinetic model (GKM). In a series of healthy volunteers (n = 12), we first investigated the effects of introducing an inversion pulse on the measured value of CBF and on the temporal stability of the perfusion signal. Next we investigated the effect of flip angle modulation on the spatial blurring of the perfusion signal. Finally, we evaluated the repeatability of the CBF measurements, including the influence of the measurement of arterial blood magnetization (a calibration factor for the GKM). RESULTS: The sequence provides sufficient perfusion signal to achieve whole brain coverage in ≈ 5 minutes. Introduction of the inversion pulse for background suppression did not significantly affect computed CBF values, but did reduce the fluctuation in the perfusion signal. Flip angle modulation reduced blurring, resulting in higher estimates of gray matter (GM) CBF and lower estimates of white matter (WM) CBF. The repeatability study showed that measurement of arterial blood signal did not result in significantly higher error in the perfusion measurement. CONCLUSION: Improvements in acquisition and sequence preparation presented here allow for better quantification and localization of perfusion signal, allowing for accurate whole-brain CBF measurements in 5 minutes.

Modulated repetition time look-locker (MORTLL): a method for rapid high resolution three-dimensional T1 mapping.

PURPOSE: To demonstrate a modification of the Look-Locker (LL) technique that enables rapid high resolution T1 mapping over the physiologic range of intracranial T1 values, ranging from white matter to cerebrospinal fluid (CSF). This is achieved by use of a three-dimensional (3D) balanced steady-state free precession (b-SSFP) acquisition (for high signal-to-noise and resolution) along with variable repetition time to allow effective full recovery of longitudinal magnetization. MATERIALS AND METHODS: Two modifications to the Look-Locker technique were made to realize high resolution imaging in a clinically reasonable scan time. The 3D b-SSFP acquisition after an initial inversion pulse was followed by a variable repetition time. This technique makes it possible to image a volume of thin contiguous slices with high resolution and accuracy using a simple fitting procedure and is particularly useful for imaging long T1 species such as CSF. The total scan time is directly proportional to the number of slices to be acquired. The scan time was reduced by almost half when the repetition time was modified using a predesigned smooth function. Phantoms and volunteers were imaged at different resolutions on a 3 Tesla scanner. Results were compared with other accepted techniques. RESULTS: T1 values in the brain corresponded well with full repetition time imaging as well as inversion recovery spin echo imaging. T1 values for white matter, gray matter, and CSF were measured to be 755 +/- 10 ms, 1202 +/- 9 ms, and 4482 +/- 71 ms, respectively. Scan times were reduced by approximately half over full repetition time measurements. CONCLUSION: High resolution T1 maps can be obtained rapidly and with a relatively simple postprocessing method. The technique is particularly well suited for long T1 species. For example, changes in the composition of proteins in CSF are linked to various pathologies. The T1 values showed excellent agreement with values obtained from inversion recovery spin-echo imaging.

Design and optimization for variable rate selective excitation using an analytic RF scaling function.

At higher B(0) fields, specific absorption rate (SAR) deposition increases. Due to maximum SAR limitation, slice coverage decreases and/or scan time increases. Conventional selective RF pulses are played out in conjunction with a time independent field gradient. Variable rate selective excitation (VERSE) is a technique that modifies the original RF and gradient waveforms such that slice profile is unchanged. The drawback is that the slice profile for off-resonance spins is distorted. A new VERSE algorithm based on modeling the scaled waveforms as a Fermi function is introduced. It ensures that system related constraints of maximum gradient amplitude and slew rate are not exceeded. The algorithm can be used to preserve the original RF pulse duration while minimizing SAR and peak b1 or to minimize the RF pulse duration. The design is general and can be applied to any symmetrical or asymmetrical RF waveform. The algorithm is demonstrated by using it to (a) minimize the SAR of a linear phase RF pulse, (b) minimize SAR of a hyperbolic secant RF pulse, and (c) minimize the duration of a linear phase RF pulse. Images with a T1-FLAIR (T1 FLuid Attenuated Inversion Recovery) sequence using a conventional and VERSE adiabatic inversion RF pulse are presented. Comparison of images and scan parameters for different anatomies and coils shows increased scan coverage and decreased SAR with the VERSE inversion RF pulse, while image quality is preserved.