Functional Multiparametric Magnetic Resonance Imaging of the Kidneys Using Blood Oxygen Level-Dependent and Diffusion-Weighted Sequences: a Reliable Tool for Monitoring Acute Upper Urinary Tract Obstruction.

PURPOSE Little data is available on noninvasive MRI-based assessment of renal function during upper urinary tract (UUT) obstruction. In this study, we determined whether functional multiparametric kidney MRI is able to monitor treatment response in acute unilateral UUT obstruction. MATERIAL AND METHODS Between 01/2008 and 01/2010, 18 patients with acute unilateral UUT obstruction due to calculi were prospectively enrolled to undergo kidney MRI with conventional, blood oxygen level-dependent (BOLD) and diffusion-weighted (DW) sequences on emergency admission and after release of obstruction. Functional imaging parameters of the obstructed and contralateral unobstructed kidneys derived from BOLD (apparent spin relaxation rate [R2*]) and DW (total apparent diffusion coefficient [ADCT], pure diffusion coefficient [ADCD] and perfusion fraction [FP]) sequences were assessed during acute UUT obstruction and after its release. RESULTS During acute obstruction, R2* and FP values were lower in the cortex (p=0.020 and p=0.031, respectively) and medulla (p=0.012 and p=0.190, respectively) of the obstructed compared to the contralateral unobstructed kidneys. After release of obstruction, R2* and FP values increased both in the cortex (p=0.016 and p=0.004, respectively) and medulla (p=0.071 and p=0.044, respectively) of the formerly obstructed kidneys to values similar to those found in the contralateral kidneys. ADCT and ADCD values did not significantly differ between obstructed and contralateral unobstructed kidneys during or after obstruction. CONCLUSIONS In our patients with acute unilateral UUT obstruction due to calculi, functional kidney MRI using BOLD and DW sequences allowed for the monitoring of pathophysiologic changes of obstructed kidneys during obstruction and after its release.

Quantitative MRI of spinal cord injury in a rat model: Correlative studies

Serial quantitative and correlative studies of experimental spinal cord injury (SCI) in rats were conducted using three-dimensional magnetic resonance imaging (MRI). Correlative measures included morphological histopathology, neurobehavioral measures of functional deficit, and biochemical assays for N-acetyl-aspartate (NAA), lactate, pyruvate, and ATP. A spinal cord injury device was characterized and provided a reproducible injury severity. Injuries were moderate and consistent to within $\pm$20% (standard deviation). For MRI, a three-dimensional implementation of the single spin-echo FATE (Fast optimum angle, short TE) pulse sequence was used for rapid acquisition, with a 128 x 128 x 32 (x,y,z) matrix size and a 0.21 x 0.21 x 1.5 mm resolution. These serial studies revealed a bimodal characteristic in the evolution in MRI pathology with time. Early and late phases of SCI pathology were clearly visualized in $T\sb2$-weighted MRI, and these corresponded to specific histopathological changes in the spinal cord. Centralized hypointense MRI regions correlated with evidence of hemorrhagic and necrotic tissue, while surrounding hyperintense regions represented edema or myelomalacia. Unexpectedly, $T\sb2$-weighted MRI pathology contrast at 24 hours after injury appeared to subside before peaking at 72 hours after injury. This change is likely attributable to ongoing secondary injury processes, which may alter local $T\sb2$ values or reduce the natural anisotropy of the spinal cord. MRI, functional, and histological measures all indicated that 72 hours after injury was the temporal maximum for quantitative measures of spinal cord pathology. Thereafter, significant improvement was seen only in neurobehavioral scores. Significant correlations were found between quantitated MRI pathology and histopathology. Also, NAA and lactate levels correlated with behavioral measures of the level of function deficit. Asymmetric (rostral/caudal) changes in NAA and lactate due to injury indicate that rostral and caudal segments from the injury site are affected differently by the injury. These studies indicate that volumetric quantitation of MRI pathology from $T\sb2$-weighted images may play an important role in early prediction of neurologic deficit and spinal cord pathology. The loss of $T\sb2$ contrast at 24 hours suggests MR may be able to detect certain delayed mechanisms of secondary injury which are not resolved by histopathology or other radiological modalities. Furthermore, in vivo proton magnetic resonance spectroscopy (MRS) studies of SCI may provide a valuable addition source of information about changes in regional spinal cord lactate and NAA levels, which are indicative of local metabolic and pathological changes. ^

Functional connectivity in BOLD and CBF data: similarity and reliability of resting brain networks

Resting-state functional connectivity (FC) fMRI (rs-fcMRI) offers an appealing approach to mapping the brain's intrinsic functional organization. Blood oxygen level dependent (BOLD) and arterial spin labeling (ASL) are the two main rs-fcMRI approaches to assess alterations in brain networks associated with individual differences, behavior and psychopathology. While the BOLD signal is stronger with a higher temporal resolution, ASL provides quantitative, direct measures of the physiology and metabolism of specific networks. This study systematically investigated the similarity and reliability of resting brain networks (RBNs) in BOLD and ASL. A 2×2×2 factorial design was employed where each subject underwent repeated BOLD and ASL rs-fcMRI scans on two occasions on two MRI scanners respectively. Both independent and joint FC analyses revealed common RBNs in ASL and BOLD rs-fcMRI with a moderate to high level of spatial overlap, verified by Dice Similarity Coefficients. Test-retest analyses indicated more reliable spatial network patterns in BOLD (average modal Intraclass Correlation Coefficients: 0.905±0.033 between-sessions; 0.885±0.052 between-scanners) than ASL (0.545±0.048; 0.575±0.059). Nevertheless, ASL provided highly reproducible (0.955±0.021; 0.970±0.011) network-specific CBF measurements. Moreover, we observed positive correlations between regional CBF and FC in core areas of all RBNs indicating a relationship between network connectivity and its baseline metabolism. Taken together, the combination of ASL and BOLD rs-fcMRI provides a powerful tool for characterizing the spatiotemporal and quantitative properties of RBNs. These findings pave the way for future BOLD and ASL rs-fcMRI studies in clinical populations that are carried out across time and scanners.

Comparison between magnetic resonance imaging and computed tomography of the lung in patients with cystic fibrosis with regard to clinical, laboratory, and pulmonary functional parameters

OBJECTIVE To evaluate whether magnetic resonance imaging (MRI) is effective as computed tomography (CT) in determining morphologic and functional pulmonary changes in patients with cystic fibrosis (CF) in association with multiple clinical parameters. MATERIALS AND METHODS Institutional review board approval and patient written informed consent were obtained. In this prospective study, 30 patients with CF (17 men and 13 women; mean (SD) age, 30.2 (9.2) years; range, 19-52 years) were included. Chest CT was acquired by unenhanced low-dose technique for clinical purposes. Lung MRI (1.5 T) comprised T2- and T1-weighted sequences before and after the application of 0.1-mmol·kg gadobutrol, also considering lung perfusion imaging. All CT and MR images were visually evaluated by using 2 different scoring systems: the modified Helbich and the Eichinger scores. Signal intensity of the peribronchial walls and detected mucus on T2-weighted images as well as signal enhancement of the peribronchial walls on contrast-enhanced T1-weighted sequences were additionally assessed on MRI. For the clinical evaluation, the pulmonary exacerbation rate, laboratory, and pulmonary functional parameters were determined. RESULTS The overall modified Helbich CT score had a mean (SD) of 15.3 (4.8) (range, 3-21) and median of 16.0 (interquartile range [IQR], 6.3). The overall modified Helbich MR score showed slightly, not significantly, lower values (Wilcoxon rank sum test and Student t test; P > 0.05): mean (SD) of 14.3 (4.7) (range, 3-20) and median of 15.0 (IQR, 7.3). Without assessment of perfusion, the overall Eichinger score resulted in the following values for CT vs MR examinations: mean (SD), 20.3 (7.2) (range, 4-31); and median, 21.0 (IQR, 9.5) vs mean (SD), 19.5 (7.1) (range, 4-33); and median, 20.0 (IQR, 9.0). All differences between CT and MR examinations were not significant (Wilcoxon rank sum tests and Student t tests; P > 0.05). In general, the correlations of the CT scores (overall and different imaging parameters) to the clinical parameters were slightly higher compared to the MRI scores. However, if all additional MRI parameters were integrated into the scoring systems, the correlations reached the values of the CT scores. The overall image quality was significantly higher for the CT examinations compared to the MRI sequences. CONCLUSIONS One major diagnostic benefit of lung MRI in CF is the possible acquisition of several different morphologic and functional imaging features without the use of any radiation exposure. Lung MRI shows reliable associations with CT and clinical parameters, which suggests its implementation in CF for routine diagnosis, which would be particularly important in follow-up imaging over the long term.

Investigation of arterial spin labeling MRI for quantitative cerebral blood flow measurement

Arterial spin labeling (ASL) is a technique for noninvasively measuring cerebral perfusion using magnetic resonance imaging. Clinical applications of ASL include functional activation studies, evaluation of the effect of pharmaceuticals on perfusion, and assessment of cerebrovascular disease, stroke, and brain tumor. The use of ASL in the clinic has been limited by poor image quality when large anatomic coverage is required and the time required for data acquisition and processing. This research sought to address these difficulties by optimizing the ASL acquisition and processing schemes. To improve data acquisition, optimal acquisition parameters were determined through simulations, phantom studies and in vivo measurements. The scan time for ASL data acquisition was limited to fifteen minutes to reduce potential subject motion. A processing scheme was implemented that rapidly produced regional cerebral blood flow (rCBF) maps with minimal user input. To provide a measure of the precision of the rCBF values produced by ASL, bootstrap analysis was performed on a representative data set. The bootstrap analysis of single gray and white matter voxels yielded a coefficient of variation of 6.7% and 29% respectively, implying that the calculated rCBF value is far more precise for gray matter than white matter. Additionally, bootstrap analysis was performed to investigate the sensitivity of the rCBF data to the input parameters and provide a quantitative comparison of several existing perfusion models. This study guided the selection of the optimum perfusion quantification model for further experiments. The optimized ASL acquisition and processing schemes were evaluated with two ASL acquisitions on each of five normal subjects. The gray-to-white matter rCBF ratios for nine of the ten acquisitions were within ±10% of 2.6 and none were statistically different from 2.6, the typical ratio produced by a variety of quantitative perfusion techniques. Overall, this work produced an ASL data acquisition and processing technique for quantitative perfusion and functional activation studies, while revealing the limitations of the technique through bootstrap analysis. ^

Objective Assessment of Olfactory Function Using Functional Magnetic Resonance Imaging=Estudio objetivo del olfato mediante resonancia magnética funcional

Objective: To show the results of a device that generates automated olfactory stimuli suitable for functional magnetic resonance imaging (fMRI) experiments. Material and methods: Te n normal volunteers, 5 women and 5 men, were studied. The system allows the programming of several sequences, providing the capability to synchronise the onset of odour presentation with acquisition by a trigger signal of the MRI scanner. The olfactometer is a device that allows selection of the odour, the event paradigm, the time of stimuli and the odour concentration. The paradigm used during fMRI scanning consisted of 15-s blocks. The odorant event took 2 s with butanol, mint and coffee. Results: We observed olfactory activity in the olfactory bulb, entorhinal cortex (4%), amygdala (2.5%) and temporo-parietal cortex, especially in the areas related to emotional integration. Conclusions: The device has demonstrated its effectiveness in stimulating olfactory areas and its capacity to adapt to fMRI equipment.RESUMEN Objetivo: Mostrar los resultados del olfatómetro capaz de generar tareas olfativas en un equipo de resonancia magnética funcional (fMRI). Material y métodos: Estudiamos 10 sujetos normales: 5 varones y 5 mujeres. El olfatómetro está dise ̃ nado para que el estímulo que produce se sincronice con el equipo de fMRI mediante la se ̃ nal desencadenante que suministra el propio equipo. El olfatómetro es capaz de: selec- cionar el olor, secuenciar los distintos olores, programar la frecuencia y duración de los olores y controlar la intensidad del olor. El paradigma utilizado responde a un dise ̃ no de activación asociada a eventos, en el que la duración del bloque de activación y de reposo es de 15 s. La duración del estímulo olfativo (butanol, menta o café) es de 2 segundos, durante toda la serie que consta de 9 ciclos. Resultados: Se ha observado reactividad (contraste BOLD) en las diferentes áreas cerebrales involucradas en las tareas olfativas: bulbo olfatorio, córtex entorrinal (4%), amigdala (2,5%) y córtex temporoparietal. Las áreas relacionadas con integración de las emociones tienen una reactividad mayor. Conclusiones: El dispositivo propuesto nos permite controlar de forma automática y sincronizada los olores necesarios para estudiar la actividad de las áreas olfatorias cerebrales mediante fMRI.

Disparate connectivity for structural and functional networks is revealed when physical location of the connected nodes is considered

Macroscopic brain networks have been widely described with the manifold of metrics available using graph theory. However, most analyses do not incorporate information about the physical position of network nodes. Here, we provide a multimodal macroscopic network characterization while considering the physical positions of nodes. To do so, we examined anatomical and functional macroscopic brain networks in a sample of twenty healthy subjects. Anatomical networks are obtained with a graph based tractography algorithm from diffusion-weighted magnetic resonance images (DW-MRI). Anatomical con- nections identified via DW-MRI provided probabilistic constraints for determining the connectedness of 90 dif- ferent brain areas. Functional networks are derived from temporal linear correlations between blood-oxygenation level-dependent signals derived from the same brain areas. Rentian Scaling analysis, a technique adapted from very- large-scale integration circuits analyses, shows that func- tional networks are more random and less optimized than the anatomical networks. We also provide a new metric that allows quantifying the global connectivity arrange- ments for both structural and functional networks. While the functional networks show a higher contribution of inter-hemispheric connections, the anatomical networks highest connections are identified in a dorsal?ventral arrangement. These results indicate that anatomical and functional networks present different connectivity organi- zations that can only be identified when the physical locations of the nodes are included in the analysis.

Guiding functional connectivity estimation by structural connectivity in MEG: an application to discrimination of mild cognitive impaired conditions

Whole brain resting state connectivity is a promising biomarker that might help to obtain an early diagnosis in many neurological diseases, such as dementia. Inferring resting-state connectivity is often based on correlations, which are sensitive to indirect connections, leading to an inaccurate representation of the real backbone of the network. The precision matrix is a better representation for whole brain connectivity, as it considers only direct connections. The network structure can be estimated using the graphical lasso (GL), which achieves sparsity through l1-regularization on the precision matrix. In this paper, we propose a structural connectivity adaptive version of the GL, where weaker anatomical connections are represented as stronger penalties on the corre- sponding functional connections. We applied beamformer source reconstruction to the resting state MEG record- ings of 81 subjects, where 29 were healthy controls, 22 were single-domain amnestic Mild Cognitive Impaired (MCI), and 30 were multiple-domain amnestic MCI. An atlas-based anatomical parcellation of 66 regions was ob- tained for each subject, and time series were assigned to each of the regions. The fiber densities between the re- gions, obtained with deterministic tractography from diffusion-weighted MRI, were used to define the anatomical connectivity. Precision matrices were obtained with the region specific time series in five different frequency bands. We compared our method with the traditional GL and a functional adaptive version of the GL, in terms of log-likelihood and classification accuracies between the three groups. We conclude that introduc- ing an anatomical prior improves the expressivity of the model and, in most cases, leads to a better classification between groups.

Imaging of H217O distribution in the brain of a live rat by using proton-detected 17O MRI

Imaging of H217O has a number of important applications. Mapping the distribution of H217O produced by oxidative metabolism of 17O-enriched oxygen gas may lead to a new method of metabolic functional imaging; regional cerebral blood flow also can be measured by measuring the H217O distribution after the injection of 17O-enriched physiological saline solution. Previous studies have proposed a method for indirect detection of 17O. The method is based on the shortening of the proton T2 in H217O solutions, caused by the residual 17O-1H scalar coupling and transferred to the bulk water via fast chemical exchange. It has been shown that the proton T2 of H217O solutions can be restored to that of H216O by irradiating the resonance frequency of the 17O nucleus. The indirect 17O image thus is obtained by taking the difference between two T2-weighted spin-echo images: one acquired after irradiation of the 17O resonance and one acquired without irradiation. It also has been established that, at relatively low concentrations of H217O, the indirect method yields an image that quantitatively reflects the H217O distribution in the sample. The method is referred to as PRIMO (proton imaging of oxygen). In this work, we show in vivo proton images of the H217O distribution in a rat brain after an i.v. injection of H217O-enriched physiological saline solution. Implementing the indirect detection method in an echo-planar imaging sequence enabled obtaining H217O images with good spatial and temporal resolution of few seconds.

Myopia: structural and functional correlates

Ocular dimensions are widely recognised as key variants of refractive error. Previously, accurate depiction of eye shape in vivo was largely restricted by limitations in the imaging techniques available. This thesis describes unique applications of the recently introduced 3-dimensional magnetic resonance imaging (MRI) approach to evaluate human eye shape in a group of young adult subjects (n=76) with a range of ametropia (MSE= -19.76 to +4.38D). Specific MRI derived parameters of ocular shape are then correlated with measures of visual function. Key findings include the significant homogeneity of ocular volume in the anterior eye for a range of refractive errors, whilst significant volume changes occur in the posterior eye as a function of ametropia. Anterior vs. posterior eye differences have also been shown through evaluations of equivalent spherical radius; the posterior 25% cap of the eye was shown to be relatively steeper in myopes compared to emmetropes. Further analyses showed differences in retinal quadrant profiles; assessments of the maximum distance from the retinal surface to the presumed visual axes showed exaggerated growth of the temporal quadrant in myopic eyes. Comparisons of retinal contour values derived from transformation of peripheral refraction data were made with MRI; flatter retinal curvature values were noted when using the MRI technique. A distinctive feature of this work is the evaluation of the relationship between ocular structure and visual function. Multiple aspects of visual function were evaluated through several vehicles: multifocal electroretinogram testing, visual field sensitivity testing, and the use of psychophysical methods to determine ganglion cell density. The results show that many quadrantic structural and functional variations exist. In general, the data could not demonstrate a significant correlation between visual function and associated measures of ocular conformation either within or between myopic and emmetropic groups.

Functional localisation of human sensory-motor cortex using magnetoencephalography

The 19 channel Neuromagnetometer system in the Clinical Neurophysiology Unit at Aston University is a multi-channel system, unique in the United Kingdom. A bite bar head localisation and MRI co-registration strategy which enabled accurate and reproducible localisation of MEG data into cortical space was developed. This afforded the opportunity to study magnetic fields of the human cortex generated by stimulation of peripheral nerve, by stimulation of visceral sensory receptors and by those evoked through voluntary finger movement. Initially, a study of sensory-motor evoked data was performed in a healthy control population. The techniques developed were then applied to patients who were to undergo neurosurgical intervention for the treatment of epilepsy and I or space occupying lesions. This enabled both validation of the effective accuracy of source localisation using MEG as well as to determine the clinical value of MEG in presurgical assessment of functional localisation in human cortex. The studies in this thesis have demonstrated that MEG can repeatedly and reliably locate sources contained within a single gyrus and thus potentially differentiate between disparate gyral activation. This ability is critical in the clinical application of any functional imaging technique; which is yet to be fully validated by any other 'non-invasive' functional imaging methodology. The technique was also applied to the study of visceral sensory representation in the cortex which yielded important data about the multiple cortical representation of visceral sensory function.

Universality in the merging dynamics of parametric active contours:a study in MRI based lung segmentation

Measurement of lung ventilation is one of the most reliable techniques in diagnosing pulmonary diseases. The time-consuming and bias-prone traditional methods using hyperpolarized H 3He and 1H magnetic resonance imageries have recently been improved by an automated technique based on 'multiple active contour evolution'. This method involves a simultaneous evolution of multiple initial conditions, called 'snakes', eventually leading to their 'merging' and is entirely independent of the shapes and sizes of snakes or other parametric details. The objective of this paper is to show, through a theoretical analysis, that the functional dynamics of merging as depicted in the active contour method has a direct analogue in statistical physics and this explains its 'universality'. We show that the multiple active contour method has an universal scaling behaviour akin to that of classical nucleation in two spatial dimensions. We prove our point by comparing the numerically evaluated exponents with an equivalent thermodynamic model. © IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.

An integrated functional mapping and source localization platform for brain research

The premise of this dissertation is to create a highly integrated platform that combines the most current recording technologies for brain research through the development of new algorithms for three-dimensional (3D) functional mapping and 3D source localization. The recording modalities that were integrated include: Electroencephalography (EEG), Optical Topographic Maps (OTM), Magnetic Resonance Imaging (MRI), and Diffusion Tensor Imaging (DTI). This work can be divided into two parts: The first part involves the integration of OTM with MRI, where the topographic maps are mapped to both the skull and cortical surface of the brain. This integration process is made possible through the development of new algorithms that determine the probes location on the MRI head model and warping the 2D topographic maps onto the 3D MRI head/brain model. Dynamic changes of the brain activation can be visualized on the MRI head model through a graphical user interface. The second part of this research involves augmenting a fiber tracking system, by adding the ability to integrate the source localization results generated by commercial software named Curry. This task involved registering the EEG electrodes and the dipole results to the MRI data. Such Integration will allow the visualization of fiber tracts, along with the source of the EEG, in a 3D transparent brain structure. The research findings of this dissertation were tested and validated through the participation of patients from Miami Children Hospital (MCH). Such an integrated platform presented to the medical professionals in the form of a user-friendly graphical interface is viewed as a major contribution of this dissertation. It should be emphasized that there are two main aspects to this research endeavor: (1) if a dipole could be situated in time at its different positions, its trajectory may reveal additional information on the extent and nature of the brain malfunction; (2) situating such a dipole trajectory with respect to the fiber tracks could ensure the preservation of these fiber tracks (axons) during surgical interventions, preserving as a consequence these parts of the brain that are responsible for information transmission.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.