TP53 as a Biomarker in Head and Neck Squamous Cell Carcinoma

Currently, there are no molecular biomarkers that guide treatment decisions for patients with head and neck squamous cell carcinoma (HNSCC). Several retrospective studies have evaluated TP53 in HNSCC, and results have suggested that specific mutations are associated with poor outcome. However, there exists heterogeneity among these studies in the site and stage of disease of the patients reviewed, the treatments rendered, and methods of evaluating TP53 mutation. Thus, it remains unclear as to which patients and in which clinical settings TP53 mutation is most useful in predicting treatment failure. In the current study, we reviewed the records of a cohort of patients with advanced, resectable HNSCC who received surgery and post-operative radiation (PORT) and had DNA isolated from fresh tumor tissue obtained at the time of surgery. TP53 mutations were identified using Sanger sequencing of exons 2-11 and the associated splice regions of the TP53 gene. We have found that the group of patients with either non-disruptive or disruptive TP53 mutations had decreased overall survival, disease-free survival, and an increased rate of distant metastasis. When examined as an independent factor, disruptive mutation was strongly associated with the development of distant metastasis. As a second aim of this project, we performed a pilot study examining the utility of the AmpliChip® p53 test as a practical method for TP53 sequencing in the clinical setting. AmpliChip® testing and Sanger sequencing was performed on a separate cohort of patients with HNSCC. Our study demonstrated the ablity of the AmpliChip® to call TP53 mutation from a single formalin-fixed paraffin-embedded slide. The results from AmpliChip® testing were identical with the Sanger method in 11 of 19 cases, with a higher rate of mutation calls using the AmpliChip® test. TP53 mutation is a potential prognostic biomarker among patients with advanced, resectable HNSCC treated with surgery and PORT. Whether this subgroup of patients could benefit from the addition of concurrent or induction chemotherapy remains to be evaluated in prospective clinical trials. Our pilot study of the p53 AmpliChip® suggests this could be a practical and reliable method of TP53 analysis in the clinical setting.

Glutamate receptors expressed in {\it Xenopus\/} oocytes: Studies in function and regulation

The amino acid glutamate is the primary excitatory neurotransmitter for the CNS and is responsible for the majority of fast synaptic transmission. Glutamate receptors have been shown to be involved in multiple forms of synaptic plasticity such as LTP, LTD, and the formation of specific synaptic connections during development. In addition to contributing to the plasticity of the CNS, glutamate receptors also are involved in, at least in part, various pathological conditions such as epilepsy, ischemic damage due to stroke, and Huntington's chorea. The regulation of glutamate receptors, particularly the ionotropic NMDA and AMPA/KA receptors is therefore of great interest. In this body of work, glutamate receptor function and regulation by kinase activity was examined using the Xenopus oocyte which is a convenient and faithful expression system for exogenous proteins. Glutamate receptor responses were measured using the two-electrode voltage clamp technique in oocytes injected with rat total forebrain RNA. NMDA elicited currents that were glycine-dependent, subject to block by Mg$\sp{2+}$ in a voltage-dependent manner and sensitive to the specific NMDA antagonist APV in a manner consistent with those types of responses found in neural tissue. Similarly, KA-evoked currents were sensitive to the specific AMPA/KA antagonist CNQX and exhibited current voltage relationships consistent with the calcium permeable type II KA receptors found in the hippocampus. There is evidence to indicate that NMDA and AMPA/KA receptors are regulated by protein kinase A (PKA). We explored this by examining the effects of activators of PKA (forskolin, 1-isobutyl-3-methylxanthine (IBMX) and 8-Br-cAMP) on NMDA and KA currents in the oocyte. In buffer where Ca$\sp{2+}$ was replaced by 2 mM Ba$\sp{2+},$ forskolin plus IBMX and 8-Br-cAMP augmented currents due to NMDA application but not KA. This augmentation was abolished by pretreating the oocytes in the kinase inhibitor K252A. The use of chloride channel blockers resulted in attenuation of this effect indicating that Ba$\sp{2+}$ influx through the NMDA channel was activating the endogenous calcium-activated chloride current and that the cAMP mediated augmentation was at the level of the chloride channel and not the NMDA channel. This was confirmed by (1) the finding that 8-Br-cAMP increased chloride currents elicited via calcium channel activation while having no effect on the calcium channels themselves and (2) the fact that lowering the Ba$\sp{2+}$ concentration to 200 $\mu$M abolished the augmentation NMDA currents by 8-Br-cAMP. Thus PKA does not appear to modulate ionotropic glutamate receptors in our preparation. Another kinase also implicated in the regulation of NMDA receptors, calcium/phospholipid-dependent protein kinase (PKC), was examined for its effects on the NMDA receptor under low Ba$\sp{2+}$ (200 $\mu$M) conditions. Phorbol esters, activators of PKC, induced a robust potentiation of NMDA currents that was blockable by the kinase inhibitor K252A. Furthermore activation of metabotropic receptors by the selective agonist trans-ACPD, also potentiated NMDA albeit more modestly. These results indicate that neither NMDA nor KA-activated glutamate receptors are modulated by PKA in Xenopus oocytes whereas NMDA receptors appear to be augmented by PKC. Furthermore, the endogenous chloride current of the oocyte was found to be responsive to Ba$\sp{2+}$ and in addition is enhanced by PKA. Both of these latter findings are novel. In conclusion, the Xenopus oocyte is a useful expression system for the analysis of ligand-gated channel activity and the regulation of those channels by phosphorylation. ^

Conditional expression of Sry, and Wnt4 by gene-targeting: Studies in sexual and skeletal development

Sry and Wnt4 cDNAs were individually introduced into the ubiquitously-expressed Rosa26 ( R26) locus by gene targeting in embryonic stem (ES) cells to create a conditional gene expression system in mice. In the targeted alleles, expression of these cDNAs should be blocked by a neomycin resistance selection cassette that is flanked by loxP sites. Transgene expression should be activated after the blocking cassette is deleted by Cre recombinase. ^ To test this conditional expression system, I have bred R26-stop- Sry and R26-stop-Wnt4 heterozygotes with a MisRII-Cre mouse line that expresses Cre in the gonads of both sexes. Analysis of these two types of bigenic heterozygotes indicated that their gonads developed normally like those of wild types. However, one XX R26-Sry/R26-Sry; MisR2-Cre/+ showed epididymis-like structures resembling those of males. In contrast, only normal phenotypes were observed in XY R26-Wnt4/R26-Wnt4; MisR2-Cre /+ mice. To interpret these results, I have tested for Cre recombinase activity by Southern blot and transcription of the Sry and Wnt4 transgenes by RT-PCR. Results showed that bigenic mutants had insufficient activation of the transgenes in their gonads at E12.5 and E13.5. Therefore, the failure to observe mutant phenotypes may have resulted from low activity of MisR2-Cre recombination at the appropriate time. ^ Col2a1-Cre transgenic mice express Cre in differentiating chondrocytes. R26-Wnt4; Col2a1-Cre bigenic heterozygous mice were found to exhibit a dramatic alteration in growth presumably caused by Wnt4 overexpression during chondrogenesis. R26-Wnt4; Col2a1-Cre mice exhibited dwarfism beginning approximately 10 days after birth. In addition, they also had craniofacial abnormalities, and had delayed ossification of the lumbar vertebrate and pelvic bones. Histological analysis of the growth plates of R26-Wnt4; Col2a1-Cre mice revealed less structural organization and a delay in onset of the primary and secondary ossification centers. Molecular studies confirmed that overexpression of Wnt4 causes decreased proliferation and early maturation of chondrocytes. In addition, R26-Wnt4; Col2a1-Cre mice had decreased expression of vascular endothelial growth factor (VEGF), suggesting that defects in vascularization may contribute to the dwarf phenotype. Finally, 9-month-old R26-Wnt4; Col2a1-Cre mice had significantly more fat cells in the marrow cavities of their metaphysis long bones, implying that long-term overexpression of Wnt4may cause bone marrow pathologies. In conclusion, Wnt4 was activated by Col2a1-Cre recombinase and was overexpressed in the growth plate, resulting in aberrant proliferation and differentiation of chondrocytes, and ultimately leads to dwarfism in mice. ^

Comparative studies of Hox genes in mammalian limb and vertebral pattern formation

Divergence of anterior-posterior (AP) limb pattern and differences in vertebral column morphology are the two main examples of mammalian evolution. The Hox genes (homeobox containing gene) have been implicated in driving evolution of these structures. However, regarding Hox genes, how they contribute to the generation of mammalian morphological diversities, is still unclear. Implementing comparative gene expression and phenotypic rescue studies for different mammalian Hox genes could aid in unraveling this mystery. In the first part of this thesis, the expression pattern of Hoxd13 gene, a key Hox gene in the establishment of the limb AP pattern, was examined in developing limbs of bats and mice. Bat forelimbs exhibit a pronounced asymmetric AP pattern and offer a good model to study the molecular mechanisms that contribute to the variety of mammalian limbs. The data showed that the expression domain of bat Hoxd13 was shifted prior to the asymmetric limb plate expansion, whereas its domain in mice was much more symmetric. This finding reveals a correlation between the divergence of Hoxd13 expression and the AP patterning difference in limb development. The second part of this thesis details a phenotypic rescue approach by human HOXB1-9 transgenes in mice with Hoxb1-9 deletion, The mouse mutants displayed homeosis in cervical and anterior thoracic vertebrae. The human transgenes entirely rescued the mouse mutants, suggesting that these human HOX genes have similar functions to their mouse orthologues in anterior axial skeletal patterning. The anterior expressing human HOXB transgenes such as HOXB1-3 were expressed in the mouse embryonic trunk in a similar manner as their murine orthologues. However, the anterior boundary of human HOXB9 expression domain was more posterior than that of the mouse Hoxb9 by 2-3 somites. These data provide the molecular support for the hypothesis that Hox genes are responsible for maintaining similar anterior axial skeletal architectures cervical and anterior thoracic regions, but different architectures in lumbar and posterior thoracic regions between humans and mice. ^

STUDIES ON CHOLINERGIC FUNCTION AND MEMORY IN MICE (CHOLINOTOXINS, WORKING, REFERENCE MEMORY, RADIAL ARM MAZE)

Considerable evidence suggests that central cholinergic neurons participate in either acquisition, storage or retrieval of information. Experiments were designed to evaluate information processing in mice following either reversible or irreversible impairment in central cholinergic activity. The cholinergic receptor antagonists, atropine and methylatropine were used to reversibly inhibit cholinergic transmission. Irreversible impairment in central cholinergic function was achieved by central administration of the cholinergic-specific neurotoxins, N-ethyl-choline aziridinium (ECA) and N-ethyl-acetylcholine aziridinium (EACA).^ ECA and EACA appear to act by irreversible inhibition of high affinity choline uptake (proposed rate-limiting step in acetylcholine synthesis). Intraventricular administration of ECA or EACA produced persistent reduction in hippocampal choline acetyltransferase activity. Other neuronal systems and brain regions showed no evidence of toxicity.^ Mice treated with either ECA or EACA showed behavioral deficits associated with cholinergic dysfunction. Passive avoidance behavior was significantly impaired by cholinotoxin treatment. Radial arm maze performance was also significantly impaired in cholinotoxin-treated animals. Deficits in radial arm maze performance were transient, however, such that rapid and apparent complete behavioral recovery was seen during retention testing. The centrally active cholinergic receptor antagonist atropine also caused significant impairment in radial arm maze behavior, while equivalent doses of methylatropine were without effect.^ The relative effects of cholinotoxin and receptor antagonist treatment on short-term (working) memory and long-term (reference) memory in radial arm maze behavior were examined. Maze rotation studies indicated that there were at least two different response strategies which could result in accurate maze performance. One strategy involved the use of response algorithms and was considered to be a function of reference memory. Another strategy appeared to be primarily dependent on spatial working memory. However, all behavioral paradigms with multiple trails have reference memory requirements (i.e. information useful over all trials). Performance was similarly affected following either cholinotoxin or anticholinergic treatment, regardless of the response strategy utilized. In addition, rates of behavioral recovery following cholinotoxin treatment were similar between response groups. It was concluded that both cholinotoxin and anticholinergic treatment primarily resulted in impaired reference memory processes. ^

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Functional studies of *Ras and Bcl-2 oncoproteins in keratinocyte homeostasis and multistep skin carcinogenesis

To assess the effect of deregulated Ha-ras and bcl-2, individually and in combination on epidermal keratinocyte homeostasis and during multistep skin carcinogenesis, we generated skin-specific transgenic mice and keratinocyte transfectants constitutively expressing oncogenic Ha-ras and bcl-2 proteins. The deregulated Ha-ras and bcl-2 expression contributing to homeostatic imbalances in the skin had an additive effect on the probability of tumor development. They were also cooperative in incidence, growth, and latency of tumor formation, and they exhibited synergistic cooperation in malignant transformation of benign papillomas. To explain the homeostatic imbalances by Ha-ras and bcl-2 overexpression in the skin, we investigated the three major cellular processes of proliferation, cell death, and differentiation. Epidermal expression of Bcl-2 retarded keratinocyte proliferation in the epidermis of neonatal mice compared with results for control littermates. Constitutive expression of Ha-ras increased keratinocyte proliferation, and co-expression of bcl-2 modestly suppressed the ras-mediated abnormal proliferation of neonatal keratinocytes. Bcl-2 proteins in keratinocytes protected UV-treated cells from apoptotic cell death regardless of oncogenic ras expression in both non-neoplastic neonatal epidermis and human keratinocyte cell lines. The spontaneous apoptotic index (AI) was also lower in papillomas constitutively expressing bcl-2 compared with the ones that developed in control mice. Ras-overexpressing epidermis, including that in ras/bcl-2 double transgenic mice, had abnormal differentiation patterns compared with controls. The oncogenic ras protein had alterations in both epidermal distribution and the extent of cytokeratin 14 and involucrin expression. Abnormal expression of the hyperproliferation marker cytokeratin 6 and modest down regulation of cytokeratin 1 were also detected. Late appearance of filaggrin was another abnormal phenotype of the ras-expressing epidermis. Overexpression of bcl-2 had no effect on epidermal differentiation. Together, these findings suggest that constitutive expression of oncogenic Ha-ras and bcl-2 are important determinants of epidermal proliferation, viability and differentiation. In summary, our results demonstrated that the disruption of epidermal homeostasis by overexpressed ras and bcl-2 predisposes to hyperplastic growth of the epidermis and to papilloma development and that these proteins with distinct mechanisms for oncogenesis are functionally synergistic for malignant transformation of chemically induced skin carcinogenesis. ^