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.

ABERRATIONS OF A PUTATIVE TUMOR SUPPRESSOR GENE SEL1L IN PANCREATIC DUCTAL ADENOCARCINOMA

Introduction: Pancreatic cancer is the fourth leading cause of cancer-related death among males and females in the United States. Sel-1-like (SEL1L) is a putative tumor suppressor gene that is downregulated in a significant proportion of human pancreatic ductal adenocarcinoma (PDAC). It was hypothesized that SEL1L expression could be down-modulated by somatic mutation, loss of heterozygosity (LOH), CpG island hypermethylation and/or aberrantly expressed microRNAs (miRNAs). Material and methods: In 42 PDAC tumors, the SEL1L coding region was amplified using reverse transcription polymerase chain reaction (RT-PCR), and analyzed by agarose gel electrophoresis and sequenced to search for mutations. Using fluorescent fragment analysis, two intragenic microsatellites in the SEL1L gene region were examined to detect LOH in a total of 73 pairs of PDAC tumors and normal-appearing adjacent tissues. Bisulfite DNA sequencing was performed to determine the methylation status of the SEL1L promoter in 41 PDAC tumors and 6 PDAC cell lines. Using real-time quantitative PCR, the expression levels of SEL1L mRNA and 7 aberrantly upregulated miRNAs that potentially target SEL1L were assessed in 42 PDAC tumor and normal pairs. Statistical methods were applied to evaluate the correlation between SEL1L mRNA and the miRNAs. Further the interaction was determined by functional analysis using a molecular biological approach. Results: No mutations were detected in the SEL1L coding region. More than 50% of the samples displayed abnormally alternate or aberrant spliced transcripts of SEL1L. About 14.5% of the tumors displayed LOH at the CAR/CAL microsatellite locus and 10.7% at the RepIN20 microsatellite locus. However, the presence of LOH did not show significant association with SEL1L downregulation. No methylation was observed in the SEL1L promoter. Statistical analysis showed that SEL1L mRNA expression levels significantly and inversely correlated with the expression of hsa-mir-143, hsa-mir-155, and hsa-mir-223. Functional analysis indicated that hsa-mir-155 acted as a suppressor of SEL1L in PL18 and MDAPanc3 PDAC cell lines. Discussion: Evidence from these studies suggested that SEL1L was possibly downregulated by aberrantly upregulated miRNAs in PDAC. Future studies should be directed towards developing a better understanding of the mechanisms for generation of aberrant SEL1L transcripts, and further analysis of miRNAs that may downregulate SEL1L.

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.

Human TNF receptor p75/80 gene structure and promoter characterization

Tumor necrosis factor receptor p75/80 ((TNF-R p75/80) is a 75 kDa type 1 transmembrane protein expressed predominately on cells of hematopoietic lineage. TNF-R p75/80 belongs to the TNF receptor superfamily characterized by cysteine-rich extracellular regions composed of three to six disulfide-linked domains. In the present report, we have characterized, for the first time, the complete gene structure for human TNF-R p75/80 which spans approximately 43 kbp. The gene consists of 10 exons (ranging from 34 bp to 2.5 kbp) and 9 introns (343 bp to 19 kbp). Consensus elements for transcription factors involved in T cell development and activation were noted in the 5$\sp\prime$ flanking region including TCF-1, Ikaros, AP-1, CK-2, IL-6RE, ISRE, GAS, NF-$\kappa$B and SP1, as well as an unusually high GC content and CpG frequency that appears characteristic of some TNF-R family members. The unusual (GATA)$\sb{\rm n}$ and (GAA)(GGA) repeats found within intron 1 may prove useful for further genome analysis within the 1p36 chromosomal locus. The human TNF-R p75/80 gene structure will permit further assessment of its involvement in normal hematopoietic cell development and function, autoimmune disease, and non-random translocations in hematopoietic malignancies. The region 1.8 kb 5$\sp\prime$ of the ATG was able to drive luciferase expression when transfected into cell lines expressing TNF-R p75/80. Further characterization of the 5$\sp\prime$-regulatory region will aid in determining factors and signal transduction pathways involved in regulating TNF-R p75/80 expression. ^

GLUCAGON GENE EXPRESSION: ANALYSIS OF PREPROGLUCAGON

Glucagon is a 29 amino acid polypeptide hormone produced in the (alpha) cells of the pancreatic islets. The purpose of this research was to understand better the role of glucagon in the regulation of metabolic processes. As with other polypeptide hormones, the synthesis of glucagon is thought to involve a larger precursor, which is then enzymatically cleaved to the functional form. The specific research objectives were to obtain cloned copies of the messenger RNA (mRNA) for pancreatic glucagon, to determine their primary sequences, and from this coding information to deduce the amino acid sequence of the initial glucagon precursor. From this suggested preproglucagon sequence and prior information on possible proglucagon intermediate processing products, the overall objective of this research is to propose a possible pathway for the biosynthesis of pancreatic glucagon.^ Synthetic oligodeoxynucleotide probes of 14-nucleotides (14-mer) and 17-nucleotides (a 17-mer) complementary to codons specifying a unique sequence of mature glucagon were synthesized. The ('32)P-labeled-14-mer was hybridized with size-fractionated fetal bovine pancreatic poly(A('+))RNA bound to nitrocellulose. RNA fractions of (TURN)14S were found to hybridize specifically, resulting in an (TURN)10-fold enrichment for these sequences. These poly(A('+))RNAs were translated in a cell-free system and the products analyzed by gel electrophoresis. The translation products were found to be enriched for a protein of the putative size of mammalian preproglucagon ((TURN)21 kd). These enriched RNA fractions were used to construct a complementary DNA (cDNA) library is plasmid pBR322.^ Screening of duplicate colony filters with the ('32)P-labeled-17-mer and a ('32)P-labeled-17-mer-primed cDNA probe indicated 25 possible glucagon clones from 3100 colonies screened. Restriction mapping of 6 of these clones suggested that they represented a single mRNA species. Primary sequence analysis of one clone containing a 1200 base pair DNA insert revealed that it contained essentially a full-length copy of glucagon cDNA.^ Analaysis of the cDNA suggested that it encoded an initial translation product of 180 amino acids with an M(,r) = 21 kd. The first initiation codon (ATG, methionine) followed by the longest open reading frame of 540 nucleotides was preceded by a 5'-untranslated region of 90 nucleotides, and was followed by a longer 3'-untranslated region of 471 nucleotides, resulting in a total of 1101 nucleotides. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

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

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

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

The role of plasmacytoid dendritic cells in the regulation of adaptive immunity

Plasmacytoid dendritic cells (pDCs) selectively express TLR7 which allows them to respond to RNA viruses and TLR9 which allows them to respond to DNA viruses and CpG oligonucleotides. Upon exposure to virus pDCs produce vast amounts of type I interferon (IFN) directly inhibiting viral replication and contributing to the activation of other immune cells. The ability of pDCs to promote B and T cell differentiation through type I IFN has been well documented although the role of additional factors including tumor necrosis factor (TNF) family members has not been thoroughly addressed. Here the expression of selected TNF family members in pDCs was examined and the role of TNF receptor-ligand interactions in the regulation of B and T lymphocyte growth and differentiation by pDCs was investigated. Upon stimulation with CpG-B, pDCs exhibit strong and stable expression of CD70, a TNF family ligand that binds to its receptor CD27 on memory B cells and promotes plasma cell differentiation and Ig secretion. Using an in vitro pDC/B cell co-culture system, it was determined that CpG-B-stimulated pDCs induce the proliferation of CD40L-activated human peripheral B cells and Ig secretion. This occurs independently of IFN and residual CpG, and requires physical contact between pDCs and B cells. CpG-stimulated pDCs induce the proliferation of both naive and memory B cells although Ig secretion is restricted to the memory subset. Blocking the interaction of CD70 with CD27 using an antagonist anti-CD70 antibody reduces the induction of B cell proliferation and IgG secretion by CpG-B-stimulated pDCs. Published studies have also indicated an important role for CD70 in promoting the expansion of CD4+ and CD8+ T cells and the development of effector function. CpG-B-stimulated pDCs induce naïve CD4+ T cell proliferation and production of multiple cytokines including IFN-γ, TNF-α, IL-10, IL-4, IL-5 and IL-13. Blocking the function of CD70 with an antagonist anti-CD70 antibody significantly reduced the induction of naïve CD4+ T cell proliferation by CpG-B-stimulated pDCs and the production of IL-4 and IL-13. Collectively these data indicate an important role for CD70 in the regulation of B and T lymphocyte growth and differentiation by pDCs. ^

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.

ALLERGEN SENSITIZATION AND IMMUNOMODULATION BY SYNTHETIC LIGANDS, PUL-042, IN AN ALLERGEN-INDUCED ASTHMA MURINE MODEL

Allergen-induced asthma is the leading form of asthma and a chronic condition worldwide. Common allergens are known to contribute to the pathogenesis of this disease. Murine models of allergic asthma have mostly used an intraperitoneal route of sensitization (not airway) to study this disease. Allergic asthma pathophysiology involves the activation of TH2-specific cells, which triggers production of IgE antibodies, the up-regulation of TH2-specific cytokines (i.e. IL-4, IL-5, IL-9 and IL-13), increased airway eosinophilia, and mucin hypersecretion. Although there are several therapeutics currently treating asthmatic patients, some of these treatments can result in drug tolerance and may be linked to increased mortality. CpG oligodeoxynucleotides (ODNs) is a synthetic ligand that targets Toll-like Receptor (TLR) 9. It has been evaluated as a therapeutic agent for the treatment of cancer, infectious diseases, and for treating allergy and asthma. PUL-042 is also a synthetic TLR ligand and is composed of two agonists against TLR2/6 heterodimer and TLR9. Previous studies have evaluated PUL-042 for its ability to confer resistance against bacterial and viral lung infection. These findings, combined with studies performed using CpG ODNs, led to speculation that PUL-042 dampens the immune response in allergen-induced asthma. My thesis research investigated airway route sensitization and airway delivery of PUL-042 to evaluate its effects in reducing an allergen-induced asthma phenotype in a murine model. The results of this study contribute to the foundation for future investigations to evaluate the efficacy of PUL-042 as a novel therapy in allergic-asthma disease.

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