Automatic model-based evaluation of magnetic resonance-guided radio frequency ablation lesions with histological correlation

PURPOSE: To develop a model-based method for automatic evaluation of radio frequency (RF) ablation treatment using magnetic resonance (MR) images. MATERIALS AND METHODS: RF current lesions were generated in a rabbit thigh model using MR imaging (MRI) guidance. We created a 12-parameter, three-dimensional, globally deformable model with quadric surfaces that delineates lesion boundaries and is automatically fitted to MR grayscale data. We applied this method to in vivo T2- and contrast-enhanced (CE) T1-weighted MR images acquired immediately post-ablation and four days later. We then compared results to manually segmented MR and three-dimensional registered corresponding histological boundaries of cellular damage. RESULTS: Resulting lesions featured a two-boundary appearance with an inner region and an outer hyperintense margin on MR images. For automated vs. manual MR boundaries, the mean errors over all specimens were 0.19 +/- 0.51 mm and 0.27 +/- 0.52 mm for the inner surface, and -0.29 +/- 0.40 mm and -0.12 +/- 0.17 mm for the outer surface, for T2- and CE T1-weighted images, respectively. For automated vs. histological boundaries, mean errors over all specimens were 0.07 +/- 0.64 mm and 0.33 +/- 0.71 mm for the inner surface, and -0.27 +/- 0.69 mm and 0.02 +/- 0.43 mm for the outer surface, for T2- and CE T1-weighted images, respectively. All boundary errors compared favorably to MR voxel dimensions, which were 0.7 mm in-plane and 3.0 mm thick. CONCLUSION: The method is accurate both in describing MR-apparent boundaries and in predicting histological response and has applications in lesion visualization, volume estimation, and treatment evaluation.     

Three-dimensional method for comparing in vivo interventional MR images of thermally ablated tissue with tissue response

PURPOSE: To investigate the ability of magnetic resonance (MR) to monitor radio-frequency (RF) ablation treatments by comparing MR images of thermal lesions to histologically assayed cellular damage. We developed a new methodology using three-dimensional registration for making spatial correlations. MATERIALS AND METHODS: A low-field, open MRI system was used to guide an ablation probe into rabbit thigh muscle and acquire MR volumes after ablation. After fixation, we sliced and photographed the tissue at 3-mm intervals, using a specially designed apparatus, to obtain a volume of tissue images. Histologic samples were digitized using a video microscopy system. For our three-dimensional registration method, we used the tissue images as the reference, and registered histology and MR images to them using two different computer alignment steps. First, the MR volume was aligned to the volume of tissue images by registering needle fiducials placed near the tissue of interest. Second, we registered the histology images with the tissue images using a two-dimensional warping technique that aligned internal features and the outside boundary of histology and tissue images. RESULTS: The MR and histology images were very well aligned, and registration accuracy, determined from displacement of needle fiducials, was 1.32 +/- 0.39 mm (mean +/- SD), which compared favorably to the MR voxel dimensions (0.70 mm in-plane and 3.0 mm thick). A preliminary comparison of MR and tissue response showed that the region inside the elliptical hyperintense rim in MR closely corresponds to the region of necrosis as established by histology, with a mean absolute distance between MR and histology boundaries of 1.17 mm, slightly smaller than the mean registration error. The MR region slightly overestimated the region of necrosis, with a mean signed distance between boundaries of 0.85 mm. CONCLUSION: Our results suggest that our methodology can be used to achieve three-dimensional registration of histology and in vivo MR images. In MR lesion images, the inner border of the hyperintense region corresponds to the border of irreversible cell damage. This is good evidence that during RF ablation treatments, iMRI lesion images can be used for real-time feedback.     

Fiducial markers for MR histological correlation in ex vivo or short-term in vivo animal experiments: a screening study

PURPOSE: To test injectable fiducial markers for magnetic resonance (MR) histological correlation in ex vivo or in vivo animal experiments. MATERIALS AND METHODS: A total of 35 potential markers were tested ex vivo in pork muscle. The end-points were: 1) visibility, size, and shape on MR images and at macroscopic examination; 2) 24-hour stability; and 3) microscopic appearance. Selected markers were injected in vivo (rabbit's muscle and breast tumor tissue) to test their three-hour in vivo stability and their potential toxicity. Finally, different dilutions of the two best markers were assessed again through the same screening tests to determine whether their size on MR images could be customized by dilution. RESULTS: Two fluid acrylic paints containing inorganic pigments were found to be potentially interesting markers. On MR images, they created well-defined susceptibility artifacts. The markers made with iridescent bronze paint (iron oxide coated mica particles) were readily visible on microscopy and their size on MR images could be customized by dilution. The iridescent stainless steel paint (iron, chromium, nickel) created ex vivo the smallest markers in tissue but needed colloidal iron staining to be visible on microscopy and could not be easily diluted. CONCLUSION: Fluid acrylic paints are potentially interesting markers for MR histological correlation. Further studies are needed to assess their long-term properties.     

Radiofrequency thermal ablation: correlation of hyperacute MR lesion images with tissue response

PURPOSE: To investigate the hypothesis that the outer boundary of the hyperintense region observed in hyperacute (several minutes post-ablation) T2 and gadolinium contrast-enhanced (CE) T1-weighted magnetic resonance (MR) lesion images is an accurate predictor of eventual cell death from radiofrequency (RF) thermal ablation. MATERIALS AND METHODS: A low-field, open MR imaging system was used to guide an ablation electrode into a thigh muscle of five rabbits and acquire in vivo T2 and CE T1-weighted MR volumes. Ablation occurred by applying RF current for two minutes with the electrode's temperature maintained at 90 degrees +/- 2 degrees C. After fixation, we sliced and photographed the tissue at 3 mm intervals, using a specially designed apparatus, to obtain a volume of tissue images. Digital images of hematoxylin and eosin (H&E) and Masson trichrome-stained histologic samples were obtained, and distinct regions of tissue damage were labeled using a video microscopy system. After the MR and histology images were aligned using a three-dimensional registration method, we compared tissue damage boundaries identified in histology with boundaries marked in MR images. RESULTS: The lesions have distinct zones of tissue damage histologically: a central zone of necrotic cells surrounded by an outer zone with cells that appeared non-viable and associated with marked interstitial edema. In 14 histology images from five lesions, the inner and outer boundaries of the outer zone were compared with the boundaries of a hyperintense rim that surrounds a central hypointense region in the T2 and CE T1-weighted MR images. For T2 and CE T1-weighted MR images, respectively, the mean absolute distance was 1.04 +/- 0.30 mm (mean +/- SD) and 1.00 +/- 0.34 mm for the inner boundaries, and 0.96 +/- 0.34 mm and 0.94 +/- 0.44 mm for the outer boundaries. The mean absolute distances for T2 and CE T1-weighted MR images were not sufficiently different to achieve statistical significance (P = 0.745, 0.818, for the inner and outer boundary, respectively). CONCLUSION: In hyperacute T2 and CE T1-weighted MR lesion images, observations strongly suggest that the outer boundary of the hyperintense rim corresponds to the region of eventual cell necrosis within a distance comparable to our ability to measure. This is good evidence that during RF ablation procedures, MR lesion images can be used to accurately localize the zone of irreversible tissue damage at the lesion margin.     

Methods for registration of magnetic resonance images of ex vivo prostate specimens with histology

Purpose:

To evaluate two methods of scanning and tissue processing to achieve accurate magnetic resonance (MR)‐histologic correlation in human prostate specimens.

Materials and Methods:

Two prostates had acrylic paint markers injected to define the plane of imaging and serve as internal fiducials. Each was placed on a polycarbonate plane‐finder device (PFD), which was adjusted to align the imaging and cutting planes. Three prostates were aligned by use of a plane finder key (PFK), a polycarbonate plate that locks the specimen in a cylindrical carrier. Markers were injected for registration analysis. Prostates were imaged, then sectioned. Imaging software was used to create registration maps of the MR and histology images. Measurements between control points were made and compared.

Results:

Accurate correlation was achieved between MR and histologic images. The mean displacement (MD) between the corresponding registration points using the PFD technique ranged from 1.11–1.38 mm for each section. The MD for all sections was 1.24 mm. The MD using the PFK technique ranged from 0.79–1.01 mm for each section, and the MD across all sections for the PFK was 0.92 mm.

Conclusion:

We describe two methods that can achieve accurate, reproducible correlation between MR imaging and histologic sections in human prostatectomy specimens. J. Magn. Reson. Imaging 2012;36:206–212. © 2012 Wiley Periodicals, Inc.     

Three-dimensional magnetic resonance imaging technique for myocardial-delayed hyperenhancement: a comparison with the two-dimensional technique

PURPOSE: To compare two-dimensional and three-dimensional techniques in the detection of myocardial infarction (MI) and in the grading transmural extent (TE). MATERIALS AND METHODS: Twelve patients with clinically proven MI were examined using two-dimensional and three-dimensional techniques with cardiac-gated, breath-hold, T1-weighted gradient echo sequence with an inversion recovery pulse following gadopentetate dimeglumine (Gd-DTPA) at 0.2 mmol/kg. Contrast-to-noise, signal-to-noise, and signal intensity ratios (CNR, SNR, and SIR, respectively) were derived and compared for each technique. RESULTS: From two-dimensional to three-dimensional, statistical significant difference was found in the mean CNR (11.65 vs. 56.59; P = 0.002), SNR (18.03 vs. 76.90; P < 0.001), and SIR (3.6 vs. 6.36; P = 0.05). Intraobserver agreement (kappa) between two-dimensional and three-dimensional were R1 = 74% and R2 = 90%. Interobserver agreements between the readers were two-dimensional = 77% and three-dimensional = 79%. CONCLUSION: Mean CNR, SNR, and SIR are significantly increased in the three-dimensional technique compared to the conventional two-dimensional technique.     

MR histological correlation: a method for cutting specimens along the imaging plane in animal or ex vivo experiments

PURPOSE: To assess a method aimed at cutting histological specimens along the magnetic resonance (MR) imaging plane. MATERIAL AND METHODS: The method is performed in two steps: the imaging plane (defined by three acrylic paint markers) is made horizontal under MR guidance by using a mobile platform that can be rotated in three directions (PlaneFinder device [PFD]); then, the specimen is embedded in wax and cut horizontally. Three-dimensional images parallel to the markers' plane were obtained on 31 pork muscles containing a central hole with a pyramidal shape, with a technique of reference (RT images) and with PFD (PF images), before and after fixation. The last 17 fixed specimens were cut in the markers' plane (tissue section [TS] images). The central hole area (CHA) in the markers' plane was used to compare RT, PF, and TS images. Using a workstation, PF images were rotated and translated to estimate the shift along each direction that could explain the entire CHA difference between RT, PF, and TS images (maximum error, worst-case scenario). RESULTS: Excellent correlation was found between RT and PF images (r = 0.989, slope = 1.0175), PF and TS images (r = 0.991, slope = 1.0058), and RT images on fresh specimens and TS images (r = 0.979, slope = 1.0732). For each step, the maximum angle error was < or = 3 degrees in 88-95% of the specimens. CONCLUSION: Our methodology can be used to cut specimens along the imaging plane with high accuracy.     

Interleaved echo-planar imaging for fast multiplanar magnetic resonance temperature imaging of ultrasound thermal ablation therapy

PURPOSE: To develop a multiplanar magnetic resonance temperature imaging (MRTI) technique based on interleaved gradient-echo echo-planar imaging (EPI), verify in phantom, develop software tools to process and display data on a clinical scanner in near real-time, and demonstrate feasibility to monitor ultrasound thermal ablation therapy in vivo. MATERIALS AND METHODS: Temperature estimation used complex phase-difference subtraction of the EPI MRTI data to indirectly measure the temperature-dependent water proton-resonance-frequency shift. Software tools were developed to run on a clinical 1.5-T MR scanner that processed and displayed relevant temperature and thermal dosimetry data during the course of thermal ablation treatments in canine brain and prostate in vivo. RESULTS: EPI MRTI provided multi-planar acquisitions and increased temperature sensitivity and lipid suppression. Relative to a single-plane fast gradient-echo MRTI sequence at comparable spatial and temporal resolutions in phantom, EPI MRTI demonstrated a three-fold increase in sensitivity and slice coverage per TR. In vivo monitoring of ultrasound thermal ablation therapy in canine brain and prostate demonstrated the usefulness of the temperature and thermal dose information. CONCLUSION: Multi-planar MRTI allowed progression of thermal damage to be monitored and treatment parameters adjusted in near real-time (less than five second delay). EPI MRTI is an effective multi-planar monitoring method during ultrasound thermal ablation procedures.     

Reactive sclerosis: hyperintense appearance on T2-weighted magnetic resonance imaging

Osteosclerosis is defined as increased density of bone on X-ray imaging studies. It is known that osteosclerosis appears hypointense on both T1- and T2-weighted magnetic resonance imaging sequences. In this review, we present our experience in various sclerotic skeletal pathologies that appear hyperintense on T2-weighted magnetic resonance imaging. We emphasize the possible pathophysiological mechanisms that may cause this appearance, such as bone marrow edema and/or composition of newly formed bone.     

Sub-acute changes in lesion conspicuity and geometry following MR-guided radiofrequency ablation

PURPOSE: To evaluate MR signal and lesion zone volume evolution through the sub-acute phase following image-guided radiofrequency (RF) thermal ablation. MATERIALS AND METHODS: For many tissues, including muscle and liver, thermal lesions that result from RF heating have a characteristic two-boundary appearance featuring an inner core (zone I) surrounded by a hyper-intense margin (zone II) and normal tissue (zone III), found in both T(2) and contrast enhanced (CE) T(1)-weighted MR images, both immediately post-ablation and four days later. First, we compared corresponding points between manually segmented zone boundaries apparent in T(2)- and CE T(1)-weighted images. Second, we examined the contrast-to-noise ratio (CNR) between all zone combinations. Third, we quantified the volume of zone I, zone II, and the entire lesion using a three-dimensional lesion geometry model fitted to segmented images. RESULTS: On a slice-by-slice basis, no statistically significant differences were found between zone boundaries apparent in T(2) and CE T(1)-weighted images. The contrast to noise ratio (CNR) of zone I vs. zone II, zone I vs. background muscle, and zone II vs. background muscle was always equal or greater for T(2)-weighted images than for CE T(1)-weighted images. In addition, by day four, zone II significantly increased in intensity compared to background muscle. The median Zone I volume increase was 44.2% (42.6%) using T(2) weighted images and 55.5% (68.7% interquartile range) using CE T(1)- weighted images. This expansion likely corresponds to an enlargement of the ablated, coagulative necrosis, region. The median Zone II volume increase was 15.0% (42.6%) using T(2)- weighted images 1.5% (38.8%) using CE T(1)-weighted images. CONCLUSIONS: 1) There are no significant differences between the apparent zone boundaries in T(2)- and CE T(1)-weighted images; 2) CNR is equal or greater for T(2)-weighted images as compared to CE T(1)-images; and 3) both the inner and outer lesion zone volumes typically increase several days post-ablation.     

Characterization of regional pulmonary mechanics from serial magnetic resonance imaging data

RATIONALE AND OBJECTIVES: The aim of this study was to investigate a method for quantifying lung motion from the registration of successive images in serial magnetic resonance imaging acquisitions during normal respiration. MATERIALS AND METHODS: Estimates of pulmonary motion were obtained by summing the normalized cross-correlation over serially acquired lung images to identify corresponding locations between the images. The estimated motions were modeled as deformations of an elastic body and thus reflect to a first order approximation the true physical behavior of lung parenchyma. The Lagrangian strain, derived from the calculated motion fields, were used to quantify the tissue deformation induced in the lung over the serial acquisition. RESULTS: The method was validated on a magnetic resonance imaging study, for which breath-hold images were acquired of a healthy volunteer at different phases of the respiratory cycle. Regional parenchymal strain was observed to be oriented toward the pulmonary hilum, with strain magnitude maximal at the midcycle of the expiratory phase. CONCLUSION: In vivo magnetic resonance imaging quantification of lung motion holds the potential of providing a new diagnostic dimension in the assessment of pulmonary function, augmenting the information provided by studies of ventilation and perfusion.     

Central nervous system aspergillosis: magnetic resonance imaging, diffusion-weighted imaging, and magnetic resonance spectroscopy features

Aspergillus infection is invasive in nature in the immunosuppressed population and disseminates throughout the body, with the brain being a common site. Conventional magnetic resonance imaging (MRI) combined with diffusion-weighted imaging (DWI) and magnetic resonance spectroscopy (MRS) play a life-saving role in the early diagnosis and treatment monitoring of this potentially fatal infection. We present MRI, DWI, and MRS findings of a case of central nervous system aspergillosis with treatment follow-up.     

Contrast-enhanced magnetic resonance cholangiography using mangafodipir compared with standard T2W MRC sequences: a pictorial essay

Mangafodipir, a manganese-containing hepatobiliary contrast agent, is excreted in bile. We review the principles and practice of a contrast-enhanced MRC technique using mangafodipir and compare it with standard T2-weighted magnetic resonance cholangiography (MRC) sequences. Potential applications include the evaluation of leaks and strictures; the assessment of drainage in normal, surgically by-passed, stented and obstructed biliary systems; the diagnosis of cholecystitis; and the evaluation of normal and variant biliary anatomy.     

Sturge-Weber syndrome: diffusion magnetic resonance imaging and proton magnetic resonance spectroscopy findings

We report on the diffusion magnetic resonance imaging (MRI) and proton MR spectroscopy findings of a 26-year-old female patient with Sturge-Weber syndrome. Echo-planar trace diffusion MRI revealed mildly high signal intensity changes at parieto-occipital lobes on b = 1000 s/mm² images, suggesting restricted diffusion. On corresponding apparent diffusion coefficient maps, those areas had moderately high signal intensity and high apparent diffusion coefficient values (around 0.9×10(-3) mm²/s) compared with the contralateral symmetrical normal side of the brain (0.776×10(-3) mm²/s). This finding was consistent with increased motion of water molecules (disintegration of the neural tissue) in these regions. Proton MR spectroscopy revealed decreased N-acetyl aspartate and increased choline peaks, indicating disintegration of neural tissue associated with neuronal loss as well.     

Noise suppression digital filter for functional magnetic resonance imaging based on image reference data

The central decision in every functional magnetic resonance imaging (fMRI) experiment is whether pixels in brain tissues are showing activation in response to neural stimulus or as a result of noise. Images are degraded not only by random (e.g., thermal) noise, but also by structured noise due to MR system characteristics, cardiac and respiratory pulsations, and patient motion. A novel digital filter has been developed to suppress cardiac and respiratory structured noise in fMRI images, using estimates of structured and random noise power spectra obtained directly from the images. It is an adaptive filter based on stationary noise statistics, and is equivalent in form to a Wiener filter. A mathematical model of the filtering process was developed to understand how the strength and distribution of structured and random noise power influenced filter performance. The filter was tested using images from an auditory activation study in ten subjects. In subjects whose structured noise power was localized to a relatively narrow frequency range, a strong relationship was found, both experimentally (R = 0.975, P < 0.0004 for H_o: R = 0) and using the model, between filter performance and the level of structured noise power contaminating the experiment frequency. The filter significantly reduced the rate of false-positive activations in the subset of subjects whose experiment frequency was relatively heavily contaminated by structured noise. Notch filters, that simply eliminate unwanted frequencies, performed poorly in all subjects. Unlike the proposed Wiener filter, these filters did not suppress structured noise power at the experiment frequency that contributes to false-positive activations.     

Whole-heart steady-state free precession coronary artery magnetic resonance angiography.

Current implementations of coronary artery magnetic resonance angiography (MRA) suffer from limited coverage of the coronary arterial system. Whole-heart coronary MRA was implemented based on a free-breathing steady-state free-precession (SSFP) technique with magnetization preparation. The technique was compared to a similar implementation of conventional, thin-slab coronary MRA in 12 normal volunteers. Three thin-slab volumes were prescribed: 1) a transverse slab, covering the left main (LM) artery and proximal segments of the left anterior ascending (LAD) and left circumflex (LCX) coronary arteries; 2) a double-oblique slab covering the right coronary artery (RCA); and 3) a double-oblique slab covering the proximal and distal segments of the LCX. The whole-heart data set was reformatted in identical orientations. Visible vessel length, vessel sharpness, and vessel diameter were determined and compared separately for each vessel. Whole-heart coronary MRA visualized LM/LAD (11.7 +/- 3.4 cm) and LCX (6.9 +/- 3.6 cm) over a significantly longer distance than the transverse volume (LM/LAD, 6.1 +/- 1.1 cm, P < 0.001; LCX, 4.2 +/- 1.2 cm, P < 0.05). Improvements in visible vessel length for RCA and LCX in the whole-heart approach vs. their respective targeted volumes were not significant. It is concluded that the whole-heart coronary MRA technique improves visible vessel length and facilitates high-quality coronary MRA of the complete coronary artery tree in a single measurement.     

Functional magnetic resonance imaging for neurosurgical planning in neurooncology

Functional magnetic resonance imaging (fMRI) is a non-invasive technique that is widely available and can be used to determine the spatial relationships between tumor tissue and eloquent brain areas. Within certain limits, this functional information can be applied in the field of neurosurgery as a pre-operative mapping tool to minimize damage to eloquent brain areas. In this article, we review the literature on the use of fMRI for neurosurgical planning. The issues addressed are: (1) stimulation paradigms, (2) the influence of tumors on the blood oxygenation level-dependent (BOLD) signal, (3) post-processing the fMRI time course, (4) integration of fMRI results into neuronavigation systems, (5) the accuracy of fMRI and (6) fMRI compared to intra-operative mapping (IOM).