Improved Techniques for Acquisition and Analysis of Dynamic Contrast-Enhanced Magnetic Resonance Imaging for Detecting Vascular Permeability in the Central Nervous System

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is a noninvasive technique for quantitative assessment of the integrity of blood-brain barrier and blood-spinal cord barrier (BSCB) in the presence of central nervous system pathologies. However, the results of DCE-MRI show substantial variability. The high variability can be caused by a number of factors including inaccurate T1 estimation, insufficient temporal resolution and poor contrast-to-noise ratio. My thesis work is to develop improved methods to reduce the variability of DCE-MRI results. To obtain fast and accurate T1 map, the Look-Locker acquisition technique was implemented with a novel and truly centric k-space segmentation scheme. In addition, an original multi-step curve fitting procedure was developed to increase the accuracy of T1 estimation. A view sharing acquisition method was implemented to increase temporal resolution, and a novel normalization method was introduced to reduce image artifacts. Finally, a new clustering algorithm was developed to reduce apparent noise in the DCE-MRI data. The performance of these proposed methods was verified by simulations and phantom studies. As part of this work, the proposed techniques were applied to an in vivo DCE-MRI study of experimental spinal cord injury (SCI). These methods have shown robust results and allow quantitative assessment of regions with very low vascular permeability. In conclusion, applications of the improved DCE-MRI acquisition and analysis methods developed in this thesis work can improve the accuracy of the DCE-MRI results.

Signal changes in the spinal cord of the rat after injection of formalin into the hindpaw: Characterization using functional magnetic resonance imaging

Changes in metabolism and local circulation occur in the spinal cord during peripheral noxious stimulation. Evidence is presented that this stimulation also causes signal intensity alterations in functional magnetic resonance images of the spinal cord during formalin-induced pain. These results indicate the potential of functional magnetic resonance imaging in assessing noninvasively the extent and intensity of spinal cord excitation in this well characterized pain model. Therefore, the aim of this study was to establish functional magnetic resonance imaging as a noninvasive method to characterize temporal changes in the spinal cord after a single injection of 50 μl of formalin subcutaneously into the hindpaw of the anesthetized rat. This challenge produced a biphasic licking activity in the freely moving conscious animal. Images of the spinal cord were acquired within 2 min, enabling monitoring of the site and the temporal evolution of the signal changes during the development of formalin-induced hyperalgesia without the need of any surgical procedure. The time course of changes in the spinal cord functional image in the isoflurane-anesthetized animal was similar to that obtained from behavioral experiments. Also, comparable physiological data, control experiments, and the inhibition of a response through application of the local anesthetic agent lidocaine indicate that the signal changes observed after formalin injection were specifically related to excitability changes in the relevant segments of the lumbar spinal cord. This approach could be useful to characterize different models of pain and hyperalgesia and, more importantly, to evaluate effects of analgesic drugs.

Dynamic contrast-enhanced magnetic resonance imaging reveals stress-induced angiogenesis in MCF7 human breast tumors.

The mechanism of contrast enhancement of tumors using magnetic resonance imaging was investigated in MCF7 human breast cancer implanted in nude mice. Dynamic contrast-enhanced images recorded at high spatial resolution were analyzed by an image analysis method based on a physiological model, which included the blood circulation, the tumor, the remaining tissues, and clearance via the kidneys. This analysis enabled us to map in rapidly enhancing regions within the tumor, the capillary permeability factor (capillary permeability times surface area per voxel volume) and the fraction of leakage space. Correlation of these maps with T2-weighted spin echo images, with histopathology, and with immunohistochemical staining of endothelial cells demonstrated the presence of dense permeable microcapillaries in the tumor periphery and in intratumoral regions that surrounded necrotic loci. The high leakage from the intratumoral permeable capillaries indicated an induction of a specific angiogenic process associated with stress conditions that cause necrosis. This induction was augmented in tumors responding to tamoxifen treatment. Determination of the distribution and extent of this stress-induced angiogenic activity by contrast-enhanced MRI might be of diagnostic and of prognostic value.

Magnetic resonance imaging (MRI) detection of the murine brain response to light: temporal differentiation and negative functional MRI changes.

Using a 9.4 T MRI instrument, we have obtained images of the mouse brain response to photic stimulation during a period between deep anesthesia and the early stages of arousal. The large image enhancements we observe (often >30%) are consistent with literature results extrapolated to 9.4 T. However, there are also two unusual aspects to our findings. (i) The visual area of the brain responds only to changes in stimulus intensity, suggesting that we directly detect operations of the M visual system pathway. Such a channel has been observed in mice by invasive electrophysiology, and described in detail for primates. (ii) Along with the typical positive response in the area of the occipital portion of the brain containing the visual cortex, another area displays decreased signal intensity upon stimulation.

Activation of single whisker barrel in rat brain localized by functional magnetic resonance imaging.

The previously established cortical representation of rat whiskers in layer IV of the cortex contains distinct cylindrical columns of cellular aggregates, which are termed barrels and correlate in a one-to-one relation to whiskers on the contralateral rat face. In the present study, functional magnetic resonance imaging (fMRI) of the rat brain was used to map whisker barrel activation during mechanical up-down movement (+/- 2.5 mm amplitude at 8 Hz) of single/multiple whisker(s). Multislice gradient echo fMRI experiments were performed at 7 T with in-plane image resolution of 220 x 220 microns, slice thickness of 1 mm, and echo time of 16 ms. Highly significant (P < 0.001) and localized contralateral regions of activation were observed upon stimulation of single/multiple whisker(s). In all experiments (n = 10), the locations of activation relative to bregma and midline were highly correlated with the neuroanatomical position of the corresponding whisker barrels, and the results were reproducible intra- and interanimal. Our results indicate that fMRI based on blood oxygenation level-dependent image contrast has the sensitivity to depict activation of a single whisker barrel in the rat brain. This noninvasive technique will supplement existing methods in the study of rat barrel cortex and should be particularly useful for the long-term investigations of central nervous system in the same animal.

Identification of human brain regions underlying responses to resistive inspiratory loading with functional magnetic resonance imaging.

Compensatory ventilatory responses to increased inspiratory loading are essential for adequate breathing regulation in a number of pulmonary diseases; however, the human brain sites mediating such responses are unknown. Midsagittal and axial images were acquired in 11 healthy volunteers during unloaded and loaded (30 cmH2O; 1 cmH2O = 98 Pa) inspiratory breathing, by using functional magnetic resonance imaging (fMRI) strategies (1.5-tesla MR; repetition time, 72 msec; echo time, 45 msec; flip angle, 30 degrees; field of view, 26 cm; slice thickness, 5 mm; number of excitations, 1; matrix, 128 x 256). Digital image subtractions and region of interest analyses revealed significantly increased fMRI signal intensity in discrete areas of the ventral and dorsal pons, interpeduncular nucleus, basal forebrain, putamen, and cerebellar regions. Upon load withdrawal, certain regions displayed a rapid fMRI signal off-transient, while in others, a slower fMRI signal decay emerged. Sustained loading elicited slow decreases in fMRI signal across activated regions, while second application of an identical load resulted in smaller signal increases compared to initial signal responses (P < 0.001). A moderate inspiratory load is associated with consistent regional activation of discrete brain locations; certain of these regions have been implicated in mediation of loaded breathing in animal models. We speculate that temporal changes in fMRI signal may indicate respiratory after-discharge and/or habituation phenomena.

Speeding up the EM algorithm for mixture model-based segmentation of magnetic resonance images

Mixture models implemented via the expectation-maximization (EM) algorithm are being increasingly used in a wide range of problems in pattern recognition such as image segmentation. However, the EM algorithm requires considerable computational time in its application to huge data sets such as a three-dimensional magnetic resonance (MR) image of over 10 million voxels. Recently, it was shown that a sparse, incremental version of the EM algorithm could improve its rate of convergence. In this paper, we show how this modified EM algorithm can be speeded up further by adopting a multiresolution kd-tree structure in performing the E-step. The proposed algorithm outperforms some other variants of the EM algorithm for segmenting MR images of the human brain. (C) 2004 Pattern Recognition Society. Published by Elsevier Ltd. All rights reserved.

Impact of scar thickness on the assessment of viability using dobutamine echocardiography and thallium single-photon emission computed tomography - A comparison with contrast-enhanced magnetic resonance imaging

OBJECTIVES We sought to determine whether the transmural extent of scar (TES) explains discordances between dobutamine echocardiography (DbE) and thallium single-photon emission computed tomography (Tl-SPECT) in the detection of viable myocardium (VM). BACKGROUND Discrepancies between DbE and Tl-SPECT are often attributed to differences between contractile reserve and membrane integrity, but may also reflect a disproportionate influence of nontransmural scar on thickening at DbE. METHODS Sixty patients (age 62 +/- 12 years; 10 women and 50 men) with postinfarction left ventricular dysfunction underwent standard rest-late redistribution Tl-SPECT and DbE. Viable myocardium was identified when dysfunctional segments showed Tl activity >60% on the late-redistribution image or by low-dose augmentation at DbE. Contrast-enhanced magnetic resonance imaging (ceMRI) was used to divide TES into five groups: 0%, 75% of the wall thickness replaced by scar. RESULTS As TES increased, both the mean Tl uptake and change in wall motion score decreased significantly (both p < 0.001). However, the presence of subendocardial scar was insufficient to prevent thickening; >50% of segments still showed contractile function with TES of 25% to 75%, although residual function was uncommon with TES >75%. The relationship of both tests to increasing TES was similar, but Tl-SPECT identified VM more frequently than DbE in all groups. Among segments without scar or with small amounts of scar (50% were viable by SPECT. CONCLUSIONS Both contractile reserve and perfusion are sensitive to the extent of scar. However, contractile reserve may be impaired in the face of no or minor scar, and thickening may still occur with extensive scar. (C) 2004 by the American College of Cardiology Foundation.

A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging

A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings. especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging. (C) 2004 Elsevier Inc. All rights reserved.

3-D magnetic resonance angiography versus conventional angiography in peripheral arterial disease: Pilot study

Angiography is usually performed as the preoperative road map for those requiring revascularization for lower extremity peripheral arterial disease (PAD). The alternative investigations are ultrasound, 3-D magnetic resonance angiography (3-D MRA) and computed tomography angiography. This pilot study aimed to assess whether 3-D MRA could replace the gold standard angiography in preoperative planning. Eight patients considered for aortoiliac or infrainguinal arterial bypass surgery were recruited. All underwent both imaging modalities within 7 days. A vascular surgeon and a radiologist each reported on the images from both the 3-D MRA and the angiography, with blinding to patient details and each others reports. Comparisons were made between the reports for the angiographic and the 3-D MRA images, and between the reports of the vascular surgeon and the radiologist. Compared to the gold standard angiogram, 3-D MRA had a sensitivity of 77% and specificity of 94% in detecting occlusion, and a sensitivity of 72% and specificity of 90% in differentiating high grade (> 50%) versus low grade (< 50%) stenoses. There was an overall concordance of 78% between the two investigations with a range of 62% in the peroneal artery to 94% in the aorta. 3-D MRA showed flow in 23% of cases where conventional angiography showed no flow. In the present pilot study, 3-D MRA had reasonable concordance with the gold standard angiography, depending on the level of the lesion. At times it showed vessel flow where occlusion was shown on conventional angiogram. 3-D MRA in peripheral vascular disease is challenging the gold standard, but is inconsistent at present.

Noninvasive localization of petroleum-derived spray oil in plants with chemical shift selective magnetic resonance imaging

Nuclear magnetic resonance (NMR) spectroscopy and magnetic resonance imaging (MRI) were used to detect petroleum-derived spray oils (PDSOs) in citrus seedlings and trees. The NMR spectrum of the phantom containing 10% (v/v) of a nC24 agricultural mineral oil (AMO) showed the resonance of the water protons at delta = 5 ppm, while the resonance of the oil protons at delta = 1.3 to 1.7 ppm. The peak resolution and the chemical shift difference of more than 3.3 ppm between water and oil protons effectively differentiated water and the oil. Chemical shift selective imaging (CSSI) was performed to localize the AMO within the stems of Citrus trifoliata L. seedlings after the application of a 4% (v/v) spray. The chemical shift selective images of the oil were acquired by excitation at delta = 1.5 ppm by averaging over 400 transients in each phase-encoding step. Oil was mainly detected in the outer cortex of stems within 10 d of spray application; some oil was also observed in the inner vascular bundle and pith of the stems at this point. CSSI was also applied to investigate the persistence of oil deposits in sprayed mature Washington navel orange (Citrus x aurantium L.) trees in an orchard. The trees were treated with either fourteen 0.25%, fourteen 0.5%, four 1.75%, or single 7% sprays of a nC23 horticultural mineral oil (HMO) 12 to 16 months before examination of plant tissues by CSSI, and were still showing symptoms of chronic phytotoxicity largely manifested as reduced yield. The oil deposits were detected in stems of sprayed flushes and unsprayed flushes produced 4 to 5 months after the last spray was applied, suggesting a potential movement of the oil via phloem and a correlation of the persistence of oil deposit in plants and the phytotoxicity. The results demonstrate that MRI is an effective method to probe the uptake and localization of PDSOs and other xenobiotics in vivo in plants noninvasively and nondestructively.

Magnetic resonance microscopy of the equine hoof wall: a study of resolution and potential

Reasons for performing study: Obtaining magnetic resonance images of the inner hoof wall tissue at the microscopic level would enable early accurate diagnosis of laminitis and therefore more effective therapy. Objectives: To optimise magnetic resonance imaging (MRI) parameters in order to obtain the highest possible resolution of the structures beneath the equine hoof wall. Methods: Magnetic resonance microscopy (MRM) was performed in front feet from 6 cadaver horses using T-2-weighted fast spin echo (FSE-T-2), and T-1-weighted gradient echo (GRE-T-1) sequences. Results: In T-2 weighted FSE images most of the stratum medium showed no signal, however the coronary, terminal and sole papillae were visible. The stratum lamellatum was clearly visible and primary epidermal lamellae could be differentiated from dermal lamellae. Conclusion: Most structures beneath the hoof wall were differentiated. Conventional scanners for diagnostic MRI in horses are low or high field. However this study used ultra-high field scanners currently not available for clinical use. Signal-to-noise ratio (SIN) increases as a function of field strength. An increase of spatial resolution of the image results in a decreased SIN. SIN can also be improved with better coils and the resolution of high field MRI scanners will increase as technology develops and surface array coils become more readily available. Potential relevance: Although MR images with microscopic resolution were obtained ex vivo, this study demonstrates the potential for detection of lamellar pathology as it occurs. Early recognition of the development of laminitis to instigate effective therapy at an earlier stage and may improve the outcome for laminitic horses. Clinical MR is now readily available at 3 T, while 4 T, 7 T and 9 T systems are being used for human whole body applications.

High-field magnetic resonance imaging with reduced field/tissue RF artefacts - A modeling study using hybrid MoM/FEM and FDTD technique

Due to complex field/tissue interactions, high-field magnetic resonance (MR) images suffer significant image distortions that result in compromised diagnostic quality. A new method that attempts to remove these distortions is proposed in this paper and is based on the use of transceiver-phased arrays. The proposed system uses, in the examples presented herein, a shielded four-element transceive-phased array head coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both the images together, the image distortion can be reduced several fold. A combined hybrid method of moments (MoM)/finite element method (FEM) and finite-difference time-domain (FDTD) technique is proposed and used to elucidate the concept of the new method and to accurately evaluate the electromagnetic field (EMF) in a human head model. In addition, the proposed method is used in conjunction with the generalized auto-calibrating partially parallel acquisitions (GRAPPA) reconstruction technique to enable rapid imaging of the two scans. Simulation results reported herein for 11-T (470-MHz) brain imaging applications show that the new method with GRAPPA reconstruction theoretically results in improved image quality and that the proposed combined hybrid MoM/FEM and FDTD technique is. suitable for high-field magnetic resonance imaging (MRI) numerical analysis.

Probing Tissue Microstructure Using Susceptibility Contrast Magnetic Resonance Imaging

Magnetic resonance imaging is a research and clinical tool that has been applied in a wide variety of sciences. One area of magnetic resonance imaging that has exhibited terrific promise and growth in the past decade is magnetic susceptibility imaging. Imaging tissue susceptibility provides insight into the microstructural organization and chemical properties of biological tissues, but this image contrast is not well understood. The purpose of this work is to develop effective approaches to image, assess, and model the mechanisms that generate both isotropic and anisotropic magnetic susceptibility contrast in biological tissues, including myocardium and central nervous system white matter.

This document contains the first report of MRI-measured susceptibility anisotropy in myocardium. Intact mouse heart specimens were scanned using MRI at 9.4 T to ascertain both the magnetic susceptibility and myofiber orientation of the tissue. The susceptibility anisotropy of myocardium was observed and measured by relating the apparent tissue susceptibility as a function of the myofiber angle with respect to the applied magnetic field. A multi-filament model of myocardial tissue revealed that the diamagnetically anisotropy α-helix peptide bonds in myofilament proteins are capable of producing bulk susceptibility anisotropy on a scale measurable by MRI, and are potentially the chief sources of the experimentally observed anisotropy.

The growing use of paramagnetic contrast agents in magnetic susceptibility imaging motivated a series of investigations regarding the effect of these exogenous agents on susceptibility imaging in the brain, heart, and kidney. In each of these organs, gadolinium increases susceptibility contrast and anisotropy, though the enhancements depend on the tissue type, compartmentalization of contrast agent, and complex multi-pool relaxation. In the brain, the introduction of paramagnetic contrast agents actually makes white matter tissue regions appear more diamagnetic relative to the reference susceptibility. Gadolinium-enhanced MRI yields tensor-valued susceptibility images with eigenvectors that more accurately reflect the underlying tissue orientation.

Despite the boost gadolinium provides, tensor-valued susceptibility image reconstruction is prone to image artifacts. A novel algorithm was developed to mitigate these artifacts by incorporating orientation-dependent tissue relaxation information into susceptibility tensor estimation. The technique was verified using a numerical phantom simulation, and improves susceptibility-based tractography in the brain, kidney, and heart. This work represents the first successful application of susceptibility-based tractography to a whole, intact heart.

The knowledge and tools developed throughout the course of this research were then applied to studying mouse models of Alzheimer’s disease in vivo, and studying hypertrophic human myocardium specimens ex vivo. Though a preliminary study using contrast-enhanced quantitative susceptibility mapping has revealed diamagnetic amyloid plaques associated with Alzheimer’s disease in the mouse brain ex vivo, non-contrast susceptibility imaging was unable to precisely identify these plaques in vivo. Susceptibility tensor imaging of human myocardium specimens at 9.4 T shows that susceptibility anisotropy is larger and mean susceptibility is more diamagnetic in hypertrophic tissue than in normal tissue. These findings support the hypothesis that myofilament proteins are a source of susceptibility contrast and anisotropy in myocardium. This collection of preclinical studies provides new tools and context for analyzing tissue structure, chemistry, and health in a variety of organs throughout the body.

Development and Optimization of Four-dimensional Magnetic Resonance Imaging (4D-MRI) for Radiation Therapy

A tenet of modern radiotherapy (RT) is to identify the treatment target accurately, following which the high-dose treatment volume may be expanded into the surrounding tissues in order to create the clinical and planning target volumes. Respiratory motion can induce errors in target volume delineation and dose delivery in radiation therapy for thoracic and abdominal cancers. Historically, radiotherapy treatment planning in the thoracic and abdominal regions has used 2D or 3D images acquired under uncoached free-breathing conditions, irrespective of whether the target tumor is moving or not. Once the gross target volume has been delineated, standard margins are commonly added in order to account for motion. However, the generic margins do not usually take the target motion trajectory into consideration. That may lead to under- or over-estimate motion with subsequent risk of missing the target during treatment or irradiating excessive normal tissue. That introduces systematic errors into treatment planning and delivery. In clinical practice, four-dimensional (4D) imaging has been popular in For RT motion management. It provides temporal information about tumor and organ at risk motion, and it permits patient-specific treatment planning. The most common contemporary imaging technique for identifying tumor motion is 4D computed tomography (4D-CT). However, CT has poor soft tissue contrast and it induce ionizing radiation hazard. In the last decade, 4D magnetic resonance imaging (4D-MRI) has become an emerging tool to image respiratory motion, especially in the abdomen, because of the superior soft-tissue contrast. Recently, several 4D-MRI techniques have been proposed, including prospective and retrospective approaches. Nevertheless, 4D-MRI techniques are faced with several challenges: 1) suboptimal and inconsistent tumor contrast with large inter-patient variation; 2) relatively low temporal-spatial resolution; 3) it lacks a reliable respiratory surrogate. In this research work, novel 4D-MRI techniques applying MRI weightings that was not used in existing 4D-MRI techniques, including T2/T1-weighted, T2-weighted and Diffusion-weighted MRI were investigated. A result-driven phase retrospective sorting method was proposed, and it was applied to image space as well as k-space of MR imaging. Novel image-based respiratory surrogates were developed, improved and evaluated.