Involvement of microRNAs in Androgen Receptor-dependent Breast Cancers

Triple negative breast cancer (TNBC) is a very aggressive tumor subtype characterized by the lack of expression of estrogen receptor 1 (ESR1), due in the most of cases to an increased expression of DNA methyltransferases (DNMTs) and hypermethylation in CpG islands, resulting in gene silencing. Furthermore, in ESR1- negative breast cancers, androgen receptor (AR) is highly expressed and some studies suggest that it can drive tumor progression and might represent a therapeutic target. A correlation between microRNAs, small non-coding RNAs that regulate gene expression, and DNMTs was investigated in a TNBC cell line to restore a normal methylation pattern of ESR1, leading to its re-expression and conferring again sensitivity to selective estrogen receptor modulators (SERMs). miR-148A and miR-29B were found to be involved in the reduction of the expression of DNMT1 and DNMT3A and in a slight increase of ESR1 expression, but not at protein level. Then, we found a down-regulation of AR by miRs-7, -9, -27a, -27b, -29a, -29b, -29c, -127-3p, -127-5p and -376 at 48h post transfection and an up-regulation by miR-15a and miR-16 at every time considered. We concomitantly investigated a possible increase of Tamoxifen, Herceptin and Metformin sensitivity after AR silencing in MDA-MB 453 and T-47D cell lines. Cells seemed more sensitive when silenced for AR only in MDA-MB-453 at 24h post Tamoxifen treatment. Studies on Metformin have basically confirmed an increase of drug sensitivity due to AR silencing in both cell lines. Analysis of Herceptin showed how MDA-MB 453 samples silenced for AR have a slight decrease in the percentage of proliferating cells, demonstrating a possible increase in the response to treatment. These preliminary data provide the basis for further study of the modulation of the expression of AR by microRNAs and it will be interesting to understand the molecular mechanisms underlying these interactions.

Natural compounds Camptothecin and Triptolide: highly specific enzyme inhibitors and tools to dissect transcriptional functions

With this work I elucidated new and unexpected mechanisms of two strong and highly specific transcription inhibitors: Triptolide and Campthotecin. Triptolide (TPL) is a diterpene epoxide derived from the Chinese plant Trypterigium Wilfoordii Hook F. TPL inhibits the ATPase activity of XPB, a subunit of the general transcription factor TFIIH. In this thesis I found that degradation of Rbp1 (the largest subunit of RNA Polymerase II) caused by TPL treatments, is preceded by an hyperphosphorylation event at serine 5 of the carboxy-terminal domain (CTD) of Rbp1. This event is concomitant with a block of RNA Polymerase II at promoters of active genes. The enzyme responsible for Ser5 hyperphosphorylation event is CDK7. Notably, CDK7 downregulation rescued both Ser5 hyperphosphorylation and Rbp1 degradation triggered by TPL. Camptothecin (CPT), derived from the plant Camptotheca acuminata, specifically inhibits topoisomerase 1 (Top1). We first found that CPT induced antisense transcription at divergent CpG islands promoter. Interestingly, by immunofluorescence experiments, CPT was found to induce a burst of R loop structures (DNA/RNA hybrids) at nucleoli and mitochondria. We then decided to investigate the role of Top1 in R loop homeostasis through a short interfering RNA approach (RNAi). Using DNA/RNA immunoprecipitation techniques coupled to NGS I found that Top1 depletion induces an increase of R loops at a genome-wide level. We found that such increase occurs on the entire gene body. At a subset of loci R loops resulted particularly stressed after Top1 depletion: some of these genes showed the formation of new R loops structures, whereas other loci showed a reduction of R loops. Interestingly we found that new peaks usually appear at tandem or divergent genes in the entire gene body, while losses of R loop peaks seems to be a feature specific of 3’ end regions of convergent genes.

Statistical methods for the analysis of DNA sequences: application to dinucleotide distribution in the human genome

Questa tesi si inserisce nell'ambito delle analisi statistiche e dei metodi stocastici applicati all'analisi delle sequenze di DNA. Nello specifico il nostro lavoro è incentrato sullo studio del dinucleotide CG (CpG) all'interno del genoma umano, che si trova raggruppato in zone specifiche denominate CpG islands. Queste sono legate alla metilazione del DNA, un processo che riveste un ruolo fondamentale nella regolazione genica. La prima parte dello studio è dedicata a una caratterizzazione globale del contenuto e della distribuzione dei 16 diversi dinucleotidi all'interno del genoma umano: in particolare viene studiata la distribuzione delle distanze tra occorrenze successive dello stesso dinucleotide lungo la sequenza. I risultati vengono confrontati con diversi modelli nulli: sequenze random generate con catene di Markov di ordine zero (basate sulle frequenze relative dei nucleotidi) e uno (basate sulle probabilità di transizione tra diversi nucleotidi) e la distribuzione geometrica per le distanze. Da questa analisi le proprietà caratteristiche del dinucleotide CpG emergono chiaramente, sia dal confronto con gli altri dinucleotidi che con i modelli random. A seguito di questa prima parte abbiamo scelto di concentrare le successive analisi in zone di interesse biologico, studiando l’abbondanza e la distribuzione di CpG al loro interno (CpG islands, promotori e Lamina Associated Domains). Nei primi due casi si osserva un forte arricchimento nel contenuto di CpG, e la distribuzione delle distanze è spostata verso valori inferiori, indicando che questo dinucleotide è clusterizzato. All’interno delle LADs si trovano mediamente meno CpG e questi presentano distanze maggiori. Infine abbiamo adottato una rappresentazione a random walk del DNA, costruita in base al posizionamento dei dinucleotidi: il walk ottenuto presenta caratteristiche drasticamente diverse all’interno e all’esterno di zone annotate come CpG island. Riteniamo pertanto che metodi basati su questo approccio potrebbero essere sfruttati per migliorare l’individuazione di queste aree di interesse nel genoma umano e di altri organismi.

Epigenetic regulation of somatic angiotensin-converting enzyme by DNA methylation and histone acetylation

Somatic angiotensin-converting enzyme (sACE) is crucial in cardiovascular homeostasis and displays a tissue-specific profile. Epigenetic patterns modulate genes expression and their alterations were implied in pathologies including hypertension. However, the influence of DNA methylation and chromatin condensation state on the expression of sACE is unknown. We examined whether such epigenetic mechanisms could participate in the control of sACE expression in vitro and in vivo. We identified two CpG islands in the human ace-1 gene 3 kb proximal promoter region. Their methylation abolished the luciferase activity of ace-1 promoter/reporter constructs transfected into human liver (HepG2), colon (HT29), microvascular endothelial (HMEC-1) and lung (SUT) cell lines (p < 0.001). Bisulphite sequencing revealed a cell-type specific basal methylation pattern of the ace-1 gene -1,466/+25 region. As assessed by RT-qPCR, inhibition of DNA methylation by 5-aza-2'-deoxycytidine and/or of histone deacetylation by trichostatin A highly stimulated sACE mRNA expression cell-type specifically (p < 0.001 vs. vehicle treated cells). In the rat, in vivo 5-aza-cytidine injections demethylated the ace-1 promoter and increased sACE mRNA expression in the lungs and liver (p = 0.05), but not in the kidney. In conclusion, the expression level of somatic ACE is modulated by CpG-methylation and histone deacetylases inhibition. The basal methylation pattern of the promoter of the ace-1 gene is cell-type specific and correlates to sACE transcription. DNMT inhibition is associated with altered methylation of the ace-1 promoter and a cell-type and tissue-specific increase of sACE mRNA levels. This study indicates a strong influence of epigenetic mechanisms on sACE expression.

Distinct DNA methylation patterns in cirrhotic liver and hepatocellular carcinoma

Abberrant DNA methylation is one of the hallmarks of cancerogenesis. Our study aims to delineate differential DNA methylation in cirrhosis and hepatic cancerogenesis. Patterns of methylation of 27,578 individual CpG loci in 12 hepatocellular carcinomas (HCCs), 15 cirrhotic controls and 12 normal liver samples were investigated using an array-based technology. A supervised principal component analysis (PCA) revealed 167 hypomethylated loci and 100 hypermethylated loci in cirrhosis and HCC as compared to normal controls. Thus, these loci show a "cirrhotic" methylation pattern that is maintained in HCC. In pairwise supervised PCAs between normal liver, cirrhosis and HCC, eight loci were significantly changed in all analyses differentiating the three groups (p < 0.0001). Of these, five loci showed highest methylation levels in HCC and lowest in control tissue (LOC55908, CELSR1, CRMP1, GNRH2, ALOX12 and ANGPTL7), whereas two loci showed the opposite direction of change (SPRR3 and TNFSF15). Genes hypermethylated between normal liver to cirrhosis, which maintain this methylation pattern during the development of HCC, are depleted for CpG islands, high CpG content promoters and polycomb repressive complex 2 (PRC2) targets in embryonic stem cells. In contrast, genes selectively hypermethylated in HCC as compared to nonmalignant samples showed an enrichment of CpG islands, high CpG content promoters and PRC2 target genes (p < 0.0001). Cirrhosis and HCC show distinct patterns of differential methylation with regards to promoter structure, PRC2 targets and CpG islands.

HIC1 tumour suppressor gene is suppressed in acute myeloid leukaemia and induced during granulocytic differentiation

A hallmark of acute myeloid leukaemia (AML) is a block in differentiation caused by deregulated gene expression. The tumour suppressor Hypermethylated In Cancer 1 (HIC1) is a transcriptional repressor, which is epigenetically silenced in solid cancers. HIC1 mRNA expression was found to be low in 128 patient samples of AML and CD34+ progenitor cells when compared with terminally differentiated granulocytes. HIC1 mRNA was induced in a patient with t(15;17)-positive acute promyelocytic leukaemia receiving all-trans retinoic acid (ATRA) therapy. We therefore investigated whether HIC1 plays a role in granulocytic differentiation and whether loss of function of this gene might contribute to the differentiation block in AML. We evaluated HIC1 mRNA levels in HL-60 and U-937 cells upon ATRA-induced differentiation and in CD34+ progenitor cells after granulocyte colony-stimulating factor-induced differentiation. In both models of granulocytic differentiation, we observed significant HIC1 induction. When HIC1 mRNA was suppressed in HL-60 cells using stably expressed short hairpin RNA targeting HIC1, granulocytic differentiation was altered as assessed by CD11b expression. Bisulphite sequencing of GC-rich regions (CpG islands) in the HIC1 promoter provided evidence that the observed suppression in HL-60 cells was not because of promoter hypermethylation. Our findings indicate a role for the tumour suppressor gene HIC1 in granulocytic differentiation. Low expression of HIC1 may very well contribute to pathogenic events in leukaemogenesis.

Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

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.

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

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

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.

Construction of a 2.8-megabase yeast artificial chromosome contig and cloning of the human methylthioadenosine phosphorylase gene from the tumor suppressor region on 9p21.

Many human malignant cells lack methylthioadenosine phosphorylase (MTAP) enzyme activity. The gene (MTAP) encoding this enzyme was previously mapped to the short arm of chromosome 9, band p21-22, a region that is frequently deleted in multiple tumor types. To clone candidate tumor suppressor genes from the deleted region on 9p21-22, we have constructed a long-range physical map of 2.8 megabases for 9p21 by using overlapping yeast artificial chromosome and cosmid clones. This map includes the type IIFN gene cluster, the recently identified candidate tumor suppressor genes CDKN2 (p16INK4A) and CDKN2B (p15INK4B), and several CpG islands. In addition, we have identified other transcription units within the yeast artificial chromosome contig. Sequence analysis of a 2.5-kb cDNA clone isolated from a CpG island that maps between the IFN genes and CDKN2 reveals a predicted open reading frame of 283 amino acids followed by 1302 nucleotides of 3' untranslated sequence. This gene is evolutionarily conserved and shows significant amino acid homologies to mouse and human purine nucleoside phosphorylases and to a hypothetical 25.8-kDa protein in the pet gene (coding for cytochrome bc1 complex) region of Rhodospirillum rubrum. The location, expression pattern, and nucleotide sequence of this gene suggest that it codes for the MTAP enzyme.

Identification of putative noncoding RNAs among the RIKEN mouse full-length cDNA collection

With the sequencing and annotation of genomes and transcriptomes of several eukaryotes, the importance of noncoding RNA (ncRNA)-RNA molecules that are not translated to protein products-has become more evident. A subclass of ncRNA transcripts are encoded by highly regulated, multi-exon, transcriptional units, are processed like typical protein-coding mRNAs and are increasingly implicated in regulation of many cellular functions in eukaryotes. This study describes the identification of candidate functional ncRNAs from among the RIKEN mouse full-length cDNA collection, which contains 60,770 sequences, by using a systematic computational filtering approach. We initially searched for previously reported ncRNAs and found nine murine ncRNAs and homologs of several previously described nonmouse ncRNAs. Through our computational approach to filter artifact-free clones that lack protein coding potential, we extracted 4280 transcripts as the largest-candidate set. Many clones in the set had EST hits, potential CpG islands surrounding the transcription start sites, and homologies with the human genome. This implies that many candidates are indeed transcribed in a regulated manner. Our results demonstrate that ncRNAs are a major functional subclass of processed transcripts in mammals.

Evolutionary turnover of mammalian transcription start sites

Alignments of homologous genomic sequences are widely used to identify functional genetic elements and study their evolution. Most studies tacitly equate homology of functional elements with sequence homology. This assumption is violated by the phenomenon of turnover, in which functionally equivalent elements reside at locations that are nonorthologous at the sequence level. Turnover has been demonstrated previously for transcription-factor-binding sites. Here, we show that transcription start sites of equivalent genes do not always reside at equivalent locations in the human and mouse genomes. We also identify two types of partial turnover, illustrating evolutionary pathways that could lead to complete turnover. These findings suggest that the signals encoding transcription start sites are highly flexible and evolvable, and have cautionary implications for the use of sequence-level conservation to detect gene regulatory elements.

GONOME: measuring correlations between GO terms and genomic positions

Background: Current methods to find significantly under- and over-represented gene ontology (GO) terms in a set of genes consider the genes as equally probable balls in a bag, as may be appropriate for transcripts in micro-array data. However, due to the varying length of genes and intergenic regions, that approach is inappropriate for deciding if any GO terms are correlated with a set of genomic positions. Results: We present an algorithm - GONOME - that can determine which GO terms are significantly associated with a set of genomic positions given a genome annotated with (at least) the starts and ends of genes. We show that certain GO terms may appear to be significantly associated with a set of randomly chosen positions in the human genome if gene lengths are not considered, and that these same terms have been reported as significantly over-represented in a number of recent papers. This apparent over-representation disappears when gene lengths are considered, as GONOME does. For example, we show that, when gene length is taken into account, the term development is not significantly enriched in genes associated with human CpG islands, in contradiction to a previous report. We further demonstrate the efficacy of GONOME by showing that occurrences of the proteosome-associated control element (PACE) upstream activating sequence in the S. cerevisiae genome associate significantly to appropriate GO terms. An extension of this approach yields a whole-genome motif discovery algorithm that allows identification of many other promoter sequences linked to different types of genes, including a large group of previously unknown motifs significantly associated with the terms 'translation' and 'translational elongation'. Conclusion: GONOME is an algorithm that correctly extracts over-represented GO terms from a set of genomic positions. By explicitly considering gene size, GONOME avoids a systematic bias toward GO terms linked to large genes. Inappropriate use of existing algorithms that do not take gene size into account has led to erroneous or suspect conclusions. Reciprocally GONOME may be used to identify new features in genomes that are significantly associated with particular categories of genes.