PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

IDENTIFICATION OF FACTORS INVOLVED IN DNA METHYLATION OF CPG-ISLAND-PROMOTERS

Repression of many tumor suppressor genes (TSGs) in cancer is mediated by aberrantly increased DNA methylation levels at promoter CpG islands (CGI). About one-fourth of empirically defined human promoters are surrounded by or contain clustered repetitive elements. It was previously observed that a sharp transition of methylation occurs between highly methylated repetitive elements (SINE or LINE) and unmethylated CGI-promoters (e.g. P16, VHL, CDH and RIL) in normal tissues. The functions that lead to increased CGI methylation in cancer remain poorly understood. We propose that CGI-promoters contain cis-elements for triggering de novo DNA methylation. In the first part of our project, we established a site-specific integration system with enforced local transcriptional repression in colorectal cancer cells and monitored the occurrence of de novo DNA methylation in exogenous fragments containing a CGI-promoter and repetitive elements. Initial de novo methylation was seeded at specific CG sites in a repetitive element, and accelerated by persistent binding of a KRAB-containing transcriptional repressor. Furthermore, additional repetitive elements (LINE and SINE) located adjacent to the promoter could confer DNA methylation spreading into the CGI particularly in the setting of KRAB-factor binding. However, a repressive chromatin alone was not sufficient to initiate DNA methylation, which required specific DNA sequences and was integration-site (and/or cell-line) specific. In addition, all the methylation observed showed slow and gradual accumulation over several months of culture. Overall, these results demonstrate a requirement for specific DNA sequences to trigger de novo DNA methylation, and repetitive elements as cis-regulatory factors to cooperate with strengthened transcriptional repression in promoting methylation spreading. In the second part, we re-introduced disrupted DNMT3B or DNMT1 into HCT116 DKO cells and mapped the remethylation pattern through a profiling method (DREAM). Moderate remethylation occurred when DNMT3B was re-expressed with a preference toward non-CGI and non-promoter regions. Hence, there exists a set of genomic regions with priority to be targets for DNMT3B in somatic cells.

The effect of DNA methylation in carcinogenesis

CpG island methylation within single gene promoters can silence expression of associated genes. We first extended these studies to bidirectional gene pairs controlled by single promoters. We showed that hypermethylation of bidirectional promoter-associated CpG island silences gene pairs (WNT9A/CD558500, CTDSPL/BC040563, and KCNK15/BF 195580) simultaneously. Hypomethylation of these promoters by 5-aza-2'-deoxycytidine treatment reactivated or enhanced gene expression bidirectionally. These results were further confirmed by luciferase assays. Methylation of WNT9A/CD558500 and CTDSPL/BC040563 promoters occurs frequently in primary colon cancers and acute lymphoid leukemia, respectively. ^ Next we sought to understand the origins of hypermethylation in cancer. CpG islands associated with tumor suppressor genes are normally free from methylation, but can be hypermethylated in cancer. It remains poorly understood how these genes are protected from methylation in normal tissues. In our studies, we aimed to determine if cis-acting elements in these genes are responsible for this protection, using the tumor suppressor gene p16 as a model. We found that Alu repeats located both upstream and downstream of the p16 promoter become hypermethylated with age. In colon cancer samples, the methylation level is particularly high, and the promoter can also be affected. Therefore, the protection in the promoter against methylation spreading could fail during tumorigenesis. This methylation pattern in p16 was also observed in cell lines of different tissue origins, and their methylation levels were found to be inversely correlated with that of active histone modification markers (H3K4-3me and H3K9-Ac). To identify the mechanism of protection against methylation spreading, we constructed serial deletions of the p16 protected region and used silencing of a neomycin reporter gene to evaluate the protective effects of these fragments. A 126 bp element was identified within the region which exerts bidirectional protection against DNA methylation, independently of its transcriptional activity. The protective strength of this element is comparable to that of the HS4 insulator. During long-term culture, the presence of this element significantly slowed methylation spreading. In conclusion, we have found that an element located in the p16 promoter is responsible for protection against DNA methylation spreading in normal tissues. The failure of protective cis-elements may be a general feature of tumor-suppressor gene silencing during tumorigenesis. ^

DNA methylation inhibits p53 mediated survivin repression in cancer cells

Survivin (BIRC5) is a member of the Inhibitor of Apoptosis (IAP) gene family and functions as a chromosomal passenger protein as well as a mediator of cell survival. Survivin is widely expressed during embryonic development then becomes transcriptionally silent in most highly differentiated adult tissues. It is also overexpressed in virtually every type of tumor. The survivin promoter contains a canonical CpG island that has been described as epigenetically regulated by DNA methylation. We observed that survivin is overexpressed in high grade, poorly differentiated endometrial tumors, and we hypothesized that DNA hypomethylation could explain this expression pattern. Surprisingly, methylation specific PCR and bisulfite pyrosequencing analysis showed that survivin was hypermethylated in endometrial tumors and that this hypermethylation correlated with increased survivin expression. We proposed that methylation could activate survivin expression by inhibit the binding of a transcriptional repressor. ^ The tumor suppressor protein p53 is a well documented transcriptional repressor of survivin and examination of the survivin promoter showed that the p53 binding site contains 3 CpG sites which often become methylated in endometrial tumors. To determine if methylation regulates survivin expression, we treated HCT116 cells with decitabine, a demethylation agent, and observed that survivin transcript and protein levels were significantly repressed following demethylation in a p53 dependent manner. Subsequent binding studies confirmed that DNA methylation inhibited the binding of p53 protein to its binding site in the survivin promoter. ^ We are the first to report this novel mechanism of epigenetic regulation of survivin. We also conducted microarray analysis which showed that many other cancer relevant genes may also be regulated in this manner. While demethylation agents are traditionally thought to inhibit cancer cell growth by reactivating tumor suppressors, our results indicate that an additional important mechanism is to decrease the expression of oncogenes. ^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

The effects of DNA cytosine methylation on H -ras promoter structure and function

The effect of DNA cytosine methylation on H-ras promoter activity was assessed using a transient expression system employing the plasmid H-rasCAT (NaeI H-ras promoter linked to the chloramphenicol acetyltransferase (CAT) gene). This 551 bp promoter is 80% GC rich, enriched with 168 CpG dinucleotides, and contains six functional GC box elements which represent major DNA methylation target sites. Prokaryotic methyltransferases HhaI (CGm$\sp5$CG) and HpaII (Cm$\sp5$CGG) alone or in combination with a human placental methyltransferase (HP MTase) were used to introduce methyl groups at different CpG sites within the promoter. To test for functional promoter activity, the methylated plasmids were introduced into CV-1 cells and CAT activity assessed 48 h post-transfection. Methylation at specific HhaI and HpaII sites reduced CAT expression by 70%, whereas more extensive methylation at generalized CpG sites with HP MTase inactivated the promoter $>$95%. The inhibition of H-ras promoter activity was not attributable to methylation-induced differences in DNA uptake or stability in the cell, topological form of the plasmid, or methylation effects in nonpromoter regions. We also observed that DNA cytosine methylation of a 360 bp promoter fragment by HP MTase induced a local change in DNA conformation. Using three independent methodologies (nitrocellulose filter binding assays, gel mobility shifts, and Southwestern blots), we determined that this change in promoter conformation affected the interaction of nuclear proteins with cis-regulatory sequences residing in the promoter region. The results provide evidence to suggest that DNA methylation may regulate gene expression by inducing changes in local promoter conformation which in turn alters the interactions between DNA and protein factors required for transcription. The results provide supportive evidence for the hypothesis of Cedar and Riggs, who postulated that DNA methylation may regulate gene expression by altering the binding affinities of proteins for DNA. ^

An Sp1/Sp3 binding polymorphism confers methylation protection.

Hundreds of genes show aberrant DNA hypermethylation in cancer, yet little is known about the causes of this hypermethylation. We identified RIL as a frequent methylation target in cancer. In search for factors that influence RIL hypermethylation, we found a 12-bp polymorphic sequence around its transcription start site that creates a long allele. Pyrosequencing of homozygous tumors revealed a 2.1-fold higher methylation for the short alleles (P<0.001). Bisulfite sequencing of cancers heterozygous for RIL showed that the short alleles are 3.1-fold more methylated than the long (P<0.001). The comparison of expression levels between unmethylated long and short EBV-transformed cell lines showed no difference in expression in vivo. Electrophorectic mobility shift assay showed that the inserted region of the long allele binds Sp1 and Sp3 transcription factors, a binding that is absent in the short allele. Transient transfection of RIL allele-specific transgenes showed no effects of the additional Sp1 site on transcription early on. However, stable transfection of methylation-seeded constructs showed gradually decreasing transcription levels from the short allele with eventual spreading of de novo methylation. In contrast, the long allele showed stable levels of expression over time as measured by luciferase and approximately 2-3-fold lower levels of methylation by bisulfite sequencing (P<0.001), suggesting that the polymorphic Sp1 site protects against time-dependent silencing. Our finding demonstrates that, in some genes, hypermethylation in cancer is dictated by protein-DNA interactions at the promoters and provides a novel mechanism by which genetic polymorphisms can influence an epigenetic state.

Evolutionary conservation of residues in vertebrate DNA polymerase N conferring low fidelity and bypass activity.

POLN is a nuclear A-family DNA polymerase encoded in vertebrate genomes. POLN has unusual fidelity and DNA lesion bypass properties, including strong strand displacement activity, low fidelity favoring incorporation of T for template G and accurate translesion synthesis past a 5S-thymine glycol (5S-Tg). We searched for conserved features of the polymerase domain that distinguish it from prokaryotic pol I-type DNA polymerases. A Lys residue (679 in human POLN) of particular interest was identified in the conserved 'O-helix' of motif 4 in the fingers sub-domain. The corresponding residue is one of the most important for controlling fidelity of prokaryotic pol I and is a nonpolar Ala or Thr in those enzymes. Kinetic measurements show that K679A or K679T POLN mutant DNA polymerases have full activity on nondamaged templates, but poorly incorporate T opposite template G and do not bypass 5S-Tg efficiently. We also found that a conserved Tyr residue in the same motif not only affects sensitivity to dideoxynucleotides, but also greatly influences enzyme activity, fidelity and bypass. Protein sequence alignment reveals that POLN has three specific insertions in the DNA polymerase domain. The results demonstrate that residues have been strictly retained during evolution that confer unique bypass and fidelity properties on POLN.

Loss of DNA polymerase zeta enhances spontaneous tumorigenesis.

Mammalian genomes encode at least 15 distinct DNA polymerases, functioning as specialists in DNA replication, DNA repair, recombination, or bypass of DNA damage. Although the DNA polymerase zeta (polzeta) catalytic subunit REV3L is important in defense against genotoxins, little is known of its biological function. This is because REV3L is essential during embryogenesis, unlike other translesion DNA polymerases. Outstanding questions include whether any adult cells are viable in the absence of polzeta and whether polzeta status influences tumorigenesis. REV3L-deficient cells have properties that could influence the development of neoplasia in opposing ways: markedly reduced damage-induced point mutagenesis and extensive chromosome instability. To answer these questions, Rev3L was conditionally deleted from tissues of adult mice using MMTV-Cre. Loss of REV3L was tolerated in epithelial tissues but not in the hematopoietic lineage. Thymic lymphomas in Tp53(-/-) Rev3L conditional mice occurred with decreased latency and higher incidence. The lymphomas were populated predominantly by Rev3L-null T cells, showing that loss of Rev3L can promote tumorigenesis. Remarkably, the tumors were frequently oligoclonal, consistent with accelerated genetic changes in the absence of Rev3L. Mammary tumors could also arise from Rev3L-deleted cells in both Tp53(+/+) and Tp53(+/-) backgrounds. Mammary tumors in Tp53(+/-) mice deleting Rev3L formed months earlier than mammary tumors in Tp53(+/-) control mice. Prominent preneoplastic changes in glandular tissue adjacent to these tumors occurred only in mice deleting Rev3L and were associated with increased tumor multiplicity. Polzeta is the only specialized DNA polymerase yet identified that inhibits spontaneous tumor development.

POLG DNA testing as an emerging standard of care before instituting valproic acid therapy for pediatric seizure disorders.

PURPOSE: To review our clinical experience and determine if there are appropriate signs and symptoms to consider POLG sequencing prior to valproic acid (VPA) dosing in patients with seizures. METHODS: Four patients who developed VPA-induced hepatotoxicity were examined for POLG sequence variations. A subsequent chart review was used to describe clinical course prior to and after VPA dosing. RESULTS: Four patients of multiple different ethnicities, age 3-18 years, developed VPA-induced hepatotoxicity. All were given VPA due to intractable partial seizures. Three of the patients had developed epilepsia partialis continua. The time from VPA exposure to liver failure was between 2 and 3 months. Liver failure was reversible in one patient. Molecular studies revealed homozygous p.R597W or p.A467T mutations in two patients. The other two patients showed compound heterozygous mutations, p.A467T/p.Q68X and p.L83P/p.G888S. Clinical findings and POLG mutations were diagnostic of Alpers-Huttenlocher syndrome. CONCLUSION: Our cases underscore several important findings: POLG mutations have been observed in every ethnic group studied to date; early predominance of epileptiform discharges over the occipital region is common in POLG-induced epilepsy; the EEG and MRI findings varying between patients and stages of the disease; and VPA dosing at any stage of Alpers-Huttenlocher syndrome can precipitate liver failure. Our data support an emerging proposal that POLG gene testing should be considered in any child or adolescent who presents or develops intractable seizures with or without status epilepticus or epilepsia partialis continua, particularly when there is a history of psychomotor regression.

CREB trans-activation of disruptor of telomeric silencing-1 mediates forskolin inhibition of CTGF transcription in mesangial cells.

Connective tissue growth factor (CTGF) participates in diverse fibrotic processes including glomerulosclerosis. The adenylyl cyclase agonist forskolin inhibits CTGF expression in mesangial cells by unclear mechanisms. We recently reported that the histone H3K79 methyltransferase disruptor of telomeric silencing-1 (Dot1) suppresses CTGF gene expression in collecting duct cells (J Clin Invest 117: 773-783, 2007) and HEK 293 cells (J Biol Chem In press). In the present study, we characterized the involvement of Dot1 in mediating the inhibitory effect of forskolin on CTGF transcription in mouse mesangial cells. Overexpression of Dot1 or treatment with forskolin dramatically suppressed basal CTGF mRNA levels and CTGF promoter-luciferase activity, while hypermethylating H3K79 in chromatin associated with the CTGF promoter. siRNA knockdown of Dot1 abrogated the inhibitory effect of forskolin on CTGF mRNA expression. Analysis of the Dot1 promoter sequence identified a CREB response element (CRE) at -384/-380. Overexpression of CREB enhanced forskolin-stimulated Dot1 promoter activity. A constitutively active CREB mutant (CREB-VP16) strongly induced Dot1 promoter-luciferase activity, whereas overexpression of CREBdLZ-VP16, which lacks the CREB DNA-binding domain, abolished this activation. Mutation of the -384/-380 CRE resulted in 70% lower levels of Dot1 promoter activity. ChIP assays confirmed CREB binding to the Dot1 promoter in chromatin. We conclude that forskolin stimulates CREB-mediated trans-activation of the Dot1 gene, which leads to hypermethylation of histone H3K79 at the CTGF promoter, and inhibition of CTGF transcription. These data are the first to describe regulation of the Dot1 gene, and disclose a complex network of genetic and epigenetic controls on CTGF transcription.

Genetic predictors of bladder cancer with an emphasis on methylation-related genes

Although tobacco exposure remains the prevailing risk factor for bladder cancer (BC), only a small percentage of exposed individuals develop cancer, suggesting that tobacco-related carcinogenesis is modulated by genetic susceptibility and possibly by DNA methylation-related events. Methylation patterns established by DNA methyltransferases (DNMTs) are influenced by dietary folate and genetic polymorphisms in the methylene-tetrahydrofolate reductase gene (MTHFR). Therefore, we hypothesized that DNA methylation-related genes, such as DNMT3B and MTHFR, might modulate BC risk. ^ In a study of 514 Caucasian BC cases and 498 healthy Caucasian controls examining the DNMT3B C46359T polymorphism, CC genotype was found to be a risk factor in women (Odds Ratio (OR) = 1.79), but not in men. This risk was further increased among women who were never smokers, consumed low dietary folate, and had adverse variants of MTHFR. In addition, higher DNMT3B expression among smokers was a risk factor (OR = 4.27) and correlated with genetic variants of the DNMT3B C46359T polymorphism, providing salient evidence for the risk associated with the CC variant. This suggests that the DNMT3B CC variant may confer a predisposition toward aberrant de novo methylation of CpG islands in critical tumor suppressor genes. ^ The convergence of alterations in DNMT3B, associated with promoter methylation, and reduced dietary folate consumption, accompanying global hypomethylation and genetic instability, may act synergistically to promote bladder carcinogenesis, especially in women. The results of this study unveiled new gender-specific paradigms of BC risk for women and demonstrated that this risk can be modified by folate consumption as well as polymorphisms in the folate pathway. ^

Mechanisms of transcription inhibition by 8-amino-adenosine

Nucleoside analogs are a class of chemotherapeutic agents with tremendous utility in treating viral infections and cancers. Traditional nucleoside analogs are DNA-directed. However, there is a new group of nucleoside analogs that induce cell death by a direct effect on RNA synthesis. The adenosine analog, 8-chloroadenosine, is incorporated into RNA and is currently in clinical trials. Another congener, 8-amino-adenosine has demonstrated toxicity in multiple myeloma cell lines. Like other nucleoside analogs, 8-amino-adenosine must be metabolized to its triphosphate to elicit a cytotoxic effect. Furthermore, 8-amino-adenosine causes a decline of the intracellular ATP pool and inhibits mRNA poly(A) adenylation. ^ Because of the previously known adenosine analog mechanism as well as the scope of the RNA directed nucleoside analog field, I hypothesized there are multiple mechanisms of transcription inhibition mediating 8-amino-adenosine-induced cell death. Prior to investigating these mechanisms, cell death by 8-amino-adenosine was characterized. 8-Amino-adenosine activates PARP cleavage and induces the caspase cascade. 8-Amino-adenosine increases Annexin V binding and the mitochondrial membrane permeability in wild-type MEF cells. In BAX/BAK deficient MEF cells, 8-amino-adenosine decreases the mitochondrial membrane permeability and induces autophagy. ^ Once cell death was characterized, the mechanisms of 8-amino-adenosine transcription inhibition were assessed. It was established that 8-aminoadenosine treatment causes 8-amino-ATP accumulation and decreases the intracellular ATP concentration, resulting in RNA synthesis inhibition. Several other mechanisms are identified. First, a relationship between ATP decline by 8-amino-adenosine or other known ATP synthesis inhibitors and RNA synthesis is established indicating that effects on cellular bioenergy, regardless of the mechanism of ATP decline, can decrease RNA synthesis. Second, 8-aminoadenosine treatment decreases the phosphorylation of serine residues on the RNA polymerase II C-terminal domain which regulates transcription initiation and elongation. Third, evidence is provided to demonstrate 8-amino-ATP is a substrate for RNA synthesis. Fourth, 8-amino-ATP is incorporated at the 3'-terminal position leading to chain termination. Finally, in vitro transcription assays show that 8-amino-ATP may compete with ATP to decrease de novo mRNA synthesis. Overall, this work demonstrates 8-amino-adenosine is a cytotoxic nucleoside analog that functions to inhibit RNA transcription through multiple mechanisms. ^

Mechanism of methylation gene reactivation in human cancer cells

Epigenetic silencing of tumor suppressor genes by DNA hypermethylation at promoter regions is a common event in carcinogenesis and tumor progression. Abrogation of methylation and reversal of epigenetic silencing is a very potent way in cancer treatment. However, the reactivation mechanisms are poorly understood. In this study, we first developed a cell line model system named YB5, derived from SW48 cancer cell line, which bears one copy of stably integrated EGFP gene on Chromosome 1p31.1 region. The GFP gene expression is transcriptionally silenced due to the hypermethylated promoter CMV. However, the GFP expression can be restored using demethylating agent 5-aza-2' deoxycytidine (DAC), and detected by FACS and fluorescent microscopy. Using this system, we observed the heterogeneous reactivation induced by DAC treatment. After flow sorting, GFP negative cells exhibited similar level of incomplete demethylation compared to GFP positive cells on repetitive LINE1 element, tumor suppressor genes such as P16, CDH13, and RASSF1a, and CMV promoter as well. However, the local chromatin of CMV-GFP locus altered to an open structure marked by high H3 lysine 9 acetylation and low H3 lysine 27 tri-methylation in GFP positive cells, while the GFP negative cells retained mostly the original repressive marks. Thus, we concluded that DAC induced DNA hypomethylation alone does not directly determine the level of re-expression, and the resetting of the local chromatin structure under hypomethylation environment is required for gene reactivation. Besides, a lentivirus vector-based shRNA screening was performed using the YB5 system. Although it is the rare chance that vector lands in the neighboring region of GFP, we found that the exogenous vector DNA inserted into the upstream region of GFP gene locus led to the promoter demethylation and reactivated the silenced GFP gene. Thus, epigenetic state can be affected by changing of the adjacent nucleic acid sequences. Further, this hypermethylation silenced system was utilized for epigenetic drug screening. We have found that DAC combined with carboplatin would enhance the GFP% yield and increase expression of other tumor suppressor genes than DAC alone, and this synergistic effect may be related to DNA repair process. In summary, these studies reveal that reversing of methylation silencing requires coordinated alterations of DNA methylation, chromatin structure, and local microenvironment. ^

A BAYESIAN APPROACH TO DOSE-RESPONSE ASSESSMENT AND DRUG-DRUG INTERACTION ANALYSIS: APPLICATION TO IN VITRO STUDIES

The considerable search for synergistic agents in cancer research is motivated by the therapeutic benefits achieved by combining anti-cancer agents. Synergistic agents make it possible to reduce dosage while maintaining or enhancing a desired effect. Other favorable outcomes of synergistic agents include reduction in toxicity and minimizing or delaying drug resistance. Dose-response assessment and drug-drug interaction analysis play an important part in the drug discovery process, however analysis are often poorly done. This dissertation is an effort to notably improve dose-response assessment and drug-drug interaction analysis. The most commonly used method in published analysis is the Median-Effect Principle/Combination Index method (Chou and Talalay, 1984). The Median-Effect Principle/Combination Index method leads to inefficiency by ignoring important sources of variation inherent in dose-response data and discarding data points that do not fit the Median-Effect Principle. Previous work has shown that the conventional method yields a high rate of false positives (Boik, Boik, Newman, 2008; Hennessey, Rosner, Bast, Chen, 2010) and, in some cases, low power to detect synergy. There is a great need for improving the current methodology. We developed a Bayesian framework for dose-response modeling and drug-drug interaction analysis. First, we developed a hierarchical meta-regression dose-response model that accounts for various sources of variation and uncertainty and allows one to incorporate knowledge from prior studies into the current analysis, thus offering a more efficient and reliable inference. Second, in the case that parametric dose-response models do not fit the data, we developed a practical and flexible nonparametric regression method for meta-analysis of independently repeated dose-response experiments. Third, and lastly, we developed a method, based on Loewe additivity that allows one to quantitatively assess interaction between two agents combined at a fixed dose ratio. The proposed method makes a comprehensive and honest account of uncertainty within drug interaction assessment. Extensive simulation studies show that the novel methodology improves the screening process of effective/synergistic agents and reduces the incidence of type I error. We consider an ovarian cancer cell line study that investigates the combined effect of DNA methylation inhibitors and histone deacetylation inhibitors in human ovarian cancer cell lines. The hypothesis is that the combination of DNA methylation inhibitors and histone deacetylation inhibitors will enhance antiproliferative activity in human ovarian cancer cell lines compared to treatment with each inhibitor alone. By applying the proposed Bayesian methodology, in vitro synergy was declared for DNA methylation inhibitor, 5-AZA-2'-deoxycytidine combined with one histone deacetylation inhibitor, suberoylanilide hydroxamic acid or trichostatin A in the cell lines HEY and SKOV3. This suggests potential new epigenetic therapies in cell growth inhibition of ovarian cancer cells.