Multiple gradient echo propeller (MGREP) MRI: Technical development and potential applications

The PROPELLER (Periodically Rotated Overlapping Parallel Lines with Enhanced Reconstruction) magnetic resonance imaging (MRI) technique has inherent advantages over other fast imaging methods, including robust motion correction, reduced image distortion, and resistance to off-resonance effects. These features make PROPELLER highly desirable for T2*-sensitive imaging, high-resolution diffusion imaging, and many other applications. However, PROPELLER has been predominantly implemented as a fast spin-echo (FSE) technique, which is insensitive to T2* contrast, and requires time-inefficient signal averaging to achieve adequate signal-to-noise ratio (SNR) for many applications. These issues presently constrain the potential clinical utility of FSE-based PROPELLER. ^ In this research, our aim was to extend and enhance the potential applications of PROPELLER MRI by developing a novel multiple gradient echo PROPELLER (MGREP) technique that can overcome the aforementioned limitations. The MGREP pulse sequence was designed to acquire multiple gradient-echo images simultaneously, without any increase in total scan time or RF energy deposition relative to FSE-based PROPELLER. A new parameter was also introduced for direct user-control over gradient echo spacing, to allow variable sensitivity to T2* contrast. In parallel to pulse sequence development, an improved algorithm for motion correction was also developed and evaluated against the established method through extensive simulations. The potential advantages of MGREP over FSE-based PROPELLER were illustrated via three specific applications: (1) quantitative T2* measurement, (2) time-efficient signal averaging, and (3) high-resolution diffusion imaging. Relative to the FSE-PROPELLER method, the MGREP sequence was found to yield quantitative T2* values, increase SNR by ∼40% without any increase in acquisition time or RF energy deposition, and noticeably improve image quality in high-resolution diffusion maps. In addition, the new motion algorithm was found to improve the performance considerably in motion-artifact reduction. ^ Overall, this work demonstrated a number of enhancements and extensions to existing PROPELLER techniques. The new technical capabilities of PROPELLER imaging, developed in this thesis research, are expected to serve as the foundation for further expanding the scope of PROPELLER applications. ^

Verification of intensity modulated stereotactic radiotherapy using Monte Carlo calculations and EPID dosimetry

The purpose of this work was to develop a comprehensive IMSRT QA procedure that examined, using EPID dosimetry and Monte Carlo (MC) calculations, each step in the treatment planning and delivery process. These steps included verification of the field shaping, treatment planning system (RTPS) dose calculations, and patient dose delivery. Verification of each step in the treatment process is assumed to result in correct dose delivery to the patient. ^ The accelerator MC model was verified against commissioning data for field sizes from 0.8 × 0.8 cm 2 to 10 × 10 cm 2. Depth doses were within 2% local percent difference (LPD) in low gradient regions and 1 mm distance to agreement (DTA) in high gradient regions. Lateral profiles were within 2% LPD in low gradient regions and 1 mm DTA in high gradient regions. Calculated output factors were within 1% of measurement for field sizes ≥1 × 1 cm2. ^ The measured and calculated pretreatment EPID dose patterns were compared using criteria of 5% LPD, 1 mm DTA, or 2% of central axis pixel value with ≥95% of compared points required to pass for successful verification. Pretreatment field verification resulted in 97% percent of the points passing. ^ The RTPS and Monte Carlo phantom dose calculations were compared using 5% LPD, 2 mm DTA, or 2% of the maximum dose with ≥95% of compared points required passing for successful verification. RTPS calculation verification resulted in 97% percent of the points passing. ^ The measured and calculated EPID exit dose patterns were compared using criteria of 5% LPD, 1 mm DTA, or 2% of central axis pixel value with ≥95% of compared points required to pass for successful verification. Exit dose verification resulted in 97% percent of the points passing. ^ Each of the processes above verified an individual step in the treatment planning and delivery process. The combination of these verification steps ensures accurate treatment delivery to the patient. This work shows that Monte Carlo calculations and EPID dosimetry can be used to quantitatively verify IMSRT treatments resulting in improved patient care and, potentially, improved clinical outcome. ^

Assessment of cone beam computed tomography techniques for imaging lung damage in mice in vivo

Lung damage is a common side effect of chemotherapeutic drugs such as bleomycin. This study used a bleomycin mouse model which simulates the lung damage observed in humans. Noninvasive, in vivo cone-beam computed tomography (CBCT) was used to visualize and quantify fibrotic and inflammatory damage over the entire lung volume of mice. Bleomycin was used to induce pulmonary damage in vivo and the results from two CBCT systems, a micro-CT and flat panel CT (fpCT), were compared to histologic measurements, the standard method of murine lung damage quantification. Twenty C57BL/6 mice were given either 3 U/kg of bleomycin or saline intratracheally. The mice were scanned at baseline, before the administration of bleomycin, and then 10, 14, and 21 days afterward. At each time point, a subset of mice was sacrificed for histologic analysis. The resulting CT images were used to assess lung volume. Percent lung damage (PLD) was calculated for each mouse on both the fpCT (PLDfpcT) and the micro-CT (PLDμCT). Histologic PLD (PLDH) was calculated for each histologic section at each time point (day 10, n = 4; day 14, n = 4; day 21, n = 5; control group, n = 5). A linear regression was applied to the PLDfpCT vs. PLDH, PLDμCT vs. PLDH and PLDfpCT vs. PLDμCT distributions. This study did not demonstrate strong correlations between PLDCT and PLDH. The coefficient of determination, R, was 0.68 for PLDμCT vs. PLDH and 0.75 for the PLD fpCT vs. PLDH. The experimental issues identified from this study were: (1) inconsistent inflation of the lungs from scan to scan, (2) variable distribution of damage (one histologic section not representative of overall lung damage), (3) control mice not scanned with each group of bleomycin mice, (4) two CT systems caused long anesthesia time for the mice, and (5) respiratory gating did not hold the volume of lung constant throughout the scan. Addressing these issues might allow for further improvement of the correlation between PLDCT and PLDH. ^

Development of automated image analysis tools for verification of radiotherapy field accuracy with an electronic portal imaging device

The successful management of cancer with radiation relies on the accurate deposition of a prescribed dose to a prescribed anatomical volume within the patient. Treatment set-up errors are inevitable because the alignment of field shaping devices with the patient must be repeated daily up to eighty times during the course of a fractionated radiotherapy treatment. With the invention of electronic portal imaging devices (EPIDs), patient's portal images can be visualized daily in real-time after only a small fraction of the radiation dose has been delivered to each treatment field. However, the accuracy of human visual evaluation of low-contrast portal images has been found to be inadequate. The goal of this research is to develop automated image analysis tools to detect both treatment field shape errors and patient anatomy placement errors with an EPID. A moments method has been developed to align treatment field images to compensate for lack of repositioning precision of the image detector. A figure of merit has also been established to verify the shape and rotation of the treatment fields. Following proper alignment of treatment field boundaries, a cross-correlation method has been developed to detect shifts of the patient's anatomy relative to the treatment field boundary. Phantom studies showed that the moments method aligned the radiation fields to within 0.5mm of translation and 0.5$\sp\circ$ of rotation and that the cross-correlation method aligned anatomical structures inside the radiation field to within 1 mm of translation and 1$\sp\circ$ of rotation. A new procedure of generating and using digitally reconstructed radiographs (DRRs) at megavoltage energies as reference images was also investigated. The procedure allowed a direct comparison between a designed treatment portal and the actual patient setup positions detected by an EPID. Phantom studies confirmed the feasibility of the methodology. Both the moments method and the cross-correlation technique were implemented within an experimental radiotherapy picture archival and communication system (RT-PACS) and were used clinically to evaluate the setup variability of two groups of cancer patients treated with and without an alpha-cradle immobilization aid. The tools developed in this project have proven to be very effective and have played an important role in detecting patient alignment errors and field-shape errors in treatment fields formed by a multileaf collimator (MLC). ^

Pencil -beam redefinition algorithm dose calculations for electron therapy treatment planning

The electron pencil-beam redefinition algorithm (PBRA) of Shiu and Hogstrom has been developed for use in radiotherapy treatment planning (RTP). Earlier studies of Boyd and Hogstrom showed that the PBRA lacked an adequate incident beam model, that PBRA might require improved electron physics, and that no data existed which allowed adequate assessment of the PBRA-calculated dose accuracy in a heterogeneous medium such as one presented by patient anatomy. The hypothesis of this research was that by addressing the above issues the PBRA-calculated dose would be accurate to within 4% or 2 mm in regions of high dose gradients. A secondary electron source was added to the PBRA to account for collimation-scattered electrons in the incident beam. Parameters of the dual-source model were determined from a minimal data set to allow ease of beam commissioning. Comparisons with measured data showed 3% or better dose accuracy in water within the field for cases where 4% accuracy was not previously achievable. A measured data set was developed that allowed an evaluation of PBRA in regions distal to localized heterogeneities. Geometries in the data set included irregular surfaces and high- and low-density internal heterogeneities. The data was estimated to have 1% precision and 2% agreement with accurate, benchmarked Monte Carlo (MC) code. PBRA electron transport was enhanced by modeling local pencil beam divergence. This required fundamental changes to the mathematics of electron transport (divPBRA). Evaluation of divPBRA with the measured data set showed marginal improvement in dose accuracy when compared to PBRA; however, 4% or 2mm accuracy was not achieved by either PBRA version for all data points. Finally, PBRA was evaluated clinically by comparing PBRA- and MC-calculated dose distributions using site-specific patient RTP data. Results show PBRA did not agree with MC to within 4% or 2mm in a small fraction (<3%) of the irradiated volume. Although the hypothesis of the research was shown to be false, the minor dose inaccuracies should have little or no impact on RTP decisions or patient outcome. Therefore, given ease of beam commissioning, documentation of accuracy, and calculational speed, the PBRA should be considered a practical tool for clinical use. ^

Studies on regulation of skeletal muscle differentiation by growth factors

Skeletal muscle differentiation involves sequential events in which proliferating undifferentiated myoblasts withdraw from the cell cycle and fuse to form multinucleated myotubes. The process of fusion is accompanied by the disappearance of proteins associated with cell proliferation and the coordinate induction of a battery of muscle-specific gene products, which includes the muscle isoenzyme of creatine kinase, nicotinic acetylcholine receptor, and contractile proteins such as alpha-actin. The molecular events associated with myogenesis are particularly amenable to experimental analysis because the events which occur in vivo can be recapitulated in vitro using established muscle cell lines. Initiation of myogenic differentiation in vitro can be achieved by removing serum from the culture medium. Myogenesis, therefore, can be considered to be regulated through a repression-type of mechanism by components in serum. The objectives of this project were to identify the components involved in regulation of myogenesis and approach the mechanism(s) whereby these components achieve their regulatory function. Initially, the effects of a series of polypeptide growth factors on myogenesis were examined. Among them TGF$\beta$ and FGF were found to be potent inhibitors of myogenic differentiation which did not affect cell proliferation. The inhibitory effects of these growth factors on differentiation requires their persistent presence in the culture medium. After myoblasts have undergone fusion, they become refractory to the inhibitory effects of TGF$\beta$, FGF, and serum. When fusion is inhibited by the presence of EGTA, a Ca$\sp{2+}$ chelator, muscle-specific genes are expressed reversibly upon removal of inhibitory growth factors. Subsequent exposure of biochemically differentiated cells to serum or TGF$\beta$ leads to down-regulation of muscle-specific genes. Stimulation with serum also leads to reentry of myocytes into the cell cycle, whereas fused myotubes are irreversibly and terminally differentiated. Measurement of levels of TGF$\beta$ receptors reveals that under non-fusing conditions, TGF$\beta$ receptor levels in biochemically differentiated myocytes remained as high as in undifferentiated myoblasts, while during terminal differentiation, TGF$\beta$ receptors decreased at least five-fold. Thus, down-regulation of TGF$\beta$ receptors is coupled to irreversible differentiation, but not reversible differentiation in the absence of fusion. The possible involvement of second messenger systems, such as cAMP and protein kinase C, in the pathway(s) by which TGF$\beta$, FGF, or serum factors transduce their signals from the cell surface to the nucleus was also examined. The results showed that myogenic differentiation is subject to negative regulation through cAMP elevation-dependent and cAMP elevation-independent pathways and that serum mitogens, TGF$\beta$ and FGF inhibit differentiation through a mechanism independent of cAMP-elevation or protein kinase C activation. ^

Studies of the mechanism of daunorubicin carcinogenesis: Modulation by vitamin E

Daunorubicin (DNR) is an anthracycline antibiotic used as a cancer chemotherapeutic agent. However, it causes mammary adenocarcinomas in female Sprague-Dawley (SD) rats. Vitamin E (E) has been found to reduce DNR carcinogenicity. I investigated the mechanism of DNR carcinogenicity and its interaction with E in SD rats by studying DNR-DNA adduct formation and the influence of E status on DNR clearance and free radical producing and detoxifying enzymes.^ The hypothesis was that DNR exerts its tumorigenic effect via free radicals generated during redox cycling and production of reactive intermediates capable of forming DNA adducts. E was postulated to act as a protective agent through a combination of its antioxidant property, modulation of drug clearance and levels of free radical producing and detoxifying enzymes.^ DNA adduct formation was measured by the nuclease P1 $\sp{32}$P-post labeling assay. In vitro, DNR was activated by rat liver microsomes and either NADPH or cumene hydrogen peroxide (CuOOH). Rat liver DNA incubated with this mixture formed two adducts when the cofactor was NADPH and three adducts when CuOOH was used. In vivo, SD rats were treated with i.v. doses of DNR. No detectable DNR-DNA adducts were formed in liver or mammary DNA in vivo, although there was an intensification of endogenous DNA adducts.^ Groups, 1, 2, 3 and 4 of weanling female SD rats were fed 0, 100, 1,000 and 10,000 mg $\alpha$-tocopheryl acetate/kg diet respectively. A comparison of Groups 1 and 4 showed no effect of E status on clearance of 10 mg tritiated DNR/kg body weight over 72 hours. However, liver cleared DNR at a faster rate than mammary epithelial cells (MEC).^ Xanthine oxidase, which catalyzes DNR redox cycling, was significantly decreased in liver and MEC of rats in group 4 compared to groups 1, 2, and 3. Detoxifying enzymes were not dramatically affected by E supplementation. Quinone reductase in MEC was significantly increased in group 4 compared to other groups. Overall, the liver had higher levels of free radical detoxifying enzymes compared to MEC.^ These data support a role of free radicals in DNR carcinogenicity because (1) endogenous DNA adducts formed due to free radical insult are further intensified by DNR treatment in vivo, (2) MEC, the specific target of DNR carcinogenicity, cannot rapidly clear DNR and have a lower free radical detoxifying capability than liver, (3) E supplementation caused lowering of free radical generating potential via xanthine oxidase, and increased DNR detoxification due to elevation of quinone reductase in MEC. ^

A genetically engineered, adherent smooth muscle cell lining for left ventricular assist devices

Due to the clinical success of left ventricular assist devices (LVADs) used for short term "bridge to transplant" and the limited availability of donor organs, heart assist devices are being considered for long term implantation as an alternative to heart transplantation. In an effort to improve biocompatibility, a nonthrombogenic cellular lining was developed from genetically engineered smooth muscle cells (GE-SMC) for the Thermocardiosystems Heartmate$\sp{\rm TM}$ LVAD. SMCs have been transduced with the genes for endothelial nitric oxide synthase (NOS III) and GTP cyclohydrolase (GTPCH) with subsequent stable expression of the NOS III protein via an Epstein Barr based DNA expression vector. Transduced SMCs produce nitric oxide at concentrations that reduce platelet deposition and smooth muscle cell proliferation when tested in vitro. In addition, the adhesive capabilities of GE-SMC linings were also examined, and optimized in physical environments mimicking typical in vivo LVAD operation. Preliminary investigations examining cell adhesion during constant shear stress exposure demonstrated an acute phase of cell loss corresponding to cytoskeletal F-actin rearrangement. Subsequently, an in vitro circulatory loop was designed to expose cell lined LVADs to in vivo operating conditions. Cumulative cell loss from cell lined LVADs was less than 10% after 24 hours of flow. Using a protocol for "preconditioning" the cell lining within the mock circulatory loop, the first implantation of an LVAD containing a genetically engineered SMC lining was successfully implemented in a bovine model. Results from this 24 hour study indicate that the flow-conditioned cellular lining remained intact with no evidence of thromboembolization and only minimal changes in coagulation studies. ^

Diffusion tensor imaging of rat spinal cord in vivo

Magnetic resonance imaging, with its exquisite soft tissue contrast, is an ideal modality for investigating spinal cord pathology. While conventional MRI techniques are very sensitive for spinal cord pathology, their specificity is somewhat limited. Diffusion MRI is an advanced technique which is a very sensitive and specific indicator of the integrity of white matter tracts. Diffusion imaging has been shown to detect early ischemic changes in white matter, while conventional imaging demonstrates no change. By acquiring the complete apparent diffusion tensor (ADT), tissue diffusion properties can be expressed in terms of quantitative and rotationally invariant parameters. ^ Systematic study of SCI in vivo requires controlled animal models such as the popular rat model. To date, studies of spinal cord using ADT imaging have been performed exclusively in fixed, excised spinal cords, introducing inevitable artifacts and losing the benefits of MRI's noninvasive nature. In vivo imaging reflects the actual in vivo tissue properties, and allows each animal to be imaged at multiple time points, greatly reducing the number of animals required to achieve statistical significance. Because the spinal cord is very small, the available signal-to-noise ratio (SNR) is very low. Prior spin-echo based ADT studies of rat spinal cord have relied on high magnetic field strengths and long imaging times—on the order of 10 hours—for adequate SNR. Such long imaging times are incompatible with in vivo imaging, and are not relevant for imaging the early phases following SCI. Echo planar imaging (EPI) is one of the fastest imaging methods, and is popular for diffusion imaging. However, EPI further lowers the image SNR, and is very sensitive to small imperfections in the magnetic field, such as those introduced by the bony spine. Additionally, The small field-of-view (FOV) needed for spinal cord imaging requires large imaging gradients which generate EPI artifacts. The addition of diffusion gradients introduces yet further artifacts. ^ This work develops a method for rapid EPI-based in vivo diffusion imaging of rat spinal cord. The method involves improving the SNR using an implantable coil; reducing magnetic field inhomogeneities by means of an autoshim, and correcting EPI artifacts by post-processing. New EPI artifacts due to diffusion gradients described, and post-processing correction techniques are developed. ^ These techniques were used to obtain rotationally invariant diffusion parameters from 9 animals in vivo, and were validated using the gold-standard, but slow, spinecho based diffusion sequence. These are the first reported measurements of the ADT in spinal cord in vivo . ^ Many of the techniques described are equally applicable toward imaging of human spinal cord. We anticipate that these techniques will aid in evaluating and optimizing potential therapies, and will lead to improved patient care. ^

STUDIES ON ENERGY METABOLISM IN THE HIPPOCAMPUS AFTER ISCHEMIA

The gerbil model of ischemia was used to determine the effect of carotid occlusion on energy metabolites in cellular layers of discrete regions of the hippocampus and dentate gyrus. Levels of glucose, glycogen, ATP and phosphocreatine (PCr) were unchanged after 1 minute of ischemia. However, 3 minutes of ischemia produced a dramatic decrease in net levels of all metabolites. No additional decrease was observed after 15 minutes of ischemia. Re-establishment of the blood flow for 5 minutes after a 15 minute ischemic episode returned all metabolites to pre-ischemia levels. Concentrations of glucose and glycogen were elevated in sham-operated animals as a function of the pentobarbital anesthetic employed. In other studies, elevated GABA levels (produced by inhibiting GABA-transaminase with (gamma)-vinyl-GABA (GVG)) were found to decrease the rate of utilization of the high-energy phosphate metabolites ATP and PCr in the mouse cortex. In addition, glucose and glycogen levels were increased. Thus, tonic inhibition by GABA produced decreased cellular activity. Additional experiments demonstrated the attenuation of ischemia-induced metabolite depletion in cellular layers of regions of the hippocampus, dentate gyrus and cortex after GVG administration. Under ether, 1 minute of bilateral carotid occlusion produced a dramatic decrease in metabolite levels. After GVG treatment, the decrease was blocked completely for glucose, glycogen and ATP, and partially for PCr. Therefore, GABA-transaminase inhibition produced increased levels of GABA which subsequently decreased cellular activity. The protection against ischemia may have been due to (a)decreased metabolic rate; the available energy stores were utilized at a slower rate, and (b)increased levels of energy substrates; additional supplies available to maintain viability. These data suggest that the functional state of neural tissue can determine the response to metabolic stress. ^

Enhanced in vivo angiogenesis within a model tissue engineered construct using biodegradable microspheres containing encapsulated VEGF

Introduction. Tissue engineering techniques offer a potential means to develop a tissue engineered construct (TEC) for the treatment of tissue and organ deficiencies. However, a lack of adequate vascularization is a limiting factor in the development of most viable engineered tissues. Vascular endothelial growth factor (VEGF) could aid in the development of a viable vascular network within TECs. The long-term goals of this research are to develop clinically relevant, appropriately vascularized TECs for use in humans. This project tested the hypothesis that the delivery of VEGF via controlled release from biodegradable microspheres would increase the vascular density and rate of angiogenesis within a model TEC. ^ Materials and methods. Biodegradable VEGF-encapsulated microspheres were manufactured using a novel method entitled the Solid Encapsulation/Single Emulsion/Solvent Extraction technique. Using a PLGA/PEG polymer blend, microspheres were manufactured and characterized in vitro. A model TEC using fibrin was designed for in vivo tissue engineering experimentation. At the appropriate timepoint, the TECs were explanted, and stained and quantified for CD31 using a novel semi-automated thresholding technique. ^ Results. In vitro results show the microspheres could be manufactured, stored, degrade, and release biologically active VEGF. The in vivo investigations revealed that skeletal muscle was the optimal implantation site as compared to dermis. In addition, the TECs containing fibrin with VEGF demonstrated significantly more angiogenesis than the controls. The TECs containing VEGF microspheres displayed a significant increase in vascular density by day 10. Furthermore, TECs containing VEGF microspheres had a significantly increased relative rate of angiogenesis from implantation day 5 to day 10. ^ Conclusions. A novel technique for producing microspheres loaded with biologically active proteins was developed. A defined concentration of microspheres can deliver a quantifiable level of VEGF with known release kinetics. A novel model TEC for in vivo tissue engineering investigations was developed. VEGF and VEGF microspheres stimulate angiogenesis within the model TEC. This investigation determined that biodegradable rhVEGF 165-encapsulated microspheres increased the vascular density and relative rate of angiogenesis within a model TEC. Future applications could include the incorporation of microvascular fragments into the model TEC and the incorporation of specific tissues, such as fat or bone. ^

4HPR-X radiation mechanisms of interactions in non-small cell lung cancer (NSCLC) cells in vitro

Lung cancer is the leading cause of cancer death. However, poor survival using conventional therapies fuel the search for more rational interventions. The objective of this study was to design and implement a 4HPR-radiation interaction model in NSCLC, employing a traditional clinical modality (radiation), a relatively new, therapeutically unexplored agent (4HPR) and rationally combining them based on molecular mechanistic findings pertaining to their interactions. To test the hypothesis that 4HPR sensitizes cells to radiation-induced cell death via G2+M accumulation, we designed a working model consisting of H522 adenocarcinoma cells (p53, K-ras mutated) derived from an NSCLC patient; 4HPR at concentrations up to 10 μM; and X radiation up to 6 Gy generated by a patient-dedicated Phillips RT-250 X ray unit at 250 KV, 15 mA, 1.85 Gy/min. We found that 4HPR produced time- and dose-dependent morphological changes, growth inhibition, and DNA damage-inducing enhancement of reactive oxygen species. A transient G2+M accumulation of cells maximal at 24 h of continuous 4HPR exposure was used for irradiation time scheduling. Our data demonstrated enhanced cell death (both apoptotic and necrotic) in irradiated cells pre-treated with 4HPR versus those with either stressor alone. 4HPR's effect of increased NSCLC cells' radioresponse was confirmed by clonogenic assay. To explore these practical findings from a molecular mechanistic perspective, we further investigated and showed that levels of cyclin B1 and p34cdc2 kinase—both components of the mitosis promoting factor (MPF) regulating the G2/M transition—did not change following 4HPR treatment. Likewise, cdc25C phosphatase was not altered. However, enhanced p34cdc2 phosphorylation on its Thr14Tyr15 residues—indicative of its inactivation and increased expression of MPF negative regulators chk1 and wee1 kinases—were supportive of explaining 4HPR-treated cells' accumulation. Hence, p34cdc2 phosphorylation, chk1, and wee1 warrant further evaluation as potential molecular targets for 4HPR-X radiation combination. In summary, we (1) demonstrated that 4HPR not only induces cell death by itself, but also increases NSCLC cells' subsequent radioresponse, indicative of potential clinical applicability, and (2) for the first time, shed light on deciphering 4HPR-X radiation molecular mechanisms of interaction, including the finding of 4HPR's role as a p34cdc2 inactivator via Thr14Tyr15 phosphorylation. ^

Correction of aberrant FGFR1 RNA splicing through blocking intronic regulatory elements

Analysis of the human genome has revealed that more than 74% of human genes undergo alternative RNA splicing. Aberrations in alternative RNA splicing have been associated with several human disorders, including cancer. ^ We studied the aberrant expression of alternative RNA splicing isoforms of the Fibroblast Growth Factor Receptor 1 (FGFR1) gene in a human glioblastoma cancer model. Normal glial cells express the FGFR1α, which contains three extracellular domains. In tumors the most abundant isoform is the FGFR1β, which lacks the first extracellular domain due to the skipping of a single exon, termed alpha. The skipping of the α-exon is regulated by two intronic silencing sequences within the precursor mRNA. Since we observed no mutations on these elements in tumor cells, we hypothesized that the over-expression of regulatory proteins that recognize these sequences is responsible for the aberrant expression of splicing isoforms. Hence, we blocked the formation of protein complexes on the ISS using antisense RNA oligonucleotides in vitro. We also evaluated the impact of the ISS antisense oligonucleotides on the endogenous FGFR1 splicing, in a glioblastoma cell model. By targeting intronic regulatory elements we were able to increase the level of alpha exon inclusion up to 90% in glioblastoma cells. The effect was dose dependent, sequence specific and reproducible in glioblastoma and other cancer cells, which also exhibit an alpha exon skipping phenotype. Targeting FGFR1 endogenous ISS1 and ISS2 sequences did not have an additive or synergistic effect, which suggest a regulatory splicing mechanism that requires the interaction of complexes formed on these elements. An increase in the levels of the FGFR1α isoform resulted in a reduction in cell invasiveness. Also, a significant increase in the levels of caspase 3/7 activities, which is indicative of an elevation in apoptosis levels, suggests that expression of FGFR1β might be relevant for tumor survival. These studies demonstrate that it is possible to prevent aberrant expression of exon skipping events through the targeting of intronic regulatory elements, providing an important new therapeutic tool for the correction of human disease caused by alternative RNA splicing. ^

Application of stellar photometry to the analysis of microarray images

Improvements in the analysis of microarray images are critical for accurately quantifying gene expression levels. The acquisition of accurate spot intensities directly influences the results and interpretation of statistical analyses. This dissertation discusses the implementation of a novel approach to the analysis of cDNA microarray images. We use a stellar photometric model, the Moffat function, to quantify microarray spots from nylon microarray images. The inherent flexibility of the Moffat shape model makes it ideal for quantifying microarray spots. We apply our novel approach to a Wilms' tumor microarray study and compare our results with a fixed-circle segmentation approach for spot quantification. Our results suggest that different spot feature extraction methods can have an impact on the ability of statistical methods to identify differentially expressed genes. We also used the Moffat function to simulate a series of microarray images under various experimental conditions. These simulations were used to validate the performance of various statistical methods for identifying differentially expressed genes. Our simulation results indicate that tests taking into account the dependency between mean spot intensity and variance estimation, such as the smoothened t-test, can better identify differentially expressed genes, especially when the number of replicates and mean fold change are low. The analysis of the simulations also showed that overall, a rank sum test (Mann-Whitney) performed well at identifying differentially expressed genes. Previous work has suggested the strengths of nonparametric approaches for identifying differentially expressed genes. We also show that multivariate approaches, such as hierarchical and k-means cluster analysis along with principal components analysis, are only effective at classifying samples when replicate numbers and mean fold change are high. Finally, we show how our stellar shape model approach can be extended to the analysis of 2D-gel images by adapting the Moffat function to take into account the elliptical nature of spots in such images. Our results indicate that stellar shape models offer a previously unexplored approach for the quantification of 2D-gel spots. ^

Respiratory rhythm generation: A modeling study of the respiratory central pattern generator and the phrenic motor neuron

The respiratory central pattern generator is a collection of medullary neurons that generates the rhythm of respiration. The respiratory central pattern generator feeds phrenic motor neurons, which, in turn, drive the main muscle of respiration, the diaphragm. The purpose of this thesis is to understand the neural control of respiration through mathematical models of the respiratory central pattern generator and phrenic motor neurons. ^ We first designed and validated a Hodgkin-Huxley type model that mimics the behavior of phrenic motor neurons under a wide range of electrical and pharmacological perturbations. This model was constrained physiological data from the literature. Next, we designed and validated a model of the respiratory central pattern generator by connecting four Hodgkin-Huxley type models of medullary respiratory neurons in a mutually inhibitory network. This network was in turn driven by a simple model of an endogenously bursting neuron, which acted as the pacemaker for the respiratory central pattern generator. Finally, the respiratory central pattern generator and phrenic motor neuron models were connected and their interactions studied. ^ Our study of the models has provided a number of insights into the behavior of the respiratory central pattern generator and phrenic motor neurons. These include the suggestion of a role for the T-type and N-type calcium channels during single spikes and repetitive firing in phrenic motor neurons, as well as a better understanding of network properties underlying respiratory rhythm generation. We also utilized an existing model of lung mechanics to study the interactions between the respiratory central pattern generator and ventilation. ^