<p>On-board image guidance, such as cone-beam CT (CBCT) and kV/MV 2D imaging, is essential in many radiation therapy procedures, such as intensity modulated radiotherapy (IMRT) and stereotactic body radiation therapy (SBRT). These imaging techniques provide predominantly anatomical information for treatment planning and target localization. Recently, studies have shown that treatment planning based on functional and molecular information about the tumor and surrounding tissue could potentially improve the effectiveness of radiation therapy. However, current on-board imaging systems are limited in their functional and molecular imaging capability. Single Photon Emission Computed Tomography (SPECT) is a candidate to achieve on-board functional and molecular imaging. Traditional SPECT systems typically take 20 minutes or more for a scan, which is too long for on-board imaging. A robotic multi-pinhole SPECT system was proposed in this dissertation to provide shorter imaging time by using a robotic arm to maneuver the multi-pinhole SPECT system around the patient in position for radiation therapy. </p><p>A 49-pinhole collimated SPECT detector and its shielding were designed and simulated in this work using the computer-aided design (CAD) software. The trajectories of robotic arm about the patient, treatment table and gantry in the radiation therapy room and several detector assemblies such as parallel holes, single pinhole and 49 pinholes collimated detector were investigated. The rail mounted system was designed to enable a full range of detector positions and orientations to various crucial treatment sites including head and torso, while avoiding collision with linear accelerator (LINAC), patient table and patient.</p><p>An alignment method was developed in this work to calibrate the on-board robotic SPECT to the LINAC coordinate frame and to the coordinate frames of other on-board imaging systems such as CBCT. This alignment method utilizes line sources and one pinhole projection of these line sources. The model consists of multiple alignment parameters which maps line sources in 3-dimensional (3D) space to their 2-dimensional (2D) projections on the SPECT detector. Computer-simulation studies and experimental evaluations were performed as a function of number of line sources, Radon transform accuracy, finite line-source width, intrinsic camera resolution, Poisson noise and acquisition geometry. In computer-simulation studies, when there was no error in determining angles (&alpha;) and offsets (&rho;) of the measured projections, the six alignment parameters (3 translational and 3 rotational) were estimated perfectly using three line sources. When angles (&alpha;) and offsets (&rho;) were provided by Radon transform, the estimation accuracy was reduced. The estimation error was associated with rounding errors of Radon transform, finite line-source width, Poisson noise, number of line sources, intrinsic camera resolution and detector acquisition geometry. The estimation accuracy was significantly improved by using 4 line sources rather than 3 and also by using thinner line-source projections (obtained by better intrinsic detector resolution). With 5 line sources, median errors were 0.2 mm for the detector translations, 0.7 mm for the detector radius of rotation, and less than 0.5° for detector rotation, tilt and twist. In experimental evaluations, average errors relative to a different, independent registration technique were about 1.8 mm for detector translations, 1.1 mm for the detector radius of rotation (ROR), 0.5° and 0.4° for detector rotation and tilt, respectively, and 1.2° for detector twist. </p><p>Simulation studies were performed to investigate the improvement of imaging sensitivity and accuracy of hot sphere localization for breast imaging of patients in prone position. A 3D XCAT phantom was simulated in the prone position with nine hot spheres of 10 mm diameter added in the left breast. A no-treatment-table case and two commercial prone breast boards, 7 and 24 cm thick, were simulated. Different pinhole focal lengths were assessed for root-mean-square-error (RMSE). The pinhole focal lengths resulting in the lowest RMSE values were 12 cm, 18 cm and 21 cm for no table, thin board, and thick board, respectively. In both no table and thin board cases, all 9 hot spheres were easily visualized above background with 4-minute scans utilizing the 49-pinhole SPECT system while seven of nine hot spheres were visible with the thick board. In comparison with parallel-hole system, our 49-pinhole system shows reduction in noise and bias under these simulation cases. These results correspond to smaller radii of rotation for no-table case and thinner prone board. Similarly, localization accuracy with the 49-pinhole system was significantly better than with the parallel-hole system for both the thin and thick prone boards. Median localization errors for the 49-pinhole system with the thin board were less than 3 mm for 5 of 9 hot spheres, and less than 6 mm for the other 4 hot spheres. Median localization errors of 49-pinhole system with the thick board were less than 4 mm for 5 of 9 hot spheres, and less than 8 mm for the other 4 hot spheres. </p><p>Besides prone breast imaging, respiratory-gated region-of-interest (ROI) imaging of lung tumor was also investigated. A simulation study was conducted on the potential of multi-pinhole, region-of-interest (ROI) SPECT to alleviate noise effects associated with respiratory-gated SPECT imaging of the thorax. Two 4D XCAT digital phantoms were constructed, with either a 10 mm or 20 mm diameter tumor added in the right lung. The maximum diaphragm motion was 2 cm (for 10 mm tumor) or 4 cm (for 20 mm tumor) in superior-inferior direction and 1.2 cm in anterior-posterior direction. Projections were simulated with a 4-minute acquisition time (40 seconds per each of 6 gates) using either the ROI SPECT system (49-pinhole) or reference single and dual conventional broad cross-section, parallel-hole collimated SPECT. The SPECT images were reconstructed using OSEM with up to 6 iterations. Images were evaluated as a function of gate by profiles, noise versus bias curves, and a numerical observer performing a forced-choice localization task. Even for the 20 mm tumor, the 49-pinhole imaging ROI was found sufficient to encompass fully usual clinical ranges of diaphragm motion. Averaged over the 6 gates, noise at iteration 6 of 49-pinhole ROI imaging (10.9 µCi/ml) was approximately comparable to noise at iteration 2 of the two dual and single parallel-hole, broad cross-section systems (12.4 µCi/ml and 13.8 µCi/ml, respectively). Corresponding biases were much lower for the 49-pinhole ROI system (3.8 µCi/ml), versus 6.2 µCi/ml and 6.5 µCi/ml for the dual and single parallel-hole systems, respectively. Median localization errors averaged over 6 gates, for the 10 mm and 20 mm tumors respectively, were 1.6 mm and 0.5 mm using the ROI imaging system and 6.6 mm and 2.3 mm using the dual parallel-hole, broad cross-section system. The results demonstrate substantially improved imaging via ROI methods. One important application may be gated imaging of patients in position for radiation therapy.</p><p>A robotic SPECT imaging system was constructed utilizing a gamma camera detector (Digirad 2020tc) and a robot (KUKA KR150-L110 robot). An imaging study was performed with a phantom (PET CT Phantom<super>TM</super>), which includes 5 spheres of 10, 13, 17, 22 and 28 mm in diameter. The phantom was placed on a flat-top couch. SPECT projections were acquired with a parallel-hole collimator and a single-pinhole collimator both without background in the phantom, and with background at 1/10th the sphere activity concentration. The imaging trajectories of parallel-hole and pinhole collimated detectors spanned 180 degrees and 228 degrees respectively. The pinhole detector viewed a 14.7 cm-diameter common volume which encompassed the 28 mm and 22 mm spheres. The common volume for parallel-hole was a 20.8-cm-diameter cylinder which encompassed all five spheres in the phantom. The maneuverability of the robotic system was tested by navigating the detector to trace the flat-top table while avoiding collision with the table and maintaining the closest possible proximity to the common volume. For image reconstruction, detector trajectories were described by radius-of-rotation and detector rotation angle &#952;. These reconstruction parameters were obtained from the robot base and tool coordinates. The robotic SPECT system was able to maneuver the parallel-hole and pinhole collimated SPECT detectors in close proximity to the phantom, minimizing impact of the flat-top couch on detector to center-of-rotation (COR) distance. In no background case, all five spheres were visible in the reconstructed parallel-hole and pinhole images. In with background case, three spheres of 17, 22 and 28 mm diameter were readily observed with the parallel-hole imaging, and the targeted spheres (22 and 28 mm diameter) were readily observed in the pinhole ROI imaging.</p><p>In conclusion, the proposed on-board robotic SPECT can be aligned to LINAC/CBCT with a single pinhole projection of the line-source phantom. Alignment parameters can be estimated using one pinhole projection of line sources. This alignment method may be important for multi-pinhole SPECT, where relative pinhole alignment may vary during rotation. For single pinhole and multi-pinhole SPECT imaging onboard radiation therapy machines, the method could provide alignment of SPECT coordinates with those of CBCT and the LINAC. In simulation studies of prone breast imaging and respiratory-gated lung imaging, the 49-pinhole detector showed better tumor contrast recovery and localization in a 4-minute scan compared to parallel-hole detector. On-board SPECT could be achieved by a robot maneuvering a SPECT detector about patients in position for radiation therapy on a flat-top couch. The robot inherent coordinate frames could be an effective means to estimate detector pose for use in SPECT image reconstruction.</p> / Dissertation

<p>Dynamic contrast-enhanced (DCE) MRI has been widely used as a quantitative imaging method for monitoring tumor response to therapy. The pharmacokinetic parameters derived from this technique have been used in more than 100 phase I trials and investigator led studies. The simultaneous challenges of increasing the temporal and spatial resolution, in a setting where the signal from the much smaller voxel is weaker, have made this MR technique difficult to implement in small-animal imaging. Existing preclinical DCE-MRI protocols acquire a limited number of slices resulting in potentially lost information in the third dimension. Furthermore, drug efficacy studies measuring the effect of an anti-angiogenic treatment, often compare the derived biomarkers on manually selected tumor regions or over the entire volume. These measurements include domains where the interpretation of the biomarkers may be unclear (such as in necrotic areas).</p><p>This dissertation describes and compares a family of four-dimensional (3D spatial + time), projection acquisition, keyhole-sampling strategies that support high spatial and temporal resolution. An interleaved 3D radial trajectory with a quasi-uniform distribution of points in k-space was used for sampling temporally resolved datasets. These volumes were reconstructed with three different k-space filters encompassing a range of possible keyhole strategies. The effect of k-space filtering on spatial and temporal resolution was studied in phantoms and in vivo. The statistical variation of the DCE-MRI measurement is analyzed by considering the fundamental sources of error in the MR signal intensity acquired with the spoiled gradient-echo (SPGR) pulse sequence. Finally, the technique was applied for measuring the extent of the opening of the blood-brain barrier in a mouse model of pediatric glioma and for identifying regions of therapeutic effect in a model of colorectal adenocarcinoma. </p><p>It is shown that 4D radial keyhole imaging does not degrade the system spatial and temporal resolution at a cost of 20-40% decrease in SNR. The time-dependent concentration of the contrast agent measured in vivo is within the theoretically predicted limits. The uncertainty in measuring the pharmacokinetic parameters with the sequences is of the same order, but always higher than, the uncertainty in measuring the pre-injection longitudinal relaxation time. The histogram of the time-to-peak provides useful knowledge about the spatial distribution of K^trans and microvascular density. Two regions with distinct kinetic parameters were identified when the TTP map from DCE-MRM was thresholded at 1000 sec. The effect of bevacizumab, as measured by a decrease in K^trans, was confined to one of these regions. DCE-MRI studies may contribute unique insights into the response of the tumor microenvironment to therapy.</p> / Dissertation

Magnetic resonance based diffusion imaging has been gaining more utility and clinical relevance over the past decade. Using conventional echo planar techniques it is possible to acquire and characterize water diffusion within the central nervous system (CNS); namely in the form of Diffusion Weighted Imaging (DWI) and Diffusion Tensor Imaging (DTI). While each modality provides valuable clinical information in terms of the presence of diffusion, DWI, and its directionality, DTI, the techniques used for analysis are limited to assuming an ideal Gaussian distribution for water displacement with no intermolecular interactions. This assumption reduces the amount of relevant information that can be interpreted in a clinical setting. By measuring the excess kurtosis, or peakedness, of the Gaussian distribution it is possible to get a better understanding of the underlying cellular structure. The objective of this work is to provide mathematical and experimental evidence that Diffusion Kurtosis Imaging (DKI) can provide additional information about the micromolecular environment of the pediatric spinal cord by more completely characterizing the probabilistic nature of random water displacement. A novel DKI imaging sequence based on a 2D spatially selective radio frequency pulse providing reduced FOV imaging with view angle tilting (VAT) was implemented, optimized on a 3Tesla MRI scanner, and tested on pediatric subjects (normal:15; patients with spinal cord injury:5). Software was developed and validated in-house for post processing of the DKI images and estimation of the tensor parameters. The results show statistically significant differences in kurtosis parameters (mean kurtosis, axial kurtosis) between normal and patients. DKI provides incremental and new information over conventional diffusion acquisitions that can be integrated into clinical protocols when coupled with higher order estimation algorithms. / Electrical and Computer Engineering

<p>The work presented in this dissertation focuses on evaluation of absolute quantification accuracy in dedicated breast SPECT-CT. The overall goal was to investigate through simulations and measurements the impact and utilization of various correction methods for scattered and attenuated photons, characterization of incomplete charge collection in Cadmium Zinc Telluride detectors as a surrogate means of improving scatter correction, and resolution recovery methods for modeling collimator blur during image reconstruction. The quantification accuracy of attenuation coefficients in CT reconstructions was evaluated in geometric phantoms, and a slice-by-slice breast segmentation algorithm was developed to separate adipose and glandular tissue. All correction and segmentation methods were then applied to a pilot study imaging parathyroid patients to determine the average uptake of Tc-99m Sestamibi in healthy breast tissue, including tissue specific uptake in adipose and glandular tissue. </p><p>Monte Carlo methods were utilized to examine the changes in incident scatter energy distribution on the SPECT detector as a function of 3D detector position about a pendant breast geometry. A simulated prone breast geometry with torso, heart, and liver was designed. An ideal detector was positioned at various azimuthal and tilted positions to mimic the capabilities of the breast SPECT subsystem. The limited near-photopeak scatter energy range in simulated spectra was linearly fit and the slope used to characterize changes in scatter distribution as a function of detector position. Results show that the detected scatter distribution changes with detector tilt, with increasing incidence of high energy scattered photons at larger detector tilts. However, reconstructions of various simulated trajectories show minimal impact on quantification (<5%) compared to a primary-only reconstruction.</p><p>Two scatter compensation methods were investigated and compared to a narrow photopeak-only windowing for quantification accuracy in large uniform regions and small, regional uptake areas: 1) a narrow ±4% photopeak energy window to minimize scatter in the photopeak window, 2) the previously calibrated dual-energy window scatter correction method, and 3) a modified dual-energy window correction method that attempts to account for the effects of incomplete charge collection in Cadmium Zinc Telluride detectors. Various cylindrical phantoms, including those with imbedded hot and cold regions, were evaluated. Results show that the Photopeak-only and DEW methods yield reasonable quantification accuracy (within 10%) for a wide range of activity concentrations and phantom configurations. The mDEW demonstrated highly accurate quantification measurements in large, uniform regions with improved uniformity compared to the DEW method. However, the mDEW method is susceptible to the calibration parameters and the activity concentration of the scanned phantom. The sensitivity of the mDEW to these factors makes it a poor choice for robust quantification applications. Thus, the DEW method using a high-performance CZT gamma camera is still a better choice for quantification purposes</p><p>Phantoms studies were performed to investigate the application of SPECT vs CT attenuation correction. Minor differences were observed between SPECT and CT maps when assuming a uniformly filled phantom with the attenuation coefficient of water, except when the SPECT attenuation map volume was significantly larger than the CT volume. Material specific attenuation coefficients reduce the corresponding measured activity concentrations compared to a water-only correction, but the results do not appear more accurate than a water-only attenuation map. Investigations on the impact of image registration show that accurate registration is necessary for absolute quantification, with errors up to 14% observed for 1.5cm shifts. </p><p>A method of modeling collimator resolution within the SPECT reconstruction algorithm was investigated for its impact on contrast and quantification accuracy. Three levels of resolution modeling, each with increasing ray-sampling, were investigated. The resolution model was applied to both cylindrical and anthropomorphic breast phantoms with hot and cold regions. Large volume quantification results (background measurements) are unaffected by the application of resolution modeling. For smaller chambers and simulated lesions, contrast generally increases with resolution modeling. Edges of lesions also appear sharper with resolution modeling. No significant differences were seen between the various levels of resolution modeling. However, Gibbs artifacts are amplified at the boundaries of high contrast regions, which can significantly affect absolute quantification measurements. Convergence with resolution modeling is also notably slower, requiring more iterations with OSEM to reach a stable mean activity concentration. Additionally, reconstructions require far more computing time with resolution modeling due to the increase in number of sampling rays. Thus while the edge enhancement and contrast improvements may benefit lesion detection, the artifacts, slower convergence, and increased reconstruction time limit the utility of resolution modeling for both absolute quantification and clinical imaging studies. </p><p>Finally, a clinical pilot study was initiated to measure the average uptake of Tc-99m Sestamibi in healthy breast tissue. Subjects were consented from those undergoing diagnostic parathyroid studies at Duke. Each subject was injected with 25mCi of Sestamibi as part of their pre-surgical parathyroid SPECT imaging studies and scanned with the dedicated breast SPECT-CT system before their diagnostic parathyroid SPECT scan. Based on phantom studies of CT reconstructed attenuation coefficient accuracy, a slice-by-slice segmentation algorithm was developed to separate breast CT data into adipose and glandular tissue. SPECT data were scatter, attenuation, and decay corrected to the time of injection. Segmented CT images were used to measure average radiotracer concentration in the whole breast, as well as adipose and glandular tissue. With 8 subjects scanned, the average measured whole breast activity concentration was found to be 0.10µCi/mL. No significant differences were seen between adipose and glandular tissue uptake. </p><p>In conclusion, the application of various characterization and correct methods for quantitative SPECT imaging were investigated. Changes in detected scatter distribution appear to have minimal impact on quantification, and characterization of low-energy tailing for a modified scatter subtraction method yields inferior overall quantification results. Comparable quantification accuracy is seen with SPECT and CT-based attenuation maps, assuming the SPECT-based volume is fairly accurate. In general, resolution recovery within OSEM yields higher contrast, but quantification accuracy appears more susceptible to measurement location. Finally, scatter, attenuation, and resolution recovery methods, along with a breast segmentation algorithm, were implemented in a clinical imaging study for quantifying Tc-99m Sestamibi uptake. While the average whole breast uptake was measured to be 0. 10µCi/mL, no significant differences were seen between adipose and glandular tissue or when implementing resolution recovery. Thus, for future clinical imaging, it's recommended that the application of the investigated correction methods should be limited to the traditional DEW method and CT-based attenuation maps for quantification studies.</p> / Dissertation

<p>The behavior of light after interacting with a biological medium reveals a wealth of information that may be used to distinguish between normal and disease states. This may be achieved by simply imaging the morphology of tissues or individual cells, and/or by more sophisticated methods that quantify specific surrogate biomarkers of disease. To this end, the work presented in this dissertation demonstrates novel tools derived from low coherence interferometry (LCI) that quantitatively measure wavelength-dependent scattering and absorption properties of biological samples, with high spectral resolution and micrometer spatial resolution, to provide insight into disease states. </p><p>The presented work first describes a dual window (DW) method, which decomposes a signal sampled in a single domain (in this case the frequency domain) to a distribution that simultaneously contains information from both the original domain and the conjugate domain (here, the temporal or spatial domain). As the name suggests, the DW method utilizes two independently adjustable windows, each with different spatial and spectral properties to overcome limitations found in other processing methods that seek to obtain the same information. A theoretical treatment is provided, and the method is validated through simulations and experiments. With this tool, the spatially dependent spectral behavior of light after interacting with a biological medium may be analyzed to extract parameters of interest, such as the scattering and absorption properties. </p><p>The DW method is employed to investigate scattering properties of samples using Fourier domain LCI (fLCI). In this method, induced temporal coherence effects provide insight into structural changes in dominant scatterers, such as cell nuclei within tissue, which can reveal the early stages of cancerous development. fLCI is demonstrated in complex, three-dimensional samples using a scattering phantom and an ex-vivo animal model. The results from the latter study show that fLCI is able to detect changes in the morphology of tissues undergoing precancerous development. </p><p>The DW method is also employed to enable a novel form of optical coherence tomography (OCT), an imaging modality that uses coherence gating to obtain micrometer-scale, cross-sectional information of tissues. The novel method, named molecular imaging true color spectroscopic OCT (METRiCS OCT), analyses the depth dependent absorption of light to ascertain quantitative information of chromophore concentration, such as hemoglobin. The molecular information is also processed to yield a true color representation of the sample, a unique capability of this approach. A number of experiments, including hemoglobin absorbing phantoms and in-vivo imaging of a chick embryo model and dorsal skinfold window chamber model, demonstrate the power of the method. </p><p>The final method presented in this dissertation, consists of a spectroscopic approach that interrogates the dispersive biochemical properties of samples to independently probe the scattering and absorption coefficients. To demonstrate this method, named non-linear phase dispersion spectroscopy (NLDS), a careful analysis of LCI signals is presented. The method is verified using measurements from samples that scatter and absorb light. Lastly, NLDS is combined with phase microscopy to achieve molecular imaging with sub-micron spatial resolution. Imaging of red blood cells (RBCs) shows that the method enables highly sensitive measurements that can quantify hemoglobin content from single RBCs.</p> / Dissertation

<p>In the United States, heart failure is a leading cause of hospitalization. The medical industry places great emphasis on diagnosing heart disease through cardiac metrics like ejection fraction. Left ventricular ejection fraction is a commonly used diagnostic indicator for heart efficiency and is measured with echocardiography through different volume calculation techniques. However, ejection fraction results can drastically vary from one examiner to another. Generally cardiologists still give ejection fraction measurements a plus or minus 10 percent error range. </p><p>A program developed at Duke called 4DViz is robust enough for users to process 3D ultrasound data. 4DViz allows examiners to determine heart chamber volumes by constructing a surface model over an imaged heart chamber with many mouse click inputs. Through 4DViz programming, a viable approach for calculating ejection fraction is attempted in this thesis. Using feature tracking, surface drawing, and voxel filling, the new approach aims to reduce examiner input and improve ejection fraction consistency. Water filled balloons were used to calibrate the algorithm's parameters. In testing, several volunteers were asked to use the 4DViz. Their results are compared to volume measurements where user input was standard. The results show promise and may remove some of the inconsistency behind ejection fraction measurements.</p> / Thesis

<p>The diagnosis and treatment of pulmonary diseases still rely on pulmonary function tests that offer archaic or insensitive biomarkers of lung structure and function. As a consequence, chronic obstructive pulmonary disease is the third leading cause of death in the US, and the hospitalization costs for asthma are on the order of $29 Billion. Pulmonary diseases have created a large and unsustainable economic burden, and hence there is still a dire need for biomarkers that can predict early changes in lung function. The work presented in this thesis looks to address this very issue, by taking advantage of the unique properties of hyperpolarized (HP) <super>129</super>Xe in conjunction with magnetic resonance imaging (MRI), to probe the fundamental function of the lung - gas-exchange. </p><p>While a bulk of the inhaled HP <super>129</super>Xe stays in the alveolar spaces, its moderate solubility in the pulmonary tissues causes a small fraction of this xenon in the alveolar spaces to diffuse into the pulmonary barrier tissue and plasma, and further into the red blood cells (RBC). Additionally, when in either of these compartments, xenon experiences a unique shift in its resonance frequency from the gas-phase (barrier - 198 ppm, RBC - 217 ppm). These unique resonances are collectively called the dissolved-phase of xenon. As the pathway taken by xenon to reach the RBCs is identical to that of oxygen, this dissolved-phase offers a non-invasive probe to study the oxygen transfer pathway, and imaging its distribution, to first order, would give us an image of gas-exchange in the lung.</p><p>Gas-exchange is controlled by ventilation, perfusion, and lastly diffusion of gases across the capillary membrane. This process of diffusion is dictated by Fick's first law of diffusion, and hence the volume of gas taken up by the capillary blood stream depends on the alveolar surface area, and the interstitial thickness. Interestingly, changes in these factors can be measured using the resonances of xenon. Changes in the alveolar surface area brought on by diseases like emphysema will increase the diffusion of xenon within the alveolus. Thus, by using diffusion-weighted imaging of the gas-phase of <super>129</super>Xe, which is the focus of chapter 3, one can extract the `apparent diffusion coefficient' (ADC) of xenon, that is sensitive to the changes in the alveolar surface area. The dissolved-phase on the other hand, while sensitive to the surface area, is also sensitive to subtle changes in the interstitial thickness. In fact, after the application of an RF pulse on the dissolved-phase, the recovery time for the xenon signal in the RBCs is significantly delayed by micron scale thickening of the interstitium. This delayed signal recovery can be used as a sensitive marker for diffusion impairment in the lung. </p><p>While direct imaging of the dissolved-phase was shown to be feasible, truly quantifying gas-exchange in the lung will require two additional technical advances - 1) As the gas-phase is the source magnetization for the dissolved-phase signal, it is imperative to acquire both the gas and dissolved-phase images in a single breath. The technical details of this achievement are discussed in chapters 4 and 5. 2) As the dissolved-phase consists of both the barrier and the RBC components, obtaining a fundamental image of gas-exchange in the lung will require creating independent images of <super>129</super>Xe in the barrier and <super>129</super>Xe in the RBCs. This goal first required creating a global metric of gas-transfer in the lung (chapter 6), which aided the implementation of the 1-point Dixon acquisition strategy to separate the components of the dissolved-phase. In conjunction with aim 1, it was finally possible to image all three resonances of <super>129</super>Xe in a single breath (chapter 7). These <super>129</super>Xe-RBC images were acquired in healthy volunteers and their efficacy was tested in subjects with idiopathic pulmonary fibrosis (IPF). These IPF subjects are known for their characteristic diffusion limitation, and in regions of fibrosis shown on their CT scans, the <super>129</super>Xe-RBC images showed gas-transfer defects. </p><p>Hyperpolarized <super>129</super>Xe MRI thus provides a non-invasive, ionizing radiation free method to probe ventilation, microstructural changes and most importantly, gas-exchange. These preliminary results indicate that xenon MRI has potential as a sensitive tool in therapeutic clinical trials to evaluate longitudinal changes in lung function.</p> / Dissertation

<p>Optical coherence tomography (OCT) is a non-invasive imaging modality that uses low coherence interferometry to generate three-dimensional datasets of a sample's structure. OCT has found tremendous clinical applications in imaging the retina and has demonstrated great utility in the diagnosis of various retinal diseases. However, such diagnoses rely upon the ability to observe abnormalities in the structure of the retina caused by pathology. By the time an ocular disease has progressed to the point of affecting the morphology of the retina, irreversible vision loss in the eye may already occur. Changes in the functionality of the tissue often precede changes to the structure. Thus, if imaging methods are developed to provide additional functional information about the behavior and response of the retinal tissue and vasculature, earlier treatment for disease may be prescribed, thus preserving vision for the patient. </p><p>Within the last decade, significant technological advances in OCT systems have enabled high-speed and high sensitivity image acquisition using either spectral domain OCT (SDOCT) or swept-source OCT (SSOCT) configurations. Such systems use Fourier processing to extract structural information of a sample from interferometric principles. But such systems also have access to the optical phase information, which allows for functional analysis of sample dynamics. This dissertation details the development and application of methods using both intensity and phase information as a tool for studying interesting biological phenomena. The goal of this work is an extension of techniques to image the vasculature in the retina and enhance the clinical utility of OCT.</p><p>I first outline basic theory necessary for understanding the principles of OCT. I then describe OCT phase imaging in cellular applications as a demonstration of the ability of OCT to provide functional information on biological dynamics. Phase imaging methods suffer from an artifact known as phase wrapping, and I have developed a software technique to overcome this problem in OCT, thus extending its usefulness in providing quantitative information. I characterize the limitations in measuring moving scatterers with Doppler OCT in both SDOCT and SSOCT system. I also show the ability to image the vasculature in the retina using variance imaging with a high-speed retinal imaging system and software based methods to correct for patient motion and create a widefield mosaic in an automated manner. Finally, future directions for this work are discussed.</p> / Dissertation

<p>Onboard cone beam computed tomography (CBCT) has become a widespread means of three-dimensional target localization for radiation therapy; however, it is susceptible to metal artifacts and beam-hardening artifacts that can hinder visualization of low contrast anatomy. Dual-CBCT provides easy access to techniques that may reduces such artifacts. Additionally, dual-CBCT can decrease imaging time and provide simultaneous orthogonal projections which may also be useful for fast target localization. However, dual-CBCT will suffer from large increases in scattered radiation due to the addition of the second source.</p><p>An experimental bench top dual CBCT system was constructed so that each imaging chain in the dual CBCT system mimics the geometry of gantry-mounted CBCT systems commonly used in the radiation therapy room. The two systems share a common axis of rotation and are mounted orthogonally. Custom control software was developed to ensure reproducible exposure and rotation timings. This software allows the implementation of the acquisition sequences required for the cross scatter avoidance and correction strategies studied.</p><p>Utilizing the experimental dual CBCT system cross scatter was characterized from 70-145 kVp in projections and reconstructed images using this system and three cylindrical phantoms (15cm, 20cm, and 30cm) with a common Catphan core. A novel strategy for avoiding cross-scatter in dual-CBCT was developed that utilized interleaved data acquisition on each imaging chain. Contrast and contrast-to-noise-ratio were measured in reconstructions to evaluate the effectiveness of this strategy to avoid the effects of cross scatter.</p><p>A novel correction strategy for cross scatter was developed wherein the cross scatter was regularly sampled during the course of data acquisition and these samples were used as the basis for low- and high- frequency corrections for the cross-scatter in projections. The cross scatter sampling interval was determined for an anthropomorphic phantom at three different sites relevant to radiation therapy by estimating the angular Nyquist frequency. The low frequency portion of the cross scatter distribution is interpolated between samples to provide an estimate of the cross scatter distribution at every projection angle and was then subtracted from the projections.</p><p>The high-frequency portion of the correction was applied after the low-frequency correction was applied. The novel high-frequency correction utilizes the fact that a direct estimate of the high-frequency components was obtained in the cross scatter samples. The high-frequency components of the measured cross scatter were subtracted from the projections in the Fourier domain, a process referred to as spectral subtraction. Each projection is corrected using the cross scatter sample taken at the closest projection angle. In order to apply this correction in the Fourier domain the high-frequency component of the cross scatter must be approximately stationary. To improve the stationarity of the high-frequency cross scatter component a novel two-dimensional, overlapping window was developed. The spectral subtraction was then applied in each window and the results added to form the final image.</p><p>The effectiveness of the correction techniques were evaluated by measuring the contrast and contrast-to-noise-ratio in an image quality phantom. Additionally, the effect of the high-frequency correction on resolution was measured using a line pair phantom.</p><p>Cross scatter in dual CBCT was shown for large phantoms to be much higher than forward scatter which has long been known to be one of the largest degrading factors of image quality in CBCT. This results in large losses of contrast and CNR in reconstructed images. The interleaving strategy for avoiding cross scatter during projection acquisition showed similar performance to cross scatter free acquisitions, however, does not acquire projections at the maximum possible rate. For those applications in which maximizing the acquisition rate of projections is important, the low- and high-frequency corrections effectively mitigated the effects of cross scatter in the dual CBCT system.</p> / Dissertation

<p>There is a critical need to target tumor hypoxia as patients with hypoxic tumors have worse prognosis due to aggressive phenotypes and resistance to radiotherapy and chemotherapy. The overall goal of this work is to improve response to conventional cancer therapies by targeting tumor hypoxia. This has been carried out and evaluated through the use of polymersome-encapsulated myoglobin (PEMs) with the hypothesis that O2-releasing PEMs will increase tumor oxygenation, and thereby improve response to radiotherapy. Mb was chosen as an O2 carrying protein to deliver to tumors because it has a strong association to O2, providing a mechanism to deliver O2 only within the hypoxic regions of the tumor. Mb was loaded within nanoscale polymeric vesicles that were expected to accumulate within solid tumors due to the enhanced permeability and retention (EPR) effect. This hypothesis has been tested through the following aims:</p><p>1. Develop NIR imaging techniques for studying the biodistribution and pharmacokinetics of polymersomes</p><p>2. Establish the effects of Mb-containing polymersomes on tumor physiology</p><p>3. Modify tumor growth through delivery of Mb polymersomes in combination with a cytotoxic therapy specific to aerobic tumors</p><p>These aims have been evaluated through numerous in vivo studies. First, polymersomes of various polymer formulations and diameters ranging from 110-550 nm were prepared with a near-infrared (NIR) -emissive fluorophore. Using live animal fluorescence imaging, I was able to study the biodistribution of the polymersomes following i.v. administration, demonstrating significant polymersome accumulation in orthotopic 4T1 mammary carcinomas. In addition, a novel method for measuring pharmacokinetics was developed, using serial small volume blood draws from individual mice. The plasma fluorescence in microcapillary tubes was used to quantify polymersome concentrations, demonstrating long circulation half-lives that varied from 6-23 h. Toxicity of various polymersome formulations were also studied in vitro and in vivo, revealing negligible toxicities.</p><p>For the second aim, PEMs were administered i.v. in tumor-bearing mice. Unexpectedly, we observed a dramatic gross tumor effect within hours of treatment in both orthotopic 4T1 tumors and flank Renca renal cell carcinomas. Histological analysis revealed endothelial cell apoptosis as early as 1 h following treatment, with scattered tumor cell death throughout the tumor by 4 h. Hematoxylin and eosin staining showed significant necrosis 24 h following PEM treatment. Vascular effects and polymersome distribution were studied in 4T1 window chamber tumors. Following i.v. treatment with PEMs, intravital microscopy was used to image polymersome fluorescence, brightfield transmission was imaged for vessel morphology and blood flow, and a tunable filter was used for determining hemoglobin (Hb) oxygen saturation. Tumor hemorrhaging was observed within hours of PEM treatment, which was not seen with empty polymersomes. This was consistent with the gross tumor effects observed initially. Hb saturation decreased in both the PEM and empty polymersome groups, but not in saline-treated mice. While we expected to observe an increase in tumor oxygenation by using Mb as an oxygen carrier, we actually observed hemorrhage, decreased oxygenation, and central tumor necrosis. In vitro studies using human endothelial cells demonstrated dramatic changes in cell morphology and increased permeability due to Mb and PEM treatments, which appear to be enhanced in an oxidative environment. These in vitro and in vivo observations are similar to what is seen with tumor vascular disrupting agents.</p><p>For the third aim, I combined radiotherapy (RT) and PEM treatment with a new hypothesis. I originally expected the PEMs to increase tumor oxygenation, thus making the tumor more susceptible to RT. However, considering the results from the second aim, this hypothesis was modified: the PEMs would result in necrosis of the tumor core, while RT would target the more oxygenated rim of the tumor, thus leading to improved tumor growth delay compared with PEM or RT alone. This hypothesis was tested in both orthotopic, syngeneic 4T1 tumors as well as flank FaDu xenografts. 4T1 tumor cells were surgically implanted within the dorsal mammary fat pad of mice and grown until ~200 mm3. A CT microirradiator with a square collimator was used in order to locate and specifically irradiate the tumor. Within 1 h following RT, the PEMs were administered i.v.. Mice receiving PEMs with no RT showed a significant decrease in tumor growth compared with saline-treated mice (p = 0.0001 for time to 3x original tumor volume). In addition, the combination of RT plus PEMs reduced tumor growth compared with RT alone (p = 0.0144 for time to 3x original tumor volume). However, this effect was not seen with FaDu tumors. This may have been due to excessive radiation dose or other compounding factors: the timing between RT and PEM treatment was not optimized, and the number of mice per group was small (3-4). </p><p>Thus, the conclusions for each aim are as follows:</p><p>1. Develop NIR imaging techniques for studying the biodistribution and pharmacokinetics of polymersomes</p><p>NIR imaging techniques were optimized for studying polymersomes, demonstrating long plasma circulation times and accumulation within tumors.</p><p>2. Establish the effects of Mb-containing polymersomes on tumor physiology</p><p>While the hypothesis was that PEMs would accumulate within hypoxic tumors and subsequently increase O2 tension, we observed a rapid decrease in tumor oxygenation followed by a dramatic hemorrhagic effect of Mb polymersomes, which appear to be due to both endothelial cell apoptosis and morphological changes, resulting in central tumor necrosis.</p><p>3. Modify tumor growth through delivery of Mb polymersomes in combination with a cytotoxic therapy specific to aerobic tumors</p><p>Combination therapy of PEMs with RT results in enhanced tumor growth delay in aggressive 4T1 mammary carcinomas compared with RT or PEMs alone.</p><p>These studies have led to a proposed mechanism for the PEM anti-tumor effect in combination with RT. Prior to PEM administration, RT is administered, resulting in tumor cell kill of the well-oxygenated tumor periphery. Mb polymersomes are then injected i.v. and begin to accumulate within tumors due to the EPR effect. As shown in Aim 1, this accumulation occurs over a short time scale. Within 30 min of PEM treatment, the Mb is believed to act on tumor vessels, resulting in morphological changes and apoptosis of endothelial cells. These effects are expected to increase permeability of the vessels and expose the basement membrane, which leads to clotting and decreased blood flow. Both decreased perfusion and increased permeability are believed to have a catastrophic effect on interior tumor vessels. Hemorrhage results as the endothelial cells die, resulting in tumor core necrosis. Therefore, the result is tumor cell kill at the periphery due to RT and central tumor necrosis due to PEM treatment.</p><p>PEMs have potential in cancer therapy as a new class of VDAs. While the mechanism requires further investigation, this work has demonstrated that PEM treatment results in tumor vessel destruction and central necrosis. PEMs accumulate within tumors, thus minimizing the systemic toxicity of treatment commonly seen with VDAs. By combining PEMs with a therapy that kills the better perfused tumor periphery, PEMs show promise in improving tumor response. Future mechanistic studies will be needed in order to maximize vessel damage and optimize combination dosing schedules to improve outcome.</p> / Dissertation