Using spatial linear models with SAR and CAR structure to examine Texas lung cancer incidence rates

Scholars have found that socioeconomic status was one of the key factors that influenced early-stage lung cancer incidence rates in a variety of regions. This thesis examined the association between median household income and lung cancer incidence rates in Texas counties. A total of 254 individual counties in Texas with corresponding lung cancer incidence rates from 2004 to 2008 and median household incomes in 2006 were collected from the National Cancer Institute Surveillance System. A simple linear model and spatial linear models with two structures, Simultaneous Autoregressive Structure (SAR) and Conditional Autoregressive Structure (CAR), were used to link median household income and lung cancer incidence rates in Texas. The residuals of the spatial linear models were analyzed with Moran's I and Geary's C statistics, and the statistical results were used to detect similar lung cancer incidence rate clusters and disease patterns in Texas.^

Risk of Second Malignant Neoplasms Following Proton Arc Therapy and Volumetric Modulated Arc Therapy for Prostate Cancer

The risk of second malignant neoplasms (SMNs) following prostate radiotherapy is a concern due to the large population of survivors and decreasing age at diagnosis. It is known that parallel-opposed beam proton therapy carries a lower risk than photon IMRT. However, a comparison of SMN risk following proton and photon arc therapies has not previously been reported. The purpose of this study was to predict the ratio of excess relative risk (RRR) of SMN incidence following proton arc therapy to that after volumetric modulated arc therapy (VMAT). Additionally, we investigated the impact of margin size and the effect of risk-minimized proton beam weighting on predicted RRR. Physician-approved treatment plans were created for both modalities for three patients. Therapeutic dose was obtained with differential dose-volume histograms from the treatment planning system, and stray dose was estimated from the literature or calculated with Monte Carlo simulations. Then, various risk models were applied to the total dose. Additional treatment plans were also investigated with varying margin size and risk-minimized proton beam weighting. The mean RRR ranged from 0.74 to 0.99, depending on risk model. The additional treatment plans revealed that the RRR remained approximately constant with varying margin size, and that the predicted RRR was reduced by 12% using a risk-minimized proton arc therapy planning technique. In conclusion, proton arc therapy was found to provide an advantage over VMAT in regard to predicted risk of SMN following prostate radiotherapy. This advantage was independent of margin size and was amplified with risk-optimized proton beam weighting.

Sulforaphane improves cognitive function administered following traumatic brain injury.

Recent studies have shown that sulforaphane, a naturally occurring compound that is found in cruciferous vegetables, offers cellular protection in several models of brain injury. When administered following traumatic brain injury (TBI), sulforaphane has been demonstrated to attenuate blood-brain barrier permeability and reduce cerebral edema. These beneficial effects of sulforaphane have been shown to involve induction of a group of cytoprotective, Nrf2-driven genes, whose protein products include free radical scavenging and detoxifying enzymes. However, the influence of sulforaphane on post-injury cognitive deficits has not been examined. In this study, we examined if sulforaphane, when administered following cortical impact injury, can improve the performance of rats tested in hippocampal- and prefrontal cortex-dependent tasks. Our results indicate that sulforaphane treatment improves performance in the Morris water maze task (as indicated by decreased latencies during learning and platform localization during a probe trial) and reduces working memory dysfunction (tested using the delayed match-to-place task). These behavioral improvements were only observed when the treatment was initiated 1h, but not 6h, post-injury. These studies support the use of sulforaphane in the treatment of TBI, and extend the previously observed protective effects to include enhanced cognition.

The Role of Cell Sterilization in Population Based Studies of Radiogenic Second Cancers Following Radiation Therapy

Advances in radiotherapy have generated increased interest in comparative studies of treatment techniques and their effectiveness. In this respect, pediatric patients are of specific interest because of their sensitivity to radiation induced second cancers. However, due to the rarity of childhood cancers and the long latency of second cancers, large sample sizes are unavailable for the epidemiological study of contemporary radiotherapy treatments. Additionally, when specific treatments are considered, such as proton therapy, sample sizes are further reduced due to the rareness of such treatments. We propose a method to improve statistical power in micro clinical trials. Specifically, we use a more biologically relevant quantity, cancer equivalent dose (DCE), to estimate risk instead of mean absorbed dose (DMA). Our objective was to demonstrate that when DCE is used fewer subjects are needed for clinical trials. Thus, we compared the impact of DCE vs. DMA on sample size in a virtual clinical trial that estimated risk for second cancer (SC) in the thyroid following craniospinal irradiation (CSI) of pediatric patients using protons vs. photons. Dose reconstruction, risk models, and statistical analysis were used to evaluate SC risk from therapeutic and stray radiation from CSI for 18 patients. Absorbed dose was calculated in two ways: with (1) traditional DMA and (2) with DCE. DCE and DMA values were used to estimate relative risk of SC incidence (RRCE and RRMA, respectively) after proton vs. photon CSI. Ratios of RR for proton vs. photon CSI (RRRCE and RRRMA) were then used in comparative estimations of sample size to determine the minimal number of patients needed to maintain 80% statistical power when using DCE vs. DMA. For all patients, we found that protons substantially reduced the risk of developing a second thyroid cancer when compared to photon therapy. Mean RRR values were 0.052±0.014 and 0.087±0.021 for RRRMA and RRRCE, respectively. However, we did not find that use of DCE reduced the number of patents needed for acceptable statistical power (i.e, 80%). In fact, when considerations were made for RRR values that met equipoise requirements and the need for descriptive statistics, the minimum number of patients needed for a micro-clinical trial increased from 17 using DMA to 37 using DCE. Subsequent analyses revealed that for our sample, the most influential factor in determining variations in sample size was the experimental standard deviation of estimates for RRR across the patient sample. Additionally, because the relative uncertainty in dose from proton CSI was so much larger (on the order of 2000 times larger) than the other uncertainty terms, it dominated the uncertainty in RRR. Thus, we found that use of corrections for cell sterilization, in the form of DCE, may be an important and underappreciated consideration in the design of clinical trials and radio-epidemiological studies. In addition, the accurate application of cell sterilization to thyroid dose was sensitive to variations in absorbed dose, especially for proton CSI, which may stem from errors in patient positioning, range calculation, and other aspects of treatment planning and delivery.

Mouse models and quantitative trait loci of susceptibility to bleomycin- and radiation-induced pulmonary fibrosis

Radiotherapy involving the thoracic cavity and chemotherapy with the drug bleomycin are both dose limited by the development of pulmonary fibrosis. From evidence that there is variation in the population in susceptibility to pulmonary fibrosis, and animal data, it was hypothesized that individual variation in susceptibility to bleomycin-induced, or radiation-induced, pulmonary fibrosis is, in part, genetically controlled. In this thesis a three generation mouse genetic model of C57BL/6J (fibrosis prone) and C3Hf/Kam (fibrosis resistant) mouse strains and F1 and F2 (F1 intercross) progeny derived from the parental strains was developed to investigate the genetic basis of susceptibility to fibrosis. In the bleomycin studies the mice received 100 mg/kg (125 for females) of bleomycin, via mini osmotic pump. The animals were sacrificed at eight weeks following treatment or when their breathing rate indicated respiratory distress. In the radiation studies the mice were given a single dose of 14 or 16 Gy (Co$\sp{60})$ to the whole thorax and were sacrificed when moribund. The phenotype was defined as the percent of fibrosis area in the left lung as quantified with image analysis of histological sections. Quantitative trait loci (QTL) mapping was used to identify the chromosomal location of genes which contribute to susceptibility to bleomycin-induced pulmonary fibrosis in C57BL/6J mice compared to C3Hf/Kam mice and to determine if the QTL's which influence susceptibility to bleomycin-induced lung fibrosis in these progenitor strains could be implicated in susceptibility to radiation-induced lung fibrosis. For bleomycin, a genome wide scan revealed QTL's on chromosome 17, at the MHC, (LOD = 11.7 for males and 7.2 for females) accounting for approximately 21% of the phenotypic variance, and on chromosome 11 (LOD = 4.9), in male mice only, adding 8% of phenotypic variance. The bleomycin QTL on chromosome 17 was also implicated for susceptibility to radiation-induced fibrosis (LOD = 5.0) and contributes 7% of the phenotypic variance in the radiation study. In conclusion, susceptibility to both bleomycin-induced and radiation-induced pulmonary fibrosis are heritable traits, and are influenced by a genetic factor which maps to a genomic region containing the MHC. ^

MEASUREMENT OF TUMOR BLOOD FLOW FOLLOWING NEUTRON IRRADIATION (ACTIVATION)

Clinical oncologists and cancer researchers benefit from information on the vascularization or non-vascularization of solid tumors because of blood flow's influence on three popular treatment types: hyperthermia therapy, radiotherapy, and chemotherapy. The objective of this research is the development of a clinically useful tumor blood flow measurement technique. The designed technique is sensitive, has good spatial resolution, in non-invasive and presents no risk to the patient beyond his usual treatment (measurements will be subsequent only to normal patient treatment).^ Tumor blood flow was determined by measuring the washout of positron emitting isotopes created through neutron therapy treatment. In order to do this, several technical and scientific questions were addressed first. These questions were: (1) What isotopes are created in tumor tissue when it is irradiated in a neutron therapy beam and how much of each isotope is expected? (2) What are the chemical states of the isotopes that are potentially useful for blood flow measurements and will those chemical states allow these or other isotopes to be washed out of the tumor? (3) How should isotope washout by blood flow be modeled in order to most effectively use the data? These questions have been answered through both theoretical calculation and measurement.^ The first question was answered through the measurement of macroscopic cross sections for the predominant nuclear reactions in the body. These results correlate well with an independent mathematical prediction of tissue activation and measurements of mouse spleen neutron activation. The second question was addressed by performing cell suspension and protein precipitation techniques on neutron activated mouse spleens. The third and final question was answered by using first physical principles to develop a model mimicking the blood flow system and measurement technique.^ In a final set of experiments, the above were applied to flow models and animals. The ultimate aim of this project is to apply its methodology to neutron therapy patients. ^

Tumor heterogeneity dictates sensitivity to gefitinib (Iressa(TM)/ZD1839) in solid tumor models

The epidermal growth factor receptor (EGFR) and its ligands are overexpressed in many human tumors, including bladder and pancreas, correlating with a more aggressive tumor phenotype and poor patient prognosis. We initiated the present study to characterize the heterogeneity of gefitinib responsiveness in a panel of human bladder and pancreatic cancer cell lines in order to identify the biological characteristics of EGFR-dependent proliferation that could be used to prospectively identify drug-sensitive tumors. A second objective was to elucidate how to best exploit these results by utilizing gefitinib in combination therapy. To these ends, we examined the effects of the EGFR antagonist gefitinib on proliferation and apoptosis in a panel of 18 human bladder cancer cell lines and 9 human pancreatic cancer cell lines. Our data confirmed the existence of marked heterogeneity in Iressa responsiveness with less than half of the cell lines displaying significant growth inhibition by clinically relevant concentrations of the drug. Gefitinib responsiveness was found to be p27 kip1 dependent as DNA synthesis was restored following exposure to p27siRNA. Unfortunately, Iressa responsiveness was not closely linked to surface EGFR or TGF-α expression in the bladder cancer cells, however, cellular TGF-α expression correlated directly with Iressa sensitivity in the pancreatic cancer cell lines. These findings provide the potential for prospectively identifying patients with drug-sensitive tumors. ^ Further studies aimed at exploiting gefitinib-mediated cell cycle effects led us to investigate if gefitinib-mediated TRAIL sensitization correlated with increased p27kip1 accumulation. We observed that increased TRAIL sensitivity following gefitinib exposure was not dependent on p27 kip1 expression. Additional studies initiated to examine the role(s) of Akt and Erk signaling demonstrated that exposure to PI3K or MEK inhibitors significantly enhanced TRAIL-induced apoptosis at concentrations that block target phosphorylation. Furthermore, combinations of TRAIL and the PI3K or MEK inhibitors increased procaspase-8 processing above levels observed with TRAIL alone, indicating that the effects were exerted at the level of caspase-8 activation, considered the earliest step in the TRAIL pathway. ^

Bayesian generalized linear models for meta-analysis of diagnostic tests

With the recognition of the importance of evidence-based medicine, there is an emerging need for methods to systematically synthesize available data. Specifically, methods to provide accurate estimates of test characteristics for diagnostic tests are needed to help physicians make better clinical decisions. To provide more flexible approaches for meta-analysis of diagnostic tests, we developed three Bayesian generalized linear models. Two of these models, a bivariate normal and a binomial model, analyzed pairs of sensitivity and specificity values while incorporating the correlation between these two outcome variables. Noninformative independent uniform priors were used for the variance of sensitivity, specificity and correlation. We also applied an inverse Wishart prior to check the sensitivity of the results. The third model was a multinomial model where the test results were modeled as multinomial random variables. All three models can include specific imaging techniques as covariates in order to compare performance. Vague normal priors were assigned to the coefficients of the covariates. The computations were carried out using the 'Bayesian inference using Gibbs sampling' implementation of Markov chain Monte Carlo techniques. We investigated the properties of the three proposed models through extensive simulation studies. We also applied these models to a previously published meta-analysis dataset on cervical cancer as well as to an unpublished melanoma dataset. In general, our findings show that the point estimates of sensitivity and specificity were consistent among Bayesian and frequentist bivariate normal and binomial models. However, in the simulation studies, the estimates of the correlation coefficient from Bayesian bivariate models are not as good as those obtained from frequentist estimation regardless of which prior distribution was used for the covariance matrix. The Bayesian multinomial model consistently underestimated the sensitivity and specificity regardless of the sample size and correlation coefficient. In conclusion, the Bayesian bivariate binomial model provides the most flexible framework for future applications because of its following strengths: (1) it facilitates direct comparison between different tests; (2) it captures the variability in both sensitivity and specificity simultaneously as well as the intercorrelation between the two; and (3) it can be directly applied to sparse data without ad hoc correction. ^

Quantitative comparison of late effects following photon versus proton external-beam radiation therapies: Toward an evidence-based approach to selecting a treatment modality

Radiation therapy has been used as an effective treatment for malignancies in pediatric patients. However, in many cases, the side effects of radiation diminish these patients’ quality of life. In order to develop strategies to minimize radiogenic complications, one must first quantitatively estimate pediatric patients’ relative risk for radiogenic late effects, which has not become feasible till recently because of the calculational complexity. The goals of this work were to calculate the dose delivered to tissues and organs in pediatric patients during contemporary photon and proton radiotherapies; to estimate the corresponding risk of radiogenic second cancer and cardiac toxicity based on the calculated doses and on dose-risk models from the literature; to test for the statistical significance of the difference between predicted risks after photon versus proton radiotherapies; and to provide a prototype of an evidence-based approach to selecting treatment modalities for pediatric patients, taking second cancer and cardiac toxicity into account. The results showed that proton therapy confers a lower predicted risk of radiogenic second cancer, and lower risks of radiogenic cardiac toxicities, compared to photon therapy. An uncertainty analysis revealed that the qualitative findings of this study are insensitive to changes in a wide variety of host and treatment related factors.

Novel Use of Dual Anti-Inflammatory Therapy to Overcome Drug Resistance and Improve Functional Recovery Following Spinal Cord Injury

Over 1.2 million Americans are currently living with a traumatic spinal cord injury (SCI). Despite the need for effective therapies, there are currently no proven effective treatments that can improve recovery of function in SCI patients. Many therapeutic compounds have shown promise in preclinical models of SCI, but all of these have fallen short in clinical trials. P-glycoprotein (Pgp) is an active transporter expressed on capillary endothelial cell membranes at the blood-spinal cord barrier (BSCB). Pgp limits passive diffusion of blood-borne drugs into the CNS, by actively extruding drugs from the endothelial cell membrane. Pgp can become pathologically up-regulated, thus greatly impeding therapeutic drug delivery (‘multidrug resistance’). Importantly, many drugs that have been evaluated for the treatment of SCI are Pgp substrates. We hypothesized that Pgp-mediated drug resistance diminishes the delivery and efficacy of neuroprotective drugs following SCI. We observed a progressive, spatial spread of Pgp overexpression within the injured spinal cord. To assess Pgp function, we examined spinal cord uptake of systemically-delivered riluzole, a drug that is currently being evaluated in clinical trials as an SCI intervention. Blood-to-spinal cord riluzole penetration was reduced following SCI in wild-type but not Pgp-null rats, highlighting a critical role for Pgp in mediating spinal cord drug resistance after injury. Others have shown that pro-inflammatory signaling drives Pgp up-regulation in cancer and epilepsy. We have detected inflammation in both acutely- and chronically-injured spinal cord tissue. We therefore evaluated the ability of the dual COX-/5-LOX inhibitor licofelone to attenuate Pgp-mediated drug resistance following SCI. Licofelone treatment both reduced spinal cord Pgp levels and enhanced spinal cord riluzole bioavailability following SCI. Thus, we propose that licofelone may offer a new combinatorial treatment strategy to enhance spinal cord drug delivery following SCI. Additionally, we assessed the ability of licofelone, riluzole, or both to enhance recovery of locomotor function following SCI. We found that licofelone treatment conferred a significant improvement in hindlimb function that was sustained through the end of the study. In contrast, riluzole did not improve functional outcome. We therefore conclude that licofelone holds promise as a potential neuroprotective intervention for SCI.

Analysis of life expectancy using linear and nonlinear models

Life expectancy has consistently increased over the last 150 years due to improvements in nutrition, medicine, and public health. Several studies found that in many developed countries, life expectancy continued to rise following a nearly linear trend, which was contrary to a common belief that the rate of improvement in life expectancy would decelerate and was fit with an S-shaped curve. Using samples of countries that exhibited a wide range of economic development levels, we explored the change in life expectancy over time by employing both nonlinear and linear models. We then observed if there were any significant differences in estimates between linear models, assuming an auto-correlated error structure. When data did not have a sigmoidal shape, nonlinear growth models sometimes failed to provide meaningful parameter estimates. The existence of an inflection point and asymptotes in the growth models made them inflexible with life expectancy data. In linear models, there was no significant difference in the life expectancy growth rate and future estimates between ordinary least squares (OLS) and generalized least squares (GLS). However, the generalized least squares model was more robust because the data involved time-series variables and residuals were positively correlated. ^