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.

Regulation of Set1-mediated methylation of Dam1

Eukaryotic genomes exist within a dynamic structure named chromatin in which DNA is wrapped around an octamer of histones forming the nucleosome. Histones are modified by a range of posttranslational modifications including methylation, phosphorylation, and ubiquitination, which are integral to a range of DNA-templated processes including transcriptional regulation. A hallmark for transcriptional activity is methylation of histone H3 on lysine (K) 4 within active gene promoters. In S. cerevisiae, H3K4 methylation is mediated by Set1 within the COMPASS complex. Methylation requires prior ubiquitination of histone H2BK123 by the E2-E3 ligases Rad6 and Bre1, as well as the Paf1 transcriptional elongation complex. This regulatory pathway exemplifies cross-talk in trans between posttranslational modifications on distinct histone molecules. Set1 has an additional substrate in the kinetochore protein Dam1, which is methylated on K233. This methylation antagonizes phosphorylation of adjacent serines by the Ipl1 Aurora kinase. The discovery of a second Set1 substrate raised the question of how Set1 function is regulated at the kinetochore. I hypothesized that transcriptional regulatory factors essential for H3K4 methylation at gene promoters might also regulate Set1-mediated methylation of Dam1K233. Here I show that the regulatory factors essential for COMPASS activity at gene promoters is also indispensable for the methylation of Dam1K233. Deletion of members of the COMPASS complex leads to loss of Dam1K233 methylation. In addition, deletion of Rad6, Bre1, or members of the Paf1 complex abolishes Dam1 methylation. The role of Rad6 and Bre1 in Dam1 methylation is dependent on H2BK123 ubiquitination, as mutation of K123 within H2B results in complete loss of Dam1 methylation. Importantly, methylation of Dam1K233 is independent of transcription and occurs at the kinetochore. My results demonstrate that Set1-mediated methylation is regulated by a general pathway regardless of substrate that is composed of transcriptional regulatory factors functioning independently of transcription at the kinetochore. My data provide the first example of cross-talk in trans between modifications on a histone and a non-histone protein. Additionally, my results indicate that several factors previously thought to be required for Set1 function at gene promoters are more generally required for the catalytic activity of the COMPASS complex regardless of substrate or cellular process.

CROSSTALK BETWEEN R1175 METHYLATION AND Y1173 PHOSPHORYLATION NEGATIVELY MODULATES EGFR-MEDIATED ERK ACTIVATION

Post-translational protein modifications are critical regulators of protein functions as they expand the signaling potentials of the modified proteins, leading to diverse physiological consequences. Currently, increasing evidence suggests that protein methylation is as important as other post-translational modifications in the regulation of various biological processes. This drives us to ask whether methylation is involved in the EGFR (epidermal growth factor receptor) signaling, a biological process extensively regulated by multiple post-translational modifications including phosphorylation, glycosylation and ubiquitination. We found that EGFR R1175 is methylated by a protein arginine methyltransferase named PRMT5. During EGFR activation, PRMT5-mediated R1175 methylation specifically enhances EGF-induced EGFR autophosphorylation at Y1173 residue. This novel modification crosstalk increases SHP1 recruitment to EGFR and suppresses EGFR-mediated ERK activation, resulting in inhibition of cell proliferation, migration, and invasion of EGFR-expressing cells. Based on these findings, we provide the first link between arginine methylation and tyrosine phosphorylation and identify R1175 methylation as an inhibitory modification specifically against EGFR-mediated ERK activation.

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.

An Sp1/Sp3 site polymorphism associated with hypermethylation of the candidate tumor suppressor gene RIL in cancer

Gene silencing due to promoter methylation is an alternative to mutations and deletions, which inactivate tumor suppressor genes (TSG) in cancer. We identified RIL by Methylated CpG Island Amplification technique as a novel aberrantly methylated gene. RIL is expressed in normal tissues and maps to the 5q31 region, frequently deleted in leukemias. We found methylation of RIL in 55/80 (69%) cancer cell lines, with highest methylation in leukemia and colon. We also observed methylation in 46/80 (58%) primary tumors, whereas normal tissues showed substantially lower degrees of methylation. RIL expression was lost in 13/16 cancer cell lines and was restored by demethylating agent. Screening of 38 cell lines and 13 primary cancers by SSCP revealed no mutations in RIL, suggesting that methylation and LOH are the primary inactivation mechanisms. Stable transfection of RIL into colorectal cancer cells resulted in reduction in cell growth, clonogenicity, and increased apoptosis upon UVC treatment, suggesting that RIL is a good candidate TSG. ^ In searching for a cause of RIL hypermethylation, we identified a 12-bp polymorphic sequence around the transcription start site of the gene that creates a long allele containing 3CTC repeat. Evolutionary studies suggested that the long allele appeared late in evolution due to insertion. Using bisulfite sequencing, in cancers heterozygous for RIL, we found that the short allele is 4.4-fold more methylated than the long allele (P = 0.003). EMSA results suggested binding of factor(s) to the inserted region of the long allele, but not to the short. EMSA mutagenesis and competition studies, as well as supershifts using nuclear extracts or recombinant Sp1 strongly indicated that those DNA binding proteins are Sp1 and Sp3. Transient transfections of RIL allele-specific expression constructs showed less than 2-fold differences in luciferase activity, suggesting no major effects of the additional Sp1 site on transcription. However, stable transfection resulted in 3-fold lower levels of transcription from the short allele 60 days post-transfection, consistent with the concept that the polymorphic Sp1 site protects against time-dependent silencing. Thus, an insertional polymorphism in the RIL promoter creates an additional Sp1/Sp3 site, which appears to protect it from silencing and methylation in cancer. ^

Emerging roles for the double strand break repair protein 53BP1 in transcriptional regulation

Disruption of the mechanisms that regulate cell-cycle checkpoints, DNA repair, and apoptosis results in genomic instability and often leads to the development of cancer. In response to double stranded breaks (DSBs) as induced by ionizing radiation (IR), generated during DNA replication, or through immunoglobulin heavy chain (IgH) rearrangements in T and B cells of lymphoid origin, the protein kinases ATM and ATR are central players that activate signaling pathways leading to DSB repair. p53 binding protein 1 (53BP1) participates in the repair of DNA double stranded breaks (DSBs) where it is recruited to or near sites of DNA damage. In addition to its well established role in DSB repair, multiple lines of evidence implicate 53BP1 in transcription which stem from its initial discovery as a p53 binding protein in a yeast two-hybrid screen. However, the mechanisms behind the role of 53BP1 in these processes are not well understood. ^ 53BP1 possesses several motifs that are likely important for its role in DSB repair including two BRCA1 C-terminal repeats, tandem Tudor domains, and a variety of phosphorylation sites. In addition to these motifs, we identified a glycine and arginine rich region (GAR) upstream of the Tudor domains, a sequence that is oftentimes serves as a site for protein arginine methylation. The focus of this project was to characterize the methylation of 53BP1 and to evaluate how methylation influenced the role of 53BP1 as a tumor suppressor. ^ Using a variety of biochemical techniques, we demonstrated that 53BP1 is methylated by the PRMT1 methyltransferase in vivo. Moreover, GAR methylation occurs on arginine residues in an asymmetric manner. We further show that sequences upstream of the Tudor domains that do not include the GAR stretch are sufficient for 53BP1 oligomerization in vivo. While investigating the role of arginine methylation in 53BP1 function, we discovered that 53BP1 associates with proteins of the general transcription apparatus as well as to other factors implicated in coordinating transcription with chromatin function. Collectively, these data support a role for 53BP1 in regulating transcription and provide insight into the possible mechanisms by which this occurs. ^

Molecular mechanisms by which p53/LSD1 and ERalpha/Trim24 complexes mediate gene regulation

In this dissertation, I identify two molecular mechanisms by which transcription factors cooperate with their co-regulators to mediate gene regulation. In the first part, I demonstrate that p53 directly recruits LSD1, a histone demethylase, to AFP chromatin to demethylate methylated H3K4 and actively mediate transcription repression. Loss of p53 and LSD1 interaction at chromatin leads to derepression of AFP in hepatic cells. In the second part, I reveal that Trim24 functions as an important co-activator in ERα-mediated gene activation in response to estrogen stimulation. Trim24 is recruited by ligand-bound ERα to chromatin and stabilizes ERα-chromatin interactions by binding to histone H3 via its PHD finger, which preferentially recognizes unmethylated H3K4. ^

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. ^

TAZ AS A REGULATOR OF MESENCHYMAL TRANSFORMATION AND CLINICAL AGGRESSIVENESS IN GLIOMAS

Glioblastoma multiforme (GBM) is an aggressive, high grade brain tumor. Microarray studies have shown a subset of GBMs with a mesenchymal gene signature. This subset is associated with poor clinical outcome and resistance to treatment. To establish the molecular drivers of this mesenchymal transition, we correlated transcription factor expression to the mesenchymal signature and identified transcriptional co-activator with PDZ-binding motif (TAZ) to be highly associated with the mesenchymal shift. High TAZ expression correlated with worse clinical outcome and higher grade. These data led to the hypothesis that TAZ is critical to the mesenchymal transition and aggressive clinical behavior seen in GBM. We investigated the expression of TAZ, its binding partner TEAD, and the mesenchymal marker FN1 in human gliomas. Western analyses demonstrated increased expression of TAZ, TEAD4, and FN1 in GBM relative to lower grade gliomas. We also identified CpG islands in the TAZ promoter that are methylated in most lower grade gliomas, but not in GBMs. TAZ-methylated glioma stem cell (GSC) lines treated with a demethylation agent showed an increase in mRNA and protein TAZ expression; therefore, methylation may be another novel way TAZ is regulated since TAZ is epigenetically silenced in tumors with a better clinical outcome. To further characterize the role of TAZ in gliomagenesis, we stably silenced or over-expressed TAZ in GSCs. Silencing of TAZ decreased invasion, self-renewal, mesenchymal protein expression, and tumor-initiating capacity. Over-expression of TAZ led to an increase in invasion, mesenchymal protein expression, mesenchymal differentiation, and tumor-initiating ability. These actions are dependent on TAZ interacting with TEAD since all these effects were abrogated with TAZ could not bind to TEAD. We also show that TAZ and TEAD directly bind to mesenchymal gene promoters. Thus, TAZ-TEAD interaction is critically important in the mesenchymal shift and in the aggressive clinical behavior of GBM. We identified TAZ as a regulator of the mesenchymal transition in gliomas. TAZ could be used as a biomarker to both estimate prognosis and stratify patients into clinically relevant subgroups. Since mesenchymal transition is correlated to tumor aggressiveness, strategies to target and inhibit TAZ-TEAD and the downstream gene targets may be warranted in alternative treatment.

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. ^