IMPROVING QUANTITATIVE TREATMENT RESPONSE MONITORING WITH DEFORMABLE IMAGE REGISTRATION

Quantitative imaging with 18F-FDG PET/CT has the potential to provide an in vivo assessment of response to radiotherapy (RT). However, comparing tissue tracer uptake in longitudinal studies is often confounded by variations in patient setup and potential treatment induced gross anatomic changes. These variations make true response monitoring for the same anatomic volume a challenge, not only for tumors, but also for normal organs-at-risk (OAR). The central hypothesis of this study is that more accurate image registration will lead to improved quantitation of tissue response to RT with 18F-FDG PET/CT. Employing an in-house developed “demons” based deformable image registration algorithm, pre-RT tumor and parotid gland volumes can be more accurately mapped to serial functional images. To test the hypothesis, specific aim 1 was designed to analyze whether deformably mapping tumor volumes rather than aligning to bony structures leads to superior tumor response assessment. We found that deformable mapping of the most metabolically avid regions improved response prediction (P<0.05). The positive predictive power for residual disease was 63% compared to 50% for contrast enhanced post-RT CT. Specific aim 2 was designed to use parotid gland standardized uptake value (SUV) as an objective imaging biomarker for salivary toxicity. We found that relative change in parotid gland SUV correlated strongly with salivary toxicity as defined by the RTOG/EORTC late effects analytic scale (Spearman’s ρ = -0.96, P<0.01). Finally, the goal of specific aim 3 was to create a phenomenological dose-SUV response model for the human parotid glands. Utilizing only baseline metabolic function and the planned dose distribution, predicting parotid SUV change or salivary toxicity, based upon specific aim 2, became possible. We found that the predicted and observed parotid SUV relative changes were significantly correlated (Spearman’s ρ = 0.94, P<0.01). The application of deformable image registration to quantitative treatment response monitoring with 18F-FDG PET/CT could have a profound impact on patient management. Accurate and early identification of residual disease may allow for more timely intervention, while the ability to quantify and predict toxicity of normal OAR might permit individualized refinement of radiation treatment plan designs.

The immune response to the VlsE antigenic variation protein of Borrelia burgdorferi

The spirochete Borrelia burgdorferi (Bb) is the causative agent of Lyme disease. During infection, a strong immune response is elicited towards Bb by its host; however, the organism is able to persist and to disseminate to many different tissues. The vls locus is located on the linear plasmid lp28-1, a plasmid shown to be important for virulence in the mouse model. During infection, vlsE undergoes antigenic variation through a series of gene conversions, which results in the insertion of sequences from the silent, unexpressed cassettes into the vlsE cassette. We hypothesize that this antigenic variation is important in the spirochete's ability to persist within mammals by allowing it to evade the immune system. To define the role of vls in immune evasion, the immune response against VlsE was determined by using a recombinant form of VlsE (VlsE1-His) as an antigen to screen patient sera. Lyme patients produce antibodies that recognize VlsE, and these antibodies are present throughout the course of disease. Immunization with the VlsE1-His protein provided protection against infection with Bb expressing the same variant of VlsE (VlsE1), but was only partially protective when mice were infected with organisms expressing VlsE variants; however, subsequent VlsE immunization studies yielded inconsistent protection. Successful immunizations produced different antibody reactivities to VlsE epitopes than non-protective immunizations, but the reason for this variable response is unclear. In the process of developing genetic approaches to transform infectious Bb, it was determined that the transformation barrier posed by plasmids lp25 and lp56 could be circumvented by replacing the required lp25 gene pncA. To characterize the role of vlsE in infectivity, Bb lacking lp28-1 were complemented with a shuttle plasmid containing the lp25 encoded virulence determinant pncA and vlsE. Complemented spirochetes express VlsE, but the gene does not undergo antigenic variation and infectivity in the mouse model was not restored, indicating that either antigenic variation of vlsE is necessary for survival in the mouse model or that other genes on lp28-1 are important for virulence. ^

Differential stress response and susceptibility to oxidative injury in alveolar macrophages from C57BL/6J and C3H/HeJ mice following oxidant exposure

Studies have demonstrated a variable response to ozone among individuals and animal species and strains. For instance, C57BL/6J mice have a greater inflammatory response to ozone exposure than C3H/HeJ mice. In these studies, I utilized these strain differences in an effort to derive a mechanistic explanation to the variable strain sensitivity to ozone exposure. Therefore, alveolar macrophages (AM) from C57BL/6J and C3H/HeJ mice were exposed in vitro to hydrogen peroxide ($\rm H\sb2O\sb2$), heat and acetyl ceramide or in vivo to ozone. Necrosis and DNA fragmentation in macrophages from the two murine strains were determined to assess cytotoxicity following these treatments. In addition, synthesis and expression of the stress proteins, stress protein 72 (SP72) and heme oxygenase (HO-1), were examined following treatments. The in vitro experiments were conducted to eliminate the possibility of in vivo confounders (i.e., differences in breathing rates in the two strains) and thus directly implicate some inherent difference between cells from the two murine strains. $\rm H\sb2O\sb2$ and heat caused greater cytotoxicity in AM from C57BL/6J than C3H/HeJ mice and DNA fragmentation was a particularly sensitive indicator of cell injury. Similarly, AM from C57BL/6J mice were more sensitive to ozone exposure than cells from C3H/HeJ mice. Exposure to either 1 or 0.4 ppm ozone caused greater cytotoxicity in macrophages from C57BL/6J mice compared to macrophages from C3H/HeJ mice. The increased sensitivity of AM to injury was associated with decreased synthesis and expression of stress proteins. AM from C57BL/6J mice synthesized and expressed significantly less stress proteins in response to heat and ozone than AM from C3H/HeJ mice. Heat treatment resulted in greater synthesis and expression of SP72. In addition, macrophages from C57BL/6J mice expressed lower amounts of HO-1 than macrophages from C3H/HeJ mice following 0.4 ppm ozone exposure. Therefore, AM from C57BL/6J mice are more susceptible to oxidative injury than AM from C3H/HeJ mice which might be due to differential expression of stress proteins in these cells. ^

Effects of variable Gs expression on $\beta\sb2$-adrenergic receptor stimulated adenylyl cyclase response

An exact knowledge of the kinetic nature of the interaction between the stimulatory G protein (G$\sb{\rm s}$) and the adenylyl cyclase catalytic unit (C) is essential for interpreting the effects of Gs mutations and expression levels on cellular response to a wide variety of hormones, drugs, and neurotransmitters. In particular, insight as to the association of these proteins could lead to progress in tumor biology where single spontaneous mutations in G proteins have been associated with the formation of tumors (118). The question this work attempts to answer is whether the adenylyl cyclase activation by epinephrine stimulated $\beta\sb2$-adrenergic receptors occurs via G$\sb{\rm s}$ proteins by a G$\sb{\rm s}$ to C shuttle or G$\sb{\rm s}$-C precoupled mechanism. The two forms of activation are distinguishable by the effect of G$\sb{\rm s}$ levels on epinephrine stimulated EC50 values for cyclase activation.^ We have made stable transfectants of S49 cyc$\sp-$ cells with the gene for the $\alpha$ protein of G$\sb{\rm s}$ $(\alpha\sb{\rm s})$ which is under the control of the mouse mammary tumor virus LTR promoter (110). Expression of G$\sb{\rm s}\alpha$ was then controlled by incubation of the cells for various times with 5 $\mu$M dexamethasone. Expression of G$\sb{\rm s}\alpha$ led to the appearance of GTP shifts in the competitive binding of epinephrine with $\sp{125}$ICYP to the $\beta$-adrenergic receptors and to agonist dependent adenylyl cyclase activity. High expression of G$\sb{\rm s}\alpha$ resulted in lower EC50's for the adenylyl cyclase activity in response to epinephrine than did low expression. By kinetic modelling, this result is consistent with the existence of a shuttle mechanism for adenylyl cyclase activation by hormones.^ One item of concern that remains to be addressed is the extent to which activation of adenylyl cyclase occurs by a "pure" shuttle mechanism. Kinetic and biochemical experiments by other investigators have revealed that adenylyl cyclase activation, by hormones, may occur via a Gs-C precoupled mechanism (80, 94, 97). Activation of adenylyl cyclase, therefore, probably does not occur by either a pure "'Shuttle" or "Gs-C Precoupled" mechanism, but rather by a "Hybrid" mechanism. The extent to which either the shuttle or precoupled mechanism contributes to hormone stimulated adenylyl cyclase activity is the subject of on-going research. ^

Prediction of sepsis in the burn intensive care unit patient

Sepsis is a significant cause for multiple organ failure and death in the burn patient, yet identification in this population is confounded by chronic hypermetabolism and impaired immune function. The purpose of this study was twofold: 1) determine the ability of the systemic inflammatory response syndrome (SIRS) and American Burn Association (ABA) criteria to predict sepsis in the burn patient; and 2) develop a model representing the best combination of clinical predictors associated with sepsis in the same population. A retrospective, case-controlled, within-patient comparison of burn patients admitted to a single intensive care unit (ICU) was conducted for the period January 2005 to September 2010. Blood culture results were paired with clinical condition: "positive-sick"; "negative-sick", and "screening-not sick". Data were collected for the 72 hours prior to each blood culture. The most significant predictors were evaluated using logistic regression, Generalized Estimating Equations (GEE) and ROC area under the curve (AUC) analyses to assess model predictive ability. Bootstrapping methods were employed to evaluate potential model over-fitting. Fifty-nine subjects were included, representing 177 culture periods. SIRS criteria were not found to be associated with culture type, with an average of 98% of subjects meeting criteria in the 3 days prior. ABA sepsis criteria were significantly different among culture type only on the day prior (p = 0.004). The variables identified for the model included: heart rate>130 beats/min, mean blood pressure<60 mmHg, base deficit<-6 mEq/L, temperature>36°C, use of vasoactive medications, and glucose>150 mg/d1. The model was significant in predicting "positive culture-sick" and sepsis state, with AUC of 0.775 (p < 0.001) and 0.714 (p < .001), respectively; comparatively, the ABA criteria AUC was 0.619 (p = 0.028) and 0.597 (p = .035), respectively. SIRS criteria are not appropriate for identifying sepsis in the burn population. The ABA criteria perform better, but only for the day prior to positive blood culture results. The time period useful to diagnose sepsis using clinical criteria may be limited to 24 hours. A combination of predictors is superior to individual variable trends, yet algorithms or computer support will be necessary for the clinician to find such models useful. ^

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

Transcriptional regulation of the Borrelia burgdorferi antigenically variable VlsE surface protein.

The Lyme disease agent Borrelia burgdorferi can persistently infect humans and other animals despite host active immune responses. This is facilitated, in part, by the vls locus, a complex system consisting of the vlsE expression site and an adjacent set of 11 to 15 silent vls cassettes. Segments of nonexpressed cassettes recombine with the vlsE region during infection of mammalian hosts, resulting in combinatorial antigenic variation of the VlsE outer surface protein. We now demonstrate that synthesis of VlsE is regulated during the natural mammal-tick infectious cycle, being activated in mammals but repressed during tick colonization. Examination of cultured B. burgdorferi cells indicated that the spirochete controls vlsE transcription levels in response to environmental cues. Analysis of PvlsE::gfp fusions in B. burgdorferi indicated that VlsE production is controlled at the level of transcriptional initiation, and regions of 5' DNA involved in the regulation were identified. Electrophoretic mobility shift assays detected qualitative and quantitative changes in patterns of protein-DNA complexes formed between the vlsE promoter and cytoplasmic proteins, suggesting the involvement of DNA-binding proteins in the regulation of vlsE, with at least one protein acting as a transcriptional activator.

Design and analysis for three level hierarchical structures with count response

In numerous intervention studies and education field trials, random assignment to treatment occurs in clusters rather than at the level of observation. This departure of random assignment of units may be due to logistics, political feasibility, or ecological validity. Data within the same cluster or grouping are often correlated. Application of traditional regression techniques, which assume independence between observations, to clustered data produce consistent parameter estimates. However such estimators are often inefficient as compared to methods which incorporate the clustered nature of the data into the estimation procedure (Neuhaus 1993).1 Multilevel models, also known as random effects or random components models, can be used to account for the clustering of data by estimating higher level, or group, as well as lower level, or individual variation. Designing a study, in which the unit of observation is nested within higher level groupings, requires the determination of sample sizes at each level. This study investigates the design and analysis of various sampling strategies for a 3-level repeated measures design on the parameter estimates when the outcome variable of interest follows a Poisson distribution. ^ Results study suggest that second order PQL estimation produces the least biased estimates in the 3-level multilevel Poisson model followed by first order PQL and then second and first order MQL. The MQL estimates of both fixed and random parameters are generally satisfactory when the level 2 and level 3 variation is less than 0.10. However, as the higher level error variance increases, the MQL estimates become increasingly biased. If convergence of the estimation algorithm is not obtained by PQL procedure and higher level error variance is large, the estimates may be significantly biased. In this case bias correction techniques such as bootstrapping should be considered as an alternative procedure. For larger sample sizes, those structures with 20 or more units sampled at levels with normally distributed random errors produced more stable estimates with less sampling variance than structures with an increased number of level 1 units. For small sample sizes, sampling fewer units at the level with Poisson variation produces less sampling variation, however this criterion is no longer important when sample sizes are large. ^ 1Neuhaus J (1993). “Estimation efficiency and Tests of Covariate Effects with Clustered Binary Data”. Biometrics , 49, 989–996^

A novel treatment planning methodology for high dose 166Ho-DOTMP therapy in patients with multiple myeloma

Bone marrow ablation, i.e., the complete sterilization of the active bone marrow, followed by bone marrow transplantation (BMT) is a comment treatment of hematological malignancies. The use of targeted bone-seeking radiopharmaceuticals to selectively deliver radiation to the adjacent bone marrow cavities while sparing normal tissues is a promising technique. Current radiopharmaceutical treatment planning methods do not properly compensate for the patient-specific variable distribution of radioactive material within the skeleton. To improve the current method of internal dosimetry, novel methods for measuring the radiopharmaceutical distribution within the skeleton were developed. 99mTc-MDP was proven as an adequate surrogate for measuring 166Ho-DOTMP skeletal uptake and biodistribution, allowing these measures to be obtained faster, safer, and with higher spatial resolution. This translates directly into better measurements of the radiation dose distribution within the bone marrow. The resulting bone marrow dose-volume histograms allow prediction of the patient disease response where conventional organ scale dosimetry failed. They indicate that complete remission is only achieved when greater than 90% of the bone marrow receives at least 30 Gy. ^ Comprehensive treatment planning requires combining target and non-target organ dosimetry. Organs in the urinary tract were of special concern. The kidney dose is primarily dependent upon the mean transit time of 166 Ho-DOTMP through the kidney. Deconvolution analysis of renograms predicted a mean transit time of 2.6 minutes for 166Ho-DOTMP. The radiation dose to the urinary bladder wall is dependent upon numerous factors including patient hydration and void schedule. For beta-emitting isotopes such as 166Ho, reduction of the bladder wall dose is best accomplished through good patient hydration and ensuring a partially full bladder at the time of injection. Encouraging the patient to void frequently, or catheterizing the patient without irrigation, will not significantly reduce the bladder wall dose. ^ The results from this work will produce the most advanced treatment planning methodology for bone marrow ablation therapy using radioisotopes currently available. Treatments can be tailored specifically for each patient, including the addition of concomitant total body irradiation for patients with unfavorable dose distributions, to deliver a desired patient disease response, while minimizing the dose or toxicity to non-target organs. ^

Logistic regression models for ordinal response: A study of self-efficacy in colorectal cancer screening

The ordinal logistic regression models are used to analyze the dependant variable with multiple outcomes that can be ranked, but have been underutilized. In this study, we describe four logistic regression models for analyzing the ordinal response variable. ^ In this methodological study, the four regression models are proposed. The first model uses the multinomial logistic model. The second is adjacent-category logit model. The third is the proportional odds model and the fourth model is the continuation-ratio model. We illustrate and compare the fit of these models using data from the survey designed by the University of Texas, School of Public Health research project PCCaSO (Promoting Colon Cancer Screening in people 50 and Over), to study the patient’s confidence in the completion colorectal cancer screening (CRCS). ^ The purpose of this study is two fold: first, to provide a synthesized review of models for analyzing data with ordinal response, and second, to evaluate their usefulness in epidemiological research, with particular emphasis on model formulation, interpretation of model coefficients, and their implications. Four ordinal logistic models that are used in this study include (1) Multinomial logistic model, (2) Adjacent-category logistic model [9], (3) Continuation-ratio logistic model [10], (4) Proportional logistic model [11]. We recommend that the analyst performs (1) goodness-of-fit tests, (2) sensitivity analysis by fitting and comparing different models.^

Molecular mechanisms of response in gastrointestinal stromal tumors treated with imatinib mesylate

Gastrointestinal stromal tumors (GIST) represent 80% of sarcoma arising from the GI tract. The inciting event in tumor progression is mutation of the kit or, rarely, platelet derived growth factor receptor-α (PDGFR) gene. These mutations encode ligand independent, constitutively active proteins: Kit or PDGFR. ^ These tumors are notoriously chemo and radio resistant. Historically, patients with advanced disease realized a median overall survival of 9 months. However, with modern management of GIST with imatinib mesylate (Novartis), a small molecule inhibitor of the Kit, PDGFR, and Abl tyrosine kinases, patients now realize a median overall survival greater than 30 months. However, almost half of patients present with surgically resectable GIST and the utility of imatinib in this context has not been prospectively studied. Also, therapeutic benefit of imatinib is variable from patient to patient and alternative targeted therapy is emerging as potential alternatives to imatinib. Thus, elucidating prognostic factors for patients with GIST in the imatinib-era is crucial to providing optimal care to each particular patient. Moreover, the exact mechanism of action of imatinib in GIST is not fully understood. Therefore, physicians find difficulty in accurately predicting which patient will benefit from imatinib, how to assess response to therapy, and the time at which to assess response. ^ I have hypothesized that imatinib is tolerable and clinically beneficial in the context of surgery, VEGF expression and kit non-exon 11 genotypes portend poor survival on imatinib therapy, and imatinib's mechanism of action is in part due to anti-vascular effects and inhibition of the Kit/SCF signaling axis of tumor-associated endothelial cells. ^ Results herein demonstrate that imatinib is safe and increases the duration of disease-free survival when combined with surgery. Radiographic and molecular (namely, apoptosis) changes occur within 3 days of imatinib initiation. I illustrate that non-exon 11 mutant genotypes and VEGF are poor prognostic factors for patients treated with imatinib. These findings may allow for patient stratification to emerging therapies rather than imatinib. I show that imatinib has anti-vascular effects via inducing tumor endothelial cell apoptosis perhaps by abrogation of the Kit/SCF signaling axis. ^

The response of Texas public health agencies and workers to the emergency preparedness initiatives after 9/11

The purpose of this research was to better understand the impact of the terrorist attacks in 2001 on public health, particularly for Texas public health. This study employed mixed methods to examine changes to public health culture within Texas local public health agencies, important attitudes of public health workers toward responding to a disaster, and the funding policies that might ensure our investment in public health emergency preparedness is protected. ^ A qualitative analysis of interviews conducted with a large sample of public health officials in Texas found that all the constituent parts of a peculiar culture for public health preparedness existed that spanned the state's local health departments regardless of size, or funding level. The new preparedness culture in Texas had the hallmarks necessary for a robust public health preparedness and emergency response system. ^ The willingness of public health workers, necessary to make these kinds of changes and mount a disaster response was examined in one of Texas' most experienced disaster response teams—the public health workers for the City of Houston. A hypothesized latent variable model showed that willingness mediated all other factors in the model (self-efficacy, knowledge, barriers, and risk perception) for self-reported likelihood of reporting to work for a disaster. The RMSEA for the final model was 0.042 with a confidence interval of 0.036—0.049 and the chi-squared difference test was P=0.08, indicating a well-fitted model that suggests willingness is an important factor for consideration by preparedness planners and researchers alike. ^ Finally, with disasters on the rise and federal funding for preparedness dwindling, a review of states' policies for the distribution of these funds and their advantages and disadvantages were examined through a review of current literature and public documents, and a survey of state-level public health officials, emergency management professionals and researchers. Although the base plus per-capita method is the most common, it is not necessarily perceived to be the most effective. No clear "optimal" method emerged from the study, but recommendations for a strategic combination of three methods were made that has the potential to maximize the benefits of each method, while minimizing the weaknesses.^

Radiomics of NSCLC: Quantitative CT Image Feature Characterization and Tumor Shrinkage Prediction

Radiomics is the high-throughput extraction and analysis of quantitative image features. For non-small cell lung cancer (NSCLC) patients, radiomics can be applied to standard of care computed tomography (CT) images to improve tumor diagnosis, staging, and response assessment. The first objective of this work was to show that CT image features extracted from pre-treatment NSCLC tumors could be used to predict tumor shrinkage in response to therapy. This is important since tumor shrinkage is an important cancer treatment endpoint that is correlated with probability of disease progression and overall survival. Accurate prediction of tumor shrinkage could also lead to individually customized treatment plans. To accomplish this objective, 64 stage NSCLC patients with similar treatments were all imaged using the same CT scanner and protocol. Quantitative image features were extracted and principal component regression with simulated annealing subset selection was used to predict shrinkage. Cross validation and permutation tests were used to validate the results. The optimal model gave a strong correlation between the observed and predicted shrinkages with . The second objective of this work was to identify sets of NSCLC CT image features that are reproducible, non-redundant, and informative across multiple machines. Feature sets with these qualities are needed for NSCLC radiomics models to be robust to machine variation and spurious correlation. To accomplish this objective, test-retest CT image pairs were obtained from 56 NSCLC patients imaged on three CT machines from two institutions. For each machine, quantitative image features with concordance correlation coefficient values greater than 0.90 were considered reproducible. Multi-machine reproducible feature sets were created by taking the intersection of individual machine reproducible feature sets. Redundant features were removed through hierarchical clustering. The findings showed that image feature reproducibility and redundancy depended on both the CT machine and the CT image type (average cine 4D-CT imaging vs. end-exhale cine 4D-CT imaging vs. helical inspiratory breath-hold 3D CT). For each image type, a set of cross-machine reproducible, non-redundant, and informative image features was identified. Compared to end-exhale 4D-CT and breath-hold 3D-CT, average 4D-CT derived image features showed superior multi-machine reproducibility and are the best candidates for clinical correlation.