Application of cell line based genomic predictors to predict response to targeted therapies in breast cancer

Cancer cell lines can be treated with a drug and the molecular comparison of responders and non-responders may yield potential predictors that could be tested in the clinic. It is a bioinformatics challenge to apply the cell line-derived multivariable response predictors to patients who respond to therapy. Using the gene expression data from 23 breast cancer cell lines, I developed three predictors of dasatinib sensitivity by selecting differentially expressed genes and applying different classification algorithms. The performance of these predictors on independent cell lines with known dasatinib response was tested. The predictor based on weighted voting method has the best overall performance. It correctly predicted dasatinib sensitivity in 11 out of 12 (92%) breast and 17 out of 23 (74%) lung cancer cell lines. These predictors were then applied to the gene expression data from 133 breast cancer patients in an attempt to predict how the patients might respond to dasatinib therapy. Two predictors identified 13 patients in common to be dasatinib sensitive. Sixty two percent of these cases are triple negative (ER-negative, HER2-negative and PR-negative) and 76% are double negative. The result is consistent with the findings from other studies, which identified a target population for dasatinib treatment to be triple negative or basal breast cancer subtype. In conclusion, we think that the cell line-derived dasatinib classifiers can be applied to the human patients. ^

Transcriptional regulation of lung tumor suppressor GPRC5A in malignant and normal cells

Chronic exposure of the airways to cigarette smoke induces inflammatory response and genomic instability that play important roles in lung cancer development. Nuclear factor kappa B (NF-κB), the major intracellular mediator of inflammatory signals, is frequently activated in preneoplastic and malignant lung lesions. ^ Previously, we had shown that a lung tumor suppressor GPRC5A is frequently repressed in human non-small cell lung cancers (NSCLC) cells and lung tumor specimens. Recently, other groups have shown that human GPRC5A transcript levels are higher in bronchial samples of former than of current smokers. These results suggested that smoking represses GPRC5A expression and thus promotes the occurrence of lung cancer. We hypothesized that cigarette smoking or associated inflammatory response repressed GPRC5A expression through NF-κB signaling. ^ To determine the effect of inflammation, we examined GPRC5A protein expression in several lung cell lines following by TNF-α treatment. TNF-α significantly suppressed GPRC5A expression in normal small airway epithelial cells (SAEC) as well as in Calu-1 cells. Real-time PCR analysis indicated that TNF-α inhibits GPRC5A expression at the transcriptional level. NF-κB, the major downstream effectors of TNF-α signaling, mediates TNF-α-induced repression of GPRC5A because over-expression of NF-κB suppressed GPRC5A. To determine the region in the GPRC5A promoter through which NF-κB acts, we examined the ability of TNF-α to inhibit a series of reporter constructs with different deletions of GPRC5A promoter. The luciferase assay showed that the potential NF-κB binding sites containing region are irresponsible for TNF-α-induced suppression. Further analysis using constructs with different deletions in p65 revealed that NF-κB-mediated repression of GPRC5A is transcription-independent. Co-immunoprecipitation assays revealed that NF-κB could form a complex with RAR/RXR heterodimer. Moreover, the inhibitory effect of NF-κB has been found to be proportional to NF-κB/RAR ratio in luciferase assay. Finally, Chromatin IP demonstrated that NF-κB/p65 bound to GPRC5A promoter as well as RAR/RXR and suppressed transcription. Taken together, we propose that inflammation-induced NF-κB activation disrupts the RA signaling and suppresses GPRC5A expression and thus contributes to the oncogenesis of lung cancer. Our studies shed new light on the pathogenesis of lung cancer and potentially provide novel interventions for preventing and treating this disease. ^

AN ASSESSMENT OF THE RELATION OF CANCER MORTALITY RATES AND SOURCES OF DRINKING WATER IN TEXAS

The geographic distribution of average annual age-adjusted mortality rates (1964-1976) for four types of cancer (all cancer sites combined, gastrointestinal, urinary, and lung cancer) were compared by sources of drinking water for 254 Texas counties and county rural areas and 301 Texas cities. Exposure variables considered were surface versus ground water, public water supplies versus individuals wells, and trihalomethane levels in municipal water supplies. Each general source of "surface" and "ground" water was further divided by aggregating ground water using areas by aquifers and surface water using study areas by river basins. Potential confounding variables taken into account included median education, employment in cancer risk industries, population mobility, ethnicity, and urbanicity. A pattern of higher and lower cancer mortality rates was found for populations using some aquifers and river basins. Further study is required to determine whether the differences in cancer mortality rates that were observed are related to drinking water content or are coincidental with differences in personal characteristics which could not be taken into account in this ecologic study design. ^

Targeting TRAF6 for Cancer Therapeutical Development

Tumor necrosis factor (TNF)-Receptor Associated Factors (TRAFs) are a family of signal transducer proteins. TRAF6 is a unique member of this family in that it is involved in not only the TNF superfamily, but the toll-like receptor (TLR)/IL-1R (TIR) superfamily. The formation of the complex consisting of Receptor Activator of Nuclear Factor κ B (RANK), with its ligand (RANKL) results in the recruitment of TRAF6, which activates NF-κB, JNK and MAP kinase pathways. TRAF6 is critical in signaling with leading to release of various growth factors in bone, and promotes osteoclastogenesis. TRAF6 has also been implicated as an oncogene in lung cancer and as a target in multiple myeloma. In the hopes of developing small molecule inhibitors of the TRAF6-RANK interaction, multiple steps were carried out. Computational prediction of hot spot residues on the protein-protein interaction of TRAF6 and RANK were examined. Three methods were used: Robetta, KFC2, and HotPoint, each of which uses a different methodology to determine if a residue is a hot spot. These hot spot predictions were considered the basis for resolving the binding site for in silico high-throughput screening using GOLD and the MyriaScreen database of drug/lead-like compounds. Computationally intensive molecular dynamics simulations highlighted the binding mechanism and TRAF6 structural changes upon hit binding. Compounds identified as hits were verified using a GST-pull down assay, comparing inhibition to a RANK decoy peptide. Since many drugs fail due to lack of efficacy and toxicity, predictive models for the evaluation of the LD50 and bioavailability of our TRAF6 hits, and these models can be used towards other drugs and small molecule therapeutics as well. Datasets of compounds and their corresponding bioavailability and LD50 values were curated based, and QSAR models were built using molecular descriptors of these compounds using the k-nearest neighbor (k-NN) method, and quality of these models were cross-validated.

Fibulin-2 stabilizes tumor extracellular matrix and drives malignant progression of lung adenocarcinoma

The ECM of epithelial carcinomas undergoes structural remodeling during periods of uncontrolled growth, creating regional heterogeneity and torsional stress. How tumors maintain ECM integrity in the face of dynamic biophysical forces is still largely unclear. This study addresses these deficiencies using mouse models of human lung adenocarcinoma. Spontaneous lung tumors were marked by disorganized basement membranes, dense collagen networks, and increased tissue stiffness. Metastasis-prone lung adenocarcinoma cells secreted fibulin-2 (Fbln2), a matrix glycoprotein involved in ECM supra-molecular assembly. Fibulin-2 depletion in tumor cells decreased the intra-tumoral abundance of matrix metalloproteinases and reduced collagen cross-linking and tumor compressive properties resulting in inhibited tumor growth and metastasis. Fbln2 deposition within intra-tumoral fibrotic bands was a predictor of poor clinical outcome in patients. Collectively, these findings support a feed-forward model in which tumor cells secrete matrix-stabilizing factors required for the assembly of ECM that preferentially favors malignant progression. To our knowledge, this is the first evidence that tumor cells directly regulate the integrity of their surrounding matrix through the secretion of matrix-stabilizing factors such as fibulin-2. These findings open a new avenue of research into matrix assembly molecules as potential therapeutic targets in cancer patients.

Comparison of Tumor Shrinkage and Cumulative Dose Distribution for Lung Cancers

Background: The physical characteristic of protons is that they deliver most of their radiation dose to the target volume and deliver no dose to the normal tissue distal to the tumor. Previously, numerous studies have shown unique advantages of proton therapy over intensity-modulated radiation therapy (IMRT) in conforming dose to the tumor and sparing dose to the surrounding normal tissues and the critical structures in many clinical sites. However, proton therapy is known to be more sensitive to treatment uncertainties such as inter- and intra-fractional variations in patient anatomy. To date, no study has clearly demonstrated the effectiveness of proton therapy compared with the conventional IMRT under the consideration of both respiratory motion and tumor shrinkage in non-small cell lung cancer (NSCLC) patients. Purpose: This thesis investigated two questions for establishing a clinically relevant comparison of the two different modalities (IMRT and proton therapy). The first question was whether or not there are any differences in tumor shrinkage between patients randomized to IMRT versus passively scattered proton therapy (PSPT). Tumor shrinkage is considered a standard measure of radiation therapy response that has been widely used to gauge a short-term progression of radiation therapy. The second question was whether or not there are any differences between the planned dose and 5D dose under the influence of inter- and intra-fractional variations in the patient anatomy for both modalities. Methods: A total of 45 patients (25 IMRT patients and 20 PSPT patients) were used to quantify the tumor shrinkage in terms of the change of the primary gross tumor volume (GTVp). All patients were randomized to receive either IMRT or PSPT for NSCLC. Treatment planning goals were identical for both groups. All patients received 5 to 8 weekly repeated 4-dimensional computed tomography (4DCT) scans during the course of radiation treatments. The original GTVp contours were propagated to T50 of weekly 4DCT images using deformable image registration and their absolute volumes were measured. Statistical analysis was performed to compare the distribution of tumor shrinkage between the two population groups. In order to investigate the difference between the planned dose and the 5D dose with consideration of both breathing motion and anatomical change, we re-calculated new dose distributions at every phase of the breathing cycle for all available weekly 4DCT data sets which resulted 50 to 80 individual dose calculations for each of the 7 patients presented in this thesis. The newly calculated dose distributions were then deformed and accumulated to T50 of the planning 4DCT for comparison with the planned dose distribution. Results: At the end of the treatment, both IMRT and PSPT groups showed mean tumor volume reductions of 23.6% ( 19.2%) and 20.9% ( 17.0 %) respectively. Moreover, the mean difference in tumor shrinkage between two groups is 3% along with the corresponding 95% confidence interval, [-8%, 14%]. The rate of tumor shrinkage was highly correlated with the initial tumor volume size. For the planning dose and 5D dose comparison study, all 7 patients showed a mean difference of 1 % in terms of target coverage for both IMRT and PSPT treatment plans. Conclusions: The results of the tumor shrinkage investigation showed no statistically significant difference in tumor shrinkage between the IMRT and PSPT patients, and the tumor shrinkage between the two modalities is similar based on the 95% confidence interval. From the pilot study of comparing the planned dose with the 5D dose, we found the difference to be only 1%. Overall impression of the two modalities in terms of treatment response as measured by the tumor shrinkage and 5D dose under the influence of anatomical change that were designed under the same protocol (i.e. randomized trial) showed similar result.

Induction of synthetic lethality in mutant KRAS cells for non-small cell lung cancers chemoprevention and therapy

Lung cancer is the leading cause of cancer death in both men and women in the United States and worldwide. Despite improvement in treatment strategies, the 5-year survival rate of lung cancer patients remains low. Thus, effective chemoprevention and treatment approaches are sorely needed. Mutations and activation of KRAS occur frequently in tobacco users and the early stage of development of non-small cell lung cancers (NSCLC). So they are thought to be the primary driver for lung carcinogenesis. My work showed that KRAS mutations and activations modulated the expression of TNF-related apoptosis-inducing ligand (TRAIL) receptors by up-regulating death receptors and down-regulating decoy receptors. In addition, we showed that KRAS suppresses cellular FADD-like IL-1β-converting enzyme (FLICE)-like inhibitory protein (c-FLIP) expression through activation of ERK/MAPK-mediated activation of c-MYC which means the mutant KRAS cells could be specifically targeted via TRAIL induced apoptosis. The expression level of Inhibitors of Apoptosis Proteins (IAPs) in mutant KRAS cells is usually high which could be overcome by the second mitochondria-derived activator of caspases (Smac) mimetic. So the combination of TRAIL and Smac mimetic induced the synthetic lethal reaction specifically in the mutant-KRAS cells but not in normal lung cells and wild-type KRAS lung cancer cells. Therefore, a synthetic lethal interaction among TRAIL, Smac mimetic and KRAS mutations could be used as an approach for chemoprevention and treatment of NSCLC with KRAS mutations. Further data in animal experiments showed that short-term, intermittent treatment with TRAIL and Smac mimetic induced apoptosis in mutant KRAS cells and reduced tumor burden in a KRAS-induced pre-malignancy model and mutant KRAS NSCLC xenograft models. These results show the great potential benefit of a selective therapeutic approach for the chemoprevention and treatment of NSCLC with KRAS mutations.

Activation of matrix metalloproteinase -9 (MMP-9) by heregulin in human breast cancer cells involves multiple signaling pathways

Matrix metalloproteinase-9 (MMP-9) plays an important role in tumor invasion and angiogenesis. Secretion of MMP-9 has been reported in various cancer types including lung cancer, brain cancer, colon cancer, and breast cancer. Heregulin is a growth factor that regulates growth and differentiation of normal breast cells as well as mammary tumor cells. To study the role of heregulin in breast cancer metastasis, we tested whether heregulin may regulate MMP-9 secretion. By screening a panel of breast cancer cell line for their ability to respond to heregulin and produce MMP-9, we have found that MMP-9 secretion can be induced by heregulin-β1 in two breast cancer cell lines, SKBr3 and MCF-7. In both cell lines, increase of MMP-9 activity as shown by zymography was accompanied by increased protein level as well as mRNA level of MMP-9. Using a reporter luciferase assay, we have identified that proximal −670bp promoter of MMP-9 had similar activity to a 2.2kb MMP-9 promoter in response to heregulin stimulation. Heregulin treatment of SKBr3 and MCF-7 activated multiple signaling pathways inside cells. These include the Erk pathway, the p38 kinase pathway, PKC pathway, and PI-3K pathway. To examine which pathways are involved in MMP-9 activation by heregulin, we have used a panel of chemical inhibitors to specifically inhibit each one of these pathways. Ro-31-8220 (PKC inhibitor) and SB203580 (p38 kinase inhibitor) completely blocked heregulin activation of MMP-9. On the other hand, PD098059 (MEK-1 inhibitor) partially blocked MMP-9 activation, whereas PI-3K inhibitor, wortmannin, had no effect. Therefore, at least three signaling pathways are involved in activation of MMP-9 by heregulin. Since MMP-9 is tightly associated with metastatic potential, our study also suggests that heregulin may enhance breast tumor metastasis through induction of MMP-9 expression. ^

Thoracic radiotherapy treatment planning with cine PET/CT

Purpose: Respiratory motion causes substantial uncertainty in radiotherapy treatment planning. Four-dimensional computed tomography (4D-CT) is a useful tool to image tumor motion during normal respiration. Treatment margins can be reduced by targeting the motion path of the tumor. The expense and complexity of 4D-CT, however, may be cost-prohibitive at some facilities. We developed an image processing technique to produce images from cine CT that contain significant motion information without 4D-CT. The purpose of this work was to compare cine CT and 4D-CT for the purposes of target delineation and dose calculation, and to explore the role of PET in target delineation of lung cancer. Methods: To determine whether cine CT could substitute 4D-CT for small mobile lung tumors, we compared target volumes delineated by a physician on cine CT and 4D-CT for 27 tumors with intrafractional motion greater than 1 cm. We assessed dose calculation by comparing dose distributions calculated on respiratory-averaged cine CT and respiratory-averaged 4D-CT using the gamma index. A threshold-based PET segmentation model of size, motion, and source-to-background was developed from phantom scans and validated with 24 lung tumors. Finally, feasibility of integrating cine CT and PET for contouring was assessed on a small group of larger tumors. Results: Cine CT to 4D-CT target volume ratios were (1.05±0.14) and (0.97±0.13) for high-contrast and low-contrast tumors respectively which was within intraobserver variation. Dose distributions on cine CT produced good agreement (< 2%/1 mm) with 4D-CT for 71 of 73 patients. The segmentation model fit the phantom data with R2 = 0.96 and produced PET target volumes that matched CT better than 6 published methods (-5.15%). Application of the model to more complex tumors produced mixed results and further research is necessary to adequately integrate PET and cine CT for delineation. Conclusions: Cine CT can be used for target delineation of small mobile lesions with minimal differences to 4D-CT. PET, utilizing the segmentation model, can provide additional contrast. Additional research is required to assess the efficacy of complex tumor delineation with cine CT and PET. Respiratory-averaged cine CT can substitute respiratory-averaged 4D-CT for dose calculation with negligible differences.

THE ROLE OF TYROSINE PHOSPHORYLATION IN THE FUNCTIONS OF THE TUMOR SUPPRESSOR GPRC5A

The retinoic acid inducible G protein coupled receptor family C group 5 type A (GPRC5A) is expressed preferentially in normal lung tissue but its expression is suppressed in the majority of human non-small cell lung cancer cell lines and tissues. This differential expression has led to the idea that GPRC5A is a potential tumor suppressor. This notion was supported by the finding that mice with a deletion of the Gprc5a gene develop spontaneous lung tumors. However, there are various tumor cell lines and tissue samples, including lung, that exhibit higher GPRC5A expression than normal tissues and some reports by other groups that GPRC5A transfection increased cell growth and colony formation. Obviously, GPRC5A has failed to suppress the development of the tumors and the growth of the cell lines where its expression is not suppressed. Since no mutations were detected in the coding sequence of GPRC5A in 20 NSCLC cell lines, it’s possible that GPRC5A acts as a tumor suppressor in the context of some cells but not in others. Alternatively, we raised the hypothesis that the GPRC5A protein may be inactivated by posttranslational modification(s) such as phosphorylation. It is well established that Serine/Threonine phosphorylation of G protein coupled receptors leads to their desensitization and in a few cases Tyrosine phosphorylation of GPCRs has been linked to internalization. Others reported that GPRC5A can undergo tyrosine phosphorylation in the cytoplasmic domain after treatment of normal human mammary epithelial cells (HMECs) with epidermal growth factor (EGF) or Heregulin. This suggested that GPRC5A is a substrate of EGFR. Therefore, we hypothesized that tyrosine phosphorylation of GPRC5A by activation of EGFR signaling may lead to its inactivation. To test this hypothesis, we transfected human embryo kidney (HEK) 293 cells with GPRC5A and EGFR expression vectors and confirmed that GPRC5A can be tyrosine phosphorylated after activation of EGFR by EGF. Further, we found that EGFR and GPRC5A can interact either directly or through other proteins and that inhibition of the EGFR kinase activity decreased the phosphorylation of GPRA5A and the interaction between GPRC5A and EGFR. In c-terminal of GPRC5A, There are four tyrosine residues Y317, Y320, Y347, Y350. We prepared GPRC5A mutants in which all four tyrosine residues had been replaced by phenylalanine (mutant 4F) or each individual Tyr residue was replaced by Phe and found that Y317 is the major site for EGFR mediated phosphorylation in the HEK293T cell line. We also found that EGF can induce GPRC5A internalization both in H1792 transient and stable cell lines. EGF also partially inactivates the suppressive function of GPRC5A on cell invasion activity and anchorage-independent growth ability of H1792 stable cell lines. These finding support our hypothesis that GPRC5A may be inactivated by posttranslational modification- tyrosine phosphorylation.

MODIFICATION AND IMPLEMENTATION OF THE RPC HETEROGENEOUS THORAX PHANTOM FOR VERIFICATION OF PROTON THERAPY TREATMENT PROCEDURES

The Radiological Physics Center (RPC) provides heterogeneous phantoms that are used to evaluate radiation treatment procedures as part of a comprehensive quality assurance program for institutions participating in clinical trials. It was hypothesized that the existing RPC heterogeneous thorax phantom can be modified to assess lung tumor proton beam therapy procedures involving patient simulation, treatment planning, and treatment delivery, and could confirm agreement between the measured dose and calculated dose within 5%/3mm with a reproducibility of 5%. The Hounsfield Units (HU) for lung equivalent materials (balsa wood and cork) was measured using a CT scanner. The relative linear stopping power (RLSP) of these materials was measured. The linear energy transfer (LET) of Gafchromic EBT2 film was analyzed utilizing parallel and perpendicular orientations in a water tank and compared to ion chamber readings. Both parallel and perpendicular orientations displayed a quenching effect underperforming the ion chamber, with the parallel orientation showing an average 31 % difference and the perpendicular showing an average of 15% difference. Two treatment plans were created that delivered the prescribed dose to the target volume, while achieving low entrance doses. Both treatment plans were designed using smeared compensators and expanded apertures, as would be utilized for a patient in the clinic. Plan 1a contained two beams that were set to orthogonal angles and a zero degree couch kick. Plan 1b utilized two beams set to 10 and 80 degrees with a 15 degree couch kick. EBT2 film and TLD were inserted and the phantom was irradiated 3 times for each plan. Both plans passed the criteria for the TLD measurements where the TLD values were within 7% of the dose calculated by Eclipse. Utilizing the 5%/3mm criteria, the 3 trial average of overall pass rate was 71% for Plan 1a. The 3 trial average for the overall pass rate was 76% for Plan 1b. The trials were then analyzed using RPC conventional lung treatment guidelines set forth by the RTOG: 5%/5mm, and an overall pass rate of 85%. Utilizing these criteria, only Plan 1b passed for all 3 trials, with an average overall pass rate of 89%.

Evaluation of deformable image registration for improved 4D CT-derived ventilation for image guided radiotherapy

Recent treatment planning studies have demonstrated the use of physiologic images in radiation therapy treatment planning to identify regions for functional avoidance. This image-guided radiotherapy (IGRT) strategy may reduce the injury and/or functional loss following thoracic radiotherapy. 4D computed tomography (CT), developed for radiotherapy treatment planning, is a relatively new imaging technique that allows the acquisition of a time-varying sequence of 3D CT images of the patient's lungs through the respiratory cycle. Guerrero et al. developed a method to calculate ventilation imaging from 4D CT, which is potentially better suited and more broadly available for IGRT than the current standard imaging methods. The key to extracting function information from 4D CT is the construction of a volumetric deformation field that accurately tracks the motion of the patient's lungs during the respiratory cycle. The spatial accuracy of the displacement field directly impacts the ventilation images; higher spatial registration accuracy will result in less ventilation image artifacts and physiologic inaccuracies. Presently, a consistent methodology for spatial accuracy evaluation of the DIR transformation is lacking. Evaluation of the 4D CT-derived ventilation images will be performed to assess correlation with global measurements of lung ventilation, as well as regional correlation of the distribution of ventilation with the current clinical standard SPECT. This requires a novel framework for both the detailed assessment of an image registration algorithm's performance characteristics as well as quality assurance for spatial accuracy assessment in routine application. Finally, we hypothesize that hypo-ventilated regions, identified on 4D CT ventilation images, will correlate with hypo-perfused regions in lung cancer patients who have obstructive lesions. A prospective imaging trial of patients with locally advanced non-small-cell lung cancer will allow this hypothesis to be tested. These advances are intended to contribute to the validation and clinical implementation of CT-based ventilation imaging in prospective clinical trials, in which the impact of this imaging method on patient outcomes may be tested.

Improving the accuracy of radiation pneumonitis dose response models

The prognosis for lung cancer patients remains poor. Five year survival rates have been reported to be 15%. Studies have shown that dose escalation to the tumor can lead to better local control and subsequently better overall survival. However, dose to lung tumor is limited by normal tissue toxicity. The most prevalent thoracic toxicity is radiation pneumonitis. In order to determine a safe dose that can be delivered to the healthy lung, researchers have turned to mathematical models predicting the rate of radiation pneumonitis. However, these models rely on simple metrics based on the dose-volume histogram and are not yet accurate enough to be used for dose escalation trials. The purpose of this work was to improve the fit of predictive risk models for radiation pneumonitis and to show the dosimetric benefit of using the models to guide patient treatment planning. The study was divided into 3 specific aims. The first two specifics aims were focused on improving the fit of the predictive model. In Specific Aim 1 we incorporated information about the spatial location of the lung dose distribution into a predictive model. In Specific Aim 2 we incorporated ventilation-based functional information into a predictive pneumonitis model. In the third specific aim a proof of principle virtual simulation was performed where a model-determined limit was used to scale the prescription dose. The data showed that for our patient cohort, the fit of the model to the data was not improved by incorporating spatial information. Although we were not able to achieve a significant improvement in model fit using pre-treatment ventilation, we show some promising results indicating that ventilation imaging can provide useful information about lung function in lung cancer patients. The virtual simulation trial demonstrated that using a personalized lung dose limit derived from a predictive model will result in a different prescription than what was achieved with the clinically used plan; thus demonstrating the utility of a normal tissue toxicity model in personalizing the prescription dose.

Development of a Bayesian Joint Logistic Model to Better Study the Association between Haplotypes and Disease

In 2011, there will be an estimated 1,596,670 new cancer cases and 571,950 cancer-related deaths in the US. With the ever-increasing applications of cancer genetics in epidemiology, there is great potential to identify genetic risk factors that would help identify individuals with increased genetic susceptibility to cancer, which could be used to develop interventions or targeted therapies that could hopefully reduce cancer risk and mortality. In this dissertation, I propose to develop a new statistical method to evaluate the role of haplotypes in cancer susceptibility and development. This model will be flexible enough to handle not only haplotypes of any size, but also a variety of covariates. I will then apply this method to three cancer-related data sets (Hodgkin Disease, Glioma, and Lung Cancer). I hypothesize that there is substantial improvement in the estimation of association between haplotypes and disease, with the use of a Bayesian mathematical method to infer haplotypes that uses prior information from known genetics sources. Analysis based on haplotypes using information from publically available genetic sources generally show increased odds ratios and smaller p-values in both the Hodgkin, Glioma, and Lung data sets. For instance, the Bayesian Joint Logistic Model (BJLM) inferred haplotype TC had a substantially higher estimated effect size (OR=12.16, 95% CI = 2.47-90.1 vs. 9.24, 95% CI = 1.81-47.2) and more significant p-value (0.00044 vs. 0.008) for Hodgkin Disease compared to a traditional logistic regression approach. Also, the effect sizes of haplotypes modeled with recessive genetic effects were higher (and had more significant p-values) when analyzed with the BJLM. Full genetic models with haplotype information developed with the BJLM resulted in significantly higher discriminatory power and a significantly higher Net Reclassification Index compared to those developed with haplo.stats for lung cancer. Future analysis for this work could be to incorporate the 1000 Genomes project, which offers a larger selection of SNPs can be incorporated into the information from known genetic sources as well. Other future analysis include testing non-binary outcomes, like the levels of biomarkers that are present in lung cancer (NNK), and extending this analysis to full GWAS studies.

Targeting of {\it HER\/}-{\it 2\/}/{\it neu\/} with E1A gene therapy

Amplification or overexpression of HER-2/neu has been demonstrated in human cancers of the ovary, breast, lung and correlated with chemoresistance and poor clinic prognosis. We have previously found that the adenovirus type 5 early region 1A (E1A) gene product can repress the overexpression and suppress the tumorigenic potential of HER-2/neu-overexpressing cancer cells. In addition, E1A has been reported to induce apoptosis and inhibit the metastatic potential of tumor cells. Therefore, E1A could be considered as a tumor suppressor gene in HER-2/neu-overexpressing cancer cells. To develop an efficient HER-2/neu-targeting gene therapy with E1A, adenoviral vector or cationic liposome was used to introduce E1A into human ovarian, breast and lung cancer cells. Successful therapeutic effects were achieved.^ A replication-deficient adenovirus containing the E1A gene, Ad.E1A(+), was used to infect HER-2/neu-overexpressing human ovarian cancer cell line. Ovarian cancer growth in vitro and colony formation in soft agarose were greatly inhibited.^ To examine tumor suppressor function of E1A in breast cancer, we introduced E1A in vitro by adenovirus into both HER-2/neu-overexpressing and low-expressing human breast cancer cell lines. In HER-2/neu-overexpressing cells, E1A greatly inhibited tumor cell growth in vitro and colony formation in soft agarose. However, in low HER-2/neu expressing cancer cell lines, E1A could only reduce colony formation in soft agarose but had no significant effect on cell growth in monolayer, indicating different effects of E1A in these two types of cancer cells. To test the local therapeutic efficacy of E1A, we used either adenovirus- or liposome-mediated E1A gene delivery systems in an orthotopic breast cancer animal model.^ To test the therapeutic efficacy of systemically-delivered E1A in vivo lung cancer, we treated mice bearing intratracheal lung cancer by i.v. tail injections of Ad.E1A(+). As a result, Ad.E1A(+) suppressed HER-2/neu overexpression and inhibited intratracheal lung cancer growth. However, no significant tumor suppression effect of Ad.E1A(+) was observed in mice bearing HER-2/neu low expressing cell line when the same therapeutic procedure was followed. (Abstract shortened by UMI.) ^