RNA binding characteristics and overall topology of the vaccinia virus polyadenylation complex

mRNA 3′ polyadenylation is central to mRNA biogenesis in prokaryotes and eukaryotes, and is implicated in numerous aspects of mRNA metabolism, including efficiency of mRNA export from the nucleus, message stability, and initiation of translation. However, due to the great complexity of the eukaryotic polyadenylation apparatus, the mechanisms of RNA 3 ′ end processing have remained elusive. Although the RNA processing reactions leading to polyadenylated messenger RNA have been studied in many systems, and much progress has been made, a complete understanding of the biochemistry of the poly(A) polymerase enzyme is still lacking. My research uses Vaccinia virus as a model system to gain a better understanding of this complicated polyadenylation process, which consist of RNA binding, catalysis and polymerase translocation. ^ Vaccinia virus replicates in the cytoplasm of its host cell, so it must employ its own poly(A) polymerase (PAP), a heterodimer of two virus encoded proteins, VP55 and VP39. VP55 is the catalytic subunit, adding 30 adenylates to a non-polyadenylated RNA in a rapid processive manner before abruptly changing to a slow, non-processive mode of adenylate addition and dissociating from the RNA. VP39 is the stimulatory subunit. It has no polyadenylation catalytic activity by itself, but when associated with VP55 it facilitates the semi-processive synthesis of tails several hundred adenylates in length. ^ Oligonucleotide selection and competition studies have shown that the heterodimer binds a minimal motif of (rU)2 (N)25 U, the “heterodimer binding motif”, within an oligonucleotide, and its primer selection for polyadenylation is base-type specific. ^ Crosslinking studies using photosensitive uridylate analogs show that within a VP55-VP39-primer ternary complex, VP55 comes into contact with all three required uridylates, while VP39 only contacts the downstream uridylate. Further studies, using a backbone-anchored photosensitive crosslinker show that both PAP subunits are in close proximity to the downstream −10 to −21 region of 50mer model primers containing the heterodimer binding motif. This equal crosslinking to both subunits suggests that the dimerization of VP55 and VP39 creates either a cleft or a channel between the two subunits through which this region of RNA passes. ^ Peptide mapping studies of VP39 covalently crosslinked to the oligonucleotide have identified residue R107 as the amino acid in close proximity to the −10 uridylate. This helps us project a conceptual model onto the known physical surface of this subunit. In the absence of any tertiary structural data for VP55, we have used a series of oligonucleotide selection assays, as well as crosslinking, nucleotide transfer assays, and gel shift assays to gain insight into the requirements for binding, polyadenylation and translocation. Collectively, these data allow us to put together a comprehensive model of the structure and function of the polyadenylation ternary complex consisting of VP39, VP55 and RNA. ^

Characterization of hLodestar: A CDC5L interacting protein

Lodestar, a Drosophila maternal-effect gene, is essential for proper chromosome segregation during embryonic mitosis. Mutations in lodestar cause chromatin bridging in anaphase, preventing the sister chromatids from fully separating and leaving chromatin tangled at the metaphase plate. Drosophila lodestar protein was originally identified, in purified fractions of Drosophila Kc cell nuclear extracts, by its ability to suppress the generation of long RNA polymerase II transcripts. The human homolog of this protein (hLodestar) was cloned and studied in comparison to the Drosophila lodestar activities. The results of these studies show, similar to the Drosophila protein, hLodestar has dsDNA-dependent ATPase and transcription termination activity in vitro. hLodestar has also been shown to release RNA polymerase I and II stalled at a cyclobutane thymine dimer. Lodestar belongs to the SNF2 family of proteins, which are members of the DExH/D helicase super-family. The SNF2 family of proteins are believed to play a critical role in altering protein-DNA interactions in a variety of cellular contexts. We have recently isolated a human cDNA (hLodestar) that shares significant homology to the Drosophila lodestar gene. The 4.6 kb clone contains an open reading frame of 1162 amino acids, and shares 55% similarity and 46% identity to the Drosophila Lodestar protein sequence. Our studies looking for hLodestar interacting proteins revealed an association with CDC5L in the yeast two-hybrid system and co-immunoprecipitation experiments. CDC5L has been well documented to be a component of the spliceosome. Our data suggests hLodestar is involved in splicing through in vitro assembly and splicing reactions, in addition to its association with spliceosomes purified from HeLa nuclear extract. Although many other members of the DExH/D helicase super-family have been linked to splicing, this is the first SNF2 family member to be implicated in the splicing reaction. ^

RNA polymerase II large subunit degradation induced by ecteinascidin 743: Molecular characterization and subsequent rational investigation of antitumor mechanisms in cancers with clinical response

Ecteinascidin 743 (Et-743), which is a novel DNA minor groove alkylator with a unique spectrum of antitumor activity, is currently being evaluated in phase II/III clinical trials. Although the precise molecular mechanisms responsible for the observed antitumor activity are poorly understood, recent data suggests that post-translational modifications of RNA polymerase II Large Subunit (RNAPII LS) may play a central role in the cellular response to this promising anticancer agent. The stalling of an actively transcribing RNAPII LS at Et-743-DNA adducts is the initial cellular signal for transcription-coupled nucleotide excision repair (TC-NER). In this manner, Et-743 poisons TC-NER and produces DNA single strand breaks. Et-743 also inhibits the transcription and RNAPII LS-mediated expression of selected genes. Because the poisoning of TC-NER and transcription inhibition are critical components of the molecular response to Et-743 treatment, we have investigated if changes in RNAPII LS contribute to the disruption of these two cellular pathways. In addition, we have studied changes in RNAPII LS in two tumors for which clinical responses were reported in phase I/II clinical trials: renal cell carcinoma and Ewing's sarcoma. Our results demonstrate that Et-743 induces degradation of the RNAPII LS that is dependent on active transcription, a functional 26S proteasome, and requires functional TC-NER, but not global genome repair. Additionally, we have provided the first experimental data indicating that degradation of RNAPII LS might lead to the inhibition of activated gene transcription. A set of studies performed in isogenic renal carcinoma cells deficient in von Hippel-Lindau protein, which is a ubiquitin-E3-ligase for RNAPII LS, confirmed the central role of RNAPII LS degradation in the sensitivity to Et-743. Finally, we have shown that RNAPII LS is also degraded in Ewing's sarcoma tumors following Et-743 treatment and provide data to suggest that this event plays a role in decreased expression of the Ewing's sarcoma oncoprotein, EWS-Fli1. Altogether, these data implicate degradation of RNAPII LS as a critical event following Et-743 exposure and suggest that the clinical activity observed in renal carcinoma and Ewing's sarcoma may be mediated by disruption of molecular pathways requiring a fully functional RNAPII LS. ^

Reactive oxygen species-dependent mechanisms of N-(4-hydroxyphenyl)retinamide (4HPR)-induced apoptosis in human carcinoma cells

4HPR is a synthetic retinoid that has shown chemopreventive and therapeutic efficacy against premalignant and malignant lesions including oral leukoplakia, ovarian and breast cancer, and neuroblastoma. 4HPR induces apoptosis in various cancer cells and production of reactive oxygen species (ROS) has been suggested as a possible cause underlying these effects. However, the mechanisms governing these effects by 4HPR are not fully elucidated. In this study, we explored the mechanisms of 4HPR-induced ROS increase and apoptosis in human cancer cells. ^ First, we identified genes modulated by 4HPR using oligonucleotide gene expression arrays and found that they fall into specific functional canonical pathways and gene networks using Ingenuity Pathways Analysis®. Further analysis has shown that 4HPR induced up-regulation of Endoplasmic Reticulum (ER)-related genes such as Heat shock proteins 70 and 90 and the transcriptional factor, GADD153. These findings were validated using quantitative real-time PCR. ^ Second, we found that 4HPR induced extensive ER stress evidenced by dilation of the ER and endoribonuclease-mediated splicing and activation of the transcriptional factor, XBP-1. In addition, 4HPR induced the up-regulation of various ER stress-related genes and their protein products, as well as cleavage and activation of the ER specific Caspase-4. Concomitantly with XBP-1 splicing, all of these effects were dependent on ROS generation by 4HPR. Furthermore, chemical inhibition and RNA interference studies revealed a novel pro-apoptotic role for HSP70/A1A in 4HPR-mediated apoptosis. ^ Third, we observed rapid activation of the small GTPase Rac by 4HPR which was upstream of ROS generation. Inhibition of Rac activity or silencing of its expression by RNA interference inhibited ROS generation and apoptosis induction by 4HPR. siRNA targeting PAK1 and expression of a dominant negative Rac, decreased 4HPR-mediated ROS generation, while expression of a constitutive active Rac increased basal and 4HPR-induced ROS generation and PARP cleavage. Furthermore, metastatic cancer cells exhibited higher Rac activation, ROS generation, and cell growth inhibition due to 4HPR exposure compared to their primary cancer cell counterparts. ^ These findings provide novel insights into 4HPR-mediated ROS generation and apoptosis induction and support the use of ROS inducing agents such as 4HPR against metastatic cancer cells. ^

Micro RNA expression alterations in patients with glioblastoma

Glioblastoma, also known as glioblastoma multiform or GBM, is the most common and most malignant primary brain tumor. The clinical history of patients with glioblastoma is short, usually less than 3 months in more than 50% of cases after diagnosis. Currently, the methods of glioblastoma treatment are chemotherapy, radiotherapy and surgery. Even with the more effective treatment options, patients with glioblastoma most likely have a median survival time of 10 to 12 months. It is necessary to seek other treatment methods, including gene-targeted treatment. The success of gene-targeted treatment depends critically on the knowledge of genes that may be the cause of, or contribute to disease. To establish a correlate between glioblastoma survival timeline and micro RNA expression alteration, a study of 91 glioblastoma patients was conducted at the University of Texas M. D. Anderson Cancer Center. These 91 glioblastoma patients were newly diagnosed from 2002 to 2007. Statistical analysis was conducted to test the association of miRNA expression alteration between long-term survival and short-term survival glioblastoma. The completion of this proposed study will provide a better understanding of the regulatory role of miRNA in glioblastoma progression.^

Identification of germ cell tumor modifier genes from the 129.MOLF-Chr19 consomic mouse strain

Naturally occurring genetic variants confer susceptibility to disease in the human population, including in testicular germ cell tumor development. Disease susceptibility loci for testicular germ cell tumors have been identified by genetic mapping in humans and mice. However, the identity of many of the susceptibility genes remains unclear. My study utilized a chromosome substitution strain, the 129.MOLF-Chr 19 (or M19 strain), to identify candidate testicular germ cell tumor susceptibility genes. Males of this strain have a high incidence of germ cell tumors in the testes. By forward genetic approaches, five susceptibility loci were fine-mapped and the genetic interactions were dissected. In addition, I identified three protein-coding genes and one micro-RNA as testicular tumor susceptibility genes by genomic screening. Using reverse genetic approaches, I verified one of the candidates, Splicing factor 1, as a modifier of testicular tumor. Deficiency of SF1 significantly reduces the incidence of testicular tumors in mice. This study highlights the advantage of the 129.MOLF-Chr 19 consomic strain in disease gene identification and validation. It also sets the stage to elucidate the molecular mechanisms of tumorigenesis in the testis. ^

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

Structure-function analysis of human Integrator subunit-4

Structure-function analysis of human Integrator subunit 4 Anupama Sataluri Advisor: Eric. J. Wagner, Ph.D. Uridine-rich small nuclear RNAs (U snRNA) are RNA Polymerase-II (RNAPII) transcripts that are ubiquitously expressed and are known to be essential for gene expression. snRNAs play a key role in mRNA splicing and in histone mRNA expression. Inaccurate snRNA biosynthesis can lead to diseases related to defective splicing and histone mRNA expression. Although the 3′ end formation mechanism and processing machinery of other RNAPII transcripts such as mRNA has been well studied, the mechanism of snRNA 3′ end processing has remained a mystery until the recent discovery of the machinery that mediates this process. In 2005, a complex of 14 subunits (the Integrator complex) associated with RNA Polymerase-II was discovered. The 14subunits were annotated Integrator 1-14 based on their size. The subunits of this complex together were found to facilitate 3′ end processing of snRNA. Identification of the Integrator complex propelled research in the direction of understanding the events of snRNA 3’end processing. Recent studies from our lab confirmed that Integrator subunit (IntS) 9 and 11 together perform the endonucleolytic cleavage of the nascent snRNA 3′ end to generate mature snRNA. However, the role of other members of the Integrator complex remains elusive. Current research in our lab is focused on deciphering the role of each subunit within the Integrator complex This work specifically focuses on elucidating the role of human Integrator subunit 4 (IntS4) and understanding how it facilitates the overall function of the complex. IntS4 has structural similarity with a protein called “Symplekin”, which is part of the mRNA 3’end processing machinery. Symplekin has been thoroughly researched in recent years and structure-function correlation studies in the context of mRNA 3’end processing have reported a scaffold function for Symplekin due to the presence of HEAT repeat motifs in its N-terminus. Based upon the structural similarity between IntS4 and Symplekin, we hypothesized that Integrator subunit 4 may be behaving as a Symplekin-like scaffold molecule that facilitates the interaction between other members of the Integrator Complex. To answer this question, the two important goals of this study were to: 1) identify the region of IntS4, which is important for snRNA 3′ end processing and 2) determine binding partners of IntS4 which promote its function as a scaffold. IntS4 structurally consists of a highly conserved N-terminus with 8 HEAT repeats, followed by a nonconserved C- terminus. A series of siRNA resistant N and C-terminus deletion constructs as well as specific point mutants within its N-terminal HEAT repeats were generated for human IntS4 and, utilizing a snRNA transcriptional readthrough GFP-reporter assay, we tested their ability to rescue misprocessing. This assay revealed a possible scaffold like property of IntS4. To probe IntS4 for interaction partners, we performed co-immunoprecipitation on nuclear extracts of IntS4 expressing stable cell lines and identified IntS3 and IntS5 among other Integrator subunits to be binding partners which facilitate the scaffold like function of hIntS4. These findings have established a critical role for IntS4 in snRNA 3′ end processing, identified that both its N and C termini are essential for its function, and mapped putative interaction domains with other Integrator subunits.

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.

Regulation of mRNA splicing in Moloney murine sarcoma virus mutants

Cells infected with the conditionally defective MuSVts110 mutant of Moloney murine sarcoma virus are transformed at 33$\sp\circ$C but appear morphologically normal at 39$\sp\circ$C. The molecular basis for this phenotype is as follows: MuSVts110 contains a 1487 nucleotide central deletion that has truncated the 3$\sp\prime$ end to the gag gene and the 5$\sp\prime$ end of the mos gene. The resulting gag-mos junction is out-of-frame and the v-mos protein is not expressed. At 33$\sp\circ$C or lower, a splicing event is activated such that a 431 base intron is removed to realign the gag and mos gene in-frame, allowing the expression of a transforming protein P85$\sp{gag-mos}$. Temperature-dependent splicing appeared to be an intrinsic property of MuSVts110 transcripts and not a general feature of pre-mRNA splicing in 6m2 cells since splicing activity of a heterologous transcript in the same cells did not vary with temperature. The possibility that the splice event was not temperature-sensitive, but that the accumulation of spliced transcript at the lower growth temperatures was due to its selective thermolability was ruled out as stability studies revealed that the relative turnover rates of the unspliced and spliced MuSVts110 transcripts were not affected by temperature.^ The consensus sequences containing the splice sites activated in the MuSVts110 mutant (5$\sp\prime$ gag and 3$\sp\prime$ mos) are present, but not utilized, in wild-type MuSV-124. To test the hypothesis that it was the reduction of the 1919 base intervening sequence in MuSV-124 to 431 bases in MuSVts110 which activated splicing, the identical 1487 base deletion was introduced into cloned wild-type MuSV-124 DNA to create the MuSVts110 equivalent, ts32.^ To examine conditions permissive for splicing, we assayed splice site activation in a series of MuSV-124 "intron-modification" mutants. Data suggest that splicing in wild-type MuSV-124 may be blocked due to the lack of a proximal branchpoint sequence, but can be activated by those intron mutations which reposition a branch site closer to the 3$\sp\prime$ splice site. (Abstract shortened with permission of author.) ^

Estrogen induced desensitization of c -{\it fos\/} messenger RNA expression {\it in vivo\/}

In this thesis, we investigated the regulation of the nuclear proto-oncogene, c-fos by estrogen in vivo. In the uterus, estrogen causes a rapid, dramatic and transient induction of c-fos mRNA and this occurs by transcriptional activation. We have discovered a previously unrecognized regulatory mechanism by which fos becomes desensitized to estrogen following the transient induction. We investigated three aspects of this desensitization: (1) the kinetics and general characteristics of the phenomenon; (2) the molecular mechanism of the desensitization; and (3) the relationship of desensitization to estrogen stimulated DNA synthesis. The desensitization occurs between 3-24 hours after initial hormonal stimulation and is reversible within 72 hours. The desensitization is not species specific, in that it occurs in both the rat and mouse. The desensitization also occurs in at least two estrogen responsive tissues, the uterus and vagina. The desensitization is not unique to c-fos, since both c-myc and c-jun show similar patterns of desensitization. However, the desensitization is not observed with creatine kinase B (CKB), indicating that not all estrogen inducible genes become desensitized. In the second general area, we determined the desensitization is at the transcriptional level. The desensitization is homologous, but not heterologous, since estrogen induction does not desensitize c-fos to other agents. Other studies show that the desensitization is not due to the lack of functional estrogen receptors. Taken together, these findings suggest that the desensitization occurs at the level of the estrogen responsive element. In the third major area, we demonstrated that the desensitization appears to be related to estrogen induced DNA synthesis. Support for this suggestion comes from the observation that short acting estrogens which induce fos, but not DNA synthesis, do not produce desensitization. ^

Inhibitor of cytochrome c messenger RNA -protein interaction in stimulated rat skeletal muscle

The purpose of this work was to examine the possible mechanisms for the regulation of cytochrome c gene expression in response to increased contractile activity in rat skeletal muscle. The working hypothesis was that increased contractile activity enhances cytochrome c gene expression through a cis-element. A 110% increase in cytochrome c mRNA concentration was observed in tibialis anterior (TA) muscle after 9 days of chronic stimulation. Similar difference (120%) exists between soleus (SO) muscle of higher contractile activity and white vastus lateralis (WV) muscle of lower contractile activity. These results suggest that the endogenous cytochrome c gene expression is regulated by contractile activity. Cytochrome c-reporter genes were injected into skeletal muscles to identify the cis-element that is responsible for the regulation. Although the data was inconclusive, part of it suggested the importance of the 3$\sp\prime$-untranslated region (3$\sp\prime$-UTR) in mediating the response to increased contractile activity.^ RNA gel mobility shift (GMSA) and ultraviolet (UV) cross-linking assays revealed specific RNA-protein interaction in a 50-nucleotide region of the 3$\sp\prime$-UTR in unstimulated TA muscle. Computer analysis predicted a stem-loop structure of 17 nucleotides, which provides a structural basis for RNA-protein interaction. These 17 nucleotides are 100% conserved among rat, mouse and human cytochrome c genes and their 13 pseudogenes, suggesting a functional role for this region. The RNA-protein interaction was significantly less in highly active SO muscle than in inactive WV muscle and was dramatically decreased in stimulated TA muscle due to a protein inhibitor(s) associated with ribosome. It is possible that cytochrome c mRNAs undergoing translation are subject to a compartmentalized regulatory influence.^ The conclusion from these results is that increases in contractile activity induce or activate a protein inhibitor(s) associated with ribosome in rat skeletal muscle. The inhibitor decreases RNA-protein interaction in the 3$\sp\prime$-UTR of cytochrome c mRNA, which may result in increased mRNA stability and/or translation. ^