High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.

Storage of bovine reproductive tissues and RNA extracts on ice for 24h or repeated freeze–thaw cycles do not affect RNA integrity

The aims of this study were to test (i) the effect of time of tissue and RNA extracts storage on ice and (ii) the effect of repeated freeze–thaw cycles on RNA integrity and gene expression of bovine reproductive tissues. Fragments of endometrium (ENDO), corpus luteum (CL) and ampulla (AMP) were subdivided and incubated for 0, 1, 3, 6, 12 or 24 h on ice. RNA extracts were incubated on ice for 0, 3, 12 or 24 h, or exposed to 1, 2, 4 or 6 freeze–thaw cycles. RNA integrity number (RIN) was estimated. Expression of progesterone receptor (PGR) and cyclophilin genes from RNA extracts stored on ice for 0 or 24 h, and 1 or 6 freeze–thaw cycles was measured by qPCR. Tissue and RNA extract incubation on ice, and repeated freeze–thaw cycles did not affect RIN values of RNA from ENDO, CL or AMP. Storage on ice or exposure to freeze–thaw cycles did not affect Cq values for PGR or cyclophilin genes. In conclusion, neither generalized RNA degradation nor specific RNA degradation was affected by storage of tissue or RNA extracts on ice for up to 24 h, or by up to 6 freeze–thaw cycles of RNA extracts obtained from bovine ENDO, CL and AMP.

The intronic long noncoding RNA ANRASSF1 recruits PRC2 to the RASSF1A promoter, reducing the expression of RASSF1A and increasing cell proliferation

The down-regulation of the tumor-suppressor gene RASSF1A has been shown to increase cell proliferation in several tumors. RASSF1A expression is regulated through epigenetic events involving the polycomb repressive complex 2 (PRC2); however, the molecular mechanisms modulating the recruitment of this epigenetic modifier to the RASSF1 locus remain largely unknown. Here, we identify and characterize ANRASSF1, an endogenous unspliced long noncoding RNA (lncRNA) that is transcribed from the opposite strand on the RASSF1 gene locus in several cell lines and tissues and binds PRC2. ANRASSF1 is transcribed through RNA polymerase II and is 5'-capped and polyadenylated; it exhibits nuclear localization and has a shorter half-life compared with other lncRNAs that bind PRC2. ANRASSF1 endogenous expression is higher in breast and prostate tumor cell lines compared with non-tumor, and an opposite pattern is observed for RASSF1A. ANRASSF1 ectopic overexpression reduces RASSF1A abundance and increases the proliferation of HeLa cells, whereas ANRASSF1 silencing causes the opposite effects. These changes in ANRASSF1 levels do not affect the RASSF1C isoform abundance. ANRASSF1 overexpression causes a marked increase in both PRC2 occupancy and histone H3K27me3 repressive marks, specifically at the RASSF1A promoter region. No effect of ANRASSF1 overexpression was detected on PRC2 occupancy and histone H3K27me3 at the promoter regions of RASSF1C and the four other neighboring genes, including two well-characterized tumor suppressor genes. Additionally, we demonstrated that ANRASSF1 forms an RNA/DNA hybrid and recruits PRC2 to the RASSF1A promoter. Together, these results demonstrate a novel mechanism of epigenetic repression of the RASSF1A tumor suppressor gene involving antisense unspliced lncRNA, in which ANRASSF1 selectively represses the expression of the RASSF1 isoform overlapping the antisense transcript in a location-specific manner. In a broader perspective, our findings suggest that other non-characterized unspliced intronic lncRNAs transcribed in the human genome might contribute to a location-specific epigenetic modulation of genes.

RNA Interference and cyclooxygenase-2 (COX-2) regulation in colon cancer cells

Despite new methods and combined strategies, conventional cancer chemotherapy still lacks specificity and induces drug resistance. Gene therapy can offer the potential to obtain the success in the clinical treatment of cancer and this can be achieved by replacing mutated tumour suppressor genes, inhibiting gene transcription, introducing new genes encoding for therapeutic products, or specifically silencing any given target gene. Concerning gene silencing, attention has recently shifted onto the RNA interference (RNAi) phenomenon. Gene silencing mediated by RNAi machinery is based on short RNA molecules, small interfering RNAs (siRNAs) and microRNAs (miRNAs), that are fully o partially homologous to the mRNA of the genes being silenced, respectively. On one hand, synthetic siRNAs appear as an important research tool to understand the function of a gene and the prospect of using siRNAs as potent and specific inhibitors of any target gene provides a new therapeutical approach for many untreatable diseases, particularly cancer. On the other hand, the discovery of the gene regulatory pathways mediated by miRNAs, offered to the research community new important perspectives for the comprehension of the physiological and, above all, the pathological mechanisms underlying the gene regulation. Indeed, changes in miRNAs expression have been identified in several types of neoplasia and it has also been proposed that the overexpression of genes in cancer cells may be due to the disruption of a control network in which relevant miRNA are implicated. For these reasons, I focused my research on a possible link between RNAi and the enzyme cyclooxygenase-2 (COX-2) in the field of colorectal cancer (CRC), since it has been established that the transition adenoma-adenocarcinoma and the progression of CRC depend on aberrant constitutive expression of COX-2 gene. In fact, overexpressed COX-2 is involved in the block of apoptosis, the stimulation of tumor-angiogenesis and promotes cell invasion, tumour growth and metastatization. On the basis of data reported in the literature, the first aim of my research was to develop an innovative and effective tool, based on the RNAi mechanism, able to silence strongly and specifically COX-2 expression in human colorectal cancer cell lines. In this study, I firstly show that an siRNA sequence directed against COX-2 mRNA (siCOX-2), potently downregulated COX-2 gene expression in human umbilical vein endothelial cells (HUVEC) and inhibited PMA-induced angiogenesis in vitro in a specific, non-toxic manner. Moreover, I found that the insertion of a specific cassette carrying anti-COX-2 shRNA sequence (shCOX-2, the precursor of siCOX-2 previously tested) into a viral vector (pSUPER.retro) greatly increased silencing potency in a colon cancer cell line (HT-29) without activating any interferon response. Phenotypically, COX-2 deficient HT-29 cells showed a significant impairment of their in vitro malignant behaviour. Thus, results reported here indicate an easy-to-use, powerful and high selective virus-based method to knockdown COX-2 gene in a stable and long-lasting manner, in colon cancer cells. Furthermore, they open up the possibility of an in vivo application of this anti-COX-2 retroviral vector, as therapeutic agent for human cancers overexpressing COX-2. In order to improve the tumour selectivity, pSUPER.retro vector was modified for the shCOX-2 expression cassette. The aim was to obtain a strong, specific transcription of shCOX-2 followed by COX-2 silencing mediated by siCOX-2 only in cancer cells. For this reason, H1 promoter in basic pSUPER.retro vector [pS(H1)] was substituted with the human Cox-2 promoter [pS(COX2)] and with a promoter containing repeated copies of the TCF binding element (TBE) [pS(TBE)]. These promoters were choosen because they are partculary activated in colon cancer cells. COX-2 was effectively silenced in HT-29 and HCA-7 colon cancer cells by using enhanced pS(COX2) and pS(TBE) vectors. In particular, an higher siCOX-2 production followed by a stronger inhibition of Cox-2 gene were achieved by using pS(TBE) vector, that represents not only the most effective, but also the most specific system to downregulate COX-2 in colon cancer cells. Because of the many limits that a retroviral therapy could have in a possible in vivo treatment of CRC, the next goal was to render the enhanced RNAi-mediate COX-2 silencing more suitable for this kind of application. Xiang and et al. (2006) demonstrated that it is possible to induce RNAi in mammalian cells after infection with engineered E. Coli strains expressing Inv and HlyA genes, which encode for two bacterial factors needed for successful transfer of shRNA in mammalian cells. This system, called “trans-kingdom” RNAi (tkRNAi) could represent an optimal approach for the treatment of colorectal cancer, since E. Coli in normally resident in human intestinal flora and could easily vehicled to the tumor tissue. For this reason, I tested the improved COX-2 silencing mediated by pS(COX2) and pS(TBE) vectors by using tkRNAi system. Results obtained in HT-29 and HCA-7 cell lines were in high agreement with data previously collected after the transfection of pS(COX2) and pS(TBE) vectors in the same cell lines. These findings suggest that tkRNAi system for COX-2 silencing, in particular mediated by pS(TBE) vector, could represent a promising tool for the treatment of colorectal cancer. Flanking the studies addressed to the setting-up of a RNAi-mediated therapeutical strategy, I proposed to get ahead with the comprehension of new molecular basis of human colorectal cancer. In particular, it is known that components of the miRNA/RNAi pathway may be altered during the progressive development of colorectal cancer (CRC), and it has been already demonstrated that some miRNAs work as tumor suppressors or oncomiRs in colon cancer. Thus, my hypothesis was that overexpressed COX-2 protein in colon cancer could be the result of decreased levels of one or more tumor suppressor miRNAs. In this thesis, I clearly show an inverse correlation between COX-2 expression and the human miR- 101(1) levels in colon cancer cell lines, tissues and metastases. I also demonstrate that the in vitro modulating of miR-101(1) expression in colon cancer cell lines leads to significant variations in COX-2 expression, and this phenomenon is based on a direct interaction between miR-101(1) and COX-2 mRNA. Moreover, I started to investigate miR-101(1) regulation in the hypoxic environment since adaptation to hypoxia is critical for tumor cell growth and survival and it is known that COX-2 can be induced directly by hypoxia-inducible factor 1 (HIF-1). Surprisingly, I observed that COX-2 overexpression induced by hypoxia is always coupled to a significant decrease of miR-101(1) levels in colon cancer cell lines, suggesting that miR-101(1) regulation could be involved in the adaption of cancer cells to the hypoxic environment that strongly characterize CRC tissues.

Mechanism of resistance to tyrosine kinase inhibitors in philadelphia-positive acute lymphblastic leukaemia (all): from genetic alterations to impaired RNA editing

The Ph chromosome is the most frequent cytogenetic aberration associated with adult ALL and it represents the single most significant adverse prognostic marker. Despite imatinib has led to significant improvements in the treatment of patients with Ph+ ALL, in the majority of cases resistance developed quickly and disease progressed. Some mechanisms of resistance have been widely described but the full knowledge of contributing factors, driving both the disease and resistance, remains to be defined. The observation of rapid development of lymphoblastic leukemia in mice expressing altered Ikaros (Ik) isoforms represented the background of this study. Ikaros is a zinc finger transcription factor required for normal hemopoietic differentiation and proliferation, particularly in the lymphoid lineages. By means of alternative splicing, Ikaros encodes several proteins that differ in their abilities to bind to a consensus DNA-binding site. Shorter, DNA nonbinding isoforms exert a dominant negative effect, inhibiting the ability of longer heterodimer partners to bind DNA. The differential expression pattern of Ik isoforms in Ph+ ALL patients was analyzed in order to determine if molecular abnormalities involving the Ik gene could associate with resistance to imatinib and dasatinib. Bone marrow and peripheral blood samples from 46 adult patients (median age 55 yrs, 18-76) with Ph+ ALL at diagnosis and during treatment with imatinib (16 pts) or dasatinib (30 pts) were collected. We set up a fast, high-throughput method based on capillary electrophoresis technology to detect and quantify splice variants. 41% Ph+ ALL patients expressed high levels of the non DNA-binding dominant negative Ik6 isoform lacking critical N-terminal zinc-fingers which display abnormal subcellular compartmentalization pattern. Nuclear extracts from patients expressed Ik6 failed to bind DNA in mobility shift assay using a DNA probe containing an Ikaros-specific DNA binding sequence. In 59% Ph+ ALL patients there was the coexistence in the same PCR sample and at the same time of many splice variants corresponded to Ik1, Ik2, Ik4, Ik4A, Ik5A, Ik6, Ik6 and Ik8 isoforms. In these patients aberrant full-length Ikaros isoforms in Ph+ ALL characterized by a 60-bp insertion immediately downstream of exon 3 and a recurring 30-bp in-frame deletion at the end of exon 7 involving most frequently the Ik2, Ik4 isoforms were also identified. Both the insertion and deletion were due to the selection of alternative splice donor and acceptor sites. The molecular monitoring of minimal residual disease showed for the first time in vivo that the Ik6 expression strongly correlated with the BCR-ABL transcript levels suggesting that this alteration could depend on the Bcr-Abl activity. Patient-derived leukaemia cells expressed dominant-negative Ik6 at diagnosis and at the time of relapse, but never during remission. In order to mechanistically demonstrated whether in vitro the overexpression of Ik6 impairs the response to tyrosine kinase inhibitors (TKIs) and contributes to resistance, an imatinib-sensitive Ik6-negative Ph+ ALL cell line (SUP-B15) was transfected with the complete Ik6 DNA coding sequence. The expression of Ik6 strongly increased proliferation and inhibited apoptosis in TKI sensitive cells establishing a previously unknown link between specific molecular defects that involve the Ikaros gene and the resistance to TKIs in Ph+ ALL patients. Amplification and genomic sequence analysis of the exon splice junction regions showed the presence of 2 single nucleotide polymorphisms (SNPs): rs10251980 [A/G] in the exon2/3 splice junction and of rs10262731 [A/G] in the exon 7/8 splice junction in 50% and 36% of patients, respectively. A variant of the rs11329346 [-/C], in 16% of patients was also found. Other two different single nucleotide substitutions not recognized as SNP were observed. Some mutations were predicted by computational analyses (RESCUE approach) to alter cis-splicing elements. In conclusion, these findings demonstrated that the post-transcriptional regulation of alternative splicing of Ikaros gene is defective in the majority of Ph+ ALL patients treated with TKIs. The overexpression of Ik6 blocking B-cell differentiation could contribute to resistance opening a time frame, during which leukaemia cells acquire secondary transforming events that confer definitive resistance to imatinib and dasatinib.

Variation in peach fruit susceptibility to Monilinia laxa and gene expression

Brown rot caused by Monilinia laxa and Monilinia fructigena is considered one of the most important diseases affecting Prunus species. Although some losses can result from the rotten fruits in the orchard, most of the damage is caused to fruits during the post-harvest phase. Several studies reported that brown rot incidence during fruit development highly varies; it was found that at a period corresponding to the the pit hardening stage, fruit susceptibility drastically decreases, to be quickly restored afterwards. However the molecular basis of this phenomenon is still not well understood. Furthermore, no difference in the rot incidence was found between wound and un-wound fruits, suggesting that resistance associated more to a specifc biochemical response of the fruit, rather than to a higher mechanical resistance. So far, the interaction Monilinia-peach was analyzed through chemical approaches. In this study, a bio-molecular approach was undertaken in order to reveal alteration in gene expression associated to the variation of susceptibility. In this thesis three different methods for gene expression analysis were used to analyze the alterations in gene expression occurring in peach fruits during the pit hardening stage, in a period encompassing the temporary change in Monilinia susceptibility: real time PCR, microarray and cDNA AFLP techniques. In 2005, peach fruits (cv.K2) were weekly harvested during a 19-week long-period, starting from the fourth week after full bloom, until full maturity. At each sampling time, three replicates of 5 fruits each were dipped in the M.laxa conidial suspension or in distilled water, as negative control. The fruits were maintained at room temperature for 3 hours; afterwards, they were peeled with a scalpel; the peel was immediately frozen in liquid nitrogen and transferred to -80 °C until use. The degree of susceptibility of peach fruit to the pathogen was determined on 3 replicates of 20 fruits each, as percentage of infected fruits, after one week at 20 °C. Real time PCR analysis was performed to study the variation in expression of those genes encoding for the enzymes of the phenylpropanoid pathway (phenylalanine ammonia lyase (PAL), chalcone synthase (CHS), cinnamate 4-hydroxylase (C4H), leucoanthocyanidine reductase (LAR), hydroxycinnamoyl CoA quinate hydroxycinnamoyl transferase (HQT) and of the jasmonate pathway, such as lipoxygenase (LOX), both involved in the production of important defense compounds. Alteration in gene expression was monitored on fruit samples of a period encompassing the pit hardening stage and the corresponding temporary resistance to M.laxa infections, weekly, from the 6thto the 12th week after full bloom (AFB) inoculated with M. laxa or mock-inoculated. The data suggest a critical change in the expression level of the phenylpropanoid pathway from the 7th to the 8th week AFB; such change could be directly physiologically associated to the peach growth and it could indirectly determine the decrease of susceptibility of peach fruit to Monilinia rot during the subsequent weeks. To investigate on the transcriptome variation underneath the temporary loss of susceptibility of peach fruits to Monilinia rot, the microarray and the cDNA AFLP techniques were used. The samples harvested on the 8th week AFB (named S, for susceptible ones) and on the 12th week AFB (named R, for resistant ones) were compared, both inoculated or mock-inoculated. The microarray experiments were carried out at the University of Padua (Dept. of Environmental Agronomy and Crop Science), using the μPEACH1.0 microarray together with the suited protocols. The analysis showed that 30 genes (corresponding to the 0.6% of the total sequences (4806) contained in the μPeach1.0 microarray) were found up-regulated and 31 ( 0.6%) down regulated in RH vs. SH fruits. On the other hand, 20 genes (0.4%) were shown to be up-regulated and 13 (0.3%) down-regulated in the RI vs. SI fruit. No genes were found differentially expressed in the mock-inoculated resistant fruits (RH) vs. the inoculated resistant ones (RI). Among the up-regulated genes an ATP sulfurylase, an heat shock protein 70, the major allergen Pru P1, an harpin inducing protein and S-adenosylmethionine decarboxylase were found, conversely among the down-regulated ones, cinnamyl alcohol dehydrogenase, an histidine- containing phosphotransfer protein and the ferritin were found. The microarray experimental results and the data indirectly derived, were tested by Real Time PCR analysis. cDNA AFLP analysis was also performed on the same samples. 339 transcript derived fragments considered significant for Monilinia resistance, were selected, sequenced and classified. Genes potentially involved in cell rescue and defence were well represented (8%); several genes (12.1%) involved in the protein folding, post-transductional modification and genes (9.2%) involved in cellular transport were also found. A further 10.3% of genes were classified as involved in the metabolism of aminoacid, carbohydrate and fatty acid. On the other hand, genes involved in the protein synthesis (5.7%) and in signal transduction and communication (5.7%) were found. Among the most interesting genes found differentially expressed between susceptible and resistant fruits, genes encoding for pathogenesis related (PR) proteins were found. To investigate on the association of Monilinia resistance and PR biological function, the major allergen Pru P1 (GenBank accession AM493970) and its isoform (here named Pru P2), were expressed in heterologous system and in vitro assayed for their anti-microbial activity. The ribonuclease activity of the recombinant Pru P1 and Pru P2 proteins was assayed against peach total RNA. As the other PR10 proteins, they showed a ribonucleolytic activity, that could be important to contrast pathogen penetration. Moreover Pru P1 and Pru P2 recombinant proteins were checked for direct antimicrobial activity. No inhibitory effect of Pru P1 or Pru P2 was detected against the selected fungi.

Transcriptional functions of DNA Topoisomerases at a genome-wide scale and a single gene levels

The DNA topology is an important modifier of DNA functions. Torsional stress is generated when right handed DNA is either over- or underwound, producing structural deformations which drive or are driven by processes such as replication, transcription, recombination and repair. DNA topoisomerases are molecular machines that regulate the topological state of the DNA in the cell. These enzymes accomplish this task by either passing one strand of the DNA through a break in the opposing strand or by passing a region of the duplex from the same or a different molecule through a double-stranded cut generated in the DNA. Because of their ability to cut one or two strands of DNA they are also target for some of the most successful anticancer drugs used in standard combination therapies of human cancers. An effective anticancer drug is Camptothecin (CPT) that specifically targets DNA topoisomerase 1 (TOP 1). The research project of the present thesis has been focused on the role of human TOP 1 during transcription and on the transcriptional consequences associated with TOP 1 inhibition by CPT in human cell lines. Previous findings demonstrate that TOP 1 inhibition by CPT perturbs RNA polymerase (RNAP II) density at promoters and along transcribed genes suggesting an involvement of TOP 1 in RNAP II promoter proximal pausing site. Within the transcription cycle, promoter pausing is a fundamental step the importance of which has been well established as a means of coupling elongation to RNA maturation. By measuring nascent RNA transcripts bound to chromatin, we demonstrated that TOP 1 inhibition by CPT can enhance RNAP II escape from promoter proximal pausing site of the human Hypoxia Inducible Factor 1 (HIF-1) and c-MYC genes in a dose dependent manner. This effect is dependent from Cdk7/Cdk9 activities since it can be reversed by the kinases inhibitor DRB. Since CPT affects RNAP II by promoting the hyperphosphorylation of its Rpb1 subunit the findings suggest that TOP 1inhibition by CPT may increase the activity of Cdks which in turn phosphorylate the Rpb1 subunit of RNAP II enhancing its escape from pausing. Interestingly, the transcriptional consequences of CPT induced topological stress are wider than expected. CPT increased co-transcriptional splicing of exon1 and 2 and markedly affected alternative splicing at exon 11. Surprisingly despite its well-established transcription inhibitory activity, CPT can trigger the production of a novel long RNA (5’aHIF-1) antisense to the human HIF-1 mRNA and a known antisense RNA at the 3’ end of the gene, while decreasing mRNA levels. The effects require TOP 1 and are independent from CPT induced DNA damage. Thus, when the supercoiling imbalance promoted by CPT occurs at promoter, it may trigger deregulation of the RNAP II pausing, increased chromatin accessibility and activation/derepression of antisense transcripts in a Cdks dependent manner. A changed balance of antisense transcripts and mRNAs may regulate the activity of HIF-1 and contribute to the control of tumor progression After focusing our TOP 1 investigations at a single gene level, we have extended the study to the whole genome by developing the “Topo-Seq” approach which generates a map of genome-wide distribution of sites of TOP 1 activity sites in human cells. The preliminary data revealed that TOP 1 preferentially localizes at intragenic regions and in particular at 5’ and 3’ ends of genes. Surprisingly upon TOP 1 downregulation, which impairs protein expression by 80%, TOP 1 molecules are mostly localized around 3’ ends of genes, thus suggesting that its activity is essential at these regions and can be compensate at 5’ ends. The developed procedure is a pioneer tool for the detection of TOP 1 cleavage sites across the genome and can open the way to further investigations of the enzyme roles in different nuclear processes.

Comparative genomic sequencing analysis of a region in human chromosome 11p15.3, mouse chromosome 7 and analysis of the novel STK33, Stk33 gene

The comparative genomic sequence analysis of a region in human chromosome 11p15.3 and its homologous segment in mouse chromosome 7 between ST5 and LMO1 genes has been performed. 158,201 bases were sequenced in the mouse and compared with the syntenic region in human, partially available in the public databases. The analysed region exhibits the typical eukaryotic genomic structure and compared with the close neighbouring regions, strikingly reflexes the mosaic pattern distribution of (G+C) and repeats content despites its relative short size. Within this region the novel gene STK33 was discovered (Stk33 in the mouse), that codes for a serine/threonine kinase. The finding of this gene constitutes an excellent example of the strength of the comparative sequencing approach. Poor gene-predictions in the mouse genomic sequence were corrected and improved by the comparison with the unordered data from the human genomic sequence publicly available. Phylogenetical analysis suggests that STK33 belongs to the calcium/calmodulin-dependent protein kinases group and seems to be a novelty in the chordate lineage. The gene, as a whole, seems to evolve under purifying selection whereas some regions appear to be under strong positive selection. Both human and mouse versions of serine/threonine kinase 33, consists of seventeen exons highly conserved in the coding regions, particularly in those coding for the core protein kinase domain. Also the exon/intron structure in the coding regions of the gene is conserved between human and mouse. The existence and functionality of the gene is supported by the presence of entries in the EST databases and was in vivo fully confirmed by isolating specific transcripts from human uterus total RNA and from several mouse tissues. Strong evidence for alternative splicing was found, which may result in tissue-specific starting points of transcription and in some extent, different protein N-termini. RT-PCR and hybridisation experiments suggest that STK33/Stk33 is differentially expressed in a few tissues and in relative low levels. STK33 has been shown to be reproducibly down-regulated in tumor tissues, particularly in ovarian tumors. RNA in-situ hybridisation experiments using mouse Stk33-specific probes showed expression in dividing cells from lung and germinal epithelium and possibly also in macrophages from kidney and lungs. Preliminary experimentation with antibodies designed in this work, performed in parallel to the preparation of this manuscript, seems to confirm this expression pattern. The fact that the chromosomal region 11p15 in which STK33 is located may be associated with several human diseases including tumor development, suggest further investigation is necessary to establish the role of STK33 in human health.

Identifikation transkriptioneller Alterationen menschlicher Tumoren und Entwicklung eines neuen auf RNA-Transfer basierenden Verfahrens zur Bestimmung der Genfunktion

Zielvorgaben der vorliegenden Arbeit war die Identifikation neuer selektiv in Tumoren aktivierter Gene sowie die Entwicklung eines methodischen Prozesses, um die molekularen Effekte der fehlerhaften Aktivierung solcher Gene zu untersuchen. Für die erste Fragestellung haben wir zwei komplementäre Methoden entwickelt. Zum einen haben wir nach neuen Mitglieder der Cancer/Germline (CG) Familie von Genen gesucht, die bereits attraktive Zielstrukturen laufender Phase I/IIa Studien sind. Zu diesem Zweck wurde ein bioinformatischer Data Mining Ansatz generiert. Dieser führte zur erfolgreichen in silico Klonierung neuer CG Gene. Zur Identifikation von in Tumorzellen überexprimierten Genen nutzten wir einen cDNA Mikroarray mit 1152 ausgewählten Genen mit direkter oder indirekter tumorimmunologischer oder tumorbiologischer Relevanz. Die komparative transkriptionelle Untersuchung von humanen Tumor- und Normalgeweben mit diesem Array führte zur Wiederentdeckung bereits bekannter, aber auch zur Aufdeckung bisher nicht beschriebener tumor-assoziierter Transkriptionsveränderungen. Der zweite große Schwerpunkt dieser Arbeit war die Technologieentwicklung eines versatilen Prozesses zur Untersuchung von molekularen Effekten eines aberrant in Zellen exprimierten Gens. Zur Simulation dieser Situation stellten wir in vitro transkribierte RNA dieses Gens her und elektroporierten diese in Zielzellen. Transkriptionsanalysen solcher Transfektanden mit Affymetrix Oligonukleotid Mikroarray deckten auf gesamt-genomischer Ebene ganze Kaskaden konsekutiver, transkriptioneller Alterationen auf.

RNA non-codificanti indotti da danno al DNA regolano il fattore di trascrizione HIF-1α

Recenti analisi sull’intero trascrittoma hanno rivelato una estensiva trascrizione di RNA non codificanti (ncRNA), le quali funzioni sono tuttavia in gran parte sconosciute. In questo lavoro è stato dimostrato che alte dosi di camptotecina (CPT), un farmaco antitumorale inibitore della Top1, aumentano la trascrizione di due ncRNA antisenso in 5’ e 3’ (5'aHIF-1α e 3'aHIF-1α rispettivamente) al locus genico di HIF-1α e diminuiscono i livelli dell’mRNA di HIF-1α stesso. Gli effetti del trattamento sono Top1-dipendenti, mentre non dipendono dal danno al DNA alla forca di replicazione o dai checkpoint attivati dal danno al DNA. I ncRNA vengono attivati in risposta a diversi tipi di stress, il 5'aHIF-1α è lungo circa 10 kb e possiede sia il CAP in 5’ sia poliadenilazione in 3’ (in letteratura è noto che il 3'aHIF-1α è un trascritto di 1,7 kb, senza 5’CAP né poliadenilazione). Analisi di localizzazione intracellulare hanno dimostrato che entrambi sono trascritti nucleari. In particolare 5'aHIF-1α co-localizza con proteine del complesso del poro nucleare, suggerendo un suo possibile ruolo come mediatore degli scambi della membrana nucleare. È stata dimostrata inoltre la trascrizione dei due ncRNA in tessuti di tumore umano del rene, evidenziandone possibili ruoli nello sviluppo del cancro. È anche noto in letteratura che basse dosi di CPT in condizioni di ipossia diminuiscono i livelli di proteina di HIF-1α. Dopo aver dimostrato su diverse linee cellulari che i due ncRNA sopracitati non potessero essere implicati in tale effetto, abbiamo studiato le variazioni dell’intero miRnoma alle nuove condizioni sperimentali. In tal modo abbiamo scoperto che il miR-X sembra essere il mediatore molecolare dell’abbattimento di HIF-1α dopo trattamento con basse dosi di CPT in ipossia. Complessivamente, questi risultati suggeriscono che il fattore di trascrizione HIF-1α venga finemente regolato da RNA non-codificanti indotti da danno al DNA.

Ribosome Biogenesis and cell cycle regulation: Effect of RNA Polymerase III inhibition

In cycling cells positive stimuli like nutrient, growth factors and mitogens increase ribosome biogenesis rate and protein synthesis to ensure both growth and proliferation. In contrast, under stress situation, proliferating cells negatively modulate ribosome production to reduce protein synthesis and block cell cycle progression. The main strategy used by cycling cell to coordinate cell proliferation and ribosome biogenesis is to share regulatory elements, which participate directly in ribosome production and in cell cycle regulation. In fact, there is evidence that stimulation or inhibition of cell proliferation exerts direct effect on activity of the RNA polymerases controlling the ribosome biogenesis, while several alterations in normal ribosome biogenesis cause changes of the expression and the activity of the tumor suppressor p53, the main effector of cell cycle progression inhibition. The available data on the cross-talk between ribosome biogenesis and cell proliferation have been until now obtained in experimental model in which changes in ribosome biogenesis were obtained either by reducing the activity of the RNA polymerase I or by down-regulating the expression of the ribosomal proteins. The molecular pathways involved in the relationship between the effect of the inhibition of RNA polymerase III (Pol III) activity and cell cycle progression have been not yet investigated. In eukaryotes, RNA Polymerase III is responsible for transcription of factors involved both in ribosome assembly (5S rRNA) and rRNA processing (RNAse P and MRP).Thus, the aim of this study is characterize the effects of the down-regulation of RNA Polymerase III activity, or the specific depletion of 5S rRNA. The results that will be obtained might lead to a deeper understanding of the molecular pathway that controls the coordination between ribosome biogenesis and cell cycle, and might give useful information about the possibility to target RNA Polymerase III for cancer treatment.

Phylogenomische Analysen bei Metazoen - Zur Stellung der Xenoturbellida und Syndermata

Im Rahmen der vorliegenden Dissertation wurde die phylogenetischen Stellungen der Xenoturbellida (Deuterostomia) und der Syndermata (Protostomia) mit phylogenomischen Techniken untersucht. Auf methodischer Ebene konnte gezeigt werden, dass ribosomale Proteine aufgrund ihres mittleren bis hohen Konservierungsgrades, ihrer Häufigkeit in kleineren EST-Projekten, damit verbunden ihrer Häufigkeit in Datenbanken und ihres phylogenetischen Informationsgehalts nützliche Werkzeuge für phylogenetische Fragestellungen sind. Es konnte durch phylogenetische Rekonstruktionen und Hypothesentests auf Basis eines 11.912 Aminosäuren langen Datensatzes gezeigt werden, dass die Xenoturbellida innerhalb der Deuterostomia eine Schwestergruppenbeziehung zu den Ambulacraria eingehen. Diese Arbeit zeigt im Vergleich aller bisher durchgeführten Arbeiten die beste statistische Unterstützung für diese Topologie. Weiterhin konnte untermauert werden, dass die Urochordata vermutlich anstelle der Cephalochordata die Schwestergruppe der Vertebrata sind. Der Vergleich der publizierten Xenoturbella EST-Datensätze mit dem eigenen Datensatz ließ den Rückschluß zu, dass ESTs offenbar klar weniger anfällig gegen Kontaminationen mit Erbmaterial (DNA+RNA) anderer Spezies sind als PCR-Amplifikate genomischer oder mitochondrialer Gene. Allerdings bestimmt anscheinend der physiologische Zustand der Tiere die Repräsentation von Transkriptklassen wie Stressproteine und mitochondriale Transkripte. Die bakteriellen Transkripte in einem der EST-Datensätze stammen vermutlich von Chlamydien, die möglicherweise symbiontisch in Xenoturbella bocki leben. Im Bereich der Protostomia wurden drei EST-Projekte für Vertreter der Syndermata durchgeführt. Basierend auf drei verschiedenen Proteinalignment-Datensätzen von ca. 11.000 Aminosäuren Länge konnte gezeigt werden, dass die Syndermata innerhalb der Spiralia einzugruppieren sind und dass sie mit den Gnathostomulida das monophyletische Supertaxon Gnathifera bilden. Die genaue phylogenetische Position der Syndermata innerhalb der Spiralia konnte hingegen noch nicht eindeutig geklärt werden, ebenso wie kein kongruenter Beweis für die Existenz des Supertaxons Platyzoa gefunden werden konnte. Im Rahmen der Untersuchung der internen Phylogenie der Syndermata konnten drei der fünf konkurrierenden Hypothesen aufgrund der Paraphylie der Eurotatoria ausgeschlossen werden. Da keine Daten der Seisonidea in den Analysen implementiert waren, bleibt die Frage der internen Phylogenie der Syndermata letztlich offen. Klar ist jedoch, dass die Eurotatoria nicht wie bislang angenommen monophyletisch sind, da die räderorgantragenden Bdelloidea keinesfalls den morphologisch diesbezüglich ähnlichen Monogononta ähnlich sind, sondern den räderorganlosen Acanthocephala näher stehen. Die Abbildung der molekularen Phylogenie auf die morphologischen Verhältnisse zeigt, dass das Räderorgan (partiell oder komplett) offenbar kurz nach der Aufspaltung der Syndermata in Monogononta und Acanthocephala + Bdelloidea in der Acanthocephala + Bdelloidea-Linie reduziert wurde. Die Entstehung des einziehbaren hinteren Körperteils (Rostrum bei Bdelloidea bzw. Proboscis bei Acanthocephala) in der Acanthocephala + Bdelloidea-Linie könnte das Schlüsselereignis zur Entstehung des Endoparasitismus der Acanthocephala gewesen sein.

Genexpressionsanalyse der fötalen Wachstumsfuge des Rindes

Erkrankungen des Skelettapparats wie beispielsweise die Osteoporose oder Arthrose gehören neben den Herz-Kreislauferkrankungen und Tumoren zu den Häufigsten Erkrankungen des Menschen. Ein besseres Verständnis der Bildung und des Erhalts von Knochen- oder Knorpelgewebe ist deshalb von besonderer Bedeutung. Viele bisherige Ansätze zur Identifizierung hierfür relevanter Gene, deren Produkte und Interaktionen beruhen auf der Untersuchung pathologischer Situationen. Daher ist die Funktion vieler Gene nur im Zusammenhang mit Krankheiten beschrieben. Untersuchungen, die die Genaktivität bei der Normalentwicklung von knochen- und knorpelbildenden Geweben zum Ziel haben, sind dagegen weit weniger oft durchgeführt worden. rnEines der entwicklungsphysiologisch interessantesten Gewebe ist die Epiphysenfuge der Röhrenknochen. In dieser sogenannten Wachstumsfuge ist insbesondere beim fötalen Gewebe eine sehr hohe Aktivität derjenigen Gene zu erwarten, die an der Knochen- und Knorpelbildung beteiligt sind. In der vorliegenden Arbeit wurde daher aus der Epiphysenfuge von Kälberknochen RNA isoliert und eine cDNA-Bibliothek konstruiert. Von dieser wurden ca. 4000 Klone im Rahmen eines klassischen EST-Projekts sequenziert. Durch die Analyse konnte ein ungefähr 900 Gene umfassendes Expressionsprofil erstellt werden und viele Transkripte für Komponenten der regulatorischen und strukturbildenden Bestandteile der Knochen- und Knorpelentwicklung identifiziert werden. Neben den typischen Genen für Komponenten der Knochenentwicklung sind auch deutlich Bestandteile für embryonale Entwicklungsprozesse vertreten. Zu ersten gehören in erster Linie die Kollagene, allen voran Kollagen II alpha 1, das mit Abstand höchst exprimierte Gen in der fötalen Wachstumsfuge. Nach den ribosomalen Proteinen stellen die Kollagene mit ca. 10 % aller auswertbaren Sequenzen die zweitgrößte Gengruppe im erstellten Expressionsprofil dar. Proteoglykane und andere niedrig exprimierte regulatorische Elemente, wie Transkriptionsfaktoren, konnten im EST-Projekt aufgrund der geringen Abdeckung nur in sehr geringer Kopienzahl gefunden werden. Allerdings förderte die EST-Analyse mehrere interessante, bisher nicht bekannte Transkripte zutage, die detaillierter untersucht wurden. Dazu gehören Transkripte die, die dem LOC618319 zugeordnet werden konnten. Neben den bisher beschriebenen drei Exonbereichen konnte ein weiteres Exon im 3‘-UTR identifiziert werden. Im abgeleiteten Protein, das mindestens 121 AS lang ist, wurden ein Signalpeptid und eine Transmembrandomäne nachgewiesen. In Verbindung mit einer möglichen Glykosylierung ist das Genprodukt in die Gruppe der Proteoglykane einzuordnen. Leicht abweichend von den typischen Strukturen knochen- und knorpelspezifischer Proteoglykane ist eine mögliche Funktion dieses Genprodukts bei der Interaktion mit Integrinen und der Zell-Zellinteraktion, aber auch bei der Signaltransduktion denkbar. rnDie EST-Sequenzierungen von ca. 4000 cDNA-Klonen können aber in der Regel nur einen Bruchteil der möglichen Transkripte des untersuchten Gewebes abdecken. Mit den neuen Sequenziertechnologien des „Next Generation Sequencing“ bestehen völlig neue Möglichkeiten, komplette Transkriptome mit sehr hoher Abdeckung zu sequenzieren und zu analysieren. Zur Unterstützung der EST-Daten und zur deutlichen Verbreiterung der Datenbasis wurde das Transkriptom der bovinen fötalen Wachstumsfuge sowohl mit Hilfe der Roche-454/FLX- als auch der Illumina-Solexa-Technologie sequenziert. Bei der Auswertung der ca. 40000 454- und 75 Millionen Illumina-Sequenzen wurden Verfahren zur allgemeinen Handhabung, der Qualitätskontrolle, dem „Clustern“, der Annotation und quantitativen Auswertung von großen Mengen an Sequenzdaten etabliert. Beim Vergleich der Hochdurchsatz Blast-Analysen im klassischen „Read-Count“-Ansatz mit dem erstellten EST-Expressionsprofil konnten gute Überstimmungen gezeigt werden. Abweichungen zwischen den einzelnen Methoden konnten nicht in allen Fällen methodisch erklärt werden. In einigen Fällen sind Korrelationen zwischen Transkriptlänge und „Read“-Verteilung zu erkennen. Obwohl schon simple Methoden wie die Normierung auf RPKM („reads per kilo base transkript per million mappable reads“) eine Verbesserung der Interpretation ermöglichen, konnten messtechnisch durch die Art der Sequenzierung bedingte systematische Fehler nicht immer ausgeräumt werden. Besonders wichtig ist daher die geeignete Normalisierung der Daten beim Vergleich verschieden generierter Datensätze. rnDie hier diskutierten Ergebnisse aus den verschiedenen Analysen zeigen die neuen Sequenziertechnologien als gute Ergänzung und potentiellen Ersatz für etablierte Methoden zur Genexpressionsanalyse.rn

Mechanistic study of the effects of 8-oxo-7,8-dihydroguanine (8-oxoG) on reporter gene expression in mammalian cells

DNA damage causes replication errors, leading to genetic instability or cell death. Besides that, many types of DNA base modifications have been shown to interfere with transcriptional elongation if they are located in the transcribed DNA strand of active genes, acting as roadblocks for RNA polymerases. It is widely assumed that transcription blockage by endogenous DNA damage is responsible for the early cell senescence in organs and accelerated ageing observed in individuals with compromised nucleotide excision repair.rnThe aims of this work were to design new experimental systems for testing transcription blocking potentials of DNA base modifications in an individual gene and to apply these test systems to the investigation of the effects of a frequent endogenously generated base modification, namely 8-oxo-7,8-hydroxyguanine (8-oxoG), on the gene transcription in cells. Several experimental strategies were employed for this purpose. First, I constructed an episomal vector encoding for a short-lived EGFP-ODC fusion protein and measured expression of the reporter gene in permanently transfected clonal cell lines exposed to DNA damaging agents. Second, the expression of plasmid-borne EGFP gene damaged with photosensitisers to obtain one or several oxidative purine modifications per plasmid molecule was determined in transiently transfected human and mouse host cells in an approach known as “host cell reactivation”. As a prerequisite for these experiments, a robust method of precise quantitative measurement of the EGFP gene expression in transiently transfected cells by flow cytometry was developed and validated. Third, I elaborated a very efficient procedure for insertion of synthetic oligonucleotides carrying 8-oxoG into plasmid DNA, avoiding any unwanted base damage and strand breaks. The consequences of 8-oxoG placed in defined positions in opposing DNA strands of the EGFP gene for transcription were measured by host cell reactivation in cells with functional 8-oxoguanine DNA glycosylase (OGG1) gene and in OGG1 null cells.rnThe results obtained in Ogg1-/- cells demonstrated that unrepaired 8-oxoG, even if situated in the transcribed DNA strand, does not have any negative effect on the reporter gene transcription. On the other hand, as few as one 8-oxoG was sufficient to cause a significant decrease of the gene expression in OGG1-proficient cell lines, i.e. in the presence of base excision repair. For two analysed positions of 8-oxoG in the plasmid DNA, the inhibition of gene transcription by the base modification correlated with the efficiency of its excision by purified OGG1 protein under cell-free conditions. Based on these findings, it has to be concluded that the observed decrease of transcription is mediated by excision of the base modification by OGG1 and probably caused by the repair-induced single-strand breaks. The mechanism of transcription inhibition by 8-oxoG is therefore clearly distinct from stalling of elongating RNA polymerase II complexes at the modified base.

Investigations into the effect of nucleoside modifications on the physicochemical properties and biological function of DNA and RNA

This thesis explores the effect of chemical nucleoside modification on the physicochemical and biological properties of nucleic acids. Positional alteration on the Watson-Crick edge of purines and pyrimidines, the “C-H” edge of pyrimidines, as well as both the Hoogsteen and sugar edges of purines were attempted by means of copper catalyzed azide-alkyne cycloaddition. For this purpose, nucleic acid building blocks carrying terminal alkynes were synthesized and introduced into oligonucleotides by solid-phase oligonucleotide chemistry. rnOf particular interest was the effect of nucleoside modification on hydrogen bond formation with complementary nucleosides. The attachment of propargyl functionalities onto the N2 of guanosine and the N4 of 5-methylcytosine, respectively, followed by incorporation of the modified analogs into oligonucleotides, was successfully achieved. Temperature dependent UV-absorption melting measurements with duplexes formed between modified oligonucleotides and a variety of complementary strands resulted in melting temperatures for the respective duplexes. As a result, the effect that both the nature and the site of nucleoside modification have on base pairing properties could thus be assisted. rnTo further explore the enzymatic recognition of chemically modified nucleosides, the oligonucleotide containing the N2-modified guanosine derivative on the 5’-end, which was clicked to a fluorescent dye, was subjected to knockdown analyses of the eGFP reporter gene in the presence of increasing concentrations of siRNA duplexes. From these dose-dependent experiments, a clear effect of 5’-labeling on the knockdown efficiency could be seen. In contrast, 3’-labeling was found to be relatively insignificant.rn